Hydrogeochemistry and isotope analyses used to determine groundwater recharge and flow in low-gradient catchments of the province of Buenos Aires,Argentina

The region known as Pampa Plain, in Argentina, is a vast area characterized by slopes of less than 0.05%. The surface sediment is silty sand, mainly Aeolian, referred to as Pampean Loess, which is a phreatic aquifer unit of utmost importance for the water supply of the region. On account of the slight gradient, hydrogeological analyses using only hydraulic measurements are difficult to perform, often leading to confusing results. Thus, the study presented relies on hydrochemical modeling and isotopic determinations as well. The study area comprises three catchments in the inter-mountainous area corresponding to El Moro, Tamangueyú and Seco creeks, and covering 2,570 km² in the province of Buenos Aires. Measurements of piezometric levels and samples of groundwater, surface water and rainwater were carried out between January and March 2005. Major ion results were analyzed by means of hydrochemical graphs and hydrogeochemical modelling using NETPATH. The resulting data show that it is possible to identify local changes in recharge, flow direction and stream/groundwater relationship by using hydrochemical and isotopic information, which may become a useful and more precise tool for the study of particularly flat landscapes.     

The Meso-Cenozoic thermo-tectonic evolution of the Eastern Pyrenees: an <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar fission track and (U–Th)/He thermochronological study of the Canigou and Mont-Louis massifs

Vertical displacements on the SW–NE Têt fault (Eastern Pyrenees Axial Zone, France), which separates the Variscan Canigou-Carança and Mont-Louis massifs, were constrained using a thermochronologic multi-method approach. ⁴⁰Ar/³⁹Ar data from the granitic Mont-Louis massif record its Variscan cooling history and reveal no ages younger than Early Cretaceous, while the Canigou-Carança gneiss massif records systematically younger ⁴⁰Ar/³⁹Ar ages. These younger ⁴⁰Ar/³⁹Ar ages in the Canigou-Carança gneiss massif are the result of partial to total rejuvenation of argon isotopic systems related to a thermal flow coeval with the Cretaceous HT-BP metamorphism in the North Pyrenean Zone. Only the deepest rocks from the Canigou-Carança suffered this extensive Mid-Cretaceous thermal overprint probably due to differential burial around 4 km at that time. The post Mid-Cretaceous vertical displacements along the Têt fault are recorded by “low” temperature thermochronology using K-feldspar ⁴⁰Ar/³⁹Ar, zircon and apatite fission track and (U–Th)/He datings. The Mont-Louis granite samples experienced a long period of protracted cooling reflecting a lack of thermo-tectonic activity in this area from Late Palaeozoic to Early Cenozoic, followed by cooling from 55–60 Ma to Late Eocene at a mean rate of 15–20°C/Ma in the final stage. This cooling stage corresponds to Têt fault reactivation with a reversed component, promoting exhumation of the Mont-Louis roof zone contemporaneously with the south-vergent Pyrenean thrusting. In the Canigou-Carança massif, the main cooling event occurred from 32 to 18 Ma at a maximum rate of 30°C/Ma during Early Oligocene followed by a more moderate rate of 3°C/Ma from Late Oligocene to Early Burdigalian, coeval with the normal reactivation of the Têt fault in brittle conditions that accommodated the final exhumation of the massif during the opening of the Gulf of Lion.     

Provenance of late Palaeozoic metasediments of the Patagonian proto-Pacific margin (southernmost Chile and Argentina)

In this provenance study of late Palaeozoic metasediments of the Eastern Andean Metamorphic Complex (EAMC) along the south Patagonian proto-Pacific margin of Gondwana, the palaeogeological setting of the continental margin in Devonian–Carboniferous and Permian times is reconstructed. The study is based on detrital heavy mineral contents, chemical compositions of tourmaline grains, and whole rock element and Nd-Sr isotopic compositions. Element and isotopic compositions reveal that Devonian–Carboniferous metaturbidites deposited before the development of a Late Carboniferous–Permian magmatic arc along the margin were mainly fed from felsic, recycled, old continental rocks. The last recycling phase involved erosion of metasediments that were exposed in Patagonia. Feeder systems to the basin cut either through epidote-rich or garnet-rich metasediments. In Permian time, EAMC metaturbidites were deposited next to the evolving magmatic arc and were derived from felsic, crustal rocks. Two provenance domains are recognised. The metasediments of the northern one are chemically similar to those of the Devonian–Carboniferous metasediments. This domain was fed from the metasedimentary host rocks of the magmatic arc. The southern domain probably was fed from the arc proper, as indicated mainly by the dominance of metaplutonic lithic fragments, abundant detrital biotite, and the major element composition of the metasediments.     

High heat flow anomalies on an old oceanic plate observed seaward of the Japan Trench

Thirty-three new measurements on the seaward slope and outer rise of the Japan Trench along a parallel of 38°45′N revealed the existence of high heat flow anomalies on the subducting Pacific plate, where the seafloor age is about 135 m.y.. The most prominent anomaly with the highest value of 114 mW/m² is associated with a small mound on the outer rise, which was reported to be a kind of mud volcano. On the seaward slope of the trench, heat flow is variable: high (70–90 mW/m²) at some locations and normal for the seafloor age (about 50 mW/m²) at others. The spatial variation of heat flow may be related to development of normal faults and horst/graben structures due to bending of the Pacific plate before subduction, with fluid flow along the fault zones enhancing the vertical heat transfer. Possible heat sources of the high heat flow anomalies are intra-plate volcanism in the last several million years like that discovered recently on the Pacific plate east of the Japan Trench.     

Estimates of heat flow and heat production and a thermal model of the São Francisco craton

An updated analysis of geothermal data from the highland area of eastern Brazil has been carried out and the characteristics of regional variations in geothermal gradients and heat flow examined. The database employed includes results of geothermal measurements at 45 localities. The results indicate that the Salvador craton and the adjacent metamorphic fold belts northeastern parts of the study area are characterized by geothermal gradients in the range of 6–17°C/km. The estimated heat flow values fall in the range of 28–53 mW/m², with low values in the cratonic area relative to the fold belts. On the other hand, the São Francisco craton and the intracratonic São Francisco sedimentary basin in the southwestern parts are characterized by relatively higher gradient values, in the range of 14–42°C/km, with the corresponding heat flow values falling in the range of 36–89 mW/m². Maps of regional variations indicate that high heat flow anomaly in the São Francisco craton is limited to areas of sedimentary cover, to the west of the Espinhaço mountain belt. Crustal thermal models have been developed to examine the implications of the observed intracratonic variations in heat flow. The thermal models take into consideration variation of thermal conductivity with temperature as well as change of radiogenic heat generation with depth. Vertical distributions of seismic velocities were used in obtaining estimates of radiogenic heat production in crustal layers. Crustal temperatures are calculated based on a procedure that makes simultaneous use of the Kirchoff and Generalized Integral Transforms, providing thereby analytical solutions in 2D and 3D geometry. The results point to temperature variations of up to 300°C at the Moho depth, between the northern Salvador and southern São Francisco cratons. There are indications that differences in rheological properties, related to thermal field, are responsible for the contrasting styles of deformation patterns in the adjacent metamorphic fold belts.     

Geochemistry of amphibolitized eclogites and cross-cutting tonalitic–trondhjemitic dykes in the Metamorphic Kimi Complex in East Rhodope (N.E. Greece): implications for partial melting at the base of a thickened crust

In the ultra-high pressure Metamorphic Kimi Complex widespread tonalitic–trondhjemitic dykes, with an intrusion age ca. 65–63 Ma, cross-cut boudins and layers of amphibolitized eclogites. Geochemical investigation proclaims the tied genetic relationship of the amphibolitized eclogites and the associated tonalitic–trondhjemitic dykes. The major and trace element contents and rare earth element patterns of the amphibolitized eclogites indicate formation of their protoliths by fractional crystallization of tholeiitic magmas in a back-arc environment. The tonalites and trondhjemites are characterized by moderate to high Sr contents (>130 ppm), and low Y (<8.2 ppm) and heavy rare earth element contents (Yb content of 0.19–0.88 ppm). The chemical composition of the tonalitic and trondhjemitic dykes are best explained by partial melting of a tholeiitic source like the amphibolitized eclogites with residual garnet and amphibole, at the base of a thickened crust during Early Tertiary subduction/accretion at the southern margins of the European continent.     

Evaluation of change in subsurface thermal environment due to groundwater flow in the Tokyo Lowland,Japan

Information on the distribution of subsurface temperature and hydraulic heads at 24 observation wells in and around the Tokyo Lowland, the eastern part of the Tokyo Metropolitan area, were examined to make clear the relationship between groundwater and the subsurface thermal environment in the urban area. Minimums in temperature–depth logs due to subsurface temperature increasing at shallow parts were recognized in 21 wells. This fact shows subsurface temperature is affected by ground surface warming in almost all of this area. Deeper than minimums, where the effects of surface warming became relatively small, regional variation is observed as follows: high temperatures are shown in the central part to the southern part, and low temperatures shown in the inland to eastern part. The high temperature area corresponds to an area where the lower boundary of groundwater flow is relatively shallow. This area corresponds also to an area with severe land subsidence resulting from excessive groundwater pumping. It is considered that this high temperature area is formed by the effects of upward groundwater flow affected by hydrogeological conditions and pumping. On the other hand, a comparison between past data (1956–1967) and present data (2001–2003) revealed widespread decreasing temperature in the inland area. This is explained by downward groundwater flow based on an analysis of temperature–depth logs. This fact suggests that subsurface temperature is not only increasing from the effects of surface warming but also decreasing from the effects of groundwater environment change due to pumping.     

Surface water–groundwater interaction and chemistry in a mineral-armored hydrothermal outflow channel, Yellowstone National Park, USA

Small quantities of groundwater interact with hydrothermal surface water to drive in-stream geochemical processes in a silica-armored hot-spring outflow channel in Yellowstone National Park, USA. The objective of this study was to characterize the hydrology and geochemistry of this unique system in order to (1) learn more about the Yellowstone Plateau’s subsurface water mixing between meteoric and hydrothermal waters and (2) learn more about the chemical and physical processes that lead to accumulation of streambed cements, i.e., streambed armor. A combination of hydrological, geochemical, mineralogical, microscopic, and petrographic techniques were used to identify groundwater and surface-water exchange. Interaction could be identified in winter because of differences in surface water and groundwater composition but interaction at other times of the year cannot be ruled out. Dissolved constituents originating from groundwater (e.g., Fe(II) and Mg) were traced downstream until oxidation and/or subsequent precipitation with silica removed them, particularly where high affinity substrates like cyanobacterial surfaces were present. Because the stream lies in a relatively flat drainage basin and is fed mainly by a seasonally relatively stable hot spring, this system allowed study of the chemical processes along a stream without the obscuring effects of sedimentation.     

Field scale heterogeneity of redox conditions in till-upscaling to a catchment nitrate model

Point scale studies in different settings of glacial geology show a large local variation of redox conditions. There is a need to develop an upscaling methodology for catchment scale models. This paper describes a study of field-scale heterogeneity of redox-interfaces in a till aquitard within an area of 600?×?600 m. The results showed significant variation of the depths to the redox-interface and thickness of the aquitard. Nitrate was present above the redox-interface but reduced to non-detectable levels a few metres below the interface. An upscaling approach for an area of 92 km² is proposed. Two models are proposed to predict the depth to the redox-interface in the aquitard and the resulting nitrate recharge concentrations to an underlying aquifer. The first model assumes that the depth to the redox-interface reflects the hydraulic head in the aquitard, and the second model assumes that the depth of the redox-interface is randomly distributed according to a log-normal probability distribution function. The upscaling approach using the random redox model estimated recharge concentrations comparable to the observed concentration in the underlying aquifer. The presented upscaling approach is applicable in distributed catchment models where sub-grid variability cannot be represented by the large grids.