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.     

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.     

Geochemistry of the Jurassic Mirdita Ophiolite (Albania) and the MORB to SSZ evolution of a marginal basin oceanic crust

The Middle Jurassic Mirdita Ophiolite in northern Albania is part of an ophiolite belt occurring between the Apulian and Pelagonian subcontinents in the Balkan Peninsula. The upper mantle and crustal units of the Mirdita Ophiolite show major changes in thickness, rock types, and chemical compositions from west to east as a result of its complex evolution in a suprasubduction zone (SSZ) environment. The  3–4-km-thick Western Mirdita Ophiolite (WMO) includes lherzolite–harzburgite, plagioclase–lherzolite, plagioclase–dunite in its upper mantle units and a plutonic complex composed of olivine gabbro, troctolite, ferrogabbro, and gabbro. These peridotites and gabbroic rocks are overlain directly by a  600-m-thick extrusive sequence containing basaltic pillow lavas and hyaloclastites. Sheeted dikes are rare in the WMO. The  12-km-thick Eastern Mirdita Ophiolite (EMO) includes tectonized harzburgite and dunite with extensive chromite deposits, as well as ultramafic cumulates including olivine clinopyroxenite, wehrlite, olivine websterite, and dunite forming a transitional Moho with the overlying lower crustal section. The plutonic rocks are made of pyroxenite, gabbronorite, gabbro, amphibole gabbro, diorite, quartz diorite, and plagiogranite. A well-developed sheeted dike complex has mutually intrusive relations with the underlying isotropic gabbros and plagiogranites and feeds into the overlying pillow lavas. Dike compositions change from older basalt to basaltic andesite, andesite, dacite, quartz diorite, to late-stage andesitic and boninitic dikes as constrained by crosscutting relations. The  1.1-km-thick extrusive sequence comprises basaltic and basaltic andesitic pillow lavas in the lower 700 m, and andesitic, dacitic and rhyodacitic massive sheet flows in the upper 400 m. Rare boninitic dikes and lavas occur as the youngest igneous products within the EMO. The basaltic and basaltic andesitic rocks of the WMO extrusive sequence display MORB affinities with Ti and Zr contents decreasing upsection (TiO₂ = 3.5–0.5%, Zr = 300–50 ppm), while _Nd(T) (+ 8 to + 6.5) varies little. These magmas were derived from partial melting of fertile MORB-type mantle. Fractional crystallization was important in the evolution of WMO magmas. The low Ti and HREE abundances and Cs and Ba enrichments in the uppermost basaltic andesites may indicate an increased subduction influence in the evolution of the late-stage WMO magmas. Basaltic andesites in the lower 700 m of the EMO volcanic sequence have lower TiO₂ ( 0.5%) and Zr ( 50 ppm) contents but _Nd(T) values (+ 7 to + 6.5) are similar to those of the WMO lavas. These rocks show variable enrichment in subduction-enriched incompatible elements (Cs, Ba, Th, U, LREE). The basaltic andesites through dacites and boninites within the upper 400 meters of EMO lavas show low TiO₂ ( 0.8–0.3%) and _Nd(T) (+ 6.5 to + 3.0). The mantle source of these rocks was variably enriched in Th by melts derived from subducted sediments as indicated by the large variations in Ba, K, and Pb contents. EMO boninitic dikes and lavas and some gabbroic intrusions with negative _{Nd (T)} values (− 1.4 and − 4.0, respectively) suggest that these magmas were produced from partial melting of previously depleted, ultra-refractory mantle. The MORB to SSZ transition (from west to east and stratigraphically upwards in the Mirdita Ophiolite and the progression of the _Nd(T) values from + 8.0 to − 4.0 towards the east resulted from an eastward shift in protoarc–forearc magmatism, keeping pace with slab rollback in this direction. The mantle flow above the retreating slab and in the arc-wedge corner played a major role in the evolution of the melting column, in which melt generation, aggregation/mixing and differentiation occurred at all levels of the sub-arc/forearc mantle. The SSZ Mirdita Ophiolite evolved during the intra-oceanic collapse and closure of the Pindos marginal basin, which had a protracted tectonic history involving seafloor spreading, protoarc rifting, and trench-continent collision.     

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²).     

A late Quaternary pollen sequence from Mfabeni Peatland, South Africa: Reconstructing forest history in Maputaland

This paper documents a continuous  44,000-yr pollen record derived from the Mfabeni Peatland on the Maputaland Coastal Plain. A detailed fossil pollen analysis indicates the existence of extensive Podocarpus-abundant coastal forests before  33,000 cal yr BP. The onset of wetter local conditions after this time is inferred from forest retreat and the development of swampy conditions. Conditions during the last glacial maximum ( 21,000 cal yr BP) are inferred to have been colder and drier than the present, as evidenced by forest retreat and replacement of swampy reed/sedge communities by dry grassland. Forest growth and expansion during the Holocene Altithermal ( 8000–6000 cal yr BP) indicates warm, relatively moist conditions. Previous records from Maputaland have suggested a northward migration of Podocarpus forest during the late Holocene. However, we interpret a mid-Holocene decline in Podocarpus at Mfabeni as evidence of deforestation. Forest clearance during the mid-Holocene is supported by the appearance of Morella serrata, suggesting a shift towards more open grassland/savanna, possibly due to burning. These signals of human impact are coupled with an increase in Acacia, indicative of the development of secondary forest and hence disturbance.     

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.     

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.     

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.     

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.     

Cirque glacier activity in arctic Norway during the last deglaciation

Numerous cirques of the Lofoten–Vesterålen archipelago in northern Norway have distinct moraine sequences that previously have been assigned to the Allerød-Younger Dryas ( 13,400 to 11,700 yr BP) interval, constraining the regional distribution of the equilibrium-line altitude (ELA) of cirque and valley glaciers. Here we present evidence from a once glacier-fed lake on southern Andøya that contests this view. Analyses of radiocarbon dated lacustrine sediments including rock magnetic parameters, grain size, organic matter, dry bulk density and visual interpretation suggest that no glacier was present in the low-lying cirque during the Younger Dryas-Allerød. The initiation of the glacial retreat commenced with the onset of the Bølling warming ( 14,700 yr BP) and was completed by the onset of Allerød Interstade ( 13,400 yr BP). The reconstructed glacier stages of the investigated cirque coincide with a cool and dry period from  17,500 to 14,700 yr BP and a somewhat larger Last Glacial Maximum (LGM) advance possibly occurring between  21,050 and 19,100 yr BP.     

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.     

A late-glacial transition from Picea glauca to Picea mariana in southern New England

Picea is an important taxon in late-glacial pollen records from eastern North America, but little is known about which species of Picea were present. We apply a recently developed palynological method for discriminating the three Picea species in eastern North America to three records from New England. Picea glauca was dominant at  14,500–14,000 cal yr BP, followed by a transition to Picea mariana between  14,000 and 13,500 cal yr BP. Comparison of the pollen data with hydrogen isotope data shows clearly that this transition began before the beginning of the Younger Dryas Chronozone. The ecological changes of the late-glacial interval were not a simple oscillation in the position of a single species' range, but rather major changes in vegetation structure and composition occurring during an interval of variations in several environmental factors, including climate, edaphic conditions, and atmospheric CO₂ levels.     

Age estimates of coastal terraces in the Andaman and Nicobar Islands and their tectonic implications

The great Indian Ocean earthquake of December 26, 2004 caused significant vertical changes in its rupture zone. About 800 km of the rupture is along the Andaman and Nicobar Islands, which forms the outer arc ridge of the subduction zone. Coseismic deformation along the exposed land could be observed as uplift/subsidence. Here we analyze the morphological features along the coast of the Andaman and Nicobar Islands, in an effort to reconstruct the past tectonics, taking cues from the coseismic effects. We obtained radiocarbon dates from coastal terraces of the island belt and used them to compute uplift rates, which vary from 1.33 mm yr^− 1 in the Little Andaman to 2.80 mm yr^− 1 in South Andaman and 2.45 mm yr^− 1 in the North Andaman. Our radiocarbon dates converge on  600 yr and  1000 yr old coastal uplifts, which we attribute to the level changes due to two major previous subduction earthquakes in the region.     

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.     

Analysis of crustal deformation in Luzon, Philippines using geodetic observations and earthquake focal mechanisms

We utilize regional GPS velocities from Luzon, Philippines, with focal mechanism data from the Harvard Centroid Moment Tensor (CMT) Catalog, to constrain tectonic deformation in the complex plate boundary zone between the Philippine Sea Plate and Eurasia (the Sundaland block). Processed satellite imagery and digital elevation models are used with existing gravity anomaly, seismicity, and geologic maps to define a suite of six elastic blocks. Geodetic and focal mechanism data are inverted simultaneously to estimate plate rotations and fault-locking parameters for each of the tectonic blocks and faults comprising Luzon. Major tectonic structures that were found to absorb the plate convergence include the Manila Trench (20–100 mm yr^− 1) and East Luzon Trough ( 9–15 mm yr^− 1)/Philippine Trench ( 29–34 mm yr^− 1), which accommodate eastward and westward subduction beneath Luzon, respectively; the left-lateral strike-slip Philippine Fault ( 20–40 mm yr^− 1), and its northward extensions, the Northern Cordillera Fault ( 17–37 mm yr^− 1 transtension), and the Digdig Fault ( 17–27 mm yr^− 1 transpression). The Macolod Corridor, a zone of active volcanism, crustal thinning, extension, and extensive normal and strike-slip faulting in southwestern Luzon, is associated with left-lateral, transtensional slip of  5–10 mm yr^− 1. The Marikina Fault, which separates the Central Luzon block from the Southwestern Luzon block, reveals  10–12 mm yr^− 1 of left-lateral transpression. Our analysis suggests that much of the Philippine Fault and associated splays are locked to partly coupled, while the Manila and Philippine trenches appear to be poorly coupled. Luzon is best characterized as a tectonically active plate boundary zone, comprising six mobile elastic tectonic blocks between two active subduction zones. The Philippine Fault and associated intra-arc faults accommodate much of the trench-parallel component of relative plate motion.     

Pre-seismic electromagnetic phenomena in the framework of percolation and fractal theories

Understanding the electromagnetic response to geodynamic processes occurring in the earth's upper crust, in particular pre-seismic and seismic processes, is a challenging task in modern geophysics. There is increasing evidence that seismo-electromagnetic (SEM) phenomena are difficult to describe quantitatively by “linear” models using “averaged” parameters of the medium, such as electrical and hydraulic conductivities. Because the upper crust is highly inhomogeneous (at all scales), porous, and can be fully or partially water-saturated, the most natural way to describe its parameters is via fractal-theoretic and percolation-theoretic models.Recent studies indicate that the electrokinetic effect is the most likely driving mechanism for the various types of SEM signals. Here we considered the hydraulic, electric, and electrokinetic conductances of a porous water-saturated medium as a function of porosity () and moisture content (θ), utilizing a percolation/fractal approach. We show that the electric conductivity and electrokinetic current in such a medium are both proportional to (–_c)² and (θ–θ_c)², where _c and θ_c are the critical values of porosity and moisture content, respectively. This behavior admits the possibility of a relatively large change in the respective electric and electrokinetic parameters due to a small change in the mechanical strain field. This is significant because it may account for the appearance of some types of SEM signals at large distances from the earthquake origin, which is the main deficiency of most models. Indeed, the anticipated strain changes related to pre-earthquake processes is usually very small except near the focal area. The expected ‘averaged’ electromagnetic response also would be very small, unless a local underground water system exists, not necessarily near the focal area, but which is close to critical point(s).We discuss the conditions under which electrotelluric and geomagnetic variations can accompany mechanical disturbances in the earth's crust.     

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.     

Seismic properties of the upper crust in the central Friuli area (northeastern Italy) based on petrophysical data

The compressional and shear wave velocities have been measured at room temperature and pressure up to 450 MPa on 5 sedimentary rock samples, representative of the most common lithologies of the upper crust in the central Friuli area (northeastern Italy). At 400 MPa confining pressure the Triassic dolomitic rock shows the highest velocities (Vp  7 km/s, Vs  3.6 km/s), the Jurassic and Triassic limestones samples intermediate velocities (Vp  6.3 /s, Vs  3.5 km/s) and the Cenozoic and Paleozoic sandstones the lowest velocities (Vp  6.15 km/s, Vs  3.35 km/s). The Paleozoic sandstone sample is characterized by the strongest anisotropy (10%) and significant birefringence (0.2 km/s) is found only on the Cenozoic sandstone sample. We elaborated the synthetic profiles of seismic velocities, density, elastic parameters and reflection coefficient, related to 4 one-dimensional geological models extended up to 22 km depth. The synthetic profiles evidence high rheological contrasts between Triassic dolomitic rocks and the soft sandstones and the Jurassic limestones. The Vp profiles obtained from laboratory measurements match very well the in-situ Vp profile measured by sonic log for the limestones and dolomitic rocks, supporting our one-dimensional modelling of the calcareous-carbonatic stratigraphic series. The Vp and Vs values of the synthetic profiles are compared with the corresponding ones obtained from the 3-D tomographic inversion of local earthquakes. The laboratory Vp are generally higher than the tomographic ones with major discrepancies for the dolomitic lithology. The comparison with the depth location of seismicity reveals that the seismic energy is mainly released in correspondence of high-contrast rheological boundaries.     

Climate variability in the Aral Sea basin (Central Asia) during the late Holocene based on vegetation changes

High-resolution pollen analyses ( 50 yr) from sediment cores retrieved at Chernyshov Bay in the NW Large Aral Sea record shifts in vegetational development from subdesertic to steppe vegetation in the Aral Sea basin during the late Holocene. Using pollen data to quantify climatic parameters, we reconstruct and date for the first time significant changes in moisture conditions in Central Asia during the past 2000 yr. Cold and arid conditions prevailed between ca. AD 0 and 400, AD 900 and 1150, and AD 1500 and 1650 with the extension of xeric vegetation dominated by steppe elements. These intervals are characterized by low winter and summer mean temperatures and low mean annual precipitation (P_mm < 250 mm/yr). Conversely, the most suitable climate conditions occurred between ca. AD 400 and 900, and AD 1150 and 1450, when steppe vegetation was enriched by plants requiring moister conditions (P_mm  250–500 mm/yr) and some trees developed. Our results are fairly consistent with other late Holocene records from the eastern Mediterranean region and the Middle East, showing that regional rainfall in Central Asia is predominantly controlled by the eastern Mediterranean cyclonic system when the North Atlantic Oscillation (NAO) is in a negative phase.