Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

Chemical alteration of tephra in the depositional environment: theoretical stability modelling

The study of the chemical stability of vitreous material in aqueous media is well‐established. There has to date been little consideration of the implications of variations in the chemical durability of tephra in Quaternary tephrochronology. Chemical alteration can take the form of cationic leaching from the matrix, or complete destruction of the silica network, either of which could constrain the ability to chemically identify distal tephra. Here we apply established models of vitreous durability to the published chemical analyses of a large number of Icelandic tephras in order to predict their relative durabilities under equivalent conditions. This suggests that some important tephras have relatively poor chemical stability, and that rhyolitic tephras are, in general, more stable than basaltic. We conclude that tephras should be expected to show predictable differential chemical stability in the post‐depositional environment. Copyright © 2003 John Wiley & Sons, Ltd.     

A Simulating Experiment in the Process of Soil Erosion on Bare Land in Mt. Tanakami

In order to understand the process of surface erosion and acquire basic data of conditions on hillslope without vege tation, a sprinkling experiment is conducted on a bare slope in Mt. Tanakami in the central part of Japan. Based on the mea surements of runoff, mean soil erosion depth, and sediment yield, etc. , the results suggest the following characteristics in the process of surface erosion in the experimental area. (1) The occurrence of sediment discharge is interrupted; (2) Surface runoff is a saturated overland flow; (3) The mean soil erosion depth is thick compared with other areas in Mt. Tanakami;(4) Sediment discharge process is detachment- limited.     

运输系统的区域效应研究

运输系统改进—运距缩短—产业布局调整这一空间过程对不同的区域可能具有完全不同的意义。本文通过对一个两区域—两部门经济模型的研究得出结论:当较落后的资源产区与经济中心之间的距离由于运输系统的改进而缩短时,资源产区的经济结构向单一化转变,初级产品的生产和输出在经济建设中逐渐占主导地位,加工工业逐渐萎缩,地方经济发展受到阻碍,而且这种变化是在距离的某一特定值处突然出现的,带有跃变的性质。这一研究结果表明,对区域运输系统的任何改进措施都应该慎而行之,以免造成与本来目标相悖的区域效应。     

THE CONSTRUCTION AND ITS DEVELOPMENT OFTHE OVERSEAS TRANSPORT SYSTEM IN NORTHEAST CHINA

ＴＨＥＣＯＮＳＴＲＵＣＴＩＯＮＡＮＤＩＴＳＤＥＶＥＬＯＰＭＥＮＴＯＦＴＨＥＯＶＥＲＳＥＡＳＴＲＡＮＳＰＯＲＴＳＹＳＴＥＭＩＮＮＯＲＴＨＥＡＳＴＣＨＩＮＡＧａｏＳｈａｌｉ（高莎丽）（ＤｅｐａｒｔｍｅｎｔｏｆＧｅｏｇｒａｐｈｙ，ＮｏｒｔｈｅａｓｔＮｏｒｍ...     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. Copyright © 2003 John Wiley & Sons, Ltd.     

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.     

TEM evidence for intracellular accumulation of lead by bacteria in subsurface environments

We monitored near-surface atmospheric fallout (15-cm above ground) and soil solution (at 15, 35 and 55 cm below ground) derived nanoparticles over an 8-month period by collecting the particles directly onto TEM grids in anthropogenically-influenced (vineyard) and pristine (native forest) sites in France. Particle clusters trapped on the grid were selected randomly and individual particles were binned into eight different groups (euhedral clays, weathered clays, salts, oxi-hydroxides, bacteria, non-living organic matter, aggregates and undetermined). Bacteria represent 9–23% of the collected nanoparticle area (ave. 9.4% and 18% for two atmospheric collection sites and ave. 23% for soil infiltration samples). Bacteria were most often associated with non-living organic matter and comprised a variety of morpho-types. Interestingly, 45% of all the bacteria analyzed by transmission electron microscopy and electron dispersive spectroscopy (TEM-EDX) showed the presence of intracellular grains significantly enriched in lead and phosphorus. Intracellular sequestration of Pb into polyphosphate bodies has been observed in the laboratory, but this is the first observation of this phenomenon in a natural environment. Furthermore, this suggests that microbial-bound Pb may be an important transport mechanism in subsurface environments.     

Landslide susceptibility mapping for a landslide-prone area (Findikli,NE of Turkey) by likelihood-frequency ratio and weighted linear combination models

Landslides are very common natural problems in the Black Sea Region of Turkey due to the steep topography, improper use of land cover and adverse climatic conditions for landslides. In the western part of region, many studies have been carried out especially in the last decade for landslide susceptibility mapping using different evaluation methods such as deterministic approach, landslide distribution, qualitative, statistical and distribution-free analyses. The purpose of this study is to produce landslide susceptibility maps of a landslide-prone area (Findikli district, Rize) located at the eastern part of the Black Sea Region of Turkey by likelihood frequency ratio (LRM) model and weighted linear combination (WLC) model and to compare the results obtained. For this purpose, landslide inventory map of the area were prepared for the years of 1983 and 1995 by detailed field surveys and aerial-photography studies. Slope angle, slope aspect, lithology, distance from drainage lines, distance from roads and the land-cover of the study area are considered as the landslide-conditioning parameters. The differences between the susceptibility maps derived by the LRM and the WLC models are relatively minor when broad-based classifications are taken into account. However, the WLC map showed more details but the other map produced by LRM model produced weak results. The reason for this result is considered to be the fact that the majority of pixels in the LRM map have high values than the WLC-derived susceptibility map. In order to validate the two susceptibility maps, both of them were compared with the landslide inventory map. Although the landslides do not exist in the very high susceptibility class of the both maps, 79% of the landslides fall into the high and very high susceptibility zones of the WLC map while this is 49% for the LRM map. This shows that the WLC model exhibited higher performance than the LRM model.