Style and evolution of salt pillows and related structures in the northern part of the Northeast German Basin

The northern part of the Northeast German Basin contains a large number of Late Permian (Zechstein) salt pillows, whereas diapiric structures are almost completely absent. This lack of diapirs facilitated the study of early stages of salt movement in the basin. Salt pillows and related structures were investigated in terms of distribution, geometry and time of initiation of salt flow within the regional geological context. The primary Zechstein thickness in the study area was reconstructed to gain more insight into the relationship between the geometry of the salt layer and the style of the salt-related structures. In this study, no clear spatial relationship between the salt structures and basement faults has been found and the location of the salt structures in this area appears to be highly independent of the underlying structural grain. The overburden is affected by minor faulting. We propose that buckling of the overburden due to regional compression significantly contributed to the initiation of the Late Jurassic to Early Cretaceous salt structures in the basin. Reverse faulting of the Gardelegen and Haldensleben Faults is related to inversion tectonics and exerted a compression on the basin fill. During the deformation, the Late Permian salt layer acted as an efficient detachment and led to a marked decoupling of the Mesozoic overburden from the underlying pre-Zechstein rocks.     

Towards the establishment of a reliable proxy for the reactive surface area of smectite

Reactive surface area is one of the key parameters for studying the kinetics of mineral dissolution. The common practice in experimental kinetics is to normalize the dissolution rate to the surface area measured by the BET method. The relationship between BET surface area and the reactive surface area is not trivial in minerals such as smectites, which possess both internal and external surface areas, and in which the dissolution is controlled by the chemical attack on the edge surface. The present study examines two proxies for the reactive surface area of the Clay Mineral Society reference smectite SAz-1: BET surface area and the edge surface area measured using AFM.Since smectites are very microporous, their BET surface area is strongly influenced by the degassing procedure. It is demonstrated that outgassing the smectite powder at 135°C in a 15 mL min⁻¹ N₂ gas flow for at least 24 hours minimizes contribution from micropores to less than 11% of the BET surface area.Following dissolution experiments in solutions with a low electrolyte concentration, the BET surface area increased from 34 ± 2 m² g⁻¹ in raw SAz-1 to 127 ± 13 m² g⁻¹ in SAz-1 sample recovered from dissolution experiments. This increase in BET surface area is explained by a decrease in the average size of the smectite aggregates, and by an increase in microporosity due to the depletion in the major interlayer cation, i.e., Ca²⁺. As the BET surface area of the raw smectite sample includes considerably less microporosity compared to the BET surface area of the smectite recovered from dissolution experiments, the former is a better approximation of the external surface area of the dried sample powder.AFM measurements show that there is no correlation between the specific external surface area of the sample and its specific edge surface area. This observation is explained by the platy morphology of the smectite particle in which the specific external surface area depends linearly on the height reciprocal, whereas the specific edge surface area is independent of the particles height and depends linearly on the sum of the reciprocals of the length of the axes. Therefore, there is no reason to expect a correlation between the BET and the edge surface area. Our results show that the edge surface area (4.9 ± 0.7 m² g⁻¹) of the smectite particles cannot be predicted based on its external surface area (136 ± 20 m² g⁻¹). Therefore, the BET surface area cannot serve as a proxy for the reactive surface area. We suggest using AFM measurements of the specific edge surface area as an alternative proxy for the reactive surface area of smectite.     

Correlation between hydrogen isotope ratios of lipid biomarkers and sediment maturity

We investigated the influence of thermal maturity on the hydrogen isotope ratios of sedimentary hydrocarbons to prove that the isotope ratio of hydrocarbons mirrors paleoclimate signatures. δD values from n-alkanes and acyclic isoprenoids of two sediment sections (Kupferschiefer [KS], 258 Ma, and Posidonienschiefer [PS], 184 Ma) with different maturation history were investigated. Both covered thermal maturity from 0.48 to 1.3 R_c (vitrinite reflectance and reflectance calculated from MPI1). Sediment burial up to 4500 m caused thermal maturation of organic matter in the KS horizon from the Early Zechstein basin of Poland, whereas contact metamorphic thermal maturation originated in the Early Toarcian PS (Posidonienschiefer) of the North German Vlotho Massif. The δD values of the extracted n-alkanes positively correlate with thermal maturity in the KS (y = 56‰ × MPI1[x] − 160‰ [VSMOW]) and in the PS (y = 104‰ × MPI1[x] − 200‰ [VSMOW]). The δD values of isoprenoids (i.e., pristane, phytane) were even more enriched with increasing maturity (y = 179‰ × MPI1[x] − 341‰ [VSMOW] in the KS; y = 300‰ × MPI1[x] − 415‰ [VSMOW] in PS).These results explain why isotope ratios of n-alkanes and isoprenoids in mature sediments are generally enriched in D and do not have the expected isotopic difference between n-alkanes and isoprenoids of ∼190‰. Moreover, the correlation between sediment maturity parameters and δD values suggests that after correction the δD values of n-alkanes can be used to reconstruct climate and environment in the geological past.     

Numerical simulation of a tunnel surrounded by sand under earthquake using a hypoplastic model

A numerical model of a centrifuge experiment on tunnel located in sand is being presented. The experiment was carried out under seismic loading using a dynamic actuator. The responses of the tunnel and of the sand were measured. The numerical model is based on a hypoplastic constitutive model with intergranular strains implemented in the FE-code TOCHNOG. The calculated accelerations in the sand match the measured results, while the surface settlement and the bending moments in the tunnel lining are only qualitatively reproduced by the numerical model.     

Methane and sulfide fluxes in permanent anoxia: In situ studies at the Dvurechenskii mud volcano (Sorokin Trough, Black Sea)

The Dvurechenskii mud volcano (DMV) is located in permanently anoxic waters at 2060 m depth (Sorokin Trough, Black Sea). The DMV was studied during the RV Meteor expedition M72/2 as an example of an active mud volcano system, to investigate the significance of submarine mud volcanism for the methane and sulfide budget of the anoxic Black Sea hydrosphere. Our studies included benthic fluxes of methane and sulfide, as well as the factors controlling transport, consumption and production of both compounds within the sediment. The pie-shaped mud volcano showed temperature anomalies as well as solute and gas fluxes indicating high fluid flow at its summit north of the geographical center. The anaerobic oxidation of methane (AOM) coupled to sulfate reduction (SR) was repressed in this zone due to the upward flow of sulfate-depleted fluids through recently deposited subsurface muds, apparently limiting microbial methanotrophic activity. Consequently, the emission of dissolved methane into the water column was high, with an estimated rate of 0.46 mol m⁻² d⁻¹. On the wide plateau and edge of the mud volcano surrounding the summit, fluid flow and total methane flux were lower, allowing higher SR and AOM rates correlated with an increase in sulfate penetration into the sediment. Here, between 50% and 70% of the methane flux (0.07-0.1 mol m⁻² d⁻¹) was consumed within the upper 10 cm of the sediment. The overall amount of dissolved methane released from the entire mud volcano structure into the water column was significant with a discharge of 1.3 × 10⁷ mol yr⁻¹. The DMV maintains also high areal rates of methane-fueled sulfide production and emission of on average 0.05 mol m⁻² d⁻¹. This is a difference to mud volcanoes in oxic waters, which emit similar amounts of methane, but not sulfide. However, based on a comparison of this and other mud volcanoes of the Black Sea, we conclude that sulfide and methane emission into the hydrosphere from deep-water mud volcanoes does not significantly contribute to the sulfide and methane inventory of the Black Sea.     

New age constraints on the western Betic intramontane basins: A late Tortonian closure of the Guadalhorce Corridor?

Several gateways connected the Mediterranean with the Atlantic during the late Miocene but the timing of closure and therefore their role prior to and during the Messinian Salinity Crisis (5.97–5.33 Ma) is still under debate. The timing of closure of the Guadalhorce Corridor is especially disputed as the common lack of marine microfossils hampers precise age determination. Here we present new biostratigraphic age constraints on the sediments of the Ronda, Antequera and Arcos regions, which formed the northern part of the Guadalhorce Corridor. The general presence of Globorotalia menardii 4 in the youngest deep‐marine sediments of all three regions indicates a late Tortonian age, older than 7.51 Ma. We conclude that the Guadalhorce Corridor closed during the late Tortonian, well before the onset of the Messinian Salinity Crisis and that the late Tortonian tectonic uplift of the eastern Betics extended into the western Betics.     

Experimental and petrological constraints on local-scale interaction of biotite-amphibole gneiss with H2O-CO2-（K, Na）Cl fluids at middle-crustal conditions： Example from the Limpopo Complex, South Africa

Reaction textures and fluid inclusions in the～2.0 Ga pyroxene-bearing dehydration zones within the Sand River biotite-hornblende orthogneisses（Central Zone of the Limpopo Complex） suggest that the formation of these zones is a result of close interplay between dehydration process along ductile shear zones triggered by H₂O-CO₂-salt fluids at 750—800℃and 5.5—6.2 kbar.partial melting,and later exsolution of residual brine and H₂O-CO₂ fluids during melt crystallization at 650—700℃.These processes caused local variations of water and alkali activity in the fluids,resulting in various mineral assemblages within the dehydration zone.The petrological observations are substantiated by experiments on the interaction of the Sand River gneiss with the H₂O-CO-2-（K,Na）Cl fluids at 750 and 800℃and 5.5 kbar.It follows that the interaction of biotite-amphibole gneiss with H₂O-CO₂-（K.Na）Cl fluids is accompanied by partial melting at 750—800℃.Orthopyroxene-bearing assemblages are characteristic for temperature 800℃and are stable in equilibrium with fluids with low salt concentrations,while salt-rich fluids produce clinopyroxene-bearing assemblages.These observations arc in good agreement with the petrological data on the dehydration zones within the Sand River orthoeneisses.     

Adjoint formulation and constraint handling for gradient-based optimization of compositional reservoir flow

An adjoint formulation for the gradient-based optimization of oil–gas compositional reservoir simulation problems is presented. The method is implemented within an automatic differentiation-based compositional flow simulator (Stanford’s Automatic Differentiation-based General Purpose Research Simulator, AD-GPRS). The development of adjoint procedures for general compositional problems is much more challenging than for oil–water problems due to the increased complexity of the code and the underlying physics. The treatment of nonlinear constraints, an example of which is a maximum gas rate specification in injection or production wells, when the control variables are well bottom-hole pressures, poses a particular challenge. Two approaches for handling these constraints are presented—a formal treatment within the optimizer and a simpler heuristic treatment in the forward model. The relationship between discrete and continuous adjoint formulations is also elucidated. Results for four example cases of increasing complexity are presented. Improvements in the objective function (cumulative oil produced) relative to reference solutions range from 4.2 to 11.6 %. The heuristic treatment of nonlinear constraints is shown to offer a cost-effective means for obtaining feasible solutions, which are, in some cases, better than those obtained using the formal constraint handling procedure.