Effect of Matric Suction on Resilient Modulus of Compacted Aggregate Base Courses

This research was conducted to investigate the effect of matric suction on resilient modulus of unbound aggregate base courses. The study characterized the water characteristic curves and resilient modulus versus matric suction relationships of aggregate base courses that were compacted at different water contents and between 98 and 103 % of the modified Proctor density. The soil–water characteristic curve (SWCC) and the relationship between resilient modulus (M _r ) and matric suction (ψ) were established for different unbound granular and recycled asphalt pavement materials. This relationship is important for predicting changes in modulus due to changes in moisture of unbound pavement materials. Resilient modulus tests were conducted according to the National Cooperative Highway Research Program (NCHRP) 1-28A procedure at varying water contents, and the measured SWCC was used to determine the corresponding matric suction. Three reference summary resilient moduli (SRM) were considered: at optimum water content, optimum water content +2 % and optimum water content ?2 %. The Bandia and Bargny limestones are characterized by a higher water-holding capacity which explains why the modulus of limestone was more sensitive to water content than for basalt or quartzite. Limestones tend to be more sensitive to changes in water content and thus to matric suction. The shape of the SWCC depends on the particle size distribution and the cementation properties from dehydration of the aggregates. Material properties required as input to the Mechanistic-Empirical Pavement Design Guide (M-EPDG) to predict changes in resilient modulus in response to changes in moisture contents in the field were determined for implementation in the M-EPDG process. Results show that the SRM was more correlated with matric suction than with compaction water content (for resilient modulus testing). The empirical models commonly used to predict the SWCC such as the Perera et al. (Prediction of the SWCC based on grain-size-distribution and index properties. GSP 130 Advances in Pavement Engineering, ASCE, 2005) and the M-EPDG (NCHRP in Guide for mechanistic-empirical design of pavement structures. National cooperative highway research program. ARA, Inc., ERES Consultants Division, Champaign, IL, 2004) models tend to underestimate the SWCC and cannot provide reasonable estimation. SRM normalized with respect to the SRM at the optimum water content varied linearly with the logarithm of matric suction. Empirical relationships between SRM and matric suction on semi-logarithmic scale were established and are reported.     

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.     

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.     

Design of Raft Foundations for High-Rise Buildings on Jointed Rock

This paper presents calculations of displacements and bending moments in a 2-m-thick reinforced-concrete foundation slab using three-dimensional finite-element software. A preliminary paper was presented by Justo et al. (Rock Mech Rock Eng 43:287–304, 2010). The slab is the base of a tower of 137 m height above foundation, supported on jointed and partly weathered basalt and scoria. Installation of rod extensometers at different depths below foundation allowed comparison between measured displacements and displacements calculated using moduli obtained from rock classification systems and three material models: elastic, Mohr–Coulomb and hardening (H). Although all three material models can provide acceptable results, the H model is preferable when there are unloading processes. Acceptable values of settlement may be achieved with medium meshing and an approximate distribution of loads. The absolute values of negative bending moments (tensions below) increase as the rock mass modulus decreases or when the mesh is refined. The paper stresses the importance of adequately representing the details of the distribution of loads and the necessity for fine meshing to obtain acceptable values of bending moments.     

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.     

Geology and geochemistry of Tertiary basalt in south Wadi Hodein area, South Eastern Desert, Egypt

Tertiary basalt is widespread in the area south of Wadi Hodein, south Eastern Desert, Egypt. It is the youngest unit in the basement rocks of the Central Eastern Desert classification of El Shazly (Proc 22nd Intl Geol Congr, New Delhi 10:88–101, 1964) and El Ramly (Ann Geol Surv Egypt II:1–17, 1972), traversed all the previous succession of the basement rocks as well as the Nubia Sandstone of Cretaceous age, forming sheets, small hills, ridges, and dikes. This Tertiary basalt is strongly associated with the opening of the Red Sea. Geologic, petrographic, and petrochemical studies as well as microprobe and X-ray analyses were performed on samples from Wadi Hodein Tertiary basalt. Field and petrographic studies classified the Tertiary basalt in south Wadi Hodein into porphyritic olivine basalt, plagiophyric basalt, and doleritic basalt. Opaque minerals (magnetite and ilmenite) constitute 6–7.5% of this basalt. Petrochemical studies and microprobe analyses reveal that they are low-TiO₂ basalt with low uranium and thorium contents, classified as being basaltic andesite to andesite, originated from calc-alkaline magma, and developed in within-plate tectonic environment. Scanning electron microscopy shows that magnetite and ilmenite are the prevalent opaque minerals in this Tertiary basalt. Field radiometric measurements of the Tertiary basalt in south Wadi Hodein reveals low uranium and thorium contents. Uranium contents range from 0.5 to 0.9 ppm, while thorium contents range from 1.2 to 3.2 ppm. Fractional crystallization and mass balance modeling indicate that the most-silica low-TiO₂ Tertiary basalt in south Wadi Hodein can be derived from the relatively less-silica low-TiO₂ Tertiary basalt of south Quseir and Gabal Qatrani through fractional crystallization of plagioclase, olivine, augite, and titanomagnetite oxides. Tertiary basalts in south Wadi Hodein and south Quseir have nearly the same age, 25 Ma (Sherif, The Fifth International Conference on the Geology of Africa, 2007), 24 Ma (Meneisy and Abdel Aal, Ain Shams Sci Bull 25(24B): 163–176, 1984), and 27 Ma (El Shazly et al., Egypt J Geol 1975), respectively. Finally, the fractionation modeling and geochemical characteristics of these basalts suggested their origination from one basaltic magma emplaced in late Oligocene.     

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}}}}$ .     

考虑吸力效应的非饱和黏土回弹模量预估模型

回弹模量是反映路基弹性支撑能力的重要指标,并受运营期间基质吸力变化的显著影响。通过3种压实黏土试样的土-水特征曲线试验和重复动三轴试验,研究了回弹模量与基质吸力的相关性,并基于非饱和土有效应力原理,提出了考虑吸力效应的回弹模量预估模型。结果表明：基质吸力与回弹模量之间具有较强的非线性相关性,可采用指数形式进行两者的拟合;含水率由最佳含水率-4%提升至最佳含水率+4%时,受基质吸力减小的影响,回弹模量降幅达到29.1%～39.0%;在 土-水特征曲线过渡区,采用以饱和度表征的有效应力参数,可较好地反映基质吸力对于回弹模量的贡献,预估效果优于采用残余含水率或进气值相关的有效应力参数;在相同饱和度的情况下,黏粒含量越多、塑性指数越高,基质吸力对有效应力的贡献比例越小。     

Extreme Learning Machine Based Modeling of Resilient Modulus of Subgrade Soils

This paper investigates the potential of two variants of extreme learning machine based regression approaches in predicting the resilient modulus of cohesive soils. Support vector regression was used to compare the performance of the proposed extreme learning machine based regression approaches. The dataset used in this study was derived from literature and consists of 9 input parameters with a total of 891 cases. For testing, two methods i.e. train/test and tenfold cross validation was used. In case of train and test methods, a total of 594 randomly selected cases were used to train different algorithms and the remaining 297 data were used to test the created models. Correlation coefficient value of 0.991 (root mean square error = 3.47 MPa) was achieved by polynomial kernel based extreme learning machine in comparison to 0.990 and 0.990 (root mean square error = 4.790 and 4.290 MPa) by simple extreme learning machine and radial basis kernel function based support vector regression respectively with test dataset. Comparisons of results with tenfold cross validation also suggest that polynomial kernel based extreme learning machine works well in terms of root mean square error and computational cost with the used dataset. Sensitivity analysis suggests the importance of confining stress and deviator stress in predicting the resilient modulus when using with polynomial kernel based extreme learning machine modeling approach.     

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.     

Effects of water and fines contents on the resilient modulus of the interlayer soil of railway substructure

This paper deals with the resilient behavior of the interlayer soil which is created mainly by the interpenetration of ballast and subgrade soils. The interlayer soil studied was taken from a site in the southeast of France. Large-scale cyclic triaxial tests were carried out at three water contents (w = 4, 6 and 12 %) and three fines contents corresponding to 5, 10 % subgrade added to the natural interlayer soil and 10 % fine particles (<80 μm) removed from the natural interlayer soil. Soil specimens underwent various deviator stresses, and for each deviator stress, a large number of cycles was applied. The effects of deviator stress, number of cycles, water content and fines content on the resilient modulus (M _r) were analyzed. It appears that the effects of water content and fines content must be analyzed together because the two effects are closely linked. Under unsaturated conditions, the soil containing high fines content has higher resilient modulus due to the contribution of suction. When the soil approaches the saturated state, it loses its mechanical enhancement with a sharp decrease in resilient modulus.     

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.     

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.     

Use of site amplification and anelastic attenuation for the determination of source parameters of the Sikkim earthquake of September 18, 2011, using far-field strong-motion data

An earthquake of magnitude 6.9 (M _w) occurred in the Sikkim region of India on September 18, 2011. This earthquake is recorded on strong-motion network in Uttarakhand Himalaya located about 900 km away from the epicenter of this earthquake. In this paper acceleration record from six far-field stations has been used to compute the source parameters of this earthquake. The acceleration spectra of ground motion at these far-field stations are strongly affected by both local site effects and near-site anelastic attenuation. In the present work the spectrum of S-phase recorded at these far-field stations has been corrected for anelastic attenuation at both source and site and the site amplification terms. Site amplifications at different stations and near-site shear wave attenuation factor have been computed by the technique of inversion of acceleration spectra given by Joshi et al. (Pure Appl Geophys 169:1821–1845, 2012a). For estimation of site amplification and shear wave quality factor [Q _β (f)] at the recording sites, ten local events recorded at various stations between July 2011 and December 2011 have been used. The obtained source spectrum from acceleration records is compared with the theoretical source spectrum defined by Brune (J Geophys Res 76:5002, 1970) at each station for both horizontal components of the records. Iterative forward modeling of theoretical source spectrum gives the average estimate of seismic moment (M _o), source radius (r _o) and stress drop (Δσ) as (3.2 ± 0.8) × 10²⁶ dyne cm, 13.3 ± 0.8 km and 59.2 ± 8.8 bars, respectively, for the Sikkim earthquake of September 18, 2011.     

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.     

Crystallization of oxidized, moderately hydrous arc basalt at mid- to lower-crustal pressures: implications for andesite genesis

This study focuses on the production of convergent margin calc-alkaline andesites by crystallization–differentiation of basaltic magmas in the lower to middle crust. Previous experimental studies show that dry, reduced, subalkaline basalts differentiate to tholeiitic (high Fe/Mg) daughter liquids, but the influences of H₂O and oxidation on differentiation are less well established. Accordingly, we performed crystallization experiments at controlled oxidized fO₂ (Re–ReO₂ ≈ ΔNi–NiO + 2) on a relatively magnesian basalt (8.7 wt% MgO) typical of mafic magmas erupted in the Cascades near Mount Rainier, Washington. The basalt was synthesized with 2 wt% H₂O and run at 900, 700, and 400 MPa and 1,200 to 950 °C. A broadly clinopyroxenitic crystallization interval dominates near the liquidus at 900 and 700 MPa, consisting of augite + olivine + orthopyroxene + Cr-spinel (in decreasing abundance). With decreasing temperature, plagioclase crystallizes, Fe–Ti-oxide replaces spinel, olivine dissolves, and finally amphibole appears, producing gabbroic and then amphibole gabbroic crystallization stages. Enhanced plagioclase stability at lower pressure narrows the clinopyroxenitic interval and brings the gabbroic interval toward the liquidus. Liquids at 900 MPa track along Miyashiro’s (Am J Sci 274(4):321–355, 1974) tholeiitic versus calc-alkaline boundary, whereas those at 700 and 400 MPa become calc-alkaline at silica contents ≥56 wt%. This difference is chiefly due to higher temperature appearance of magnetite (versus spinel) at lower pressures. Although the evolved liquids are similar in many respects to common calc-alkaline andesites, the 900 and 700 MPa liquids differ in having low CaO concentrations due to early and abundant crystallization of augite, with the result that those liquids become peraluminous (ASI: molar Al/(Na + K + 2Ca) > 1) at ≥61 wt% SiO₂, similar to liquids reported in other studies of the high-pressure crystallization of hydrous basalts (Müntener and Ulmer in Geophys Res Lett 33(21):L21308, 2006). The lower-pressure liquids (400 MPa) have this same trait, but to a lesser extent due to more abundant near-liquidus plagioclase crystallization. A compilation of >6,500 analyses of igneous rocks from the Cascades and the Sierra Nevada batholith, representative of convergent margin (arc) magmas, shows that ASI increases continuously and linearly with SiO₂ from basalts to rhyolites or granites and that arc magmas are not commonly peraluminous until SiO₂ exceeds 69 wt%. These relations are consistent with plagioclase accompanying mafic silicates over nearly all the range of crystallization (or remelting). The scarcity of natural peraluminous andesites shows that progressive crystallization–differentiation of primitive basalts in the deep crust, producing early clinopyroxenitic cumulates and evolved liquids, does not dominate the creation of intermediate arc magmas or of the continental crust. Instead, mid- to upper-crustal differentiation and/or open-system processes are critical to the production of intermediate arc magmas. Primary among the open-system processes may be extraction of highly evolved (granitic, rhyolitic) liquids at advanced degrees of basalt solidification (or incipient partial melting of predecessor gabbroic intrusions) and mixing of such liquids into replenishing basalts. Furthermore, if the andesitic-composition continents derived from basaltic sources, the arc ASI–SiO₂ relation shows that the mafic component returned to the mantle was gabbroic in composition, not pyroxenitic.     

Effect of freeze—thaw cycles on resilient properties of fine-grained soils

Stress-deformation data for silt and clay subgrade soils were obtained from in-situ tests and laboratory tests, for use in mechanistic models for design of pavements that will experience freezing and thawing. Plate-bearing tests were run on in-service allbituminous-concrete (ABC) pavements constructed directly on silt subgrade, and on an experimental ABC pavement constructed on clay subgrade, applying repeated loads to the pavement surfaces while the subgrade was frozen, thawing, thawed, and fully recovered. The in-service pavement had experienced several seasons of natural cyclic freezing and thawing, while the experimental pavement was artificially frozen and thawed twice. Repeated-load laboratory triaxial compression tests were performed on the same soils in the frozen and thawed states.

Analysis of deflection data from the in-situ tests showed resilient moduli of the subgrade soils up to more than 10 GPa when frozen, as low as 2 MPa during the thawing period, and up to more than 100 MPa when fully recovered. Analysis of the laboratory tests, which gave moduli comparable to the latter values, showed that resilient modulus and Poisson's ratio in the thawed and recovering conditions can be expressed as a function of the stress state, the moisture content, and the dry density.     

Solution of the MoMaS reactive transport benchmark with MIN3P—model formulation and simulation results

This paper summarizes the governing equations as implemented in the MIN3P multicomponent flow and reactive transport code (Mayer et al., Water Resour Res 38:1174, 2002) and introduces the equations in discretized form. Linearization and solution methods are presented including adaptive time stepping and update modification schemes. Code-specific details for the implementation of the GdR MoMaS benchmark simulations (Carrayrou et al., Comput Geosci, 2009) are presented. The standard version of the MIN3P code was used to solve the Easy, Medium, and Hard Test Cases, in one and two spatial dimensions, for both advection- and diffusion-dominated conditions. An analysis of the sensitivity of the solution in relation to spatial and temporal discretization parameters is provided for the Easy Test Case, selected results are presented for the Medium and Hard Test Cases, and the performance of the code as a function of discretization parameters is evaluated for all test cases.     

Principal component analysis method in the detection of total electron content anomalies in the 24 hrs prior to large earthquakes

The goal of this research is to examine one-dimensional total electron content (TEC) data using principal component analysis (PCA) to search for total electron content (TEC) anomalies associated with large earthquakes in 24 h prior to nucleation. The characteristics of principal eigenvalues generated for TEC prior to 24 earthquakes of magnitude scale M?≥?5.0 and 6 lesser earthquakes of magnitude scale M?<?5 that occurred in Taiwan from 01 January 2000 to 31 December 2001 are examined. In an earlier paper, I was able to confirm the statistical findings of Liu et al. (J Geophys Res 111, 2006) that sparse earthquake-associated TEC anomalies existed in 5 days prior to the 12 large earthquakes they examined (Lin, Terr Atm Ocean Sci, 2010). In this paper, I wish to examine the subtlety of principal component analysis in detecting earthquake-associated TEC anomalies by examining if such precursors can be detected in 24 h prior to large earthquakes. Of the earthquakes examined, TEC anomalies given by clear extreme principal eigenvalues were evident within 24 h of nucleation for 21 of the 24 earthquakes of M?≥?5.0. After making allowance for the general status of background TEC, it is clear that these extreme principal eigenvalues are representative of earthquake-associated anomalies. For the smaller earthquakes (M?<?5), it was not possible to differentiate earthquake-associated anomalies from background effects on TEC status. These new findings confirm the validity of PCA in searching for earthquake-associated TEC anomalies and show that it is subtle enough to detect TEC anomalies within 24 h leading to a large earthquake. If this approach continues to prove successful, it could theoretically be used in real-time prediction of large earthquakes through early detection of earthquake-associated TEC anomalies.