The distribution of radiogenic heat production as a function of depth in the Sierra Nevada Batholith, California

Geochemical analyses and geobarometric determinations have been combined to create a depth vs. radiogenic heat production database for the Sierra Nevada batholith, California. This database shows that mean heat production values first increase, then decrease, with increasing depth. Heat production is 2 μW/m³ within the 3-km-thick volcanic pile at the top of the batholith, below which it increases to an average value of 3.5 μW/m³ at 5.5 km depth, then decreases to 0.5–1 μW/m³ at 15 km depth and remains at these values through the entire crust below 15 km. Below the crust, from depths of 40–125 km, the batholith's root and mantle wedge that coevolved beneath the batholith appears to have an average radiogenic heat production rate of 0.14 μW/m³. This is higher than the rates from most published xenolith studies, but reasonable given the presence of crustal components in the arc root assemblages. The pattern of radiogenic heat production interpreted from the depth vs. heat production database is not consistent with the downward-decreasing exponential distribution predicted from modeling of surface heat flow data. The interpreted distribution predicts a reasonable range of geothermal gradients and shows that essentially all of the present day surface heat flow from the Sierra Nevada could be generated within the 35 km thick crust. This requires a very low heat flux from the mantle, which is consistent with a model of cessation of Sierran magmatism during Laramide flat-slab subduction, followed by conductive cooling of the upper mantle for 70 m.y. The heat production variation with depth is principally due to large variations in uranium and thorium concentration; potassium is less variable in concentration within the Sierran crust, and produces relatively little of the heat in high heat production rocks. Because silica content is relatively constant through the upper 30 km of the Sierran batholith, while U, Th, and K concentrations are highly variable, radiogenic heat production does not vary directly with silica content.     

A reassessment of models for hydrocarbon generation in the Khibiny nepheline syenite complex, Kola Peninsula, Russia

Although hydrocarbon-bearing fluids have been known from the alkaline igneous rocks of the Khibiny intrusion for many years, their origin remains enigmatic. A recently proposed model of post-magmatic hydrocarbon (HC) generation through Fischer-Tropsch (FT) type reactions suggests the hydration of Fe-bearing phases and release of H₂ which reacts with magmatically derived CO₂ to form CH₄ and higher HCs. However, new petrographic, microthermometric, laser Raman, bulk gas and isotope data are presented and discussed in the context of previously published work in order to reassess models of HC generation. The gas phase is dominated by CH₄ with only minor proportions of higher hydrocarbons. No remnants of the proposed primary CO₂-rich fluid are found in the complex. The majority of the fluid inclusions are of secondary nature and trapped in healed microfractures. This indicates a high fluid flux after magma crystallisation. Entrapment conditions for fluid inclusions are 450–550 °C at 2.8–4.5 kbar. These temperatures are too high for hydrocarbon gas generation through the FT reaction. Chemical analyses of rims of Fe-rich phases suggest that they are not the result of alteration but instead represent changes in magma composition during crystallisation. Furthermore, there is no clear relationship between the presence of Fe-rich minerals and the abundance of fluid inclusion planes (FIPs) as reported elsewhere. δ¹³C values for methane range from − 22.4‰ to − 5.4‰, confirming a largely abiogenic origin for the gas. The presence of primary CH₄-dominated fluid inclusions and melt inclusions, which contain a methane-rich gas phase, indicates a magmatic origin of the HCs. An increase in methane content, together with a decrease in δ¹³C isotope values towards the intrusion margin suggests that magmatically derived abiogenic hydrocarbons may have mixed with biogenic hydrocarbons derived from the surrounding country rocks.     

A Decision Support System for the Evaluation of Eco-engineering Strategies for Slope Protection

A decision support system (DSS) has been developed to assist expert and non-expert users in the evaluation and selection of eco-engineering strategies for slope protection. This DSS combines a qualitative hazard assessment of erosion and mass movements with a detailed catalogue of eco-engineering strategies for slope protection of which the suitability is evaluated in relation to the data entered. The slope decision support system (SDSS) is a knowledge based DSS in which knowledge is stored in frames containing rules that can evaluate the available information for a project, stored as project specific information (PSI) in a data file. The advantages of such a system are that it accepts incomplete information and that the qualitative nature of the information does not instil the user with a sense of unjustified exactitude. By its multidisciplinary and progressive nature, the DSS will be of value during the initial stages of an eco-engineering project when data collection and the potential of different eco-engineering strategies are considered. The accent of the output of the DSS is on the application of eco-engineering strategies for slope protection as an environmentally-friendly solution aiding sustainable development. For its acceptance within the engineering community, the DSS needs to prove its predictive capacity. Therefore, its performance has been benchmarked against successful and unsuccessful cases of slope stabilisation using eco-engineering. The target audience and the areas of application of this DSS are reviewed and the strategies for further development in this area suggested.     

Microvariability of line profiles in the spectrum of the star ι Her

We present the results of a search for and analysis of line-profile variations in the spectrum of the star ι Her. The observations were acquired with the 1.8 m telescope of the Bohyunsan Optical Astronomy Observatory (Republic of Korea) in May–June 2004. We obtained 69 spectra of the star with signal-to-noise ratios ≈300 and a time resolution of 5–7 min. Profile variability was revealed for six lines of HI, HeI, and SiIII, in the central parts of the lines. The variability amplitude is ≈(1–2)% in units of the intensity of the adjacent continuum. Evidence was found for cyclic variations of the lines, with periods from ≈7^h to ≈2.9^d. We conclude that ι Her belongs to the group of slowly pulsating stars.     

Deep seismic reflection study over the Vindhyans of Rajasthan: Implications for geophysical setting of the basin

This paper presents results of high-resolution deep seismic reflection profiling of the Proterozoic Vindhyan basin of the Rajasthan area along the Chandli-Bundi-Kota-Kunjer profile. Seismic images have been used to estimate the thickness of Vindhyan strata as well as to understand the tectonic framework of the basin. The results are constrained by gravity, magnetic and magnetotelluric data. The study reveals gentle SE-dipping reflection bands representing the Vindhyan strata. The seismic sections depict gradual thickening of the Vindhyan succession towards southeast from Bundi. The velocities of the upper and lower Vindhyans are identified as 4.6-4.8 km/s and 5.1-5.3 km/s. The NW limit of the Vindhyan basin is demarcated by the Great Boundary Fault (GBF) that manifests as a 30 km wide NW dipping thrust fault extending to a depth of 30 km.     

Geochemistry of silicic magmas in the Macolod Corridor, SW Luzon, Philippines: evidence of distinct, mantle-derived, crustal sources for silicic magmas

Silicic volcanic deposits (>65 wt% SiO₂), which occur as domes, lavas and pyroclastic deposits, are relatively abundant in the Macolod Corridor, SW Luzon, Philippines. At Makiling stratovolcano, silicic domes occur along the margins of the volcano and are chemically similar to the silicic lavas that comprise part of the volcano. Pyroclastic flows are associated with the Laguna de Bay Caldera and these are chemically distinct from the domes and lavas at Makiling stratovolcano. As a whole, samples from the Laguna de Bay Caldera contain lower concentrations of MgO and higher concentrations of Fe₂O_3(t) than the samples from domes and lavas. The Laguna de Bay samples are more enriched in incompatible trace elements. The silicic rocks from the domes, Makiling Volcano and Laguna de Bay Caldera all contain high alkalis and high K₂O/Na₂O ratios. Melting experiments of primitive basalts and andesites demonstrate that it is difficult to produce high K₂O/Na₂O silicic magmas by fractional crystallization or partial melting of a low K₂O/Na₂O source. However, recent melting experiments (Sisson et al., Contrib Mineral Petrol 148:635–661, 2005) demonstrate that extreme fractional crystallization or partial melting of K-rich basalts can produce these silicic magmas. Our model for the generation of the silicic magmas in the Macolod Corridor requires partial melting of mantle-derived, evolved, moderate to K-rich, crystallized calc-alkaline magmas that ponded and crystallized in the mid-crust. Major and trace element variations, along with oxygen isotopes and ages of the deposits, are consistent with this model. Electronic Supplementary Material Supplementary material is available for this article at     

Atomistic calculations of structural and elastic properties of serpentine minerals: the case of lizardite

The physical properties of the hydrous phyllosilicate lizardite have been investigated by atomistic simulation using the GULP code based on transferable semi-empirical interatomic potentials. Lizardite behavior was first investigated during structure relaxation at room temperature. The Helmholtz free energy is minimum for an equilibrium structure that is in agreement with experiment. The bulk, shear, and Young modulii for lizardite were calculated along with the Poisson ratio. From the shear and bulk modulii, we also calculated translational and longitudinal acoustic wave velocities that are important quantities for tectonophysics models. As expected, lizardite is stiffer in the a direction parallel to the layers than in the c perpendicular direction; the variation of the unit cell parameters with pressure is in good agreement with experiment. The cohesive energy between two successive layers along c direction was calculated at 0.33 eV (i.e., 0.11 eV per OH bond) in good agreement with recent ab initio calculations. Upon pressure and temperature variations, we evidenced that structural changes are mainly pressure induced; pressure being accommodated by a decrease of the c parameter up to 10 GPa. We also found that the change of slope in the derivative of the c cell parameter with respect to pressure occurring around 2 GPa originates from the bending of the interlayer hydroxyl groups with respect to the layer normal direction.