CyberShake: A Physics-Based Seismic Hazard Model for Southern California

CyberShake, as part of the Southern California Earthquake Center??s (SCEC) Community Modeling Environment, is developing a methodology that explicitly incorporates deterministic source and wave propagation effects within seismic hazard calculations through the use of physics-based 3D ground motion simulations. To calculate a waveform-based seismic hazard estimate for a site of interest, we begin with Uniform California Earthquake Rupture Forecast, Version 2.0 (UCERF2.0) and identify all ruptures within 200?km of the site of interest. We convert the UCERF2.0 rupture definition into multiple rupture variations with differing hypocenter locations and slip distributions, resulting in about 415,000 rupture variations per site. Strain Green Tensors are calculated for the site of interest using the SCEC Community Velocity Model, Version 4 (CVM4), and then, using reciprocity, we calculate synthetic seismograms for each rupture variation. Peak intensity measures are then extracted from these synthetics and combined with the original rupture probabilities to produce probabilistic seismic hazard curves for the site. Being explicitly site-based, CyberShake directly samples the ground motion variability at that site over many earthquake cycles (i.e., rupture scenarios) and alleviates the need for the ergodic assumption that is implicitly included in traditional empirically based calculations. Thus far, we have simulated ruptures at over 200 sites in the Los Angeles region for ground shaking periods of 2?s and longer, providing the basis for the first generation CyberShake hazard maps. Our results indicate that the combination of rupture directivity and basin response effects can lead to an increase in the hazard level for some sites, relative to that given by a conventional Ground Motion Prediction Equation (GMPE). Additionally, and perhaps more importantly, we find that the physics-based hazard results are much more sensitive to the assumed magnitude-area relations and magnitude uncertainty estimates used in the definition of the ruptures than is found in the traditional GMPE approach. This reinforces the need for continued development of a better understanding of earthquake source characterization and the constitutive relations that govern the earthquake rupture process.     

Adsorption of Rb and Sr at the orthoclase (0 0 1)-solution interface

Adsorption of Rb⁺ and Sr²⁺ at the orthoclase (0 0 1)-solution interface is probed with high-resolution X-ray reflectivity and resonant anomalous X-ray reflectivity. Specular X-ray reflectivity data for orthoclase in contact with 0.01 m RbCl solution at pH 5.5 reveal a systematic increase in electron density adjacent to the mineral surface with respect to that observed in contact with de-ionized water (DIW). Quantitative analysis indicates that Rb⁺ adsorbs at a height of 0.83 ± 0.03 Å with respect to the bulk K⁺ site with a nominal coverage of 0.72 ± 0.10 ions per surface unit mesh (55.7 Å²). These results are consistent with an ion-exchange reaction in which Rb⁺ occupies an inner-sphere adsorption (IS) site. In contrast, X-ray reflectivity data for orthoclase in contact with 0.01 m Sr(NO₃)₂ solution at pH 5.3 reveal few significant changes with respect to DIW. Resonant anomalous X-ray reflectivity was used to probe Sr²⁺ adsorption and to image its vertical distribution. This element-specific measurement reveals that Sr²⁺ adsorbs with a total coverage of 0.37 ± 0.02 ions per surface unit mesh, at a substantially larger height (3.28 ± 0.05 Å) than found for Rb⁺, and with a relatively broad density distribution (having a root-mean-square width of 1.88 ± 0.08 Å for a single-peak model), implying that Sr²⁺ adsorbs primarily as a fully-hydrated outer-sphere (OS), species. Comparison to a two-height model suggests that 13 ± 5% of the adsorbed Sr²⁺ may be present as an IS species. This partitioning implies a ∼5 kJ/mol difference in free energy between the IS and OS Sr²⁺ on orthoclase. Differences in the partitioning of Sr²⁺ between IS and OS species for orthoclase (0 0 1) and muscovite (0 0 1) suggest control by the geometry of the IS adsorption site. Results for the OS distribution are compared to predictions of the Poisson-Boltzmann equation in the strong coupling regime, which predicts an intrinsically narrow vertical diffuse ion distribution; the OS distribution might thus be thought of as the diffuse ion profile in the limit of high surface charge.     

The speciation of water in silicate melts

Previous models of water solubility in silicate melts generally assume essentially complete reaction of water molecules to hydroxyl groups. In this paper a new model is proposed that is based on the hypothesis that the observed concentrations of molecular water and hydroxyl groups in hydrous silicate glasses reflect those of the melts from which they were quenched. The new model relates the proportions of molecular water and hydroxyl groups in melts via the following reaction describing the homogeneous equilibrium between melt species: H₂O_molecular (melt) + oxygen (melt) = 2OH (melt). An equilibrium constant has been formulated for this reaction and species are assumed to mix ideally. Given an equilibrium constant for this reaction of 0.1–0.3, the proposed model can account for variations in the concentrations of molecular water and hydroxyl groups in melts as functions of the total dissolved water content that are similar to those observed in glasses. The solubility of molecular water in melt is described by the following reaction: H₂O (vapor) = H₂O_molecular (melt).These reactions describing the homogeneous and heterogeneous equilibria of hydrous silicate melts can account for the following observations: the linearity between f_H2O and the square of the mole fraction of dissolved water at low total water contents and deviations from linearity at high total water contents; the difference between the partial molar volume of water in melts at low total water contents and at high total water contents; the similarity between water contents of vapor-saturated melts of significantly different compositions at high pressures versus the dependence on melt composition of water solubility in silicate melts at low pressures; and the variations of viscosity, electrical conductivity, the diffusivity of “water,” the diffusivity of cesium, and phase relationships with the total dissolved water contents of melts.This model is thus consistent with available observations on hydrous melt systems and available data on the species concentrations of hydrous glasses and is easily tested, since measurements of the concentrations of molecular water and hydroxyl groups in silicate glasses quenched from melts equilibrated over a range of conditions and total dissolved water contents are readily obtainable.     

Galileiite: A new meteoritic phosphate minera

Abstract— A new mineral named galileiite, NaFe₄(PO₄)₃, has been found within troilite nodules in iron meteorites of the IIIA and IIIB groups. the mineral is optically positive (ω = 1.72, ω = 1.75), colorless in transmitted light and pale amber in reflected light. Grains of galileiite are very small, generally 10 μm or less; rarely, grains are up to 30 μm. It is associated with Ca-free graftonite (or Ca-free sarcopside), chromite and, occasionally, schreibersite. Johnsomervilleite may occur within troilite nodules in the same meteorite as galileiite, but they have never been observed together in the same troilite nodule. Because of the small sample size, single crystal x-ray work was not successful; however, Gandolfi diffraction measurements were made. The three strongest diffraction peaks are 2.71 Å, 3.01 Å and 4.13 Å. On the basis of its composition and similar diffraction pattern, it is considered to be related to johnsomervilleite, fillowite and chladniite, all of which are rhombohedral and isostructural. Galileiite may also be rhombohedral, but that is yet to be demonstrated.     

Geology,fluid inclusion and sulphur isotope characteristics of the El Cobre VHMS deposit,Southern Cuba

The El Cobre deposit is located in eastern Cuba within the volcanosedimentary sequence of the Sierra Maestra Paleogene arc. The deposit is hosted by tholeiitic basalts, andesites and tuffs and comprises thick stratiform barite and anhydrite bodies, three stratabound disseminated up to massive sulphide bodies produced by silicification and sulphidation of limestones or sulphates, an anhydrite stockwork and a siliceous stockwork, grading downwards to quartz veins. Sulphides are mainly pyrite, chalcopyrite and sphalerite; gold occurs in the stratabound ores. Fluid inclusions measured in sphalerite, quartz, anhydrite and calcite show salinities between 2.3 and 5.7 wt% NaCl eq. and homogenisation temperatures between 177 and 300°C. Sulphides from the stratabound mineralisation display δ ³⁴S values of 0‰ to +6.0‰, whilst those from the feeder zone lie between −1.4‰ and +7.3‰. Sulphides show an intra-grain sulphur isotope zonation of about 2‰; usually, δ ³⁴S values increase towards the rims. Sulphate sulphur has δ ³⁴S in the range of +17‰ to +21‰, except two samples with values of +5.9‰ and +7.7‰. Sulphur isotope data indicate that the thermochemical reduction of sulphate from a hydrothermal fluid of seawater origin was the main source of sulphide sulphur and that most of the sulphates precipitated by heating of seawater. The structure of the deposit, mineralogy, fluid inclusion and isotope data suggest that the deposit formed from seawater-derived fluids with probably minor supply of magmatic fluids.     

Hydrochemistry of urban groundwater in Seoul, South Korea: effects of land-use and pollutant recharge

The ionic and isotopic compositions (δD, δ¹⁸O, and ³H) of urban groundwaters have been monitored in Seoul to examine the water quality in relation to land-use. High tritium contents (6.1–12.0 TU) and the absence of spatial/seasonal change of O–H isotope data indicate that groundwaters are well mixed within aquifers with recently recharged waters of high contamination susceptibility. Statistical analyses show a spatial variation of major ions in relation to land-use type. The major ion concentrations tend to increase with anthropogenic contamination, due to the local pollutants recharge. The TDS concentration appears to be a useful contamination indicator, as it generally increases by the order of forested green zone (average 151 mg/l), agricultural area, residential area, traffic area, and industrialized area (average 585 mg/l). With the increased anthropogenic contamination, the groundwater chemistry changes from a Ca–HCO₃ type toward a Ca–Cl(+NO₃) type. The source and behavior of major ions are discussed and the hydrochemical backgrounds are proposed as the basis of a groundwater management plan.     

中国有壳变形虫Argynnia(Vucetich,1974)属的修订及若干环境重建问题探讨

有壳变形虫是一类具外壳的陆相淡水根足纲(Rhizopoda)原生动物,广泛栖息于湖泊、泥炭、沼泽、土壤等各种淡水潮湿环境,由于其壳体具有很好的抗腐蚀性而在沉积物中保存下来,近年来广泛应用于泥炭和湖泊湿地的古环境重建研究中.大部分有壳变形虫是世界广泛分布的,但是也有一些种类具有一定的地带性分布格局.形态学和分子生物学的研...     

Practically Applicable Central Limit Theorem for Spatial Statistics

Let {Z(s):s∈D⊆ℝ ^d } be a zero mean stationary random field observed at a finite number of locations. Lahiri (Sankhya Ser. A 65:356–388, 2003) proved spatial central limit theorems (CLT) for ∑ _i=1 ⁿ Z(s _i ) assuming a ‘nearly infill domain sampling’. Applications of his results depended on the underlying spatial sampling region and the design in a complicated fashion. The main objective of this paper is to provide CLTs that could be applied easily in practice. We present two main results assuming a ‘nearly infill domain sampling’ defined mainly in terms of dependence. Theorem 1 establishes a CLT for ∑ _i=1 ⁿ Z(s _i ) and Theorem 2 is obtained mainly for applications to density estimates. We report on a simulation study for illustrating a way of applying our results in practice.