Assessing Natural Gas Versus CO2 Potential Underground Storage Sites in Greece: A Pragmatic Approach †

Although subsurface traps have been regularly explored for hydrocarbon exploration, natural gas and CO2 storage has drawn industrial attention over the past few decades, thanks to the increasing demand for energy resources and the need for greenhouse gas mitigation. With only one depleted hydrocarbon field in Greece, saline aquifers, salt caverns and sedimentary basins ought to be evaluated in furtherance of the latter. Within this study the potential of the Greek subsurface for underground storage is discussed. An overview and re-evaluation of the so-far studied areas is implemented based on the available data. Lastly, a pragmatic approach for the storage potential in Greece was created, delineating gaps and risks in the already proposed sites. Based on the above details, a case study for CO2 storage is presented, which is relevant to the West Katakolo field saline aquifer.

Environmental implications of fluoride presence in drinking water along the southeastern of El Bajío Guanajuatese, Guanajuato, Mexico: sources and health effects

Chronic exposure from drinking water with naturally high concentrations of fluoride (F⁻) has serious health consequences in several regions across the world including north-central Mexico as Guanajuato State, where the rural population is particularly dependent on untreated groundwater pumped from wells that have natural F- concentrations higher than those allowed by national and international regulations. The contaminated aquifers in the area are usually located in fractured volcanic environment that interacts with sedimentary basins and have a carbonate basement. Few studies focused on identifying the origin and hydrogeochemical processes related to fluoride release and mobilization have been developed, and even fewer that quantify the natural content of F in the geological environment. In this study, an evaluation of fluoride in volcanic rocks collected from 11 sampling sites along the Sierra de Codornices (Guanajuato State, Central Mexico) was carried out. The fluoride content is disseminated in volcanic rocks and the highest contents were obtained in felsic rocks. According to results obtained of a sampling campaign of 32 wells in 2019 their statistical and hydrogeochemical evaluation suggest that F- mobilization in groundwater from Juventino Rosas and Villagran municipalities, is a product of volcanic glass dissolution, a process involved in alkaline desorption occurring on the surfaces of F-containing minerals, and possibly on ion exchange occurring in minerals and some clays or even in deep fluids enriched in F. All these processes may be accelerated by the geothermal characteristics of the groundwater in the study area. The hydrogeochemical results and the epidemiological survey conducted indicate that children and older adults of Praderas de la Venta are at risk of exposure to F- due to the high concentrations ingested over a long period of time, to the toxicity of the element and its ability to accumulate in the bones. Prolonged exposure to high concentrations increases the risk.

Two-End-Member Mixing in the Fluids Emitted From Mud Volcano Lei-Gong-Huo, Eastern Taiwan: Evidence From Sr Isotopes

Mud volcano is one of the most important conduits for deep seated materials to migrate upward in sedimentary basins, convergent margins, and subduction zones. Understanding their temporal and spatial characteristics and variations provides us the important information on fluid sources and chemical compositions at depth. Mud volcano Lei-Gong-Huo (MV LGH) is a unique mud volcano, which is located on the mélange formation lying on the andesitic volcanic arc. Fluids emitted from 46 mud pools in MV LGH in eastern Taiwan were sampled and their major trace constitutes as well as H, O, and Sr isotopes (87Sr/86Sr and δ88Sr) were measured. Major constitutes of the fluids are Cl−, Na, and Ca. Compared with seawater, LGH fluids have lower Cl−, δD, δ18O, Na/Cl, K/Cl, and Mg/Cl but higher Ca/Cl ratios, indicating water–rock interaction of igneous rock and the ancient seawater at the source region. This interpretation is further supported by Sr isotopes, which show low value of 87Sr/86Sr ratio down to 0.70708. The result of spatial distribution showing strong negative correlation between Na and Ca concentration as well as Ca and 87Sr/86Sr ratios indicates that two end-member mixing is the major chemical characteristic. The fluids interacting with igneous rock carry high Ca, high δ88Sr, low Na, and low 87Sr/86Sr ratio, while those interacting with sedimentary rock carry low Ca, low δ88Sr, high Na, and high 87Sr/86Sr ratio. The source from the igneous region dominates the eastern and southeastern parts of the mud pools while sedimentary source dominates the western and northwestern parts. Most mud pools show mixing behavior between the two sources. Some of the sedimentary-dominated mud pools reveal existence of residual ancient water as indicated by 87Sr/86Sr. The major factor to fractionate the stable Sr isotopes in LGH waters is the source lithology. In summary, fluids emitted by mud pools in LGH originate from two sources, which are water–rock interactions of igneous rock with the ancient seawater from the east and sedimentary rock from the west at depth, resulting from the complex geologic background of mélange formation.

Hydrocarbon potential and prospects for exploration of Eastern Arctic oil and gas deposits

Conditions for the formation of hydrocarbon systems and prospects for searching for accumulations of oil and gas in the waters of the Eastern Arctic are considered. Significant hydrocarbon potential is predicted in the sedimentary basins of this region. All known manifestations of oil hydrocarbons are installed on land adjacent to the south, as well as on the east of the shelf. The East Arctic waters are included in a single model in order to perform an adequate comparative analysis of the evolution of hydrocarbon systems. The purpose of the research was to build space-time digital models of sedimentary basins and hydrocarbon systems, and to quantify the volume of generation, migration, and accumulation of hydrocarbons for the main horizons of source rocks. To achieve this goal, a spatiotemporal numerical basin simulation was carried out, based on which the distribution of probable hydrocarbon systems was determined and further analyzed. Following to the data obtained the most probable HC accumulation zones and types of fluids contained in potential traps were predicted. Keywords: numerical space-time basin modeling; modeling of hydrocarbon systems; evidence of oil and gas presence; Eastern Arctic; elements of hydrocarbon systems; oil and gas reservoirs; migration; accumulation; perspective objects

Tracing fluid evolution in sedimentary basins with calcite geochemical, isotopic and U-Pb geochronological data: Implications for petroleum and mineral resource accumulation in the Nanpanjiang Basin, South China

Fluid migration in sedimentary basins enable mass and energy transport and play critical roles in geochemical and geodynamical evolution of sedimentary basins. Moreover, reconstructing sedimentary basin fluid evolution from the geological record aids in constraining the evolution of associated petroleum and mineralization systems. As a relict of fluid flow activity, calcite is often a record of fluid flow and therefore can be used to characterize the fluids responsible for its precipitation. Here we study the Nanpanjiang Basin in South China where petroleum reservoirs and Carlin-type gold deposits spatially coincide. Through in situ U-Pb dating and geochemical analysis (87Sr/86Sr, δ18OVienna standard mean ocean water, δ13CVienna Peedee belemnite, rare earth elements) of calcite, this work constrains the key times related to petroleum migration/accumulation and Carlin-type gold mineralization, defines the basin fluid evolution, and proposes a genetic model for petroleum accumulation and gold mineralization within the Nanpanjiang Basin. The U-Pb age (ca. 241.4 Ma) for the gray/black calcite related to bitumen indicates the petroleum migration/accumulation occurred during the Triassic. The U-Pb date (ca. 106−121 Ma) of the white calcite associated with the gold-bearing pyrite, realgar, and fluorite record the lower timing limit of the Carlin-type gold systems. The geochemical data suggest both calcite types are cogenetic but suffered complex evolution with the gray/black calcite precipitating under low temperatures related to the continuous basin burial and the white calcite affected by post formation alteration related to both hydrothermal and meteoric fluids. Combined with the regional tectonic history, the Early Triassic petroleum migration/accumulation and the Early Cretaceous secondary Carlin-type gold mineralization events are considered to be related to the collision between the Indo-China and South China blocks, and the subduction between the Paleo-Pacific and Eurasian plates, respectively.

The westernmost Late Miocene–Pliocene volcanic activity in the Vardar zone (North Macedonia)

Late Miocene to Pleistocene volcanism within the Vardar zone (North Macedonia) covers a large area, has a wide range in composition, and is largely connected to the tectonic evolution of the South Balkan extensional system, the northern part of the Aegean extensional regime. The onset of the scattered potassic to ultrapotassic volcanism south from the Scutari-Peć transverse zone occurred at ca. 8.0 Ma based on this study. Here, we focused on three volcanic centers located on deep structures or thrust faults along the western part of the Vardar zone, for which there is none to very little geochronological and geochemical data available. Pakoševo and Debrište localities are represented as small remnants of lava flows cropping out at the southern edge of Skopje basin and at the western edge of Tikveš basin, respectively. Šumovit Greben center is considered as part of the Kožuf-Voras volcanic system, and it is located on its westernmost side, at the southern edge of Mariovo basin, which is largely composed of volcaniclastic sediments. We present new eruption ages applying the unspiked Cassignol-Gillot K–Ar technique on groundmass, as well as petrological and geochemical data, supplemented with Sr and Nd isotopes to complement and better understand the Neogene-Pleistocene volcanism in the region. Eruption ages on these rocks interlayered between sedimentary formations allow to better constrain the evolution of those sedimentary basins. Rocks from the three volcanic centers belong to the high-K calc-alkaline–shoshonitic series based on their elevated K content. The oldest center amongst these three localities, as well as other Late Miocene centers within the region, is the trachyandesitic Debrište, which formed at ca. 8.0 Ma, and exhibits the highest Nd and lowest Sr isotopic ratios (0.512441–0.512535 and 0.706759–0.706753, respectively). The basaltic trachyandesite Pakoševo center formed at ca. 3.8 Ma and its Nd and Sr isotopic ratios (0.512260 and 0.709593, respectively) bear the strongest signature of crustal contamination. The rhyolitic Šumovit Greben center is a composite volcanic structure formed at ca. 3.0–2.7 Ma. Its youngest eruption unit has a slightly higher Nd and lower Sr isotopic ratios (0.512382 and 0.709208, respectively) representing a magma with a lesser extent of crustal assimilation than the other samples from this center. The overall trend through time in the Sr and Nd isotopic ratios of the Late Miocene to Pleistocene mafic volcanic centers in the region implies an increasing rate of metasomatism of the lithospheric mantle.

Contrasting transform and passive margin subsidence history and heat flow evolution: insights from 3D thermo-mechanical modelling

Transform and passive margins developed during the continental rifting and opening of oceanic basins are fundamental elements of plate tectonics. It has been suggested that inherited structures, plate divergence velocities and surface processes exert a first order control on the topographic and bathymetric evolution and thermal history of these margins and related sedimentary basins. Their complex spatial-temporal dynamics have remained controversial. Here, we conducted 3D magmatic-thermo-mechanical numerical experiments coupled with surface processes modelling to simulate the dynamics of continental rifting, continental transform fault zone formation and persistent oceanic transform faulting and zero-offset oceanic fracture zones development. Numerical modelling results allow to explain the first order observations from passive and transform margins, such as diachronous rifting, heat flow rise and cooling in individual depocenters and contrasting basin tectonics of extensional and transtensional origin. In addition, the models reproduce the rise of both marginal ridges and transform marginal plateaus and their interaction with erosion and sedimentation. Comparison of model results with observations from natural examples yield new insights into the tectono-sedimentary and thermal evolution of several key passive and transform continental margins worldwide.Supplementary material at https://doi.org/10.6084/m9.figshare.c.5756555

Evolution of rift systems and their fault networks in response to surface processes

Continental rifting is responsible for the generation of major sedimentary basins, both during rift inception and during the formation of rifted continental margins. Geophysical and field studies revealed that rifts feature complex networks of normal faults but the factors controlling fault network properties and their evolution are still matter of debate. Here, we employ high-resolution 2D geodynamic models (ASPECT) including two-way coupling to a surface processes code (FastScape) to conduct 12 models of major rift types that are exposed to various degrees of erosion and sedimentation. We further present a novel quantitative fault analysis toolbox (Fatbox), which allows us to isolate fault growth patterns, the number of faults, and their length and displacement throughout rift history. Our analysis reveals that rift fault networks may evolve through five major phases: 1) distributed deformation and coalescence, 2) fault system growth, 3) fault system decline and basinward localization, 4) rift migration, and 5) breakup. These phases can be correlated to distinct rifted margin domains. Models of asymmetric rifting suggest rift migration is facilitated through both ductile and brittle deformation within a weak exhumation channel that rotates subhorizontally and remains active at low angles. In sedimentation-starved settings, this channel satisfies the conditions for serpentinization. We find that surface processes are not only able to enhance strain localization and to increase fault longevity but that they also reduce the total length of the fault system, prolong rift phases and delay continental breakup.