Hydrogeological investigations in the upper Ouémé catchment in Benin, West Africa

The scope of this PhD thesis was the hydrogeological conceptualisation of the Upper Ouémé river catchment in Benin. The study area exceeds 14,500 km**2 and is underlain by a crystalline basement. At this setting the typical sequence of aquifers - a regolith aquifer at the top and a fractured bedrock aquifer at the bottom - is encountered, which is found in basement areas all over Africa and elsewhere in the world. The chosen regional approach revealed important information about the hydrochemistry and hydrogeology of this catchment. Based on the regional conceptual model a numerical groundwater flow model was designed. The numerical model was used to estimate the impact of climate change on the regional groundwater resources. This study was realised within the framework of the German interdisciplinary research project IMPETUS (English translation: "Integrated approach to the efficient management of scarce water resources in West Africa"), which is jointly managed by the German universities of Bonn and Cologne. Since the year 2000 the Upper Ouémé catchment was the principal target for investigations into the relevant processes of the regional water cycle. A first study from 2000 to 2003 (Fass, 2004, http://nbn-resolving.de/urn:nbn:de:hbz:5n-03849) focused on the hydrogeology of a small local catchment (~30 km**2). In the course of this thesis five field campaigns were underdone from the year 2004 to 2006. In the beginning of 2004 a groundwater monitoring net was installed based on 12 automatic data loggers. Manual piezometric measurements and the sampling of groundwater and surface water were realised for each campaign throughout the whole study area. Water samples were analysed for major ions, for a choice of heavy metals and for their composition by deuterium, oxygen-18 and tritium. The numerical model was performed with FEFLOW. The hydraulic and hydrochemical characteristics were described for the regolith aquifer and the bedrock aquifer. The regolith aquifer plays the role of the groundwater stock with low conductivity while the fractures of the bedrock may conduct water relatively fast towards extraction points. Flow in fractures of the bedrock depends on the connectivity of the fracture network which might be of local to subregional importance. Stable isotopes in combination with hydrochemistry proved that recharge occurs on catchment scale and exclusively by precipitation. Influx of groundwater from distant areas along dominant structures like the Kandi fault or from the Atacora mountain chain is excluded. The analysis of tritium in groundwater from different depths revealed the interesting fact of the strongly rising groundwater ages. Bedrock groundwater may possibly be much older than 50 years. Equilibrium phases of the silicate weathering products kaolinite and montmorillonite showed that the deeper part of the regolith aquifer and the bedrock aquifer feature either stagnant or less mobile groundwater while the shallow aquifer level is influenced by seasonal groundwater table fluctuations. The hydrochemical data characterised this zone by the progressive change of the hydrochemical facies of recently infiltrated rainwater on its flow path into deeper parts of the aquifers. Surprisingly it was found out that seasonal influences on groundwater hydrochemistry are minor, mainly because they affect only the groundwater levels close to the surface. The transfer of the hydrogeological features of the Upper Ouémé catchment into a regional numerical model demanded a strong simplification. Groundwater tables are a reprint of the general surface morphology. Pumping or other types of groundwater extraction would have only very local impact on the available groundwater resources. It was possible to integrate IMPETUS scenario data into the groundwater model. As a result it was shown that the impact of climate change on the groundwater resources until the year 2025 under the given conditions will be negligible due to the little share of precipitation needed for recharge and the low water needs for domestic use. Reason for concern is the groundwater quality on water points in the vicinity of settlements because of contamination by human activities as shown for the village of Dogué. Nitrate concentrations achieved in many places already alerting levels. Health risks from fluoride or heavy metals were excluded for the Upper Ouémé area.

Interstitial water chemistry at DSDP Leg 57 Holes

Analyses of water samples taken by means of an in-hole sampler generally show good agreement with analyses of samples collected by routine shipboard squeezing techniques. At Sites 438 and 439, a decrease in salinity with depth is related to former freshwater flow from an aquifer that crops out at an anticline on a deep sea terrace between Japan and the top of the trench slope of the Japan Trench. This former subaerial recharge suggests significant late Cenozoic subsidence of the terrace, because it now lies at a water depth of 1500 meters. Samples from the trench slope at Site 440 have extremely high values of alkalinity and ammonia, presumably because of a favorable combination of high sedimentation rate and organic carbon content. Diagenetic conditions on the trench slope favor formation of the Fe-Mg carbonate mineral, ankerite; at Site 440 it first occurs at a depth below the sea floor of only 29 meters in late Pleistocene strata. Undissolved diatoms persist to relatively great depth at the sites of Leg 57 because of a low geothermal gradient caused by subduction. Secondary silica lepispheres first appear at 851 meters at the most landward and warmest site, Site 438, in strata 16 million years old with an ambient temperature of 31 °C.

Hydrogen and oxygen stable isotope analyses of precipitation sampled at Erlangen (southern Germany) between 2010 and 2013

Shallow groundwater aquifers are often influenced by anthropogenic contaminants or increased nutrient levels. In contrast, deeper aquifers hold potentially pristine paleo-waters that are not influenced by modern recharge. They thus represent important water resources, but their recharge history is often unknown. In this study groundwater from two aquifers in southern Germany were analyzed for their hydrogen and oxygen stable isotope compositions. One sampling campaign targeted the upper aquifer that is actively recharged by modern precipitation, whereas the second campaign sampled the confined, deep Benkersandstein aquifer. The groundwater samples from both aquifers were compared to the local meteoric water line to investigate sources and conditions of groundwater recharge. In addition, the deep groundwater was dated by tritium and radiocarbon analyses. Stable and radiogenic isotope data indicate that the deep-aquifer groundwater was not part of the hydrological water cycle in the recent human history. The results show that the groundwater is older than ~20,000 years and most likely originates from isotopically depleted melt waters of the Pleistocene ice age. Today, the use of this aquifer is strictly regulated to preserve the pristine water. Clear identification of such non-renewable paleo-waters by means of isotope geochemistry will help local water authorities to enact and justify measures for conservation of these valuable resources for future generations in the context of a sustainable water management.

Biogeochemistry and microbiology of groundwater during acetate and bicarbonate amendment to an alluvial aquifer

These data are from a field experiment conducted in a shallow alluvial aquifer along the Colorado River in Rifle, Colorado, USA. In this experiment, bicarbonate-promoted uranium desorption and acetate amendment were combined and compared to an acetate amendment-only experiment in the same experimental plot. Data include names and location data for boreholes, geochemical data for all the boreholes between June 1, 2010 and January 1, 2011, microarray data provided as signal to noise ratio (SNR) for individual microarray probes, microarray data provided as signal to noise ratio (SNR) by Genus.

Análisis de la técnica de recarga artificial en la Cubeta de Sant Andreu de la Barca (Barcelona)

El presente proyecto pretende mostrar las posibilidades de la recarga artificial como elemento de gestión de los recursos de agua subterránea del acuífero aluvial del río Llobregat, en su sector de la Cubeta de Sant Andreu de la Barca (Barcelona), el cual es fuente de abastecimiento urbano, agrícola e industrial del entorno de Sant Andreu de la Barca. Igualmente, se revisa la efectividad de la aplicación de dicha técnica que se practica actualmente. Hace años, las condiciones naturales del río hacían posible la recarga natural de la Cubeta, pero la implantación de industrias en la zona y la regulación mediante la presa de La Baells desde el año 1976, han originado cambios en el comportamiento natural del acuífero, que se han manifestado básicamente en una disminución de la recarga natural como consecuencia de la impermeabilización de los suelos originada por la urbanización de los mismos. Se ha producido también un fenómeno de colmatación del lecho del río, al dificultarse la removilización de los materiales finos depositados en el fondo del cauce, produciendo también una disminución de la recarga al acuífero a través del lecho. Por último, la mayor demanda de suministro por parte del sector industrial ha influido negativamente en el almacenamiento del acuífero. Desde hace décadas, se lleva efectuando una recarga artificial en el lecho del río para tratar de paliar en lo posible estos efectos perjudiciales. La misma se efectúa mediante dos sistemas distintos: escarificando el lecho del río para aumentar su capacidad de infiltración, y mediante la utilización de balsas de recarga, a través de las cuales se infiltra el agua en el acuífero. El objetivo de este proyecto es analizar el efecto actual que la recarga mediante balsas está teniendo sobre el acuífero. Para ello, se utilizarán, entre otras técnicas, la elaboración de un modelo matemático con el que poder simular el efecto de la recarga. VIII ABSTRACT This project aims to show the possibilities of artificial recharge as groundwater resources management element in the alluvial aquifer of the Llobregat River, in the sector of the Basin of Sant Andreu de la Barca (Barcelona), which is a source of urban, agricultural and industrial supply of the Sant Andreu de la Barca area. Too, the effectiveness of the current implementation of this technique is reviewed. Years ago, natural river conditions made possible the natural recharge of the aquifer, but the establishment of industries in the area and the river regulation by the la Baells dam since 1976, have led to changes in the natural behavior of the aquifer, which basically are a decrease in natural recharge, because the soil has become impervious, due to land urbanization. There has been also a phenomenon of clogging of the river bed, caused by the difficult in the remobilization of fine materials deposited on the river bed, which has produced a decrease in the groundwater recharge through the river. Finally, the increase in water demand by the industrial sector, has diminished aquifer storage. For decades, artificial recharge in the river bed has been practiced to alleviate negative effects. It is done by two different systems: scratching the river bed to increase its infiltration capacity, and using recharge ponds, through which water infiltrates into the aquifer. The objective of this project is to analyze the current effect that recharge by ponds is having in the aquifer. Among other techniques, mathematical modeling has been used for simulating the effect of artificial recharge in the aquifer.

Cálculo de la recarga natural en grandes áreas, en función de la litología y la precipitación = Calculation of natural recharge in large areas, as a function of lithology and rainfall

En una región amplia como España se demuestra —mediante inferencias estadísticas sobre una muestra completa de 875 manantiales en los que se conoce su caudal medio y la litología de su área de alimentación y que han sido agrupados en regiones de distinta pluviometría— que la recarga media anual es una fracción fija de la precipitación media para cada litología. Se han establecido así unas tasas de recarga respecto de la precipitación para seis grupos litológicos de diferente permeabilidad: arenas, gravas y formaciones aluviales en general, 8.3%; conglomerados, 5.6%; areniscas, 7.3%; calizas y dolomías, 34.3%; margas, margocalizas, limos y arcillas, 3.3%; otras rocas, 1.3%. Teniendo en cuenta la representatividad de España, la cual tiene una gran variabilidad de litología, pluviometría, topografía, etcétera, estas tasas de recarga respecto de la precipita-ción son probablemente valores cuasi universales que pueden ser utilizados para estimar la recarga media o los recursos hídricos subterráneos medios de regiones amplias en cualquier parte del mundo, salvo en regiones especiales, como las que tienen permafrost, por ejemplo. En todo caso, estas tasas de recarga podrían ser retocadas para cada región según sus particulares características. Los datos de precipitación y litología son muy corrientes, por lo que el método puede ser ampliamente utilizado para completar balances hidráulicos.In a region as large as Spain, annual mean recharge is shown to be a fixed proportion of the mean rainfall for each lithology. This determination is based on statistical inferences from a complete sample of 875 springs for which mean flow and catchment areas are known and which have been grouped into distinct rainfall regions. Recharge rates have thus been established with respect to rainfall for six lithological groups with different permeability: sands, gravels and generally alluvial formations, 8.3%; conglomerates, 5.6%; sandstones, 7.3%; limestone and dolomite 34.3%; marls, marly limestones, silts and clays, 3.3%; and hard rocks, 1.3%. Considering the representativeness of Spain, which is large in size and has a highly varied lithology, topography and rainfall, these recharge rates for rainfall are probably quasi-universal values that can be used to estimate average recharge or average groundwater resources of large regions in any part of the world (except in special cases such as areas with permafrost, for example). For any case, these recharge rates can be adapted to each region according to its particular characteristics. Rainfall and lithology data are very common, and so the method can be widely used to calculate hydraulic balances.

Comportamiento hidrológico de cuencas de media montaña españolas: efectos de los procesos de acumulación / fusión de nieve en terrenos permeables sobre la infiltración y régimen de caudales Alto Tajo

La infiltración de agua en el suelo y la recarga profunda del agua subterránea contenida en los acuíferos es un proceso lento en relación con otros fenómenos hidrológicos. La redacción de esta tesis ha pretendido contribuir al estudio de la influencia que el almacenamiento de la precipitación sólida en forma de manto de nieve y su eventual fusión puedan tener sobre dicho proceso en áreas de media montaña (1.000 – 2.000 m.) en las que con gran frecuencia se sitúan las cabeceras de los ríos peninsulares. Para ello se ha partido del análisis de las diferentes variables intervinientes durante un determinado periodo temporal y sobre un espacio geográfico concreto, por lo que su metodología es de naturaleza empírica. La extensión del periodo (2002/03 a 2010/11) ha venido condicionada por la disponibilidad de los valores de algunas de sus principales variables, como han sido el equivalente en agua de la nieve acumulada y los caudales procedentes de su fusión. Éstos se han obtenido como resultado de la aplicación del modelo ASTER, desarrollado en el programa de Evaluación de los Recursos Hídricos procedentes de la Innivación (ERHIN), calibrado – entre otros- con datos de precipitaciones, temperatura y caudales provenientes a su vez del Sistema Automático de Información Hidrológica (SAIH). Ambos programas fueron implantados por la Administración en las diferentes Confederaciones Hidrográficas y en determinados Organismos de cuenca actuales, en cuyo desarrollo participó el autor de esta tesis. En cuanto a la zona de estudio se ha procedido a su elección considerando las posibles áreas de media montaña en las que la presencia de la nieve fuera hidrológicamente significativa y estuvieran constituidas litológicamente por afloramientos permeables que no impidieran la infiltración en el terreno y la formación de acuíferos de cierta relevancia. El interés se centró discrecionalmente en la cuenca del Tajo, tanto por el carácter estratégico de la misma -como suministradora en la actualidad de excedentes a otras cuencas deficitarias- como por el valor representativo de sus condiciones climáticas y orográficas en relación con otras cuencas hidrográficas peninsulares. Para ello se partió de las cabeceras de ríos identificadas por el programa ERHIN por su interés nivológico para la implantación del modelo ASTER y de las Masas de Agua Subterráneas MASb (antes Unidades Hidrogeológicas UUHH) definidas en los planes hidrológicos. La intersección en el territorio de ambos criterios condujo, finalmente, a la zona del Alto Tajo, en la que se cumplen ambos requisitos. El tramo quedó concretado en el comprendido entre las cabeceras de los ríos Tajo y Guadiela y la cola de los embalses de Entrepeñas y Buendía respectivamente, puntos de cierre para la calibración llevada a cabo en la modelización ASTER. Gran parte de éste discurre, en su parte alta, sobre rocas carbonatadas (calizas y dolomías del Jurásico y Cretácico), relacionados con las MASb de Tajuña-Montes Universales, Molina de Aragón y Sigüenza-Maranchón. Los valores diarios de las reservas de agua en forma de nieve, evapotranspiración y caudales procedentes de la fusión se han obtenido a partir de los resultados del mencionado modelo, procediéndose al cálculo de la infiltración por balance hídrico durante el periodo de estudio considerado, teniendo en cuenta los valores de precipitación, evapotranspiración y aportaciones de caudales. Esto ha requerido el estudio previo de las condiciones hidrogeológicas de la zona seleccionada con objeto de conocer las posibles interconexiones subterráneas que pudieran alterar los saldos entre las variables intervinientes anteriormente citadas. Para ello se ha llevado a cabo la recopilación y análisis de la información hidrogeológica correspondiente a la documentación de los planes hidrológicos del Tajo (Plan Hidrológico de la cuenca del Tajo RD 1664/1998 y el actual Plan Hidrológico de la parte española de la Demarcación Hidrográfica del Tajo RD 270/2014) y de los estudios previos realizados por el organismo de cuenca y el Instituto Geológico y Minero de España (lGME) fundamentalmente. En relación con la MASb Tajuña-Montes Universales -cuya extensión supera la zona seleccionada- dichos estudios consideran su estructura geológica y distribución litológica, con intercalaciones impermeables que actúan como barreras, dividiendo a éstas en Subunidades e identificando las zonas de drenaje de sus respectivos acuíferos. También se ha considerado la documentación y estudios previos del Plan Hidrológico Nacional sobre las Unidades Hidrogeológicas compartidas entre ámbitos geográficos de diferentes planes hidrológicos. Se concluye que las divisorias hidrográficas de las cabeceras son sensiblemente coincidentes o abarcan las Subunidades Montes Universales meridionales, Priego, Cifuentes, Zaorejas, u Montes Universales septentrionales, que drenan hacia el Tajo/Guadiela (bien directamente, bien a través de afluentes como el Gallo, Ablanquejo, Cabrillas, Cuervo…), MASb Molina de Aragón, que drena al Tajo a través del río Gallo y MASb Sigüenza—Maranchón, que drena su parte correspondiente hacia el Tajo a través del Ablanquejo. Se descartan – salvo la pequeña salvedad del manantial de Cifuentes- las conexiones hidrogeológicas con otras MASb o Subunidades por lo que las cabeceras del Tajo y del Guadiela pueden considerarse como un Sistema independiente donde las precipitaciones no evaporadas escurren superficialmente o se infiltran y descargan hacia los embalses de Entrepeñas y Buendía. La cuantificación diaria y acumulada de los balances hídricos ha permitido calcular la evolución aproximada de las reservas de agua subterránea desde la fecha inicial. Originalmente los balances se realizaron de forma separada en las cabeceras del Tajo y del Guadiela, cuyos valores acumulados manifestaron una tendencia creciente en la primera y decreciente en la segunda. Dicha situación se equilibra cuando el balance se practica conjuntamente en ambas, apreciándose en la variación del volumen de agua subterránea una evolución acorde hidrológicamente con los ciclos de verano/invierno y periodos de sequía, manteniéndose sus valores medios a largo/medio plazo, poniendo en evidencia la existencia de interconexiones subterráneas entre ambas cuencas. El balance conjunto, agregando la cabecera del Tajuña (que también comparte los materiales permeables de la MASb Tajuña-Montes Universales) no reveló la existencia de nuevas interrelaciones hidrogeológicas que influyeran en los balances hídricos realizados Tajo/Guadiela, confirmando las conclusiones de los estudios hidrogeológicos anteriormente analizados. Se ha procedido a confrontar y validar los resultados obtenidos de la evolución de las reservas de agua subterránea mediante los siguientes procedimientos alternativos: - Cálculo de los parámetros de desagüe de la curva de agotamiento correspondiente al volumen de agua subterránea drenante hacia el Tajo/Guadiela. Éste se ha realizado a partir de las aportaciones mensuales entrantes en los embalses de Entrepeñas y Buendía durante los meses de junio, julio, agosto y septiembre, cuyos valores responden al perfil típico de descargas de un acuífero. A partir de éstos se ha determinado el volumen drenante correspondiente al primero de junio de cada año de la serie histórica considerada. - Determinación del caudal base por el método Wallingford y deducción de los volúmenes drenantes. Estimación de las recarga anuales - Cuantificación de la recarga anual por el método Sanz, Menéndez Pidal de Navascués y Távara. Se obtuvieron valores de recarga muy aproximados entre los calculados por los dos últimos procedimientos citados. Respecto a las reservas de agua subterránea almacenadas siguen una evolución semejante en todos los casos, lo que ha permitido considerar válidos los resultados conseguidos mediante balance hídrico. Confirmada su solidez, se han buscado correlaciones simples entre el volumen de las reservas subterráneas (como indicador estimativo del efecto de la infiltración) y los volúmenes procedentes de la fusión. La conclusión es que estos últimos no tienen un efecto determinante a escala anual sobre la infiltración,recarga y variación de los volúmenes de agua subterránea, frente al peso de otras variables (precipitación y evapotranspiración). No obstante se ha encontrado una buena correlación múltiple entre la recarga estimada y la precipitación eficaz (precipitación menos evapotranspiración) y fusión, que ha permitido cuantificar la contribución de esta última. Posteriormente se ha recurrido a la selección de los episodios más intensos de acumulación /fusión en las cabeceras del Tajo y Guadiela. Y se procedió a la comparación entre los resultados obtenidos por aplicación del modelo de simulación en los mismos periodos (normalmente de varios días de duración) con datos reales y con datos ficticios de temperatura que anularan o disminuyeran la presencia de nieve, apreciándose una gran sensibilidad del efecto de la temperatura sobre la evapotranspiración y estableciéndose nuevamente correlaciones lineales entre los volúmenes de fusión y el incremento de reservas subterráneas. Las mismas confirman el efecto “favorecedor” de la acumulación de agua en forma de nieve y su posterior licuación, sobre sobre la infiltración de agua en el suelo y almacenamiento subterráneo. Finalmente se establecieron varios escenarios climáticos (+1ºC; +3ºC; +1ºC y – 10% precipitación; y 3ºC – 10% precipitación) compatibles con las previsiones del IPCC para mediados y finales del presente siglo, determinándose mediante simulación ASTER los correspondientes valores de fusión. La correlación establecida a escala anual ha permitido evaluar el efecto de la disminución del volumen de fusión - en los diferentes escenarios – sobre la recarga, pronosticando un descenso de los caudales de estiaje y la desaparición del “efecto nieve” sobre la infiltración y recarga con un aumento de 3ºC de temperatura. Teniendo en cuenta las condiciones de representatividad de la zona elegida, resulta verosímil la extensión de las anteriores conclusiones a otras cabeceras fluviales enclavadas en áreas de media montaña situadas entre 1000 a 2000m y sus efectos aguas abajo.Water infiltration into the soil and groundwater recharge deep water in aquifers is slow relative to other hydrological phenomena. The wording of this thesis aims to contribute to the study of the influence that the storage of solid precipitation as snow cover and its eventual melting may have on this process in mid-mountain areas (1000 - 2,000 m) where very often the headwaters of the peninsular rivers are located. For this party analysis of the different variables involved has over a given time period and a particular geographical area, so that their methodology is empirical in nature. The extension of the period (2002/03 to 2010/11) has been conditioned by the availability of the values of some of its key variables, as were the water equivalent of the snow and flows from melting. These have been obtained as a result of the application of ASTER model, developed in the program Evaluation of Water Resources from the Innivation (ERHIN), calibrated - among others data of rainfall, temperature and flow from turn System Automatic Hydrological Information (SAIH). Both programs were implemented by the Administration in the different Water Boards and to undertakings for current basin, in which the author participated development of this thesis. As for the study area has proceeded at its option considering the possible areas of midmountain in the presence of snow outside hydrological meaningful and they were lithology consisting of permeable outcrops that did not prevent infiltration into the ground and forming aquifers of some significance. We were interested discretion in the Tagus basin, therefore the strategic nature of it, as currently supplying surplus to other basins deficit- as the representative value of its climate and terrain conditions in relation to other peninsular river basins . To do this we started from the headwaters identified by the ERHIN program for its implementation snow interest to the ASTER model and Ground Water Bodies MASb (before UUHH Hydrogeological Units) defined in hydrological plans. The intersection in the territory of both criteria led eventually to the Alto Tajo, in which both requirements are met. The section was finalized in the period between the headwaters of the Tagus and Guadiela rivers and reservoirs end Entrepeñas and Buendia respectively checking points for calibration performed in ASTER modeling. Much of it runs on carbonate rocks (limestones and dolomites of Jurassic and Cretaceous) related MASb of Tajuña -Montes Universal, Molina de Aragón and Sigüenza-Maranchón. The daily values of water reserves in the form of snow, evapotranspiration and flow from melting were obtained from the results of this model, proceeding to the calculation of infiltration water balance during the study period considered, taking into account values of precipitation, evapotranspiration and input flow. This has required the prior examination of the hydrogeological conditions of your required in order to know the possible underground interconnections that could alter the balance between the intervening variables aforementioned area. For this we have carried out the collection and analysis of hydrogeological information relevant documentation Tagus river management plans (Hydrological Plan Tajo Basin RD 1664/1998 and the current Hydrological Plan of the Spanish part of the River Basin Tagus RD 270/2014) and previous studies by the basin organization and the Geological Survey of Spain (IGME) mainly. Regarding the MASb Tajuña- Montes Universal - whose length exceeds the area selected - these studies consider its geological structure and lithology distribution with waterproof collations that act as barriers, dividing it into subunits and identifying areas draining their respective aquifers. It has also considered the documentation and previous studies of the National Hydrological Plan on shared among different geographical areas management plans Hydrogeological Units. We conclude that river dividing the headers are substantially coincident or covering Subunits southern Universal Montes, Priego Cifuentes, Zaorejas and northern Universal Mounts, which drain into the Tagus / Guadiela (either directly or through tributaries such as Gallo, Ablanquejo , whitecaps , Raven ...), MASb Molina de Aragón which drains through the Tajo del Gallo and MASb Sigüenza- Maranchón river that drains into the Tagus using the Ablanquejo . Discarded - except the small exception of spring Cifuentes -hydrogeological connections with other MASb or Subunits so the headwaters of the Tagus and Guadiela be considered as a separate system, where rainfall not evaporated runs on surface or infiltrates and eventually discharged into reservoirs Entrepeñas and Buendia. The daily and cumulative quantification of water balances allowed us to compute the approximate evolution of groundwater reserves from its initial date. Initially balances were performed separately in the headwaters of the Tagus and Guadiela, whose cumulative values showed an increasing trend in the first and decreasing in the second. This situation is balanced when the balance is practiced together in both , appreciating the change in volume of groundwater hydrological evolution commensurate with the cycles of summer / winter and drought periods , keeping their average long / medium term values and putting in shows the existence of underground interconnections between the two basins. The overall balance, adding header Tajuña (which also shares the permeable materials MASb Tajuña -Montes Universal ) did not reveal the existence of new hydrogeological interrelationships that influenced water balances made Tajo / Guadiela, confirming the findings of the hydrogeological studies previously analyzed. We proceeded to confront and validate the results of the evolution of groundwater reserves by the following alternative procedures: - Calculate the parameters drain depletion curve corresponding to the volume of groundwater draining into the Tajo / Guadiela. This has been made from monthly inflows in the reservoirs of Entrepeñas and Buendia during the months of June, July, August and September, whose values match the typical profile of an aquifer discharges. From these has been determined for the first of June each year of the time series considered drainage volume - Determination of base flow by Wallingford method and deduction of drainage volumes. Estimate of annual recharge - Quantification of the annual recharge by the method Sanz Menéndez Pidal of Navascués and Távara. Very approximate values recharge between calculated for the last two mentioned methods were obtained. Concerning groundwater reserves stored follow a similar pattern in all cases, allowing consider valid the results achieved through water balance. Confirmed its robustness, simple correlations were sought between the volume of groundwater reserves (as estimated indicator of the effect of infiltration) and volumes from the melting. The conclusion is that the latter do not have a decisive effect on the annual scale infiltration, recharge and variation in volumes of groundwater, against the weight of other variables (precipitation and evapotranspiration). However found a good multiple correlation between the estimated recharge and effective precipitation (precipitation minus evapotranspiration) and fusion, which allowed quantify the contribution of the latter. Subsequently it has resorted to the selection of the most intense episodes of accumulation / melting in the headwaters of the Tagus and Guadiela. And we proceeded to the comparison between the results obtained by application of the simulation model in the same periods (usually several days) with real data and fictitious temperature data to annul or decrease the presence of snow, appreciating a great sensitivity of the effect of temperature on evapotranspiration and establishing linear correlations between the volumes of melting and increased groundwater reserves again. They confirm the “flattering " effect of water accumulation as snow and subsequent liquefaction of the infiltration of water into the soil and underground storage. Finally various climate scenarios (+1ºC; +3ºC; +1ºC y – 10% precipitation; y 3ºC – 10% precipitation) were established consistent with IPCC projections for mid - to late - century, determined through simulation ASTER corresponding values of melting. The correlation established on an annual scale has allowed to evaluate the effect of decreasing the volume of melt - in different scenarios - on recharge, predicting a decline in low flows and the disappearance of "snow effect" on infiltration and recharge with an increase of 3°C temperature. Given the conditions of representativeness of the chosen area, plausible extension of the above findings to other landlocked headwaters in mid-mountain areas located between 1000 to 2000m and its downstream effects.

Desarrollo de métodos numérico-interpretativos para la realización de ensayos de bombeo

Los ensayos de bombeo son, sin lugar a dudas, una de las pruebas más fiables y de mayor interés que se hacen en el medio físico. No son pruebas estrictamente puntuales, dado que el bombeo atrae flujo desde distancias lejanas al pozo, la prueba tiene una excelente representatividad espacial. Los métodos de interpretación mediante ensayos de bombeo se empezaron a plantear en la primera mitad del pasado siglo. Con los ensayos de bombeo se puede calcular la transmisividad y coeficiente de almacenamiento de las formaciones acuíferas y suministran información sobre el tipo de acuífero, la calidad constructiva del pozo de extracción, la existencia de barreras impermeable o bordes de recarga próximos, e incluso en algunas circunstancias permiten el cálculo del área de embalse subterráneo. Desde mediados del siglo 20 existe una eficaz y abundante gama de métodos analítico-interpretativos de ensayos de bombeo, tanto en régimen permanente como transitorio. Estos métodos son ampliamente conocidos y están muy experimentados a lo largo de muchos países, sin embargo, hoy día, podrían utilizarse modelos de flujo para la interpretación, logrando la misma fiabilidad e incluso mejores posibilidades de análisis. Muchos ensayos que no pueden interpretarse porque las configuraciones del medio son demasiado complejas y no están disponibles, o no es posible, el desarrollo de métodos analíticos, tienen buena adaptación y en ocasiones muy fácil solución haciendo uso de los métodos numéricos de simulación del flujo. En esta tesis se ha buscado una vía de interpretar ensayos de bombeo haciendo uso de modelos de simulación del flujo. Se utiliza el modelo universal MODFLOW del United States Geological Survey, en el cual se configura una celda de simulación y mallado particularmente adecuados para el problema a tratar, se valida con los métodos analíticos existentes. Con la célula convenientemente validada se simulan otros casos en los que no existen métodos analíticos desarrollados dada la complejidad del medio físico a tratar y se sacan las oportunas conclusiones. Por último se desarrolla un modelo específico y la correspondiente aplicación de uso general para la interpretación numérica de ensayos de bombeo tanto con las configuraciones normales como con configuraciones complejas del medio físico. ABSTRACT Pumping tests are, without doubt, one of the most reliable and most interesting tests done in the physical environment. They are not strictly anecdotal evidence, since pumping flow attracts from far distances to the well, the test has excellent spatial representation. Methods of interpretation by pumping tests began to arise in the first half of last century. With pumping tests, can be calculated transmissivity and storage coefficient of the aquifer formations, and provide information on the type of aquifer, the construction quality of the well, the existence of waterproof barriers or borders next recharge, and even in some circumstances allow calculating the area of underground reservoir. Since the mid-20th century there is effective and abundant range of analytical interpretative pumping tests, both in steady state and transient methods. These methods are very widely known and experienced over many countries, however, nowadays, may flow models used for interpretation, obtaining equally reliable or even better possibilities for analysis. Many trials cannot be interpreted as environmental settings are too complex and are not available, or not possible, the development of analytical methods, have good adaptation and sometimes very easily solved using numerical flow simulation methods. This thesis has sought a way to interpret pumping tests using flow simulation models. MODFLOW universal model of United States Geological Survey, in which a simulation cell and meshing particularly suitable for the problem to be treated, is validated with existing analytical methods used is set. With suitably validated cell other cases where there are no analytical methods developed given the complexity of the physical environment to try and draw appropriate conclusions are simulated. Finally, a specific model and the corresponding application commonly used for numerical interpretation of pumping tests both with normal settings as complex configurations of the physical environment is developed.

DInSAR monitoring of aquifer compaction due to water withdrawal: Vegas Baja and Media of the Segura river (SE, Spain) case study

This study was partially financed by the Spanish Ministry of Education and Science and EU FEDER under project TEC2005-06863, by the Valencia Regional Government under projects GV006/179 and ACOMP07/087, and by the University of Alicante under projects VIGROB2004102, VIGROB-053, and VIGROB-114.

A quasi-elastic aquifer deformational behavior: Madrid aquifer case study

The purpose of this paper is to analyze the quasi-elastic deformational behavior that has been induced by groundwater withdrawal of the Tertiary detrital aquifer of Madrid (Spain). The spatial and temporal evolution of ground surface displacement was estimated by processing two datasets of radar satellite images (SAR) using Persistent Scatterer Interferometry (PSI). The first SAR dataset was acquired between April 1992 and November 2000 by ERS-1 and ERS-2 satellites, and the second one by the ENVISAT satellite between August 2002 and September 2010. The spatial distribution of PSI measurements reveals that the magnitude of the displacement increases gradually towards the center of the well field area, where approximately 80 mm of maximum cumulated displacement is registered. The correlation analysis made between displacement and piezometric time series provides a correlation coefficient greater than 85% for all the wells. The elastic and inelastic components of measured displacements were separated, observing that the elastic component is, on average, more than 4 times the inelastic component for the studied period. Moreover, the hysteresis loops on the stress–strain plots indicate that the response is in the elastic range. These results demonstrate the quasi-elastic behavior of the aquifer. During the aquifer recovery phase ground surface uplift almost recovers from the subsidence experienced during the preceding extraction phase. Taking into account this unique aquifer system, a one dimensional elastic model was calibrated in the period 1997–2000. Subsequently, the model was used to predict the ground surface movements during the period 1992–2010. Modeled displacements were validated with PSI displacement measurements, exhibiting an error of 13% on average, related with the inelastic component of deformation occurring as a long-term trend in low permeability fine-grained units. This result further demonstrates the quasi-elastic deformational behavior of this unique aquifer system.

Brackish groundwater desalination by reverse osmosis in southeastern Spain. Presence of emerging contaminants and potential impacts on soil-aquifer media

Desalinated brackish groundwater is becoming a new source of water supply to comply with growing water demands, especially in (semi-) arid countries. Recent publications show that some chemical compounds may persist in an unaltered form after the desalination processes and that there is an associated risk of mixing waters with different salinity for irrigation. At the university of Alicante campus (Spain), a mix of desalinated brackish groundwater and water from the existing aquifer is currently applied for landscape irrigation. The presence of 209 emerging compounds, surfactants, priority substances according to the 2008/105/EC Directive, 11 heavy metals and microbiological organisms in blended water and aquifer samples was investigated. Thirty-five compounds were detected (pesticides, pharmaceuticals and surfactants) among them two priority substances α-endosulfan and Ni were found above the permitted maximum concentration. Blended water used for landscape irrigation during the summer period is supersaturated with respect to carbonates, which may ultimately lead to mineral precipitation in the soil-aquifer media and changes in hydraulic parameters.

Efectos del cambio climático y usos del suelo sobre los recursos hídricos de la cuenca del río Tordera (Barcelona, España)

Se presentan los efectos del cambio global en la cuenca del río Tordera (España) para el periodo 2000-2050, escenarios climáticos A2 (medio-alto) definidos por el Panel Intergubernamental del Cambio Climático (IPCC, 200) y escenarios socioeconómicos (cambios previstos en la cuenca) denominados estable y tendencial. Los efectos sobre los recursos hídricos se han analizado de forma conjunta superficial-subterránea mediante una metodológica de tipo acoplado. Para establecer los impactos futuros sobre los recursos hídricos se ha seleccionado el Modelo de Circulación Global ECHAM5 (Max Planck Institute). Los resultados obtenidos indican una disminución de la precipitación del 11.3% y un aumento de la temperatura de 1ºC, respecto a los valores históricos de la zona. De acuerdo a la proyección futura (2050) sobre cambios en los recursos hídricos, la escorrentía superficial obtenida mediante simulación con el código HEC-HMS 3.4 experimenta una reducción del 31.8% respecto al valor histórico y la recarga natural, estimada mediante VISUAL-Balan, se reduce en un 11.7%. El balance en el acuífero deltaico simulado mediante MODFLOW 2009.1 Pro muestra igualmente una disminución de los parámetros del balance. Los cambios del uso del suelo previstos de acuerdo a la legislación vigente (escenarios socioeconómicos) no conducen a la generación de un impacto apreciable en los recursos hídricos; según los escenarios definidos la variación de precipitación y temperatura son los parámetros fundamentales del cambio previsto.