Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.     

Morphological analysis of the drainage system in the Eastern Alps

We study the morphology of the major rivers draining the Eastern Alps to test whether the active tectonics of this part of the orogen is reflected in the shape of channel profiles of the river network. In our approach we compare channel profiles measured from digital elevation models with numerically modelled channel profiles using a stream power approach. It is shown that regions of high stream power coincide largely with regions of highest topography and largest uplift rates, while the forelands and the Pannonian Basin are characterised by a significantly lower stream power. From stream power modelling we conclude that there is young uplift at the very east of the Eastern Alps, in the Bohemian Massif and in the Pohorje Range. The impact of the Pleistocene glaciations is explored by comparing properties of rivers that drain in proximal and distal positions relative to the ice sheet during the last glacial maximum. Our analysis shows that most knick points, wind gaps and other non-equilibrium features of catchments covered by ice during the last glaciations (Salzach, Enns) can be correlated with glacial processes. In contrast the ice free catchments of the Mur and Drava are characterized by channels in morphological equilibrium at the first approximation and are showing only weak evidence of the strong tectonic activity within these catchments. Finally, the channel profiles of the Adige and the divide between the upper Rhine and Danube catchments differ significantly from the other catchments. We relate this to the fact that the Adige and the Rhine respond to different base levels from the remainder of the Eastern Alps: The Adige may preserve a record from the Messininan base level change and the Rhine is subject to the base level lowering in the Rhine Graben.     

辽河断陷大民屯凹陷静北潜山孢粉化石新发现

静北潜山的油气勘探始于上世纪70年代初,并将钻井揭示的未见化石的碳酸盐岩和石英砂岩沉积划归中元古界长城系,沿用至今.2002年在对该区的地层研究中,首次发现了晚古生代孢粉化石,从而证实静北潜山存在晚古生代地层.本文以孢粉化石为主要依据,对潜山地层进行了重新划分,并与邻区进行对比.在此基础上对本区的构造和沉积作一讨论.     

汶川大地震（MS 8.0）同震变形作用及其与地质灾害的关系

2008年5月12日发生于四川盆地西部龙门山断裂带的汶川大地震（MS 8.0）波及半个亚洲,震撼整个中国。本文通过地震后的实地调查,对发育在龙门山断裂带上的同震地表破裂带的分布、产状、继承性复活与变形特征,以及同震变形与地震地质灾害的关系等进行了初步总结,分析表明这次汶川大地震(MS 8.0)沿北川-映秀逆冲断裂和安县-灌县逆冲断裂同时发生地表破裂,前者产生以高角度逆冲兼右旋走滑为特征的地表破裂带长约275 km,后者产生以缓倾角逆冲作用为特征的地表破裂带长约80 km。汶川大地震的同震地表破裂带分布具有分段性特征,并与地表破坏程度的分带性有着一定的内在联系,详细研究表明,同震地表破裂带的产状直接影响地表破坏程度和地震地质灾害的强度,汶川大地震（MS 8.0）沿呈高角度陡倾的北川-映秀逆冲断裂发育的同震地表变形所产生的地表破坏程度和地震地质灾害的强度比沿缓倾角的安县-灌县逆冲断裂要强。从各种类型的地震断裂来看,具有垂直运动的逆冲型地震断裂所造成的地表破坏程度和地质灾害强度比具水平运动的走滑型地震断裂要强。因此,汶川大地震发生的破裂过程和同震地表变形与地震地质灾害的关系值得深入研究。     

Development of a 3D GIS and its application to karst areas

There is a growing interest in modeling and analyzing karst phenomena in three dimensions. This paper integrates geology, groundwater hydrology, geographic information system (GIS), database management system (DBMS), visualization and data mining to study karst features in Huaibei, China. The 3D geo-objects retrieved from the karst area are analyzed and mapped into different abstract levels. The spatial relationships among the objects are constructed by a dual-linker. The shapes of the 3D objects and the topological models with attributes are stored and maintained in the DBMS. Spatial analysis was then used to integrate the data in the DBMS and the 3D model to form a virtual reality (VR) to provide analytical functions such as distribution analysis, correlation query, and probability assessment. The research successfully implements 3D modeling and analyses in the karst area, and meanwhile provides an efficient tool for government policy-makers to set out restrictions on water resource development in the area.     

The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

Groundwater flow and contaminant transport modelling at an air weapons range

Numerical modelling was done at the Cold Lake Air Weapons Range, Canada, to test whether the dissolved RDX and nitrate detected in groundwater come from the same sources, and to predict whether contamination poses a threat to the surface water receptors near the site. Military live fire training activities may indeed pose a risk of contamination to groundwater resources, however field investigations on military bases are quite recent, and little information is available on the long-term behaviour of munition residues related contaminants. Very limited information was available about the contaminant source zones, which were assigned based on our knowledge of current training activities. The RDX plume was well represented with the model, but the heterogeneous distribution of nitrate concentrations was more difficult to reproduce. It was nonetheless determined that both contaminants originate from the same areas. According to the model, both contaminants should reach the nearby river, but concentrations in the river should remain very low if the source zone concentration does not change. Finally, the model allowed the recommendation of a new location for the main bombing target, which would offer added protection to the river and the lake into which it flows.     

EPR investigation of the methyl radical,the hydrogen atom and carbon oxide radicals in Maxixe-type beryl

The EPR spectra of Maxixe-type beryl contain a large number of overlapping signals. The angular dependence of the 1:3:3:1 signal typical for the CH₃ radical shows that this radical is located at the center of the channel cavity with its symmetry axis parallel to the crystal c-axis and is rotating around this axis. Its EPR spectrum is axially symmetric with g _// = 2.00263, g _⊥ = 2.00249 and A_// = 2.288 mT, A_⊥ = 2.256 mT. These anisotropies have the opposite signs of those found for surface-adsorbed methyl radicals. Hydrogen atoms are located at position 2a at the center of the beryl cavity and the EPR parameters of the narrow doublet signal are A ₀ = 1,407 MHz and g = 2.00230. Another doublet signal, which is broader and has axial symmetry with g _// = 2.00265, g _⊥ = 2.00625 and A_// = 0.895 mT, A_⊥ = 0.885 mT, could come from a HCO₃ radical. One narrow and easily saturated signal with g _// = 2.00227 and g _⊥ = 2.00386 is interpreted to arise from a carbon monoxide radical in the beryl channel, oriented with its axis parallel to the crystal c-axis. Additional weak doublet lines, which have similar g values as the carbon monoxide radical, are created by nearby hydrogens. A powder spectrum with g _// = 2.0017 and g _⊥ = 2.0004 appears upon UV irradiation of the single crystal and is easily saturated. This spectrum is interpreted to arise from a carbon dioxide radical, which rotates around its symmetry axis.     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran