Conceptual hydrodynamic model of the Pamukkale hydrothermal field, southwestern Turkey, based on hydrochemical and isotopic data

The results of study on the hydrochemical and isotope characteristics of shallow and deep waters at Pamukkale hydrothermal field Turkey are described in order to obtain a better understanding of the hydrological circulation. The field can be grouped into two groundwater sub-systems; cold water springs of Ca–HCO₃ type (10–12 °C), and CO₂-rich thermal waters of Ca–HCO₃–SO₄ type (25–58 °C). The occurrence of these water types is closely related to the morphology of the region, where intense tectonism formed horst and graben structures. Hence, two hydrogeological systems were defined: a deep geothermal system which is related to extensive and deep circulation of meteoric water in the regional flow system, and a shallow system which is related to local groundwater flow through sedimentary strata. The meteoric water falling at higher elevations percolates to the local groundwater system at a shallow level and flows to the deep geothermal system. During a deep convection cycle from a recharge to discharge area, the cold water attains heat from the asthenospheric intrusions, causing it to ascend. Variations of chemical and isotopic composition of thermal waters result from their mixing with cool groundwater in a shallow aquifer during their ascent to the surface.

---

Résumé Les résultats d'une étude portant sur les caractéristiques hydrochimiques et isotopiques de puits phréatiques et profonds situés dans le champ hydrothermal de Pamukkale, sont décrits de telle manière à éclairer le fonctionnement des circulations hydrologiques. Le champ peut étre divisé en deux sous-systèmes d'eaux souterraines, l'un avec des eaux de sources froides (10–12 °C) de type Ca–HCO3, et les eaux thermales (25–58 °C) riches en CO₂ et de type Ca–HCO₃–SO4. L'occurrence de ces types d'eaux est fermement liée à la morphologie de la région, oùne tectonique intense a engendré des structures en horsts et en grabens. Dés lors deux systèmes hydrogéologiques ont été définis : un système profond, qui est lié à la circulation extensive et profonde des eaux météoritiques dans le système régional d'écoulement, et un système phréatique lié aux écoulements locaux des eaux souterraines à travers les strates sédimentaires. Les eaux météoritiques aux altitudes élevées, percolent jusqu'aux systèmes locaux phréatiques, puis coulent jusqu'aux systèmes géothermaux plus profonds. Durant le cycle de convection profond des zones de recharge jusqu'aux zones de décharge, l'eau froide atteint les zones chaudes liées aux intrusions athenosphériques, provoquant la remontée. Les variations de la composition chimique et isotopique des eaux thermales résultent dans leurs mélanges avec des eaux souterraines froides dans les aquifères phréatiques durant leur remontée jusqu'à la surface.

---

Resumen Se describen los resultados del estudio de las características isotópicas e hidroquímicas de las aguas someras y profundas para obtener un mejor entendimiento de la circulación hidrológica del campo hidrotermal Pamukkale. El campo puede agruparse en dos sub-sistemas de agua subterránea: manantiales de agua fría del tipo Ca–HCO₃ (10–12°C) y aguas termales ricas en CO₂ del tipo Ca–HCO₃–SO₄(25–58°C). El ambiente de estos tipos de aguas se relaciona estrechamente con la morfología de la región donde el tectonismo intenso ha formado estructuras extensionales tipo graben y horst. Se definieron dos sistemas hidrogeológicos: un sistema geotermal profundo que se relaciona con la circulación profunda y extensa de agua meteórica en el sistema regional de flujo y un sistema somero el cual se relaciona con flujo local de agua subterránea a través de estratos sedimentarios. El agua meteórica que cae en altas elevaciones percola al sistema local de agua subterránea en un nivel somero y fluye hacia el sistema geotermal profundo. Durante un ciclo de convección del área de recarga hacia la zona de descarga, el agua fría se calienta a partir de los intrusivos astenosféricos lo que ocasiona que asciendan. Como resultado de la mezcla de las aguas recalentadas, con agua subterránea fría en un acuífero somero durante el ascenso hacia la superficie, se derivan variaciones en la composición química e isotópica de las aguas termales.

     

Driving forces of tropical deforestation: The role of remote sensing and spatial models

Remote sensing technologies are increasingly used to monitor landscape change in many parts of the world. While the availability of extensive and timely imagery from various satellite sensors can aid in identifying the rates and patterns of deforestation, modelling techniques can evaluate the socioeconomic and biophysical forces driving deforestation processes. This paper briefly reviews some emerging spatial methodologies aimed at identifying driving forces of land use change and applies one such methodology to understand deforestation in Mexico. Satellite image classification, change analysis and econometric modelling are used to identify the rates, hotspots and drivers of deforestation in a case study of the southern Yucatán peninsular region, an enumerated global hotspot of biodiversity and tropical deforestation. In particular, the relative roles of biophysical and socioeconomic factors in driving regional deforestation rates are evaluated. Such methodological approaches can be applied to other regions of the forested tropics and contribute insights to conservation planning and policy.     

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.