Refining the pre-industrial atmospheric Pb isotope evolution curve in Europe using an 8000 year old peat core from NW Spain

Pb pollution has existed for several millennia and remains relevant today. By using peat cores as environmental archives it is possible to reconstruct this long history on a regional scale. This is a significant contribution to the findings from ice core records, the only other archive recording purely atmospheric additions. Without information that allows linking and comparison between sites regionally, within Europe and elsewhere, our ability to make coherent global models of the natural Pb cycle, and anthropogenic forcing of this cycle, is limited. In this respect, the characteristics of the Pb pre-pollution aerosol (PPA) are important to define globally. We characterize for the first time a PPA in Southern Europe with [Pb] = 0.78 ± 0.86 μg g^− 1, net Pb accumulation rates of 0.032 ± 0.030 mg m^− 2 y^− 1 and a ²⁰⁶Pb / ²⁰⁷Pb signature of 1.25470 ± 0.02575. This PPA Pb isotope signature is more radiogenic than that found thus far in Western and Northern Europe. Spain is a historically important mining site. Using three-isotope plots and a pool of potential Pb isotope signatures, a detailed source appointment of both natural and anthropogenic Pb sources was made. We found evidence of Saharan aridification and its termination ∼4400 BP and/or agricultural signals and strong local control (from rock and soil) of the Pb PPA. Human impact is first recorded at 3210 BP but does not exceed 50% of deposited Pb until 3005 BP. Mines in SE Spain dominate early Pb pollution history at this site. During the rise of Roman rule, contributions come from mines in N, NW and SW Spain with no strong indication of other European mining activities. In Medieval and Industrial times local contributions to the peat bog are reduced.     

Provenance and Magmatic-Metamorphic Evolution of a Variscan Island-Arc Complex: Constraints from U-Pb Dating, Petrology, and Geospeedometry of the Kyffhauser Crystalline Complex, Central Germany

The Kyffhäuser Crystalline Complex, Central Germany, formspart of the Mid-German Crystalline Rise, which is assumed torepresent the Variscan collision zone between the East Avalonianterrane and the Armorican terrane assemblage. High-precisionU–Pb zircon and monazite dating indicates that sedimentaryrocks of the Kyffhäuser Crystalline Complex are youngerthan c. 470 Ma and were intruded by gabbros and diorites between345 ± 4 and 340 ± 1 Ma. These intrusions had magmatictemperatures between 850 and 900°C, and caused a contactmetamorphic overprint of the sediments at P–T conditionsof 690–750°C and 5–7 kbar, corresponding toan intrusion depth of 19–25 km. At 337 ± 1 Ma themagmatic–metamorphic suite was intruded by granites, syenitesand diorites at a shallow crustal level of some 7–11 km.This is inferred from a diorite, and conforms to P–T pathsobtained from the metasediments, indicating a nearly isothermaldecompression from 5–7 to 2–4 kbar at 690–750°C.Subsequently, the metamorphic–magmatic sequence underwentaccelerated cooling to below 400°C, as constrained by garnetgeospeedometry and a previously published K–Ar muscoviteage of 333 ± 7 Ma. With respect to P–T–D–tdata from surrounding units, rapid exhumation of the KCC canbe interpreted to result from NW-directed crustal shorteningduring the Viséan. KEY WORDS: contact metamorphism; U–Pb dating; hornblende; garnet; Mid-German Crystalline Rise; P–T pseudosection     

Assimilation and Fractional Crystallization Controlled by Transport Process of Crustal Melt: Implications from an Alkali Basalt-Dacite Suite from Rishiri Volcano, Japan

Mechanisms of fractional crystallization with simultaneous crustalassimilation (AFC) are examined for the Kutsugata and Tanetomilavas, an alkali basalt–dacite suite erupted sequentiallyfrom Rishiri Volcano, northern Japan. The major element variationswithin the suite can be explained by boundary layer fractionation;that is, mixing of a magma in the main part of the magma bodywith a fractionated interstitial melt transported from the mushyboundary layer at the floor. Systematic variations in SiO₂ correlatewith variations in the Pb, Sr and Nd isotopic compositions ofthe lavas. The geochemical variations of the lavas are explainedby a constant and relatively low ratio of assimilated mass tocrystallized mass (‘r value’). In the magma chamberin which the Kutsugata and Tanetomi magmas evolved, a strongthermal gradient was present and it is suggested that the marginalpart of the reservoir was completely solidified. The assimilantwas transported by crack flow from the partially fused floorcrust to the partially crystallized floor mush zone throughfractures in the solidified margin, formed mainly by thermalstresses resulting from cooling of the solidified margin andheating of the crust. The crustal melt was then mixed with thefractionated interstitial melt in the mushy zone, and the mixedmelt was further transported by compositional convection tothe main magma, causing its geochemical evolution to be characteristicof AFC. The volume flux of the assimilant from the crust tothe magma chamber is suggested to have decreased progressivelywith time (proportional to t^–1/2), and was about 3 x 10^–2m/year at t = 10 years and 1 x 10^–2 m/year at t = 100years. It has been commonly considered that the heat balancebetween magmas and the surrounding crust controls the couplingof assimilation and fractional crystallization processes (i.e.absolute value of r). However, it is inferred from this studythat the ratio of assimilated mass to crystallized mass canbe controlled by the transport process of the assimilant fromthe crust to magma chambers. KEY WORDS: assimilation and fractional crystallization; mass balance model; magma chamber; melt transport; Pb isotope     

江苏王港盐沼的现代沉积速率

江苏海岸湿地总面积超过5000km^2，但围垦活动的加剧使盐沼湿地面积迅速减少；为了达到海岸防护、保滩促淤的目的，先后引种了大米草和互花米草。通过采集江苏王港潮滩盐沼湿地9处柱状样和挖取探槽剖面，以及对沉积剖面样品进行^210Pb和^137Cs测年分析，探讨了^210Pb的本底及影响^210Pb富集的因素，估算了王港盐沼湿地的现代沉积速率。王港潮滩盐沼沉积物对^210Pb的吸附作用较弱；受风暴潮、生物扰动、物源变化及实验误差等因素的影响，在^210Pbex剖面上出现了数据异常点，将这些数据点剔除后计算得到的沉积速率为33cm yr^-1，^137Cs测年显示，该地区1963年以来的平均沉积速率3．1cm yr^-1，与^210Pb法及前人研究结果相一致。根据一个典型剖面的^137Cs测年数据分析，王港潮滩盐沼的沉积过程根据地貌特征的不同可分3个阶段，大米草覆盖阶段滩面高程迅速增加，互花米草的生长提高了滩面淤积速率。     

Zircon U–Pb ages and Hf isotope compositions of migmatites from the North Qinling terrane and their geological implications

Migmatites are predominant in the North Qinling (NQ) orogen, but their formation ages are poorly constrained. This paper presents a combined study of cathodoluminescence imaging, U–Pb age, trace element and Hf isotopes of zircon in migmatites from the NQ unit. In the migmatites, most zircon grains occur as new, homogeneous crystals, while some are present as overgrowth rims around inherited cores. Morphological and trace element features suggest that the zircon crystals are metamorphic and formed during partial melting. The inherited cores have oscillatory zoning and yield U–Pb ages of c. 900 Ma, representing their protolith ages. The early Neoproterozoic protoliths probably formed in an active continental margin, being a response to the assembly of the supercontinent Rodinia. The migmatite zircon yields Hf model ages of 1911 ± 20 to 990 ± 22 Ma, indicating that the protoliths were derived from reworking of Palaeoproterozoic to Neoproterozoic crustal materials. The anatexis zircon yields formation ages ranging from 455 ± 5 to 420 ± 4 Ma, with a peak at c. 435 Ma. Combined with previous results, we suggest that the migmatization of the NQ terrane occurred at c. 455–400 Ma. The migmatization was c. 50 Ma later than the c. 490 Ma ultra‐high‐P (UHP) metamorphism, indicating that they occurred in two independent tectonic events. By contrast, the migmatization was coeval with the granulite facies metamorphism and the granitic magmatism in the NQ unit, which collectively argue for their formation due to the northward subduction of the Shangdan Ocean. UHP rocks were distributed mainly along the northern margin and occasionally in the inner part of the NQ unit, indicating that they were exhumed along the northern edge and detached from the basement by the subsequent migmatization process.     

Syn‐metamorphic folding in the Tauern Window,Austria dated by Th–Pb ages from individual allanite porphyroblasts

High‐precision ²³²Th–²⁰⁸Pb dates have been obtained from allanite porphyroblasts that show unambiguous microstructural relationships to fabrics in a major syn‐metamorphic fold in the SE Tauern Window, Austria. Three porphyroblasts were analysed from a single garnet mica schist from the Peripheral Schieferhülle in the core of the Ankogel Synform, one of a series of folds which developed shortly before the thermal peak of Alpine epidote–amphibolite facies metamorphism: allanite grain 1 provided two analyses with a combined age of 27.7 ± 0.7 Ma; grain 2, which was slightly bent and fractured during crenulation, provided two analyses with a combined age of 27.7 ± 0.4 Ma; a single analysis from grain 3, which overgrew an already crenulated fabric, gave an age of 28.0 ± 1.4 Ma. The five ²³²Th–²⁰⁸Pb ages agree within error and define an isochron with an age of 27.71 ± 0.36 Ma (95% confidence level; MSWD = 0.46). The results imply that the crenulation event was in progress in a short interval (<1 Ma) c. 28 Ma, and that the Ankogel Synform was forming at this time. The thermal peak of regional metamorphism in the SE Tauern Window was probably attained shortly after 28 Ma, only c. 5 Ma after eclogite facies metamorphism in the central Tauern Window. Metasediment may contain allanite porphyroblasts with clear‐cut microstructural relationships to fabric development and metamorphic crystallization; for such rocks, ²³²Th–²⁰⁸Pb dating on microsamples offers a powerful geochronological tool.     

Intensity of in‐lake processes in floodplain lakes within the Bug River zone of fluvial activity

The paper presents a new approach to calculating the erosion and deposition values of floodplain lake basins, the erosion–deposition index (EDI). The EDI is a sum of the basin geometry indices (BGIs), which can be calculated for a separate cross section of the lake. The distribution of processes within the basin was investigated in two selected floodplain lakes with the use of BGIs. Field research was carried out in the Bug River valley from 1 November 2006 to 31 October 2011. The highest erosion was observed in the lakes located close to the parent river. Deposition processes were observed in lakes with high inflow of groundwater. The results showed that EDI values of 48 out of the 71 floodplain lakes ranged from ?0.2 to 0.2. Spatial distribution of erosion and deposition processes within the lake basins resulted from a velocity of water inflowing or flowing through the basin. This was observed especially in contrafluent–confluent lake. Inflow of rivers water via upstream crevasse occurred later than via downstream one, but energy of flowing water was higher, which favoured erosion of this part of the lake basin. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling sediment deposition and phosphorus retention in a river floodplain

Phosphorus (P) is one of the major limiting nutrient in many freshwater ecosystems. During the last decade, attention has been focused on the fluxes of suspended sediment and particulate P through freshwater drainage systems because of severe eutrophication effects in aquatic ecosystems. Hence, the analysis and prediction of phosphorus and sediment dynamics constitute an important element for ecological conservation and restoration of freshwater ecosystems. In that sense, the development of a suitable prediction model is justified, and the present work is devoted to the validation and application of a predictive soluble reactive phosphorus (SRP) uptake and sedimentation models, to a real riparian system of the middle Ebro river floodplain. Both models are coupled to a fully distributed two‐dimensional shallow‐water flow numerical model. The SRP uptake model is validated using data from three field experiments. The model predictions show a good accuracy for SRP concentration, where the linear regressions between measured and calculated values of the three experiments were significant (r² ≥ 0.62; p ≤ 0.05), and a Nash–Sutcliffe coefficient (E) that ranged from 0.54 to 0.62. The sedimentation model is validated using field data collected during two real flooding events within the same river reach. The comparison between calculated and measured sediment depositions showed a significant linear regression (p ≤ 0.05; r² = 0.97) and an E that ranged from 0.63 to 0.78. Subsequently, the complete model that includes flow dynamics, solute transport, SRP uptake and sedimentation is used to simulate and analyse floodplain sediment deposition, river nutrient contribution and SRP uptake. According to this analysis, the main SRP uptake process appears to be the sediment sorption. The analysis also reveals the presence of a lateral gradient of hydrological connectivity that decreases with distance from the river and controls the river matter contribution to the floodplain. Copyright © 2014 John Wiley & Sons, Ltd.     

Variability in the vertical hyporheic water exchange affected by hydraulic conductivity and river morphology at a natural confluent meander bend

River confluences and their associated tributaries are key morphodynamic nodes that play important roles in controlling hydraulic geometry and hyporheic water exchange in fluvial networks. However, the existing knowledge regarding hyporheic water exchange associated with river confluence morphology is relatively scarce. On January 14 and 15, 2016, the general hydraulic and morphological characteristics of the confluent meander bend (CMB) between the Juehe River and the Haohe River in the southern region of Xi'an City, Shaanxi Province, China, were investigated. The patterns and magnitudes of vertical hyporheic water exchange (VHWE) were estimated based on a one‐dimensional heat steady‐state model, whereas the sediment vertical hydraulic conductivity (K_v) was calculated via in situ permeameter tests. The results demonstrated that 6 hydrodynamic zones and their extensions were observed at the CMB during the test period. These zones were likely controlled by the obtuse junction angle and low momentum flux ratio, influencing the sediment grain size distribution of the CMB. The VHWE patterns at the test site during the test period mostly showed upwelling flow dominated by regional groundwater discharging into the river. The occurrence of longitudinal downwelling and upwelling patterns along the meander bend at the CMB was likely subjected to the comprehensive influences of the local sinuosity of the meander bend and regional groundwater discharge and finally formed regional and local flow paths. Additionally, in dominated upwelling areas, the change in VHWE magnitudes was nearly consistent with that in K_v values, and higher values of both variables generally occurred in erosional zones near the thalweg paths of the CMB, which were mostly made up of sand and gravel. This was potentially caused by the erosional and depositional processes subjected to confluence morphology. Furthermore, lower K_v values observed in downwelling areas at the CMB were attributed to sediment clogging caused by local downwelling flow. The confluence morphology and sediment K_v are thus likely the driving factors that cause local variations in the VHWE of fluvial systems.     

Significance of post‐peak metamorphic reaction microstructures in the ultrahigh temperature Eastern Ghats Province,India

Ultrahigh temperature (UHT) granulites in the Eastern Ghats Province (EGP) have a complex P–T–t history. We review the P–T histories of UHT metamorphism in the EGP and use that as a framework for investigating the P–T–t history of Mg–Al‐rich granulites from Anakapalle, with the express purpose of trying to reconcile the down‐pressure‐dominated P–T path with other UHT localities in the EGP. Mafic granulite that is host to Mg–Al‐rich metasedimentary granulites at Anakapalle has a protolith age of c. 1,580 Ma. Mg–Al‐rich metasedimentary granulites within the mafic granulite at Anakapalle were metamorphosed at UHT conditions during tectonism at 960–875 Ma, meaning that the UHT metamorphism was not the result of contact metamorphism from emplacement of the host mafic rock. Reworking occurred during the Pan‐African (c. 600–500 Ma) event, and is interpreted to have produced hydrous assemblages that overprint the post‐peak high‐T retrograde assemblages. In contrast to rocks elsewhere in the EGP that developed post‐peak cordierite, the metasedimentary granulites at Anakapalle developed post‐peak, generation ‘2’ reaction products that are cordierite‐absent and nominally anhydrous. Therefore, rocks at Anakapalle offer the unique opportunity to quantify the pressure drop that occurred during so‐called M2 that affected the EGP. We argue that M2 is either a continuation of M1 and that the overall P–T path shape is a complex counter‐clockwise loop, or that M1 is an up‐temperature counter‐clockwise deviation superimposed on the M2 path. Therefore, rather than the rocks at Anakapalle having a metamorphic history that is apparently anomalous from the rest of the EGP, we interpret that other previously studied localities in the EGP record a different part of the same P–T path history as Anakapalle, but do not preserve a significant record of pressure decrease. This is due either to the inability of refractory rocks to extensively react to produce a rich mineralogical record of pressure decrease, or because the earlier high‐P part of the rocks history was erased by the M1 loop. Irrespective of the specific scenario, models for the tectonic evolution of the EGP must take the substantial pressure decrease during M2 into account, as it is probable the P–T record at Anakapalle is a reflection of tectonics affecting the entire province.