Effects on the shallow aquifers by CO2 leakages in a tertiary basin (province of Murcia, Spain)

Geological storage of CO2 is nowadays internationally considered as the most effective method for greenhouse gas emission mitigation, in order to minimize its effects on the global climatology. One of the main options is to store CO2 in deep saline aquifers at more than 800m depth, because it reaches its supercritical state. Study of the CO2 natural accumulations as natural analogues of an artificial CO2 storage is very useful in order to understand the CO2 long term behaviour and thus to predict its possible impact on the surficial environment and life. Therefore the main objective of this work is to detect the affection of the CO2 leakages from a deep saline aquifer on the shallow aquifers, all of them located in the Gañuelas-Mazarrón Tertiary basin (Province of Murcia, Spain). This CO2 storage and leakage natural system can be analogous to an artificial CO2 storage with leakage phenomena. In order to reach these objectives, groundwaters from different aquifers in the site have been sampled and analysed for major elements, free and dissolved gases and stable isotopes, particularly ∂ 13 C and 3 He/ 4 He. The results obtained allow to conclude that this natural system is an interesting example of natural analogue for an artificial CO2 storage affected by leakage processes because the shallow fresh aquifers in the site are polluted by CO2 from the deep saline aquifer as a consequence of an intensive over-exploitation of these freshwater aquifers

Interaction between karstic aquifers and allogenic rivers: the aquifer of the national park of the ephemeral river Lobos Canyon (Spain)

A methodology is described for understanding the interaction of karstic aquifers with allogenic rivers, where little information is available. This methodology includes conventional hydrogen- ology methods tracer tests and measurements of flow into, out of and circulating within the karstic system. The method is designed to un- derstand the hydrogeological behaviour of a river in sufficient detail, given a short study pe- riod. The methodology is applied to a karstic system in Spain, obtaining useful, quantitative results for a hydrological year, such as an esti- mate of the water balance, differentiation be- tween autogenic and allogenic natural recharge, relationship and connection between the river and the aquifer, and measurements of infiltration capacity in watercourses under different hydro- logical situations. The paper deals with a useful example that could be applied to other rivers and aquifers where few data are available. It can be applied to aquifers under a natural regime and Mediterranean climate.

Optimización del Abonado Nitrogenado en el Melón (Cucumis Melo L.) Tipo Piel de Sapo

Spain is the fifth-largest producer of melon (Cucumis melo L.) and the second exporter in the world. To a national level, Castilla-La Mancha emphasize and, specifically, Ciudad Real, where is cultivated 27% of national area dedicated to this crop and 30% of melon national production. Melon crop is cultivating majority in Ciudad Real and it is mainly located in the Alto Guadiana, where the major aquifers of the region are located, the aquifer 23 or Mancha Occidental and the aquifer 24 or Campo de Montiel, both declared overexploited and vulnerable zones to nitrate pollution from agricultural sources. The problem is exacerbated because in this area, groundwater is the basic resource of supply to populations, and even often the only one. Given the importance of melon in the area, recent research has focused on the irrigation of melon crop. Unfortunately, scant information has been forthcoming on the effect of N fertilizer on melon piel de sapo crop, so it is very important to tackle in a serious study that lead to know the N requirements on the melon crop melon by reducing the risks of contamination by nitrate leaching without affecting productivity and crop quality. In fact, the recommended dose is often subjective and practice is a N overdose. In this situation, the taking of urgent measures to optimize the use of N fertilization is required. To do it, the effect of N in a melon crop, fertirrigated and on plastic mulch, was studied. The treatments consisted in different rates of N supply, considering N fertilizer and N content in irrigation water, so the treatment applied were: 30 (N30), 85 (N85), 112 (N112) and 139 (N139) Kg N ha-1 in 2005; 93 (N93), 243 (N243) and 393 (N393) kg ha-1 in 2006; and 11 (N11), 61 (N61), 95 (N95) and 148 (N148) kg ha-1 in 2007. A randomized complete-block design was used and each treatment was replicated four times. The results showed a significant effect of N on dry biomass and two patterns of growth were observed. On the one hand, a gradual increase in vegetative biomass of the plant, leaves and stem, with increasing N, and on the other hand, an increase of fruit biomass also with increasing N up to a maximum of biomass corresponding to the optimal dose determined in 90 kg ha-1 of N applied, corresponding to 160 kg ha-1 of N available for melon crop, since this optimum dose, the fruit biomass suffers a decline. A significant effect was observed in concentration and N uptake in leaf, steam, fruit and whole plant, increasing in all of them with increasing of N doses. Fast N uptake occurred from 30-35 to 70-80 days after transplanting, coinciding with the fruit development. The N had a clear influence on the melon yield, its components, skin thickness and flesh ratio. The melon yield increased, as the mean fruit weight and number of fruits per m2 with increasing N until achieve an above 95% of the maximum yield when the N applied is 90 kg ha-1 or 160 kg ha-1 of N available. When N exceeds the optimal amount, there is a decline in yield, reducing the mean fruit weight and number of fruits per square meter, and was also observed a decrease in fruit quality by increasing the skin thickness and decrease the flesh ratio, which means an increase in fruit hollowed with excessive N doses. There was a trend for all indexes of N use efficiency (NUE) to decline with increasing N rate. We observed two different behaviours in the calculation result of the NUE; on the one hand, all the efficiency indexes calculated with N applied and N available had an exponential trend, and on the other hand, all the efficiency indexes calculated with N uptake has a linear trend. The linear regression cuts the exponential curve, delimiting a range within which lies the optimum quantity of N. The N leaching as nitrates increased exponentially with the amount of N. The increase of N doses was affected on the N mineralization. There was a negative exponential effect of N available on the mineralization of this element that occurs in the soil during the growing season, calculated from the balances of this element. The study of N leaching for each N rate used, allowed to us to establish several environmental indices related to environmental risk that causes the use of such doses, a simple way for them to be included in the code of Best Management Practices.

Diffuse Soil Co2 Flux To Assess The Reliability Of Co2 Storage In The Mazarrón-Gañuelas Tertiary Basin (Spain)

Geological storage of CO2 is nowadays internationally considered as the most effective method for greenhouse gas emission mitigation, in order to minimize its effects on the global climatology. One of the main options is to store the CO2 in deep saline aquifers at more than 800 m depth, because it achieves its supercritical state. Among the most important aspects concerning the performance assessment of a deep CO2 geological repository is the evaluation of the CO2 leakage rate from the chosen storage geological formation. Therefore, it is absolutely necessary to increase the knowledge on the interaction among CO2, storage and sealing formations, as well as on the flow paths for CO2 and the physico-mechanical resistance of the sealing formation. Furthermore, the quantification of the CO2 leakage rate is essential to evaluate its effects on the environment. One way to achieve this objective is to study of CO2 leakage on natural analogue systems, because they can provide useful information about the natural performance of the CO2, which can be applied to an artificial CO2 geological storage. This work is focused on the retention capacity of the cap-rock by measuring the diffuse soil CO2 flux in a site selected based on: i) the presence of a natural and deep CO2 accumulation; ii) its structural geological characteristics; and iii) the nature of the cap-rocks. This site is located in the so-called Mazarrón-Gañuelas Tertiary Basin, in the Guadalentin Valley, province of Murcia (Spain) Therefore the main objective of this investigation has been to detect the possible leakages of CO2 from a deep saline aquifer to the surface in order to understand the capability of this area as a natural analogue for Carbon Capture and Sequestration (CCS). The results obtained allow to conclude that the geological sealing formation of the basin seems to be appropriate to avoid CO2 leakages from the storage formation.

Sampling strategies using the "accumulation chamber" for monitoring geological

Fundación Ciudad de la Energía (CIUDEN) is carrying out a project of geological storage of CO2, where CO2 injection tests are planned in saline aquifers at a depth of 1500 m for scientific objectives and project demonstration. Before any CO2 is stored, it is necessary to determine the baseline flux of CO2 in order to detect potential leakage during injection and post-injection monitoring. In November 2009 diffuse flux measurements of CO2 using an accumulationchamber were made in the area selected by CIUDEN for geological storage, located in Hontomin province of Burgos (Spain). This paper presents the tests carried out in order to establish the optimum sampling methodology and the geostatistical analyses performed to determine the range, with which future field campaigns will be planned.

A method to assess annual average renewable groundwater reserves for large regions in Spain

This paper proposes a method for assessing the groundwater renewable reserves of large regions for an average year, based on the integration of the recession curves for their basins springs or the natural base flow of their rivers. In this method, the hydrodynamic volume (or renewable reserves), were estimated from the baseflow equation. It was assumed that the flow was the same as the natural recharge, and that the recession coefficients were derived by the hydrogeological parameters and geometrical characteristics of aquifers, and adjusted to fit the recession curves at gauging stations. The method was applied to all the aquifers of Spain, which have a total groundwater renewable reserve of 86,895 hm3 four times the mean annual recharge. However, the distribution of these reserves is very variable; 18.6% of the country aquifers contain 94.7% of the entire reserve.

Hydrogeochemistry of surface and spring waters in the surroundings of the CO2 injection site at Hontomín–Huermeces (Burgos, Spain)

In this paper the very first geochemical and isotopic data related to surface and spring waters and dissolved gases in the area of Hontomín–Huermeces (Burgos, Spain) are presented and discussed. Hontomín–Huermeces has been selected as a pilot site for the injection of pure (>99%) CO2. Injection and monitoring wells are planned to be drilled close to 6 oil wells completed in the 1980s for which detailed stratigraphical logs are available, indicating the presence of a confined saline aquifer at the depth of about 1500 m into which less than 100,000 tons of iquid CO2 will be injected, possibly starting in 2013. The chemical and features of the spring waters suggest that they are related to a shallow hydrogeological system as the concentration of the Total Dissolved Solids approaches 800 mg/L with a Ca2+(Mg2+)-HCO3− composition, similar to that of the surface waters. This is also supported by the oxygen and hydrogen isotopic ratios that have values lying between those of the Global and the Mediterranean Meteoric Water Lines. Some spring waters close to the oil wells are haracterized by relatively high concentrations of NO3− (up to 123 mg/L), unequivocally suggesting an anthropogenic source that adds to the main water–rock interaction processes. The latter can be referred to Ca-Mg-carbonate and, at a minor extent, Al-silicate dissolution, being the outcropping sedimentary rocks characterized by Palaeozoic to Quaternary rocks. Anomalous concentrations of Cl−, SO42−, As, B and Ba were measured in two springs discharging a few hundred meters from the oil wells and in the Rio Ubierna. These contents are significantly higher than those of the whole set of the studied waters and are possibly indicative of mixing processes, although at very low extent, between deep and shallow aquifers. No evidence of deep-seated gases interacting with the Hontomín–Huermeces waters was recognized in the chemistry of the disolved gases. This is likely due to the fact that they are mainly characterized by an atmospheric source as highlighted by the high contents of N2, O2 and Ar and by N2/Ar ratios that approach that of ASW (Air Saturated Water) and possibly masking any contribution related to a deep source. Nevertheless, significant concentrations (up to 63% by vol.) of isotopically negative CO2 (<−17.7‰ V-PDB) were found in some water samples, likely related to a biogenic source. The geochemical and isotopic data of this work are of particular importance when a monitoring program will be established to verify whether CO2 leakages, induced by the injection of this greenhouse gas, may be affecting the quality of the waters in the shallow hydrological circuits at Hontomín–Huermeces. In this respect, carbonate chemistry, the isotopic carbon of dissolved CO2 and TDIC (Total Dissolved Inorganic Carbon) and selected trace elements can be considered as useful parameters to trace the migration of the injected CO2 into near-surface environments.

CO2 soil flux baseline at the technological development plant for CO2 injection at Hontomin (Burgos, Spain)

From the end of 2013 and during the following two years, 20 kt of CO2sc are planned to be injected in a saline reservoir (1500 m depth) at the Hontomín site (NE Spain). The target aquifers are Lower Jurassic limestone formations which are sealed by Lower Cretaceous clay units at the Hontomín site (NE Spain). The injection of CO2 is part of the activities committed in the Technology Development phase of the EC-funded OXYCFB300 project (European Energy Program for Recovery – EEPR, http://www.compostillaproject.eu), which include CO2 injection strategies, risk assessment, and testing and validating monitoring methodologies and techniques. Among the monitoring works, the project is intended to prove that present-day technology is able to monitor the evolution of injected CO2 in the reservoir and to detect potential leakage. One of the techniques is the measurement of CO2 flux at the soil–atmosphere interface, which includes campaigns before, during and after the injection operations. In this work soil CO2 flux measurements in the vicinity of oil borehole, drilled in the eighties and named H-1 to H-4, and injection and monitoring wells were performed using an accumulation chamber equipped with an IR sensor. Seven surveys were carried out from November 2009 to summer 2011. More than 4000 measurements were used to determine the baseline flux of CO2 and its seasonal variations. The measured values were low (from 5 to 13 g m−2 day−1) and few outliers were identified, mainly located close to the H-2 oil well. Nevertheless, these values cannot be associated to a deep source of CO2, being more likely related to biological processes, i.e. soil respiration. No anomalies were recognized close to the deep fault system (Ubierna Fault) detected by geophysical investigations. There, the CO2 flux is indeed as low as other measurement stations. CO2 fluxes appear to be controlled by the biological activity since the lowest values were recorded during autumn-winter seasons and they tend to increase in warm periods. Two reference CO2 flux values (UCL50 of 5 g m−2 d−1 for non-ploughed areas in autumn–winter seasons and 3.5 and 12 g m−2 d−1 for in ploughed and non-ploughed areas, respectively, in spring–summer time, and UCL99 of 26 g m−2 d−1 for autumn–winter in not-ploughed areas and 34 and 42 g m−2 d−1 for spring–summer in ploughed and not-ploughed areas, respectively) were calculated. Fluxes higher than these reference values could be indicative of possible leakage during the operational and post-closure stages of the storage project.

Water chemistry at Hontomín-Huermeces (Burgos, Spain): insights for a pre-, intra- and post-CO2 injection geochemical monitoring.

In this study, the very first geochemical and isotopic data related to surface and spring waters and dissolved gases in the area of Hontomín-Huermeces (Burgos, Spain) are presented and discussed. Hontomín-Huermeces was selected as a pilot site for the injection of pure (>99 %) CO2. Injection and monitoring wells are planned to be drilled close to 6 oil wells completed in the 1980’s. Stratigraphical logs indicate the presence of a confined saline aquifer at the depth of about 1,500 m into which less than 100,000 tons of liquid CO2 will be injected, possibly starting in 2013. The chemical and isotopic features of the spring waters suggest the occurrence of a shallow aquifer having a Ca2+(Mg2+)-HCO3- composition, relatively low salinity (Total Dissolved Solids _800 mg/L) and a meteoric isotopic signature. Some spring waters close to the oil wells are characterized by relatively high concentrations of NO3- (up to 123 mg/L), unequivocally indicating anthropogenic contamination that adds to the main water-rock interaction processes. The latter can be referred to Ca-Mg-carbonate and, at a minor extent, Al-silicate dissolution, being the outcropping sedimentary rocks characterized by Palaeozoic to Quaternary rocks. Anomalous concentrations of Cl-, SO42-, As, B and Ba were measured in two springs discharging a few hundreds meters from the oil wells and in the Rio Ubierna, possibly indicative of mixing processes, although at very low extent, between deep and shallow aquifers. Gases dissolved in spring waters show relatively high concentrations of atmospheric species, such as N2, O2 and Ar, and isotopically negative CO2 (<-17.7 h V-PDB), likely related to a biogenic source, possibly masking any contribution related to a deep source. The geochemical and isotopic data of this study are of particular importance when a monitoring program will be established to verify whether CO2 leakages, induced by the injection of this greenhouse gas, may affect the quality of the waters of the shallow Hontomín-Huermeces hydrological circuit. In this respect, carbonate chemistry, the isotopic carbon of dissolved CO2 and TDIC (Total Dissolved Inorganic Carbon) and selected trace elements can be considered as useful parameters to trace the migration of the injected CO2 into near-surface environments.

Identificación de los riesgos ambientales y sanitarios de la producción de gas mediante fracturación hidráulica y bases para una propuesta metodológica de vulnerabilidad de las aguas subterráneas

Hasta ahora, la gran mayoría de los recursos explotados de gas natural procedían de acumulaciones convencionales de gas aislado y de gas asociado y disuelto en el petróleo. Sin embargo, el gas natural se encuentra también en yacimientos que, debido a su baja porosidad y permeabilidad, tienen unas características que hacen que hasta muy recientemente no hayan sido económicamente rentables y que sólo puedan ser explotados mediante técnicas no convencionales, dando lugar al denominado gas no convencional. Las técnicas utilizadas para su extracción son la fracturación hidráulica o “fracking” y la perforación horizontal. Entre los diversos tipos de gas no convencional, es de prever que el gas de pizarra sea el que sufra mayor desarrollo a medio plazo en nuestro país, por lo que se están generando grandes debates, debido al riesgo de contaminación de las aguas superficiales y subterráneas del entorno, provocados por la elevada cantidad de agua utilizada, los aditivos empleados, los fluidos de retorno, por la alteración del medio físico, así como por la dificultad de monitorización de estos procesos. En este proyecto se identifican los riesgos ambientales y sanitarios asociados a la extracción de gas no convencional. El trabajo se basa principalmente en experiencias ocurridas en países donde el fracking se ha convertido en una práctica habitual. Se trata además de establecer las bases necesarias para la estimación de la vulnerabilidad de los acuíferos frente a la contaminación inducida por la fracturación hidráulica. Abstract Until now, most of the natural gas resources exploited were from isolated conventional gas accumulations and associated and dissolved gas in oil. However, the natural gas is also in reservoirs that, due to their low porosity and permeability, have characteristics that make until recently not been economically profitable and can be exploited only by unconventional techniques, leading to the so called unconventional gas. The techniques used for extraction are hydraulic fracturing or "fracking" and horizontal drilling. Among the various types of unconventional gas, it is expected that shale gas is the suffering greater medium-term development in our country, so it is generating much debate, due to the risks of contamination in surface waters and subterranean environment, caused by the high amount of water used, the additives used, the return fluid, by altering the physical environment, and the difficulty of monitoring these processes. In this project identifies the environmental and health risks associated with unconventional gas extraction. The work is mainly based on experiences that occurred in countries where fracking has become a common practice. This is for establish the necessary basis for estimating the vulnerability of aquifers from contamination induced by hydraulic fracturing.

Estrategias de monitorización de CO2 y otros gases en los estudios de análogos naturales

La Universidad Politécnica de Madrid (UPM) y la Università degli Studi di Firenze (UniFi), bajo la coordinación técnica de AMPHOS21, participan desde 2009 en el proyecto de investigación “Estrategias de Monitorización de CO2 y otros gases en el estudio de Análogos Naturales”, financiado por la Fundación Ciudad de la Energía (CIUDEN) en el marco del Proyecto Compostilla OXYCFB300 (http://www.compostillaproject.eu), del Programa “European Energy Program for Recovery - EEPR”. El objetivo principal del proyecto fue el desarrollo y puesta a punto de metodologías de monitorización superficiales para su aplicación en el seguimiento y control de los emplazamientos donde se realice el almacenamiento geológico de CO2, analizando técnicas que permitan detectar y cuantificar las posibles fugas de CO2 a la atmósfera. Los trabajos se realizaron tanto en análogos naturales (españoles e italianos) como en la Planta de Desarrollo Tecnológico de Almacenamiento de CO2 de Hontomín. Las técnicas analizadas se centran en la medición de gases y aguas superficiales (de escorrentía y manantiales). En cuanto a la medición de gases se analizó el flujo de CO2 que emana desde el suelo a la atmósfera y la aplicabilidad de trazadores naturales (como el radón) para la detección e identificación de las fugas de CO2. En cuanto al análisis químico de las aguas se analizaron los datos geoquímicos e isotópicos y los gases disueltos en las aguas de los alrededores de la PDT de Hontomín, con objeto de determinar qué parámetros son los más apropiados para la detección de una posible migración del CO2 inyectado, o de la salmuera, a los ambientes superficiales. Las medidas de flujo de CO2 se realizaron con la técnica de la cámara de acúmulo. A pesar de ser una técnica desarrollada y aplicada en diferentes ámbitos científicos se estimó necesario adaptar un protocolo de medida y de análisis de datos a las características específicas de los proyectos de captura y almacenamiento de CO2 (CAC). Donde los flujos de CO2 esperados son bajos y en caso de producirse una fuga habrá que detectar pequeñas variaciones en los valores flujo con un “ruido” en la señal alto, debido a actividad biológica en el suelo. La medida de flujo de CO2 mediante la técnica de la cámara de acúmulo se puede realizar sin limpiar la superficie donde se coloca la cámara o limpiando y esperando al reequilibrio del flujo después de la distorsión al sistema. Sin embargo, los resultados obtenidos después de limpiar y esperar muestran menor dispersión, lo que nos indica que este procedimiento es el mejor para la monitorización de los complejos de almacenamiento geológico de CO2. El protocolo de medida resultante, utilizado para la obtención de la línea base de flujo de CO2 en Hontomín, sigue los siguiente pasos: a) con una espátula se prepara el punto de medición limpiando y retirando el recubrimiento vegetal o la primera capa compacta de suelo, b) se espera un tiempo para la realización de la medida de flujo, facilitando el reequilibrio del flujo del gas tras la alteración provocada en el suelo y c) se realiza la medida de flujo de CO2. Una vez realizada la medición de flujo de CO2, y detectada si existen zonas de anomalías, se debe estimar la cantidad de CO2 que se está escapando a la atmósfera (emanación total), con el objetivo de cuantificar la posible fuga. Existen un amplio rango de metodologías para realizar dicha estimación, siendo necesario entender cuáles son las más apropiadas para obtener el valor más representativo del sistema. En esta tesis se comparan seis técnicas estadísticas: media aritmética, estimador insegado de la media (aplicando la función de Sichel), remuestreo con reemplazamiento (bootstrap), separación en diferentes poblaciones mediante métodos gráficos y métodos basados en criterios de máxima verosimilitud, y la simulación Gaussiana secuencial. Para este análisis se realizaron ocho campañas de muestreo, tanto en la Planta de Desarrollo Tecnológico de Hontomón como en análogos naturales (italianos y españoles). Los resultados muestran que la simulación Gaussiana secuencial suele ser el método más preciso para realizar el cálculo, sin embargo, existen ocasiones donde otros métodos son más apropiados. Como consecuencia, se desarrolla un procedimiento de actuación para seleccionar el método que proporcione el mejor estimador. Este procedimiento consiste, en primer lugar, en realizar un análisis variográfico. Si existe una autocorrelación entre los datos, modelizada mediante el variograma, la mejor técnica para calcular la emanación total y su intervalo de confianza es la simulación Gaussiana secuencial (sGs). Si los datos son independientes se debe comprobar la distribución muestral, aplicando la media aritmética o el estimador insesgado de la media (Sichel) para datos normales o lognormales respectivamente. Cuando los datos no son normales o corresponden a una mezcla de poblaciones la mejor técnica de estimación es la de remuestreo con reemplazamiento (bootstrap). Siguiendo este procedimiento el máximo valor del intervalo de confianza estuvo en el orden del ±20/25%, con la mayoría de valores comprendidos entre ±3,5% y ±8%. La identificación de las diferentes poblaciones muestrales en los datos de flujo de CO2 puede ayudar a interpretar los resultados obtenidos, toda vez que esta distribución se ve afectada por la presencia de varios procesos geoquímicos como, por ejemplo, una fuente geológica o biológica del CO2. Así pues, este análisis puede ser una herramienta útil en el programa de monitorización, donde el principal objetivo es demostrar que no hay fugas desde el reservorio a la atmósfera y, si ocurren, detectarlas y cuantificarlas. Los resultados obtenidos muestran que el mejor proceso para realizar la separación de poblaciones está basado en criterios de máxima verosimilitud. Los procedimientos gráficos, aunque existen pautas para realizarlos, tienen un cierto grado de subjetividad en la interpretación de manera que los resultados son menos reproducibles. Durante el desarrollo de la tesis se analizó, en análogos naturales, la relación existente entre el CO2 y los isótopos del radón (222Rn y 220Rn), detectándose en todas las zonas de emisión de CO2 una relación positiva entre los valores de concentración de 222Rn en aire del suelo y el flujo de CO2. Comparando la concentración de 220Rn con el flujo de CO2 la relación no es tan clara, mientras que en algunos casos aumenta en otros se detecta una disminución, hecho que parece estar relacionado con la profundidad de origen del radón. Estos resultados confirmarían la posible aplicación de los isótopos del radón como trazadores del origen de los gases y su aplicación en la detección de fugas. Con respecto a la determinación de la línea base de flujo CO2 en la PDT de Hontomín, se realizaron mediciones con la cámara de acúmulo en las proximidades de los sondeos petrolíferos, perforados en los ochenta y denominados H-1, H-2, H-3 y H-4, en la zona donde se instalarán el sondeo de inyección (H-I) y el de monitorización (H-A) y en las proximidades de la falla sur. Desde noviembre de 2009 a abril de 2011 se realizaron siete campañas de muestreo, adquiriéndose más de 4.000 registros de flujo de CO2 con los que se determinó la línea base y su variación estacional. Los valores obtenidos fueron bajos (valores medios entre 5 y 13 g•m-2•d-1), detectándose pocos valores anómalos, principalmente en las proximidades del sondeo H-2. Sin embargo, estos valores no se pudieron asociar a una fuente profunda del CO2 y seguramente estuvieran más relacionados con procesos biológicos, como la respiración del suelo. No se detectaron valores anómalos cerca del sistema de fracturación (falla Ubierna), toda vez que en esta zona los valores de flujo son tan bajos como en el resto de puntos de muestreo. En este sentido, los valores de flujo de CO2 aparentemente están controlados por la actividad biológica, corroborado al obtenerse los menores valores durante los meses de otoño-invierno e ir aumentando en los periodos cálidos. Se calcularon dos grupos de valores de referencia, el primer grupo (UCL50) es 5 g•m-2•d-1 en las zonas no aradas en los meses de otoño-invierno y 3,5 y 12 g•m-2•d-1 en primavera-verano para zonas aradas y no aradas, respectivamente. El segundo grupo (UCL99) corresponde a 26 g•m-2•d- 1 durante los meses de otoño-invierno en las zonas no aradas y 34 y 42 g•m-2•d-1 para los meses de primavera-verano en zonas aradas y no aradas, respectivamente. Flujos mayores a estos valores de referencia podrían ser indicativos de una posible fuga durante la inyección y posterior a la misma. Los primeros datos geoquímicos e isotópicos de las aguas superficiales (de escorrentía y de manantiales) en el área de Hontomín–Huermeces fueron analizados. Los datos sugieren que las aguas estudiadas están relacionadas con aguas meteóricas con un circuito hidrogeológico superficial, caracterizadas por valores de TDS relativamente bajos (menor a 800 mg/L) y una fácie hidrogeoquímica de Ca2+(Mg2+)-HCO3 −. Algunas aguas de manantiales se caracterizan por concentraciones elevadas de NO3 − (concentraciones de hasta 123 mg/l), lo que sugiere una contaminación antropogénica. Se obtuvieron concentraciones anómalas de of Cl−, SO4 2−, As, B y Ba en dos manantiales cercanos a los sondeos petrolíferos y en el rio Ubierna, estos componentes son probablemente indicadores de una posible mezcla entre los acuíferos profundos y superficiales. El estudio de los gases disueltos en las aguas también evidencia el circuito superficial de las aguas. Estando, por lo general, dominado por la componente atmosférica (N2, O2 y Ar). Sin embargo, en algunos casos el gas predominante fue el CO2 (con concentraciones que llegan al 63% v/v), aunque los valores isotópicos del carbono (<-17,7 ‰) muestran que lo más probable es que esté relacionado con un origen biológico. Los datos geoquímicos e isotópicos de las aguas superficiales obtenidos en la zona de Hontomín se pueden considerar como el valor de fondo con el que comparar durante la fase operacional, la clausura y posterior a la clausura. En este sentido, la composición de los elementos mayoritarios y traza, la composición isotópica del carbono del CO2 disuelto y del TDIC (Carbono inorgánico disuelto) y algunos elementos traza se pueden considerar como parámetros adecuados para detectar la migración del CO2 a los ambientes superficiales. ABSTRACT Since 2009, a group made up of Universidad Politécnica de Madrid (UPM; Spain) and Università degli Studi Firenze (UniFi; Italy) has been taking part in a joint project called “Strategies for Monitoring CO2 and other Gases in Natural analogues”. The group was coordinated by AMPHOS XXI, a private company established in Barcelona. The Project was financially supported by Fundación Ciudad de la Energía (CIUDEN; Spain) as a part of the EC-funded OXYCFB300 project (European Energy Program for Recovery -EEPR-; www.compostillaproject.eu). The main objectives of the project were aimed to develop and optimize analytical methodologies to be applied at the surface to Monitor and Verify the feasibility of geologically stored carbon dioxide. These techniques were oriented to detect and quantify possible CO2 leakages to the atmosphere. Several investigations were made in natural analogues from Spain and Italy and in the Tecnchnological Development Plant for CO2 injection al Hontomín (Burgos, Spain). The studying techniques were mainly focused on the measurements of diffuse soil gases and surface and shallow waters. The soil-gas measurements included the determination of CO2 flux and the application to natural trace gases (e.g. radon) that may help to detect any CO2 leakage. As far as the water chemistry is concerned, geochemical and isotopic data related to surface and spring waters and dissolved gases in the area of the PDT of Hontomín were analyzed to determine the most suitable parameters to trace the migration of the injected CO2 into the near-surface environments. The accumulation chamber method was used to measure the diffuse emission of CO2 at the soil-atmosphere interface. Although this technique has widely been applied in different scientific areas, it was considered of the utmost importance to adapt the optimum methodology for measuring the CO2 soil flux and estimating the total CO2 output to the specific features of the site where CO2 is to be stored shortly. During the pre-injection phase CO2 fluxes are expected to be relatively low where in the intra- and post-injection phases, if leakages are to be occurring, small variation in CO2 flux might be detected when the CO2 “noise” is overcoming the biological activity of the soil (soil respiration). CO2 flux measurements by the accumulation chamber method could be performed without vegetation clearance or after vegetation clearance. However, the results obtained after clearance show less dispersion and this suggests that this procedure appears to be more suitable for monitoring CO2 Storage sites. The measurement protocol, applied for the determination of the CO2 flux baseline at Hontomín, has included the following steps: a) cleaning and removal of both the vegetal cover and top 2 cm of soil, b) waiting to reduce flux perturbation due to the soil removal and c) measuring the CO2 flux. Once completing the CO2 flux measurements and detected whether there were anomalies zones, the total CO2 output was estimated to quantify the amount of CO2 released to the atmosphere in each of the studied areas. There is a wide range of methodologies for the estimation of the CO2 output, which were applied to understand which one was the most representative. In this study six statistical methods are presented: arithmetic mean, minimum variances unbiased estimator, bootstrap resample, partitioning of data into different populations with a graphical and a maximum likelihood procedures, and sequential Gaussian simulation. Eight campaigns were carried out in the Hontomín CO2 Storage Technology Development Plant and in natural CO2 analogues. The results show that sequential Gaussian simulation is the most accurate method to estimate the total CO2 output and the confidential interval. Nevertheless, a variety of statistic methods were also used. As a consequence, an application procedure for selecting the most realistic method was developed. The first step to estimate the total emanation rate was the variogram analysis. If the relation among the data can be explained with the variogram, the best technique to calculate the total CO2 output and its confidence interval is the sequential Gaussian simulation method (sGs). If the data are independent, their distribution is to be analyzed. For normal and log-normal distribution the proper methods are the arithmetic mean and minimum variances unbiased estimator, respectively. If the data are not normal (log-normal) or are a mixture of different populations the best approach is the bootstrap resampling. According to these steps, the maximum confidence interval was about ±20/25%, with most of values between ±3.5% and ±8%. Partitioning of CO2 flux data into different populations may help to interpret the data as their distribution can be affected by different geochemical processes, e.g. geological or biological sources of CO2. Consequently, it may be an important tool in a monitoring CCS program, where the main goal is to demonstrate that there are not leakages from the reservoir to the atmosphere and, if occurring, to be able to detect and quantify it. Results show that the partitioning of populations is better performed by maximum likelihood criteria, since graphical procedures have a degree of subjectivity in the interpretation and results may not be reproducible. The relationship between CO2 flux and radon isotopes (222Rn and 220Rn) was studied in natural analogues. In all emissions zones, a positive relation between 222Rn and CO2 was observed. However, the relationship between activity of 220Rn and CO2 flux is not clear. In some cases the 220Rn activity indeed increased with the CO2 flux in other measurements a decrease was recognized. We can speculate that this effect was possibly related to the route (deep or shallow) of the radon source. These results may confirm the possible use of the radon isotopes as tracers for the gas origin and their application in the detection of leakages. With respect to the CO2 flux baseline at the TDP of Hontomín, soil flux measurements in the vicinity of oil boreholes, drilled in the eighties and named H-1 to H-4, and injection and monitoring wells were performed using an accumulation chamber. Seven surveys were carried out from November 2009 to summer 2011. More than 4,000 measurements were used to determine the baseline flux of CO2 and its seasonal variations. The measured values were relatively low (from 5 to 13 g•m-2•day-1) and few outliers were identified, mainly located close to the H-2 oil well. Nevertheless, these values cannot be associated to a deep source of CO2, being more likely related to biological processes, i.e. soil respiration. No anomalies were recognized close to the deep fault system (Ubierna Fault) detected by geophysical investigations. There, the CO2 flux is indeed as low as other measurement stations. CO2 fluxes appear to be controlled by the biological activity since the lowest values were recorded during autumn-winter seasons and they tend to increase in warm periods. Two reference CO2 flux values (UCL50 of 5 g•m-2•d-1 for non-ploughed areas in autumn-winter seasons and 3.5 and 12 g•m-2•d-1 for in ploughed and non-ploughed areas, respectively, in spring-summer time, and UCL99 of 26 g•m-2•d-1 for autumn-winter in not-ploughed areas and 34 and 42 g•m-2•d-1 for spring-summer in ploughed and not-ploughed areas, respectively, were calculated. Fluxes higher than these reference values could be indicative of possible leakage during the operational and post-closure stages of the storage project. The first geochemical and isotopic data related to surface and spring waters and dissolved gases in the area of Hontomín–Huermeces (Burgos, Spain) are presented and discussed. The chemical and features of the spring waters suggest that they are related to a shallow hydrogeological system as the concentration of the Total Dissolved Solids approaches 800 mg/L with a Ca2+(Mg2+)-HCO3 − composition, similar to that of the surface waters. Some spring waters are characterized by relatively high concentrations of NO3 − (up to 123 mg/L), unequivocally suggesting an anthropogenic source. Anomalous concentrations of Cl−, SO4 2−, As, B and Ba were measured in two springs, discharging a few hundred meters from the oil wells, and in the Rio Ubierna. These contents are possibly indicative of mixing processes between deep and shallow aquifers. The chemistry of the dissolved gases also evidences the shallow circuits of the Hontomín– Huermeces, mainly characterized by an atmospheric source as highlighted by the contents of N2, O2, Ar and their relative ratios. Nevertheless, significant concentrations (up to 63% by vol.) of isotopically negative CO2 (<−17.7‰ V-PDB) were found in some water samples, likely related to a biogenic source. The geochemical and isotopic data of the surface and spring waters in the surroundings of Hontomín can be considered as background values when intra- and post-injection monitoring programs will be carried out. In this respect, main and minor solutes, the isotopic carbon of dissolved CO2 and TDIC (Total Dissolved Inorganic Carbon) and selected trace elements can be considered as useful parameters to trace the migration of the injected CO2 into near-surface environments.

Remediation of crude oil contaminated soils by means of bitumen emulsions

Effective remediation of contaminated sites with oils, heavy metals and other chemical agents is one of the most important environmental problems all over the world. Contaminated soils by petroleum from different sites and origins, might contaminate groundwater aquifers and then be spread by rain. Many years and the requirement of remedial techniques may be needed to remediate them. However, previous experiences show the uselessness of these methods for the solution of all problems. As any case of soil contamination is different, specific studies with relevant factors at financial aspects, legal limits and waste and soil characteristics are needed. In this work, a study of the use of bitumen asphalt emulsion to remediate contaminated soils by crude oils is exposed. By means of the soil stabilization technology, using the contaminated soils as aggregates and the tailor made emulsion as binder, the feasibility of the mix application to produce stable and resistant pavements is demonstrated

Sampling strategies using the "accumulation chamber" for monitoring geological storage of CO2

Fundación Ciudad de la Energía (CIUDEN) is carrying out a project of geological storage of CO2, where CO2 injection tests are planned in saline aquifers at a depth of 1500 m for scientific objectives and project demonstration. Before any CO2 is stored, it is necessary to determine the baseline flux of CO2 in order to detect potential leakage during injection and post-injection monitoring. In November 2009 diffuse flux measurements of CO2 using an accumulation chamber were made in the area selected by CIUDEN for geological storage, located in Hontomin province of Burgos (Spain). This paper presents the tests carried out in order to establish the optimum sampling methodology and the geostatistical analyses performed to determine the range, with which future field campaigns will be planned.

Hidrogeología e hidrodinámica del acuífero de los manantiales de Gormaz y su importancia en el caudal base del río Duero, España = Hydrogeology and hydrodynamic of the Gormaz Springs Aquifer and its importance to the base flow of the Duero River, Spain

Se cuantifican las descargas subterráneas de un acuífero a un río que lo atraviesa utilizando correlaciones estadísticas. El río Duero, España, incrementa su caudal base en varios m3/s, al atravesar unos afloramientos carbonatados mesozoicos en un pequeño tramo de su cabecera; esto es de especial importancia en época de estiaje, cuando la mayor parte del caudal base del río procede de manantiales que allí se sitúan. Dichos afloramientos corresponden a uno de los dos acuíferos calcáreos confinados, que se desarrollan en paralelo y están hidráulicamente desconectados por una capa impermeable, que forman el sistema acuífero de los manantiales de Gormaz. Este sistema se encuentra en estado de régimen natural y está apenas explotado. Se define el modelo conceptual de funcionamiento hidrogeológico, considerando el papel hidrogeológico de la falla de Gormaz, situada en la zona de descarga del sistema. Analizando información geológica antecedente y la geofísica exploratoria realizada, se obtuvo un mejor conocimiento de la geometría y los límites de los acuíferos, definiéndose un sistema acuífero con una zona de recarga en el sur, correspondiente a los afloramientos calcáreos, los cuales se confinan hacia el norte bajo el Terciario, hasta intersecar con la falla normal de Gormaz. El salto de falla genera una barrera para las formaciones permeables situadas al extremo norte (margen derecha del río Duero); a su vez, el plano de falla facilita el ascenso del agua subterránea del sistema acuífero en estudio y pone en conexión hidráulica los dos acuíferos. Se estimaron, además, los parámetros hidráulicos de los acuíferos en los alrededores de la falla. La buena correlación entre los niveles piezométricos y las descargas subterráneas al río Duero han permitido la reconstrucción del hidrograma de los manantiales de Gormaz en el periodo 1992-2006. Se calcula así que la contribución subterránea al río Duero es de 135.9 hm3/año, que supone el 18.9% de la aportación total del río. In a short stretch of its headwaters, the base flow of the River Duero increases by several m3/s as it traverses some Mesozoic carbonate outcrops. This is of special importance during the dry season, when the majority of the base flow of the river proceeds from springs in this reach. The outcrops correspond to one of two confined calcareous aquifers that developed in parallel but which are not hydraulically connected because of an impermeable layer. Together, they constitute the aquifer system of the Gormaz Springs. The system is still in its natural regime and is hardly exploited. This study defines the conceptual model of hydrogeological functioning, taking into consideration the role of the Gormaz Fault, which is situated in the discharge zone of the system. Analysis of both antecedent geological information and geophysical explorations has led to a better understanding of the geometry and boundaries of the aquifers, defining an aquifer system with a recharge zone in the south corresponding to in the calcareous outcrops. These calcareous outcrops are confined to the north below Tertiary formations, as far as their intersection with the normal fault of Gormaz. The throw of the fault forms the barrier of the permeable formations situated in the extreme north (right bank of the River Duero). In turn, the fault plane facilitates the upflow of groundwater from the aquifer system and creates hydraulic connection between the two aquifers. In addition, the study estimated the hydraulic parameters of the aquifer around the fault. The close correlation between piezometric levels and the groundwater discharges to the River Duero has enabled the reconstruction of the hydrogram of Gormaz springs over the period 1992-2006. By this means, it is calculated that the groundwater contribution to the River Duero is 135.9 hm3/year, or 18.9% of the total river inflow.

Residuos orgánicos como fuentes de fósforo

En este trabajo, se han llevado a cabo distintos experimentos en laboratorio, con el objetivo de estudiar el efecto de la aplicación de residuos orgánicos como fuentes de P en las pérdidas de este elemento, que se producen en suelo, tanto por escorrentía superficial como por lixiviación. El interés por evaluar las pérdidas de P se debe a la necesidad de conocer mejor los factores que influyen en los procesos de pérdidas de este elemento y así, poder reducir los problemas de eutrofización de aguas, tanto superficiales como subterráneas, provocadas por un exceso de este nutriente, junto con otros como el nitrógeno. Los trabajos experimentales que se han llevado a cabo se detallan a continuación: Se ha realizado el estudio de las formas de P contenidas en una serie de 14 residuos orgánicos, de distinto origen y tratamiento (compost, lodos, purines y digestato), comparando la información aportada por tres protocolos de fraccionamientos de P diferentes, seleccionados entre los principales métodos utilizados: protocolo de fraccionamiento de Ruttemberg (1992), protocolo de Normas, medidas y ensayos (Ruban et al., 2001a) y protocolo de Huang et al. (2008). Todos los métodos de fraccionamiento empleados aportaron información útil para conocer las formas de P de los residuos, a pesar de que alguno de ellos fue descrito para sedimentos o suelos. Sin embargo, resulta difícil comparar los resultados entre unos y otros, ya que cada uno emplea extractantes y tiempos de ensayos diferentes. Las cantidades de P total determinadas por cada método mantienen una relación lineal, aunque el método SMT, por ser más directo, obtiene las cantidades más elevadas para todos los residuos. Los métodos secuenciales (métodos de Huang y Ruttemberg), a pesar de ser más largos y tediosos, aportan información más detallada de la disponibilidad de las formas de P, y con ello, permiten obtener una mejor estimación de las pérdidas potenciales de este elemento tras su aplicación en suelo. Se han encontrado relaciones positivas entre las fracciones determinadas por estos dos métodos. Así mismo, se encuentra una correlación entre las formas solubles de P y la concentración de los iones [Ca + Fe + Al] de los residuos, útiles como indicadores de la disponibilidad de este elemento. Sin embargo, el protocolo SMT, no obtiene información de la solubilidad o disponibilidad de las formas de P contenidas, ni las formas determinadas mantienen relaciones directas con las de los otros métodos, con excepción del P total. Para el estudio del comportamiento de los residuos aplicados en suelos, se pusieron a punto sistemas de lluvia simulada, con el objetivo de caracterizar las pérdidas de P en la escorrentía superficial generada. Por otra parte, se emplearon columnas de suelos enmendados con residuos orgánicos, para el estudio de las pérdidas de P por lixiviación. Los ensayos de simulación de lluvia se llevaron a cabo de acuerdo al “National Phosphorus Research proyect“ (2001), que consigue simular eventos sucesivos de lluvia en unas condiciones semejantes a la realidad, empleando cajas llenas de suelo del horizonte superficial, con residuos aplicados tanto superficialmente como mediante mezcla con el propio suelo. Los ensayos se realizaron con seis residuos de diferente naturaleza y sometidos a distintos tratamientos. Se encontraron diferencias significativas en las pérdidas de las formas de P analizadas, tanto disueltas como particuladas, en las aguas de escorrentía generadas. En general, las pérdidas en el primer evento de lluvia tras la aplicación de los residuos fueron mayores a las generadas en el segundo evento, predominando las formas de P particuladas respecto a las disueltas en ambos. Se encontró una relación positiva entre las pérdidas de P en las aguas de escorrentía generadas en cada ensayo, con los contenidos de P soluble en agua y fácilmente disponible de los residuos empleados, determinados por los protocolos de fraccionamientos secuenciales. Además, se emplearon los modelos matemáticos desarrollados por Vadas et al. (2005, 2007), de evaluación de las pérdidas de P por escorrentía para fertilizantes y estiércoles. La predicción de estos modelos no se cumple en el caso de todos los residuos. Las distintas propiedades físicas de los residuos pueden afectar a las diferencias entre las pérdidas experimentales y las esperadas. Los ensayos de simulación del proceso de lixiviación se llevaron a cabo en columnas de percolación, con suelos enmendados con residuos orgánicos, de acuerdo a la norma “CEN/TS 14405–2004: Caracterización de los residuos – Test de comportamiento de lixiviación – Test de flujo ascendente”. Las pérdidas de P por procesos de lixiviación de agua, han sido despreciadas durante mucho tiempo respecto a las pérdidas por escorrentía. Sin embargo, se ha demostrado que deben tenerse en consideración, principalmente en algunos tipos de suelos o zonas cercanas a acuíferos. Se utilizaron tres suelos de distinta procedencia para los ensayos, de manera que se pudo estudiar la influencia del tipo de suelo en las pérdidas de P para cada tipo de residuo (purín, compost, digestato y lodo de EDAR). Los índices de adsorción de P determinados para cada suelo permiten evaluar aquellos que presentarán más riesgo de producir pérdidas de este elemento al aplicarse fuentes externas de P, encontrando una relación positiva entre ambos. Las pérdidas de P en los lixiviados varían en función tanto del residuo como del suelo empleado. Para el compost, el purín y el lodo, se encontró una relación entre las pérdidas generadas en el agua lixiviada de las columnas y las formas de P soluble contenidas en los residuos. Sin embargo, en el caso del digestato, no existía esta correlación. Las pérdidas para este residuo fueron en todos los casos menores a las estimadas, considerando las formas de P contenido. El estudio de la mojabilidad, propiedad física del residuo que evalúa la capacidad de interacción residuo-agua, permitió explicar el comportamiento anómalo de este residuo, con una mayor resistencia a que el agua entrara en su estructura y por tanto, una mayor dificultad de solubilizar el P contenido en el propio residuo, que en el caso de otros residuos. En general, podemos considerar que el estudio de las formas de P más disponibles o solubles en agua, aporta información útil de las pérdidas potenciales de P. Sin embargo, es necesario estudiar las propiedades físicas de los residuos orgánicos aplicados y la capacidad de adsorción de P de los suelos, para estimar las pérdidas de P y con ello, colaborar a controlar los procesos de eutrofización en aguas. ABSTRACT This dissertation explores the effect of organic wastes application as sources of P in losses of this element that occur by both surface runoff and leaching in soil. To do so, diverse laboratory experiments are conducted and presented here. Evaluating P losses is necessary to better understand the factors that influence the processes behind the loss of this element. Reducing P losses reduces eutrophication problems of both surface water and groundwater caused by an excess of this nutrient, along with other as nitrogen. Details of the experiments are presented below: The first experiment studies the forms of P contained in a series of 14 organic wastes of different origin and treatment (compost, sludge, slurry and digestate), comparing the information provided by three methods of P fractionation. The methods selected were: Ruttemberg protocol (1992); Standards, Measurements and Testing protocol (Ruban et al., 2001a); and Huang protocol (Huang et al., 2008). All fractionation methods employed successfully contribute to our knowledge of P forms in wastes, even though one of them was originally described for sediments or soils information. However, it is difficult to compare results among each other, as each protocol employs different extractants and time in the trials. Total amounts of P obtained by each method show a linear relationship, although the SMT method, which is more direct, obtains the highest amounts for all residues. Sequential methods (Huang and Ruttemberg’s protocols), despite being longer and more tedious, provide more detailed information on the availability of the forms of P. Therefore, allow the estimation of the potential losses of P after application in soil. Furthermore, positive relationships have been found among fractions obtained by these methods. Positive relationship has been found also among soluble forms of P and the concentration of ions Fe + Ca + Al, which is useful as an indicator of the availability of this element. However, the SMT protocol does not collect information about solubility or availability of forms of P contained; neither do certain forms maintain direct relations with the forms from other methods, with the exception of total P methods. To study the behavior of wastes applied to soils two experiments were conducted. Simulated rain systems were prepared to characterize P losses in the surface runoff generated. In addition, columns of soils amended with organic waste were developed for the study of P leaching losses. Simulated rain systems were carried out according to the ’National Phosphorus Research Project’ (2001), which manages to simulate successive rainfall events in conditions resembling reality. The experiment uses boxes filled with soil from the surface horizon amended with residues, both superficially and by mixing with the soil. Tests were conducted with six residues of different type and subjected to diverse treatments. Findings show significant differences in losses of the P forms analyzed in the generated runoff water, in both solution and particulate forms. In general, losses in the first rainfall event after application of waste were higher than the losses generated in the second event, predominating particulate forms of P over dissolved forms in both events. In all trials, a positive relationship was found between various P forms determined by sequential fractionation protocols (water soluble P and readily available P forms) and P losses in runoff. Furthermore, results from Vadas´s mathematical models (Vadas et al., 2005; 2007) to assess P losses by runoff fertilizers and manures indicate that the prediction of this model is not fulfilled in the case of all residues. The diverse physical properties of wastes may affect the differences between experimental and expected losses. Finally, leaching simulation processes were carried out in percolation columns, filled with soils amended with organic wastes, following the ‘CEN/TS 14405-2004 standard: Characterization of waste - Leaching behavior test - Test Flow ascending ’. P losses by leaching have been neglected for a long time with respect to runoff losses. However, findings corroborate previous studies showing that these P losses have to be taken into account, especially in certain types of soils and in zones near aquifers. To study the influence of soil type on P losses, experiments were carried out with three different soils and for each type of waste (manure, compost, digestate and sludge WWTP). Each soil’s P adsorption rates allow assessing which soils imply a higher risk of P losses when external sources of P are applied. P losses in leachate vary according to the type of soil employed and according to the specific residue. In the case of compost, manure and sludge, there is a relationship between leaching losses and residues’ soluble forms of P. The exception being the digestate, where there was no such correlation. Digestate P losses by leaching were lower than expected in all cases considering the forms of P contained. Moreover, examining digestate wettability -- that is, the physical property of the residue that assesses the capacity of waste-water interaction -- allowed explaining the anomalous behavior of this residue. Digestate has a high resistance to water entering its structure and thus higher difficulty to solubilize the P contained. Overall, studying the more available or soluble P forms provides useful information about the potential loss of P. However, this dissertation shows that it is necessary to examine the physical properties of organic residues applied as well as the P adsorption capacity of soils to estimate P losses, and thus to control eutrophication in water.