Vegetation and peat accumulation steer Holocene tidal–fluvial basin filling and overbank sedimentation along the Old Rhine River,The Netherlands

In the transformation from tidal systems to freshwater coastal landscapes, plants act as eco-engineering species that reduce hydrodynamics and trap sediment, but nature and timing of the mechanisms of land creation along estuaries remains unclear. This article focuses on the Old Rhine estuary (The Netherlands) to show the importance of vegetation in coastal landscape evolution, predominantly regarding tidal basin filling and overbank morphology. This estuary hosted the main outflow channel of the river Rhine between ca 6500 to 2000 cal bp , and was constrained by peat during most of its existence. This study reconstructs its geological evolution, by correlating newly integrated geological data and new field records to varying conditions. Numerical modelling was performed to test the inferred mechanisms. It was found that floodbasin vegetation and resulting organic accumulation strongly accelerated back-barrier infill, by minimizing tidal influence. After tidal and wave transport had already sufficiently filled the back-barrier basin, reed rapidly expanded from its edges under brackish conditions, as shown by diatom analysis and datings. Reed growth provided a positive infilling feedback by reducing tidal flow and tidal prism, accelerating basin infilling. New radiocarbon dates show that large-scale crevassing along the Old Rhine River – driven by tidal backwater effect – only started as nutrient-rich river water transformed the floodbasin into an Alder carr in a next phase of estuary evolution. Such less dense vegetation promotes crevassing as sediments are more easily transported into the floodbasin. As river discharge increased and estuary mouth infilling progressed, crevasse activity diminished around 3800 to 3000 cal bp , likely due to a reduced tidal backwater effect. The insights from this data-rich Holocene study showcase the dominant role that vegetation may have in the long-term evolution of coastal wetlands. It provides clues for effective use of vegetation in vulnerable wetland landscapes to steer sedimentation patterns to strategically adapt to rising water levels.     

A roadmap for a dedicated Earth Science Grid platform

Due to its intensive data processing and highly distributed organization, the multidisciplinary Earth Science applications community is uniquely positioned for the uptake and exploitation of Grid technologies. Currently Enabling Grids for E-sciencE, and other large Grid infrastructures are already deployed and capable of operational services. So far however, the adoption and exploitation of Grid technology throughout the Earth Science community has been slower than expected. The Dissemination and Exploitation of GRids in Earth sciencE project, proposed by the European Commission to assist and accelerate this process in a number of different ways, had between its main goals the creation of a roadmap towards Earth Science Grid platform. This paper presents the resulting roadmap.     

Effect of Fe on garnet–melt trace element partitioning: experiments in FCMAS and quantification of crystal-chemical controls in natural systems

Garnet–melt trace element partitioning experiments were performed in the system FeO–CaO–MgO–Al₂O₃–SiO₂ (FCMAS) at 3 GPa and 1540°C, aimed specifically at studying the effect of garnet Fe²⁺ content on partition coefficients (D^Grt/Melt). D^Grt/Melt, measured by SIMS, for trivalent elements entering the garnet X-site show a small but significant dependence on garnet almandine content. This dependence is rationalised using the lattice strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal–melt partition coefficients from elastic moduli. Nature 372, 452–454], which describes partitioning of an element i with radius r_i and valency Z in terms of three parameters: the effective radius of the site r₀(Z), the strain-free partition coefficient D₀(Z) for a cation with radius r₀(Z), and the apparent compressibility of the garnet X-site given by its Young's modulus E_X(Z). Combination of these results with data in Fe-free systems [Van Westrenen, W., Blundy, J.D., Wood, B.J., 1999. Crystal-chemical controls on trace element partitioning between garnet and anhydrous silicate melt. Am. Mineral. 84, 838–847] and crystal structure data for spessartine, andradite, and uvarovite, leads to the following equations for r₀(3+) and E_X(3+) as a function of garnet composition (X) and pressure (P):

r₀(3+) [Å]=0.930X_Py+0.993X_Gr+0.916X_Alm+0.946X_Spes+1.05(X_And+X_Uv)−0.005(P [GPa]−3.0)(±0.005 Å)

E_X(3+) [GPa]=3.5×10¹²(1.38+r₀(3+) [Å])^−26.7(±30 GPa)

Accuracy of these equations is shown by application to the existing garnet–melt partitioning database, covering a wide range of P and T conditions (1.8 GPa<P<5.0 GPa; 975°C<T<1640°C). D^Grt/Melt for all 3+ elements entering the X-site (REE, Sc and Y) are predicted to within 10–40% at given P, T, and X, when D^Grt/Melt for just one of these elements is known. In the absence of such knowledge, relative element fractionation (e.g. D_Sm^Grt/Melt/D_Nd^Grt/Melt) can be predicted. As an example, we predict that during partial melting of garnet peridotite, group A eclogite, and garnet pyroxenite, r₀(3+) for garnets ranges from 0.939±0.005 to 0.953±0.009 Å. These values are consistently smaller than the ionic radius of the heaviest REE, Lu. The above equations quantify the crystal-chemical controls on garnet–melt partitioning for the REE, Y and Sc. As such, they represent a major advance en route to predicting D^Grt/Melt for these elements as a function of P, T and X.     

Molecular modelling of rare earth element complexation in subduction zone fluids

Complexation of (trace) elements in fluids plays a critical role in determining element mobility in subduction zones, but to date, the atomic-scale processes controlling elemental solubilities are poorly understood. As a first step towards computer simulation of element complexation in subduction zone fluids, a thermodynamic cycle was developed to investigate the hydration environment and energetics of lanthanide complexes using density functional theory. The first solvation shell is explicitly defined and the remaining part of the aqueous fluid is modelled using a polarisable continuum model, which allows extrapolation to a broad pressure and temperature range.We illustrate our method by comparing solvation of lanthanide series elements in H₂O in the presence of fluoride or chloride for conditions relevant to subduction zones. The energetics of lanthanide- and lanthanide-fluoride/chloride hydration complexes were determined computationally. Calculated hydration free energies for trivalent lanthanides with explicit eight- and nine-fold coordinated first hydration shells show good agreement with literature data at room pressure and temperature. The hydration free energy is more negative for smaller complexes (heavy lanthanides) relative to larger complexes (light lanthanides), with the difference between La and Lu in water amounting to 361 kJ mol⁻¹. The hydration free energy of all lanthanide ions becomes less negative with increasing pressure (p) and temperature (T) up to 2.5 GPa and 1000 K (typical conditions in the upper part of subducting slabs). The free energy difference between light- and heavy-lanthanides remains essentially unchanged at elevated (p, T) conditions. There are minor geometrical differences in local hydration environment between light lanthanide-chloride (La-Nd) and heavy lanthanide-chloride (Pm-Lu) hydrated complexes, without a distinguishable energy difference. Complexation with fluoride is energetically more favourable than with chloride by 206 ± 4 kJ mol⁻¹ across the entire lanthanide series. The association of fluoride-water and chloride-water fragments with lanthanide-water complexes is energetically more favourable for aqueous lanthanide complexes surrounded by fewer first hydration shell water molecules.The methods developed in this study, in conjunction with simulation of the energetics of trace element incorporation into minerals, open the possibility to use molecular modelling to constrain elemental behaviour in subduction zones.     

Mantle resistance against Gibraltar slab dragging as a key cause of the Messinian Salinity Crisis

The Messinian Salinity Crisis (5.97–5.33 Ma) was caused by the closure of the Atlantic‐Mediterranean gateways that cut through the Gibraltar orogenic system. The geodynamic drivers underlying gateway closure and re‐opening are still debated. Here, we interrogate the gateway successions to find the imprints of surface deformation, infer the timing and nature of associated geodynamic drivers, and test such inferences against numerical simulations of slab dynamics. We find that since the latest Miocene, a tectonic framework was established in the gateway region dominated simultaneously by (a) relative plate convergence, (b) slab tearing under the eastern Betic Cordillera and (c) mantle resistance against north‐northeastward dragging of the Gibraltar slab by the African plate's absolute motion. We propose that mantle‐resisted slab dragging and slab tearing operated in concert closing the gateways that caused the Messinian Salinity Crisis, whereas sinking of heavy oceanic lithosphere located between buoyant continental plates re‐opened the Strait of Gibraltar at 5.33 Ma.     

Effect of attachment configuration on the trapping efficiency of Vaseline-coated slide catchers for windblown particles

There are various types of the windblown sediment traps developed for wind tunnel and field studies. One of the main supports expected from these traps is in measuring surface dust concentrations to appropriately derive flux equations. The measurement performance and accuracy of a trap is very important and depends strictly upon the physical characteristics and the behaviors of dust grains with air flows. This paper presents the measurement results of static pressure distribution (SPD) of wind flow around Vaseline-coated slide (VCS) catchers with an aim of finding out whether or not particle trapping efficiency (η) of the VCS is related to the SPD. The SPD was evaluated by a wind reduction coefficient (R _c) in a series of wind tunnel experiments with different VCS settings which have different attachment configurations on a pole. Three VCS configurations were considered: a configuration on a circular plastic pole (CPP) and two configurations on wooden square poles (WSP₁ and WSP₂, respectively). Thus, the primary contribution of this work was to experimentally analyze the effect of the different attachment configurations on the SPD, and the secondary objective was to determine the effect of the SPD on the η. It was shown that spatial correlation and spatial pattern of the R _c were different in the surrounding area of each configuration, and ANOVA and DUNCAN tests indicated that η(s) of WSP₁, WSP₂, and CPP were different at the significant level of P ≤ 0.05 with the mean of 0.94 ± 0.09, 0.63 ± 0.14, and 1.13 ± 0.07, respectively. Additionally, the amount of PM₂₀, PM₄₀, PM₆₀, PM₈₀, and PM₁₀₀ trapped by the configurations of WSP₁, WSP₂, and CPP considerably varied depending upon the particular aerodynamic circumstances associated with every configuration.     

A comparison of the evolution of arc complexes in Paleozoic interior and peripheral orogens: Speculations on geodynamic correlations

We discuss the potential geodynamic connections between Paleozoic arc development along the flanks of the interior (e.g. the Iapetus and Rheic) oceans and the exterior Paleopacific Ocean. Paleozoic arcs in the Iapetus and Rheic oceanic realms are preserved in the Appalachian–Caledonide and Variscan orogens, and in the Paleopacific Ocean realm they are preserved in the Terra Australis Orogen. Potential geodynamic connections are suggested by paleocontinental reconstructions showing Cambrian–Early Ordovician contraction of the exterior ocean as the interior oceans expanded, and subsequent Paleozoic expansion of the exterior oceans while the interior oceans contracted. Subduction initiated in the eastern segment of Iapetus at ca. 515 Ma and Early to Middle Ordovician orogenesis along the flanks of this ocean is highlighted by arc–continent collisions and ophiolite obductions. Over a similar time interval, subduction and orogenesis took place in the exterior ocean and included formation of the Macquarie arc in the Tasmanides of Eastern Australia and the Famatina arc and correlatives in the periphery of the proto-Andean margin of Gondwana. Major changes in the style of subduction (from retreating to advancing) in interior oceans occurred during the Silurian, following accretion of the peri-Gondwanan terranes and Baltica, and closure of the northeastern segment of Iapetus. During the same time interval, subduction in the Paleopacific Ocean was predominantly in a retreating mode, although intermittent episodes of contraction closed major marginal basins. In addition, however, there were major disturbances in the Earth tectonic systems during the Ordovician, including an unprecedented rise in marine life diversity, as well as significant fluctuations in sea level, atmospheric CO₂, and ⁸⁷Sr/⁸⁶Sr and ¹³C in marine strata carbonates. Stable and radiogenic isotopic data provide evidence for the addition of abundant mantle-derived magma, fluids and large mineral deposits that have a significant mantle-derived component. When considered together, the coeval, profound changes in the style of tectonic activity and the disturbances recorded in Earth Systems are consistent with the emergence of a superplume during the Ordovician. We speculate that the emergence of a superplume triggered by slab avalanche events within the Iapetus and Paleopacific oceans was associated with the establishment of a new geoid high within the Paleopacific regime, the closure of the interior Rheic Ocean and the amalgamation of Laurussia and Gondwana, which was a key event in the Late Carboniferous amalgamation of Pangea.     

Coastal flood risks and seasonal tourism: analysing the effects of tourism dynamics on casualty calculations

Since coastal tourism is one of the fastest growing sectors of tourism industry, coastal areas have become increasingly vulnerable in the case of flooding. While in recent years a number of different methods have been put forward to map coastal flood risks, the implications of tourism dynamics for the assessment of human casualties has remained largely overlooked in these models. This chapter examines to what extent the ignorance of (residential) coastal tourism may bias the calculations of human casualties. To this end, a case study has been conducted on the Belgian coast. Both the dynamic nature of coastal tourism and the behaviour of residential tourists in storm surge scenarios are considered. The results of this study show that including tourism dynamics in flood risk management is justified and appropriate, depending on the tourist attractiveness of the flood-prone area and its temporal fluctuations.