CO2 fluxes of alpine shrubland ecosystem on the north-eastern Tibetan plateau

We measured ecosystem CO2 fluxes for an alpine shrubland on the north-eastern Tibetan Plateau, Qinghai, China. The study is to understand (1) the seasonal variation of CO2 flux and (2) how environmental factors affect the seasonality of CO2 exchange in the alpine ecosystem. Daytime ecosystem respiration was extrapolated from the relationship between temperature and nighttime CO2 fluxes under high turbulent conditions.Seasonal patterns of gross ecosystem production, ecosystem respiration and net ecosystem CO2 exchange followed highly the seasonal change of aboveground biomass in the alpine shrubland. The net ecosystem CO2 exchange was mainly controlled by the variation of photosynthetic photon flux density, while the ecosystem respiration was closely correlated to the soil temperature at 5-cm depth. Integrated values of gross ecosystem production, ecosystem respiration and net ecosystem CO2 exchange for the period from November 1, 2002 to October 31 2003 were estimated to be 1418, 1155 and 222 g CO2 m(-2) yr(-1), respectively.

Diurnal, seasonal and annual variation in net ecosystem CO2 exchange of an alpine shrubland on Qinghai-Tibetan plateau

Thus far, grassland ecosystem research has mainly been focused on low-lying grassland areas, whereas research on high-altitude grassland areas, especially on the carbon budget of remote areas like the Qinghai-Tibetan plateau is insufficient. To address this issue, flux of CO2 were measured over an alpine shrubland ecosystem (37 degrees 36'N, 101 degrees 18'E; 325 above sea level [a. s. l.]) on the Qinghai-Tibetan Plateau, China, for 2 years (2003 and 2004) with the eddy covariance method. The vegetation is dominated by formation Potentilla fruticosa L. The soil is Mol-Cryic Cambisols. To interpret the biotic and abiotic factors that modulate CO2 flux over the course of a year we decomposed net ecosystem CO2 exchange (NEE) into its constituent components, and ecosystem respiration (R-eco). Results showed that seasonal trends of annual total biomass and NEE followed closely the change in leaf area index. Integrated NEE were -58.5 and -75.5 g C m(-2), respectively, for the 2003 and 2004 years. Carbon uptake was mainly attributed from June, July, August, and September of the growing season. In July, NEE reached seasonal peaks of similar magnitude (4-5 g C m(-2) day(-1)) each of the 2 years. Also, the integrated night-time NEE reached comparable peak values (1.5-2 g C m(-2) day(-1)) in the 2 years of study. Despite the large difference in time between carbon uptake and release (carbon uptake time < release time), the alpine shrubland was carbon sink. This is probably because the ecosystem respiration at our site was confined significantly by low temperature and small biomass and large day/night temperature difference and usually soil moisture was not limiting factor for carbon uptake. In general, R-eco was an exponential function of soil temperature, but with season-dependent values of Q(10). The temperature-dependent respiration model failed immediately after rain events, when large pulses of R-eco were observed. Thus, for this alpine shrubland in Qinghai-Tibetan plateau, the timing of rain events had more impact than the total amount of precipitation on ecosystem R-eco and NEE.

Characterizing CO2 fluxes for growing and non-growing seasons in a shrub ecosystem on the Qinghai-Tibet Plateau

To assess carbon budget for shrub ecosystems on the Qinghai-Tibet Plateau, CO2 flux was measured with an open-path eddy covariance system for an alpine shrub ecosystem during growing and non-growing seasons. CO2 flux dynamics was distinct between the two seasons. During the growing season from May to September, the ecosystem exhibited net CO2 uptake from 08:00 to 19:00 (Beijing Standard Time), but net CO2 emission from 19:00 to 08:00. Maximum CO2 uptake appeared around 12:00 with values of 0.71, 1,19, 1.46 and 0.67 g CO2 m(-2) h(-1) for June, July, August and September, respectively. Diurnal fluctuation Of CO2 flux showed higher correlation with photosynthetic photon flux density than temperature. The maximum net CO2 influx occurred in August with a value of 247 g CO2 m(-2). The total CO2 uptake by the ecosystem was up to 583 g CO2 m(-2) for the growing season. During the non-growing season from January to April and from October to December, CO2 flux showed small fluctuation with the largest net CO2 efflux of 0.30 g CO2 m(-2) h(-1) in April. The diurnal CO2 flux was close to zero during most time of the day, but showed a small net CO2 eff lux from 11:00 to 18:00. Diurnal CO2 flux, is significantly correlated to diurnal temperature in the non-growing season. The maximum monthly net CO2 eff lux appeared in April, with a value of 105 g CO2 m(-2). The total net CO2 eff lux for the whole non-growing season was 356 g CO2 m(-2).

Coccolithophore surface distributions in the North Atlantic and their modulation of the air-sea flux of CO2 from 10 years of satellite Earth observation data

Coccolithophores are the primary oceanic phytoplankton responsible for the production of calcium carbonate (CaCO3). These climatically important plankton play a key role in the oceanic carbon cycle as a major contributor of carbon to the open ocean 5 carbonate pump (�50%) and their formation can affect the atmosphere-to-ocean (airsea) uptake of carbon dioxide (CO2) through increasing the seawater partial pressure of CO2 (pCO2). Here we document variations in the areal extent of surface blooms of the globally important coccolithophore, Emiliania huxleyi, in the North Atlantic over a 10-year period (1998–2007), using Earth observation data from the Sea-viewing Wide 10 Field of view Sensor (SeaWiFS).We calculate the annual mean surface areal coverage of E. huxleyi in the North Atlantic to be 474 000±119 000km2 yr−1, which results in a net CaCO3 production of 0.62±0.15 Tg CaCO3 carbon per year. However, this surface coverage and net production can fluctuate by −54/+81% about these mean values and are strongly correlated with the El Ni˜no/Southern Oscillation (ENSO) climate os15 cillation index (r =0.75, p<0.02). Our analysis evaluates the spatial extent over which the E. huxleyi blooms in the North Atlantic can increase the pCO2 and thus decrease the localised sink of atmospheric CO2. In regions where the blooms are prevalent, the average reduction in the monthly CO2 sink can reach 12 %. The maximum reduction of the monthly CO2 sink in the time series is 32 %. This work suggests that the high 20 variability, frequency and distribution of these calcifying plankton and their impact on pCO2 should be considered within modelling studies of the North Atlantic if we are to fully understand the variability of its air-to-sea CO2 flux.

Air-sea CO2 exchange in the Kuroshio and its importance to the global CO2 uptake

Senior thesis written for Oceanography 444

Estimating surface CO2 fluxes from space-borne dry air mole fraction observations using an Ensemble Kalman filter

We have developed an ensemble Kalman Filter (EnKF) to estimate 8-day regional surface fluxes of CO2 from space-borne CO2 dry-air mole fraction observations (XCO2) and evaluate the approach using a series of synthetic experiments, in preparation for data from the NASA Orbiting Carbon Observatory (OCO). The 32-day duty cycle of OCO alternates every 16 days between nadir and glint measurements of backscattered solar radiation at short-wave infrared wavelengths. The EnKF uses an ensemble of states to represent the error covariances to estimate 8-day CO2 surface fluxes over 144 geographical regions. We use a 12×8-day lag window, recognising that XCO2 measurements include surface flux information from prior time windows. The observation operator that relates surface CO2 fluxes to atmospheric distributions of XCO2 includes: a) the GEOS-Chem transport model that relates surface fluxes to global 3-D distributions of CO2 concentrations, which are sampled at the time and location of OCO measurements that are cloud-free and have aerosol optical depths <0.3; and b) scene-dependent averaging kernels that relate the CO2 profiles to XCO2, accounting for differences between nadir and glint measurements, and the associated scene-dependent observation errors. We show that OCO XCO2 measurements significantly reduce the uncertainties of surface CO2 flux estimates. Glint measurements are generally better at constraining ocean CO2 flux estimates. Nadir XCO2 measurements over the terrestrial tropics are sparse throughout the year because of either clouds or smoke. Glint measurements provide the most effective constraint for estimating tropical terrestrial CO2 fluxes by accurately sampling fresh continental outflow over neighbouring oceans. We also present results from sensitivity experiments that investigate how flux estimates change with 1) bias and unbiased errors, 2) alternative duty cycles, 3) measurement density and correlations, 4) the spatial resolution of estimated flux estimates, and 5) reducing the length of the lag window and the size of the ensemble. At the revision stage of this manuscript, the OCO instrument failed to reach its orbit after it was launched on 24 February 2009. The EnKF formulation presented here is also applicable to GOSAT measurements of CO2 and CH4.

An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.

Fractal dimension and anisotropy of soil CO2 emission in an agricultural field during fallow

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Variabilidade espacial da emissão de CO2, temperatura e umidade em latossolo sob cultivo de cana-de-açúcar em sistemas de colheita manual com queima e mecanizada

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Variabilidade espacial da emissão de CO2 e sua relação com propriedades do solo em área de cana-de-açúcar no Sudeste do Brasil

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Iron oxides as proxies for characterizing anisotropy in soil CO2 emission in sugarcane areas under green harvest

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

On the spatial and temporal dependence of CO2 emission on soil properties in sugarcane (Saccharum spp.) production

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

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.

Seasonal differences in the CO2 exchange of a short-hydroperiod Florida Everglades marsh

Although wetlands are among the world's most productive ecosystems, little is known of long-term CO2 exchange in tropical and subtropical wetlands. The Everglades is a highly managed wetlands complex occupying >6000 km2 in south Florida. This ecosystem is oligotrophic, but extremely high rates of productivity have been previously reported. To evaluate CO2 exchange and its response to seasonality (dry vs. wet season) in the Everglades, an eddy covariance tower was established in a short-hydroperiod marl marsh. Rates of net ecosystem exchange and ecosystem respiration were small year-round and declined in the wet season relative to the dry season. Inundation reduced macrophyte CO2 uptake, substantially limiting gross ecosystem production. While light and air temperature exerted the primary controls on net ecosystem exchange and ecosystem respiration in the dry season, inundation weakened these relationships. The ecosystem shifted from a CO2 sink in the dry season to a CO2 source in the wet season; however, the marsh was a small carbon sink on an annual basis. Net ecosystem production, ecosystem respiration, and gross ecosystem production were −49.9, 446.1 and 496.0 g C m−2 year−1, respectively. Unexpectedly low CO2 flux rates and annual production distinguish the Everglades from many other wetlands. Nonetheless, impending changes in water management are likely to alter the CO2 balance of this wetland and may increase the source strength of these extensive short-hydroperiod wetlands.

IMPACTS OF ENHANCED TEMPERATURE AND SNOW DEPOSITION ON SENESCENCE DATE, VEGETATION COVER, AND CO2 EXCHANGE IN A CANADIAN HIGH ARCTIC MESIC ECOSYSTEM

Arctic regions are expected to experience an increase in both temperature and precipitation over the coming decades, which is likely to impact vegetation dynamics and greenhouse gas exchange. To test this response, an experiment was installed at the Cape Bounty Arctic Watershed Observatory, on Melville Island, NU, in 2008 as part of the International Tundra Experiment (ITEX). Snow fences and open top chambers (OTCs) were used to manipulate snow depth and air temperature, respectively. Unlike most ITEX sites to date, enhanced temperature and snowfall were combined here in a factorial design with eight replicates. As an added control, four plots were established well outside the enhanced snow area. Senescence date was recorded at the end of the season, and at the peak of the growing season a vegetation survey was conducted within each plot in order to determine the total percent cover of each plot, as well as the percent cover of individual species. Carbon dioxide (CO2) exchange was also measured within each plot throughout the growing season. The date of senescence occurred significantly earlier in plots which had not been manipulated in any way, compared to all other treatments for all species. Salix arctica showed the greatest increase in cover over time at the species level. Lichen cover increased significantly in the deepened snow plots, and in general there were significant increases in percent cover in some functional groups over time. During June and into July the net CO2 flux was to the atmosphere. It was not until July 27 that these ecosystems became net carbon sinks. However, warming alone resulted in the ecosystem acting as a significant net carbon sink for the entire growing season. Plots exposed to warming alone were estimated to have removed approximately 19.94 g C m-2 from the atmosphere, whereas all other treatments were very similar to one another and estimated to have added approximately 3.12 g C m-2 to the atmosphere. Active layer depth and soil temperatures suggest that plots within the ambient snow zone may be receiving some additional snow due to their proximity to the fences. CO2 fluxes measured within the outer control plots suggest that the effect of warming alone could lead to this ecosystem being an even stronger net C sink under truly ambient snow conditions.