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.     

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.     

Discriminating between pore‐filling load and bed‐structure load: a new porosity‐based method,exemplified for the river Rhine

Sediments contained in the river bed do not necessarily contribute to morphological change. The finest part of the sediment mixture often fills the pores between the larger grains and can be removed without causing a drop in bed level. The discrimination between pore‐filling load and bed‐structure load, therefore, is of practical importance for morphological predictions. In this study, a new method is proposed to estimate the cut‐off grain size that forms the boundary between pore‐filling load and bed‐structure load. The method evaluates the pore structure of the river bed geometrically. Only detailed grain‐size distributions of the river bed are required as input to the method. A preliminary validation shows that the calculated porosity and cut‐off size values agree well with experimental data. Application of the new cut‐off size method to the river Rhine demonstrates that the estimated cut‐off size decreases in a downstream direction from about 2 to 0·05 mm, covariant with the downstream fining of bed sediments. Grain size fractions that are pore‐filling load in the upstream part of the river thus gradually become bed‐structure load in the downstream part. The estimated (mass) percentage of pore‐filling load in the river bed ranges from 0% in areas with a unimodal river bed, to about 22% in reaches with a bimodal sand‐gravel bed. The estimated bed porosity varies between 0·15 and 0·35, which is considerably less than the often‐used standard value of 0·40. The predicted cut‐off size between pore‐filling load and bed‐structure load (D_c,p) is fundamentally different from the cut‐off size between wash‐load and bed‐material load (D_c,w), irrespective of the method used to determine D_c,p or D_c,w. D_c,w values are in the order of 10^?1 mm and mainly dependent on the flow characteristics, whereas D_c,p values are generally much larger (about 10⁰ mm in gravel‐bed rivers) and dependent on the bed composition. Knowledge of D_c,w is important for the prediction of the total sediment transport in a river (including suspended fines that do not interact with the bed), whereas knowledge of D_c,p helps to improve morphological predictions, especially if spatial variations in D_c,p are taken into account. An alternative to using a spatially variable value of D_c,p in morphological models is to use a spatially variable bed porosity, which can also be predicted with the new method. In addition to the morphological benefits, the new method also has sedimentological applications. The possibility to determine quickly whether a sediment mixture is clast‐supported or matrix‐supported may help to better understand downstream fining trends, sediment entrainment thresholds and variations in hydraulic conductivity.     

Trace element distribution among rock-forming minerals from metamorphosed to partially molten basic igneous rocks in a contact aureole (Fuerteventura, Canaries)

Low pressure partial melting of basanitic and ankaramitic dykes gave rise to unusual, zebra-like migmatites, in the contact aureole of a layered pyroxenite–gabbro intrusion, in the root zone of an ocean island (Basal Complex, Fuerteventura, Canary Islands). These migmatites are characterised by a dense network of closely spaced, millimetre-wide leucocratic segregations. Their mineralogy consists of plagioclase (An_32–36), diopside, biotite, oxides (magnetite, ilmenite), +/− amphibole, dominated by plagioclase in the leucosome and diopside in the melanosome. The melanosome is almost completely recrystallised, with the preservation of large, relict igneous diopside phenocrysts in dyke centres. Comparison of whole-rock and mineral major- and trace-element data allowed us to assess the redistribution of elements between different mineral phases and generations during contact metamorphism and partial melting.

Dykes within and outside the thermal aureole behaved like closed chemical systems. Nevertheless, Zr, Hf, Y and REEs were internally redistributed, as deduced by comparing the trace element contents of the various diopside generations. Neocrystallised diopside – in the melanosome, leucosome and as epitaxial phenocryst rims – from the migmatite zone, are all enriched in Zr, Hf, Y and REEs compared to relict phenocrysts. This has been assigned to the liberation of trace elements on the breakdown of enriched primary minerals, kaersutite and sphene, on entering the thermal aureole. Major and trace element compositions of minerals in migmatite melanosomes and leucosomes are almost identical, pointing to a syn- or post-solidus reequilibration on the cooling of the migmatite terrain i.e. mineral–melt equilibria were reset to mineral–mineral equilibria.     

Partial crustal melting beneath the Betic Cordillera (SE Spain): The case study of Mar Menor volcanic suite

The Neogene Volcanic Province (NVP) within the Betic Cordillera (SE Spain) consists of three main metapelitic enclave suites (from SW to NE: El Hoyazo, Mazarrón and Mar Menor). Since the NVP represents a singular place in the world where crustal enclaves were immediately quenched after melting, their microstructures provide a “photograph” of the conditions at depth just after the moment of the melting.

The thermobarometric information provided by the different microstructural assemblages has been integrated with the geophysical and geodynamical published data into a model of the petrologic evolution of the Mar Menor enclaves. They were equilibrated at 2–3 kbar, 850–900 °C, and followed a sequence of heating melt producing reactions. A local cooling event evidenced by minor melt crystallization preceded the eruption.

The lower crustal studies presented in this work contribute to the knowledge of: (i) the partial melting event beneath the Mar Menor volcanic suite through a petrologic detailed study of the enclaves; (ii) how the microstructures of fast cooled anatectic rocks play an important role in tracing the magma evolution in a chamber up to the eruption, and how they can be used as pseudothermobarometers; (iii) the past and current evolution of the Alborán Domain (Betic Cordillera) and Mediterranean Sea, and how the base of a metapelitic crust has melted within an active geodynamic setting.     

Cenozoic Eucla Basin and associated palaeovalleys,southern Australia — Climatic and tectonic influences on landscape evolution,sedimentation and heavy mineral accumulation

The Eucla Basin including the vast Nullarbor Plain lies on the margins of the Yilgarn, Musgrave and Gawler cratons in southern Australia and owes its distinctive landscape to a unique set of interactions between eustatic, climatic and tectonic processes over the last ~ 50 Ma. Understanding of the history of the basin and the palaeovalleys that drained from the surrounding cratons are important because they contain major mineral deposits, and the sediments derived from them contain remobilised gold, uranium, and heavy minerals. In particular, a remarkably preserved palaeoshoreline sequence along the north-eastern margin of the Eucla Basin is highly prospective for heavy mineral placer deposits. The record of marine, marginal marine, estuarine, fluvial and lacustrine environments, as constrained mainly by an extensive borehole dataset, reflects major depositional events during the Palaeocene–Early Eocene, Middle–Late Eocene, Oligocene–Early Miocene, Middle Miocene–Early Pliocene and Pliocene–Quaternary. These events reflect the key role of eustatic sea-level variation which, during highstands, inundated the craton margins, flooding palaeovalleys to up to 400 km inboard of the present coastline. However, a systematic eastward migration of the depocentre across the Eucla Basin during the Neogene, together with apparent flow reversals in a number of palaeovalley systems draining the Gawler Craton, suggest that the Eucla Basin has also been subject to differential vertical movements, expressed as a west-side up, east-side down tilting of ~ 100–200 m. This differential movement forms part of a broader north-down–southwest-up dynamic topographic tilting of the Australian continent associated with relatively fast (6–7 cm/yr) northward plate motion since fast spreading commenced in the Southern Ocean at ~ 43 Ma. We suggest that the evolving dynamic topography field has played a key role in facilitating development of placer deposits, largely through multistage, eastward reworking of near-shore sequences during highstand transgressive cycles on a progressively tilting platform under the influence of persistent westerly weather systems.     

The Anarak, Jandaq and Posht-e-Badam metamorphic complexes in central Iran: New geological data, relationships and tectonic implications

The Anarak, Jandaq and Posht-e-Badam metamorphic complexes occupy the NW part of the Central-East Iranian Microcontinent and are juxtaposed with the Great Kavir block and Sanandaj-Sirjan zone. Our recent findings redefine the origin of these complexes, so far attributed to the Precambrian–Early Paleozoic orogenic episodes, and now directly related to the tectonic evolution of the Paleo-Tethys Ocean. This tectonic evolution was initiated by Late Ordovician–Early Devonian rifting events and terminated in the Triassic by the Eocimmerian collision event due to the docking of the Cimmerian blocks with the Asiatic Turan block.

The “Variscan accretionary complex” is a new name we proposed for the most widely distributed metamorphic rocks connected to the Anarak and Jandaq complexes. This accretionary complex exposed from SW of Jandaq to the Anarak and Kabudan areas is a thick and fine grain siliciclastic sequence accompanied by marginal-sea ophiolitic remnants, including gabbro-basalts with a supra-subduction-geochemical signature. New ⁴⁰Ar/³⁹Ar ages are obtained as 333–320 Ma for the metamorphism of this sequence under greenschist to amphibolite facies. Moreover, the limy intercalations in the volcano-sedimentary part of this complex in Godar-e-Siah yielded Upper Devonian–Tournaisian conodonts. The northeastern part of this complex in the Jandaq area was intruded by 215 ± 15 Ma arc to collisional granite and pegmatites dated by ID-TIMS and its metamorphic rocks are characterized by some ⁴⁰Ar/³⁹Ar radiometric ages of 163–156 Ma.

The “Variscan” accretionary complex was northwardly accreted to the Airekan granitic terrane dated at 549 ± 15 Ma. Later, from the Late Carboniferous to Triassic, huge amounts of oceanic material were accreted to its southern side and penetrated by several seamounts such as the Anarak and Kabudan. This new period of accretion is supported by the 280–230 Ma ⁴⁰Ar/³⁹Ar ages for the Anarak mild high-pressure metamorphic rocks and a 262 Ma U–Pb age for the trondhjemite–rhyolite association of that area. The Triassic Bayazeh flysch filled the foreland basin during the final closure of the Paleo-Tethys Ocean and was partly deposited and/or thrusted onto the Cimmerian Yazd block.

The Paleo-Tethys magmatic arc products have been well-preserved in the Late Devonian–Carboniferous Godar-e-Siah intra-arc deposits and the Triassic Nakhlak fore-arc succession. On the passive margin of the Cimmerian block, in the Yazd region, the nearly continuous Upper Paleozoic platform-type deposition was totally interrupted during the Middle to Late Triassic. Local erosion, down to Lower Paleozoic levels, may be related to flexural bulge erosion. The platform was finally unconformably covered by Liassic continental molassic deposits of the Shemshak.

One of the extensional periods related to Neo-Tethyan back-arc rifting in Late Cretaceous time finally separated parts of the Eocimmerian collisional domain from the Eurasian Turan domain. The opening and closing of this new ocean, characterized by the Nain and Sabzevar ophiolitic mélanges, finally transported the Anarak–Jandaq composite terrane to Central Iran, accompanied by large scale rotation of the Central-East Iranian Microcontinent (CEIM). Due to many similarities between the Posht-e-Badam metamorphic complex and the Anarak–Jandaq composite terrane, the former could be part of the latter, if it was transported further south during Tertiary time.     

Petrography and textural development of inorganic and biogenic lithotypes in a relict barite tufa deposit at Flybye Springs, NT, Canada

A relict mound of Holocene barite (BaSO₄) tufa underlies the Flybye Springs, a small, barium‐rich, cold sulphur spring system in the Northwest Territories of Canada. The tufa is composed of relatively pure barite with ≤0·34 wt% Ca²⁺ and ≤0·77 wt% Sr²⁺. The mound is made up of coated bubble, raft, undulatory sheet, stromatolitic, coated grain and detrital conglomerate barite tufa. Although previously unreported in barite, these lithotypes are akin to facies found in many carbonate spring deposits. Raft and ooid‐coated grain tufa was formed via ‘inorganic’ barite precipitation in spring water ponds and tributaries where rapid oxidation of sulphide to sulphate established barite supersaturation. Undulatory sheet tufa may have formed by the reaction of dissolved barium with sulphate derived from the oxidation of extracellular polysaccharide‐rich colloidal sulphur films floating in oxygenated, barite‐saturated spring water ponds. Coated bubble, oncoid‐coated grain and stromatolitic tufa with filamentous microfossils was formed in close association with sulphur‐tolerant microbes inhabiting dysoxic and oxygenated spring water tributaries and ponds. Adsorption of dissolved barium to microbial extracellular polysaccharide probably facilitated the development of these ‘biogenic’ lithotypes. Detrital conglomerate tufa was formed by barite cementation of microdetrital tufa, allochthonous lithoclasts and organic detritus, including caribou hair. Biogenic textures, organic artefacts and microfossils in the Flybye barite tufa have survived diagenetic aggradational recrystallization and precipitation of secondary cements, indicating the potential for palaeoecological information to be preserved in barite in the geological record. Similarities between the Flybye barite tufa and carbonate spring deposits demonstrate that analogous textures can develop in chemical sedimentary systems with distinct mineralogy, biology and physiochemistry.     

Late Quaternary contourites and glaciomarine sedimentation in the Fram Strait

Two sites in the eastern Fram Strait, the Vestnesa Ridge and the Yermak Plateau, have been surveyed and sampled providing a depositional record over the last glacial‐interglacial cycle. The Fram Strait is the only deep‐water connection from the Arctic Ocean to the North Atlantic and contains a marine sediment record of both high latitude thermohaline flow and ice sheet interaction. On the Vestnesa Ridge, the western Svalbard margin, a sediment drift was identified in 1226 m of water. Gravity and multicores from the crest of the drift recovered turbidites and contourites. ¹⁴C dating indicates an age range of 8287 to 26 900 years BP (Early Holocene to Late Weichselian). The Yermak Plateau is characterized by slope sediments in 961 m of water. Gravity and multicores recovered contourites and hemipelagites. ¹⁴C ages were between 8615 and 46 437 years BP (Early Holocene to mid‐Weichselian). Downcore dinoflagellate cyst analyses from both sites provide a record of changing surface water conditions since the mid‐Weichselian, suggesting variable sea ice extent, productivity and polynyas present even during the Last Glacial Maximum. Four layers of ice‐rafted debris were also identified and correlated within the cores. These events occurred ca at 9, 24 to 25, 26 to 27 and 43 ka, asynchronous with Heinrich layers in the wider north‐east Atlantic and here interpreted as reflecting instability in the Svalbard/Barents Ice sheet and the northward advection of warm Atlantic water during the Late Weichselian. The activity of the ancestral West Spitsbergen Current is interpreted using mean sortable silt records from the cores. On the Vestnesa Ridge drift the modern mass accumulation rate, calculated using excess ²¹⁰Pb, is 0·076 g cm^?2 year^?1. On the Yermak Plateau slope the modern mass accumulation rate is 0·053 g cm^?2 year^?1.     

Fluid and source magma evolution of the Questa porphyry Mo deposit,New Mexico,USA

Combined fluid inclusion microthermometry and microanalysis by laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) are used to constrain the hydrothermal processes forming a typical Climax-type porphyry Mo deposit. Molybdenum mineralisation at Questa occurred in two superimposed hydrothermal stages, a magmatic-hydrothermal breccia and later stockwork veining. In both stages, texturally earliest fluids were single-phase, of low salinity (~7 wt.% NaCl_equiv.) and intermediate-density. Upon decompression to ~300 bar, they boiled off a vapour phase, leaving behind a residual brine (up to 45 wt.% NaCl_equiv) at temperatures of ~420°C. The highest average Mo concentrations in this hot brine were ~500 μg/g, exceeding the Mo content of the intermediate-density input fluid by about an order of magnitude and reflecting pre-concentration of Mo by fluid phase separation prior to MoS₂ deposition from the brine. Molybdenum concentrations in brine inclusions, then, decrease down to 5 μg/g, recording Mo precipitation in response to cooling of the saline liquid to ~360°C. Molybdenite precipitation from a dense, residual and probably sulphide-depleted brine is proposed to explain the tabular shape of the ore body and the absence of Cu-Fe sulphides in contrast to the more common Cu-Mo deposits related to porphyry stocks. Cesium and Rb concentrations in the single-phase fluids of the breccia range from 2 to 8 and from 40 to 65 μg/g, respectively. In the stockwork veins, Cs and Rb concentrations are significantly higher (45–90 and 110–230 μg/g, respectively). Because Cs and Rb are incompatible and hydrothermally non-reactive elements, the systematic increase in their concentration requires two distinct pulses of fluid exsolution from a progressively more fractionated magma. By contrast, major element and ore metal concentrations of these two fluid pulses remain essentially constant. Mass balance calculations using fluid chemical data from LA-ICPMS suggest that at least 25 km³ of melt and 7 Gt of deep input fluid were necessary to provide the amount of Mo contained in the stockwork vein stage alone. While the absolute amounts of fluid and melt are uncertain, the well-constrained element ratios in the fluids together with empirical fluid/melt partition coefficients derived from the inclusion analyses suggest a high water content of the source melt of ~10%. In line with other circumstantial evidence, these results suggest that initial fluid exsolution may have occurred at a confining pressure exceeding 5 kbar. The source of the molybdenum-mineralising fluids probably was a particularly large magma chamber that crystallised and fractionated in the lower crust or at mid-crustal level, well below the shallow intrusions immediately underlying Questa and other porphyry molybdenum deposits. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.