A multi-method approach to study the stability of natural slopes and landslide susceptibility mapping

In this paper, a multi-method approach for the assessment of the stability of natural slopes and landslide hazard mapping applied to the Dakar coastal region is presented. This approach is based on the effective combination of geotechnical field and laboratory works, of GIS, and of mechanical (deterministic and numerical) stability analysis. By using this approach, valuable results were gained regarding instability factors, landslide kinematics, simulation of slope failure and coastal erosion. This led to a thorough assessment and strong reduction in the subjectivity of the slope stability and hazard assessment and to the development of an objective landslide danger map of the SW coast of Dakar. Analysis of the results shows that the slides were influenced by the geotechnical properties of the soil, the weathering, the hydrogeological situation, and the erosion by waves. The landslide susceptibility assessment based on this methodological approach has allowed for an appropriate and adequate consideration of the multiple factors affecting the stability and the optimization of planning and investment for land development in the city.     

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.     

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     

The role of diapycnal mixing for the equilibrium response of thermohaline circulation

Using a zonally averaged, one-hemispheric numerical model of the thermohaline circulation, the dependence of the overturning strength on the surface equator-to-pole density difference is investigated. It is found that the qualitative behavior of the thermohaline circulation depends crucially on the nature of the small-scale vertical mixing in the interior of the ocean. Two different representations of this process are considered: constant vertical diffusivity and the case where the rate of mixing energy supply is taken to be a fixed quantity, implying that the vertical diffusivity decreases with increasing stability of the water column. When the stability-dependent diffusivity parameterization is applied, a weaker density difference is associated with a stronger circulation, contrary to the results for a fixed diffusivity. A counterintuitive consequence of the stability-dependent mixing is that the poleward atmospheric freshwater flux, which acts to reduce the thermally imposed density contrast, strengthens the thermally dominated circulation and its attendant poleward heat transport. However, for a critical value of the freshwater forcing, the thermally dominated branch of steady states becomes unstable, and is succeeded by strongly time-dependent states that oscillate between phases of forward and partly reversed circulation. When a constant vertical diffusivity is employed, on the other hand, the thermally dominated circulation is replaced by a steady salinity-dominated state with reversed flow. Thus in this model, the features of the vertical mixing are essential for the steady-state response to freshwater forcing as well as for the character of flow that is attained when the thermally dominated circulation becomes unstable.Responsible Editor: Jin-Song von Storch     

The role of the seismic trigger in the Calitri landslide (Italy): historical reconstruction and dynamic analysis

The study covered by this paper was focused on the historical case of the Calitri landslide, which was repeatedly reactivated by earthquakes, as reported since 1694. The town of Calitri (Southern Italy) is located on a ridge whose southern slope, from its top to the Ofanto river valley floor, has been historically affected by major landsliding. The last record of recurrence of the Calitri landslide leads back to the 1980 Irpinia earthquake, which caused significant damage to the town and had pervasive and visible ground effects. Based on a detailed historical reconstruction of landsliding and seismicity at Calitri, the study analysed the current static and dynamic stability of the landsliding slope by means of a finite-difference numerical analysis, taking into account the various factors of landslide initiation and kinematics. The results of the numerical analysis: (i) were consistent with the roto-translational mechanism observed upon the latest reactivation of the landslide; (ii) demonstrated that excess pore pressure redistribution caused a lag between the seismic trigger and the initiation of landsliding; and (iii) showed the impact of seismic input frequency on propagation and depth of slope instability.     

Western Alpine back-thrusting as subsidence mechanism in the Tertiary Piedmont Basin (Western Po Plain, NW Italy)

Basin formation dynamics of the Tertiary Piedmont Basin (TPB) are here investigated by means of cross-section numerical modelling. Previous works hypothesised that basin subsidence occurred due first to extension (Oligocene) and then to subsequent loading due to back-thrusting (Miocene). However, structural evidence shows that the TPB was mainly under contraction from Oligocene until post Pliocene time while extension played a minor role. Furthermore, thermal indicators strongly call for a cold (flexure-induced) mechanism but are strictly inconsistent with a hot (thermally induced) mechanism. Our new modelling shows that the TPB stratigraphic features can be reproduced by flexure of a visco-elastic plate loaded by back-thrusts active in the Western Alps in Oligo-Miocene times. Far-field compression contributed to the TPB subsidence and controlled the basin infill geometry by enhancing basin tilting, forebulge uplift and erosion of the southern margin of the basin. These results suggest that the TPB subsidence is the result of a combination of mechanisms including thrust loading and far-field compressional stresses.     

不同对流参数化方案和初值场组合对四川两次暴雨影响的数值模拟

针对2004年9月2~6日四川盆地东北部一次持续性暴雨过程和2004年6月30日成都地区暴雨过程,利用MM5中尺度数值模式,进行了Grell、Kuo、KF及BM四种积云对流参数化方案与T213模式分析值及NCEP再分析值两种初始值组合的数值试验,分析了不同积云方案和初值组合的模拟性能。比较结果表明,MM5模式较好地再现了这两次强降水过程。在同一模式中,不同对流参数化方案和不同初值组合在降水落区和强度模拟上存在一定程度的差异。初步模拟表明NCEP资料模拟的降水强度较T213资料模拟的降水强度偏弱,但降水落区较T213资料模拟的更接近实况,T213资料模拟的降水空报现象较严重。以T213资料模拟得到的高度、温度场都较NCEP资料更有利于影响系统的加强,模拟的湿度场、散度场和垂直运动场都较NCEP资料更有利于强降水的发生。相对而言,NCEP资料模拟结果更接近实况,Kuo和Grell方案对初值表现出更为敏感。     

Phase Equilibria of the Lyngdal Granodiorite (Norway): Implications for the Origin of Metaluminous Ferroan Granitoids

The Proterozoic (950 Ma) Lyngdal granodiorite of southern Norwaybelongs to a series of hornblende–biotite metaluminousferroan granitoids (HBG suite) coeval with the post-collisionalRogaland Anorthosite–Mangerite–Charnockite (AMC)suite. This granitoid massif shares many geochemical characteristicswith rapakivi granitoids, yet granodiorites dominate over granites.To constrain both crystallization (P, T, fO₂, H₂O in melt) andmagma generation conditions, we performed crystallization experimentson two samples of the Lyngdal granodiorite (with 60 and 65 wt% SiO₂) at 4–2 kbar, mainly at fO₂ of NNO (nickel–nickeloxide) to NNO + 1, and under fluid-saturated conditions withvarious H₂O–CO₂ ratios for each temperature. Comparisonbetween experimental phase equilibria and the mineral assemblagein the Lyngdal granodiorite indicates that it crystallized between4 and 2 kbar, from a magma with 5–6 wt % H₂O at an fO₂of NNO to NNO + 1. These oxidized and wet conditions sharplycontrast with the dry and reduced conditions inferred for thepetrogenesis of the AMC suite and many other rapakivi granitesworldwide. The high liquidus temperature and H₂O content ofthe Lyngdal granodiorite imply that it is not a primary magmaproduced by the partial melting of the crust but is derivedby the fractionation of a mafic magma. Lyngdal-type magmas appearto have volcanic equivalents in the geological record. In particular,our results show that oxidized high-silica rhyolites, such asthe Bishop Tuff, could be derived via fractionation of oxidizedintermediate magmas and do not necessarily represent primarycrustal melts. This study underlines the great variability ofcrystallization conditions (from anhydrous to hydrous and reducedto oxidized) and petrogenetic processes among the metaluminousferroan magmas of intermediate compositions (granodiorites,quartz mangerites, quartz latites), suggesting that there isnot a single model to explain these rocks. KEY WORDS: ferroan granitoids; crystallization conditions; experiments; Norway; Sveconorwegian; Bishop Tuff