A multi‐level parallelized substructuring‐frontal solution for coupled thermo/hydro/mechanical problems in unsaturated soil

This paper proposes a multi‐level parallelized substructuring–frontal combined algorithm for the analysis of the problem of thermo/hydraulic/mechanical behaviour of unsaturated soil. Temperature, displacement, pore water pressure and pore air pressure are treated as the primary variables in a non‐linear analysis. Details are given firstly of the substructuring–frontal combined approach. The incorporation of the algorithm in a multi‐level parallel strategy is then discussed. The parallel processing can thus be carried out at different substructural levels. The method thus developed impacts, in a positive way, on both computer storage requirement and execution time. Copyright © 2003 John Wiley & Sons, Ltd.     

Himalayan magmatism and porphyry copper–molybdenum mineralization in the Yulong ore belt, East Tibet

Summary ?The NW–SE-trending Yulong porphyry Cu–Mo ore belt, situated in the Sanjiang0 area of eastern Tibet, is approximately 400 km long and 35 to 70 km wide. Complex tectonic and magmatic processes during the Himalayan epoch have given rise to favorable conditions for porphyry-type Cu–Mo mineralization. Porphyry masses of the Himalayan epoch in the Yulong ore belt are distributed in groups along regional NW–SE striking tectonic lineaments. They were emplaced mainly into Triassic and Lower Permian sedimentary-volcanic rocks. K–Ar und U–Pb isotopic datings give an intrusion age range of 57–26 Ma. The porphyries are mainly of biotite monzogranitic and biotite syenogranitic compositions. Geological and geochemical data indicate that the various porphyritic intrusions in the belt had a common or similar magma source, are metaluminous to peraluminous, Nb–Y–Ba-depleted, I-type granitoids, and belong to the high-K calc-alkaline series. Within the Yulong subvolcanic belt a number of porphyry stocks bear typical porphyry type Cu–Mo alteration and mineralization. The most prominent porphyry Co–Mo deposits include Yulong, Malasongduo, Duoxiasongduo, Mangzong and Zhanaga, of which Yulong is one of the largest porphyry Cu (Mo) deposits in China with approximately 8 × 10⁶ tons of contained Cu metal. Hydrothermal alteration at Yulong developed around a biotite–monzogranitic porphyry stock that was emplaced within Upper Triassic limestone, siltstone and mudstone. The earliest alteration was due to the effects of contact metamorphism of the country rocks and alkali metasomatism (potassic alteration) within and around the porphyry body. The alteration of this stage was accompanied by a small amount of disseminated and veinlet Cu–Mo sulfide mineralization. Later alteration–mineralization zones form more or less concentric shells around the potassic zone, around which are distributed a phyllic or quartz–sericite–pyrite zone, a silicification and argillic zone, and a propylitic zone. Fluid inclusion data indicate that three types of fluids were involved in the alteration–mineralization processes: (1) early high temperature (660–420 °C) and high salinity (30–51 wt% NaCl equiv) fluids responsible for the potassic alteration and the earliest disseminated and/or veinlet Cu–Mo sulfide mineralization; (2) intermediate unmixed fluids corresponding to phyllic alteration and most Cu–Mo sulfide mineralization, with salinities of 30–50 wt% NaCl equiv and homogenization temperatures of 460–280 °C; and (3) late low to moderate temperature (300–160 °C) and low salinity (6–13 wt% NaCl equiv) fluids responsible for argillic and propylitic alteration. Hydrogen and oxygen isotopic studies show that the early hydrothermal fluids are of magmatic origin and were succeeded by increasing amounts of meteoric-derived convective waters. Sulfur isotopes also indicate a magmatic source for the sulfur in the early sulfide mineralization, with the increasing addition of sedimentary sulfur outward from the porphyry stock. Received August 29, 2001; revised version accepted May 1, 2002 Published online: November 29, 2002     

Relationship between footwall composition,crustal contamination,and fluid–rock interaction in the Platreef,Bushveld Complex,South Africa

In the northern limb of the 2.06-Ga Bushveld Complex, the Platreef is a platinum group elements (PGE)-, Cu-, and Ni-mineralized zone of pyroxenite that developed at the intrusion margin. From north to south, the footwall rocks of the Platreef change from Archaean granite to dolomite, hornfels, and quartzite. Where the footwall is granite, the Sr-isotope system is more strongly perturbed than where the footwall is Sr-poor dolomite, in which samples show an approximate isochron relationship. The Nd-isotope system for samples of pyroxenite and hanging wall norite shows an approximate isochron relationship with an implied age of 2.17 ± 0.2 Ga and initial Nd-isotope ratio of 0.5095. Assuming an age of 2.06 Ga, the ɛNd values range from −6.2 to −9.6 (ave. −7.8, n = 17) and on average are slightly more negative than the Main Zone of the Bushveld. These data are consistent with local contamination of an already contaminated magma of Main Zone composition. The similarity in isotope composition between the Platreef pyroxenites and the hanging wall norites suggests a common origin. Where the country rock is dolomite, the Platreef has generally higher plagioclase and pyroxene δ ¹⁸O values, and this indicates assimilation of the immediate footwall. Throughout the Platreef, there is considerable petrographic evidence for sub-solidus interaction with fluids, and the Δ _{plagioclase–pyroxene} values range from −2 to +6, which indicates interaction at both high and low temperatures. Whole-rock and mineral δD values suggest that the Platreef interacted with both magmatic and meteoric water, and the lack of disturbance to the Sr-isotope system suggests that fluid–rock interaction took place soon after emplacement. Where the footwall is granite, less negative δD values suggest a greater involvement of meteoric water. Consistently higher values of Δ _{plagioclase–pyroxene} in the Platreef pyroxenites and hanging wall norites in contact with dolomite suggest prolonged interaction with CO₂-rich fluid derived from decarbonation of the footwall rocks. The overprint of post crystallization fluid–rock interaction is the probable cause of the previously documented lack of correlation between PGE and sulfide content on the small scale. The Platreef in contact with dolomite is the focus of the highest PGE grades, and this suggests that dolomite contamination played a role in PGE concentration and deposition, but the exact link remains obscure. It is a possibility that the CO₂ produced by decarbonation of assimilated dolomite enhanced the process of PGE scavenging by sulfide precipitation.     

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

In situ strength measurements on natural upper-mantle minerals

Using in situ strength measurements at pressures up to 10 GPa and at room temperature, 400, 600, and 700°C, we examined rheological properties of olivine, orthopyroxene, and chromian-spinel contained in a mantle-derived xenolith. Mineral strengths were estimated using widths of X-ray diffraction peaks as a function of pressure, temperature, and time. Differential stresses of all minerals increase with increasing pressure, but they decrease with increasing temperature because of elastic strain on compression and stress relaxation during heating. During compression at room temperature, all minerals deform plastically at differential stress of 4–6 GPa. During subsequent heating, thermally induced yielding is observed in olivine at 600°C. Neither orthopyroxene nor spinel shows complete stress relaxation, but both retain some stress even at 700°C. The strength of the minerals decreases in the order of chromian-spinel ≈ orthopyroxene > olivine for these conditions. This order of strength is consistent with the residual pressure of fluid inclusions in mantle xenoliths.     

南海沉积岩芯的地球化学研究

本文对南海中部陆坡区和海盆区的5个岩芯中的 Si,Al,Fe,Mn,Ca,Mg 和 CaCO_3成分做了分析,并对其分布特征与有关影响因素进行了讨论。受沉积区环境差异的影响,元素的分布具有区域性特征;受海区特定地理位置的影响,元素的分布呈现出过渡性特征(浅海到大洋的过渡);受沉积过程中环境变化的影响,元素的分布展示了周期性变化或旋回性。根据岩芯中元素的分布特征,并参考有关古地磁,氧碳同位素和钙质超微化石的测试分析结果,将本区岩芯划分为14层地球化学层,并将本区划分为两个沉积地球化学区,揭示了早更新世以来本区沉积环境的变化和存在的两种沉积模式。     

Constraining the origin and evolution of confined turbidite systems: southern Cretan margin,Eastern Mediterranean Sea (34°30–36°N)

Bathymetric, 9.5-kHz long-range sidescan sonar (OKEAN), seismic reflection and sediment-core data are used in the analysis of two tectonic troughs south of Crete, Eastern Mediterranean Sea. Here, up to 1.2 s two-way travel time (TWTT) of strata have accumulated since the Middle Miocene in association with extension in the South Aegean region. The study area comprises >100-km- long by >25-km-wide basins filled by sediments subdivided into two seismic units: (1) an upper Unit 1 deposited in sub-basins which follow the present-day configuration of the southern Cretan margin; (2) a basal Unit 2, more than 500 ms (TWTT) thick, accumulated in deeper half-graben/grabens distinct from the present-day depocentres. Both units overlap a locally stratified Unit 3 comprising the pre-Neogene core complex of Crete and Gavdos. In this work, the interpreted seismic units are correlated with the onshore stratigraphy, demonstrating that denudation processes occurring on Crete and Gavdos in response to major tectonic events have been responsible for high sedimentation rates along the proximal southern Cretan margin. Consequently, topographically confined sedimentary units have been deposited south of Crete in the last 12 Ma, including turbidites and other mass-flow deposits fed by evolving transverse and axial channel systems. Surface processes controlling facies distribution include the direct inflow of sediment from alluvial-fan systems and incising mountain rivers onto the Cretan slope, where significant sediment instability processes occur at present. In this setting, seismic profiles reveal eight different types of stratigraphic contacts on basin-margin highs, and basinal areas show evidence of halokinesis and/or fluid escape. The acquired data also show that significant changes to the margin’s configuration occurred in association with the post-Alpine tectonic and eustatic episodes affecting the Eastern Mediterranean.     

厦门筼筜港的形成与演变

本文根据实际调查的资料,概述了筼筜港区的地质地貌特征,重点论述了第四纪沉积物的成因和时代,阐述了该区的形成与演变过程,把它划分为5个发育阶段。     

热带地区流场异常与江淮特大暴雨成因的研究

本文通过计算和对比分析，得到１９９１年春、夏江淮地区特大暴雨形成的五个基本特征：西北太平洋上副高显著偏西、偏北；较弱的西南季风与东南气流汇合成强西南气流伸向江淮地区；偏东风与强西南风形成的辐合带在江淮地区维持；江淮地区对流层中下部为强上升运动、上部为辐散；１３０°Ｅ附近南半球向北半球的较强越赤道气流持续。     

中国海岸风成沙ESR测年的研究

介绍ＥＳＲ测年的基本原理，并讨论风成沙Ｇｅ心测年的可行性、风成沙Ｅ＇心的特征、测年误差及可靠性等问题。１９９２－１９９３年在渤、黄、东、南海沿岸采集的１００余个风成沙样品ＥＳＲ测年实验表明，绝大部分海岸风成沙的年龄为１－７万年。