Intraplate deformation and 3D rheological structure of the Rhine Rift System and adjacent areas of the northern Alpine foreland

The lithosphere of the Northern Alpine foreland has undergone a polyphase evolution during which interacting stress-induced intraplate deformation and upper mantle thermal perturbations controlled folding of the thermally weakened lithosphere. In this paper we address relationships among deeper lithospheric processes, neotectonics and surface processes in the Northern Alpine foreland with special emphasis on tectonically induced topography. We focus on lithosphere memory and neotectonics, paying special attention to the thermo-mechanical structure of the Rhine Graben System and adjacent areas of the northern Alpine foreland lithosphere. We discuss implications for mechanisms of large-scale intraplate deformation and links with surface processes and topography evolution.     

Estimation of Surface-Layer Scaling Parameters in the Unstable Boundary Layer: Implications for Orographic Flow Speed-Up

A method is proposed for estimating the surface-layer depth \((z_s)\) and the friction velocity \((u_*)\) as a function of stability (here quantified by the Obukhov length, L) over the complete range of unstable flow regimes. This method extends that developed previously for stable conditions by Argaín et al. (Boundary-Layer Meteorol 130:15–28, 2009), but uses a qualitatively different approach. The method is specifically used to calculate the fractional speed-up \((\varDelta S)\) in flow over a ridge, although it is suitable for more general boundary-layer applications. The behaviour of \(z_s \left( L\right) \) and \(u_*\left( L\right) \) as a function of L is indirectly assessed via calculation of \(\varDelta S\left( L\right) \) using the linear model of Hunt et al. (Q J R Meteorol Soc 29:16–26, 1988) and its comparison with the field measurements reported in Coppin et al. (Boundary-Layer Meteorol 69:173–199, 1994) and with numerical simulations carried out using a non-linear numerical model, FLEX. The behaviour of \(\varDelta S\) estimated from the linear model is clearly improved when \(u_*\) is calculated using the method proposed here, confirming the importance of accounting for the dependences of \(z_s\left( L \right) \) and \(u_*\left( L \right) \) on L to better represent processes in the unstable boundary layer.     

Adaptation of Pressure Based CFD Solvers for Mesoscale Atmospheric Problems

General purpose Computational Fluid Dynamics (CFD) solvers are frequently used in small-scale urban pollution dispersion simulations without a large extent of ver- tical flow. Vertical flow, however, plays an important role in the formation of local breezes, such as urban heat island induced breezes that have great significance in the ventilation of large cities. The effects of atmospheric stratification, anelasticity and Coriolis force must be taken into account in such simulations. We introduce a general method for adapting pressure based CFD solvers to atmospheric flow simulations in order to take advantage of their high flexibility in geometrical modelling and meshing. Compressibility and thermal stratification effects are taken into account by utilizing a novel system of transformations of the field variables and by adding consequential source terms to the model equations of incompressible flow. Phenomena involving mesoscale to microscale coupled effects can be analyzed without model nesting, applying only local grid refinement of an arbitrary level. Elements of the method are validated against an analytical solution, results of a reference calculation, and a laboratory scale urban heat island circulation experiment. The new approach can be applied with benefits to several areas of application. Inclusion of the moisture transport phenomena and the surface energy balance are important further steps towards the practical application of the method.     

Petrogenesis of low-δ<Superscript>18</Superscript>O quartz porphyry dykes,Koegel Fontein complex,South Africa

This paper investigates the origin of low-δ¹⁸O quartz porphyry dykes associated with the 144–133 Ma Koegel Fontein Igneous Complex, which was intruded during the initial phase of breakup of Africa and South America. The 25-km diameter Rietpoort Granite is the largest and youngest phase of activity, and is roofed by a 10-km diameter pendant of gneiss. Quartz porphyry (QP) dykes, up to 15 m in width, strike NW–SE across the complex. The QP dykes that intruded outside the granite have similar quartz phenocryst δ¹⁸O values (average 8.0‰, ± 0.7, n?=?33) to the granite (average 8.3?±?1.0, n?=?7). The QP dykes that intruded the roof pendant have quartz phenocrysts with more variable δ¹⁸O values (average 1.6‰, ± 2.1, n?=?55). In some cases quartz phenocrysts have δ¹⁸O values as low as ? 2.5‰. The variation in δ¹⁸O value within the quartz crystal population of individual dykes is small relative to the overall range, and core and rim material from individual quartz phenocrysts in three samples are identical within error. There is no evidence that quartz phenocryst δ¹⁸O values have been affected by fluid–rock interaction. Based on a ?_quartz?magma value of 0.6‰, magma δ¹⁸O values must have been as low as ? 3.1‰. Samples collected along the length of the two main QP dykes that traverse the roof pendant have quartz phenocryst δ¹⁸O values that range from +?1.1 to +?4.6‰, and ? 2.3 to +?5.6‰, respectively. These δ¹⁸O values correlate negatively (r = ? 0.96) with initial ⁸⁷Sr/⁸⁶Sr, which can be explained by the event that lowered δ¹⁸O values of the source being older than the dykes. We suggest that the QP dykes were fed by magma produced by partial melting of gneiss, which had been variably altered at high temperature by ¹⁸O-depleted meteoric water during global glaciation at ~?550 Ma. The early melts had variable δ¹⁸O value but as melt pockets interconnected during melting, the δ¹⁸O values approached that of average gneiss. Variable quartz phenocryst δ¹⁸O values in the same dyke can be explained by vertical emplacement, at variable rates of ascent along the dyke. The lateral variation in quartz, and hence magma δ¹⁸O value at a particular point along a single dyke would depend on the rate of ascent of magma at that point along the dyke, and the ‘age’ of the particular magma batch.     

Composition of barian mica in multiphase solid inclusions from orogenic garnet peridotites as evidence of mantle metasomatism in a subduction zone setting

Multiphase solid inclusions in minerals formed at ultra-high-pressure (UHP) provide evidence for the presence of fluids during deep subduction. This study focuses on barian mica, which is a common phase in multiphase solid inclusions enclosed in garnet from mantle-derived UHP garnet peridotites in the Saxothuringian basement of the northern Bohemian Massif. The documented compositional variability and substitution trends provide constraints on crystallization medium of the barian mica and allow making inferences on its source. Barian mica in the multiphase solid inclusions belongs to trioctahedral micas and represents a solid solution of phlogopite KMg₃(Si₃Al)O₁₀(OH)₂, kinoshitalite BaMg₃(Al₂Si₂)O₁₀(OH)₂ and ferrokinoshitalite BaFe₃(Al₂Si₂)O₁₀(OH)₂. In addition to Ba (0.24–0.67 apfu), mica is significantly enriched in Mg (X_Mg ~ 0.85 to 0.95), Cr (0.03–0.43 apfu) and Cl (0.04–0.34 apfu). The substitution vector involving Ba in the I-site which describes the observed chemical variability can be expressed as BaFe^IVAlClK_?1Mg_?1Si_?1(OH)_?1. A minor amount of Cr and ^VIAl enters octahedral sites following a substitution vector ^VI(Cr,Al)₂□^VI(Mg,Fe)_?3 towards chromphyllite and muscovite. As demonstrated by variable Ba and Cl contents positively correlating with Fe, barian mica composition is partly controlled by its crystal structure. Textural evidence shows that barian mica, together with other minerals in multiphase solid inclusions, crystallized from fluids trapped during garnet growth. The unusual chemical composition of mica reflects the mixing of two distinct sources: (1) an internal source, i.e. the host peridotite and its garnet, providing Mg, Fe, Al, Cr, and (2) an external source, represented by crustal-derived subduction-zone fluids supplying Ba, K and Cl. At UHP–UHT conditions recorded by the associated diamond-bearing metasediments (c. 1100 °C and 4.5 GPa) located above the second critical point in the pelitic system, the produced subduction-zone fluids transporting the elements into the overlying mantle wedge had a solute-rich composition with properties of a hydrous melt. The occurrence of barian mica with a specific chemistry in barium-poor mantle rocks demonstrates the importance of its thorough chemical characterization.     

Carbon costs and benefits of France’s biomass energy production targets

Background

Concern about climate change has motivated France to reduce its reliance on fossil fuel by setting targets for increased biomass-based renewable energy production. This study quantifies the carbon costs and benefits for the French forestry sector in meeting these targets. A forest growth and harvest simulator was developed for French forests using recent forest inventory data, and the wood-use chain was reconstructed from national wood product statistics. We then projected wood production, bioenergy production, and carbon balance for three realistic intensification scenarios and a business-as-usual scenario. These intensification scenarios targeted either overstocked, harvest-delayed or currently actively managed stands.

Results

All three intensification strategies produced 11.6–12.4 million tonnes of oil equivalent per year of wood-based energy by 2026, which corresponds to the target assigned to French wood-energy to meet the EU 2020 renewable energy target. Sustaining this level past 2026 will be challenging, let alone further increasing it. Although energy production targets can be reached, the management intensification required will degrade the near-term carbon balance of the forestry sector, compared to continuing present-day management. Even for the best-performing intensification strategy, i.e., reducing the harvest diameter of actively managed stands, the carbon benefits would only become apparent after 2040. The carbon balance of a strategy putting abandoned forests back into production would only break even by 2055; the carbon balance from increasing thinning in managed but untended stands would not break even within the studied time periods, i.e. 2015–2045 and 2046–2100. Owing to the temporal dynamics in the components of the carbon balance, i.e., the biomass stock in the forest, the carbon stock in wood products, and substitution benefits, the merit order of the examined strategies varies over time.

Conclusions

No single solution was found to improve the carbon balance of the forestry sector by 2040 in a way that also met energy targets. We therefore searched for the intensification scenario that produces energy at the lowest carbon cost. Reducing rotation time of actively managed stands is slightly more efficient than targeting harvest-delayed stands, but in both cases, each unit of energy produced has a carbon cost that only turns into a benefit between 2060 and 2080.

     

Predicting the impact of management and climate scenarios on groundwater nitrate concentration trends in southern Portugal

Hydrogeology Journal - The dynamics related to evolution of nitrate-contaminated groundwater are analyzed with focus on the impact of intrinsic aquifer properties, agricultural activities and...