Finite-fault slip models for the 15 April 1979 (Mw 7.1) Montenegro earthquake and its strongest aftershock of 24 May 1979 (Mw 6.2)

We revisit the April 1979 Montenegro earthquake sequence to invert for finite-fault slip models for the mainshock of 15 April 1979 (M_w 7.1) and of the strongest aftershock of 24 May 1979 (M_w 6.2) using P, SH and SV waveforms, retrieved from IRIS data center. We also used body waveform modelling inversion to confirm the focal mechanism of the mainshock as a pure thrust mechanism and rule out the existence of considerable strike slip component in the motion. The mainshock occurred along a shallow (depth 7 km), low angle (14°) thrust fault, parallel to the coastline and dipping to the NE. Our preferred slip distribution model for the mainshock indicates that rupture initiated from SE and propagated towards NW, with a speed of 2.0 km/s. Moment was released in a main slip patch, confined in an area of L  50 km × W  23 km. The maximum slip ( 2.7 m) occurred  30 km to the NW of the hypocenter (location of rupture initiation). The average slip is 49 cm and the total moment release over the fault is 4.38e19 Nm. The slip model adequately fits the distribution of the M_w ≥ 4.3 aftershocks, as most of them are located in the regions of the fault plane that did not slip during the mainshock. The 24 May 1979 (M_w 6.2) strongest aftershock occurred  40 km NW of the mainshock. Our preferred slip model for this event showed a characteristic two-lobe pattern, where each lobe is  7.5 × 7.5 km². Rupture initiated in the NW lobe, where the slip obtained its maximum value of 45 cm, very close to the hypocenter, and propagated towards the south-eastern lobe where it reached another maximum value — for this lobe — of 30 cm, approximately 10 km away from the hypocenter. To indirectly validate our slip models we produced synthetic PGV maps (Shake maps) and we compared our predictions with observations of ground shaking from strong motion records. All comparisons were made for rock soil conditions and in general our slip models adequately fit the observations especially at the closest stations where the shaking was considerably stronger. Through the search of the parameter space for our inversions we obtained an optimum location for the mainshock at 42.04°N and 19.21° E and we also observed that better fit to the observations was obtained when the fault was modeled as a blind thrust fault.     

The evolution of a tropical rainforest/grassland mosaic in southeastern Brazil since 28,000 C yr BP based on carbon isotopes and pollen records

The lack of paleoecological records from the montane Atlantic Rainforest of coastal Brazil, a hotspot of biological diversity, has been a major obstacle to our understanding of the vegetational changes since the last glacial cycle. We present carbon isotope and pollen records to assess the impact of the glaciation on the native vegetation of the Serra do Mar rainforest in São Paulo, Brazil. From ca. 28,000 to  22,000 ¹⁴C yr BP, a subtropical forest with conifer trees is indicative of cool and humid conditions. In agreement carbon isotopic data on soil organic matter suggest the presence of C₃ plants and perhaps C₄ plants from  28,000 to  19,000 ¹⁴C yr BP. The significant increase in the sedimentation rate and algal spores from  19,450 to  19,000 ¹⁴C yr BP indicates increasing humidity, associated to an erosion process between  19,000 and  15,600 ¹⁴C yr BP. From  15,600 ¹⁴C yr BP to present there is a substantial increase in arboreal elements and herbs, indicating more humid and warmer climate. From  19,000 to  1000 ¹⁴C yr BP, δ¹³C values indicated the predominance of C₃ plants. These results are in agreement with studies in speleothems of caves, which suggest humid conditions during the last glacial maximum.     

Mesozoic plutons of the Yidun Arc, SW China: U/Pb geochronology and Hf isotopic signature

The Yidun Arc is a Triassic volcanic arc located between the Songpan Garzê Fold Belt and the Qiangtang Block, southwest China. To constrain the age of a number of the major granitic plutons from the Yidun Arc, laser ablation ICP-MS U/Pb analysis of zircon was conducted. Hafnium isotope data was also acquired through laser-ablation multicollector ICPMS analysis of zircon, with the aim of gaining insight into the age and nature of the source region of the plutons. Three age groups have been identified from seven granite samples: Early–Middle Triassic ( 245 to 229 Ma), Late Triassic ( 219 to 216 Ma) and Cretaceous ( 105 to 95 Ma). Hafnium analysis shows the Triassic granites to have negative and variable ε_Hf values and Mesoproterozoic ( 1.6 Ga) depleted-mantle model ages, which is interpreted to reflect derivation from an isotopically heterogeneous, largely crustal source. The Cretaceous granite shows higher and less variable ε_Hf values and slightly younger model ages ( 1.3 Ga), and is interpreted to be derived from melting of a more homogeneous crustal source. A depleted-mantle model age of  1.5 Ga is calculated from the pooled Triassic and Cretaceous samples. The source region for these magmas may be tentatively correlated with Mesoproterozoic material of the Yangtze Craton, which has been suggested to underlie the Yidun Arc; however, further work is necessary to demonstrate this suggestion.     

A late Quaternary paleotemperature record from Hanging Lake, northern Yukon Territory, eastern Beringia

The late Quaternary paleoclimate of eastern Beringia has primarily been studied by drawing qualitative inferences from vegetation shifts. To quantitatively reconstruct summer temperatures, we analyzed lake sediments for fossil chironomids, and additionally we analyzed the sediments for fossil pollen and organic carbon content. A comparison with the δ¹⁸O record from Greenland indicates that the general climatic development of the region throughout the last glaciation–Holocene transition differed from that of the North Atlantic region. Between  17 and 15 ka, mean July air temperature was on average 5°C colder than modern, albeit a period of near-modern temperature at  16.5 ka. Total pollen accumulation rates ranged between  180 and 1200 grains cm^− 2 yr^− 1. At  15 ka, approximately coeval with the Bølling interstadial, temperatures again reached modern values. At  14 ka, nearly 1000 yr after warming began, Betula pollen percentages increased substantially and mark the transition to shrub-dominated pollen contributors. Chironomid-based inferences suggest no evidence of the Younger Dryas stade and only subtle evidence of an early Holocene thermal maximum, as temperatures from  15 ka to the late Holocene were relatively stable. The most recognizable climatic oscillation of the Holocene occurred from  4.5 to 2 ka.     

Experimental frictional heating of coal gouge at seismic slip rates: Evidence for devolatilization and thermal pressurization of gouge fluids

High velocity (1 m/s) friction experiments on bituminous coal gouge display several earthquake-related phenomena, including devolatilization by frictional heating, gas pressurization, and slip weakening. Stage I is characterized by sample shortening and reduction in the coefficient of friction (μ) from  1 to 0.6. Stage II is characterized by high frequency ( 5 Hz) oscillations in stress and strain records and by gas emissions. Stage III is marked by rapid weakening (μ  0.1 to 0.35) and sample shortening, together with continued gas emissions. Stage IV produces stable stress records and continued weakness (μ  0.2), but without gas emission. Stage I shortening is due to compaction of the gouge and the weakening is attributed to mechanical or thermal effects. Stage II behavior is interpreted as due to coal gasification and fluctuations in fluid pressure, resulting in high frequency stick-slip type behavior. Dramatic reduction in shear stress in stage III is attributed to gas pressurization by pore collapse and corresponds to a frictional instability, analogous to nucleation of an earthquake. Microstructural observations indicate the deformation was brittle during stages I and II but ductile during stages III and IV. Time dependent finite element frictional heat models indicate the center of the samples became hot ( 900 °C) during stage II, whereas the edge of samples remained relatively cold (< 300 °C). Vitrinite reflectance of coal samples shows an increase in reflectance from  0.5 to  0.8% over the displacement interval 20–40 m (20–40 s), indicating that the reflectance responds to frictional heating on a short time scale. The energy expended per unit area in these low stress, large displacement experiments is similar to that of higher stress ( 50 MPa), short displacement ( 1 m) earthquakes ( 10⁷ J/m²).     

Evolution of mineralizing fluids in the Zn–Pb–Cu(–Ag ± Au) skarn and epithermal deposits of the world-class San Martín district, Zacatecas, Mexico

Mineralizing fluids at the San Martín skarn show an evolution characterized by prograde and retrograde associations. The prograde mineral associations consist of (1) a massive garnet zone, (2) a tremolite ± garnet zone, and (3) a late association of quartz, sphalerite, calcite and fluorite lining the vugs in the garnet zone. The fluids of the prograde associations exhibit decreasing temperatures of homogenization (Th) and variable salinities. The fluids of the massive garnet zone have salinities of 36 wt.% NaCl equiv. and Th of 645 to 570 °C, corresponding to pressures of 1055 bar. At the tremolite ± garnet zone, Th range from 438 to 354 °C. In the late association at the endoskarn, the following evolution can be drawn: (a) salinities of 50 to 42 wt.% NaCl equiv., and Th of 455 to 346 °C in quartz, (b) salinities of 46 wt.% NaCl equiv., and Th of 415 to 410 °C in sphalerite, (c) salinities of 50 to 37 wt.% NaCl equiv., and Th of 479 to 310 °C in calcite, (d) salinities of 33 to 28 wt.% NaCl equiv. and of 24 to 22 wt.% KCl in fluorite, and (e) two types of fluids with salinities of 2 and 39 wt.% NaCl equiv. and Th 344 and 300 °C, respectively, in later saccharoidal quartz segregations. The retrograde mineral associations comprise pervasive propylitic alteration to carbonization, and mantos with sulfides. Fluids in epidote have salinities of 7.6 wt.% NaCl equiv. and Th of 287 to 252 °C, and in calcite have salinities of 9.2 to 1 wt.% NaCl equiv. and Th of 188 to 112 °C. Fluids in the sulfide assemblages in the mantos have salinities of 8 to 3 wt.% NaCl equiv. and Th 300 °C, with corresponding pressures of 94 bar. Fluids in late epithermal veins close to the intrusive body have salinities of 10 to 5 wt.% NaCl equiv. and Th of 275 to 200 °C, and distal veins show salinities of 2 to 1 wt.% NaCl equiv. and Th of 160 °C.     

Kinematics of the Iberia–Maghreb plate contact from seismic moment tensors and GPS observations

The Iberian Peninsula and the Maghreb experience moderate earthquake activity and oblique,  NW–SE convergence between Africa and Eurasia at a rate of  5 mm/yr. Coeval extension in the Alboran Basin and a N35°E trending band of active, left-lateral shear deformation in the Alboran–Betic region are not straightforward to understand in the context of regional shortening, and evidence complexity of deformation at the plate contact. We estimate 86 seismic moment tensors (M_W 3.3 to 6.9) from time domain inversion of near-regional waveforms in an intermediate period band. Those and previous moment tensors are used to describe regional faulting style and calculate average stress tensors. The solutions associated to the Trans-Alboran shear zone show predominantly strike-slip faulting, and indicate a clockwise rotation of the largest principal stress orientation compared to the regional convergence direction (σ₁ at N350°E). At the N-Algerian and SW-Iberian margins, reverse faulting solutions dominate, corresponding to N350°E and N310°E compression, respectively. Over most of the Betic range and intraplate Iberia, we observe predominately normal faulting, and WSW–ENE extension (σ₃ at N240°E). From GPS observations we estimate that more than 3 mm/yr of African (Nubian)–Eurasian plate convergence are currently accommodated at the N-Algerian margin,  2 mm/yr in the Moroccan Atlas, and  2 mm/yr at the SW-Iberian margin. 2 mm/yr is a reasonable estimate for convergence within the Alboran region, while Alboran extension can be quantified as  2.5 mm/yr along the stretching direction (N240°E). Superposition of both motions explains the observed left-lateral transtensional regime in the Trans-Alboran shear zone. Two potential driving mechanisms of differential motion of the Alboran–Betic–Gibraltar domain may coexist in the region: a secondary stress source other than plate convergence, related to regional-scale dynamic processes in the upper mantle of the Alboran region, as well as drag from the continental-scale motion of the Nubian plate along the southern limit of the region. In the Atlantic Ocean, the  3.5 mm/yr, westward motion of the Gibraltar Arc relative to intraplate Iberia can be accommodated at the transpressive SW-Iberian margin, while available GPS observations do not support an active subduction process in this area.     

Glide systems of hematite single crystals in deformation experiments

The critical resolved shear stresses (CRSSs) of hematite crystals were determined in compression tests for r-twinning, c-twinning and {a}<m>-slip in the temperature range 25 °C to 400 °C, at 400 MPa confining pressure, and a strain rate of 10^− 5 s^− 1 by Hennig-Michaeli, Ch., Siemes, H., 1982. Experimental deformation of hematile crstals betwen 25 °C and 400 °C at 400 MPa confining pressure. In: Schreyer, W. (Ed.) High Pressure Research in Geoscience, Schweizerbart'sche Verlagsbuchhandlung, Stuttgart, p. 133–150. In the present contribution newly performed experiments on hematite single crystals at temperatures up to 800 °C at strain rates of 10^− 5 s^− 1 and 300 MPa confining pressure extends the knowledge about the CRSS of twin and slip modes. Optical observations, neutron diffraction goniometry, SEM forescatter electron images and electron backscatter diffraction are applied in order to identify the glide modes. Both twinning systems and {a}<m>-slip were confirmed by these methods. Besides the known glide systems the existence of the (c)<a>-slip system could be stated. Mechanical data establish that the CRSS of r-twinning decreases from 140 MPa at 25 °C to  5 MPa at 800 °C and for {a}<m>-slip from > 560 MPa at 25 °C to  40 MPa at 700 °C. At room temperature the CRSS for c-twinning is around 90 MPa and at 600 °C  60 MPa. The data indicate that the CRSSs above 200 °C seem to be between the values for r-twinning and {a}<m>-slip. For (c)<a>-slip only the CRSS at 600 °C could be evaluated to  60 MPa. Exact values are difficult to determine because other glide systems are always simultaneously activated.     

Development of crystallographic preferred orientations by nucleation and growth of new grains in experimentally deformed quartz single crystals

Deformation experiments have been carried out to investigate the effect of dynamic recrystallisation on crystallographic preferred orientation (CPO) development. Cylindrical samples of natural single crystals of quartz were axially deformed together with 1 vol.% of added water and 20 mg of Mn₂O₃ powder in a Griggs solid medium deformation apparatus in different crystallographic orientations with compression direction: (i) parallel to <c>, (ii) at 45° to <c> and 45° to <a> and (iii) parallel to <a>. The experiments were performed at a temperature of 800 °C, a confining pressure of 1.2 GPa, a strain rate of  10^− 6 s^− 1, to bulk finite strains of  14–36%. The deformed samples were analysed in detail using optical microscopy, electron backscatter diffraction (EBSD) and transmission electron microscopy (TEM). Two different microstructural domains were distinguished in the deformed samples: (i) domains with undulatory extinction and deformation lamellae, and (ii) domains with new recrystallised grains. Within the domains of undulatory extinction, crystal-plastic deformation caused gradual rotations of the crystal lattice up to  30° away from the host orientation. New recrystallised grains show a strong CPO with c-axis maxima at  45° to the compression direction. This is the case in all experiments, irrespective of the initial crystallographic orientation. The results show that c-axes are not continuously rotated towards the new maxima. The new grains thus developed through a mechanism different from subgrain rotation recrystallisation. New grains have a subeuhedral shape and numerous microcavities, voids, fluid channels and fluid inclusions at their grain boundaries. No host control is recorded in misorientation axes across their large angle grain boundaries. New grains might have been created by nucleation from solution in the μm-scale voids and microfractures. The CPO most likely developed due to preferred growth of the freshly precipitated grains with orientations suitable for intracrystalline deformation at the imposed experimental conditions.     

Neoproterozoic adakitic plutons in the northern margin of the Yangtze Block, China: Partial melting of a thickened lower crust and implications for secular crustal evolution

Numerous Neoproterozoic felsic and mafic–ultramafic intrusions occur in the Hannan region at the northern margin of the Yangtze Block. Among these, the Wudumen and Erliba plutons consist of granodiorites and have SHRIMP zircon U–Pb ages of  735 Ma. The rocks have high K₂O (0.8–3.6 wt.%) and Na₂O (4.4–6.4 wt.%) and low MgO (0.4–1.7 wt.%). They also have high Sr/Y (32–209) and (La/Yb)_n ratios (4.4–38.6). Their εNd values range from − 0.41 to − 0.92 and zircon initial ¹⁷⁶Hf/¹⁷⁷Hf ratios from 0.282353 to 0.282581. These geochemical features are similar to those of adakitic rocks produced by partial melting of a thickened lower crust. Our new analytical results, combined with the occurrence of voluminous arc-related mafic–ultramafic intrusions emplaced before 740 Ma, lead us to propose that the crustal evolution in the northern margin of the Yangtze Block during Neoproterozoic involved: (1) rapid crustal growth and thickening by underplating of mafic magmas from the mantle which was modified by materials coming from the subducting oceanic slab from  1.0 to  0.74 Ga, and (2) partial melting of the thickened lower crust due to a thermal anomaly induced by upwelling of asthenosphere through an oceanic slab window, producing the  735 Ma adakitic Wudumen and Erliba plutons. Our model suggests that the crustal thickness was more than 50 km at the northern margin of the Yangtze Block at  735 Ma, and rule out the possibility of a mantle plume impact causing the > 735 Ma magmatism in the region.     

Brittle–ductile–brittle deformation during cooling of tonalite (Adamello, Southern Italian Alps)

Late- to post-magmatic deformation in slightly diachronous contiguous intrusions of the north-western Adamello batholith (Southern Alps, Italy) is recorded as, from oldest to youngest: (i) joints, (ii) solid-state ductile shear zones, (iii) faults associated with epidote-K-feldspar veins and (iv) zeolite veins and faults. Structures (ii) to (iv) are localized on the pervasive precursory network of joints (i), which developed during the earliest stages of pluton cooling. High temperature ( 500 °C), ductile overprinting of joints produced lineations, defined by aligned biotite and hornblende, on the joint surfaces and highly localized mylonites. The main phase of faulting, producing cataclasites and pseudotachylytes, occurred at  250 °C and was associated with extensive fluid infiltration. Cataclasites and pseudotachylytes are clustered along different E–W-striking dextral strike-slip fault zones correlated with the activity of the Tonale fault, a major tectonic structure that bounds the Adamello batholith to the north. Ductile deformation and cataclastic/veining episodes occurred at P = 0.25–0.3 GPa during rapid cooling of the batholith to the ambient temperatures ( 250 °C) that preceded the exhumation of the batholith. Timing of the sequence of deformation can be constrained by ³⁹Ar–⁴⁰Ar ages of  30 Ma on pseudotachylytes and various existing mineral ages. In the whole composite Adamello batholith, multiple magma pulses were intruded over the time span 42–30 Ma and each intrusive body shows the same ductile-to-brittle structural sequence localized on the early joint sets. This deformation sequence of the Adamello might be typical of intrusions undergoing cooling at depths close to the brittle–ductile transition.     

Frictional sliding of gabbro gouge under hydrothermal conditions

We investigated the frictional sliding behaviour of gabbro gouge under hydrothermal conditions. Experiments were performed on 1-mm-thick gabbro gouge sandwiched between country rock pieces (with gouge inclined 35° to the sample axis) in a triaxial testing system with argon gas as the confining medium. In the first series, experiments were conducted under effective normal stresses of 200 MPa and 300 MPa respectively, with pore pressure of 10 MPa. For temperature over 400 °C, pore pressure of 30 MPa was also applied to implement supercritical water conditions. At temperatures up to 615 °C, slip rate steps ranging from 0.0488 μm/s to 1.22 μm/s were applied to obtain the rate dependence of friction.At 200 MPa effective normal stress and a pore pressure of 10 MPa, the steady state rate dependence a–b shows velocity-weakening behaviour for temperatures between  200 and  310 °C. The higher temperature limit for velocity-weakening behaviour to occur extends up to  510 °C under supercritical water conditions with a pore pressure of 30 MPa. For the limited sliding distance in our experiments, only velocity-strengthening behaviour occurred at 300 MPa effective normal stress. Considering the limited displacement (< 3.5 mm), velocity-weakening behaviour may not be excluded in the high effective normal stress case for temperature below  510 °C.The coefficient of friction shows an increasing trend with increasing temperature in the low temperature range. The cut-off temperatures for the increasing trend are  250 °C and  440 °C, respectively for the 200 MPa and 300 MPa effective normal stress cases. Above the cut-off temperatures, the coefficient of friction at 1.83 mm permanent displacement varies around an average of 0.73, which is identical to the average for the oven-dried case [He, C., Yao, W., Wang, Z., Zhou, Y., 2006. Strength and stability of frictional sliding of gabbro gouge at elevated temperatures. Tectonophysics 427, 217–229, doi:10.1016/j.tecto.2006.05.023]. Together with the small value of rate dependence (a–b < 0.0073) for the whole temperature range, these results indicate the absence of fluid-assisted creep.With the result of our experiments as a constraint on strength of frictional sliding, comparison between converted strength for strike–slip faults and creep strength of gabbro-like rocks implies fracturing and faulting behaviours in the lower crust of a cool area (Zhangbei) in North China.     

Cretaceous and Cenozoic cooling history across the ultrahigh pressure Tongbai–Dabie belt, central China, from apatite fission-track thermochronology

The crystalline terrane of the Tongbai–Dabie region, central China, comprising the Earth's largest ultrahigh-pressure (UHP) exposure was formed during Triassic collision between the Sino–Korean and Yangtze cratons. New apatite fission-track (AFT) data presented here from the UHP terrane, extends over a significantly greater area than reported in previous studies, and includes the (eastern) Dabie, the Hong'an (northwestern Dabie) and Tongbai regions. The new data yield ages ranging from 44 ± 3 to 142 ± 36 Ma and mean track lengths between 10 and 14.4 μm. Thermal history models based on the AFT data taken together with published ⁴⁰Ar/³⁹Ar, K–Ar, apatite and zircon (U–Th)/He and U–Pb data, exhibit a three-stage cooling pattern that is similar across the study region, commencing with an Early Cretaceous rapid cooling event, followed by a period of relative thermal stability during which rocks remained at temperatures within the AFT partial annealing zone (60–110 °C) and ending with a possible renewed phase of accelerated cooling during Pliocene to Recent time. The first cooling phase followed large-scale transtensional deformation between 140 and 110 Ma and is related to Early Cretaceous eastward tectonic escape and Pacific back arc extension. Between this phase and the subsequent slow cooling phase, a transition period from 120 to 80 Ma (to 70 to 45 Ma along the Tan–Lu fault) was characterised by a relatively low cooling rate (3–5 °C/Ma). This transition is likely related to a tectonic response associated with the mid-Cretaceous subduction of the Izanagi–Pacific plate as well as lithospheric extension and thinning in eastern Asia. The present regional AFT age pattern is therefore basically controlled by the Early Cretaceous rapid cooling event, but finally shaped through active Cenozoic faulting. Following the transition phase the subsequent slow cooling phase pattern implies a net reduction in horizontal compressional stress corresponding to increased extension rates along the continental margin due to the decrease in plate convergence. Modelling of the AFT data suggests a possible Pliocene–Recent cooling episode, which may be supported by increased rates of sedimentation observed in adjacent basins. This cooling phase may be interpreted as a response to the far-field effects of the frontal India–Eurasia collision to the west. Approximate estimates suggest that the total amount of post 120 Ma denudation across the UHP orogen ranged from 2.4 to 13.2 km for different tectonic blocks and ranged from 0.8 to 9.7 km during the Cretaceous to between 1.7 and 3.8 km during the Cenozoic.     

Paleomagnetic constraints on the paleogeography and oroclinal bending of the Devonian volcanic arc in Kazakhstan

Numerical modelling, incorporating coupling between surface processes and induced flow in the lower continental crust, is used to address the Quaternary evolution of the Gulf of Corinth region in central Greece. The post-Early Pleistocene marine depocentre beneath this Gulf overlies the northern margin of an older (Early Pleistocene and earlier) lacustrine basin, the Proto Gulf of Corinth Basin or PGCB. In the late Early Pleistocene, relief in this region was minimal but, subsequently, dramatic relief has developed, involving the creation of  900 m of bathymetry within the Gulf and the uplift by many hundreds of metres of the part of the PGCB, south of the modern Gulf, which forms the Gulf's main sediment supply. It is assumed that, as a result of climate change around 0.9 Ma, erosion of this sediment source region and re-deposition of this material within the Gulf began, both processes occurring at spatial average rates of  0.2 mm a^− 1. Modelling of the resulting isostatic response indicates that the local effective viscosity of the lower crust is  4 × 10¹⁹ Pa s, indicating a Moho temperature of  560 °C. It predicts that the  10 mm a^− 1 of extension across this  70 km wide model region, at an extensional strain rate of  0.15 Ma^− 1, is partitioned with  3 mm a^− 1 across the sediment source,  2 mm a^− 1 across the depocentre, and  5 mm a^− 1 across the ‘hinge zone’ in between, the latter value being an estimate of the extension rate on normal faults forming the major topographic escarpment at the southern margin of the Gulf. This modelling confirms the view, suggested previously, that coupling between this depocentre and sediment source by lower-crustal flow can explain the dramatic development in local relief since the late Early Pleistocene. The effective viscosity of the lower crust in this region is not particularly low; the strong coupling interpreted between the sediment source and depocentre results instead from their close proximity. In detail, the effective viscosity of the lower crust is expected to decrease northward across this model region, due to the northward increase in exposure of the base of the continental lithosphere to the asthenosphere; in the south the two are separated by the subducting Hellenic slab. The isostatic consequences of such a lateral variation in viscosity provide a natural explanation for why, since  0.9 Ma, the modern Gulf has developed asymmetrically over the northern part of the PGCB, leaving the rest of the PGCB to act as its sediment source.     

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.     

Guidelines for assessing the petroleum potential of coaly source rocks using Rock-Eval pyrolysis

Coaly source rocks are sufficiently different from marine and lacustrine source rocks in their organic matter characteristics to warrant separate guidelines for their assessment using Rock-Eval pyrolysis. The rank threshold for oil generation is indicated by the increase in BI (S1/TOC) at Rank(S_r)9–10 (T_max 420–430 °C, R_o 0.55–0.6%), and the threshold for oil expulsion is indicated by the peak in QI ([S1+S2]/TOC) at Rank(S_r)11–12.5 (T_max 430–440 °C, R_o 0.65–0.85%). The pronounced rank-related increase in HI (S2/TOC) prior to oil expulsion renders the use of immature samples inappropriate for source rock characterisation. A more realistic indication of the petroleum generative potential and oil expulsion efficiency of coaly source rocks can be gained from samples near the onset of expulsion. Alternatively, effective HI′ values (i.e. HIs near the onset of expulsion) can be estimated by translating the measured HIs of immature samples along the maturation pathway defined by the New Zealand (or other defined) Coal Band. Coaly source rocks comprise a continuum of coaly lithologies, including coals, shaly coals and coaly mudstones. Determination of the total genetic potential of coaly source rock sequences is best made using lithology-based samples near the onset of expulsion.     

Palaeomagnetism of the Loch Doon Granite Complex, Southern Uplands of Scotland: The Late Caledonian palaeomagnetic record and an Early Devonian episode of True Polar Wander

The Southern Uplands terrane is an Ordovician–Silurian back-arc/foreland basin emplaced at the northern margin of the Iapetus Ocean and intruded by granite complexes including Loch Doon (408.3 ± 1.5 Ma) during Early Devonian times. Protracted cooling of this 130 km³ intrusion recorded magnetic remanence comprising a predominant (‘A’) magnetisation linked to initial cooling with dual polarity and mean direction D / I = 237 / 64° (α₉₅ = 4°, palaeopole at 316°E, 21°N). Subsidiary magnetisations include Mesozoic remanence correlating with extensional tectonism in the adjoining Irish Sea Basin (‘B’, D / I = 234/− 59°) and minority populations (‘C’, D / I = 106/− 2° and ‘D’, D / I = 199/1°) recording emplacement of younger ( 395 Ma) granites in adjoining terranes and the Variscan orogenic event. The ‘A’ directions have an arcuate distribution identifying anticlockwise rotation during cooling. A comparable rotation is identified in the Orthotectonic Caledonides to the north and the Paratectonic Caledonides to the south following closure of Iapetus. Continental motion from midsoutherly latitudes ( 40°S) at 408 Ma to equatorial palaeolatitudes by  395 Ma is identified and implies minimum rates of continental movement between 430 and 390 Ma of 30–70 cm/year, more than double maximum rates induced by plate forces and interpreted as a signature of true polar wander. Silurian–Devonian palaeomagnetic data from the British–Scandinavian Caledonides define a 430–385 Ma closed loop comparable to the distributed contemporaneous palaeomagnetic poles from Gondwana. They reconcile pre-430 Ma and post-380 Ma APW from this supercontinent and show that Laurentia–Baltica–Avalonia lay to the west of South America with a relict Rheic Ocean opening to the north which closed to produce Variscan orogeny by a combination of pivotal closure and right lateral transpression.     

Insights into Wasatch fault vertical slip rates using the age of sediments in Timpanogos Cave, Utah

Timpanogos Cave, located near the Wasatch fault, is about 357 m above the American Fork River. Fluvial cave sediments and an interbedded carbonate flowstone yield a paleomagnetic and U–Th depositional age of 350 to 780 ka. Fault vertical slip rates, inferred from calculated river downcutting rates, range between 1.02 and 0.46 mm yr^− 1. These slip rates are in the range of the 0–12 Ma Wasatch Range exhumation rate ( 0.5–0.7 mm yr^− 1), suggesting that the long-term vertical slip rate remained stable through mid-Pleistocene time. However, the late Pleistocene (0–250 ka) decelerated slip rate ( 0.2–0.3 mm yr^− 1) and the accelerated Holocene slip rate ( 1.2 mm yr^− 1) are consistent with episodic fault activity. Assuming that the late Pleistocene vertical slip rate represents an episodic slowing of fault movement and the long-term (0–12 Ma) average vertical slip rate, including the late Pleistocene and Holocene, should be  0.6 mm yr^− 1, there is a net late Pleistocene vertical slip deficit of  50–75 m. The Holocene and late Pleistocene slip rates may be typical for episodes of accelerated and slowed fault movement, respectively. The calculated late Pleistocene slip deficit may mean that the current accelerated Wasatch fault slip rate will extend well into the future.     

Linking source and sink: Evaluating the balance between onshore erosion and offshore sediment accumulation since Gondwana break-up, South Africa

The geomorphic origin and evolution of the tectonically unique interior highland of southern Africa, the Kalahari Plateau, and its flanking low-lying coastal planes, remain largely unresolved because of a lack of regional quantitative analyses of its uplift and erosion history. Here we focus on the southern Cape, South Africa and link onshore denudation, based on new apatite fission track thermochronology results, to offshore sediment accumulation, using abundant well data and a seismic reflection profile. We attempt to relate source and sink in order to resolve some first order issues concerning timing of the exhumation and development of the topographic features of southern Africa. The volume of sediment accumulated off South Africa's south coast is calculated using 173 wells and a seismic reflection profile. A total, uncompacted, sediment volume of 268,500 km³ accumulated off South Africa's south coasts since  136 Ma, in the Outeniqua and Southern Outeniqua Basins. Accumulation volumes and rates were highest in the early Cretaceous (48,800 × 10⁴ km³ at  8150 km³/Ma from  136 to 130 Ma, and 57,500 × 10⁴ km³ at 5750 km³/Ma from  130 to 120 Ma) and mid–late Cretaceous (83,700 × 10⁴ km³ at 3200 km³/Ma from  93 to 67 Ma). Volumes and accumulation rates were lowest for the early–mid-Cretaceous (47,400 × 10⁴ km³ at 1750 km³/Ma from  120 to 93 Ma) and the Cenozoic (31,200 × 10⁴ km³ at 450 km³/Ma from  67 to 0 Ma). Although our analysis shows that the accumulated volume of offshore sediments does not match the calculated volume of onshore erosion, as quantified through apatite fission track thermochronology (e.g. Tinker, J.H., de Wit, M.J., Brown, R., 2008. Mesozoic exhumation of the 439 southern Cape, South Africa, quantified using apatite fission track thermochronology. Tectonophysics, doi: 10.1016/j.tecto.2007.10.009), the timing of increased sediment accumulation closely matches the timing of increased onshore denudation. This suggests that the greatest volumes of material were transported from source to sink during two distinct Cretaceous episodes, and that the processes driving onshore denudation decreased by an order of magnitude during the Cenozoic.     

Anatolian surface wave evaluated at GEOFON Station ISP Isparta, Turkey

We have studied seismic surface waves of 255 shallow regional earthquakes recently recorded at GEOFON station ISP (Isparta, Turkey) and have selected these 52 recordings with high signal-to-noise ratio for further analysis. An attempt was made by the simultaneous use of the Rayleigh and Love surface wave data to interpret the planar crust and uppermost mantle velocity structure beneath the Anatolian plate using a differential least-square inversion technique. The shear-wave velocities near the surface show a gradational change from approximately 2.2 to 3.6 km s^− 1 in the depth range 0–10 km. The mid-crustal depth range indicating a weakly developed low velocity zone has shear-wave velocities around 3.55 km s^− 1. The Moho discontinuity characterizing the crust–mantle velocity transition appears somewhat gradual between the depth range  25–45 km. The surface waves approaching from the northern Anatolia are estimated to travel a crustal thickness of  33 km whilst those from the southwestern Anatolia and part of east Mediterranean Sea indicate a thicker crust at  37 km. The eastern Anatolia events traveled even thicker crust at  41 km. A low sub-Moho velocity is estimated at  4.27 km s^− 1, although consistent with other similar studies in the region. The current velocities are considerably slower than indicated by the Preliminary Reference Earth Model (PREM) in almost all depth ranges.