Closure of a seaway: stratigraphic record and facies (Guadix basin, Southern Spain)

During the late Tortonian (upper Miocene), the Guadix Basin in S Spain formed one of the Betic corridors that connected the Mediterranean Sea with the Atlantic Ocean. The closure of this connection occurred in a series of steps, documented by three sedimentary units. A lower unit, consisting of basinal marls, shallow-water calcarenites and sands records the formation of a wide seaway. During deposition of the following unit this narrowed to a strait no more than 2 km in wide, triggering an intensification of currents that caused migration of submarine dunes preserved as giant cross-beds in bioclastic sands and conglomerates. Current flowed from the Mediterranean to the Atlantic. The third unit constitutes the youngest marine episode of the filling of the Guadix Basin. At this stage, the connection between the Mediterranean Sea and the Atlantic Ocean was broken, and a system of coastal coral reefs was established in the northern part of the Basin.     

Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions

Precise and efficient numerical simulation of transport processes in subsurface systems is a prerequisite for many site investigation or remediation studies. Random walk particle tracking (RWPT) methods have been introduced in the past to overcome numerical difficulties when simulating propagation processes in porous media such as advection-dominated mass transport. Crucial for the precision of RWPT methods is the accuracy of the numerically calculated ground water velocity field. In this paper, a global node-based method for velocity calculation is used, which was originally proposed by Yeh (Water Resour Res 7:1216–1225, 1981). This method is improved in three ways: (1) extension to unstructured grids, (2) significant enhancement of computational efficiency, and (3) extension to saturated (groundwater) as well as unsaturated systems (soil water). The novel RWPT method is tested with numerical benchmark examples from the literature and used in two field scale applications of contaminant transport in saturated and unsaturated ground water. To evaluate advective transport of the model, the accuracy of the velocity field is demonstrated by comparing several published results of particle pathlines or streamlines. Given the chosen test problem, the global node-based velocity estimation is found to be as accurate as the CK method (Cordes and Kinzelbach in Water Resour Res 28(11):2903–2911, 1992) but less accurate than the mixed or mixed-hybrid finite element methods for flow in highly heterogeneous media. To evaluate advective–diffusive transport, a transport problem studied by Hassan and Mohamed (J Hydrol 275(3–4):242–260, 2003) is investigated here and evaluated using different numbers of particles. The results indicate that the number of particles required for the given problem is decreased using the proposed method by about two orders of magnitude without losing accuracy of the concentration contours as compared to the published numbers.     

Process-Driven Direction-Dependent Asymmetry: Identification and Quantification of Directional Dependence in Spatial Fields

Natural processes generate spatial fields which reflect their specific properties. In this paper the effect of the direction of processes on the resulting spatial fields is investigated. This is done by extending the concept of reversibility used for time series to space. A novel copula based measure of asymmetry is defined which is an indicator of directional dependence. Contrary to traditional geostatistics where all points separated by a vector are considered irrespective of its sign, in this study the direction of the vector is also taken into consideration, leading to differences in the dependence corresponding to the vector \({\mathbf {h}}\) and \({-}{\mathbf {h}}\). The concept of directional dependence and the corresponding measure of asymmetry are defined using spatial copulas, and are thus independent of the scale of measurement. The result is a bivariate directional third-order moment based measure which can identify the direction in which the processes generating the spatial field acted. A statistical test to find the statistical significance of the asymmetry indicating directional dependence is presented. The methodology is tested on a number of synthetic and observed cases. Precipitation and groundwater quality parameters obtained using numerical models are first investigated. Regular dense grids obtained by numerical simulations show good correspondence between properties of the modeled processes and the new measure introduced. Measured variables observed on sparse irregular networks show similar behavior to the theoretical examples. Mean flow directions in groundwater and advection directions of precipitation fields can be detected from single snapshots. As a further example, dominant wind directions in the Sahara are found by investigating the digital terrain model.     

Architecture of tunnel valleys in the southeastern North Sea: new insights from high-resolution seismic imaging

Tunnel valleys are assumed to form near the margin of ice sheets. Hence, they can be used to reconstruct the dynamics of former ice margins. The detailed formation and infill of tunnel valleys, however, are still not well understood. Here, we present a dense grid of high-resolution 2D multi-channel reflection seismic data from the German sector of the southeastern North Sea imaging tunnel valleys in very great detail. Three tunnel valley systems were traced over distances ranging between 11 and 21 km. All tunnel valleys are completely filled and buried. They differ in incision depth, incision width and number of incisions. The tunnel valleys cut 130–380 m deep into Neogene, Palaeogene and Cretaceous sediments; they show a lower V-shaped and an upper U-shaped morphology. For individual tunnel valleys, the overall incision direction ranges from east–west to northeast–southwest. Two tunnel valleys intersect at an oblique angle without reuse of the thalweg. These valleys incise into a pre-existing glaciotectonic complex consisting of thrust sheets in the northwest of the study area. The analysis of the glaciotectonic complex and the tunnel valleys leads us to assume that we identified several marginal positions of (pre-)Elsterian ice lobes in the southeastern North Sea.     

Structural evolution of the Piqiang Fault Zone,NW Tarim Basin,China

The Piqiang Fault is a prominent strike-slip (tear) fault that laterally partitions the Keping Shan Thrust Belt in the NW Tarim Basin, China. In satellite images, the Piqiang Fault appears as a sharp, NW-trending lineament that can be traced for more than 70 km. It is oblique to the general structural trend of the thrust belt and subparallel to the thrust transport direction. This paper presents a structural analysis of the Piqiang Fault, based on satellite image interpretation and field data. A net loss of Late Paleozoic sediment across the fault zone implies that it was initiated as a major normal fault during the Early Permian, and corresponds to widespread extension and magmatism during this period. Differential erosion across the fault resulted in the subsequent removal of sediment from the east relative to the west. During the Middle to Late Cenozoic, contraction of the NW Tarim Basin and the formation of the Keping Shan Thrust Belt resulted in reactivation of the Piqiang Fault as a strike-slip (tear) fault. The fault has accommodated lateral differences in thrust density and spacing which have arisen due to the abrupt, pre-existing change in stratigraphic thickness across it. The Piqiang Fault provides an insight into the formation of oblique, strike-slip (tear) faults in contractional belts and demonstrates the importance of inherited basement structures in such settings.     

Distribution of intact and core GDGTs in marine sediments

We conducted a survey of archaeal GDGT (glycerol dibiphytanyl glycerol tetraether) distributions in marine sediments deposited in a range of depositional settings. The focus was comparison of two pools presumed to have distinct geobiological significance, i.e. intact polar GDGTs (IP GDGTs) and core GDGTs (C GDGTs). The former pool has been suggested to be related to living communities of benthic archaea in marine sediments, while the latter is commonly interpreted to consist of molecular fossils from past planktonic archaeal communities that inhabited the surface ocean. Understanding the link between these two pools is important for assessment of the validity of current molecular proxies for sedimentary archaeal biomass and past sea surface temperatures. The relative distributions of GDGTs in the two pools in a core at a CH₄ rich site in the Black Sea provide evidence for in situ production of glycosidic IP GDGTs and their subsequent degradation to corresponding C GDGTs on timescales that are short in geological terms. In addition, we monitored the relationship between the IP GDGT and C GDGT pools in a sample set from various ocean basins with subseafloor depth from a few cm to 320 m and 0 to 4 Myr in age. Notable differences between the two pools can be summarized as follows: the GDGT with acyclic biphytanes, GDGT-0, and its analogues with two and three cyclopentane moieties (GDGT-2 and -3) are generally more abundant in the pool of IP GDGTs, while crenarchaeol tends to be more abundant in the C GDGT pool. Consequently, the ring index is generally higher for the C GDGTs while TEX₈₆, a molecular proxy ratio not considering the two major GDGTs, tends to be higher in the IP GDGT pool. These differences in the proportion of individual GDGTs in the two pools are probably due to in situ production of IP GDGTs with distributions differing from those of C GDGTs. Despite these differences, we observed significant correlation of these two ratios between the two pools. Specifically, in both pools TEX₈₆ is high in sediments from warm oceanic regimes and low in cold regimes. We discuss these relationships and suggest that recycling of core GDGTs by benthic archaea is an important mechanism linking both molecular pools.     

Statistical analysis of the 2012–2013 Torreperogil–Sabiote seismic series,Spain

The study of earthquake swarms and their characteristics can improve our understanding of the transient processes that provoke seismic crises. The spatio-temporal process of the energy release is often linked with changes of statistical properties, and thus, seismicity parameters can help to reveal the underlying mechanism in time and space domains. Here, we study the Torreperogil–Sabiote 2012–2013 seismic series (southern Spain), which was relatively long lasting, and it was composed by more than 2000 events. The largest event was a magnitude 3.9 event which occurred on February 5, 2013. It caused slight damages, but it cannot explain the occurrence of the whole seismic crises which was not a typical mainshock–aftershock sequence. To shed some light on this swarm occurrence, we analyze the change of statistical properties during the evolution of the sequence, in particular, related to the magnitude and interevent time distributions. Furthermore, we fit a modified version of the epidemic type aftershock sequence (ETAS) model in order to investigate changes of the background rates and the trigger potential. Our results indicate that the sequence was driven by an aseismic transient stressing rate and that the system passes after the swarm occurrence to a new forcing regime with more typical tectonic characteristics.     

Rapid directivity detection by azimuthal amplitude spectra inversion

An early detection of the presence of rupture directivity plays a major role in the correct estimation of ground motions and risks associated to the earthquake occurrence. We present here a simple method for a fast detection of rupture directivity, which may be additionally used to discriminate fault and auxiliary planes and have first estimations of important kinematic source parameters, such as rupture length and rupture time. Our method is based on the inversion of amplitude spectra from P-wave seismograms to derive the apparent duration at each station and on the successive modelling of its azimuthal behaviour. Synthetic waveforms are built assuming a spatial point source approximation, and the finite apparent duration of the spatial point source is interpreted in terms of rupture directivity. Since synthetic seismograms for a point source are calculated very quickly, the presence of directivity may be detected within few seconds, once a focal mechanism has been derived. The method is here first tested using synthetic datasets, both for linear and planar sources, and then successfully applied to recent Mw 6.2–6.8 shallow earthquakes in Peloponnese, Greece. The method is suitable for automated application and may be used to improve kinematic waveform modelling approaches.