Global trench migration velocities and slab migration induced upper mantle volume fluxes: Constraints to find an Earth reference frame based on minimizing viscous dissipation

Since the advent of plate tectonics different global reference frames have been used to describe the motion of plates and trenches. The difference in plate motion and trench migration between different reference frames can be substantial (up to 4 cm/yr). This study presents an overview of trench migration velocities for all the mature and incipient subduction zones on Earth as calculated in eight different global reference frames. Calculations show that, irrespective of the reference frame: (1) trench retreat always dominates over trench advance, with 62–78% of the 244 trench segments retreating; (2) the mean and median trench velocity are always positive (retreating) and within the range 1.3–1.5 cm/yr and 0.9–1.3 cm/yr, respectively; (3) rapid trench retreat is only observed close to lateral slab edges (< 1500 km); and (4) trench retreat is always slow far from slab edges (> 2000 km). These calculations are predicted by geodynamic models with a varying slab width, in which plate motion, trench motion and mantle flow result from subduction of dense slabs, suggesting that trench motion is indeed primarily driven by slab buoyancy forces and that proximity to a lateral slab edge exerts a dominant control on the trench migration velocity. Despite these four general conclusions, significant differences in velocities between such reference frames remain. It is therefore important to determine which reference frame most likely describes the true absolute velocities to get an understanding of the forces driving plate tectonics and mantle convection. It is here proposed that, based on fluid dynamic considerations and predictions from geodynamic modelling, the best candidate is the one, which optimises the number of trench segments that retreat, minimizes the trench–perpendicular trench migration velocity (v_T) in the centre of wide (> 4000 km) subduction zones, maximizes the number of retreating trench segments located within 2000 km of the closest lateral slab edge, minimizes the average of the absolute of the trench–perpendicular trench migration velocity (|v_T|) for all subduction zones on Earth, and minimizes the global upper mantle toroidal volume flux (_To) that results from trench migration and associated lateral slab migration (i.e. slab rollback or slab roll-forward). Calculations show that these conditions are best met in one particular Indo-Atlantic hotspot reference frame, where 75% of the subduction zones retreat, v_T in the centre of wide subduction zones ranges between − 3.5 and 1.8 cm/yr, 83% of the trench segments located within 2000 km of the closest lateral slab edge retreat, the average of |v_T| is 2.1 cm/yr, and _To = 456 km³/yr (lower limit) and 539 km³/yr (upper limit). Inclusion of all the incipient subduction zones on Earth results in slightly greater fluxes of 465 km³/yr (lower limit) and 569 km³/yr (upper limit). It is also found that this reference frame is close to minimizing the total sub-lithospheric upper mantle volume flux (_K) associated with motion of continental keels located below the major cratons. It is stressed, however, that _K is an order of magnitude smaller than _To, and thus of subordinate importance. In conclusion, the Indo-Atlantic hotspot reference frame appears preferable for calculating plate velocities and plate boundary velocities.     

Equilibrium H/H fractionations in organic molecules. II: Linear alkanes, alkenes, ketones, carboxylic acids, esters, alcohols and ethers

Equilibrium ²H/¹H fractionation factors (α_eq) for various H positions in alkanes, alkenes, ketones, carboxylic acids, esters, alcohols, and ethers were calculated between 0 and 100 °C using vibrational frequencies from ab initio QM calculations (B3LYP/6-311G**). Results were then corrected using a temperature-dependent linear calibration curve based on experimental data for H_α in ketones (Wang et al., 2009). The total uncertainty in reported α_eq values is estimated at 10–20‰. The effects of functional groups were found to increase the value of α_eq for H next to electron-donating groups, e.g. OR, OH or O(CO)R, and to decrease the value of α_eq for H next to electron-withdrawing groups, e.g. (CO)R or (CO)OR. Smaller but significant functional group effects are also observed for H_β and sometimes H_γ. By summing over individual H positions, we estimate the equilibrium fractionation relative to water to be −90‰ to −70‰ for n-alkanes and around −100‰ for pristane and phytane. The temperature dependence of these fractionations is very weak between 0 and 100 °C. Our estimates of α_eq agree well with field data for thermally mature hydrocarbons (δ²H values between −80‰ and −110‰ relative to water). Therefore the observed δ²H increase of individual hydrocarbons and the disappearance of the biosynthetic δ²H offset between n-alkyl and linear isoprenoid lipids during maturation of organic matter can be confidently attributed to H exchange towards an equilibrium state. Our results also indicate that many n-alkyl lipids are biosynthesized with δ²H values that are close to equilibrium with water. In these cases, constant down-core δ²H values for n-alkyl lipids cannot be reliably used to infer a lack of isotopic exchange.     

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.     

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.     

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.     

The motion of Adria during the Late Jurassic and Cretaceous: New paleomagnetic results from stable Istria

The motion of Adria, the largest lithospheric fragment in the Central Mediterranean region, has played an important role in the tectonic development of the surrounding mountain chains and even of distant areas, like the Eastern Alps or the Pannonian basin. The available paleomagnetic data were insufficient to constrain this motion, except in a general way. In this paper, new paleomagnetic results are presented from one of the stable parts of Adria which emerge from the Adriatic Sea. The results were obtained on weakly magnetic platform carbonates of the mud-supported type, collected from 21 geographically distributed localities.The results, combined with mean paleomagnetic directions from selected localities from a pioneer study in Istria that were chosen using statistical criteria, were divided into three age groups (Tithonian–Aptian, Albian–Cenomanian, Turonian–Coniacian). The paleomagnetic poles calculated for each of them (Tithonian–Aptian): λ(N) = 47°, (E) = 275°, k = 67, α₉₅ = 9.4°, N = 5; Albian-Cenomanian: λ(N) = 58°, (E) = 253°, k = 145, α₉₅ = 4.3°, N = 9; Turonian–Coniacian: λ(N) = 63°, (E) = 261°, k = 50, α₉₅ = 7.3°, N = 9) reveal a moderate shift during the Cretataceous, which is comparable with that calculated from the African reference poles. However, the Istrian apparent polar wander path is slightly displaced from the African curve, as a consequence of about 10° counterclockwise rotation of Istria, with respect to Africa. This rotation angle is more that 10° smaller than the difference measured for the Mid-Late Eocene between the paleomagnetic direction of platform carbonates from Istria and the African reference direction. This difference may be the consequence of a small clockwise rotation of Istria, with respect to Africa, most probably at the end of Cretaceous.     

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.     

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.     

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

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.     

Feldspathic clasts in Yamato-86032: Remnants of the lunar crust with implications for its formation and impact history

Low concentrations of Th and Fe in the Yamato (Y)-86032 bulk meteorite support earlier suggestions that Y-86032 comes from a region of the moon far distant from the Procellarum KREEP Terrain (PKT), probably from the lunar farside. ³⁹Ar–⁴⁰Ar, Rb–Sr, Sm–Nd, and Sm-isotopic studies characterize the chronology of Y-86032 and its precursors in the mega regolith. One of the rock types present in a light gray breccia lithology is an anorthosite characterized by plagioclase with An 93, i.e., more sodic than lunar FANs, but with very low ⁸⁷Rb/⁸⁶Sr and ⁸⁷Sr/⁸⁶Sr similar to those of FANs. (FAN stands for Ferroan Anorthosite). This “An93 anorthosite” has Nd-isotopic systematics similar to those of nearside norites. A FAN-like “An97 anorthosite” is present in a second light-colored feldspathic breccia clast and has a more negative ε_Nd value consistent with residence in a LREE-enriched environment as would be provided by an early plagioclase flotation crust on the Lunar Magma Ocean (LMO). This result contrasts with generally positive values of ε_Nd for Apollo 16 FANs suggesting the possibility of assymetric development of the LMO. Other possible explanations for the dichotomy in ε_Nd values are advanced in the text. The Y-86032 protolith formed at least 4.43 ± 0.03 Ga ago as determined from a Sm–Nd isochron for mineral fragments from the breccia clast composed predominantly of An93 anorthosite and a second clast of more varied composition. We interpret the mineral fragments as being predominatly from a cogenetic rock suite. An ³⁹Ar–⁴⁰Ar age of 4.36–4.41 ± 0.035 Ga for a third clast composed predominantly of An97 anorthosite supports an old age for the protolith. Initial ¹⁴³Nd/¹⁴⁴Nd in that clast was −0.64 ± 0.13 ε-units below ¹⁴³Nd/¹⁴⁴Nd in reservoirs having chondritic Sm/Nd ratios, consistent with prior fractionation of mafic cumulates from the LMO. A maximum in the ³⁹Ar–⁴⁰Ar age spectrum of 4.23 ± 0.03 Ga for a second sample of the same feldspathic breccia clast probably reflects some diffusive ⁴⁰Ar loss. Lack of solar wind and lunar atmosphere implanted Ar in the light gray breccia clast allows determination of an ³⁹Ar/⁴⁰Ar age of 4.10 ± 0.02 Ga, which is interpreted as the time of initial brecciation of this litholgy. After correction for implanted lunar atmosphere ⁴⁰Ar, impact melt and dark regolith clasts give Ar ages of 3.8 ± 0.1 Ga implying melt formation and final breccia assembly 3.8 Ga ago. Some breccia lithologies were exposed to thermal neutron fluences of 2 × 10¹⁵ n/cm², only about 1% of the fluence experienced by some other lunar highlands meteorites. Other lithologies experienced neutron fluences of 1 × 10¹⁵ n/cm². Thus, Y-86032 spent most of the time following final brecciation deeply buried in the megaregolith. The neutron fluence data are consistent with cosmogenic ³⁸Ar_cos cosmic ray exposure ages of 10 Ma. Variations among differing lithologies in the amount of several regolith exposure indicators, including cosmogenic noble gas abundances, neutron capture induced variations in Sm isotopic abundances, and Ir contents, are consistent with a period of early (>3.8 Ga ago) lunar regolith exposure, subsequent deep burial at >5 m depth, and ejection from the moon 7–10 Ma ago.     

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.     

Vertical uplift resistance of circular plate anchors in clays under undrained condition

The vertical uplift resistance of circular plate anchors, embedded horizontally in a clayey stratum whose cohesion increases linearly with depth, has been obtained under undrained ( = 0) condition. The axi-symmetric static limit analysis formulation in combination with finite elements proposed recently by the authors has been employed. The variation of the uplift factor (F_c) with changes in the embedment ratio (H/B) has been computed for several rates of increases of soil cohesion with depth. It is noted that in all the cases, the magnitude of F_c increases continuously with depth up to a certain value of H_cr/B beyond which the uplift factor becomes essentially constant. The proposed static limit analysis formulation is seen to provide acceptable results even for the two other simple chosen axi-symmetric problems.     

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.     

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.     

Diagenetic mineralization in Pennsylvanian coals from Indiana, USA: C/C and O/O implications for cleat origin and coalbed methane generation

Cleats and fractures in southwestern Indiana coal seams are often filled with authigenic kaolinite and/or calcite. Carbon- and oxygen-stable isotope ratios of kaolinite, calcite, and coalbed CO₂ were evaluated in combination with measured values and published estimates of δ¹⁸O of coalbed paleowaters that had been present at the time of mineralization. δ¹⁸O_mineral and δ¹⁸O_water values jointly constrain the paleotemperature of mineralization. The isotopic evidence and the thermal and tectonic history of this part of the Illinois Basin led to the conclusion that maximum burial and heat-sterilization of coal seams approximately 272 Ma ago was followed by advective heat redistribution and concurrent precipitation of kaolinite in cleats at a burial depth of < 1600 m at  78 ± 5 °C. Post-Paleozoic uplift, the development of a second generation of cleats, and subsequent precipitation of calcite occurred at shallower burial depth between  500 to  1300 m at a lower temperature of 43 ± 6 °C. The available paleowater in coalbeds was likely ocean water and/or tropical meteoric water with a δ¹⁸O_water  − 1.25‰ versus VSMOW. Inoculation of coalbeds with methanogenic CO₂-reducing microbes occurred at an even later time, because modern microbially influenced ¹³C-enriched coalbed CO₂ (i.e., the isotopically fractionated residue of microbial CO₂ reduction) is out of isotopic equilibrium with ¹³C-depleted calcite in cleats.     

pH-Dependent iron oxide precipitation in a subterranean estuary

Iron-oxide-coated sediment particles in subterranean estuaries can act as a geochemical barrier (“iron curtain”) for various chemical species in groundwater (e.g. phosphate), thus limiting their discharge to coastal waters. Little is known about the factors controlling this Fe-oxide precipitation. Here, we implement a simple reaction network in a 1D reactive transport model (RTM), to investigate the effect of O₂ and pH gradients along a flow-line in the subterranean estuary of Waquoit Bay (Cape Cod, Massachusetts) on oxidative precipitation of Fe(II) and subsequent PO₄ sorption. Results show that the observed O₂ gradient is not the main factor controlling precipitation and that it is the pH gradient at the mixing zone of freshwater (pH 5.5) and seawater (pH 7.9) near the beach face that causes a  7-fold increase in the rate of oxidative precipitation of Fe(II) at  15 m. Thus, the pH gradient determines the location and magnitude of the observed iron oxide accumulation and the subsequent removal of PO₄ in this subterranean estuary.     

Changes in bulk properties and molecular compositions within New Zealand Coal Band solvent extracts from early diagenetic to catagenetic maturity levels

The bulk properties and bitumen molecular compositions of a rank-series of 38 humic coals from the New Zealand Coal Band (Cretaceous–Cenozoic) have been analysed to investigate early maturation processes affecting coaly organic matter through diagenesis to moderate catagenesis (Rank(S_r) 0.0–11.8, R_o 0.23–0.81%). The samples comprise a relatively restricted range of vitrinite rich coal types formed largely from higher land plant material under relatively oxic conditions, but with a significant contribution from microbial biomass. With increasing rank, total organic carbon contents show a general increase, whereas moisture and asphaltene contents decrease. Bitumen yields also decrease through the stages of diagenesis and early catagenesis (Rank(S_r) < 9, R_o < 0.55%), indicating partial loss of initial bitumen during early maturation. Thermal generation of hydrocarbons begins slowly at Rank(S_r)  5–6 (R_o  0.40%) as indicated by the constant occurrence and gradual increase of isoprenoids (e.g., pristane and phytane) and hopanoids in their more mature αβ configuration. This early phase of catagenesis, not previously recognised in New Zealand coals, is followed at Rank(S_r)  9 (R_o  0.55%) by the main catagenesis phase characterised by a more rapid increase in the generation of hydrocarbons, including total n-alkanes, isoprenoids and αβ-hopanes. Changes in the maturity of New Zealand coals can be traced by the Carbon Preference Index and several hopane maturity parameters, including 22S/(22S + 22R), αβ/(αβ + βα) and ββ/(αβ + βα + ββ).     

Palaeomagnetic results from the Early Permian Copacabana Group, southern Peru: Implication for Pangaea palaeogeography

Samples collected from folded carbonate rocks of the Early Permian Copacabana Group exposed in the Peruvian Subandean Zone have been subjected to detailed palaeomagnetic analysis. Thermal demagnetisation of most samples yield stable high unblocking temperature directions dominantly carried by titanomagnetite minerals. This remanence, identified in 32 samples (43 specimens), is exclusively of reverse polarity consistent with the Permian–Carboniferous Reversal Superchron (PCRS). The overall directions pass the fold test at the 99% confidence level and are considered as being a pre-folding remanence acquired in Early Permian times. The Copacabana Group yields an overall mean direction of D = 166°, I = +49° (α₉₅ = 4.5°, k = 131.5, N = 9 sites) in stratigraphic coordinates and a corresponding palaeosouth pole position situated at λ = 68°S,  = 321°E (A₉₅ = 5.2°, K = 100). Combining this pole with the coeval high quality data from South America, Africa and Australia results in a mean pole for Gondwana situated at λ = 34.4°S,  = 065.6°E (A₉₅ = 4.9°, K = 73.6, N = 13 studies) in African coordinates. This pole position supports a Pangaea B palaeogeography in Early Permian times. In contrast, the combined pole for Gondwana diverges from the coeval Laurasian mean pole when assuming the Pangaea A-type configuration. Poor quality of the Gondwana dataset and inclination shallowing in sediments seem to play no role in the misfit between the Permian–Triassic poles from Gondwana and Laurasia in Pangaea A reconstruction.