Differential aerobic and anaerobic oxidation of hydrocarbon gases discharged at mud volcanoes in the Nile deep-sea fan

The present study investigates hydrocarbon oxidation processes at Isis and Amon mud volcanoes (MV’s), in the eastern Nile deep-sea fan. In the water column, molecular and carbon isotopic signatures of light hydrocarbons indicate that gases rapidly dissolve in seawater and are partially oxidized.In the upper sediments, anaerobic oxidation of the light hydrocarbons takes place, as clearly shown by their molecular and isotopic composition. These processes lead to the presence of a distinct Sulfate-Hydrocarbon Interface at 120-145 cm and 20-50 cm below the seafloor, for Isis and Amon MV’s, respectively. In contrast to processes occurring in the water column, a clear preferential oxidation of methane, propane and n-butane over ethane and i-butane is observed in the anoxic sediments. Furthermore, for the first time, fractionation factors have been determined for the anaerobic oxidation of propane and butane, being respectively −4.80‰ and −0.7‰ for δ¹³C, and −43.3‰ for δ²H of propane.     

Emerging Technologies and Synergies for Airborne and Space-Based Measurements of Water Vapor Profiles

A deeper understanding of how clouds will respond to a warming climate is one of the outstanding challenges in climate science. Uncertainties in the response of clouds, and particularly shallow clouds, have been identified as the dominant source of the discrepancy in model estimates of equilibrium climate sensitivity. As the community gains a deeper understanding of the many processes involved, there is a growing appreciation of the critical role played by fluctuations in water vapor and the coupling of water vapor and atmospheric circulations. Reduction of uncertainties in cloud-climate feedbacks and convection initiation as well as improved understanding of processes governing these effects will result from profiling of water vapor in the lower troposphere with improved accuracy and vertical resolution compared to existing airborne and space-based measurements. This paper highlights new technologies and improved measurement approaches for measuring lower tropospheric water vapor and their expected added value to current observations. Those include differential absorption lidar and radar, microwave occultation between low-Earth orbiters, and hyperspectral microwave remote sensing. Each methodology is briefly explained, and measurement capabilities as well as the current technological readiness for aircraft and satellite implementation are specified. Potential synergies between the technologies are discussed, actual examples hereof are given, and future perspectives are explored. Based on technical maturity and the foreseen near-mid-term development path of the various discussed measurement approaches, we find that improved measurements of water vapor throughout the troposphere would greatly benefit from the combination of differential absorption lidar focusing on the lower troposphere with passive remote sensors constraining the upper-tropospheric humidity.     

The probabilistic seismic hazard assessment of Germany—version 2016, considering the range of epistemic uncertainties and aleatory variability

The basic seismic load parameters for the upcoming national design regulation for DIN EN 1998-1/NA result from the reassessment of the seismic hazard supported by the German Institution for Civil Engineering (DIBt). This 2016 version of the national seismic hazard assessment for Germany is based on a comprehensive involvement of all accessible uncertainties in models and parameters and includes the provision of a rational framework for integrating ranges of epistemic uncertainties and aleatory variabilities in a comprehensive and transparent way. The developed seismic hazard model incorporates significant improvements over previous versions. It is based on updated and extended databases, it includes robust methods to evolve sets of models representing epistemic uncertainties, and a selection of the latest generation of ground motion prediction equations. The new earthquake model is presented here, which consists of a logic tree with 4040 end branches and essential innovations employed for a realistic approach. The output specifications were designed according to the user oriented needs as suggested by two review teams supervising the entire project. Seismic load parameters, for rock conditions of \(v_{S30}\) = 800 m/s, are calculated for three hazard levels (10, 5 and 2% probability of occurrence or exceedance within 50 years) and delivered in the form of uniform hazard spectra, within the spectral period range 0.02–3 s, and seismic hazard maps for peak ground acceleration, spectral response accelerations and for macroseismic intensities. Results are supplied as the mean, the median and the 84th percentile. A broad analysis of resulting uncertainties of calculated seismic load parameters is included. The stability of the hazard maps with respect to previous versions and the cross-border comparison is emphasized.     

Low‐Resolution Modeling of Dense Drainage Networks in Confining Layers

Groundwater‐surface water (GW‐SW) interaction in numerical groundwater flow models is generally simulated using a Cauchy boundary condition, which relates the flow between the surface water and the groundwater to the product of the head difference between the node and the surface water level, and a coefficient, often referred to as the “conductance.” Previous studies have shown that in models with a low grid resolution, the resistance to GW‐SW interaction below the surface water bed should often be accounted for in the parameterization of the conductance, in addition to the resistance across the surface water bed. Three conductance expressions that take this resistance into account were investigated: two that were presented by Mehl and Hill (2010) and the one that was presented by De Lange (1999). Their accuracy in low‐resolution models regarding salt and water fluxes to a dense drainage network in a confined aquifer system was determined. For a wide range of hydrogeological conditions, the influence of (1) variable groundwater density; (2) vertical grid discretization; and (3) simulation of both ditches and tile drains in a single model cell was investigated. The results indicate that the conductance expression of De Lange (1999) should be used in similar hydrogeological conditions as considered in this paper, as it is better taking into account the resistance to flow below the surface water bed. For the cases that were considered, the influence of variable groundwater density and vertical grid discretization on the accuracy of the conductance expression of De Lange (1999) is small.     

Forward induced seismic hazard assessment: application to a synthetic seismicity catalogue from hydraulic stimulation modelling

The M _w 3.2-induced seismic event in 2006 due to fluid injection at the Basel geothermal site in Switzerland was the starting point for an ongoing discussion in Europe on the potential risk of hydraulic stimulation in general. In particular, further development of mitigation strategies of induced seismic events of economic concern became a hot topic in geosciences and geoengineering. Here, we present a workflow to assess the hazard of induced seismicity in terms of occurrence rate of induced seismic events. The workflow is called Forward Induced Seismic Hazard Assessment (FISHA) as it combines the results of forward hydromechanical-numerical models with methods of time-dependent probabilistic seismic hazard assessment. To exemplify FISHA, we use simulations of four different fluid injection types with various injection parameters, i.e. injection rate, duration and style of injection. The hydromechanical-numerical model applied in this study represents a geothermal reservoir with preexisting fractures where a routine of viscous fluid flow in porous media is implemented from which flow and pressure driven failures of rock matrix and preexisting fractures are simulated, and corresponding seismic moment magnitudes are computed. The resulting synthetic catalogues of induced seismicity, including event location, occurrence time and magnitude, are used to calibrate the magnitude completeness M _c and the parameters a and b of the frequency-magnitude relation. These are used to estimate the time-dependent occurrence rate of induced seismic events for each fluid injection scenario. In contrast to other mitigation strategies that rely on real-time data or already obtained catalogues, we can perform various synthetic experiments with the same initial conditions. Thus, the advantage of FISHA is that it can quantify hazard from numerical experiments and recommend a priori a stimulation type that lowers the occurrence rate of induced seismic events. The FISHA workflow is rather general and not limited to the hydromechanical-numerical model used in this study and can therefore be applied to other fluid injection models.     

Rapid crustal uplift in Patagonia due to enhanced ice loss

A vertical crustal uplift rate of 39 mm yr^? 1 is measured between 2003 and 2006 using Global Positioning System (GPS) measurements at the northeastern edge of the Southern Patagonia Icefield (SPI). This is the largest present-day glacial isostatic rate ever recorded. The combination of SPI's rapid melting and the unique regional slab-window tectonics that promotes a relatively low viscosity, is central to our interpretation of the observations. The two effects lead to a strong interaction of short relaxation times with ice loads that change on a comparable time scale. The profile of GPS observations link ice loss to the soft viscoelastic isostatic flow response over the time scale of the Little Ice Age (LIA), including ice loss in the period of observation. The agreement of the results with our model predictions strongly suggests the large crustal uplift in Patagonia is due an accelerated glacier wasting since the termination of the LIA and that the effective regional mantle viscosity is near 4.0–8.0 × 10¹⁸ Pa s. A century-long diminution of the icefields, at rates that are about 1/4 – 1/2 the contemporary loss rates, is consistent with multidecadal-scale temperature trends estimated for the past 50–100 years and is, in fact, a key feature in any model capable of explaining the uplift data.     

Phenocryst complexity in andesites and dacites from the Tequila volcanic field,Mexico: resolving the effects of degassing vs. magma mixing

The petrology of five phenocryst-poor (2–5%) andesites and dacites, all of which were erupted from different short-lived, monogenetic vents, is compared to that of phenocryst-rich (10–25%) andesites erupted from the adjacent stratovolcano, Volcán Tequila, in the Mexican arc. Despite differences in phenocryst abundances, these magmas have comparable phase assemblages (plagioclase + orthopyroxene + titanomagnetite + ilmenite + apatite ± augite ± hornblende), and similarly wide variations in phenocryst compositions, coupled to complex zoning patterns. For the phenocryst-poor lavas, equilibrium pairs of two Fe–Ti oxides lead to a narrow range of calculated temperatures for each sample that range from 934 (±24) to 1,073 (±6)°C and oxygen fugacities that range from +0.1 to +0.7 log units relative to the Ni–NiO buffer. Application of the plagioclase-liquid hygrometer to each sample at these calculated temperatures leads to maximum melt water concentrations of 4.6–3.1 wt% during plagioclase crystallization, indicating that the magmas were fluid saturated at depths ≥6.4–4.5 km. There is a wide, continuous range in the composition of plagioclase (≤44 mol% An) and orthopyroxene (≤16% Mg#) phenocrysts in each sample, which is consistent with a loss of dissolved water (≤2.8 wt%) from the melt phase during degassing as the magmas ascended rapidly to the surface. Evidence is presented that shows the effect of dissolved water is to reduce the activity of MgO relative to FeO in the melt phase, which indicates that degassing will also affect the Mg# of pyroxene phenocrysts, with higher melt water concentrations favoring Fe-rich pyroxene. Both plagioclase and orthopyroxene commonly display diffusion-limited growth textures (e.g., skeletal and hopper crystals, large interior melt hollows, and swallow tails), which are consistent with large undercoolings produced by degassing-induced crystallization. Therefore, degassing is proposed as a possible cause for the phenocryst compositional diversity documented in the phenocryst-poor andesite and dacite lavas erupted from peripheral vents, including the coexistence of normally zoned plagioclase and reversely zoned orthopyroxene. Degassing-induced crystallization may also explain some of the phenocryst complexity in crystal-rich andesites erupted from large stratovolcanoes, including Volcán Tequila.     

A quantitative reconstruction of organic matter and nutrient diagenesis in Mediterranean Sea sediments over the Holocene

A multicomponent diagenetic model was developed and applied to reconstruct the conditions under which the most recent sapropel, S1, was deposited in the eastern Mediterranean Sea. Simulations demonstrate that bottom waters must have been anoxic and sulphidic during the formation of S1 and that organic matter deposition was approximately three times higher than at present. Nevertheless, most present day sediment and pore water profiles — with the exception of pyrite, iron oxyhydroxides, iron-bound phosphorus and phosphate — can be reproduced under a wide range of redox conditions during formation of S1 by varying the depositional flux of organic carbon. As a result, paleoredox indicators (e.g., C_org:S ratio, C_org:P_org ratio, trace metals) are needed when assessing the contribution of oxygen-depletion and enhanced primary production to the formation of organic-rich layers in the geological record. Furthermore, simulations show that the organic carbon concentration in sediments is a direct proxy for export production under anoxic bottom waters.The model is also used to examine the post-depositional alteration of the organic-rich layer focussing on nitrogen, phosphorus, and organic carbon dynamics. After sapropel formation, remineralisation is dominated by aerobic respiration at a rate that is inversely proportional to the time since bottom waters became oxic once again. A sensitivity analysis was undertaken to identify the most pertinent parameters in regulating the oxidation of sapropels, demonstrating that variations in sedimentation rate, depositional flux of organic carbon during sapropel formation, bottom water oxygen concentration, and porosity have the largest impact. Simulations reveal that sedimentary nutrient cycling was markedly different during the formation of S1, as well as after reoxygenation of bottom waters. Accumulation of organic nitrogen in sediments doubled during sapropel deposition, representing a significant nitrogen sink. Following reventilation of deep waters, N₂ production by denitrification was almost 12 times greater than present day values. Phosphorus cycling also exhibits a strong redox sensitivity. The benthic efflux of phosphate was up to 3.5 times higher during the formation of S1 than at present due to elevated depositional fluxes of organic matter coupled with enhanced remineralisation of organic phosphorus. Reoxygenation of bottom waters leads to a large phosphate pulse to the water column that declines rapidly with time due to rapid oxidation of organic material. The oxidation of pyrite at the redox front forms iron oxyhydroxides that bind phosphorus and, thus, attenuate the benthic phosphate efflux. These results underscore the contrasting effects of oxygen-depletion on sedimentary nitrogen and phosphorus cycling. The simulations also confirm that the current conceptual paradigm of sapropel formation and oxidation is valid and quantitatively coherent.