Variations of b‐values at the western edge of the Ryukyu Subduction Zone,north‐east Taiwan

Using earthquakes relocated in north‐east Taiwan, we estimated b‐value distribution along a cross‐section located near the Ryukyu slab edge, and four b‐value anomalous areas are evidenced: (1) a high b‐value body lying on top of a low Vp, low Vs and high Vp/Vs sausage‐like body was considered as a region of enhanced partial melt or water supply above which seismicity occurs; (2) beneath the Ilan Plain, an anomalous area characterized by b‐values slightly higher than 1.1 might give evidence to the magma conduits to the Kueishantao Island; (3) above the Ryukyu Wadati‐Benioff zone, at depths ranging from 90 to 110 km, a high b‐value anomaly might correspond to the depth where dehydration occurs in the subducting oceanic plate; and (4) a low b‐value area located within the Ryukyu slab, at depths ranging from 70 to 90 km, might be linked to the compressive mechanisms shown by focal mechanisms and the bending of the subducting plate.     

Inferring autogenically induced depositional discontinuities from observations on experimental deltaic shoreline trajectories

Palaeo shoreline is a commonly used proxy for palaeo sea level, but only if deposition is continuous and constant will shoreline trajectory T(l) completely capture sea‐level time‐series E(t). Artificial deltas were generated in the Eurotank flume facility under stepwise tectonic subsidence, periodic sea‐level fluctuation and two periodic water‐discharge scenarios, one in‐phase and the other out‐of‐phase with sea level. Independent input variables tectonic subsidence Y, sea level E and water discharge Q (controlling sediment supply S) were varied and dependent output variable shoreline trajectory T was monitored. These experiments confirm that deposition is discontinuous even for continuous sediment supply, and this hinders the inference of sea‐level curve from shoreline trajectory. These results justify the here‐developed methodology for converting shoreline trajectory from the space domain to the time domain, thereby improving the accuracy of the inferred sea‐level curve.     

An (in‐)coherent metamorphic evolution of high‐P eclogites and their host rocks in the Chinese southwest Tianshan?

In the (ultra‐)high‐P–low‐T metamorphic terrane of the Chinese South Tianshan, discontinuous mafic blocks and boudins (former upper oceanic crust) are now embedded in voluminous (mainly metasedimentary) host rocks. Two different models were proposed and relate the occurrence of both high‐P and ultra‐high‐P mafic and metasedimentary rocks to either (i) a tectonic mélange style exhumation, with no exhumation of coherent units, but different lithologies derived from different depths juxtaposed and intermingled during exhumation in the subduction channel, or (ii) the evolution of two coherent metamorphic belts: one with high‐P and the other with ultra‐high‐P conditions. In contrast to most previous studies in the Chinese South Tianshan which focused either on single eclogites or metasedimentary rocks (assumed as representative), this study concentrates on the systematic investigation of both mafic boudins and their immediate sedimentary host rocks, because the investigation of both lithologies and the comparison of their metamorphic evolution is crucial to reconstruct the geodynamical context of the whole (ultra‐)high‐P–low‐T metamorphic complex. Several sample pairs consisting of both lithologies were geochemically investigated and their respective metamorphic evolution was reconstructed using geothermobarometry and thermodynamic modelling. The latter approach considers changes in the mineral assemblage during the metamorphic evolution, as well as changes in mineral composition, which may help to determine the metamorphic history of a rock despite the preservation of critical mineral assemblages. All samples experienced a clockwise P–T path with overall maximum P–T conditions of 540–550 °C and 1.9–2.25 GPa for the host rocks, and 555–575 °C and 2.2–2.5 GPa for the eclogites. Peak‐metamorphic temperatures of ~525–540 °C of the metasedimentary host rocks were also confirmed by Raman spectroscopy of carbonaceous material. Results from thermobarometry and thermodynamic modelling are consistent with the observation that none of the samples contains mineral relicts indicating UHP conditions (like coesite in garnet) and neither conventional thermobarometry, nor thermodynamical modelling resulted in P–T conditions in the stability field of coesite. Thus, no evidence of ultra‐high‐P conditions was found. Given that the whole sampled river valley lies within the proposed ‘ultra‐high‐P sub‐belt’ and considering former studies, which showed that at other places within this ‘unit’ both ultra‐high‐P and high‐P rocks are now juxtaposed on a small scale, the formation of the whole (ultra‐)high‐P–low‐T metamorphic belt in the Chinese South Tianshan as a tectonic mélange style exhumation is more convincing than the formation and juxtaposition of two coherent metamorphic units with high‐P and ultra‐high‐P conditions respectively.     

Calculated phase equilibria in K2O-FeO-MgO-Al2O3-SiO2-H2O for silica-undersaturated sapphirine-bearing mineral assemblages

Silica‐undersaturated, sapphirine‐bearing granulites occur in a large number of localities worldwide. Such rocks have historically been under‐utilized for estimating P–T evolution histories because of limited experimental work, and a consequent poor understanding of the topology and P–T location of silica‐undersaturated mineral equilibria. Here, a calculated P–T projection for sapphirine‐bearing, silica‐undersaturated metapelitic rock compositions is constructed using THERMOCALC for the FeO‐MgO‐Al₂O₃‐SiO₂ (FMAS) and KFMASH (+K₂O + H₂O) chemical systems, allowing quantitative analysis of silica‐undersaturated mineral assemblages. This study builds on that for KFMASH sapphirine + quartz equilibria [Kelsey et al. (2004) Journal of Metamorphic Geology, vol. 22, pp. 559–578]. FMAS equilibria are significantly displaced in P–T space from silicate melt‐bearing KFMASH equilibria. The large number of univariant silica‐undersaturated KFMASH equilibria result in a P–T projection that is topologically more complex than could be established on the basis of experiments and/or natural assemblages. Coexisting sapphirine and silicate melt (with or without corundum) occur down to c. 900 °C in KFMASH, some 100 °C lower than in silica‐saturated compositions, and from pressures of c.≤1 to ≥12 kbar. Mineral compositions and composition ranges for the calculated phases are consistent with natural examples. Bulk silica has a significant effect on the stability of sapphirine‐bearing assemblages at a given P–T, resulting in a wide variety of possible granulite facies assemblages in silica‐undersaturated metapelites. Calculated pseudosections are able to reproduce many naturally occurring silica‐undersaturated assemblages, either within a single assemblage field or as the product of a P–T trajectory crossing several fields. With an understanding of the importance of bulk composition on sapphirine stability and textural development, silica‐undersaturated assemblages may be utilized in a quantitative manner in the detailed metamorphic investigation of high‐grade terranes.     

Evolution of a crustal-scale transpressive shear zone in the Albany–Fraser Orogen, SW Australia: 2. Tectonic history of the Coramup Gneiss and a kinematic framework for Mesoproterozoic collision of the West Australian and Mawson cratons

Within the Albany–Fraser Orogen of southwestern Australia, the Coramup Gneiss is a NE–SW trending zone of high‐strain rocks that preserves a detailed record of orogenesis related to Mesoproterozoic convergence of the West Australian and Mawson cratons. New structural, metamorphic and U–Pb SHRIMP zircon age data establish that the Coramup Gneiss underwent high‐grade tectonism during both Stage I (c. 1290 Ma) and Stage II (c. 1170 Ma) of the Albany–Fraser Orogeny. Stage I commenced with c. 1300 Ma high‐T, low‐P M1a metamorphism during extension, and the formation of small‐scale ptygmatic folds within a subhorizontal S1a gneissosity. High‐P M1b metamorphism at c. 1290 Ma was accompanied by the transposition and shearing of S1a into a composite, shallow SE‐dipping S1b foliation, and the development of tight recumbent F1b folds with S1‐parallel axial surfaces and asymmetries indicating NW‐directed thrusting. The preservation of a similar P–T–time record in the Fraser Complex (NE of the Coramup Gneiss) is consistent with large‐scale, NW‐directed Stage I thrusting of the Mawson Craton margin over the south‐eastern edge of the West Australian Craton. Stage II tectonism in the western Coramup Gneiss involved high‐T, low‐P M2a metamorphism and the formation of subvertical SE‐dipping D2 shear zones, shallow SW‐plunging L2 mineral stretching lineations, and NW‐verging F2 folds with S2‐parallel axial surfaces. A synkinematic pegmatite dyke emplaced into a D2 shear zone yielded a U–Pb SHRIMP zircon age of 1168 ± 12 Ma. Kinematic indicators suggest a combination of pure shear flattening perpendicular to S2, and dextral simple shear. However, contemporaneous structures elsewhere in the Albany–Fraser Orogen are consistent with continued NW–SE convergence at craton‐scale during Stage II, and oblique compression in the Coramup Gneiss is attributed to the arcuate geometry of the orogen‐scale deformation front.     

P–T–t path of the Hercynian low-pressure rocks from the Mandatoriccio complex (Sila Massif, Calabria, Italy): new insights for crustal evolution

The tectono‐metamorphic evolution of the Hercynian intermediate–upper crust outcropping in eastern Sila (Calabria, Italy) has been reconstructed, integrating microstructural analysis, P–T pseudosections, mineral isopleths and geochronological data. The studied rocks belong to a nearly complete crustal section that comprises granulite facies metamorphic rocks at the base and granitoids in the intermediate levels. Clockwise P–T paths have been constrained for metapelites of the basal level of the intermediate–upper crust (Umbriatico area). These rocks show noticeable porphyroblastic textures documenting the progressive change from medium‐P metamorphic assemblages (garnet‐ and staurolite‐bearing assemblages) towards low‐P/high‐T metamorphic assemblages (fibrolite‐ and cordierite‐bearing assemblages). Peak‐metamorphic conditions of ～590 °C and 0.35 GPa are estimated by integrating microstructural observations with P–T pseudosections calculated for bulk‐rock and reaction‐domain compositions. The top level of the intermediate–upper crust (Campana area) recorded only the major heating phase at low‐P (～550 °C and 0.25 GPa), as documented by the static growth of biotite spots and of cordierite and andalusite porphyroblasts in metapelites. In situ U–Th–Pb dating of monazite from schists containing low‐P/high‐T metamorphic assemblages gave a weighted mean U–Pb concordia age of 299 ± 3 Ma, which has been interpreted as the timing of peak metamorphism. In the framework of the whole Hercynian crustal section the peak of low‐P/high‐T metamorphism in the intermediate‐to‐upper crust took place concurrently with granulite facies metamorphism in the lower crust and with emplacement of the granitoids in the intermediate levels. In addition, decompression is a distinctive trait of the P–T evolution both in the lower and upper crust. It is proposed that post–collisional extension, together with exhumation, is the most suitable tectonic setting in which magmatic and metamorphic processes can be active simultaneously in different levels of the continental crust.     

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

A comparison of variational formats for porous media subjected to dynamic loading

This paper presents a comparison of two variational formats for fully saturated porous media subjected to dynamic loading, whereby the general situation of relative fluid acceleration is considered: (1) the classical three‐field ( u , p, w )‐format and (2) a novel two‐field ( u , p)‐format, where the seepage velocity w is a spatially ‘local’ field whose treatment resembles that of internal variables in material models. The limited numerical comparison shows that the ( u , p)‐format competes well with the ( u , p, w )‐format. Indeed, it is consistent with the general acceleration modeling in the full range of permeabilities. Moreover, in the low permeability regime (where the magnitude of w is insignificant), the new format reflects the situation pertinent to ‘added‐mass’ and is more efficient than the classical ( u , p, w )‐format. Finally, the ( u , p)‐formatcan conveniently be implemented in existing FE‐codes based on the ‘added mass’ formulation. Copyright © 2010 John Wiley & Sons, Ltd.     

Latitudinally different responses of Tethyan shoal‐water carbonate systems to the Early Aptian oceanic anoxic event (OAE 1a)

Ocean‐wide anoxic events represent intensively investigated anomalies in the global carbon cycle. Most previous research has focussed on hemipelagic and pelagic settings and on the relationship between black‐shale deposition and carbon‐isotope excursions. The study of ocean‐wide anoxic events and coeval shallow‐water settings is now increasingly seen as an interesting complementary approach, but one that is not without problems. Whereas platform drowning characterizes the Early Aptian of the northern Tethyan margin, Lithocodium–Bacinella‐rich facies and ongoing shoal‐water sedimentation at the southern Tethyan margin (Oman) bears important information on potential causes of carbon‐cycle perturbations. The present paper seeks to test the supra‐regional relevance of the Oman data by investigating coeval central Tethyan limestones. Three Lower Aptian shoal‐water sections in Istria (Croatia), deposited on the isolated Adriatic Carbonate Platform, are investigated applying chemostratigraphy (carbon and strontium) and detailed sedimentological analysis. The focus is on peritidal to lagoonal facies characterized by mass occurrences of Lithocodium–Bacinella, an enigmatic microencruster community. Lithocodium–Bacinella facies occurs predominantly in layers ranging from one to several centimetres in stratigraphic thickness, with several layers merging to metre‐thick packages. Growth fabrics within the layers include oncoidal morphotypes, lumps, interconnected patches and columns, layers and rare nodular to massive bindstone facies. These growth patterns show a remarkable regional extent and consistency over study sites distributed several kilometres apart. This widespread distribution suggests that specific Lithocodium–Bacinella morphotypes might serve as regional stratigraphic markers. The high‐resolution carbon‐isotope chemostratigraphy presented here is based on pristine rudist shells and matrix micrite samples and calibrated against strontium‐isotope data obtained from screened rudist low‐Mg calcite. The chemostratigraphic data are consistent with existing biostratigraphic data and place the studied strata at the onset of Early Aptian oceanic anoxic event 1a. Moreover, results indicate the near‐coeval nature of Lithocodium–Bacinella bloom facies in Istria and Oman. The outcomes of this study point to latitudinally different responses of Tethyan shoal‐water carbonate systems (platform drowning versus Lithocodium–Bacinella blooms) to the ocean‐wide anoxic event 1a.     

High‐precision Determination of Gold Mass Fractions in Geological Reference Materials by Internal Standardisation

Determination of gold abundances in natural rock is critical for applications, but very challenging. Here, we report a method for determining gold with a very low mass fraction (> 0.01 ng g^?1) in rocks. The method involves Carius tube digestion with reverse aqua regia, chromatographic separation to remove most of the sample matrix and measurement by high‐sensitivity ICP‐MS. The mono‐isotopic element gold was quantified by external calibration using an internal standardisation of gold to platinum that was precisely determined by isotope dilution. The method is robust and the obtained results are indistinguishable (< 5–10%, 2s) from those independently obtained by a standard addition technique on the same solution. The results from reference materials TDB‐1 and GPt‐2 are consistent with the certified values and those determined by HF‐aqua regia digestion, confirming the validity of the method. TDB‐1 (n = 20), GPt‐2 (n = 6), BHVO‐2 (n = 9) and other mafic RMs are homogenous for gold (10–20%, 2s) at the 2 g test portion level; however, sample heterogeneity affects some RMs. Gold and platinum‐group elements also display different extents of sample heterogeneity for different RMs. Given the homogeneity observed for TDB‐1, GPt‐2 and BHOV‐2, they are recommended as well‐suited RMs for inter‐laboratory comparison studies of gold.     

Anatectic record and contrasting P–T paths of aluminous gneisses from the central Grenville Province

Anatectic aluminous gneisses, some derived from sedimentary rocks of broadly pelitic composition and others from hydrothermally altered felsic volcanic rocks, are exposed in the mid‐P and high‐P segments of the hinterland in the central Grenville Province. These gneisses consist dominantly of garnet, biotite, K‐feldspar, plagioclase and quartz, with sillimanite or kyanite, and display microstructural evidence of anatexis by fluid‐absent reactions consuming muscovite and/or biotite. Melt‐related microstructures, such as inter‐granular films and/or interstitial quartz or feldspar enclosing relict phases, are most abundant in the metasedimentary samples. Despite anatexis at granulite facies conditions, the hydrothermally altered rocks preserve earlier features attributed to the circulation of hydrothermal fluids, such as sillimanite seams, dismembered quartz veins and garnet‐rich aluminous nodules in a K‐feldspar‐dominated matrix. Microstructural and mineral chemical data, integrated with P–T pseudosections calculated with thermocalc for the metasedimentary rocks, permit qualitative constraints on the P–T paths. Data from a high‐P kyanite‐bearing sample are consistent with a steep prograde P–T path up to ~14.5 kbar and 860–900 °C, followed by decompression with minor cooling to the solidus at ~11 kbar and 870 °C. This pressure‐dominated P–T path is similar to those inferred in other parts of the high‐P segment in the central Grenville Province. In contrast, the P–T path predicted from a mid‐P sillimanite‐bearing paragneiss has a strong temperature gradient with P–T of ~9.5 kbar and 850 °C at the thermal peak, and a retrograde portion down to ~8 kbar and 820 °C. In a broad sense, these two contrasting P–T patterns are consistent with predictions of thermo‐mechanical modelling of large hot orogens in which P–T paths with strong pressure gradients exhume deeper rocks in the orogenic flanks, whereas P–T paths with strong temperature gradients in the orogenic core reflect protracted lateral transport of ductile crust beneath a plateau.     

Fault‐bend folding as an end‐member solution of (double‐edge) fault‐propagation folding

Fault‐bend folding is the most commonly used kinematic mechanism to interpret the architecture and evolution of thrust‐related anticlines in thrust wedges. However, its basic requirement of an instantaneous propagation of the entire fault before hangingwall deformation, limits its kinematic effectiveness. To overcome this limitation, we used the interdependence between fold shape and fault slip vs. propagation rate (S/P ratio) implemented in double‐edge fault‐propagation folding. We show that very small S/P values produce fault‐propagation anticlines that, when transported forelandward along an upper décollement layer, closely resemble fault‐bend anticlines. Accordingly, if small geometric discrepancies between the two solutions are accepted, transported double‐edge fault‐propagation provides an effective kinematic alternative to fault‐bend folding. Even at very low S/P values, it in fact predicts a fast but finite propagation rate of the fault. We thus propose that double‐edge fault‐propagation folding provides a broadly applicable model of fault‐related folding that includes fault‐bend folding as an end‐member kinematic solution. Terra Nova, 18, 270–275, 2006     

Evidence from Fe- and Al-rich metapelites for thrust loading in the Transangarian region of the Yenisey Ridge, eastern Siberia

In the Transangarian region of the Yenisey Ridge in eastern Siberia (Russia), Fe‐ and Al‐rich metapelitic schists of the Korda plate show field and petrological evidence of two superimposed metamorphic events. An early middle Proterozoic event with age of c.1100 Ma produced low‐pressure, andalusite‐bearing assemblages at c. 3.5–4 kbar and 540–560 °C. During a subsequent late Proterozoic event at c. 850 Ma, a medium‐pressure, regional metamorphic overprint produced kyanite‐bearing mineral assemblages that replaced minerals formed in the low‐pressure event. Based on the results of geothermobarometry and P–T path calculations it can be shown that pressure increased from 4.5 to 6.7 kbar at a relatively constant temperature of 540–600 °C towards a major suture zone called the Panimba thrust. In order to produce such nearly isothermal loading of 1–7 °C km ^?1, we propose a model for the tectono‐metamorphic evolution of the study area based on crustal thickening caused by south‐westward thrusting of the 5–7 km‐thick upper‐plate metacarbonates over lower‐plate metapelites with velocity of c. 350 m Myr^?1. A small temperature increase (up to 20 ± 15 °C) of the upper part of the overlapped plate is explained by specific behaviour of steady‐state geotherms calculated using lower radioactive heat production of metacarbonates as compared with metapelites. The suggested thermal‐mechanical model corresponds well with P–T paths inferred from obtained thermobarometric data and correlates satisfactorily with P–T trajectories predicted by other two‐dimensional thermal models for different crustal thickening and exhumation histories.     

Normal Alkane Distributions in Soil Samples along a Lhasa-Bharatpur Transect

We present n-alkane distributions in the soil samples along a transect from Lhasa on the Tibetan Plateau to Bharatpur in Nepal, which covers a large geographical area and a wide range of climatic conditions. These data allow us to assess the significance of n-alkane distributions in different vegetation types and their relationships to temperature and precipitation. In the tropical rainforest and broadleaved forest zones, n-alkanes exhibit a bimodal distribution pattern with dominant homologues around n-C_(31) and n-C_(23). The bimodal distribution of tropical rainforest n-alkanes may be a general pattern because of the presence of many lianas, epiphytic plants, algae and the strong microbial activity and degradation involved in the postdepositional process. In the warm-temperate mixed forest and needle-leaved forest zones, the long-chain alkanes have a pattern of n-C_(31) n-C_(33) n-C_(27). In the alpine shrub and grassland zone, although the most abundant homologue is n-C_(31) , relatively high n-C_(23) concentrations have also been observed in some samples. The statistical results show a good correlation between n-alkane proxies and climatic factors. The average chain length(ACL) values are positively correlated with precipitation and temperature across the transect. The carbon preference index(CPI) values show a negative correlation with temperature and precipitation along a Lhasa–Bharatpur transect. Although there is a positive relationship between ACL and temperature, it is difficult to separate the two climatic variables(temperature and precipitation) because they are well coupled in the monsoon region.     

Self‐consistent subduction initiation induced by mantle flow

Mantle circulation in planets with strongly temperature‐dependent viscosity results in stagnant‐lid convection. It is fundamental to understand how this stagnant‐lid regime can change into a plate‐like convection regime as on the present‐day Earth. Here, we use 2D numerical models to study subduction initiation from an initial stagnant lid with laboratory‐consistent parameters and without pre‐existing weak zones or kinematic boundary conditions. Our results show that subduction can be initiated dynamically as a result of a thermal localization instability. The lithosphere may deform in a stagnant‐lid mode, an un‐necking mode, a symmetric‐subduction mode or an asymmetric‐subduction mode. The asymmetric‐subduction mode occurs only for relatively large friction angles and moderate thermal ages, and the presence of heterogeneities increases the parameter space of this mode. The limited parameter space might explain why only the present‐day Earth has plate tectonics, and suggests that the initiation of plate tectonics is more difficult than previously anticipated.     

Lipid Biomarkers and Their Stable Carbon Isotopes inAncient Seep Carbonates from SW Taiwan,China

Four massive brecciated, chimney-like, and slender pipe network carbonate samples(JA-4, JA-5, JX-8 and BG-12) were collected from southwestern Taiwan, which were suggested to have formed as a result of anaerobic oxidization of methane(AOM). Considering that the environmental conditions of the carbonates precipitation and the sources of carbon and organic matter need to be further declared, molecular fossils and compound-specific carbon isotopic investigations of the carbonates were conducted in this study. According to lipid biomarkers of 2,6,10,15,19-pentamethyleicosane(PMI) and squalane diagnostic to methanotrophic archaea, as well as the extremely low δ¹³C values(as low as -1¹³.4‰) detected in samples JA-4, JA-5 and JX-8, these carbonates were revealed to be a result of AOM. Based on the varied δ¹³C values of characteristic archaea biomarkers in specific samples, biogenic methane was proposed to be responsible for the formation of samples JA-4 and JA-5, whereas a mixed carbon source of ¹³C-depleted methane and ¹³C-enriched residual CO_2 from methanogenesis was suggested for the carbonate of JX-8 due to the co-occurrence of a highly positive δ¹³ C_carb value(+8‰) and a moderate ¹³C depletion of PMI. The low content of AOM-related biomarkers and the absence of indicators for ANME-2 suggested that these carbonates were formed in weak seep settings. By comparison, no typical lipid biomarkers for methanotrophic archaea was detected in carbonate BG-12. The short-chain and long-chain n-alkanes accounted for 30% and 45% of all hydrocarbons, respectively, with a CPI value of 1.2, suggesting that the n-alkanes were derived from both marine organisms and terrestrial inputs. A low thermal maturity could be revealed by the incomplete equilibrium value of the C³¹αβ 22S/(22S+22R) ratio(0.5), and the carbonate BG-12 was probably deposited in a suboxic condition indicated by a value of Pr/Ph ratio(2.5).     

Palaeoproterozoic high-T, low-P metamorphism and dehydration melting in metapelites from the Mopunga Range, Arunta Inlier, central Australia

A sequence of psammitic and pelitic metasedimentary rocks from the Mopunga Range region of the Arunta Inlier, central Australia, preserves evidence for unusually low pressure (c. 3 kbar), regional‐scale, upper amphibolite and granulite facies metamorphism and partial melting. Upper amphibolite facies metapelites of the Cackleberry Metamorphics are characterised by cordierite‐andalusite‐K‐feldspar assemblages and cordierite‐bearing leucosomes with biotite‐andalusite selvages, reflecting P–T conditions of c. 3 kbar and c. 650–680 °C. Late development of a sillimanite fabric is interpreted to reflect either an anticlockwise P–T evolution, or a later independent higher‐P thermal event. Coexistence of andalusite with sillimanite in these rocks appears to reflect the sluggish kinematics of the Al₂SiO₅ polymorphic inversion. In the Deep Bore Metamorphics, 20 km to the east, dehydration melting reactions in granulite facies metapelites have produced migmatites with quartz‐absent sillimanite‐spinel‐cordierite melanosomes, whilst in semipelitic migmatites, discontinuous leucosomes enclose cordierite‐spinel intergrowths. Metapsammitic rocks are not migmatised, and contain garnet–orthopyroxene–cordierite–biotite–quartz assemblages. Reaction textures in the Deep Bore Metamorphics are consistent with a near‐isobaric heating‐cooling path, with peak metamorphism occurring at 2.6–4.0 kbar and c. 750–800 °C. SHRIMP U–Pb dating of metamorphic zircon rims in a cordierite‐orthopyroxene migmatite from the Deep Bore Metamorphics yielded an age of 1730 ± 7 Ma, whilst detrital zircon cores define a homogeneous population at 1805 ± 7 Ma. The 1730 Ma age is interpreted to reflect the timing of high‐T, low‐P metamorphism, synchronous with the regional Late Strangways Event, whereas the 1805 Ma age provides a maximum age of deposition for the sedimentary precursor. The Mopunga Range region forms part of a more extensive low‐pressure metamorphic terrane in which lateral temperature gradients are likely to have been induced by localised advection of heat by granitic and mafic intrusions. The near‐isobaric Palaeoproterozoic P–T–t evolution of the Mopunga Range region is consistent with a relatively transient thermal event, due to advective processes that occurred synchronous with the regional Late Strangways tectonothermal event.     

The δ60/58Ni Values of Twenty‐Six Selected Geological Reference Materials

The high‐precision δ^60/58Ni values of twenty‐six geological reference materials, including igneous rocks, sedimentary rocks, stream sediments, soils and plants are reported. The δ^60/58Ni values of all samples were determined by double‐spike MC‐ICP‐MS (Nu Plasma III). Isotope standard solution (NIST SRM 986) and geological reference materials (BHVO‐2, BCR‐2, JP‐1, PCC‐1, etc.) were used to evaluate the measurement bias and intermediate precision over a period of six months. Our results show that the intermediate precision of Ni isotope determination was 0.05‰ (2s, n = 69) for spiked NIST SRM 986 and typically 0.06‰ for actual samples, and the δ^60/58Ni _{NIST SRM 986} values were in excellent agreement with previous studies. Eighteen high‐precision Ni isotope ratios of geological reference materials are first reported here, and their δ^60/58Ni values varied from ?0.27‰ to 0.52‰, with a mean of 0.13 ± 0.34‰ (2s, n = 18). Additionally, SGR‐1b (0.56 ± 0.04‰, 2s), GSS‐1 (?0.27 ± 0.06‰, 2s), GSS‐7 (?0.11 ± 0.01‰, 2s), GSD‐10 (0.46 ± 0.06‰, 2s) and GSB‐12 (0.52 ± 0.06‰, 2s) could potentially serve as candidate reference materials for Ni isotope fractionation and comparison of Ni isotopic compositions among different laboratories.     

An example of a highly bioturbated,storm‐influenced shoreface deposit: Upper Jurassic Ula Formation,Norwegian North Sea

Integrated ichnological and sedimentological analyses of core samples from the Upper Jurassic Ula Formation in the Norwegian Central Graben were undertaken to quantify the influence of storm waves on sedimentation. Two main facies associations (offshore and shoreface) that form a progradational coarsening upward succession are recognizable within the cores. The offshore deposits are characterized by massive to finely laminated mudstones and fine‐grained sandstones, within a moderately to highly bioturbated complex. The trace fossil assemblage is dominated by deposit‐feeding structures (for example, Planolites, Phycosiphon and Rosselia) and constitutes an expression of the proximal Zoophycos to distal Cruziana ichnofacies. The absence of grazing behaviours and dominance of deposit‐feeding ichnofossils is a reflection of the increased wave energies present (i.e. storm‐generated currents) within an offshore setting. The shoreface succession is represented by highly bioturbated fine‐grained to medium‐grained sandstones, with intervals of planar and trough cross‐bedding, thin pebble lags and bivalve‐rich shell layers. The ichnofossil assemblage, forming part of the Skolithos ichnofacies, is dominated by higher energy Ophiomorpha nodosa ichnofossils and lower energy Ophiomorpha irregulaire and Siphonichnus ichnofossils. The presence of sporadic wave‐generated sedimentary structures and variability in ichnofossil diversity and abundance attests to the influence of storm‐generated currents during deposition. As a whole, the Ula Formation strongly reflects the influence of storm deposits on sediment deposition; consequently, storm‐influenced shoreface most accurately describes these depositional environments.     

Effect of pedogenesis on the stable isotopic composition of calcretes and n‐alkanes: Implications for palaeoenvironmental reconstruction

In a fluvial system, depending on sub‐aerial exposure, non‐pedogenic pond calcretes can be modified into pedogenic calcretes. The present study attempts to understand the effect of sub‐aerial exposure and pedogenesis on calcretes using carbon and oxygen isotopic composition. For this purpose, two profiles (profile‐A and profile‐B) from the same stratigraphic level in Rayka from the western part of India were selected. The profiles are separated by a distance of 500 m and showed differences in calcrete characteristics. In profile‐A, the calcretes showed pedogenic features (root traces and void filling spar) whereas calcretes in profile‐B showed non‐pedogenic characteristics (fine laminations). However, some of the calcretes in profile‐A exhibited remnants of fine laminations suggesting that initially the calcretes had a non‐pedogenic origin but were modified due to pedogenesis. In profile‐A, the carbon and oxygen isotope values of pedogenic calcrete (δ¹³C_PC and δ¹⁸O_PC) showed more variation compared with non‐pedogenic pond calcretes (δ¹³C_SPC and δ¹⁸O_SPC) in profile‐B. The δ¹³C_PC and δ¹³C_SPC values exhibited a spread of 3·0‰ and 1·3‰, respectively, and δ¹⁸O_PC and δ¹⁸O_SPC values showed a spread of 2·3‰ and 1·3‰, respectively. The differences in the isotopic composition between the two profiles suggest that pedogenesis controlled the isotopic inheritance in calcretes. In addition, the carbon isotopic composition of organic matter (δ¹³C_OM) and n‐alkanes (δ¹³C_n‐alk) that forms the basis of palaeovegetational reconstruction have also been measured to understand the effect of pedogenesis on organic matter in both of the profiles. The average δ¹³C_OM values in profile‐A and profile‐B are ?23·4‰ and ?21·1‰, respectively. The disparity in δ¹³C_OM values is a result of the difference in the sources and preservation of organic matter. However, the δ¹³C_n‐alk values show a similar trend in profile‐A and profile‐B, indicating that sources of n‐alkanes are the same in both of the profiles and δ¹³C_n‐alk values are unaffected by the pedogenic modifications.