Natural and man-induced stress evolution of slopes: the Monte Mario hill in Rome

The paper deals with stress-release effects induced by man-made cuts or excavations into natural stiff clay slopes that experienced erosion in response to valley deepening. The study was focused on the Monte Mario hill in Rome (Italy), which formed part of an area of recent urban expansion. The methodology of the study relied on a reference engineering-geology model, which was developed on the basis of site and laboratory data and stress–strain analyses. The latter analyses were carried out with the finite-difference code FLAC 4.0. Numerical modelling was based on a sequential approach, taking into account the main evolutionary stages of the Tiber river valley in Romeȁ9s urban area and then making cuts at the bottom of the slope located south of the Monte Mario Astronomical Observatory. The simulation revealed the stress-release effects that fluvial erosion and excavation fronts have caused on the investigated slopes and their consequent gravitational instabilities. These processes appear with metre-scale displacements, followed by stress-release cracks (actually observed on the slopes under review). In quantifying stress-release deformations, the simulation took into account the possible role of creep in the observed retardation of stress-release effects.Research activities were carried on with the co-operation of Dȁ9Arcangelo A. (Consorzio TREESSE, andrea.darcangelo@libero.it) and Moretti S. (IMG S.r.l. Servizi Tecnici per lȁ9Ingegneria e lȁ9Ambiente, moretti@img-srl.it)     

Occurrence of salt water above fresh water in dynamic equilibrium in a coastal groundwater flow system near De Panne, Belgium

A salt water lens is found above fresh water under the shore between Dunkerque (France) and Nieuwpoort (Belgium). This inverse density distribution is in a dynamic equilibrium. It develops due to the infiltration of salt water on the back shore during high tide. Under this salt water lens, water infiltrated in the adjacent dune area flows towards the sea and discharges at the seabed. This water quality distribution differs from the classic salt water wedge under fresh water described in the literature. Here, the evolution to this water quality distribution is simulated with a density dependent numerical model. A large tidal range, shore morphology and a permeable groundwater reservoir are the main conditions for the observed water quality distribution. By altering these conditions, intermediate water quality distributions between the classic salt water wedge and the one discussed here develop. Based on these simulations, it is expected that similar kinds of inverse density distribution could be present in a number of coastal areas, which have tides, a gently sloping shore and a permeable substratum.     

3D numerical modelling of fluid flow in the Val-d’Or orogenic gold district: major crustal shear zones drain fluids from overpressured vein fields

Fluid flow patterns have been determined using oxygen isotope isopleths in the Val-d’Or orogenic gold district. 3D numerical modelling of fluid flow and oxygen isotope exchange in the vein field shows that the fluid flow patterns can be reproduced if the lower boundary of the model is permeable, which represents middle or lower crustal rocks that are infiltrated by a metamorphic fluid generated at deeper levels. This boundary condition implies that the major crustal faults so conspicuous in vein fields do not act as the only major channel for upward fluid flow. The upper model boundary is impermeable except along the trace of major crustal faults where fluids are allowed to drain out of the vein field. This upper impermeable boundary condition represents a low-permeability layer in the crust that separates the overpressured fluid from the overlying hydrostatic fluid pressure regime. We propose that the role of major crustal faults in overpressured vein fields, independent of tectonic setting, is to drain hydrothermal fluids out of the vein field along a breach across an impermeable layer higher in the crust and above the vein field. This breach is crucial to allow flow out of the vein field and accumulation of metals in the fractures, and this breach has major implications for exploration for mineral resources. We propose that tectonic events that cause episodic metamorphic dehydration create a short-lived pulse of metamorphic fluid to rise along zones of transient permeability. This results in a fluid wave that propagates upward carrying metals to the mineralized area. Earthquakes along crustal shear zones cause dilation near jogs that draw fluids and deposit metals in an interconnected network of subsidiary shear zones. Fluid flow is arrested by an impermeable barrier separating the hydrostatic and lithostatic fluid pressure regimes. Fluids flow through the evolving and interconnected network of shear zones and by advection through the rock matrix. Episodic breaches in the impermeable barrier along the crustal shear zones allow fluid flow out of the vein field.     

用逐步回归预测棉铃虫发生期和发生量

收集了1973年至2000年的气象，农作物、棉铃虫虫害等196年因子，以运城、汾阳和临汾为山西省代表站，利用最优二分割法把发生量分为10级，5级和不分级3种情况；使用逐步回归计算了219个模型，从中选出27个棉铃虫二、三代发生期，发生量最优模型，在农气服务中应用，预报准确率为95％，效果良好。     

Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

Strike-slip fault bridge fluid pumping mechanism: insights from field-based palaeostress analysis and numerical modelling

We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ₁ towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ₃ rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures.     

Hybrid vibration experiments with a bridge foundation system model

In order to improve seismic design technology of bridges, it is necessary to evaluate the vibration characteristics of a bridge–soil system that consists of soil, foundation structure, pier and superstructure. However, there have been few experimental studies on seismic behavior of bridge–soil system. In this paper, we conducted the hybrid vibration experiment on seismic behavior of bridge–soil system, and examined the applicability of hybrid vibration experiment to study seismic response of bridge–soil system. Based on the experiment results, seismic response of bridge was quantitatively studied.     

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

湖南省2001年夏秋季对流云降水潜力数值模拟结果

应用对流云数值模式模拟了2001年7～9月湖南省三个探空站0000和1200(世界时)462个对流云降水算例.模拟结果表明(1)其中53天的198个算例属于有利人工增雨的天气形势,133个算例的可播度大于零,54个算例的增雨率大于零;(2)7～9月天气形势有利于人工增雨日的算例,对流云平均含水量为575万吨,平均降水效率11.4%;(3)8月降水潜力大于7月,9月降水潜力最小.这表明即使大旱的2001年湖南省夏秋季对流云仍然有一定的增雨潜力.