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.     

Computing derivative information of sequentially coupled subsurface models

A generic framework for the computation of derivative information required for gradient-based optimization using sequentially coupled subsurface simulation models is presented. The proposed approach allows for the computation of any derivative information with no modification of the mathematical framework. It only requires the forward model Jacobians and the objective function to be appropriately defined. The flexibility of the framework is demonstrated by its application in different reservoir management studies. The performance of the gradient computation strategy is demonstrated in a synthetic water-flooding model, where the forward model is constructed based on a sequentially coupled flow-transport system. The methodology is illustrated for a synthetic model, with different types of applications of data assimilation and life-cycle optimization. Results are compared with the classical fully coupled (FIM) forward simulation. Based on the presented numerical examples, it is demonstrated how, without any modifications of the basic framework, the solution of gradient-based optimization models can be obtained for any given set of coupled equations. The sequential derivative computation methods deliver similar results compared to FIM methods, while being computationally more efficient.     

The accessibility of nodes in the Belgian road network

This paper presents an analysis of the accessibility of nodes in the Belgian road network. The intersection of roads are considered as the nodes and their connecting roads as the edges of a network. By means of the Dijkstra algorithm accessibility values are computed for each node by summing shortest travel times to a reference system of 158 nodes spread evenly over the Belgian territory. A comparison of the measured actual accessibility with that determined from an optimal situation provides an indication of the quality of directness of the Belgian road system. Accessibility and directness are shown cartographically. An isochronal map illustrates the access to Brussels.     

Correction to: Leveraging HPC accelerator architectures with modern techniques — hydrologic modeling on GPUs with ParFlow

Computational Geosciences - A Correction to this paper has been published: https://doi.org/10.1007/s10596-021-10065-y     

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.     

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.     

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.     

The Effect of Surface Heterogeneity on the Temperature–Humidity Correlation and the Relative Transport Efficiency

We derive a conceptual model of the flow over heterogeneous terrain consisting of patches with contrasting Bowen ratios. Upward moving eddies are assumed to carry heterogeneous properties, whereas downward moving eddies carry homogeneous properties. This results in a decorrelation of temperature and humidity as the contrast between the patches increases. We show that this model is able to reproduce the relationship developed by Lamaud and Irvine (Boundary-Layer Meteorol. 120:87–109, 2006). Some details differ from their expression but are in accordance with data obtained over African savannah. We extend the conceptual model to a combination of any scalars, not necessarily linked through the surface energy balance (as is the case for temperature and humidity). To this end we introduce a new parameter that describes the surface heterogeneity in surface fluxes. The results of the current model can be used to predict the discrepancy between similarity relationships for different scalars over heterogeneous terrain.     

Products and Mechanism of the Gas Phase Reaction of Ozone with β-Pinene

Gas phase ozonolysis of -pinene was performedin a 570 l static reactor at 730 Torr and 296 K insynthetic air and the products were analysed by acombination of gas phase FTIR spectroscopy, HPLC andIC analyses of gas phase and aerosol samples,respectively. The reaction mechanism was investigatedby adding HCHO, HCOOH and H₂O as Criegeeintermediate scavenger and cyclohexane as OH radicalscavenger. Main identified products (yields inparentheses) in the presence of cyclohexane as OHradical scavenger were HCHO (0.65 ± 0.04),nopinone (0.16 ± 0.04), 3-hydroxy-nopinone (0.15± 0.05), CO₂ (0.20 ± 0.04), CO (0.030± 0.002), HCOOH (0.020 ± 0.002), the secondaryozonide of -pinene (0.16 ± 0.05), andcis-pinic acid (0.02 ± 0.01). The decompositionof the primary ozonide was found to yieldpredominantly the excited C₉-Criegee intermediateand HCHO (0.84 ± 0.04) and to a minor extent theexcited CH₂OO intermediate and nopinone (0.16± 0.04). Roughly 40% of the excitedC₉-Criegee intermediate becomes stabilised andcould be shown to react with HCHO, HCOOH and H₂O. The atmospherically important reaction of thestabilised C₉-Criegee intermediate with H₂Owas found to result in a nopinone increase of (0.35± 0.05) and in the formation of H₂O₂(0.24 ± 0.03). Based on the observed products,the unimolecular decomposition/isomerisationchannels of the C₉-Criegee intermediate arediscussed in terms of the hydroperoxide and esterchannels. Subsequent reactions of the nopinonylradical, formed in the hydroperoxide channel, lead tomajor products like 3-hydroxy-nopinone but also tominor products like cis-pinic acid. A mechanismfor the formation of this dicarboxylic acid isproposed and its possible role in aerosol formationprocesses discussed.     

Predictability of large interannual Arctic sea-ice anomalies

In projections of twenty-first century climate, Arctic sea ice declines and at the same time exhibits strong interannual anomalies. Here, we investigate the potential to predict these strong sea-ice anomalies under a perfect-model assumption, using the Max-Planck-Institute Earth System Model in the same setup as in the Coupled Model Intercomparison Project Phase 5 (CMIP5). We study two cases of strong negative sea-ice anomalies: a 5-year-long anomaly for present-day conditions, and a 10-year-long anomaly for conditions projected for the middle of the twenty-first century. We treat these anomalies in the CMIP5 projections as the truth, and use exactly the same model configuration for predictions of this synthetic truth. We start ensemble predictions at different times during the anomalies, considering lagged-perfect and sea-ice-assimilated initial conditions. We find that the onset and amplitude of the interannual anomalies are not predictable. However, the further deepening of the anomaly can be predicted for typically 1 year lead time if predictions start after the onset but before the maximal amplitude of the anomaly. The magnitude of an extremely low summer sea-ice minimum is hard to predict: the skill of the prediction ensemble is not better than a damped-persistence forecast for lead times of more than a few months, and is not better than a climatology forecast for lead times of two or more years. Predictions of the present-day anomaly are more skillful than predictions of the mid-century anomaly. Predictions using sea-ice-assimilated initial conditions are competitive with those using lagged-perfect initial conditions for lead times of a year or less, but yield degraded skill for longer lead times. The results presented here suggest that there is limited prospect of predicting the large interannual sea-ice anomalies expected to occur throughout the twenty-first century.