The study of hydrodynamic behaviour of a complex karst system under low-flow conditions using natural and artificial tracers (the catchment of the Unica River,SW Slovenia)

In this study a multi-tracer test with fluorescent tracers was combined with time series analyses of natural tracers to characterize the dynamics of the solute transport through different recharge pathways and to study hydraulic behaviour of a binary karst system under low-flow conditions. Fluorescent tracer testing included the introduction of uranine, amidorhodamine G, or naphthionate at three injection points. Sampling and monitoring took place at two karst springs (Malenščica, Unica) and at two underground rivers (Pivka, Rak) recharging the Unica River at the Polje of Planina, SW Slovenia. Other monitored parameters included precipitation, spring or underground river discharge, water temperature, and electrical conductivity. Water samples were collected and analyzed for total organic carbon, Mg²⁺, SO₄ ²⁻, and NO₃ ⁻ in the laboratory. In the study area, results of the tracer test suggest that contaminant transport in karst may be retarded for several weeks during low-flow conditions followed by increases in contaminant concentrations after subsequent rainfall events. Based on interpretation of tracer concentration breakthrough curves, low apparent dominant flow velocities (i.e., between 5.8 and 22.8 m/h through the well developed karst conduits, and 3.6 m/h through the prevailing vadose zone with a dominant influence of a diffuse recharge) were obtained. Together with analyses of hydro-chemographs the artificial tracing identified different origins of water recharging the studied aquifer. During prolonged low-water conditions the Malenščica spring is mainly recharged from the karst aquifer and the Unica spring by the sinking Pivka River. After more intensive rainfall events allogenic recharge from Cerknica prevails in the Malenščica spring, while the Unica spring drains mainly the allogenic water from the Pivka Valley. These findings of alternating hydraulic connections and drainage areas due to respective hydrological conditions are important and should be considered when monitoring water quality, implementing groundwater protection measures, and optimizing future water exploitation.     

Climate controls on rainfall isotopes and their effects on cave drip water and speleothem growth: the case of Molinos cave (Teruel,NE Spain)

The interpretation of stable isotopes in speleothems in terms of past temperature variability or precipitation rates requires a comprehensive understanding of the climatic factors and processes that influence the δ¹⁸O signal in the way through the atmosphere to the cave, where carbonate precipitates acquiring its final isotopic composition. This study presents for the first time in the Iberia Peninsula an integrated analysis of the isotopic composition of rainfall (δ¹⁸O_p) during 2010–2012 years and, through a detailed monitoring survey, the transference of the primary isotopic signal throughout the soil and epikarst into the Molinos cave (Teruel, NE Spain). Both air temperature and amount of precipitation have an important effect on δ¹⁸O_p values, clearly imprinting a seasonal variability modulated by an amount effect when rainfall events are more frequent or intense. Air mass history and atmospheric circulation influences are considered through the study of weather types, synoptic-scale climate patterns and large-scale atmospheric circulation indexes (North Atlantic Oscillation and Western Mediterranean Oscillation) revealing a dominant source effect on δ¹⁸O_p values in this region where tropical North Atlantic and Western Mediterranean are the two moisture source regions. A delay of 2–3 months occurs between the dripwater oxygen isotopic composition (δ¹⁸O_d) respect to δ¹⁸O_p values as a consequence of large residence time in the epikarst. Limited calcite precipitates are found from winter to spring when δ¹⁸O_d values are less negative and dripwater rates are constant. This study suggests that NE Iberian δ¹⁸O_calcite proxy records are best interpreted as reflecting a combination of parameters, not just paleotemperature or paleorainfall and, if extending present-day situation towards the recent past, a biased signal towards winter values should be expected in Molinos speleothem records.     

Solar influence and hydrological variability during the Holocene from a speleothem annual record (Molinos Cave,NE Spain)

We present a multi‐proxy approach to reconstructing Holocene climate conditions in northeastern Spain based on an excellent correlation among the lamina thickness, colour parameters and isotope (δ¹⁸O and δ¹³C) variations recorded in a speleothem. An age model constructed from five U/Th dates and annual lamina counting suggests that the uppermost 14.7 cm of the MO‐7 stalagmite grew between 7.2 and 2.5 ka before present but experienced a growth hiatus from 4.9 to 4.3 ka. Three spectral analysis methods were applied to 11 time series. The results reveal common solar periodicities on decennial (Gleissberg cycle) and centennial (De Vries‐Suess cycle) scales. The onset of Holocene carbonate precipitation in the MO‐7 stalagmite appears to be associated with a cold, wet period, whereas the hiatus and the end of growth are related to warm, dry periods. This environmental trend fits well within the regional Holocene climate.     

Hydrochemical dynamics of TOC and NO3 − contents as natural tracers of infiltration in karst aquifers

Analysis of Total Organic Carbon (TOC) and NO₃ ^? contents in the water at six springs in the province of Malaga (Southern Spain), sampled under different hydrologic conditions, revealed two different hydrochemical behaviour patterns for these natural tracers of infiltration. TOC content increased during every recharge period, following the rapid arrival of water infiltrating through the soil. On the other hand, NO₃ ^? content only rose during the first flood episodes (normally in autumn), and fell during the winter and spring ones. This difference is consequence of the distinctive biogeochemical kinetics of nitrogen with respect to organic carbon, both in the soil and within the aquifer. Unlike the mineralisation undergone by TOC from the surface to the spring, the NO₃ ^? ion remains in the aquifer almost unaffected, due to the oxidizing conditions prevailing within the karst medium, which do not allow its denitrification and favour its preservation within the saturated zone. In non-polluted aquifers, TOC and NO₃ ^? have a common origin in the soil and can be used to determine infiltration processes and the hydrogeological functioning of karst aquifers. Their different hydrochemical evolution provides information about mineralization and degradation processes of organic matter within karst aquifers, which can be used to validate the vulnerability to contamination in this type of medium.     

Imaging the Seismic Crustal Structure of the Western Mexican Margin between 19°N and 21°N

Three thousand kilometres of multichannel (MCS) and wide-angle seismic profiles, gravity and magnetic, multibeam bathymetry and backscatter data were recorded in the offshore area of the west coast of Mexico and the Gulf of California during the spring 1996 (CORTES survey). The seismic images obtained off Puerto Vallarta, Mexico, in the Jalisco subduction zone extend from the oceanic domain up to the continental shelf, and significantly improve the knowledge of the internal crustal structure of the subduction zone between the Rivera and North American (NA) Plates. Analyzing the crustal images, we differentiate: (1) An oceanic domain with an important variation in sediment thickness ranging from 2.5 to 1 km southwards; (2) an accretionary prism comprised of highly deformed sediments, extending for a maximum width of 15 km; (3) a deformed forearc basin domain which is 25 km wide in the northern section, and is not seen towards the south where the continental slope connects directly with the accretionary prism and trench, thus suggesting a different deformational process; and (4) a continental domain consisting of a continental slope and a mid slope terrace, with a bottom simulating reflector (BSR) identified in the first second of the MCS profiles. The existence of a developed accretionary prism suggests a subduction–accretion type tectonic regime. Detailed analysis of the seismic reflection data in the oceanic domain reveals high amplitude reflections at around 6 s [two way travel time (twtt)] that clearly define the subduction plane. At 2 s (twtt) depth we identify a strong reflection which we interpret as the Moho discontinuity. We have measured a mean dip angle of 7° ± 1° at the subduction zone where the Rivera Plate begins to subduct, with the dip angle gently increasing towards the south. The oceanic crust has a mean crustal thickness of 6.0–6.5 km. We also find evidence indicating that the Rivera Plate possibly subducts at very low angles beneath the Tres Marias Islands.