Spatial measurements of stream baseflow,a relevant method for aquifer characterization and permeability evaluation. Application to a hard‐rock aquifer,the Oman ophiolite

The structure, functioning and hydrodynamic properties of aquifers can be determined from an analysis of the spatial variability of baseflow in the streams with which they are associated. Such analyses are based on simple low‐cost measurements. Through interpreting the hydrological profiles (Q = f(A)) it is possible to locate the aquifer(s) linked to the stream network and to determine the type of interrelated flow, i.e. whether the stream drains or feeds the aquifer. Using an analytical solution developed for situations with a positive linear relationship, i.e. where the baseflow increases linearly with increasing catchment size, it is also possible to estimate the permeability of the aquifer(s) concerned at catchment scale. Applied to the hard‐rock aquifers of the Oman ophiolite, this method shows that the ‘gabbro’ aquifer is more permeable than the ‘peridotite’ aquifer. As a consequence the streams drain the peridotites and ‘leak’ into the gabbro. The hydrological profiles within the peridotite are linear and positive, and indicate homogeneity in the hydrodynamic properties of these formations at the kilometre scale. The permeability of the peridotite is estimated at 5 · 10^?7 to 5 · 10^?8 m/s. Copyright © 2004 John Wiley & Sons, Ltd.     

Metamorphic phase relations in orthopyroxene-bearing granitoids: implication for high-pressure metamorphism and prograde melting in the continental crust

In this work, the factors controlling the formation and preservation of high-pressure mineral assemblages in the metamorphosed orthopyroxene-bearing metagranitoids of the Sandmata Complex, Aravalli-Delhi Mobile Belt (ADMB), northwestern India have been modelled. The rocks range in composition from farsundite through quartz mangerite to opdalite, and with varying K₂O, Ca/(Ca + Na)_rock and FeO_tot + MgO contents. A two stage metamorphic evolution has been recorded in these rocks.
An early hydration event stabilized biotite with or without epidote at the expense of magmatic orthopyroxene and plagioclase. Subsequent high-pressure granulite facies metamorphism (∼15 kbar, ∼800 °C) of these hydrated rocks produced two rock types with contrasting mineralogy and textures. In the non-migmatitic metagranitoids, spectacular garnet ± K-feldspar ± quartz corona was formed around reacting biotite, plagioclase, quartz and/or pyroxene. In contrast, biotite ± epidote melting produced migmatites, containing porphyroblastic garnet incongruent solids and leucosomes.
Applying NCKFMASHTO T–M (H₂O) and P–T pseudosection modelling techniques, it is demonstrated that the differential response of these magmatic rocks to high-pressure metamorphism is primarily controlled by the scale of initial hydration. Rocks, which were pervasively hydrated, produced garnetiferous migmatites, while for limited hydration, the same metamorphism formed sub-solidus garnet-bearing coronae. Based on the sequence of mineral assemblage evolution and the mineral compositional zoning features in the two metagranitoids, a clockwise metamorphic P–T path is constrained for the high-pressure metamorphic event. The finding has major implications in formulating geodynamic model of crustal amalgamation in the ADMB.     

Reliability assessment of granular filters in embankment dams

Empirical criteria have been used successfully to design filters of most embankment large dam projects throughout the world. However, these empirical rules are only applicable to a particular range of soils tested in laboratory and do not take into account the variability of the base material and filter particle sizes. In addition, it is widely accepted that the safety of fill dams is mainly dependent on the reliability of their filter performance. The work herein presented consists in a new general method for assessing the probability of fulfilling any empirical filter design criteria accounting for base and filter heterogeneity by means of first‐order reliability methods (FORM), so that reliability indexes and probabilities of fulfilling any particular criteria are obtained. This method will allow engineers to estimate the safety of existing filters in terms of probability of fulfilling their design criteria and might also be used as a decision tool on sampling needs and material size tolerances during construction. In addition, sensitivity analysis makes possible to analyse how reliabilities are influenced by different sources of input data. Finally, in case of a portfolio risk assessment, this method will allow engineers to compare the safety of several existing dams in order to prioritize safety investments and it is expected to be a very useful tool to evaluate probabilities of failure due to internal erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

PROPAGATION OF 30—60 DAY LOW FREQUENCY OSCILLATIONS AND THEIR INFLUENCE UPON THE SUBTROPICAL WESTERLIES JET STREAM DURING NORTHERN HEMISPHERE WINTER

Based on daily ECMWF gridpoint data of two winters during 1981—1983 including an ENSOyear,propagation of low frequency oscillations(LFO)during Northern Hemisphere winters andtheir influences upon 30—60 day oscillations of the subtropical jet stream are studied with the sta-tistical methods as complex empirical orthogonal function(CEOF)and so on.Results show that inthe winter of a normal year(1981—1982),30—60 day oscillations in the subtropical zone aremainly in the northern and southern flanks of exit region of jet stream.In the ENSO year(1982—1983),they are mainly in the vicinity of entrance and exit regions of jet stream.Intraseasonalchanges of subtropical jet stream manifested themselves as latitudinal fluctuation or longitudinalprogression or regression of about 40 day period.There are marked differences between propagat-ing passages of low frequency modes responsible for changes of subtropical jet stream in the normalyear(1981—1982)and in the ENSO year(1982—1983).Changes of oscillation amplitude showobvious phases.In general,the one in late winter is stronger than that in early winter,strongestone occurs in February.     

On salinity regimes and the vertical structure of residual flows in narrow tidal estuaries

An examination is made of the circulation in narrow estuaries subject to a predominant tidal forcing. Velocity structures are derived separately for residual flow components associated with (a) river flow, (b) wind stress, (c) a well-mixed longitudinal density gradient and (d) a fully stratified saline wedge. Dimensionless parameters are introduced to indicate the magnitude of each component and these parameters are evaluated for 9 major estuaries, thereby revealing their sensitivity to each component.For a channel of constant breadth and depth, formulae are deduced for the length of saline intrusion, L. Comparisons with observed data show that such formulae may be used with confidence to predict changes in L arising from variations in river flow, tidal range or channel depths.The level of stratification is shown to be related to a product of two parameters, one associated with velocity structure and a second involving the square of the ‘flow ratio’ $\text{u}\text{u}\text{?}$ (i.e. residual velocity/amplitude of the tidal velocity). This relationship provides a simple classification system for estuarine stratification which can be used to indicate the sensitivity of any particular estuary to changing conditions.     

全球Ms≥7￣3／4级地震的定向迁移

对全球尺度的６条大地震带内１９００～１９９０年中１８４次Ｍｓ≥７３／４级地震进行了沿地震带方向定向迁移的分析，获得了全球统一的地震定向迁移规律，总体是由西向东，迁移速度由７００ｋｍ／ａ变为１５０ｋｍ／ａ，此现象可以作多种暂态地球动力作用过程的推论，如以大西洋脊间歇式张裂引起上地幔软流物质自西向东运动，呈现纵波式的振荡传播；也可解释为非洲板块、阿拉伯板块和印度板块自西南向东北对欧亚地震带依次的推压引起向东的应变波的传播；太平洋脊两侧洋底板块向西北和东北两侧的斜向推压，可能是造成两侧地震带地震向北迁移的触发源     

Short-term Climate Characteristics at Ny-(A)lesund over the Arctic Tundra Area

Based on the Germany Koldwey Station's 1994-2003 conventional observation hourly data, this paper conducts a statistical analysis on the short-term climate characteristics for an arctic tundra region (Ny-(A)lesund island) where our first arctic expedition station (Huanghe Station) was located. Affected by the North Atlantic warming current, this area has a humid temperate climate, and the air temperature at Ny-(A)lesund rose above 0 ℃ even during deep winter season during our research period. The wind speed in this area was low and appeared most at southeast direction. We find that the temperature at Ny-(A)lesund rose in the faster rate (0.68 ℃/10 a) than those at the whole Arctic area. Compared with the floating ices where our expedition conducted in the Arctic, Ny-(A)lesund was warmer and more humid and had lower wind speed. Comparison of the near surface air temperature derived by NCEP/NCAR reanalysis to the conventional measurements conducted at the Koldwey site in Ny-(A)lesund area shows a good agreement for winter season and a significant difference for summer season.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.