Paleogeographical restoration and ramp tectonic evidence in Tunisian Tellian domain: Ain El Bey-Bou Awen area

The imbrication’s area in northern Tunisia is the most external segment of Alpine range, where several associated folds types with thrust ramps are recognized within imbricate units beneath Numidian front slope. Their presence help to understand thrusting mechanisms installation through studied area. In fact, this zone was considered as a result of Paleogene gravitary slop (Kujawski (Ann Miner Géol Tunis (24):281, 1969); Caire (Ann Min Géol Tunis 26:87–110, 1973); Rouvier 1977), which is proved to be affected by major deep decollement, given rise to various structures, some are propagation folds, specific of foreland front, limited to this area, and those in more external position: Tunisian Atlas (Creusot et al. (C R Acad Sci Paris 314(Sér II):961–965, 1992); Ouali and Mercier (PII: S0191-8141(97):00048-5, 1997); Ouali 1984; Ahmadi et al. (J Struct Geol 28:721–728, 2006)). Various categories of fold ramps could be identified: frontal folds ramp NE–SW and others as lateral or oblique ramp with NW–SE trend (Aridhi et al. (C R Geosci 343:360–369, 2011)). The relation between various structures has been used as recognition tools of thrusting sequences and to propose a new deformation chronology. Delimited outcropping of these structures between two both parallel faults strikes with regional displacement, leads to interpret these faults as cogenetic tear faults of propagation thrusts; this fault separates two domains with different deformation styles from each other side.     

Reply to ‘Comment on “The beginnings of hydrous mantle wedge melting” by Till et al.’ by Green, Rosenthal and Kovacs

The comment of Green et al. debates the interpretation of the temperature of the H₂O-saturated peridotite solidus and presence of silicate melt in the experiments of Till et al. (Contrib Mineral Petrol 163:669–688, 2012) at <1,000?°C. The criticisms presented in their comment do not invalidate any of the most compelling observations of Till et al. (Contrib Mineral Petrol 163:669–688, 2012) as discussed in the following response, including the changing minor element and Mg# composition of the solid phases with increasing temperature in our experiments with 14.5?wt% H₂O at 3.2?GPa, as well as the results of our chlorite peridotite melting experiments with 0.7?wt% H₂O. The point remains that Till et al. (Contrib Mineral Petrol 163:669–688, 2012) present data that call into question the H₂O-saturated peridotite solidus temperature preferred by Green (Tectonophysics 13(1–4):47–71, 1972; Earth Planet Sci Lett 19(1):37–53, 1973; Can Miner 14:255–268, 1976); Millhollen et al. (J Geol 82(5):575–587, 1974); Mengel and Green (Stability of amphibole and phlogopite in metasomatized peridotite under water-saturated and water-undersaturated conditions, Geological Society of Australia Special Publication, Blackwell, pp 571-581, 1989); Wallace and Green (Mineral Petrol 44:1–19, 1991) and Green et al. (Nature 467(7314):448–451, 2010).     

Higher Himalaya in the Bhagirathi section (NW Himalaya, India): its structures, backthrusts and extrusion mechanism by both channel flow and critical taper mechanisms

The Higher Himalayan Crystalline (HHC) in the Bhagirathi river section (India) on fieldwork reveals two extensional ductile top-to-N/NE shear sub-zones—the ‘South Tibetan Detachment System’ and the ‘Basal Detachment’—besides a preceding top-to-S/SW ductile shear. A top-to-N/NE brittle shear was identified as backthrusts from the HHC (except its northern portion) that occur repeatedly adjacent to numerous top-to-S/SW brittle shears as fore-thrusts. The northern portion of the HHC—the Gangotri Granite—exhibits infrequent total six extensional and compressional brittle shear senses. The backthrusts could be due to a low friction between the lower boundary of the HHC (i.e. the Main Central Thrust-Zone) and the partially molten hot rock materials of the HHC. Subduction of the Eurasian plate towards S/SW below the Indian plate more extensively in the Garhwal sector could be the second possible reason. Presence of two ductile extensional shear sub-zones may indicate channel flow (or several exhumation mechanisms) of the HHC in a shifting mode (similar to Mukherjee et al. in Int J Earth Sci 101:253–272, 2012). The top-to-S/SW extensional brittle shear exclusively within the upper (northern portion) of the HHC and a top-to-S/SW brittle shear within the remainder of it is a possible indicator of critical taper deformation mechanism. Thus, this work provides the field evidences of possibly both channel flow and critical taper conditions from a Higher Himalayan section, besides that by Larson et al. (Geol Soc Am Bull 122:1116–1134, 2010).     

Applicability of different ground-motion prediction models for northern Iran

A total of 163 free-field acceleration time histories recorded at epicentral distances of up to 200 km from 32 earthquakes with moment magnitudes ranging from M _w 4.9 to 7.4 have been used to investigate the predictive capabilities of the local, regional, and next generation attenuation (NGA) ground-motion prediction equations and determine their applicability for northern Iran. Two different statistical approaches, namely the likelihood method (LH) of Scherbaum et al. (Bull Seismol Soc Am 94:341–348, 2004) and the average log-likelihood method (LLH) of Scherbaum et al. (Bull Seismol Soc Am 99:3234–3247, 2009), have been applied for evaluation of these models. The best-fitting models (considering both the LH and LLH results) over the entire frequency range of interest are those of Ghasemi et al. (Seismol 13:499–515, 2009a) and Soghrat et al. (Geophys J Int 188:645–679, 2012) among the local models, Abrahamson and Silva (Earthq Spectra 24:67–97, 2008) and Chiou and Youngs (Earthq Spectra 24:173–215, 2008) among the NGA models, and finally Akkar and Bommer (Seism Res Lett 81:195–206, 2010) among the regional models.     

Out-of-sequence thrusts and paleogeography of the Rhenodanubian Flysch Belt (Eastern Alps) revisited

The Oberstdorf nappe of the Western and the Laab nappe of the Eastern Rhenodanubian Flysch (ERF) were independently identified as out-of-sequence thrust units by facies studies (Mattern 1999) and zircon analyses (Trautwein et al. 2001a, b, c), respectively. A new look at both areas reveals mutual similarities and new evidence for the out-of-sequence concept. Paleocurrent and heavy mineral data make it possible to reconstruct the sediment influx directions. From the Barremian to the mid-Campanian, the western and eastern basin segments were fed with south-derived garnet and north-derived zircon/”ZTR” (i.e., zircon, tourmaline, and rutile). Because both out-of-sequence units are relatively rich in zircon/ZTR they must have occupied the northernmost basin position. In the Western Rhenodanubian Flysch segment, the Sigiswang nappe occupied the central and the Üntschen nappe the southernmost basin position. In the ERF segment the central basin is represented by the Greifenstein nappe and the southernmost basin by the Kahlenberg nappe. Both out-of-sequence units do not occur in the northernmost and tectonically lowest position in their respective nappe piles as they were thrust over the other nappes. The reconstructed basin positions of the thrust units are suggested by the observation of a gradient in heavy mineral content in the thrust units. This paleogeographic arrangement is least problematic and renders models with differently positioned thrust units, requiring debris-shedding intrabasinal ridges, as unnecessarily complicated. Instead, we suggest that gradual changes in heavy mineral composition existed in across-basin direction. Garnet may stem from the Central Gneiss Complex of the Tauern window and formerly exposed lateral equivalents, all representing the southern Mid-Penninic zone. We assign the Falknis/Tasna nappe and formerly exposed lateral equivalents to the northern Mid-Penninic zone which served as the zircon/ZTR source. Interpreting Ebbing’s (Ph.D. thesis, Freie Universität Berlin, pp 1-143, 2002; Fig. 6.10) density section, we suggest that Mid-Penninic crust exists beneath the Central Gneiss Complex. During the latest Cretaceous much garnet was also N-derived. This may reflect processes related to the consumption of the North Penninic basin.     

Geomorphic analysis reveals active tectonic deformation on the eastern flank of the Pir Panjal Range,Kashmir Valley,India

There are plenty of faults that show evidence that they are active. Most of the valley’s floor is occupied by unconsolidated Karewa deposits, in particular on the south–southwest of the Kashmir Valley. In such situations, geomorphic data can reveal the location of active faults. Accordingly, we tried to identify geomorphic indices in SW of the Kashmir Valley (Veshav, Rambiara, and Romushi drainage basins), which revealed the area to be potentially tectonically active. This active faulting was further substantiated by drainage anomalies and field investigations, which provides evidence for an emergent out-of-sequence NE-dipping active reverse fault (identified first time on ground) named the Balapur Fault (BF). The BF can be traced over at least 40 km along the southwest side of the Kashmir Valley. The existence of the active Balapur Fault and of two other inferred faults north of the Panjal Thrust or Murree Thrust presents a picture of a more complex strain-partitioning regime in the Kashmir Himalayas than is usually visualized.     

Assessment of Liquefaction Potential of Guwahati City: A Case Study

The liquefaction potential of saturated cohesionless deposits in Guwahati city, Assam, was evaluated. The critical cyclic stress ratio required to cause liquefaction and the cyclic stress ratio induced by an earthquake were obtained using the simplified empirical method developed by Seed and Idriss (J soil Mech Found Eng ASCE 97(SM9):1249–1273, 1971, Ground motions and soil liquefaction during earthquakes. Earthquake Engineering Research Institute, Berkeley, CA, 1982) and Seed et al. (J Geotech Eng ASCE 109(3):458–483, 1983, J Geotech Eng ASCE 111(12):1425–1445, 1985) and the Idriss and Boulanger (2004) method. Critical cyclic stress ratio was based on the empirical relationship between standard penetration resistance and cyclic stress ratio. The liquefaction potential was evaluated by determining factor of safety against liquefaction with depth for areas in the city. A soil database from 200 boreholes covering an area of 262 km² was used for the purpose. A design peak ground acceleration of 0.36 g was used since Guwahati falls in zone V according to the seismic zoning map of India. The results show that 48 sites in Guwahati are vulnerable to liquefaction according to the Seed and Idriss method and 49 sites are vulnerable to liquefaction according to the Idriss and Boulanger method. Results are presented as maps showing zones of levels of risk of liquefaction.     

Comment on “A Comparison of Modified Fuzzy Weights of Evidence,Fuzzy Weights of Evidence,and Logistic Regression for Mapping Mineral Prospectivity” by Daojun Zhang,Frits Agterberg,Qiuming Cheng,and Renguang Zuo Math Geosci DOI 10.1007/s11004-013-9496-8

Despite a missing definition of equivalence of mathematical models or methods by Zhang et al. (Math Geosci, 2013), an “equivalence” (Zhang et al., Math Geosci, 2013, p. 6,7,8,14) of modified weights-of-evidence (Agterberg, Nat Resour Res 20:95–101, 2011) and logistic regression does not generally exist. Its alleged proof is based on a previously conjectured linear relationship between weights of evidence and logistic regression parameters (Deng, Nat Resour Res 18:249–258, 2009), which does not generally exist either (Schaeben and van den Boogaart, Nat Resour Res 20:401–406, 2011). In fact, an extremely simple linear relationship exists only if the predictor variables are conditionally independent given the target variable, in which case the contrasts, i.e., the differences of the weights, are equal to the logistic regression parameters. Thus, weights-of-evidence is the special case of logistic regression if the predictor variables are binary and conditionally independent given the target variable.     

Thermal metamorphism of the Arunachal Himalaya, India: Raman thermometry and thermochronological constraints on the tectono-thermal evolution

Determination of the peak thermal condition is vital in order to understand tectono-thermal evolution of the Himalayan belt. The Lesser Himalayan Sequence (LHS) in the Western Arunachal Pradesh, being rich in carbonaceous material (CM), facilitates the determination of peak metamorphic temperature based on Raman spectroscopy of carbonaceous material (RSCM). In this study, we have used RSCM method of Beyssac et al. (J Metamorph Geol 20:859–871, 2002a) and Rahl et al. (Earth Planet Sci Lett 240:339–354, 2005) to estimate the thermal history of LHS and Siwalik foreland from the western Arunachal Pradesh. The study indicates that the temperature of 700–800 °C in the Greater Himalayan Sequence (GHS) decreases to 650–700 °C in the main central thrust zone (MCTZ) and decreases further to <200 °C in the Mio-Pliocene sequence of Siwaliks. The work demonstrates greater reliability of Rahl et al.’s (Earth Planet Sci Lett 240:339–354, 2005) RSCM method for temperatures >600 and <340 °C. We show that the higher and lower zones of Bomdila Gneiss (BG) experienced temperature of ~600 °C and exhumed at different stages along the Bomdila Thrust (BT) and Upper Main Boundary Thrust (U.MBT). Pyrolysis analysis of the CM together with the Fission Track ages from upper Siwaliks corroborates the RSCM thermometry estimate of ~240 °C. The results indicate that the Permian sequence north of Lower MBT was deposited at greater depths (>12 km) than the upper Siwalik sediments to its south at depths <8 km before they were exhumed. The ⁴⁰Ar/³⁹Ar ages suggest that the upper zones of Se La evolved ~13–15 Ma. The middle zone exhumed at ~11 Ma and lower zone close to ~8 Ma indicating erosional unroofing of the MCT sheet. The footwall of MCTZ cooled between 6 and 8 Ma. Analyses of P–T path imply that LHS between MCT and U.MBT zone falls within the kyanite stability field with near isobaric condition. At higher structural level, the temperatures increase gradually with P–T conditions in the sillimanite stability field. The near isothermal (700–800 °C) condition in the GHS, isobaric condition in the MCTZ together with T–t path evidence of GHS that experienced relatively longer duration of near peak temperatures and rapid cooling towards MCTZ, compares the evolution of GHS and inverted metamorphic gradient closely to channel flow predictions.     

Convergence of iterative coupling for coupled flow and geomechanics

In this paper, we study solving iteratively the coupling of flow and mechanics. We demonstrate the stability and convergence of two widely used schemes: the undrained split method and the fixed stress split method. To our knowledge, this is the first time that such results have been rigorously obtained and published in the scientific literature. In addition, we propose a new stress split method, with faster convergence rate than known schemes. These results are specially important today due to the interest in hydraulic fracturing (Dean and Schmidt SPE J. 14:707–714, 2009; Ji et al. SPE J. 14:423–430, 2009; Samier and De Gennaro 2007; Settari and Maurits SPE J. 3:219–226, 1998), in oil and gas shale reservoirs.     

Melting phase relations of a mica–clinopyroxenite from the Milk River area, southern Alberta, Canada

Melting experiments were conducted on a mica–clinopyroxenite xenolith brought up in a minette dyke in southern Alberta, Canada, near Milk River. Both the minettes and mica–clinopyroxenite xenoliths were studied by Buhlmann et al. (Can J Earth Sci 37:1629–1650, 2000), who hypothesized that the minettes formed by partial melting of a mantle source containing clinopyroxene + phlogopite ± olivine, at pressures ≥1.7 GPa. In liquidus experiments performed on the most primitive minette in our previous study (Funk and Luth in Contrib Mineral Petrol 164:999–1009, 2012), we found a multiple saturation point where olivine and orthopyroxene coexisted with liquid at 1.77 GPa and 1,350 °C. We argued that the minette originally formed by partial melting of clinopyroxene + phlogopite, but had re-equilibrated with a harzburgite during ascent. In the current study, we wanted to test both the source region hypothesis of Buhlmann et al. and our re-equilibration hypothesis by studying the near-solidus phase equilibria of a mica + clinopyroxene assemblage. We found the solidus for our xenolith has a steep slope in P–T space and lies at temperatures above those of a normal cratonic geotherm, implying that this mica–clinopyroxenite is stable in the cratonic mantle. Melting could occur at greater depths, where the solidus is extrapolated to cross the geotherm or must be induced either by raising the temperatures of the surrounding rocks or by introducing hydrous fluids into the source. Our melts are in equilibrium with clinopyroxene and olivine. The compositions of the liquids derived from melting this xenolith are similar to madupitic lamproites from the Leucite Hills, Wyoming, studied by Carmichael (Contrib Mineral Petrol 15:24–66, 1967) and Barton and Hamilton (Contrib Mineral Petrol 66:41–49, 1978; Contrib Mineral Petrol 69:133–142, 1979). Barton and Hamilton (Contrib Mineral Petrol 69:133–142, 1979) proposed that the madupitic lamproites may have come from a source containing mica and pyroxene. This study supports their hypothesis. The composition of the most primitive minette from southern Alberta lies between our experimental melt and a population of representative mantle orthopyroxenes. We conclude from our study that the Milk River minettes were likely derived from a source containing phlogopite, clinopyroxene and trace amounts of apatite, which formed olivine upon melting. During ascent, the melts changed composition by reacting with orthopyroxene.     

Crystal reaming during the assembly, maturation, and waning of an eleven-million-year crustal magma cycle: thermobarometry of the Aucanquilcha Volcanic Cluster

Phenocryst assemblages of lavas from the long-lived Aucanquilcha Volcanic Cluster (AVC) have been probed to assess pressure and temperature conditions of pre-eruptive arc magmas. Andesite to dacite lavas of the AVC erupted throughout an 11-million-year, arc magmatic cycle in the central Andes in northern Chile. Phases targeted for thermobarometry include amphibole, plagioclase, pyroxenes, and Fe–Ti oxides. Overall, crystallization is documented over 1–7.5 kbar (~25 km) of pressure and ~680–1,110 °C of temperature. Pressure estimates range from ~1 to 5 kbar for amphiboles and from ~3 to 7.5 kbar for pyroxenes. Pyroxene temperatures are tightly clustered from ~1,000–1,100 °C, Fe–Ti oxide temperatures range from ~750–1,000 °C, and amphibole temperatures range from ~780–1,050 °C. Although slightly higher, these temperatures correspond well with previously published zircon temperatures ranging from ~670–900 °C. Two different Fe–Ti oxide thermometers (Andersen and Lindsley 1985; Ghiorso and Evans 2008) are compared and agree well. We also compare amphibole and amphibole–plagioclase thermobarometers (Ridolfi et al. 2010; Holland and Blundy 1994; Anderson and Smith 1995), the solutions from which do not agree well. In samples where we employ multiple thermometers, pyroxene temperature estimates are always highest, zircon temperature estimates are lowest, and Fe–Ti oxide and amphibole temperature estimates fall in between. Maximum Fe–Ti oxide and zircon temperatures are observed during the middle stage of AVC activity (~5–3 Ma), a time associated with increased eruption rates. Amphibole temperatures during this time are relatively restricted (~850–1,000 °C). The crystal record presented here offers a time-transgressive view of an evolving, multi-tiered subvolcanic reservoir. Some crystals in AVC lavas are likely to be true phenocrysts, but the diversity of crystallization temperatures and pressures recorded by phases in individual AVC lavas suggests erupting magma extensively reams and accumulates crystals from disparate levels of the middle to upper crust.     

The Use of Rock Mass Classification Systems to Estimate the Modulus and Strength of Jointed Rock

Three-dimensional, elastic and elasto-plastic finite element (FE) programs have permitted calculation of the displacements and the factor of safety (FOS) for the excavation for a tower, 132.70 m high (above foundation) on the island of Tenerife. The tower is supported by a 2 m thick reinforced concrete slab on jointed, vesicular and weathered basalt and scoria. The installation of rod extensometers at different depths below the slab has permitted comparison between measured and calculated displacements and the estimation of in situ deformation modulus. The moduli deduced from the simple empirical equations proposed by Hoek et al. (In: NARMS-TAC, 2002) and Gokceoglu et al. (Int J Rock Mech Min Sci 40:701–710, 2003) as a function of GSI, and Nicholson and Bieniawski (Int J Min Geol Eng 8:181–202, 1990) as a function of RMR, provide an acceptable fit with the measured settlements in this type of rock. Good correlation is also obtained with the empirical equation presented by Verman et al. (Rock Mech Rock Eng 30(3):121–127, 1997) that incorporates the influence of confining stress in the deformation modulus. The FOS obtained from different correlations with geomechanical classifications is within a relatively narrow range. These results increase our confidence in the use of classification schemes to estimate the deformation and stability in jointed rock.     

Little-known brachiopods from the Cretaceous of the Helvetic realm of NE Switzerland (Alpstein) and W Austria (Vorarlberg)

Three species from different stratigraphical levels of the Cretaceous of the Helvetic Alps are described. (1) Rhynchonellid specimens from the upper Öhrli-Kalk (Öhrli Formation, Late Berriasian) of NE Switzerland (Alpstein) identified as Lamellaerhynchia heimi (Sulser 2008) [Rhynchonellida, Hemithiridoidea]. Its range appears to be limited to a small area of the carbonate platform of the northern Alpstein chain. Based on internal and external morphological criteria L. heimi differs from other species of Lamellaerhynchia, as well as from Burrirhynchia cf. sayni (Jacob & Fallot 1913), occurring in the younger carbonate platform of the Schrattenkalk Formation (Early Aptian). (2) Recently collected material in various localities of the Altmann Member (Tierwis Formation, Late Hauterivian to Early Barremian) in the Alpstein area identified as Oblongarcula cf. alemannica Owen 1977 [Terebratellidina, Laqueoidea]. This species is closely related or identical to O. alemannica of the North European Boreal province and gives a reference to the occurrence of the genus Oblongarcula in the Tethyan domain of the Alps. Partially silicified specimens enable in rare cases a direct access to internal structures after that they were exposed by an acid treatment. (3) Tulipina koutaisensis (Loriol 1896) [Terebratellidina, Kingenoidea], known from Aptian deposits in the central Caucasus of Georgia, has been recorded as a rare species in the Plattenwald-Bed (Selun Member of the Garschella Formation, Albian) of W Austria (Vorarlberg). The localised occurrence and the temporal gap between the Caucasian and the Helvetic T. koutaisensis suggest an east–west directed migration along the northern margin of the Tethys Ocean during the Early Cretaceous.     

The Role of Shear Failure on Stress Characterization

Leak-off pressure and lost circulation data are generally thought to be reflective of minimum stress. We propose an alternative interpretation should be considered where the data may reflect a shear failure along zones of pre-existing weakness rather than opening of tensile fractures against the minimum stress. This mechanism has been discussed in a small number of borehole stability and hydraulic fracture papers, but has not been widely applied to leak-off test or lost circulation interpretation. In this paper, we will revisit and expand the concept introduced recently by Couzens-Schultz and Chan (J Struct Geol, doi:10.1016/j.jsg.2010.06.013, 2010) based on abnormally low leak-off tests in an active thrust belt to the analysis of lost circulation observations in modern-day deltaic environments. In the Gulf of Mexico, lost circulations historically are interpreted as a representation of the minimum horizontal stress due to initiating or reopening of a fracture in tensile mode. However, shear failure or fault reactivation can occur at pressures well below the minimum far-field stress that is typically considered a safe upper bound for mud pressure if pre-existing planes of weakness such as faults or fracture networks exist. We demonstrated a mud loss event is shown to be inconsistent with the tensile failure mode in a normal stress environment, but in good agreement with expectations for shear failure along pre-existing faults.     

Oxygen fugacity control in piston-cylinder experiments: a re-evaluation

Jakobsson (Contrib Miner Petrol 164(3):397–407, 2012) investigated a double capsule assembly for use in piston-cylinder experiments that would allow hydrous, high-temperature, and high-pressure experiments to be conducted under controlled oxygen fugacity conditions. Using a platinum outer capsule containing a metal oxide oxygen buffer (Ni–NiO or Co–CoO) and H₂O, with an inner gold–palladium capsule containing hydrous melt, this study was able to compare the oxygen fugacity imposed by the outer capsule oxygen buffer with an oxygen fugacity estimated by the AuPdFe ternary system calibrated by Barr and Grove (Contrib Miner Petrol 160(5):631–643, 2010). H₂O loss or gain, as well as iron loss to the capsule walls and carbon contamination, is often observed in piston-cylinder experiments and often go unexplained. Only a few have attempted to actually quantify various aspects of these changes (Brooker et al. in Am Miner 83(9–10):985–994, 1998; Truckenbrodt and Johannes in Am Miner 84:1333–1335, 1999). It was one of the goals of Jakobsson (Contrib Miner Petrol 164(3):397–407, 2012) to address these issues by using and testing the AuPdFe solution model of Barr and Grove (Contrib Miner Petrol 160(5):631–643, 2010), as well as to constrain the oxygen fugacity of the inner capsule. The oxygen fugacities of the analyzed melts were assumed to be equal to those of the solid Ni–NiO and Co–CoO buffers, which is incorrect since the melts are all undersaturated in H₂O and the oxygen fugacities should therefore be lower than that of the buffer by 2 log $a_{{{\text{H}}_{ 2} {\text{O}}}}$ .     

A phase-by-phase upstream scheme that converges to the vanishing capillarity solution for countercurrent two-phase flow in two-rock media

We discuss the convergence of the upstream phase-by-phase scheme (or upstream mobility scheme) towards the vanishing capillarity solution for immiscible incompressible two-phase flows in porous media made of several rock types. Troubles in the convergence were recently pointed out by Mishra and Jaffré (Comput. Geosci. 14, 105–124, 2010) and Tveit and Aavatsmark (Comput. Geosci. 16, 809–825, 2012). In this paper, we clarify the notion of vanishing capillarity solution, stressing the fact that the physically relevant notion of solution differs from the one inferred from the results of Kaasschieter (Comput. Geosci. 3, 23–48, 1999). In particular, we point out that the vanishing capillarity solution depends on the formally neglected capillary pressure curves, as it was recently proven in by Andreianov and Cancès (Comput. Geosci. 17, 551–572, 2013). Then, we propose a numerical procedure based on the hybridization of the interfaces that converges towards the vanishing capillarity solution. Numerical illustrations are provided.     

Flow rule in a high-cycle accumulation model backed by cyclic test data of 22 sands

The flow rule used in the high-cycle accumulation (HCA) model proposed by Niemunis et al. (Comput Geotech 32: 245, 2005) is examined on the basis of the data from approximately 350 drained long-term cyclic triaxial tests (N = 10⁵ cycles) performed on 22 different grain-size distribution curves of a clean quartz sand. In accordance with (Wichtmann et al. in Acta Geotechnica 1: 59, 2006), for all tested materials, the “high-cyclic flow rule (HCFR)”, i.e., the ratio of the volumetric and deviatoric strain accumulation rates \(\dot{\varepsilon}_{\rm{v}}^{{\rm acc}}/\dot{\varepsilon}_{\rm{q}}^{{\rm acc}}\) , was found dependent primarily on the average stress ratio η ^av = q ^av/p ^av and independent of amplitude, soil density and average mean pressure. The experimental HCFR can be fairly well approximated by the flow rule of the modified Cam-clay (MCC) model. Instead of the critical friction angle \(\varphi_{\rm{c}}\) which enters the flow rule for monotonic loading, the HCA model uses the MCC flow rule expression with a slightly different parameter \(\varphi_{\rm{cc}}\) . It should be determined from cyclic tests. \(\varphi_{\rm{cc}}\) and \(\varphi_{\rm{c}}\) are of similar magnitude but not always identical, because they are calibrated from different types of tests. For a simplified calibration in the absence of cyclic test data, \(\varphi_{\rm{cc}}\) may be estimated from the angle of repose \(\varphi_{\rm{r}}\) determined from a pluviated cone of sand (Wichtmann et al. in Acta Geotechnica 1: 59, 2006). However, the paper demonstrates that the MCC flow rule with \(\varphi_{\rm{r}}\) does not fit well the experimentally observed HCFR in the case of coarse or well-graded sands. For an improved simplified calibration procedure, correlations between \(\varphi_{\rm{cc}}\) and parameters of the grain-size distribution curve (d ₅₀,   C _u) have been developed on the basis of the present data set. The approximation of the experimental HCFR by the generalized flow rule equations proposed in (Wichtmann et al. in J Geotech Geoenviron Eng ASCE 136: 728, 2010), considering anisotropy, is also discussed in the paper.     

Investigation of the variability of strong ground motions from Vrancea earthquakes

A systematic investigation of the applicability of several ground motion prediction models for Vrancea intermediate-depth seismic source is conducted in this research. Two ground motion prediction models recommended by previous evaluations (Vacareanu et al. in Bull Earthq Eng 11(6):1867–1884, 2013a; Pavel et al. in Earthq Struct 6(1):1–18, 2014), as well as two new state-of-the-art ground motion prediction equations (Vacareanu et al. in J Earthq Eng, 2013b; Earthq Struct 6(2):141–161, 2014) are tested using an increased strong ground motion database consisting of 150 recordings from Vrancea subcrustal earthquakes. The evaluation is performed by using several goodness-of-fit parameters from the literature. Moreover, the applicability of the single-station sigma method is also investigated by using the same strong ground motion database recorded in 30 seismic stations from southern and eastern Romania. The influence of the soil conditions on the numerical results obtained in this study is investigated and discussed using the results provided by the analysis of variance method. The impact of the single-station standard deviation on the levels of seismic hazard is also assessed in this study, and the results show, in the analyzed cases, significant reductions of the hazard levels.     

Chromite alteration processes within Vourinos ophiolite

The renewed interest in chromite ore deposits is directly related to the increase in Cr price ruled by international market trends. Chromite, an accessory mineral in peridotites, is considered to be a petrogenetic indicator because its composition reflects the degree of partial melting that the mantle experienced while producing the chromium spinel-bearing rock (Burkhard in Geochim Cosmochim Acta 57:1297–1306, 1993). However, the understanding of chromite alteration and metamorphic modification is still controversial (e.g. Evans and Frost in Geochim Cosmochim Acta 39:959–972, 1975; Burkhard in Geochim Cosmochim Acta 57:1297–1306, 1993; Oze et al. in Am J Sci 304:67–101, 2004). Metamorphic alteration leads to major changes in chromite chemistry and to the growth of secondary phases such as ferritchromite and chlorite. In this study, we investigate the Vourinos complex chromitites (from the mines of Rizo, Aetoraches, Xerolivado and Potamia) with respect to textural and chemical analyses in order to highlight the most important trend of alteration related to chromite transformation. The present study has been partially funded by the Aliakmon project in collaboration between the Public Power Corporation of Greece and Institute of Geology and Mineral Exploration of Kozani.