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).     

Pushover and nonlinear time history analysis evaluation of a RC building collapsed during the Van (Turkey) earthquake on October 23, 2011

The nonlinear seismic behavior of a collapsed reinforced concrete (RC) residential building in the city of Van in Turkey is investigated by the static pushover and nonlinear time history analyses. The selected RC structure was designed according to the 1975 version of Turkish Earthquake Code (TEC-1975). The building had experienced heavy damage, and it was demolished in the Van earthquake on October 23, 2011. The 2007 version of Turkish Earthquake Code (TEC-2007) is considered for the assessing seismic performance evaluation of the selected RC building. The RC structure presents collapse performance level under the earthquake loads. Besides, the analytical solutions show that different performance levels for the sections are obtained from the pushover and nonlinear time history methods.     

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.     

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.     

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.     

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.     

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.     

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.     

Infiltration effects on a two-dimensional molecular dynamics model of landslides

We propose a two-dimensional computational model for deep landslides triggered by rainfall, based on interacting particles or grains. The model describes a vertical section of a fictitious granular material along a slope, in order to study the behavior of a wide-thickness landslide. The triggering of the landslide is caused by the exceeding of two conditions: a threshold speed and a condition on the static friction of the particles, the latter based on the Mohr–Coulomb failure criterion (Coulomb in Mem Acad R Div Sav 7:343–387, 1776; Mohr in Abhandlungen aus dem Gebiete der Technischen Mechanik. Ernst, Berlin, 1914). The interparticle interactions are represented as a potential that, in the absence of suitable experimental data and due to the arbitrariness of the grain dimension, is modeled similarly to the Lennard-Jones’ one (Lennard-Jones in Proc R Soc Lond A 106(738):463–477, 1924), i.e., with an attractive and a repulsive part. For the updating of the particle positions, we use a molecular dynamics method, which is quite suitable for this type of systems (Herrmann and Luding in Continuum Mech Thermodyn 10:189–231, 1998). An infiltration scheme is introduced for modeling the increasing pore pressure due to the rainfall. Finally, we also introduce the viscosity in the dynamical equations of motion. The statistical characterization and dynamical behavior of the results of simulations are quite satisfactory relative to real landslides: We obtain a power law distribution of landslide triggering times, and the velocity patterns are typical of real cases, including the acceleration progression. Therefore, we can claim that this type of modeling can represent a new method to simulate landslides triggered by rainfall.     

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.     

Effect of frequency-dependent radiation pattern in the strong motion simulation of the 2011 Tohoku earthquake,Japan, using modified semi-empirical method

We perform a strong ground motion simulation using a modified semi-empirical technique (Midorikawa in Tectonophysics 218:287–295, 1993), with frequency-dependent radiation pattern model. Joshi et al. (Nat Hazards 71:587–609, 2014) have modified the semi-empirical technique to incorporate the modeling of strong motion generation areas (SMGAs). A frequency-dependent radiation pattern model is applied to simulate high-frequency ground motion more precisely. Identified SMGAs (Kurahashi and Irikura in Earth Planets Space 63:571–576, 2011) of the 2011 off the Pacific coast of Tohoku earthquake (M _w  = 9.0) were modeled using this modified technique. We analyzed the effect of changing seismic moment values of SMGAs on the simulated acceleration time series. Final selection of the moment values of SMGAs is based on the root-mean-square error (RMSE) of waveform comparison. Records are simulated for both frequency-dependent and constant radiation pattern function. Simulated records for both cases are compared with observed records in terms of peak ground acceleration, peak ground velocity and pseudo-acceleration response spectra at different stations. Comparison of simulated and observed records in terms of RMSE suggests that the method is capable of simulating record, which matches in a wide frequency range for this earthquake and bears realistic appearance in terms of shape and strong motion parameters. The results confirm the efficacy and suitability of rupture model defined by five SMGAs for the developed modified technique.     

Seismic hazard analysis of Sinop province,Turkey using probabilistic and statistical methods

Using 4.0 and greater magnitude earthquakes which occurred between 1 January 1900 and 31 Dec 2008 in the Sinop province of Turkey this study presents a seismic hazard analysis based on the probabilistic and statistical methods. According to the earthquake zonation map, Sinop is divided into first, second, third and fourth-degree earthquake regions. Our study area covered the coordinates between 40.66°– 42.82°N and 32.20°– 36.55°E. The different magnitudes of the earthquakes during the last 108 years recorded on varied scales were converted to a common scale (Mw). The earthquake catalog was then recompiled to evaluate the potential seismic sources in the aforesaid province. Using the attenuation relationships given by Boore et al. (1997) and Kalkan and Gülkan (2004), the largest ground accelerations corresponding to a recurrence period of 475 years are found to be 0.14 g for bedrock at the central district. Comparing the seismic hazard curves, we show the spatial variations of seismic hazard potential in this province, enumerating the recurrence period in the order of 475 years.     

Impedance Functions of Slab Foundations with Rigid Piles

The earthquake resistant design of structures depends upon the soil-structure interaction during seismic excitation. The dynamic behavior of surface foundations and deep foundations has been investigated by several authors (Gazetas in Soil Dyn Earthq Eng 2(1):2–42, 1983; Pecker in Dynamique des sols. Presses de l’Ecole Nationale des Ponts et Chaussées, Paris, 1984; Sieffert and Cevaer in Manuel des fonctions d’impédance. Ouest-Editions, Nantes, 1992, etc.) but the dynamic behavior of pile-reinforced soils has not been sufficiently studied yet. The design of pile-reinforced soils comprises the design of the rigid piles and the design of the earth-platform. The foundation system studied in this article consists of an earth-platform over a soft ground reinforced by deep piles. In order to understand the dynamic behavior of a rigid pile, a series of dynamic tests is conducted on an experimental site. The vertical and horizontal impedances of the slab foundation are obtained with and without rigid piles. The numerical models are developed to interpret these dynamic tests. The numerical and experimental impedance functions are compared in both the vertical and the horizontal directions. A sensitivity analysis on the influence of the physical and geometrical properties of rigid piles on the impedance functions is presented.     

Discussion of “The Bushveld Complex,South Africa: formation of platinum–palladium,chrome- and vanadium-rich layers via hydrodynamic sorting of a mobilized cumulate slurry in a large,relatively slowly cooling,subsiding magma chamber” by Maier et al. (2013) Miner Deposita 48:1–56

We believe the hypothesis presented by Maier et al. (Miner Deposita 48:1–56, 2012) for the formation of the various ore bodies in the Bushveld Complex to be overly simplistic, and we raise concerns that some of our work, used in support of this hypothesis, has been misrepresented. The formation of both diverse metalliferous layers (platinum-group element (PGE) reefs and Ti-magnetite layers) and some discordant (pipe) ore deposits has been ascribed by Maier et al. to the single unifying process of hydrodynamic sorting. The problem faced by authors of universal hypotheses for the Bushveld Complex is the sheer size and complexity of the intrusion. We disagree with many aspects of the overall Maier et al. model and have also identified several minor errors on maps and photographs, although some of these do not have a material effect on the model. The nature and origin of the layering is, however, too complex a topic to deal with in the context of this commentary, and we restrict ourselves to noting that our preferred hypothesis, namely the incremental buildup of layering from numerous episodes of replenishment, by different magma lineages, is consistent with field relationships. Our hypothesis for the origin of the ultramafic-hosted PGE-rich reefs, i.e., lateral mixing, is applicable to economically mineralized reefs (Mitchell and Scoon, Econ Geol 102:971–1009, 2007) and poorly mineralized layers such as the Pseudoreef harzburgite (Scoon and De Klerk, Canad Mineral 25:51–77, 1987) and the chromitite layers below the UG2 (Scoon and Teigler, Econ Geol 89:1094–1121, 1994).     

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.     

Indirect economic losses of drought under future projections of climate change: a case study for Spain

This study presents the results of numerical simulations of the 2004 Indian Ocean earthquake and tsunami in the Bay of Lhok Nga (northwestern coast of Sumatra, Indonesia) integrating sediment erosion and deposition. We investigate the transport of sediment both by suspension and by bedload under different scenarii of long breaking dispersive waves through a series of numerical experiments. The earthquake source model used by Koshimura et al. (Coast Eng J 51:243–273, 2008) with a 25-m dislocation better reproduces the wave travel time, flow depth and inundation area than the other models tested. The model reproduces realistically the pronounced coastal retreat in the northern part of Lhok Nga Bay (retreat ranging between 50 and 150 m), where Paris et al. (Geomorphology 104:59–72, 2009) estimated a mean retreat of 80 m. There is also a good agreement between the simulated area of coastal retreat (195,400 m²) and the field observations (203,200 m²). The simulation may underestimate the volume of tsunami deposits (611,700 m³ vs. 500,000–1,000,000 m³ estimated by Paris et al. (2009). The model fully reproduces the observed thickness of tsunami deposits when considering both bedload and suspension, even if bedload transport dominates. Limitations are due to micro-scale topographic, anthropic features (which are not always represented by the DEM) and the amount of debris which may influence flow dynamics and sediment transport.     

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.     

Earthquake geology of Kashmir Basin and its implications for future large earthquakes

Two major traces of active thrust faults were identified in the Kashmir Basin (KB) using satellite images and by mapping active geomorphic features. The ~N130°E strike of the mapped thrust faults is consistent with the regional ~NE–SW convergence along the Indian–Eurasian collision zone. The ~NE dipping thrust faults have uplifted the young alluvial fan surfaces at the SW side of the KB. This created a major tectono-geomorphic boundary along the entire strike length of the KB that is characterised by (1) a low relief with sediment-filled sluggish streams to the SE and (2) an uplifted region, with actively flowing streams to the SW. The overall tectono-geomorphic expression suggests that recent activity along these faults has tilted the entire Kashmir valley towards NE. Further, the Mw 7.6 earthquake, which struck Northern Pakistan and Kashmir on 8 October 2005, also suggests a similar strike and NE dipping fault plane, which could indicate that the KB fault is continuous over a distance of ~210 km and connects on the west with the Balakot Bagh fault. However, the geomorphic and the structural evidences of such a structure are not very apparent on the north-west, which thus suggest that it is not a contiguous structure with the Balakot Bagh fault. Therefore, it is more likely that the KB fault is an independent thrust, a possible ramp on the Main Himalayan Thrust, which has uplifting the SW portion of the KB and drowning everything to the NE (e.g. Madden et al. 2011). Furthermore, it seems very likely that the KB fault could be a right stepping segment of the Balakot Bagh fault, similar to Riasi Thrust, as proposed by Thakur et al. (2010). The earthquake magnitude is measured by estimating the fault rupture parameters (e.g. Wells and Coppersmith in Bull Seismol Soc Am 84:974–1002, 1994). Therefore, the total strike length of the mapped KB fault is ~120 km and by assuming a dip of 29° (Avouac et al. in Earth Planet Sci Lett 249:514–528, 2006) and a down-dip limit of 20 km, a Mw of 7.6 is possible on this fault.