Erosion of continental margins in the Western Mediterranean due to sea-level stagnancy during the Messinian Salinity Crisis

High-resolution multi-channel seismic data from continental slopes with minor sediment input off southwest Mallorca Island, the Bay of Oran (Algeria) and the Alboran Ridge reveal evidence that the Messinian erosional surface is terraced at an almost constant depth interval between 320 and 380 m below present-day sea level. It is proposed that these several hundred- to 2,000-m-wide terraces were eroded contemporaneously and essentially at the same depth. Present-day differences in these depths result from subsidence or uplift in the individual realms. The terraces are thought to have evolved during one or multiple periods of sea-level stagnancy in the Western Mediterranean Basin. According to several published scenarios, a single or multiple periods of relative sea-level stillstand occurred during the Messinian desiccation event, generally known as the Messinian Salinity Crisis. Some authors suggest that the stagnancy started during the refilling phase of the Mediterranean basins. When the rising sea level reached the height of the Sicily Sill, the water spilled over this swell into the eastern basin. The stagnancy persisted until sea level in the eastern basin caught up with the western Mediterranean water level. Other authors assigned periods of sea-level stagnancy to drawdown phases, when inflowing waters from the Atlantic kept the western sea level constant at the depth of the Sicily Sill. Our findings corroborate all those Messinian sea-level reconstructions, forwarding that a single or multiple sea-level stagnancies at the depth of the Sicily Sill lasted long enough to significantly erode the upper slope. Our data also have implications for the ongoing debate of the palaeo-depth of the Sicily Sill. Since the Mallorcan plateau experienced the least vertical movement, the observed terrace depth of 380 m there is inferred to be close to the Messinian depth of this swell.     

Numerical modelling of regional faults in land subsidence prediction above gas/oil reservoirs

The stress variation induced by gas/oil production may activate pre‐existing regional faults. This may enhance the expected land subsidence due to the generation of mechanically weak points close to the producing field. A class of elasto‐plastic interface elements (IE), specifically designed to address the mechanical behaviour of faults over a regional scale, is integrated into a finite element (FE) geomechanical model and used to investigate the role exerted by active faults in anthropogenic land subsidence. The importance of regional faults depends on a variety of factors including depth of the depleted reservoir, fault number, orientation and size, geomechanical properties of porous medium, pore pressure drawdown induced by fluid production, etc. With the aid of some representative examples, a useful indication is provided as to where and how fault activation may influence both magnitude and extent of the land subsidence bowl above producing gas/oil reservoirs, pointing to a generally limited impact on the ground surface. The simulation of a real faulted gas reservoir in a complex 3‐D setting shows that the proposed IE can be simply and efficiently incorporated into a FE geomechanical model, thus improving the quality of the stress and displacement prediction. Copyright © 2007 John Wiley & Sons, Ltd.     

The geodynamic and limnological evolution of Balkan Lake Ohrid,possibly the oldest extant lake in Europe

Studies of the upper 447 m of the DEEP site sediment succession from central Lake Ohrid, Balkan Peninsula, North Macedonia and Albania provided important insights into the regional climate history and evolutionary dynamics since permanent lacustrine conditions established at 1.36 million years ago (Ma). This paper focuses on the entire 584-m-long DEEP sediment succession and a comparison to a 197-m-long sediment succession from the Pestani site ~5 km to the east in the lake, where drilling ended close to the bedrock, to unravel the earliest history of Lake Ohrid and its basin development. ²⁶Al/¹⁰Be dating of clasts from the base of the DEEP sediment succession implies that the sedimentation in the modern basin started at c. 2 Ma. Geophysical, sedimentological and micropalaeontological data allow for chronological information to be transposed from the DEEP to the Pestani succession. Fluvial conditions, slack water conditions, peat formation and/or complete desiccation prevailed at the DEEP and Pestani sites until 1.36 and 1.21 Ma, respectively, before a larger lake extended over both sites. Activation of karst aquifers to the east probably by tectonic activity and a potential existence of neighbouring Lake Prespa supported filling of Lake Ohrid. The lake deepened gradually, with a relatively constant vertical displacement rate of ~0.2 mm a⁻¹ between the central and the eastern lateral basin and with greater water depth presumably during interglacial periods. Although the dynamic environment characterized by local processes and the fragmentary chronology of the basal sediment successions from both sites hamper palaeoclimatic significance prior to the existence of a larger lake, the new data provide an unprecedented and detailed picture of the geodynamic evolution of the basin and lake that is Europe’s presumed oldest extant freshwater lake.     

Parallel solution to ill‐conditioned FE geomechanical problems

Parallel computers are potentially very attractive for the implementation of large size geomechanical models. One of the main difficulties of parallelization, however, relies on the efficient solution of the frequently ill‐conditioned algebraic system arising from the linearization of the discretized equilibrium equations. While very efficient preconditioners have been developed for sequential computers, not much work has been devoted to parallel solution algorithms in geomechanics. The present study investigates the state‐of‐the‐art performance of the factorized sparse approximate inverse (FSAI) as a preconditioner for the iterative solution of ill‐conditioned geomechanical problems. Pre‐and post‐filtration strategies are experimented with to increase the FSAI efficiency. Numerical results show that FSAI exhibits a promising potential for parallel geomechanical models mainly because of its almost ideal scalability. With the present formulation, however, at least 4 or 8 processors are required in the selected test cases to outperform one of the most efficient sequential algorithms available for FE geomechanics, i.e. the multilevel incomplete factorization (MIF). Further research is needed to improve the FSAI efficiency with a more effective selection of the preconditioner non‐zero pattern. Copyright © 2011 John Wiley & Sons, Ltd.     

Palaeointensity results from Ethiopian basalts: implications for the Oligocene geomagnetic field strength

From a large collection of Ethiopian flood basalts (~30  Myr old) sampled for magnetostratigraphy, ⁴⁰Ar/³⁹Ar geochronology and geochemical investigations, 47 samples were selected in order to test their suitability for Thellier palaeointensity experiments. Only 17 samples from eight individual flows yielded reliable palaeointensity estimates, with flow-mean virtual dipole moments ranging from 3.0 to 10.5 × 10²²  A  m² .
  A critical review of the Oligocene palaeointensity data set, including these new Ethiopian data, indicates an Oligocene mean virtual dipole moment of 5.1 ± 2.5 × 10²²  A  m² for the complete data set. After applying mild selection criteria, the reduced data set yields a mean value of only 4.6 ± 1.9 × 10²²  A  m² . This value is significantly lower than the present-day field strength but is higher than the Mesozoic dipole low mean field. This low Oligocene field might be in agreement with the high palaeosecular variation and rather high non-dipole field invoked around 30  Ma. However, the Oligocene data set is largely dependent on the palaeointensity determinations from Armenia, obtained mainly from baked contacts, which show amazingly low dispersion at both flow and between-flow levels. More data are needed to reduce the weight of these determinations on the mean value and avoid a possible bias.