Semi-analytical solutions of a contaminant transport equation with nonlinear sorption in 1D

A new method to determine semi-analytical solutions of one-dimensional contaminant transport problem with nonlinear sorption is described. This method is based on operator splitting approach where the convective transport is solved exactly and the diffusive transport by finite volume method. The exact solutions for all sorption isotherms of Freundlich and Langmuir type are presented for the case of piecewise constant initial profile and zero diffusion. Very precise numerical results for transport with small diffusion can be obtained even for larger time steps (e.g., when the Courant-Friedrichs-Lewy (CFL) condition failed).     

Hydrology and hydrogeology of Sakaryabaşı Karstic springs, Çifteler, Turkey

The Sakarya River is one of the largest rivers in Turkey and is fed mainly from Sakaryabaşı springs. The Sakaryabaşı springs are located in the Central Anatolia and issue from confined/semi-confined karst having a thermal component and therefore, having quite different hydrogeological characteristics as compared to the Taurus Karst region, a typical example of the Mediterranean type of karst. The karstic carbonate rocks that form the groundwater reservoir are overlain by a thick semi-pervious overburden of mainly clastics of Neogene age. Tectonics is the major factor controlling the occurrence of the karst springs in the area where topography is rather flat. This study aimed at explaining the occurrence and movement of the karst groundwater within the system by use of hydrogeological, chemical, and isotopic tools. Isotopic composition of the waters revealed that all waters in the region are of meteoric origin and the thermal component is due to deep circulation. The catchment area of the hydrogeological system extends to the south and groundwater movement is towards the outlets, which are in a depression along a major fault. The movement of the groundwater, based on analysis of remotely sensed images, is controlled mainly by structural elements.     

GGR Critical Review of Analytical Developments in 2004–2005

In 2005 Geostandards and Geoanalytical Research embarked upon a new initiative for its readers. Key researchers in various fields of geoanalytical technique development and their application were identified and invited to provide reviews pertinent to their expertise. As noted in the first of these publications "…instead of revisiting the historical context or decades of development in each analytical technique, the goal here has been to capture a snapshot of "hot topics" across a range of fields as represented in the… literature" (Hergt et al . 2005). Rather than prepare an annual review, a decision was taken earlier this year to provide a biennial summary of progress and accomplishments, in this case for the years 2004–2005. The principal techniques employed in Earth and environmental sciences are covered here, and include laser ablation and multicollector ICP-MS, ICP-AES, thermal ionisation and secondary ion mass spectrometry, as well as neutron activation analysis, X-ray fluorescence and atomic absorption spectrometry. A comprehensive review of the development of reference materials, often essential to these techniques, is also provided. The contributions assembled serve both to keep readers informed of advances they may be unfamiliar with, but also as a means of showcasing examples of the breadth and depth of work being conducted in these fields.     

Reference Materials in Geoanalytical Research -Review for 2004 and 2005

This review gives an overview of the use and development of reference materials of geochemical and environmental interest in the literature of the years 2004 and 2005. In these years the performance of existing methods has been improved and new geochemical applications using new techniques have been developed. Accordingly, there was an increasing need for new reference materials, especially for in situ microanalysis and for precise stable isotope measurements. In addition, there was a notable trend for further characterisation of existing reference materials, mainly for the platinum-group elements. This review focuses on five topics: reference materials for platinum-group elements, reference glasses for in situ microanalysis, zircon reference materials, isotopic reference materials, and the development and certification of reference materials.     

The structure and dynamics of 2-dimensional fluids in swelling clays

The interlayer pores of swelling 2:1 clays provide an ideal 2-dimensional environment in which to study confined fluids. In this paper we discuss our understanding of the structure and dynamics of interlayer fluid species in expanded clays, based primarily on the outcome of recent molecular modelling and neutron scattering studies. Counterion solvation is compared with that measured in bulk solutions, and at a local level the cation-oxygen coordination is found to be remarkably similar in these two environments. However, for the monovalent ions the contribution to the first coordination shell from the clay surfaces increases with counterion radius. This gives rise to inner-sphere (surface) complexes in the case of potassium and caesium. In this context, the location of the negative clay surface charge (i.e. arising from octahedral or tetrahedral substitution) is also found to be of major importance. Divalent cations, such as calcium, eagerly solvate to form outer-sphere complexes. These complexes are able to pin adjacent clay layers together, and thereby prevent colloidal swelling. Confined water molecules form hydrogen bonds to each other and to the clays' surfaces. In this way their local environment relaxes to close to the bulk water structure within two molecular layers of the clay surface. Finally, we discuss the way in which the simple organic molecules methane, methanol and ethylene glycol behave in the interlayer region of hydrated clays. Quasi-elastic neutron scattering of isotopically labelled interlayer CH₃OD and (CH₂OD)₂ in deuterated clay allows us to measure the diffusion of the CH₃- and CH₂-groups in both clay and liquid environments. We find that in both the one-layer methanol solvates and the two-layer glycol solvates the diffusion of the most mobile organic molecules is close to that in the bulk solution.     

Application of the QUILF thermobarometer to the peralkaline trachytes and pantellerites of the Eburru volcanic complex, East African Rift, Kenya

The Quaternary Eburru volcanic complex in the south-central Kenya Rift consists of pantelleritic trachytes and pantellerites. The phenocryst assemblage in the trachytes is sanidine + fayalite + ferrohedenbergite + aenigmatite ± quartz ± ilmenite ± magnetite ± pyrrhotite ± pyrite. In the pantellerites, the assemblage is sanidine + quartz + ferrohedenbergite + fayalite + aenigmatite + ferrorichterite + pyrrhotite ± apatite, although fayalite, ferrohedenbergite and ilmenite are absent from more evolved rocks (e.g. with SiO₂ > 71%). QUILF temperature calculations for the trachytes range from 709 to 793 °C and for the pantellerites 668–708 °C, the latter temperatures being among the lowest recorded for peralkaline silicic magmas. The QUILF thermobarometer demonstrates that the Eburru magmas crystallized at relatively low oxidation states (ΔFMQ + 0.5 to − 1.6) for both trachytes and pantellerites. The trachytes and pantellerites evolved along separate liquid lines of descent, the trachytes possibly deriving from a more mafic parent by fractional crystallization and the pantellerites from extreme fractionation of comenditic magmas.     

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.

     

New Generation of Probabilistic Seismic Hazard Assessment for the Area Cologne/Aachen Considering the Uncertainties of the Input Data

A new earthquake catalogue for central, northern and northwestern Europe with unified M_w magnitudes, in part derived from chi-square maximum likelihood regressions, forms the basis for seismic hazard calculations for the Lower Rhine Embayment. Uncertainties in the various input parameters are introduced, a detailed seismic zonation is performed and a recently developed technique for maximum expected magnitude estimation is adopted and quantified. Applying the logic tree algorithm, resulting hazard values with error estimates are obtained as fractile curves (median, 16% and 84% fractiles and mean) plotted for pga (peak ground acceleration; median values for Cologne 0.7 and 1.2 m/s² for probabilities of exceedence of 10% and 2%, respectively, in 50 years), 0.4 s (0.8 and 1.5 m/s²) and 1.0 s (0.3 and 0.5 m/s²) pseudoacclerations, and intensity (I₀ = 6.5 and 7.2). For the ground motion parameters, rock foundation is assumed. For the area near Cologne and Aachen, maps show the median and 84% fractile hazard for 2% probability of exceedence in 50 years based on pga (maximum median value about 1.5 m/s²), and 0.4 s (>2 m/s²) and 1.0 s (about 0.8 m/s²) pseudoaccelerations, all for rock. The pga 84% fractile map also has a maximum value above 2 m/s² and shows similarities with the median map for 0.4 s. In all maps, the maximum values fall within the area 6.2–6.3° E and 50.8–50.9° N, i.e., east of Aachen.