Crustal motion and deformation in Greece from a decade of GPS measurements, 1993–2003

The Hellenic plate boundary region, located in the collision zone between the Nubian/Arabian and Eurasian lithospheric plates, is one of the seismo-tectonically most active areas of Europe. During the last 15 years, GPS measurements have been used to determine the crustal motion in the area of Greece with the aim to better understand the geodynamical processes of this region. An extended reoccupation network covering whole Greece has been measured periodically in numerous GPS campaigns since the late eighties, and a continuous GPS network has been operated in the region of the Ionian Sea since 1995. In this paper, we present a new detailed high-quality solution of continuous and campaign-type measurements acquired between 1993 and 2003. During the GPS processing, a special effort was made to obtain consistent results with highest possible accuracies and reliabilities. Data of 54 mainly European IGS and EUREF sites were included in the GPS processing in order to obtain results which are internally consistent with the European kinematic field and order to allow for a regional interpretation. After an overview of the results of the IGS/EUREF sites, the results from more than 80 stations in Greece are presented in terms of velocities, time series, trajectories and strain rates. Previous geodetic, geological and seismological findings are generally confirmed and substantially refined. New important results include the observation of deformation zones to the north and to the south of the North Aegean Trough and in the West Hellenic arc region, arc-parallel extension of about 19 mm/yr along the Hellenic arc, and compression between the Ionian islands and the Greek mainland. Due to continuous long-term observations of 4–8 years, it was possible to extract height changes from the GPS time series. In Greece, we observe a differential subsidence of the order of 2 mm/yr between the northern and central Ionian islands across the Kefalonia fault zone. The differential subsidence of the central Ionian islands with respect to the northwestern Greek mainland amounts to 4 mm/yr.     

Quaternary fault segmentation and interaction in the epicentral area of the 1561 earthquake (Mw = 6.4), Vallo di Diano, southern Apennines, Italy

The main structural characteristics of the Caggiano and Polla faults, exposed in the epicentral area of the 1561 earthquake (Mw = 6.4), southern Italy, have been investigated in detail to assess their spatial and temporal properties, and to evaluate their seismogenic potential. These right stepping normal faults show an overlap of about 7 km and an across strike separation of about 4 km. The geometric relationships between the Caggiano and Polla faults, but also the displacement distribution along each fault, demonstrate that they have been strongly interacting throughout the Pleistocene. Nevertheless, geological evidence of Holocene tectonic activity was mainly recognized along the Caggiano Fault (faulted late glacial deposits) and in the southernmost part of the Polla Fault (faulted deposits of probably Late Pleistocene age). This suggests that the Caggiano Fault can be considered as the most tectonically active fault in the Vallo di Diano Fault System. By calculating Coulomb stress changes, we have constrained modes of mechanical interactions between the two faults in a scenario compatible with the 1561 earthquake. This approach allows us to argue that both the Caggiano and the Polla Faults are probably linked at depth, and part of the same seismogenic structure which may be potentially responsible for composite ruptures with magnitude ≥ 6.5.     

Holocene activity of the Scilla Fault, Southern Calabria: Insights from coastal morphological and structural investigations

Integration of on-land and offshore geomorphological and structural investigations coupled to extensive radiometric dating of co-seismically uplifted Holocene beaches allows characterization of the geometry, kinematics and seismotectonics of the Scilla Fault, which borders the eastern side of the Messina Strait in Calabria, Southern Italy. This region has been struck by destructive historical earthquakes, but knowledge of geologically-based source parameters for active faults is relatively poor, particularly for those running mostly offshore, as the Scilla Fault does. The  30 km-long normal fault may be divided into three segments of  10 km individual length, with the central and southern segments split in at least two strands. The central and northern segments are submerged, and in this area marine geophysical data indicate a youthful morphology and locally evidence for active faulting. The on-land strand of the western segment displaces marine terraces of the last interglacial (124 to 83 ka), but seismic reflection profiles suggest a full Quaternary activity. Structural data collected on bedrock faults exposed along the on-land segment provide evidence for normal slip and  NW-SE extension, which is consistent with focal mechanisms of large earthquakes and GPS velocity fields in the region. Detailed mapping of raised Holocene marine deposits exposed at the coastline straddling of the northern and central segments supplies evidence for two co-seismic displacements at  1.9 and  3.5 ka, and a possible previous event at  5 ka. Co-seismic displacements show a consistent site value and pattern of along-strike variation, suggestive of characteristic-type behaviour for the fault. The  1.5–2.0 m average co-seismic slips during these events document M_e  6.9–7.0 earthquakes with  1.6–1.7 ka recurrence time. Because hanging-wall subsidence cannot be included into slip magnitude computation, these slips reflect footwall uplift, and represent minimum average estimates. The palaeoseismological record based on the palaeo-shorelines suggests that the last rupture on the Scilla Fault during the February 6, 1783 M_w = 5.9–6.3 earthquake was at the expected time but it may have not entirely released the loaded stress since the last great event at  1.9 ka. Comparison of the estimated co-seismic extension rate based on the Holocene shoreline record with available GPS velocities indicates that the Scilla Fault accounts for at least  15–20% of the contemporary geodetic extension across the Messina Strait.     

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.     

The Lower Permian Wasp Head Formation,Sydney Basin: high‐latitude,shallow marine sedimentation following the late Asselian to early Sakmarian glacial event in eastern Australia

The Lower Permian Wasp Head Formation (early to middle Sakmarian) is a ～95 m thick unit that was deposited during the transition to a non‐glacial period following the late Asselian to early Sakmarian glacial event in eastern Australia. This shallow marine, sandstone‐dominated unit can be subdivided into six facies associations. (i) The marine sediment gravity flow facies association consists of breccias and conglomerates deposited in upper shoreface water depths. (ii) Upper shoreface deposits consist of cross‐stratified, conglomeratic sandstones with an impoverished expression of the Skolithos Ichnofacies. (iii) Middle shoreface deposits consist of hummocky cross‐stratified sandstones with a trace fossil assemblage that represents the Skolithos Ichnofacies. (iv) Lower shoreface deposits are similar to middle shoreface deposits, but contain more pervasive bioturbation and a distal expression of the Skolithos Ichnofacies to a proximal expression of the Cruziana Ichnofacies. (v) Delta‐influenced, lower shoreface‐offshore transition deposits are distinguished by sparsely bioturbated carbonaceous mudstone drapes within a variety of shoreface and offshore deposits. Trace fossil assemblages represent distal expressions of the Skolithos Ichnofacies to stressed, proximal expressions of the Cruziana Ichnofacies. Impoverished trace fossil assemblages record variable and episodic environmental stresses possibly caused by fluctuations in sedimentation rates, substrate consistencies, salinity, oxygen levels, turbidity and other physio‐chemical stresses characteristic of deltaic conditions. (vi) The offshore transition‐offshore facies association consists of mudstone and admixed sandstone and mudstone with pervasive bioturbation and an archetypal to distal expression of the Cruziana Ichnofacies. The lowermost ～50 m of the formation consists of a single deepening upward cycle formed as the basin transitioned from glacioisostatic rebound following the Asselian to early Sakmarian glacial to a regime dominated by regional extensional subsidence without significant glacial influence. The upper ～45 m of the formation can be subdivided into three shallowing upward cycles (parasequences) that formed in the aftermath of rapid, possibly glacioeustatic, rises in relative sea‐level or due to autocyclic progradation patterns. The shift to a parasequence‐dominated architecture and progressive decrease in ice‐rafted debris upwards through the succession records the release from glacioisostatic rebound and amelioration of climate that accompanied the transition to broadly non‐glacial conditions.     

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.     

Glauconite from Paleogene volcano-terrigenous rocks in Western Kamchatka

It is shown that glauconite-bearing interbeds are widespread in the layer-by-layer studied sections on the Sea of Okhotsk coast (Mainach section) and Kheisliveem River valley (Kavran section), the volcanoterrigenous rocks of the Kovachin, Amanin, and Gakkhin formations of the Paleogene in western Kamchatka (Upper Eocene-Lower Oligocene boundary beds). Detailed mineralogical and structural-crystallochemical characteristics of glauconite from the Amanin Formation are presented. It is suggested that such glauconite should not be used for geochronological purposes.Some specific features of glauconite formation, particularly, the preservation of specific morphological forms at high accumulation rates of volcano-terrigenous rocks, are discussed. Possibility of the formation of glauconite with the active influence of bacterial metabolism is considered.     

A Pliocene mega-tsunami deposit and associated features in the Ranquil Formation,southern Chile

An exceptionally large tsunami affected the coastline of southern Chile during the Pliocene. Its backflow eroded coarse beach and coastal dune sediments and redistributed them over the continental shelf and slope. Sandstone dykes and sills injected from the base of the resulting hyperconcentrated flow into underlying cohesive muds, assisted in plucking up large blocks of the latter and incorporating them into the flow. Locally, the rip-up intraclasts were fragmented further by smaller-scale injections to form a distinct breccia of angular to rounded mudstone clasts within a medium to coarse sandstone matrix. Sandstone sills in places mimic normal sedimentary beds, complete with structures resembling inverse gradation, planar laminae, as well as ripple and trough cross-lamination. These were probably formed by internal sediment flow and shear stress as the semi-liquefied sand was forcefully injected into cracks. In borehole cores, such sills can easily be misinterpreted as normal sedimentary beds, which can have important implications for hydrocarbon exploration.     

Anhydrite pseudomorphs and the origin of stratiform Cu–Co ores in the Katangan Copperbelt (Democratic Republic of Congo)

The stratiform Cu–Co ore mineralisation in the Katangan Copperbelt consists of dispersed sulphides and sulphides in nodules and lenses, which are often pseudomorphs after evaporites. Two types of pseudomorphs can be distinguished in the nodules and lenses. In type 1 examples, dolomite precipitated first and was subsequently replaced by Cu–Co sulphides and authigenic quartz, whereas in type 2 examples, authigenic quartz and Cu–Co sulphides precipitated prior to dolomite and are coarse-grained. The sulphur isotopic composition of the copper–cobalt sulphides in the type 1 pseudomorphs is between −10.3 and 3.1‰ relative to the Vienna Canyon Diablo Troilite, indicating that the sulphide component was derived from bacterial sulphate reduction (BSR). The generation of during this process caused the precipitation and replacement of anhydrite by dolomite. A second product of BSR is the generation of H₂S, resulting in the precipitation of Cu–Co sulphides from the mineralising fluids. Initial sulphide precipitation occurred along the rim of the pseudomorphs and continued towards the core. Precipitation of authigenic quartz was most likely induced by a pH decrease during sulphide precipitation. Fluid inclusion data from quartz indicate the presence of a high-salinity (8–18 eq. wt.% NaCl) fluid, possibly derived from evaporated seawater which migrated through the deep subsurface. ⁸⁷Sr/⁸⁶Sr ratios of dolomite in type 1 nodules range between 0.71012 and 0.73576, significantly more radiogenic than the strontium isotopic composition of Neoproterozoic marine carbonates (⁸⁷Sr/⁸⁶Sr = 0.7056–0.7087). This suggests intense interaction with siliciclastic sedimentary rocks and/or the granitic basement. The low carbon isotopic composition of the dolomite in the pseudomorphs (−7.02 and −9.93‰ relative to the Vienna Pee Dee Belemnite, V-PDB) compared to the host rock dolomite (−4.90 and +1.31‰ V-PDB) resulted from the oxidation of organic matter during BSR.     

Development of the mixed conifer forest in northern New Mexico and its relationship to Holocene environmental change

Chihuahueños Bog (2925 m) in the Jemez Mountains of northern New Mexico contains one of the few records of late-glacial and postglacial development of the mixed conifer forest in southwestern North America. The Chihuahueños Bog record extends to over 15,000 cal yr BP. An Artemisia steppe, then an open Picea woodland grew around a small pond until ca. 11,700 cal yr BP when Pinus ponderosa became established. C/N ratios, δ¹³C and δ¹⁵N values indicate both terrestrial and aquatic organic matter was incorporated into the sediment. Higher percentages of aquatic algae and elevated C/N ratios indicate higher lake levels at the opening of the Holocene, but a wetland developed subsequently as climate warmed. From ca. 8500 to 6400 cal yr BP the pond desiccated in what must have been the driest period of the Holocene there. C/N ratios declined to their lowest Holocene levels, indicating intense decomposition in the sediment. Wetter conditions returned after 6400 cal yr BP, with conversion of the site to a sedge bog as groundwater levels rose. Higher charcoal influx rates after 6400 cal yr BP probably result from greater biomass production rates. Only minor shifts in the overstory species occurred during the Holocene, suggesting that mixed conifer forest dominated throughout the record.