Numerical modelling of deformation of multi-directional random waves over a varying topography

Two numerical formulations of the breaking phenomenon were implemented in a numerical model for random wave propagation based on the elliptic formulation of the mild-slope equation. The randomness of the wave field was simulated based on a spectral component method, in which the 3-D spectrum is discretised in components of equal energy. One of the breaking process formulations is based on the concept of breaking each independent spectral component. The other is based on the distribution of the local amount of energy dissipated through the independent spectral components. The model based on the concept of breaking each independent spectral component produces the best estimates of the wave field, when the numerical results are compared with laboratory data.     

Mirdita Zone ophiolites and associated sediments in Albania reveal Neotethys Ocean origin

The Mirdita Ophiolite Zone in Albania is associated with widespread mélanges containing components of up to nappe-size. We dated matrix and components of the mélange by radiolarians, conodonts, and other taxa. The components consist of radiolarites, pelagic limestones and shallow-water limestones, all of Triassic age, as well as ophiolites. Triassic radiolarite as a primary cover of ophiolite material proves Middle Triassic onset of Mirdita ocean-floor formation. The mélange contains a turbiditic radiolarite-rich matrix (“radiolaritic flysch”), dated as Late Bajocian to Early Oxfordian. It formed as a synorogenic sediment during west-directed thrusting of ophiolite and sediment-cover nappes representing ocean floor and underplated fragments of the western continental margin. The tectonic structures formed during these orogenic events (“Younger Kimmeridian or Eohellenic Orogeny”) are sealed by Late Jurassic platform carbonates. The geological history conforms with that of the Inner Dinarides and adjoining areas; we therefore correlate the Mirdita-Pindos Ophiolite Zone with the Vardar Zone and explain its present position by far-distance west-directed thrusting.     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Dynamical evolution of NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some fast-track dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.On leave from the Department of Mathematics, University of Pisa, Via Buonarroti 2, 56127 Pisa, Italy, thanks to the G. Colombo fellowships of the European Space Agency.     

Problems in relating the optical and X-ray emissions from a solar flare

We evaluate the possibility that the short-lived Balmer line emission at H9 3835 Å of the 1972, August 2 (1839 UT) solar flare is due to heating of the chromosphere by bombarding electrons. We point out some of the problems of comparing the time behavior and spatial distribution of simultaneous hard and soft X-ray emissions. It is concluded that the present data do not justify the attribution of the short-lived optical emission to the presumed hard X-ray producing electrons.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

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     

Phénomènes karstiques fossiles et actuels au sein des formations métamorphiques silico-alumineuses de la nappe pan-africaine de Yaoundé (Sud-Cameroun)Fossil and present-time karstic phenomena in silico-aluminous metamorphic formations of the Pan-African nappe of Yaoundé (South-Cameroon)

Numerous karstic features have been recognised in the non-carbonaceous micaschists and gneisses of the Yaoundé Pan-African nappe, south of Cameroon. It is shown that their formation was controlled by the structural features of the rocks. The wells and the pipes in the bedrock outcrops point out a current karstification process, resulting from the plagioclase dissolution by the acid rain waters. Hill wall alveoli and caves, of pre-Miocene age, are exhumed features that were done by dissolution in the aquiferous underground. To cite this article: J.-P. Vicat et al., C. R. Geoscience 334 (2002) 545–550.