青藏高原北羌塘半岛湖新生代粗面玄武岩橄榄石电子探针和激光探针分析

利用电子探针和激光探针剥蚀系统（LA-ICP-MS），对北羌塘新第三纪粗面玄武岩中的橄榄石主元素和微量、稀土元素进行了系统分析。结果表明，本区橄榄石Fo平均为88，属贵橄榄石种属。相对富集Ni、Co和重稀土，而强烈亏损轻稀土及Rb、Sr、Ba、Zr等大离子亲石元素。其稀土元素配分型配分型式与粗面玄武岩全岩稀土配分型式呈互补状态。     

Use of chemical tracing to study acquisition modalities of the mineralization and behaviour of unconfined groundwaters under a semi-arid climate: the case of the Souss plain (Morocco)

The chemical content of the Souss unconfined groundwater displays spatial variations in conductivity (between 400 and 6,000 µS cm^-1). The chemical tracers (Cl^-, SO₄^2-, Sr²⁺, Br^-), which characterize the different components of the groundwater, allowed the determination of the origin of water salinity. Cl^- and SO₄^2-, reaching respectively 2,000 and 1,650 mg L^-1, display localized salinity anomalies. Br^-/Cl^- ratio distinguishes marine-influenced impoverished zones versus the oceanic domain. Thus, salinity anomalies can be attributed: (1) downstream, to a currently existing salt-encroachment (with added waste water) and sedimentary palaeosalinity, (2) in the middle-Souss, to High Atlas evaporites and to irrigation water recycling. Sr²⁺/Ca²⁺ ratio (evaporites if >1‰), confirms the evaporitic origin of the anomalies along the right bank of oued Souss. Furthermore, it facilitates the distinction between the different aquifer contributions (Cretaceous, Jurassic and Triassic), and it highlights leakage from deep Turonian limestones in the groundwater recharge system. To the south, recharge is from the Anti Atlas (evaporite-free) waters. Oxygen-18 measurements confirm the groundwater recharge from the High and Anti Atlas as piezometric maps and chemical tracers suggested, plus from leakage from the Turonian and the marine aquifers.     

Combination of Dependent Realizations Within the Gradual Deformation Method

Gradual deformation is a parameterization method that reduces considerably the unknown parameter space of stochastic models. This method can be used in an iterative optimization procedure for constraining stochastic simulations to data that are complex, nonanalytical functions of the simulated variables. This method is based on the fact that linear combinations of multi-Gaussian random functions remain multi-Gaussian random functions. During the past few years, we developed the gradual deformation method by combining independent realizations. This paper investigates another alternative: the combination of dependent realizations. One of our motivations for combining dependent realizations was to improve the numerical stability of the gradual deformation method. Because of limitations both in the size of simulation grids and in the precision of simulation algorithms, numerical realizations of a stochastic model are never perfectly independent. It was shown that the accumulation of very small dependence between realizations might result in significant structural drift from the initial stochastic model. From the combination of random functions whose covariance and cross-covariance are proportional to each other, we derived a new formulation of the gradual deformation method that can explicitly take into account the numerical dependence between realizations. This new formulation allows us to reduce the structural deterioration during the iterative optimization. The problem of combining dependent realizations also arises when deforming conditional realizations of a stochastic model. As opposed to the combination of independent realizations, combining conditional realizations avoids the additional conditioning step during the optimization process. However, this procedure is limited to global deformations with fixed structural parameters.     

The Petrology and Geochemistry of Calc-Alkaline Andesites on Shodo-Shima Island, SW Japan

Petrographical and geochemical characteristics of calc-alkalineandesites on Shodo-Shima Island, SW Japan, having bulk compositionslargely identical to the continental crust, are presented. Thefollowing petrographic observations suggest a role for magmamixing in producing such andesite magmas: (1) two types of olivinephenocrysts and spinel inclusions, one with compositions identicalto those in high-Mg andesites and the other identical to thosein basalts, are recognized in terms of Ni–Mg and Cr–Al–Fe³⁺relations, respectively; (2) the presence of orthopyroxene phenocrystswith mg-number >90 suggests the contribution of an orthopyroxene-bearinghigh-Mg andesite magma to production of calc-alkaline andesites;(3) reversely zoned pyroxene phenocrysts may not be in equilibriumwith Mg-rich olivine, suggesting the involvement of a differentiatedandesite magma as an endmember component; (4) the presence ofvery Fe-rich orthopyroxene phenocrysts indicates the associationof an orthopyroxene-bearing rhyolitic magma. Contributions fromthe above at least five endmember magmas to the calc-alkalineandesite genesis can also provide a reasonable explanation ofthe Pb–Sr–Nd isotope compositions of such andesites. KEY WORDS: calc-alkaline andesites; high-Mg andesites; magma mixing; continental crust; SW Japan     

A scanning electron microscope study of the effects of dynamic recrystallization on lattice preferred orientation in olivine

Effects of dynamic recrystallization on lattice preferred orientation (LPO) in olivine were investigated through the combination of two SEM-based techniques, electron backscattered diffraction (EBSD) technique for crystallographic orientation measurement and backscattered electron imaging (BEI) for dislocation observation. Samples are experimentally deformed olivine aggregates in simple shear geometry. In the sample deformed at T=1473 K and high stresses (480 MPa), only incipient dynamic recrystallization is observed along grain-boundaries. Orientations of these small recrystallized grains are more random than that of relict grains, suggesting an important role of grain-boundary sliding at this stage of recrystallization. In the sample deformed at T=1573 K and low stress (160 MPa), dynamic recrystallization is nearly complete and the LPO is characterized by two [100] peaks. One peak is located at the orientation subparallel to the shear direction and is dominated by grains with high Schmid factor. The other occurs at high angles to the shear direction and is due to the contribution from grains with low Schmid factor. Grains with high Schmid factor tend to have higher dislocation densities than those with low Schmid factor. Based on these observations, we identify two mechanisms by which dynamic recrystallization affects LPO: (1) enhancement of grain-boundary sliding due to grain-size reduction, leading to the modification of LPO caused by the relaxation of constraint for deformation; (2) grain-boundary migration by which grains with lower dislocation densities grow at the expense of grains with higher dislocation densities. Based on the deformation mechanism maps and stress versus recrystallized grain-size relation, we suggest that the first mechanism always plays an important role whereas the second mechanism has an important effect only under limited conditions.     

Rates of Deglaciation during the Last Glaciation and Holocene in the Cordillera Vilcanota-Quelccaya Ice Cap Region, Southeastern Perú

Moraine chronology is combined with digital topography to model deglacial rates of paleoglacier volumes in both the Huancané Valley on the west side of the Quelccaya Ice Cap and the Upismayo Valley on the northwest side of the Cordillera Vilcanota. The fastest rates of deglaciation (39×10⁻⁵ to 114×10⁻⁵ km³ yr⁻¹ and 112×10⁻⁵ to 247×10⁻⁵ km³ yr⁻¹ for each valley, respectively) were calculated for the most recent paleoglaciers, corresponding to the last few centuries. These results are consistent with observations in the Venezuelan Andes showing high rates of deglaciation since the Little Ice Age. These rates also fall within the range of 20th century rates of deglaciation measured on the Quelccaya Ice Cap (29×10⁻⁵ to 220×10⁻⁵ km³ yr⁻¹, Brecher and Thompson, 1993; Thompson, 2000). These results imply that rates of deglaciation may fluctuate significantly over time and that high rates of deglaciation may not be exclusive to the late 20th century. Equilibrium line altitude (ELA) depressions for the ice volumes of the last glaciation modeled here were computed as 230 m for the Quelccaya Ice Cap and 170 m for the Cordillera Vilcanota. Maximum ELA depressions are lower than previously published: <500 m for the Cordillera Vilcanota and <400 m for the Quelccaya Ice Cap. These lower values could imply a topographic control over paleoglacier extent.     

A three-dimensional parametric study of the use of soil nails for stabilising tunnel faces

Applying an effective nailing system at a tunnel heading, not only improves the stability of the tunnel heading and limits deformation at the tunnel face, but it also reduces volume loss during excavation and hence reduces ground surface settlement. The effectiveness of a soil nail system is affected by many factors such as the diameter and stiffness of the nails. In this paper, a systematic parametric study was conducted to study the axial rigidity of a nail, E_nA_n, for improving the stability of tunnel headings and reducing ground movements in stiff clay. The parametric study involved a series of three-dimensional elasto-plastic coupled-consolidation finite element analyses. The stability of the tunnel face is improved with increasing E_nA_n. For a given nail density applied at the tunnel face, an optimum axial rigidity of the nail (E_nA_n)_opt can be identified. The efficiency of the nailing system diminishes when (E_nA_n)_opt is reached. The use of a soil nailing system reduces the magnitude of stress relief at the tunnel heading during excavation. Thus, this reduction of stress relief minimises the amount of soil yielding and excess pore water pressure generated in the soil around the tunnel heading.     

Crack damage stress as a composite function of porosity and elastic matrix stiffness in dolomites and limestones

Twenty-five uniaxial compression tests were performed to determine stress at onset of dilation, referred to herein as “the crack damage stress,” in heterogeneous dolomites and limestones. A simplified model for crack damage stress (σ_cd) is developed here using porosity, elastic modulus, Poisson's ratio and three empirical coefficients. The model shows that when porosity decreases and elastic modulus increases, σ_cd rapidly increases and approaches its maximum value. On the other hand, when porosity increases and elastic modulus decreases, σ_cd rapidly decreases and approaches its minimum value. The proposed model is validated for six heterogeneous limestone and dolomite formations which are widely distributed in Israel.     

Hydraulic Characteristics of Rough Fractures in Linear Flow under Normal and Shear Load

Summary A hydro-mechanical testing system, which is capable of measuring both the flow rates and the normal and shear displacement of a rock fracture, was built to investigate the hydraulic behaviour of rough tension fractures. Laboratory hydraulic tests in linear flow were conducted on rough rock fractures, artificially created using a splitter under various normal and shear loading. Prior to the tests, aperture distributions were determined by measuring the topography of upper and lower fracture surfaces using a laser profilometer. Experimental variograms of the initial aperture distributions were classified into four groups of geostatistical model, though the overall experimental variograms could be well fitted to the exponential model. The permeability of the rough rock fractures decayed exponentially with respect to the normal stress increase up to 5 MPa. Hydraulic behaviours during monotonic shear loading were significantly affected by the dilation occurring until the shear stress reached the peak strength. With the further dilation, the permeability of the rough fracture specimens increased more. However, beyond shear displacement of about 7 to 8 mm, permeability gradually reached a maximum threshold value. The combined effects of both asperity degradation and gouge production, which prohibited the subsequent enlargement of mean fracture aperture, mainly caused this phenomenon. Permeability changes during cyclic shear loading showed somewhat irregular variations, especially after the first shear loading cycle, due to the complex interaction from asperity degradations and production of gouge materials. The relation between hydraulic and mechanical apertures was analyzed to investigate the valid range of mechanical apertures to be applied to the cubic law. Received June 12, 2001; accepted February 26, 2002 Published online September 2, 2002