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.     

Centrifuge modelling-based verification of response spectra design methods for some vulnerable structures

Residential RC framed structures suffered heavily during the 2001 Bhuj earthquake in Gujarat, India. These types of structures also saw severe damage in other earthquakes such as the 1999 Kocaeli earthquake in Turkey and 921 Ji-Ji earthquake in Taiwan. In this paper the seismic response of residential structures was investigated using physical modelling. Idealised soft storey and top heavy, two degrees of freedom (2DOF) portal frame structures were developed and tested on saturated and dry sand models at 25 g using the Schofield Centre 10-m Beam Centrifuge. It was possible to recreate observed field behaviour using these models. As observed in many of the recent earthquakes, soft storey structures were found to be particularly vulnerable to seismic loads. Elastic response spectra methods are often used in the design of simple portal frame structures. The seismic risk of these structures can be significantly increased due to modifications such as removal of a column or addition of heavy water tanks on the roof. The experimental data from the dynamic centrifuge tests on such soft storey or top-heavy models was used to evaluate the predictions obtained from the response spectra. Response spectra were able to predict seismic response during small to moderate intensity earthquakes, but became inaccurate during strong earthquakes and when soil structure interaction effects became important. Re-evaluation of seismic risk of such modified structures is required and time domain analyses suggested by building codes such as IBC, UBC or NEHRP may be more appropriate.     

Observations And Modelling Of The On-Ice And Off-Ice Air Flow Over The Northern Baltic Sea

Two cases of on-ice and off-ice air flow characterizing the opposite weather situations over the ice-edge zone in the northern Baltic Sea are analysed on the basis of aircraft observations, and modelled using atwo-dimensional mesoscale model. The stable boundary layer (SBL) during theon-ice flow exhibited little thermal modification, but a low-level jet (LLJ) was generated at the 250-m high top of the SBL. In the model, the LLJ was associated with inertial oscillations in space, while the baroclinicity explained the shape of the wind profile well above the SBL. Although the observed LLJ was most pronounced over the ice, the modelling suggests that it was not generated by the ice edge but by the coastline some 400 km upwind of the ice edge, where a much more drastic change in the thermal stratification and surface roughness took place. The generation, maintenance, and strength of the LLJ were very sensitive to the parameterization of turbulent mixing in the SBL. In the case of the off-ice flow, the modification of the air mass and the development of a convective boundary layer (CBL) both over the ice and open sea were reasonably well modelled. Sensitivity runs suggested that it was essential to take into account the effects of subgrid-scale leads, a forest in the archipelago (which was crossed by the air flow), and water vapour condensationinto ice crystals. The heat flux from leads was particularly important for the heatbudget of the CBL, and the observed growth of the CBL was partly due to theeffective mixing over the rough and relatively warm forest.     

An integrated hydrological model for rain‐induced landslide prediction

This paper describes an extension to the Combined Hydrology And Stability Model (CHASM) to fully include the effects of vegetation and slope plan topography on slope stability. The resultant physically based numerical model is designed to be applied to site‐specific slopes in which a detailed assessment of unsaturated and saturated hydrology is required in relation to vegetation, topography and slope stability. Applications are made to the Hawke's Bay region in New Zealand where shallow‐seated instability is strongly associated with spatial and temporal trends in vegetation cover types, and the Mid‐Levels region in Hong Kong, an area subject to a variety of landslide mechanisms, some of which may be subject to strong topographic control. An improved understanding of process mechanism, afforded by the model, is critical for reliable and appropriate design of slope stabilization and remedial measures. Copyright © 2002 John Wiley & Sons, Ltd.     

A mathematical hysteretic model for elastomeric isolation bearings

An analytical model for high damping elastomeric isolation bearings is presented in this paper. The model is used to describe mathematically the damping force and restoring force of the rubber material and bearing. Ten parameters to be identified from cyclic loading tests are included in the model. The sensitivity of the ten parameters in affecting the model is examined. These ten parameters are functions of a number of influence factors on the elastomer such as the rubber compound, Mullins effect, scragging effect, frequency, temperature and axial load. In this study, however, only the Mullins effect, scragging effect, frequency and temperature are investigated. Both material tests and shaking table tests were performed to validate the proposed model. Based on the comparison between the experimental and the analytical results, it is found that the proposed analytical model is capable of predicting the shear force–displacement hysteresis very accurately for both rubber material and bearing under cyclic loading reversals. The seismic response time histories of the bearing can also be captured, using the proposed analytical model, with a practically acceptable precision. Copyright © 2002 John Wiley & Sons, Ltd.     

Modelling evaporation using an artificial neural network algorithm

This paper investigates the prediction of Class A pan evaporation using the artificial neural network (ANN) technique. The ANN back propagation algorithm has been evaluated for its applicability for predicting evaporation from minimum climatic data. Four combinations of input data were considered and the resulting values of evaporation were analysed and compared with those of existing models. The results from this study suggest that the neural computing technique could be employed successfully in modelling the evaporation process from the available climatic data set. However, an analysis of the residuals from the ANN models developed revealed that the models showed significant error in predictions during the validation, implying loss of generalization properties of ANN models unless trained carefully. The study indicated that evaporation values could be reasonably estimated using temperature data only through the ANN technique. This would be of much use in instances where data availability is limited. Copyright © 2002 John Wiley & Sons, Ltd.     

High-T, low-P metamorphism in the Palaeoproterozoic Halls Creek Orogen, northern Australia: the middle crustal response to a mantle-related transient thermal pulse

High‐T, low‐P metamorphic rocks of the Palaeoproterozoic central Halls Creek Orogen in northern Australia are characterised by low radiogenic heat production, high upper crustal thermal gradients (locally exceeding 40 °C km^?1) sustained for over 30 Myr, and a large number of layered mafic‐ultramafic intrusions with mantle‐related geochemical signatures. In order to account for this combination of geological and thermal characteristics, we model the middle crustal response to a transient mantle‐related heat pulse resulting from a temporary reduction in the thickness of the mantle lithosphere. This mechanism has the potential to raise mid‐crustal temperatures by 150–400 °C within 10–20 Myr following initiation of the mantle temperature anomaly, via conductive dissipation through the crust. The magnitude and timing of maximum temperatures attained depend strongly on the proximity, duration and lateral extent of the thermal anomaly in the mantle lithosphere, and decrease sharply in response to anomalies that are seated deeper than 50–60 km, maintained for <5 Myr in duration and/or have half‐widths <100 km. Maximum temperatures are also intimately linked to the thermal properties of the model crust, primarily due to their influence on the steady‐state (background) thermal gradient. The amplitudes of temperature increases in the crust are principally a function of depth, and are broadly independent of crustal thermal parameters. Mid‐crustal felsic and mafic plutonism is a predictable consequence of perturbed thermal regimes in the mantle and the lowermost crust, and the advection of voluminous magmas has the potential to raise temperatures in the middle crust very quickly. Although pluton‐related thermal signatures significantly dissipate within <10 Myr (even for very large, high‐temperature intrusive bodies), the interaction of pluton‐ and mantle‐related thermal effects has the potential to maintain host rock temperatures in excess of 400–450 °C for up to 30 Myr in some parts of the mid‐crust. The numerical models presented here support the notion that transient mantle‐related heat sources have the capacity to contribute significantly to the thermal budget of metamorphism in high‐T, low‐P metamorphic belts, especially in those characterised by low surface heat flow, very high peak metamorphic geothermal gradients and abundant mafic intrusions.     

The geological and geodynamic evolution of the eastern Black Sea basin: insights from 2-D and 3-D tectonic modelling

Subsidence mechanisms that may have controlled the evolution of the eastern Black Sea have been studied and simulated using a numerical model that integrates structural, thermal, isostatic and surface processes in both two- (2-D) and three-dimensions (3-D). The model enables the forward modelling of extensional basin evolution followed by deformation due to subsequent extensional and compressional events. Seismic data show that the eastern Black Sea has evolved via a sequence of interrelated tectonic events that began with early Tertiary rifting followed by several phases of compression, mainly confined to the edges of the basin. A large magnitude (approximately 12 km) of regional subsidence also occurred in the central basin throughout the Tertiary. Models that simulate the magnitude of observed fault controlled extension (β=1.13) do not reproduce the total depth of the basin. Similarly, the modelling of compressional deformation around the edges of the basin does little to enhance subsidence in the central basin. A modelling approach that quantifies lithosphere extension according to the amount of observed crustal thinning and thickening across the basin provides the closest match to overall subsidence. The modelling also shows that deep crustal and mantle–lithosphere processes can significantly influence the rate and magnitude of syn- to post-rift subsidence and shows that such mechanisms may have played an important role in forming the anomalously thin syn-rift and thick Miocene–Quaternary sequences observed in the basin. It is also suggested that extension of a 40–45 km thick pre-rift crust is required to generate the observed magnitude of total subsidence when considering a realistic bathymetry.     

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.