Fault slip analysis in the Koralm Massif (Eastern Alps) and consequences for the final uplift of “cold spots” in Miocene times

The Paleogene and Neogene evolution of Austroalpine basement units east of the Tauern Window is characterised by the formation of two major sets of faults: (1) ESE–WNW- to E–W-trending faults, associated with ENE- and NNW-trending conjugate structures and (2) N–S to NNE-SSW striking structures, mainly acting as high-angle normal faults, often associated with E-dipping low-angle normal faults along the western margin of the Styrian Basin.Together with the stratigraphic evolution of the Styrian and Lavanttal Basins and the related subsidence histories a tectonic evolution may be reconstructed for this part of the Eastern Alps. In the southern part of the Koralm Massif, WNW-trending fractures were activated as dextral strike-slip faults, associated with the evolution of WNW-trending troughs filled up with coarse block debris. W- to WNW-trending fractures were reactivated as normal faults, indicating N–S extension. It is assumed that these phases resulted in subsidence and block debris sedimentation in Karpatian and Badenian times (ca. 17–13 Ma).In the Western Styrian Basin no Sarmatian (13–11.5 Ma) sediments are observed; Pannonian (11.5 to 7.1 Ma) sediments are restricted to the Eastern Styrian Basin. This indicates, that the Koralm basement and the Western Styrian Basin were affected by post-Sarmatian uplift, coinciding with a re-activation of N-trending normal faults along the eastern margin of the Koralm Massif. Therefore, we suggest that the final uplift of the Koralm Complex, partly together with the Western Styrian Basin, occurred during the early Pannonian (at approximately 10 Ma). The elevation of clastic deposits indicates that the Koralm Complex was elevated by approximately 800 m during this phase, associated with an additional phase of E–W-directed extension accommodated by N–S striking normal faults.     

Joint atmospheric-terrestrial water balances for East Africa: a WRF-Hydro case study for the upper Tana River basin

Theoretical and Applied Climatology - For an improved understanding of the hydrometeorological conditions of the Tana River basin of Kenya, East Africa, its joint atmospheric-terrestrial water...     

A generalized cubic equation of state with application to pure CO<Subscript>2</Subscript> injection in aquifers

A generalized cubic equation of state is given. The Peng-Robinson and the Soave-Redlich-Kwong equations are special cases of this equation. The generalized equation of state is precisely as simple and computationally efficient as these classical equations. Through comparison with the Span-Wagner equation for CO ₂, we obtain an improved density accuracy in predefined temperature-pressure domains. The generalized equation is then verified through two relevant examples of CO ₂ injection and migration. Comparisons are made with other standard cubic EOS in order to show the range of solutions obtained with less accurate EOS.     

Lateral terrestrial water flow contribution to summer precipitation at continental scale – A comparison between Europe and West Africa with WRF-Hydro-tag ensembles

It is well accepted that summer precipitation can be altered by soil moisture condition. Coupled land surface – atmospheric models have been routinely used to quantify soil moisture – precipitation feedback processes. However, most of the land surface models (LSMs) assume a vertical soil water transport and neglect lateral terrestrial water flow at the surface and in the subsurface, which potentially reduces the realism of the simulated soil moisture – precipitation feedback. In this study, the contribution of lateral terrestrial water flow to summer precipitation is assessed in two different climatic regions, Europe and West Africa, for the period June–September 2008. A version of the coupled atmospheric-hydrological model WRF-Hydro with an option to tag and trace land surface evaporation in the modelled atmosphere, named WRF-Hydro-tag, is employed. An ensemble of 30 simulations with terrestrial routing and 30 simulations without terrestrial routing is generated with random realizations of turbulent energy with the stochastic kinetic energy backscatter scheme, for both Europe and West Africa. The ensemble size allows to extract random noise from continental-scale averaged modelled precipitation. It is found that lateral terrestrial water flow increases the relative contribution of land surface evaporation to precipitation by 3.6% in Europe and 5.6% in West Africa, which enhances a positive soil moisture – precipitation feedback and generates more uncertainty in modelled precipitation, as diagnosed by a slight increase in normalized ensemble spread. This study demonstrates the small but non-negligible contribution of lateral terrestrial water flow to precipitation at continental scale.     

End member and Bayesian mixing models consistently indicate near-surface flowpath dominance in a pristine humid tropical rainforest

The impacts of forest conversion on runoff generation in the tropics have received much interest, but scientific progress is still hampered by challenging fieldwork conditions and limited knowledge about runoff mechanisms. Here, we assessed the runoff generation, flow paths and water source dynamics of a pristine rainforest catchment in Costa Rica using end member mixing analysis (EMMA) and a Bayesian mixing model (MixSIAR). Geochemical tracer data collected over a 4-week field campaign were combined with tritium data used to assess potential deeper groundwater flow pathways to the perennial stream. The streamflow composition was best captured using three end-members, namely throughfall, shallow (5–15 cm) and deeper (15–50 cm) soil water. We estimated the end-member contributions to the main stream and two tributaries using the two mixing approaches and found good agreement between results obtained from EMMA and MixSIAR. The system was overwhelmingly dominated by near-surface sources, with little evidence for deeper and older groundwater as tritium-derived baseflow mean transit time was between 2.0 and 4.4 years. The shallow soil flow pathway dominated streamflow contributions in the main stream (median 39% and 49% based on EMMA and MixSIAR, respectively), followed by the deeper soil (32% and 31%) and throughfall (25% and 19%). The two tributaries had even greater shallow soil water contributions relative to the main stream (83% and 74% for tributary A and 42% and 63% for tributary B). Tributary B had no detectable deep soil water contribution, reflecting the morphology of the hillslope (steeper slopes, shallower soils and lower vegetation density compared to hillslope A). Despite the short sampling campaign and associated uncertainties, this study allowed to thoroughly assess runoff generation mechanisms in a humid tropical catchment. Our results also provide a first comparison of two increasingly used mixing models and suggest that EMMA and MixSIAR yield comparable estimates of water source partitioning in this tropical, volcanic rainforest environment.     

The effect of alkalis and polymerization on the solubility of H<Subscript>2</Subscript>O and CO<Subscript>2</Subscript> in alkali-rich silicate melts

The effect of alkalis on the solubility of H₂O and CO₂ in alkali-rich silicate melts was investigated at 500 MPa and 1,250 °C in the systems with H₂O/(H₂O + CO₂) ratio varying from 0 to 1. Using a synthetic analog of phonotephritic magma from Alban Hills (AH1) as a base composition, the Na/(Na + K) ratio was varied from 0.28 (AH1) to 0.60 (AH2) and 0.85 (AH3) at roughly constant total alkali content. The obtained results were compared with the data for shoshonitic and latitic melts having similar total alkali content but different structural characteristics, e.g., NBO/T parameter (the ratio of non-bridging oxygens over tetrahedrally coordinated cations), as those of the AH compositions. Little variation was observed in H₂O solubility (melt equilibrated with pure H₂O fluid) for the whole compositional range in this study with values ranging between 9.7 and 10.2 wt. As previously shown, the maximum CO₂ content in melts equilibrated with CO₂-rich fluids increases strongly with the NBO/T from 0.29 wt % for latite (NBO/T = 0.17) to 0.45 wt % for shoshonite (NBO/T = 0.38) to 0.90 wt % for AH2 (NBO/T = 0.55). The highest CO₂ contents determined for AH3 and AH1 are 1.18 ± 0.05 wt % and 0.86 ± 0.12 wt %, respectively, indicating that Na is promoting carbonate incorporation stronger than potassium. At near constant NBO/T, CO₂ solubility increases from 0.86 ± 0.12 wt % in AH1 [Na/(Na + K)] = 0.28, to 1.18 ± 0.05 wt % in AH3 [Na/(Na + K)] = 0.85, suggesting that Na favors CO₂ solubility on an equimolar basis. An empirical equation is proposed to predict the maximum CO₂ solubility at 500 MPa and 1,100–1,300 °C in various silicate melts as a function of the NBO/T, (Na + K)/∑cations and Na/(Na + K) parameters: \({\text{wt}}\% \;{\text{CO}}_{2} = - 0.246 + 0.014\exp \left( {6.995 \cdot \frac{\text{NBO}}{T}} \right) + 3.150 \cdot \frac{{{\text{Na}} + {\text{K}}}}{{\varSigma {\text{cations}}}} + 0.222 \cdot \frac{\text{Na}}{{{\text{Na}} + {\text{K}}}}.\) This model is valid for melt compositions with NBO/T between 0.0 and 0.6, (Na + K)/∑cation between 0.08 and 0.36 and Na/(Na + K) ratio from 0.25 to 0.95 at oxygen fugacities around the quartz–fayalite–magnetite buffer and above.     

The Effects of Temperature and Pressure on the Porosity Evolution of Flechtinger Sandstone

A porosity change influences the transport properties and the elastic moduli of rock while circulating water in a geothermal reservoir. The static and dynamic elastic moduli can be derived from the slope of stress–strain curves and velocity measurements, respectively. Consequently, the acoustic velocities were measured while performing hydrostatic drained tests. The effect of temperature on static and dynamic elastic moduli and porosity variations of Flechtinger sandstone was investigated in a wide range of confining pressure from 2 to 55 MPa. The experiments were carried out in a conventional triaxial system whereas the pore pressure remained constant, confining pressure was cycled, and temperature was increased step wise (25, 60, 90, 120, and 140 °C). The porosity variation was calculated by employing two different theories: poroelasticity and crack closure. The porosity variation and crack porosity were determined by the first derivative of stress–strain curves and the integral of the second derivative of stress–strain curves, respectively. The crack porosity analysis confirms the creation of new cracks at high temperatures. The porosity variation was increasing with an increase in temperature at low effective pressures and was decreasing with a rise in temperature at high effective pressures. Both compressional and shear wave velocities were increasing with increasing pressure due to progressive crack closure. Furthermore, the thermomechanical behavior of Flechtinger sandstone was characterized by an inversion effect where the sign of the temperature derivative of the drained bulk modulus changes.     

Evaluating the S-isotope fractionation associated with Phanerozoic pyrite burial

This study examines the sulfur isotope record of seawater sulfate proxies using δ³⁴S and Δ³³S to place constraints on the average global fractionation (Δ³⁴S_py) associated with pyrite formation and burial and the exponent λ that relates variations of the ³⁴S/³²S to variations of the ³³S/³²S. The results presented here use an analysis of the sulfur isotope record from seawater sulfate proxies and sedimentary sulfide to extract this quantity as the arithmetic difference between δ³⁴S of seawater sulfate and contemporaneous sulfide. It also uses an independent method that draws on inferences about the Δ³³S evolution of seawater sulfate to evaluate this further. These two methods yield similar results suggesting that Δ³⁴S_py and λ changed over the course of the Phanerozoic from slightly lower values of Δ³⁴S_py (lower values of λ) in the early Phanerozoic (Cambrian-Permian) to higher values of Δ³⁴S_py (higher values of λ) starting in the Triassic. This change of Δ³⁴S_py and the exponent λ is interpreted to reflect a change in the proportion of sulfide that was reoxidized and processed by bacterial disproportionation on a global scale. The revised record of Δ³⁴S_py also yields model pyrite burial curves making them more closely resemble model evolution curves for other element systems and global sea level curves. It is suggested that possible links to sea level may occur via changes in the area of submerged continental shelves which would provide additional loci for pyrite burial.The slightly different constraints used by the two approaches to calculate this fractionation may allow for additional information to be obtained about the sulfur cycle with future studies. For instance, the correspondence of these results suggests that the inferred variation of ³⁴S/³²S of pyrite is real, and that there is no significant missing sink of fractionated sulfur at the resolution of the present study (such as might be associated with organic sulfur). Burial of organic sulfur may, however, have been important at some times in the Phanerozoic and shorter timescale deviations between results provided by these methods may be observed with higher resolution sampling. If observed, this would suggest either that the record for pyrite (or less likely sulfate) is biased, or that another sink (possibly as organic sulfur) was important during these times in the Phanerozoic.     

Viscosity of andesitic melts—new experimental data and a revised calculation model

The viscosity of a synthetic andesite-like melt was measured in the low viscosity range (10¹-10⁶ Pa s) using the falling sphere(s) method and in the high viscosity range (10⁸-10¹³ Pa s) using parallel-plate viscometry. Falling sphere experiments with melts containing 2.3 and 5.6 wt.% H₂O were carried out in an internally heated gas pressure vessel (IHPV) at 500 MPa confining pressure. The sinking velocity of Pt and Pd spheres and in one case of a corundum sphere was used to measure the melt viscosity. In addition, a creep experiment was performed at ambient pressure using a glass containing 2.73 wt.% H₂O . A more water-rich glass (5.6 wt.% H₂O ) was investigated with a high pressure parallel-plate viscometer at 400 MPa confining pressure in an IPHV. By combining our new data with previous results for a similar melt composition we derived the following expression to describe the viscosity η (in Pa s) as a function of temperature T (in K) and water content w (in wt.%)