Seismic damage analysis including inelastic shear–flexure interaction

The paper focusses on seismic damage analysis of reinforced concrete (R/C) members, accounting for shear–flexure interaction in the inelastic range. A finite element of the beam-column type recently proposed by the writers for the seismic analysis of R/C structures is first briefly described. The analytical model consists of two distributed flexibility sub-elements which interact throughout the analysis to simulate inelastic flexural and shear response. The finite element accounts for shear strength degradation with inelastic curvature demand, as well as coupling between inelastic flexural and shear deformations after flexural yielding. Based on this model, a seismic damage index is proposed taking into account both inelastic flexural and shear deformations, as well as their interaction. The finite element and the seismic damage index are used to analyse the response of R/C columns tested under cyclic loading and failing either in shear or in flexure. It is shown that the analytical model and damage index can predict and describe well the hysteretic response of R/C columns with different types of failure.     

What Can be Expected from the GRACE-FO Laser Ranging Interferometer for Earth Science Applications?

The primary objective of the gravity recovery and climate experiment follow-on (GRACE-FO) satellite mission, due for launch in August 2017, is to continue the GRACE time series of global monthly gravity field models. For this, evolved versions of the GRACE microwave instrument, GPS receiver, and accelerometer will be used. A secondary objective is to demonstrate the effectiveness of a laser ranging interferometer (LRI) in improving the satellite-to-satellite tracking measurement performance. In order to investigate the expected enhancement for Earth science applications, we have performed a full-scale simulation over the nominal mission lifetime of 5 years using a realistic orbit scenario and error assumptions both for instrument and background model errors. Unfiltered differences between the synthetic input and the finally recovered time-variable monthly gravity models show notable improvements with the LRI, on a global scale, of the order of 23 %. The gain is realized for wavelengths smaller than 240 km in case of Gaussian filtering but decreases to just a few percent when anisotropic filtering is applied. This is also confirmed for some typical regional Earth science applications which show randomly distributed patterns of small improvements but also degradations when using DDK4-filtered LRI-based models. Analysis of applied error models indicates that accelerometer noise followed by ocean tide and non-tidal mass variation errors are the main contributors to the overall GRACE-FO gravity model error. Improvements in these fields are therefore necessary, besides optimized constellations, to make use of the increased LRI accuracy and to significantly improve gravity field models from next-generation gravity missions.     

Transport of salt and freshwater in the Atlantic Subpolar Gyre

Transport of salt in the Irminger Current, the northern branch of the Atlantic Subpolar Gyre coupling the eastern and western subpolar North Atlantic, plays an important role for climate variability across a wide range of time scales. High-resolution ocean modeling and observations indicate that salinities in the eastern subpolar North Atlantic decrease with enhanced circulation of the North Atlantic subpolar gyre (SPG). This has led to the perception that a stronger SPG also transports less salt westward. In this study, we analyze a regional ocean model and a comprehensive global coupled climate model, and show that a stronger SPG transports more salt in the Irminger Current irrespective of lower salinities in its source region. The additional salt converges in the Labrador Sea and the Irminger Basin by eddy transports, increases surface salinity in the western SPG, and favors more intense deep convection. This is part of a positive feedback mechanism with potentially large implications for climate variability and predictability.     

Probability of earthquake occurrence and magnitude estimation in the post shut-in phase of geothermal projects

Induced seismicity in geothermal projects is observed to continue after shut-in of the fluid injection. Recent experiments show that the largest events tend to occur after the termination of injection. We use a probabilistic approach based on Omori??s law and the Gutenberg?CRichter magnitude?Cfrequency distribution to demonstrate that the probability of exceeding a certain maximum magnitude still increases after shut-in. This increase is governed by the exponent of Omori??s law q and the Gutenberg?CRichter b value. For a reduced b value in the post-injection phase, the probability of occurrence directly after shut-in can be even higher than the corresponding probability for an ongoing injection. For the reference case of q?=?2 and a 10% probability at shut-in time t _S to exceed a given maximum magnitude, we obtain an increase to 14.6% for t?=?2t _S at a constant Gutenberg?CRichter b value after shut-in. A reduction of the b value by one quarter leads to a probability of 20.5%. If we consider a constant probability level of occurrence for an event larger than a given magnitude at shut-in time, this maximum magnitude increases by 0.12?units for t?=?2t _S (0.26?units for a reduced b value). For the Soultz-sous-Forêts (France) injection experiment in 2000, recent studies reveal q?=?9.5 and a b value reduction by 14%. A magnitude 2.3 event 9?h after shut-in falls in the phase with a probability higher than for the continued injection. The probability of exceeding the magnitude of this post-injection event is determined to 97.1%.     

Comparison of an Electrochemical and Luminescence‐Based Oxygen Measuring System for Use in the Biodegradability Testing According to Closed Bottle Test (OECD 301D)

The effects of persistent organic pollutants on the aquatic environment still can hardly be assessed. Therefore, ready biodegradability is a desired key characteristic of chemicals. Accordingly, biodegradation testing is of utmost importance and is a vital part within the risk assessment for organic chemicals. In the closed bottle test (CBT; OECD 301D) the degree of biological degradation of a chemical is assessed by monitoring the oxygen consumption so far with a Clark electrode as standard method. However, this method is time‐consuming, laborious and led to frequent fluctuations of the test results, which limited their validity. Due to these shortcomings, an optical system that enables contact‐free measuring of the oxygen concentration in a closed bottle, based on dynamic luminescence quenching, was examined. The test results were compared to those obtained with the “classical” electrochemical technique (Clark electrode). At first, a basic standard for the handling of the new instrument in the laboratory was developed. Then possible influencing factors, reliability, comparability, and reproducibility of the test results were investigated by running the CBTs in parallel with the electrode and optode method. The findings proved the new optode method to be unambiguously superior to the electrode technique. The frequency of fluctuations of the test results and time and effort necessary for the test run could considerably be reduced. The degradation kinetics of the test substances could be followed easily in the CBT by increasing the measuring points without much additional effort.     

Sea‐level and ocean‐current control on carbonate‐platform growth,Maldives, Indian Ocean

Multichannel high‐resolution seismic and multibeam data were acquired from the Maldives‐isolated carbonate platform in the Indian Ocean for a detailed characterization of the Neogene bank architecture of this edifice. The goal of the research is to decipher the controlling factors of platform evolution, with a special emphasis on sea‐level changes and changes of the oceanic currents. The stacking pattern of Lower to Middle Miocene depositional sequences, with an evolution of a ramp geometry to a flat‐topped platform, reflects variations of accommodation, which here are proposed to be primarily governed by fluctuations of relative sea level. Easterly currents during this stage of bank growth controlled an asymmetric east‐directed progradation of the bank edge. During the late middle Miocene, this system was replaced by a twofold configuration of bank development. Bank growth continued synchronously with partial bank demise and associated sediment‐drift deposition. This turnover is attributed to the onset and/or intensification of the Indian monsoon and related upwelling and occurrence of currents, locally changing environmental conditions and impinging upon the carbonate system. Mega spill over lobes, shaped by reversing currents, formed as large‐scale prograding complexes, which have previously been interpreted as deposits formed during a forced regression. On a regional scale, a complex carbonate‐platform growth can occur, with a coexistence of bank‐margin progradation and aggradation, as well as partial drowning. It is further shown that a downward shift of clinoforms and offlapping geometries in carbonate platforms are not necessarily indicative for a sea‐level driven forced regression. Findings are expected to be applicable to other examples of Cenozoic platforms in the Indo‐Pacific region.     

Geology of the epithermal Ag–Au Huevos Verdes vein system and San José district, Deseado massif, Patagonia, Argentina

The San José district is located in the northwest part of the Deseado massif and hosts a number of epithermal Ag–Au quartz veins of intermediate sulfidation style, including the Huevos Verdes vein system. Veins are hosted by andesitic rocks of the Bajo Pobre Formation and locally by rhyodacitic pyroclastic rocks of the Chon Aike Formation. New ⁴⁰Ar/³⁹Ar constraints on the age of host rocks and mineralization define Late Jurassic ages of 151.3 ± 0.7 Ma to 144.7 ± 0.1 Ma for volcanic rocks of the Bajo Pobre Formation and of 147.6 ± 1.1 Ma for the Chon Aike Formation. Illite ages of the Huevos Verdes vein system of 140.8 ± 0.2 and 140.5 ± 0.3 Ma are 4 m.y. younger than the volcanic host rock unit. These age dates are among the youngest reported for Jurassic volcanism in the Deseado massif and correlate well with the regional context of magmatic and hydrothermal activity. The Huevos Verdes vein system has a strike length of 2,000 m, with several ore shoots along strike. The vein consists of a pre-ore stage and three main ore stages. Early barren quartz and chalcedony are followed by a mottled quartz stage of coarse saccharoidal quartz with irregular streaks and discontinuous bands of sulfide-rich material. The banded quartz–sulfide stage consists of sulfide-rich bands alternating with bands of quartz and bands of chlorite ± illite. Late-stage sulfide-rich veinlets are associated with kaolinite gangue. Ore minerals are argentite and electrum, together with pyrite, sphalerite, galena, chalcopyrite, minor bornite, covellite, and ruby silver. Wall rock alteration is characterized by narrow (< 3 m) halos of illite and illite/smectite next to veins, grading outward into propylitic alteration. Gangue minerals are dominantly massive quartz intergrown with minor to accessory adularia. Epidote, illite, illite/smectite, and, preferentially at deeper levels, Fe-chlorite gangue indicate near-neutral pH hydrothermal fluids at temperatures of >220°C. Kaolinite occurring with the late sulfide-rich veinlet stage indicates pH < 4 and a temperature of <200°C. The Huevos Verdes system has an overall strike of 325°, dipping on average 65° NE. The orientations of individual ore shoots are controlled by vein strike and intersecting north-northwest-striking faults. We propose a structural model for the time of mineralization of the San José district, consisting of a conjugate shear pair of sinistral north-northwest- and dextral west-northwest-striking faults that correspond to R and R′ in the Riedel shear model and that are related to master faults (M) of north-northeast-strike. Veins of 315° strike can be interpreted as nearly pure extensional fractures (T). Variations in vein strike predict an induced sinistral shear component for strike directions of >315°, whereas strike directions of <315° are predicted with an induced dextral strike–slip movement. The components of the structural model appear to be present on a regional scale and are not restricted to the San José district.     

Automated mapping of rock slope geometry,kinematics and stability with RSS-GIS

A GIS-implemented, deterministic approach for the automated spatial evaluation of geometrical and kinematical properties of rock slope terrains is presented. Based on spatially distributed directional information on planar geological fabrics and DEM-derived topographic attribute data, the internal geometry of rock slopes can be characterized on a grid cell basis. For such computations, different approaches for the analysis and regionalization of available structural directional information applicable in specific tectonic settings are demonstrated and implemented in a GIS environment. Simple kinematical testing procedures based on feasibility criteria can be conducted on a pixel basis to determine which failure mechanisms are likely to occur at particular terrain locations. In combination with hydraulic and strength data on geological discontinuities, scenario-based rock slope stability evaluations can be performed. For conceptual investigations on rock slope failure processes, a GIS-based specification tool for a 2-D distinct element code (UDEC) was designed to operate with the GIS-encoded spatially distributed rock slope data. The concepts of the proposed methodology for rock slope hazard assessments are demonstrated at three different test sites in Germany.     

Phase relations of phlogopite with magnesite from 4 to 8 GPa

To evaluate the stability of phlogopite in the presence of carbonate in the Earth’s mantle, we conducted a series of experiments in the KMAS–H₂O–CO₂ system. A mixture consisting of synthetic phlogopite (phl) and natural magnesite (mag) was prepared (phl₉₀-mag₁₀; wt%) and run at pressures from 4 to 8 GPa at temperatures ranging from 1,150 to 1,550°C. We bracketed the solidus between 1,200 and 1,250°C at pressures of 4, 5 and 6 GPa and between 1,150 and 1,200°C at a pressure of 7 GPa. Below the solidus, phlogopite coexists with magnesite, pyrope and a fluid. At the solidus, magnesite is the first phase to react out, and enstatite and olivine appear. Phlogopite melts over a temperature range of ~150°C. The amount of garnet increases above solidus from ~10 to ~30 modal% to higher pressures and temperatures. A dramatic change in the composition of quench phlogopite is observed with increasing pressure from similar to primary phlogopite at 4 GPa to hypersilicic at pressures ≥5 GPa. Relative to CO₂-free systems, the solidus is lowered such, that, if carbonation reactions and phlogopite metasomatism take place above a subducting slab in a very hot (Cascadia-type) subduction environment, phlogopite will melt at a pressure of ~7.5 GPa. In a cold (40 mWm⁻²) subcontinental lithospheric mantle, phlogopite is stable to a depth of 200 km in the presence of carbonate and can coexist with a fluid that becomes Si-rich with increasing pressure. Ascending kimberlitic melts that are produced at greater depths could react with peridotite at the base of the subcontinental lithospheric mantle, crystallizing phlogopite and carbonate at a depth of 180–200 km.     

Uncertainties of geochemical modeling during CO<Subscript>2</Subscript> sequestration applying batch equilibrium calculations

One of the uncertainties in the field of carbon dioxide capture and storage (CCS) is caused by the parameterization of geochemical models. The application of geochemical models contributes significantly to calculate the fate of the CO₂ after its injection. The choice of the thermodynamic database used, the selection of the secondary mineral assemblage as well as the option to calculate pressure dependent equilibrium constants influence the CO₂ trapping potential and trapping mechanism. Scenario analyses were conducted applying a geochemical batch equilibrium model for a virtual CO₂ injection into a saline Keuper aquifer. The amount of CO₂ which could be trapped in the formation water and in the form of carbonates was calculated using the model code PHREEQC. Thereby, four thermodynamic datasets were used to calculate the thermodynamic equilibria. Furthermore, the equilibrium constants were re-calculated with the code SUPCRT92, which also applied a pressure correction to the equilibrium constants. Varying the thermodynamic database caused a range of 61% in the amount of trapped CO₂ calculated. Simultaneously, the assemblage of secondary minerals was varied, and the potential secondary minerals dawsonite and K-mica were included in several scenarios. The selection of the secondary mineral assemblage caused a range of 74% in the calculated amount of trapped CO₂. Correcting the equilibrium constants with respect to a pressure of 125 bars had an influence of 11% on the amount of trapped CO₂. This illustrates the need for incorporating sensitivity analyses into reaction pathway modeling.