A conceptual hydrogeological model of ophiolite hard-rock aquifers in Oman based on a multiscale and a multidisciplinary approach

Ophiolites are found all over the world: from the Alps to the Himalayas, in Cuba, Papua-New Guinea, New Caledonia, Newfoundland, etc. They are composed of hard rocks—basalt, dolerite, gabbro and peridotite, which are formed at the mid-oceanic ridges, with specific ridge-related tectonic fracturing and intense hydrothermal alteration. Their geological and thus their hydrogeological properties differ from those of both granite or classical gabbro and classical basaltic lava. A conceptual hydrogeological model of these hard-rock aquifers was developed based on the convergent results of a multidisciplinary approach at several spatial scales, from rock-sample (centimetre) to catchment (kilometre), on well-preserved ophiolite rocks in Oman. In ophiolite rocks, groundwater circulation takes place mostly in the fissured near-surface horizon (50 m thick), and, to a lesser degree, in the tectonic fractures. Hydrograph analysis (Water Resour Res 34:233–240, 1977), interpretation of numerous pumping tests using both classical Theis and dual porosity models [Water Resour Res 32:2733–2745, 1996; Comput Geosci J (in press)], and mercury porosity and hydraulic conductivity lab-measurements support the aquifer parameter estimates. The hydraulic conductivity K of the fissured horizon is estimated at 10^–5 to 10^–6 m/s for gabbro and dolerite, and 10^–7 m/s for peridotite. The storage coefficient S of the peridotite aquifer is estimated at 10^–3 and appears to be controlled mainly by microcracks (20 to 100 m wide). Tectonic fractures in the ophiolite have similar hydrodynamic properties regardless of lithology (10^–1<T<10^–4 m²/s and 10^–1<S<10^–3) though the probability of obtaining productive wells is two to three times greater in gabbro and dolerite than in peridotite. Some of the tectonic fractures produce small hydrothermal, hyperalkaline springs in the peridotite. The water budget and hydrochemistry of the Oman ophiolite are characterized and support the conceptual hydrogeological model. Despite low annual rainfall, a relatively low hydraulic conductivity and a significant storage coefficient explain why most of the streams in peridotite are perennial.

---

Résumé Les ophiolites se trouvent partout dans le monde : des Alpes à lHymalaya, à Cuba, Papouasie Nouvelle-Guinée, Nouvelle Calédonie, etc. Les ophiolites sont composées de roches dures basaltiques, de dolérites, gabbros et péridotites, qui sont formées aux rides mi-océaniques, avec, en relation avec les rides, une fracturation tectonique et une altération hydrothermale intense. Leurs propriétés géologiques et hydrogéologiques diffèrent de celles des granites ou des gabbros classiques ou des laves basaltiques classiques. Un modèle hydrogéologique conceptuel des aquifères de hard-rock a été développé sur les résultats convergents dune approche multidisciplinaire à différentes échelles spatiales, de léchantillon de roche de quelques centimètres à la taille du bassin versant en kilomètres, des ophiolites bien conservées dOman. Dans les roches ophiolitiques, la circulation des eaux souterraines est surtout localisée dans les fissures dun horizon de surface de plus dune cinquantaine de mètres et à un degré moindre dans les fractures tectoniques. Lanalyse des hydrographes (Water Resour Res 34:233–240, 1977), linterprétation des nombreux essais de pompage utilisant des modèles classiques de Theis et des modèles à double porosité (Water Resour Res 32:2733–2745, 1996; Comput Geosci J (Soumis), la porosité au mercure et la conductivité hydraulique calculée en laboratoire supportent lestimation des paramètres de laquifère. La conductivité hydraulique K de lhorizon fissuré est estimée entre 10^–5 et 10^–6 m/s pour les gabbros et les dolérites, et 10^–7 m/s pour les péridotites. Le coefficient demmagasinement S de laquifère de péridotite est estimé à 10^–3 et semble être contrôlé par les micro-fractures (20 à 100 m de largeur). Les fractures tectoniques dans les ophiolites possèdent des propriétés hydrodynamiques en regard de la lithologie (10-1<T<10-4 m²/s et 10-1/S/10-3) bien que la probabilité dobtenir des puits productifs est trois plus grand dans le grabbro et la dolérite que la péridétite. Quelques unes des fractures tectoniques produisent de petites sources hydrothermales hyperalcalines dans la péridotite. Le bilan hydrologique et hydrochimique de lophiolite dOman caractérise le modèle hydrogéologique conceptuel. Malgré des faibles précipitations annuelles, une relativement basse conductivité hydraulique et un coefficient demmagasinement significatif expliquent pourquoi la plus part des cours deau sur les péridotites ne sont jamais taris.

---

Resumen Las rocas ofiolitas se encuentran en muchas partes del mundo desde los Alpes a los Himilaya, Newfoundland, Nueva Caledonia etc. Son rocas duras compuestas de basalto, dolerita, gabro y peridotita que son formadas por los cordones medio-oceánicos y por fracturación relacionada con los cordones tectónicos específicos y alteración hidrotermal intensa. La geología y los propiedades hidrogeológicas varian entre las de granito y gabro clásico y lava de basalto clásico. Se desarrolló un modelo conceptual de los acuíferos de las rocas duras Basado en los resultados convergentes de un enfoce multidisciplinario a varias escalas a partir de muestras de roca (centímetro) a cuenca kilómetro) enfocado en las rocas ofilíticas bien preservados en Oman. En las rocas ofiolíticas la circulación de agua ocurre principalmente en el horizonte de fisuras cercanas a la superficie (>50 m de espesor) y en un grado menos importante en fracturas tectónicas. Los estimados de los parámetros de estos acuíferos se basan en análisis de hidrógrafos (Water Resour Res 34:233–240, 1977), en la interpretación de pruebas de bombeo usando los modelos clásicos de Theis y de porosidad doble (Hamm y Bidaux 1996; Lods y Gouze, [en prensa]), y porosidad de mecurio y conductividad hidraúlica medida en el laboratorio. La conductividad hidraúlica K del horizonte fisurado se estima en 10E-5 a 10E-6 m-s para gabro y dolerita y 10E-7 para peridotita. El coeficiente de abastecimiento S del acuífero peridotita se estima en 10–3 y parece ser controlado principalmente por microfracturas (20 to 100 m de ancho). Las fracturas tectónicas en la ofiolita tienen propiedades hidrodinámicas similares independientes de la litología (10–1<T<10–4 m²/s and 10–1<S<10–3) pero la probabilidad de obtener pozos productivos es dos o tres veces mayor en gabro que en dolerita o peridotita. Algunas de las fracturas tectónicas producen pequeñas descargas de agua hidrotermal hiperalcalina en la peridotita. El presupuesto e hidroquímica de las ofilitas de Omán se caracteriza y apoya en el modelo conceptual hidrogeológico. La baja lluvia anual, una conductividad hidraúlica relativamente baja y un coeficiente de abastecimiento importante explican porqué la mejoría de los drenajes en peridotita es peremne.

     

The solubility of nantokite (CuCl(s)) and Cu speciation in low-density fluids near the critical isochore: An in-situ XAS study

The solubility of nantokite (CuCl_(s)) and the structure of the predominant copper species in supercritical water (290-400 bar at 420 °C; 350-450 °C at 290 bar; 500 °C at 350 bar; density = 0.14-0.65 g/cm³) were investigated concurrently using synchrotron X-ray absorption spectroscopy (XAS) techniques. These conditions were chosen as they represent single phase solutions near the critical isochore, where the fluid density is intermediate of typical values for vapour and brine and is highly sensitive to even small changes in pressure. X-ray absorption near edge spectroscopy (XANES) and extended X-ray absorption spectroscopy (EXAFS) analyses show that aqueous copper occurs in a slightly distorted linear coordination in the solutions studied, with an average of 1.35(±0.3) Cl and 0.65(±0.3) O neighbours. The solubility of CuCl_(s) decreases exponentially with decreasing water density (i.e., decreasing pressure at constant temperature), in a manner similar to the solubility behaviour of salts such as NaCl in water vapour. Based on this similarity, an apparent equilibrium constant for the dissolution reaction of 0.5 ± 0.4 was calculated from a regression of the data at 420 °C, and it was determined that each Cu atom is solvated by approximately three water molecules. This indicates that under these conditions, copper solubility is controlled mainly by the structure of the second-shell hydration, which is essentially invisible to the XAS techniques used in this study.These results demonstrate that for a supercritical fluid near the critical isochore, decreasing pressure may initiate precipitation of copper even before boiling or phase separation. Such a process could be responsible for near-surface ore deposition in seafloor hydrothermal systems, where supercritical fluids experience rapid pressure changes during the transition between lithostatic and hydrostatic domains.     

In situ LA-ICP-MS U–Pb titanite dating of Na–Ca metasomatism in orogenic belts: the North Pyrenean example

In the Pyrenees, in association with the rotation of the Iberian plate around Europe during the Mid-Cretaceous, a Na–Ca metasomatism is recognized as a complementary record of the hydrothermal activity that led to Na-metasomatism (albitization) and talc–chlorite mineralization. It affected metasedimentary rocks as well as Hercynian granites. In situ laser ablation ICP-MS U–Pb analyses of titanite grains formed in albitites during metasomatism date the Na–Ca metasomatism between 110 and 92 Ma. The temperature of the Na–Ca metasomatism is estimated to be approximately 550 °C. Both the time constraints and temperature estimates suggest that the Na–Ca metasomatism is related to the low-P high-T North Pyrenean metamorphism.     

The Montalet granite, Montagne Noire, France: An Early Permian syn-extensional pluton as evidenced by new U-Th-Pb data on zircon and monazite

Dating the magmatism in the Montagne Noire gneiss dome in the southern French Massif Central is a key point for understanding the Late Palaeozoic evolution of this part of the Variscan belt, which is characterised by compressive tectonics during the Carboniferous and extensional tectonics during Stephanian-Permian times. The Montalet granite crops out in the north-western part of the dome and was first considered as an early syntectonic intrusion related to compressive deformation. More recently, it has been dated at 327 Ma and considered as contemporaneous with the diapiric ascent of the Montagne Noire gneiss dome before the Stephanian-Permian extension. We show that in fact, this pluton was emplaced 294 ± 1 Ma ago and is therefore contemporaneous with the Stephanian-Permian extension. This age is consistent with the interpretation of the Montagne Noire Massif as an extensional gneiss dome.     

Miocene to Recent exhumation of the central Himalaya determined from combined detrital zircon fission-track and U/Pb analysis of Siwalik sediments, western Nepal

Fission‐track (FT) analysis of detrital zircon from synorogenic sediment is a well‐established tool to examine the cooling and exhumation history of convergent mountain belts, but has so far not been used to determine the long‐term evolution of the central Himalaya. This study presents FT analysis of detrital zircon from 22 sandstone and modern sediment samples that were collected along three stratigraphic sections within the Miocene to Pliocene Siwalik Group, and from modern rivers, in western and central Nepal. The results provide evidence for widespread cooling in the Nepalese Himalaya at about 16.0±1.4 Ma, and continuous exhumation at a rate of about 1.4±0.2 km Myr^?1 thereafter. The ～16 Ma cooling is likely related to a combination of tectonic and erosional activity, including movement on the Main Central thrust and Southern Tibetan Detachment system, as well as emplacement of the Ramgarh thrust on Lesser Himalayan sedimentary and meta‐sedimentary units. The continuous exhumation signal following the ～16 Ma cooling event is seen in connection with ongoing tectonic uplift, river incision and erosion of lower Lesser Himalayan rocks exposed below the MCT and Higher Himalayan rocks in the hanging wall of the MCT, controlled by orographic precipitation and crustal extrusion. Provenance analysis, to distinguish between Higher Himalayan and Lesser Himalayan zircon sources, is based on double dating of individual zircons with the FT and U/Pb methods. Zircons with pre‐Himalayan FT cooling ages may be derived from either nonmetamorphic parts of the Tethyan sedimentary succession or Higher Himalayan protolith that formerly covered the Dadeldhura and Ramgarh thrust sheets, but that have been removed by erosion. Both the Higher and Lesser Himalaya appear to be sources for the zircons that record either ～16 Ma cooling or the continuous exhumation afterwards.     

Impacts of greenhouse gases and aerosol direct and indirect effects on clouds and radiation in atmospheric GCM simulations of the 1930–1989 period

Among anthropogenic perturbations of the Earths atmosphere, greenhouse gases and aerosols are considered to have a major impact on the energy budget through their impact on radiative fluxes. We use three ensembles of simulations with the LMDZ general circulation model to investigate the radiative impacts of five species of greenhouse gases (CO₂, CH₄, N₂O, CFC-11 and CFC-12) and sulfate aerosols for the period 1930–1989. Since our focus is on the atmospheric changes in clouds and radiation from greenhouse gases and aerosols, we prescribed sea-surface temperatures in these simulations. Besides the direct impact on radiation through the greenhouse effect and scattering of sunlight by aerosols, strong radiative impacts of both perturbations through changes in cloudiness are analysed. The increase in greenhouse gas concentration leads to a reduction of clouds at all atmospheric levels, thus decreasing the total greenhouse effect in the longwave spectrum and increasing absorption of solar radiation by reduction of cloud albedo. Increasing anthropogenic aerosol burden results in a decrease in high-level cloud cover through a cooling of the atmosphere, and an increase in the low-level cloud cover through the second aerosol indirect effect. The trend in low-level cloud lifetime due to aerosols is quantified to 0.5 min day^–1 decade^–1 for the simulation period. The different changes in high (decrease) and low-level (increase) cloudiness due to the response of cloud processes to aerosols impact shortwave radiation in a contrariwise manner, and the net effect is slightly positive. The total aerosol effect including the aerosol direct and first indirect effects remains strongly negative.     

Anticipation of ENSO: what teach us the resonantly forced baroclinic waves

The El Niño-Southern Oscillation (ENSO) is characterized based on the date the events are mature. Their time lag defined relative to the central value of successive intervals of 4 years length, e.g. 01/1868–01/1872, 01/1872–01/1876 …, 01/1996–01/2000 … affects their evolution and, for a given amplitude, their variability. It specifies the dynamics of the quadrennial Quasi-Stationary Wave (QSW) in the tropical Pacific since ENSO always occurs at the end of the eastward phase propagation of that QSW. A third of events are unlagged with very low variability, SST anomalies being nearly concomitant between the extreme eastern and the central-eastern Pacific. A third of events are weakly lagged, in phase with the annual QSW, whose variability is much greater. Several months may elapse between the maximum SST anomalies east of the basin and along its equatorial central part. The last third of the events exhibits considerable variability, whether they are out of phase with the annual QSW or strongly lagged. The SST anomaly between 5°N and 20°N plays a key role in the maturation of the events out of phase. The events in phase (10% of the total population) are characterized by a negative SST anomaly in the central-eastern Pacific that reverses just before the maturation stage of ENSO. Sea water temperature 125 m deep in the central-eastern Pacific carries the earliest relevant information with a lead time of one year for forecasting the amplitude of unlagged ENSO while reporting how SST anomalies will develop until ENSO is fully developed. Besides, long-term forecast of the resumption of La Niña can be achieved accurately when weakly lagged events in phase with the annual QSW occur. The well differentiated typology of events vs. their time lag is the best clue to prove the leading role of the quadrennial QSW in the genesis of ENSO, while the related dynamic of the atmosphere ensues.     

Modelling the Mixing Function to Constrain Coseismic Hydrochemical Effects: An Example from the French Pyrénées

Groundwater coseismic transient anomalies are evidenced and characterized by modelling the mixing function F characteristic of the groundwater dynamics in the Ogeu (western French Pyrénées) seismic context. Investigations of water-rock interactions at Ogeu indicate that these mineral waters from sedimentary environments result from the mixing of deep waters with evaporitic signature with surficial karstic waters. A 3-year hydrochemical monitoring of Ogeu springwater evidences that using arbitrary thresholds constituted by the mean ± 1 or 2σ, as often performed in such studies, is not a suitable approach to characterize transient anomalies. Instead, we have used a mixing function F calculated with chemical elements, which display a conservative behavior not controlled by the precipitation of a mineral phase. F is processed with seismic energy release (E_s) and effective rainfalls (R). Linear impulse responses of F to E_s and R have been calculated. Rapid responses (10 days) to rainwater inputs are evidenced, consisting in the recharge of the shallow karstic reservoir by fresh water. Complex impulse response of F to microseismic activity is also evidenced. It consists in a 2-phase hydrologic signal, with an inflow of saline water in the shallow reservoir with a response delay of 10 days, followed by an inflow of karstic water with a response delay of 70 days, the amount being higher than the saline inflow. Such a process probably results from changes in volumetric strain with subsequent microfracturation transient episodes allowing short inflow of deep salted water in the aquifer. This study demonstrates that groundwater systems in such environments are unstable systems that are highly sensitive to both rainfall inputs and microseismic activity. Impulse responses calculation of F to E_s is shown to be a powerful tool to identify transient anomalies. Similar processing is suggested to be potentially efficient to detect precursors of earthquakes when long time-series (5 years at least) are available in areas with high seismicity.     

Diagenetic trends of a tertiary low-rank coal series

The Mahakam delta (Kalimantan, Indonesia) coals represent all the evolution stages between freshly-deposited plant/peat material, lignites and bituminous coals. The geochemical techniques used to study this coal series included elemental analysis, extraction of humic compounds, infrared spectroscopy and ¹³C nuclear magnetic resonance of the total coal.The main mechanisms of early maturation in this series are loss of oxygenated compounds, aromatisation and condensation of the organic matter. These changes, which have already been suggested for other coal series and partially reported for sedimentary organic matter, were confirmed and described in more detail for the Mahakam coal series.