Recharge source and hydrogeochemical evolution of shallow groundwater in a complex alluvial fan system,southwest of North China Plain

Many cities around the world are developed at alluvial fans. With economic and industrial development and increase in population, quality and quantity of groundwater are often damaged by over-exploitation in these areas. In order to realistically assess these groundwater resources and their sustainability, it is vital to understand the recharge sources and hydrogeochemical evolution of groundwater in alluvial fans. In March 2006, groundwater and surface water were sampled for major element analysis and stable isotope (oxygen-18 and deuterium) compositions in Xinxiang, which is located at a complex alluvial fan system composed of a mountainous area, Taihang Mt. alluvial fan and Yellow River alluvial fan. In the Taihang mountainous area, the groundwater was recharged by precipitation and was characterized by Ca–HCO₃ type water with depleted δ¹⁸O and δD (mean value of −8.8‰ δ¹⁸O). Along the flow path from the mountainous area to Taihang Mt. alluvial fan, the groundwater became geochemically complex (Ca–Na–Mg–HCO₃–Cl–SO₄ type), and heavier δ¹⁸O and δD were observed (around −8‰ δ¹⁸O). Before the surface water with mean δ¹⁸O of −8.7‰ recharged to groundwater, it underwent isotopic enrichment in Taihang Mt. alluvial fan. Chemical mixture and ion exchange are expected to be responsible for the chemical evolution of groundwater in Yellow River alluvial fan. Transferred water from the Yellow River is the main source of the groundwater in the Yellow River alluvial fan in the south of the study area, and stable isotopic compositions of the groundwater (mean value of −8.8‰ δ¹⁸O) were similar to those of transferred water (−8.9‰), increasing from the southern boundary of the study area to the distal end of the fan. The groundwater underwent chemical evolution from Ca–HCO₃, Na–HCO₃, to Na–SO₄. A conceptual model, integrating stiff diagrams, is used to describe the spatial variation of recharge sources, chemical evolution, and groundwater flow paths in the complex alluvial fan aquifer system.     

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Biaxial test simulations using a packing of polygonal particles

The mechanical response of cohesionless granular materials under monotonic loading is studied by performing molecular dynamic simulations. The diversity of shapes of soil grains is modelled by using randomly generated convex polygons as granular particles. Results of the biaxial test obtained for dense and loose media show that samples achieve the same void ratio at large strains independent of their initial density state. This limit state resembles the so‐called critical state of soil mechanics, except for some stress fluctuations, which remain for large deformations. These fluctuations are studied at the micro‐mechanical level, by following the evolution of the co‐ordination number, force chains and the fraction of the sliding contacts of the sample. Copyright © 2007 John Wiley & Sons, Ltd.     

华南前寒武系基底变质杂岩和中-新生代地幔岩的稀有气体同位素组成及其大地构造意义

对华南前寒武系变质岩浆杂岩稀有气体He、Ne、Ar和Xe的系统研究表明:扬子克拉通基底为含高³He的下地壳"原始岩石层",(³He/⁴He)×10^-6比值为2.8～4.6;而华夏板块基底变质岩浆杂岩则是在缺乏³He、低(³He/⁴He)×10^-8比值(3.15～17.7)的构造环境下形成的大陆中-上地壳变质岩浆杂岩层,反映出两者基底性质迥然不同。华南中-新生代爆破岩筒He同位素组成相反,相对稳定的扬子克拉通(³He/⁴He)×10^-8比值仅0.18～4.22,而郯庐-四会-吴川断裂以东,中-新生代活动地块(太平洋构造域)(³He/⁴He)×10^-8比值高达3.7～20.5。He同位素表明郯庐-四会-吴川断裂带为切割深至地幔的边界深大断裂,是扬子克拉通与华夏板块间的边界且控制了燕山期火山-侵入岩浆向西扩展。Ar同位素组成表明华南大陆中-新生代地幔形成接近"均一"的地幔组份。¹³⁶Xe/¹³⁰Xe-¹²⁹Xe/¹³⁰Xe相关组份表明它们具有地幔柱岩石同位素组成特征。     

Flexural isostatic response of the Alps to increased Quaternary erosion recorded by foreland basin remnants,SE France

We test the hypothesis that flexural isostatic compensation of the mass removed by enhanced Quaternary erosion is responsible for uplift of the Western European Alps and their forelands. We use two well‐preserved and well‐dated (1.8 Ma) abandonment surfaces of foreland basin remnants in SE France (the Chambaran and Valensole plateaux) as passive benchmarks for tilting of the foreland. Estimating their initial slope from morphometric scaling relationships, we determine bulk post‐depositional tilting of 0.5–0.8% for these surfaces. The calculated isostatic response of the Alpine lithosphere to erosional unloading, using the method recently proposed by Champagnac et al. [Geology 35 (2007) 195–198] , yields a predicted tilting of 0.3–0.4% in the considered areas, explaining approximately half of the determined post‐depositional tilting. Such long‐term deformation being insensitive to cyclic loading/unloading because of glaciations, we suspect the other half to be related to as yet undetermined long‐wavelength and long‐lived tectonic process(es).     

A three-dimentional model of the northeast East China Sea, with application to the M2 current field

A simple three-dimensional tidal model is used to examine the M2 tidal current distribution in a northeastern part of the East China Sea, especially the vertical variation of the current in the region. Computed M2 current is compared with observations available and found to be in good agreement.Main features of the calculating method in this study are: (1) Vertical variation of the tidal current is taken as a funetion of the depth-mean velocity: (2) the method is applicable to a variety of the vertical eddy viscosities; (3) it has a fine vertical resolution, especially near the sea bootom. So, this method not only enables us to get a steady state solution easily but also depicts effects of the friction on the vertical variation of the current much better.     

Relationship between biogeochemical features of biogenic elements and flocculation in the Changjiang Estuary

ＲｅｌａｔｉｏｎｓｈｉｐｂｅｔｗｅｅｎｂｉｏｇｅｏｃｈｅｍｉｃａｌｆｅａｔｕｒｅｓｏｆｂｉｏｇｅｎｉｃｅｌｅｍｅｎｔｓａｎｄｆｌｏｃｃｕｌａｔｉｏｎｉｎｔｈｅＣｈａｎｇｊｉａｎｇＥｓｔｕａｒｙ￥ＬｉｎＹｉ＇ａｎ；ＴａｎｇＲｅｎｙｏｕ；ＬｉＹａｎ；Ｄｏ...     

East-west oscillation of sea temperature in the Equatorial Pacific and its modelling

In this paper, by analyzing the data of sea temperature in the Equatorial Western Pacific (EWP), we found that the sea temperature at sensitive district (WSST) (about at 6° N, depth of 125-150 m) is correlated negatively to the sea temperature in the Equatorial Eastern Pacific (EEP). On the basis of basic characteristics of Equatorial Pacific atmosphere and ocean, we designed a simple air-sea coupled model. Using this model,an ocean circulation was simulated. The east-west oscillation of sea temperature in Equatorial Pacific is related to this ocean circulation.