Geochemistry of igneous rocks recovered from a transect across the Mariana Trough, arc, fore-arc, and trench, at DSDP Leg 60 Holes and DM17

Major and trace element analyses are presented for 110 samples from the DSDP Leg 60 basement cores drilled along a transect across the Mariana Trough, arc, fore-arc, and Trench at about 18°N. The igneous rocks forming breccias at Site 453 in the west Mariana Trough include plutonic cumulates and basalts with calc-alkaline affinities. Basalts recovered from Sites 454 and 456 in the Mariana Trough include types with compositions similar to normal MORB and types with calc-alkaline affinities within a single hole. At Site 454 the basalts show a complete compositional transition between normal MORB and calc-alkaline basalts. These basalts may be the result of mixing of the two magma types in small sub-crustal magma reservoirs or assimilation of calc-alkaline, arc-derived vitric tuffs by normal MORB magmas during eruption or intrusion. A basaltic andesite clast in the breccia recovered from Site 457 on the active Mariana arc and samples dredged from a seamount in the Mariana arc are calc-alkaline and similar in composition to the basalts recovered from the Mariana Trough and West Mariana Ridge. Primitive island arc tholeiites were recovered from all four sites (Sites 458-461) drilled on the fore-arc and arc-side wall of the trench. These basalts form a coherent compositional group distinct from the Mariana arc, West Mariana arc, and Mariana Trough calc-alkaline lavas, indicating temporal (and perhaps spatial?) chemical variations in the arc magmas erupted along the transect. Much of the 209 meters of basement cored at Site 458 consists of endiopside- and bronzite-bearing, Mg-rich andesites with compositions related to boninites. These andesites have the very low Ti, Zr, Ti/Zr, P, and rare-earthelement contents characteristic of boninites, although they are slightly light-rare-earth-depleted and have lower MgO, Cr, Ni, and higher CaO and Al2O3 contents than those reported for typical boninites. The large variations in chemistry observed in the lavas recovered from this transect suggest that diverse mantle source compositions and complex petrogenetic process are involved in forming crustal rocks at this intra-oceanic active plate margin.

C and H geochemistry of altered oceanic crus of ODP/IODP Hole 1256D

Carbon and hydrogen concentrations and isotopic compositions were measured in 19 samples from altered oceanic crust cored in ODP/IODP Hole 1256D through lavas, dikes down to the gabbroic rocks. Bulk water content varies from 0.32 to 2.14 wt% with dD values from -64per mil to -25per mil. All samples are enriched in water relative to fresh basalts. The dD values are interpreted in terms of mixing between magmatic water and another source that can be either secondary hydrous minerals and/or H contained in organic compounds such as hydrocarbons. Total CO2, extracted by step-heating technique, ranges between 564 and 2823 ppm with d13C values from -14.9per mil to -26.6per mil. As for water, these altered samples are enriched in carbon relative to fresh basalts. The carbon isotope compositions are interpreted in terms of a mixing between two components: (1) a carbonate with d13C = -4.5per mil and (2) an organic compound with d13C = -26.6per mil. A mixing model calculation indicates that, for most samples (17 of 19), more than 75% of the total C occurs as organic compounds while carbonates represent less than 25%. This result is also supported by independent estimates of carbonate content from CO2 yield after H3PO4 attack. A comparison between the carbon concentration in our samples, seawater DIC (Dissolved Inorganic Carbon) and DOC (Dissolved Organic Carbon), and hydrothermal fluids suggests that CO2 degassed from magmatic reservoirs is the main source of organic C addition to the crust during the alteration process. A reduction step of dissolved CO2 is thus required, and can be either biologically mediated or not. Abiotic processes are necessary for the deeper part of the crust (>1000 mbsf) because alteration temperatures are greater than any hyperthermophilic living organism (i.e. T > 110 °C). Even if not required, we cannot rule out the contribution of microbial activity in the low-temperature alteration zones. We propose a two-step model for carbon cycling during crustal alteration: (1) when "fresh" oceanic crust forms at or close to ridge axis, alteration starts with hot hydrothermal fluids enriched in magmatic CO2, leading to the formation of organic compounds during Fischer-Tropsch-type reactions; (2) when the crust moves away from the ridge axis, these interactions with hot hydrothermal fluids decrease and are replaced by seawater interactions with carbonate precipitation in fractures. Taking into account this organic carbon, we estimate C isotope composition of mean altered oceanic crust at ? -4.7per mil, similar to the d13C of the C degassed from the mantle at ridge axis, and discuss the global carbon budget. The total flux of C stored in the altered oceanic crust, as carbonate and organic compound, is 2.9 ± 0.4 * 10**12 molC/yr.

Grain size record of Oligocene to Pleistocene sediments from Broken Ridge

At Ocean Drilling Program Sites 752 and 754, located on Broken Ridge in the eastern Indian Ocean, we recovered a sequence of shallow-water pelagic sediments that span the past 90 m.y. The Oligocene to Pleistocene portion of these sediments are unconsolidated carbonate oozes that display a coherent variation in bulk grain size. We believe these sediments to be winnowed, and suggest that their grain size is a measure of that winnowing energy. The largest increase in grain size, interpreted to represent an enhancement in the energy of ocean currents, occurs in the earliest late Miocene. This increase occurs about 20 m upcore from the oxygen isotope indication of ice-volume increase about 13 Ma, and is about 3 m.y. younger. If this distinct temporal separation between proxy indicators of ice volume and of current intensity observed in the Broken Ridge cores is correct, the general impression of paleoclimatologists that the planetary temperature gradient and therefore atmospheric and oceanic circulation intensity varies directly with ice volume needs to be reconsidered.

Upper-ocean microstructure data from the tropical Atlantic

SST variability within the Atlantic cold tongue (ACT) region is of climatic relevance for the surrounding continents. A multi cruise data set of microstructure observations is used to infer regional as well as seasonal variability of upper ocean mixing and diapycnal heat flux within the ACT region. The variability in mixing intensity is related to the variability in large scale background conditions, which were additionally observed during the cruises. The observations indicate fundamental differences in background conditions in terms of shear and stratification below the mixed layer (ML) for the western and eastern equatorial ACT region causing critical Froude numbers (Fr) to be more frequently observed in the western equatorial ACT. The distribution of critical Fr occurrence below the ML reflects the regional and seasonal variability of mixing intensity. Turbulent dissipation rates (?) at the equator (2°N-2°S) are strongly increased in the upper thermocline compared to off-equatorial locations. In addition, ? is elevated in the western equatorial ACT compared to the east from May to November, whereas boreal summer appears as the season of highest mixing intensities throughout the equatorial ACT region, coinciding with ACT development. Diapycnal heat fluxes at the base of the ML in the western equatorial ACT region inferred from ? and stratification range from a maximum of 90 Wm-2 in boreal summer to 55 Wm-2 in September and 40 Wm-2 in November. In the eastern equatorial ACT region maximum values of about 25 Wm-2 were estimated during boreal summer reducing to about 5 Wm-2 towards the end of the year. Outside the equatorial region, inferred diapycnal heat fluxes are comparably low rarely exceeding 10 Wm-2. Integrating the obtained heat flux estimates in the ML heat budget at 10°W on the equator accentuates the diapycnal heat flux as the largest ML cooling term during boreal summer and early autumn. In the western equatorial ACT elevated meridional velocity shear in the upper thermocline contributes to the enhanced diapycnal heat flux within this region during boreal summer and autumn. The elevated meridional velocity shear appears to be associated with intra-seasonal wave activity.

Thermodynamics of zinc insertion in CuGaS2:Ti, used as a modulator agent in an intermediate-band photovoltaic material

An intermediate-bandphotovoltaicmaterial, which has an isolated metallic band located between the top of the valence band and bottom of the conduction band of some semiconductors, has been proposed as third generation solar cell to be used in photovoltaic applications. Density functional theory calculations of Zn in CuGaS2:Ti have previously shown that, the intermediate-band position can be modulated in proportion of Zn insertion in such a way that increasing Zn concentration can lead to aband-gap reduction, and an adjustment of the intermediate-band position. This could be interesting in the formation of an intermediate-bandmaterial, that has the maximum efficiency theoretically predicted for the intermediate-band solar cell. In this work, the energetics of several reaction schemes that could lead to the decomposition of the modulated intermediate-bandphotovoltaicmaterial, CuGaS2:Ti:Zn, is studied in order to assess the thermodynamic stability of this material. Calculations of the total free energy and disorder entropy have been taken into account, to get the reaction energy and free energy of the compound decomposition, which is found to be thermodynamically favorable

Review of the Performance of Residential PV systems in Belgium

The main objective of this paper is to review the state of the art of residential PV systems in Belgium by the analysis of the operational data of 993 installations. For that, three main questions are posed: how much energy do they produce? What level of performance is associated to their production? Which are the key parameters that most influence their quality? This work brings answers to these questions. A middling commercial PV system, optimally oriented, produces a mean annual energy of 892 kWh/kWp. As a whole, the orientation of PV generators causes energy productions to be some 6% inferior to optimally oriented PV systems. The mean performance ratio is 78% and the mean performance index is 85%. That is to say, the energy produced by a typical PV system in Belgium is 15% inferior to the energy produced by a very high quality PV system. Finally, on average, the real power of the PV modules falls 5% below its corresponding nominal power announced on the manufacturer's datasheet. Differences between real and nominal power of up to 16% have been detected.

Review of the Performance of Residential PV systems in France

The main objective of this paper is to review the state of the art of residential PV systems in France. This is done analyzing the operational data of 6868 installations. Three main questions are posed. How much energy do they produce? What level of performance is associated to their production? Which are the key parameters that most influence their quality? During the year 2010, the PV systems in France have produced a mean annual energy of 1163 kWh/kWp. As a whole, the orientation of PV generators causes energy productions to be some 7% inferior to optimally oriented PV systems. The mean Performance Ratio is 76% and the mean Performance Index is 85%. That is to say, the energy produced by a typical PV system in France is 15% inferior to the energy produced by a very high quality PV system. On average, the real power of the PV modules falls 4.9% below its corresponding nominal power announced on the manufacturer's datasheet. A brief analysis by PV modules technology has led to relevant observations about two technologies in particular. On the one hand, the PV systems equipped with heterojunction with intrinsic thin layer (HIT) modules show performances higher than average. On the other hand, the systems equipped with the copper indium (di)selenide (CIS) modules show a real power that is 16% lower than their nominal value.

Equation of State and Opacities for Warm Dense Matter.

We will present recent developments in the calculation of opacity and equation of state tables suitable for including in the radiation hydrodynamic code ARWEN [1] to study processes like ICF or X-ray secondary sources. For these calculations we use the code BiG BART to compute opacities in LTE conditions, with self-consistent data generated with the Flexible Atomic Code (FAC) [2]. Non-LTE effects are approximately taken into account by means of the improved RADIOM model [3], which makes use of existing LTE data tables. We use the screened-hydrogenic model [4] to derive the Equation of State using the population and energy of the levels avaliable from the atomic data

Performance Analysis of 10,000 Residential PV Systems in France and Belgium

The main objective of this paper is to review the state of the art of residential PV systems in France and Belgium. This is done analyzing the operational data of 10650 PV systems (9657 located in France and 993 in Belgium). Three main questions are posed. How much energy do they produce? What level of performance is associated to their production? Which are the key parameters that most influence their quality? During the year 2010, the PV systems in France have produced a mean annual energy of 1163 kWh/kWp in France and 852 kWh/kWp in Belgium. As a whole, the orientation of PV generators causes energy productions to be some 7% inferior to optimally oriented PV systems. The mean Performance Ratio is 76% in France and 78% in Belgium, and the mean Performance Index is 85% in both countries. On average, the real power of the PV modules falls 4.9% below its corresponding nominal power announced on the manufacturer?s datasheet. A brief analysis by PV modules technology has lead to relevant observations about two technologies in particular. On the one hand, the PV systems equipped with Heterojunction with Intrinsic. Thin layer (HIT) modules show performances higher than average. On the other hand, the systems equipped with Copper Indium (di)Selenide (CIS) modules show a real power that is 16 % lower than their nominal value.

Development of corrosion of steel bars embedded in mortar made with slag from secondary metallurgy

The aim of this work is to study the evolution of the corrosion rate of reinforcements embedded in mortar specimens that have been partly or fully replaced by the sand ladle furnace white slag. Prisms are manufactured mortar 6cm x 8cm x 2cm in which are embedded reinforcing steel bars of 6mm diameter B500SD. At the time of mixing were added varying amounts of chloride ion content by weight of cement (0%, 0.4%, 0.8%, 1.2%, 2%). The specimens were made totally or partially replacing the white slag, getting four different mixes depending on the degree of substitution. After curing the specimens for 28 days in moist chambers proceeded to dry up naturally. Here are gradually dampened by its conservation in a moist chamber, periodically measuring the corrosion rate of the bars using the technique of polarization curve. The results, in terms of corrosion current and corrosion potential, were compared with those obtained on standard samples, without replacement by slag aggregate. The analysis of results allows us to know, depending on the type of mortar used, the chloride threshold with the depassivation produced steel and the corrosion rates achieved in steels in the active state in terms of mortar moisture, obtained from qualitatively using gravimetric techniques. The results achieved to date support the conclusion that no significant differences in the behavior against corrosion induced by chloride ions, between the steel bars embedded in standard samples and the steel bars embedded in samples including with aggregates from slag. Both the chloride threshold resulting in the depassivation steel as the corrosion rate reached through the bars in an active state are very similar in both types of mortars when they have the same moisture content.

Numerical simulation of Ge solar cells using D-AMPS-1D code

A solar cell is a solid state device that converts the energy of sunlight directly into electricity by the photovoltaic effect. When light with photon energies greater than the band gap is absorbed by a semiconductor material, free electrons and free holes are generated by optical excitation in the material. The main characteristic of a photovoltaic device is the presence of internal electric field able to separate the free electrons and holes so they can pass out of the material to the external circuit before they recombine. Numerical simulation of photovoltaic devices plays a crucial role in their design, performance prediction, and comprehension of the fundamental phenomena ruling their operation. The electrical transport and the optical behavior of the solar cells discussed in this work were studied with the simulation code D-AMPS-1D. This software is an updated version of the one-dimensional (1D) simulation program Analysis of Microelectronic and Photonic Devices (AMPS) that was initially developed at The Penn State University, USA. Structures such as homojunctions, heterojunctions, multijunctions, etc., resulting from stacking layers of different materials can be studied by appropriately selecting characteristic parameters. In this work, examples of cells simulation made with D-AMPS-1D are shown. Particularly, results of Ge photovoltaic devices are presented. The role of the InGaP buffer on the device was studied. Moreover, a comparison of the simulated electrical parameters with experimental results was performed.

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.

Laser Shock Processing influence on local properties and overall tensile behavior of friction stir welded joints

Based on laser beam intensities above 109 W/cm2 with pulse energy of several Joules and duration of nanoseconds, Laser Shock Processing (LSP) is capable of inducing a surface compressive residual stress field. The paper presents experimental results showing the ability of LSP to improve the mechanical strength and cracking resistance of AA2024-T351 friction stir welded (FSW) joints. After introducing the FSW and LSP procedures, the results of microstructural analysis and micro-hardness are discussed. Video Image Correlation was used to measure the displacement and strain fields produced during tensile testing of flat specimens; the local and overall tensile behavior of native FSW joints vs. LSP treated were analyzed. Further, results of slow strain rate tensile testing of the FSW joints, native and LSP treated, performed in 3.5% NaCl solution are presented. The ability of LSP to improve the structural behavior of the FSW joints is underscored.

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.