O Flúor nas águas subterrâneas da Formação Santa Maria, na região de Santa Cruz do Sul e Venâncio Aires, RS, Brasil

Concentrações de fluoreto, até 11 mg L-1, têm sido detectadas no aqüífero da Formação Santa Maria, Sistema Aqüífero Guarani, na região central do Estado do Rio Grande do Sul, Brasil. As águas subterrâneas são utilizadas para abastecimento público e levaram a incidência de fluorose dental. Esta área é a maior produtora de tabaco no Brasil, com uma longa história de aplicação de fertilizantes. Duas hipóteses são investigadas na busca de elucidar a origem das anomalias de flúor, a primeira, se refere à contaminação desencadeada pelo uso de fertilizantes fosfatados nas plantações de tabaco, e a segunda relacionada à causa natural, geogênica, determinada por interações rocha-fluido no aqüífero confinado. A hipótese antropogênica é testada nos solos e águas do aqüífero freático. O monitoramento das concentrações de fluoreto em campo, durante o período de dois anos e os experimentos de laboratório de retenção e lixiviação de fluoreto em solos mostram que este elemento é retido pelos solos argilosos. Durante o período de monitoramento das águas do freático, concentrações maiores que 0,16 mg L-1 não foram registradas, confirmando que o fluoreto é sorvido pelos solos argilosos ou transportado para horizontes mais permeáveis para superfície ou água subterrânea, em processo lento. Análise estatística fatorial e de agrupamento associada a técnicas de modelamento geoquímico apontaram as afinidades químicas das águas subterrâneas e possíveis origens para o fluoreto. As águas subterrâneas da Formação Santa Maria são distinguidas em quatro grupos hidroquímicos, o primeiro, cloretado sódico, evolui para águas subterrâneas bicarbonatadas sódicas, e ambos grupos apresentam anomalias de fluoreto. O terceiro grupo é de águas bicarbonatadas cálcicas e o quarto tem o magnésio como parâmetro característico. As hipóteses lançadas são também discutidas com base em análises de isótopos estáveis da água, e de sulfato e nitrato dissolvidos nas águas subterrâneas, que indicam a precipitação meteórica como a principal fonte de recarga do aqüífero estudado. A composição isotópica do sulfato da água é similar a do sulfato nos fertilizantes, não sendo possível distingui-lo de outras fontes que também apresentam assinaturas isotópicas semelhantes. No entanto, os isótopos do nitrato dissolvido nas águas subterrâneas e dos fertilizantes sugerem que não existe uma relação direta com o uso de fertilizantes NPK. Os dados apresentados indicam que os elevados teores de fluoreto no aqüífero da Formação Santa Maria resultam de interações água-rocha associadas à circulação de fluídos em estruturas tectônicas regionais levando a mistura de águas subterrâneas de aqüíferos mais profundos.

Caracterização da geometria de depósitos sedimentares na borda sudoeste da Bacia`Potiguar

This study focuses on the potential of several techniques used to identify depositional geometries and paleogeographical investigation on the SW border of the Potiguar Basin. Three areas were selected for an integrated geological, geophysical and geochemistry study. The main used techniques were facies analysis, remote sensing,ground penetrating radar (GPR) and gamma-ray in outcrops, as well as petrographic microscope observations and the using of scanning eletronic microscopic (SEM), and Carbon and Oxygen Isotopic study in the carbonate tufa. These methodological approaches were very efficient in the facies analysis of 2D geometries. The GPR profiles carried out in Quixeré identified important geological reflectors which allowed to the identification of depositional geometries of tufa. However, GPR profiles were not able to identify geological reflectors in the Apodi and Olho d´Água da Bica outcrops. Gammaray profiles also presented good results, which justify their use in 1D and 2D geometric analysis. Carbon and Oxygen Isotopic analyses were also used to investigate paleoenvironmental setting of tufa deposits. It is important to remark the excellent resultsof GRP using in the identification of deposition al geometries of tufa and their contact relationships with the underlying rocks. Field analysis of faults indicate a vertical sigma-1 orientation which was associated to normal faults

Rastreabilidade de farinha de vísceras de aves por isótopos estáveis em penas de frangos de corte

O objetivo deste trabalho foi detectar a presença de farinha de vísceras de aves (FV) na alimentação de frangos de corte, por meio da análise de penas por isótopos estáveis de carbono (13C/12C) e nitrogênio (15N/14N) e espectrofotometria de massas. Setecentos e vinte pintos machos Cobb foram submetidos aos seguintes tratamentos: ração vegetal (RV) à base de milho e farelo de soja, de 1 a 42 dias de idade; ração com 8% de farinha de vísceras de frango (FV), de 1 a 42 dias de idade; ração vegetal de 1 a 21 dias e ração com FV de 22 a 42 dias; ração vegetal de 1 a 35 dias e ração com FV de 36 a 42 dias; ração com FV de 1 a 21 dias, e ração vegetal de 22 a 42 dias; ração com FV de 1 a 35 dias, e ração vegetal de 36 a 42 dias. Foram colhidas amostras de penas de quatro aves por tratamento aos 7, 14, 21, 28, 35 e 42 dias de idade, que foram submetidas à análise isotópica de carbono (13C/12C) e nitrogênio (15N/14N) por espectrometria de massas. A aplicação da técnica dos isótopos estáveis de C e N a penas de frango de corte, após a utilização de farinha de vísceras na alimentação de frangos de corte, pelo período de até 21 dias, permite detectar a inclusão da FV na dieta.

Rastreabilidade da farinha de carne e ossos bovinos em ovos de poedeiras alimentadas com ingredientes alternativos

O objetivo deste trabalho foi detectar traços de farinha de carne e ossos bovinos, em ovos de poedeiras alimentadas com dietas comerciais com inclusão de ingredientes vegetais alternativos e leveduras. A detecção foi feita pela técnica dos isótopos estáveis do carbono e do nitrogênio. Foram utilizadas 384 poedeiras, distribuídas aleatoriamente em oito tratamentos. Os tratamentos consistiram de uma dieta-controle - à base de milho e farelo de soja - e sete dietas com inclusão de farinha de carne e ossos bovinos, acrescidas ou não de outros ingredientes (farelo de trigo, quirera de arroz, farelo de algodão, glúten de milho, levedura de cana e levedura de cerveja). No 35º dia, foram tomados aleatoriamente 24 ovos por tratamento: 12 para análise de ovo e 12 para análise de gema e albúmen, em separado. Após análise isotópica de carbono e nitrogênio, os resultados foram submetidos à análise multivariada de variância. As médias dos pares isotópicos dos ovos, gema e albúmen, em todos os tratamentos, diferiram daquelas do tratamento-controle. A técnica dos isótopos estáveis permite detectar, nos ovos, gema e albúmen, a farinha de carne e ossos bovinos utilizada na dieta de poedeiras, mesmo com a inclusão de outros ingredientes vegetais e leveduras.

Evaluation of two genera of benthic foraminifera for down-core paleotemperature studies in the western south atlantic

Neste estudo, compara-se a composição de isótopos de oxigênio de dois gêneros de foraminíferos bentônicos (Uvigerina e Cibicidoides) de amostras de topo de testemunho com a composição isotópica moderna da água do mar (d¹⁸O). Baseados em uma nova relação entre d¹⁸O e salinidade para a latitude média do Atlântico Sul ocidental, estimou-se a composição isotópica da calcita em equilíbrio (d¹⁸O_eq) a partir de duas equações diferentes: (1) O'Neil et al. (1969), modificada por McCorkle et al. (1997) e (2) Kim & O'Neil (1997). Utilizando (1), a pequena diferença entre d¹⁸O_eq e d¹⁸O de Uvigerina sugere que este gênero precipita as suas testas próximo ao equilíbrio com a água. Já os dados de d¹⁸O de Cibicidoides são 0,82 ‰ menores que a composição isotópica prevista. Ao contrário, utilizando (2) os dados de d¹⁸O de Cibicidoides mostram uma concordância excelente com a composição isotópica esperada, enquanto que os dados de d¹⁸O de Uvigerina são 0,69 ‰ maiores que os valores de equilíbrio previstos. A partir das evidências apresentadas neste trabalho e em estudos prévios sugerimos a utilização do gênero Cibicidoides e a aplicação da equação de Kim & O'Neil (1997) para pesquisas de paleotemperatura. Na ausência de suficientes espécimens de Cibicidoides sugerimos a utilização de d¹⁸O deUvigerina aplicando um fator de correção de -0,69 ‰.

Mineralogia e química dos sedimentos de fundo do médio e baixo Madeira e de seus principais tributários: Amazonas - Brasil

O objetivo deste trabalho foi investigar a variação composicional e mineralógica dos sedimentos de fundo transportados pelos grandes rios da Amazônia, com nascentes nos Andes provenientes de rochas cratônicas. O estudo foi realizado com base em análises granulométricas, mineralógicas e químicas, incluindo isótopos de Pb, no rio Madeira e seus tributários. O conjunto de dados mostra que os sedimentos de fundo do rio Madeira são granulométrica, mineralógica e quimicamente distintos de seus principais tributários. Os sedimentos do rio Madeira são mais arenosos; têm maior quantidade de quartzo; menor conteúdo de Al₂O₃, Fe₂O₃, K₂O, MgO, Na₂O, PF, Ti₂O, P₂O₅, MnO e elementos-traço; e a maturidade aumenta de montante para jusante. Dentre os sedimentos dos tributários, os rios Machado e secundariamente o rio Marmelos são os que mais se aproximam da composição dos sedimentos do rio Madeira, enquanto que o Jamari com sedimentos com maior proporção de TiO_2,, Zr, Y, Nb, Ga, Hf, U, Ta e ETR, é o mais distinto. Apesar do ambiente de intenso intemperismo e erosão, em condições tropicais úmidas, a que estão submetidas as rochas drenadas pela bacia do rio Madeira, os valores das razões Th/Co, Th/Pb, Th/Yb, Al/Pb, Zr/Co e a composição isotópica de Pb indicam fontes distintas para os sedimentos de fundo estudados. Os sedimentos do rio Madeira são provenientes, principalmente, de rochas máficas, enquanto que os sedimentos de seus tributários têm como fonte essencialmente rochas félsicas.

Detecção de farinha de vísceras nas fases de alimentação em frangos de corte pela técnica de isótopos estáveis

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

Datação por luminescência opticamente estimulada (Loe) e por 14C de materiais de uma topossequência de solos na Serra do Espinhaço Meridional - MG: abordagem paleoambiental

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

Estudos de inclusões fluidas e de isótopos estáveis no depósito Moreira Gomes do campo mineralizado do Cuiú-Cuiú, Província Aurífera do Tapajós, Estado do Pará

Moreira Gomes é um dos depósitos do campo mineralizado do Cuiú-Cuiú, província Aurífera do Tapajós, com recursos de 21,7 t de ouro. A zona mineralizada, com 1200 metros de comprimento, 30-50 metros de largura e, pelo menos, 400 metros de profundidade é controlada por uma estrutura subvertical de orientação E-W, associada a um sistema de falhas transcorrentes sinistrais. As rochas hospedeiras nesse depósito são predominantemente tonalitos de 1997 ± 2 Ma (Suite Intrusiva Creporizão). O estilo da alteração hidrotermal relacionado à mineralização é predominantemente fissural e localmente pervasivo. Os tipos de alteração hidrotermal são sericitização, carbonatação, cloritização, sulfetação, silicificação e epidotização, além da formação de veios de quartzo de espessuras variadas. Pirita é principal sulfeto e contém inclusões de galena, esfalerita, calcopirita e, em menor quantidade, de hessita e bismutinita. O ouro ocorre mais comumente como inclusão em cristais de pirita e, secundariamente, na forma livre em veios de quartzo. Ag, Pb e Bi foram detectados por análise semi-quantitativa como componentes das partículas de ouro. Estudo de inclusões fluidas identificou fluidos compostos por CO₂ (Tipo 1), H₂O-C O₂-sal (Tipo 2) e H₂O-sal (Tipo 3). O volátil CO₂ é predominante na fase carbônica. O fluido do Tipo 2 apresenta densidade baixa a moderada, salinidade entre 1,6 e 11,8 % em peso equivalente de NaCl e foi aprisionado principalmente entre 280° e 350°C. No fluido do Tipo 3 o sistema químico pode conter aCl₂ e, talvez, MgCl₂, e a salinidade varia de zero a 10,1% em peso equivalente de NaCl. Apenas localmente a salinidade atingiu 25% em peso equivalente de NaCl. Esse fluido foi aprisionado principalmente entre 120° e 220°C e foi interpretado como resultado de mistura de fluido aquoso mais quente e levemente mais salino, com fluido mais frio e diluído. Globalmente, o estudo das inclusões fluidas indica estado heterogêneo durante o aprisionamento e ocorrência de separação de fases, mistura, flutuação de pressão e reequilíbrio das inclusões durante aprisionamento. A composição isotópica do fluido em equilíbrio com minerais hidrotermais (quartzo, clorita e calcita e pirita) e de inclusões fluidas apresenta valores de δ¹⁸O e δD entre +0,5 e +9,8 ‰, e -49 a -8 ‰, respectivamente. Os valores de ³⁴S de pirita (-0,29 ‰ a 3,95 ‰) são provavelmente indicativos da presença de enxofre magmático. Pares minerais forneceram temperaturas de equilíbrio isotópico em geral concordante com as temperaturas de homogeneização de inclusões fluidas e compatíveis com as relações texturais. Os resultados isotópicos, combinados com os dados mineralógicos e de inclusões fluidas são interpretados como produto da evolução de um sistema magmático hidrotermal em três estágios. (1) Exsolução de fluido magmático aquoso e portador de CO₂ entre 400°C e 320-350°C, seguido de separação de fases e precipitação principal da assembleia clorita-sericita-pirita-quartzo-ouro sob pressões menores que 2,1 kb e a 6-7 km de profundidade. (2) Resfriamento e continuação da exsolução do CO₂ do fluido magmático geraram fluido aquoso, mais pobre a desprovido de CO₂ e levemente mais salino, com aprisionamento dominantemente a 250°-280°C. A assembleia hidrotermal principal ainda precipitou, mas epidoto foi a principal fase nesse estágio. (3) Mistura do fluido aquoso do estágio 2, mais quente e mais salino, com um fluido aquoso mais frio e menos salino, de origem meteórica. Carbonatação está associada com esse estágio. A assembleia hidrotermal e os valores isotópicos indicam que fluido foi neutro a levemente alcalino e relativamente reduzido, que H₂S (ou HS-) pode ter sido a espécie de enxofre predominante, e que Au(HS) ^-2 deve ter sido o complexo transportador de ouro. A deposição do ouro em Moreira Gomes ocorreu em resposta a diversos mecanismos, envolvendo a separação de fases, mistura e reações fluido-rocha. O depósito Moreira Gomes é interpretado como o produto de um sistema magmático-hidrotermal, mas não possui feições clássicas de depósitos relacionados a intrusões graníticas, tanto oxidadas como reduzidas. A idade de deposição do minério (1,86 Ga) sugere que o sistema magmático-hidrotermal pode estar relacionado com a fase final do extenso magmatismo cálcio-alcalino da Suíte Intrusiva Parauari, embora o magmatismo transicional a alcalino da Suíte Intrusiva Maloquinha não possa ser descartado.

Estudos isotópicos e de inclusões fluidas no depósito central do campo mineralizado do Cuiú-Cuiú, província aurífera do Tapajós, estado do Pará

Central é um depósito aurífero do campo mineralizado do Cuiú-Cuiú, Província Aurífera do Tapajós, Cráton Amazônico. A zona mineralizada está hospedada em falha e compreende 800m de comprimento na direção NW-SE, seguindo o trend regional da província Tapajós, com largura entre 50 e 70m e profundidade vertical de pelo menos 450m. A mineralização está hospedada em monzogranito datado em 1984±3 Ma e atribuído à Suíte Intrusiva Parauari. Os recursos auríferos preliminarmente definidos são de 18,6t de ouro. A alteração hidrotermal é predominantemente fissural. Sericitização, cloritização, silicificação, carbonatação e sulfetação foram os tipos de alteração identificados. Pirita é o sulfeto principal e os demais sulfetos (calcopirita, esfalerita e galena) estão em fraturas ou nas bordas da pirita. O ouro preenche fraturas da pirita e análises semi-quantitativas detectaram Ag associada ao ouro. Foram identificados três tipos de inclusões fluidas hospedados em veios e vênulas de quartzo. O tipo 1 é o menos abundante e consiste em inclusões fluidas compostas por uma (CO_2vapor) ou duas fases (CO_2liq-CO_2vapor), o tipo 2 tem abundância intermediária e é formado por inclusões fluidas compostas por duas (H₂O_liq-CO_2liq) ou três fases (H₂O_liq-CO_2liq-CO_2vapor) e o tipo 3 é o mais abundante e consiste em inclusões fluidas compostas por duas fases (H₂O_liq- H₂Ovapor). O CO₂ representa o volátil nas inclusões com CO₂ e essas (tipo 1 e 2) foram geradas pelo processo de separação de fases oriundo de um fluido aquo-carbônico. A densidade global (0,33 - 0,80 g/cm³) e a salinidade (11,15 - 2,42 % em peso equivalente de NaCl) desse fluido são baixas a moderadas e a temperatura de homogeneização mostra um máximo em 340ºC. Quanto ao tipo 3, o NaCl é o principal sal, a densidade global está no intervalo de 0,65 a 1,11 g/cm³, a salinidade compreendida entre 1,16 e 13,3 % em peso equivalente de NaCl e a temperatura de homogeneização é bimodal, com picos em 120-140ºC e 180ºC. A composição isotópica das inclusões fluidas presentes no quartzo e do quartzo, calcita e clorita mostram valores de δ¹⁸O e δD de +7,8 a +13,6 ‰ e -15 a -35 ‰, respectivamente. Os valores de δ³⁴S na pirita são de +0,5 a +4,0 ‰ e δ¹³C na calcita e CO₂ de inclusões fluidas de -18 a -3,7 ‰. Os valores de δ¹⁸O_H2O e de δD_H2O no quartzo e inclusões fluidas, respectivamente, plotam no campo das águas metamórficas, com um desvio em direção à linha da água meteórica. Considerando a inexistência de evento metamórfico na região do Tapajós à época da mineralização, o sistema hidrotermal responsável pela mineralização no Central, inicialmente, deu-se a partir de fluidos aquo-carbônicos magmático-hidrotermais, exsolvidos por magma félsico relacionado com a fase mais tardia de evolução da Suíte Intrusiva Parauari. As inclusões aquo-carbônicas e carbônicas formaram-se nessa etapa, predominantemente em torno de 340°C. A contínua exsolução de fluido pelo magma levou ao empobrecimento em CO₂ nas fases mais tardias e, com o resfriamento do fluido, as inclusões aquosas passaram a predominar. A partir daí o sistema pode ter interagido com água meteórica, responsável pelo aprisionamento da maior parte das inclusões aquosas de mais baixa temperatura. É possível que parte das inclusões aquosas (as de maior temperatura) represente a mistura local dos fluidos de origens distintas. Essas observações e interpretações permitem classificar Central como um depósito de ouro magmático-hidrotermal relacionado à fase final da formação da Suíte Intrusiva Parauari.

Carbon isotope fractionation for cotton genotype selection

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

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.

Conductancia del mesófilo en especies forestales: implicaciones en la respuesta fotosintética a la disponibilidad de agua y luz

La fotosíntesis es el proceso biológico que permite la producción primaria y, por tanto, la vida en nuestro planeta. La tasa fotosintética viene determinada por la ‘maquinaria’ bioquímica y las resistencias difusivas al paso del CO2 desde la atmósfera hasta su fijación en el interior de los cloroplastos. Históricamente la mayor resistencia difusiva se ha atribuido al cierre estomático, sin embargo ahora sabemos, debido a las mejoras en las técnicas experimentales, que existe también una resistencia grande que se opone a la difusión del CO2 desde los espacios intercelulares a los lugares de carboxilación. Esta resistencia, llamada normalmente por su inversa: la conductancia del mesófilo (gm), puede ser igual o incluso superior a la resistencia debida por el cierre estomático. En la presente tesis doctoral he caracterizado la limitación que ejerce la resistencia del mesófilo a la fijación de CO2 en diversas especies forestales y en distintos momentos de su ciclo biológico. En la fase de regenerado, hemos estudiado tres situaciones ambientales relevantes en el mayor éxito de su supervivencia, que son: el déficit hídrico, su interacción con la irradiancia y el paso del crecimiento en la sombra a mayor irradiancia, como puede suceder tras la apertura de un hueco en el dosel forestal. En la fase de arbolado adulto se ha caracterizado el estado hídrico y el intercambio gaseoso en hojas desarrolladas a distinta irradiancia dentro del dosel vegetal durante tres años contrastados en pluviometría. Para cada tipo de estudio se han empleado las técnicas ecofisiológicas más pertinentes para evaluar el estado hídrico y el intercambio gaseoso. Por su complejidad y la falta de un método que permita su cuantificación directa, la gm ha sido evaluada por los métodos más usados, que son: la discriminación isotópica del carbono 13, el método de la J variable, el método de la J constante y el método de la curvatura. Los resultados más significativos permiten concluir que la limitación relativa a la fotosíntesis por la conductancia estomática, del mesófilo y bioquímica es dependiente de la localización de la hoja en el dosel forestal. Por primera vez se ha documentado que bajo estrés hídrico las hojas desarrolladas a la sombra estuvieron más limitadas por una reducción en la gm, mientras que las hojas desarrolladas a pleno sol estuvieron más limitadas por reducción mayor de la conductancia estomática (gsw). Encontramos buena conexión entre el aparato fotosintético foliar y el sistema hídrico debido al alto grado de correlación entre la conductancia hidráulica foliar aparente y la concentración de CO2 en los cloroplastos en distintas especies forestales. Además, hemos mostrado diferentes pautas de regulación del intercambio gaseoso según las particularidades ecológicas de las especies estudiadas. Tanto en brinzales crecidos de forma natural y en el arbolado adulto como en plántulas cultivadas en el invernadero la ontogenia afectó a las limitaciones de la fotosíntesis producidas por estrés hídrico, resultando que las limitaciones estomáticas fueron dominantes en hojas más jóvenes mientras que las no estomáticas en hojas más maduras. La puesta en luz supuso un gran descenso en la gm durante los días siguientes a la transferencia, siendo este efecto mayor según el grado de sombreo previo en el que se han desarrollado las hojas. La aclimatación de las hojas a la alta irradiancia estuvo ligada a las modificaciones anatómicas foliares y al estado de desarrollo de la hoja. El ratio entre la gm/gsw determinó la mayor eficiencia en el uso del agua y un menor estado oxidativo durante la fase de estrés hídrico y su posterior rehidratación, lo cual sugiere el uso de este ratio en los programas de mejora genética frente al estrés hídrico. Debido a que la mayoría de modelos de estimación de la producción primaria bruta (GPP) de un ecosistema no incluye la gm, los mismos están incurriendo en una sobreestimación del GPP particularmente bajo condiciones de estrés hídrico, porque más de la mitad de la reducción en fotosíntesis en hojas desarrolladas a la sombra se debe a la reducción en gm. Finalmente se presenta un análisis de la importancia en las estimas de la gm bajo estrés hídrico de la refijación del CO2 emitido en la mitocondria a consecuencia de la fotorrespiración y la respiración mitocondrial en luz. ABSTRACT Photosynthesis is the biological process that supports primary production and, therefore, life on our planet. Rates of photosynthesis are determined by biochemical “machinery” and the diffusive resistance to the transfer of CO2 from the atmosphere to the place of fixation within the chloroplasts. Historically the largest diffusive resistance was attributed to the stomata, although we now know via improvements in experimental techniques that there is also a large resistance from sub-stomatal cavities to sites of carboxylation. This resistance, commonly quantified as mesophyll conductance (gm), can be as large or even larger than that due to stomatal resistance. In the present PhD I have characterized the limitation exerted by the mesophyll resistance to CO2 fixation in different forest species at different stages of their life cycle. In seedlings, we studied three environmental conditions that affect plant fitness, namely, water deficit, the interaction of water deficit with irradiance, and the transfer of plants grown in the shade to higher irradiance as can occur when a gap opens in the forest canopy. At the stage of mature trees we characterized water status and gas exchange in leaves developed at different irradiance within the canopy over the course of three years that had contrasting rainfall. For each study we used the most relevant ecophysiological techniques to quantify water relations and gas exchange. Due to its complexity and the lack of a method that allows direct quantification, gm was estimated by the most commonly used methods which are: carbon isotope discrimination, the J-variable, constant J and the curvature method The most significant results suggest that the relative limitation of photosynthesis by stomata, mesophyll and biochemistry depending on the position of the leaf within the canopy. For the first time it was documented that under water stress shaded leaves were more limited by a reduction in gm, while the sun-adapted leaves were more limited by stomatal conductance (gsw). The connection between leaf photosynthetic apparatus and the hydraulic system was shown by the good correlations found between the apparent leaf hydraulic conductance and the CO2 concentration in the chloroplasts in shade- and sun-adapted leaves of several tree species. In addition, we have revealed different patterns of gas exchange regulation according to the functional ecology of the species studied. In field grown trees and greenhouse-grown seedlings ontogeny affected limitations of photosynthesis due to water stress with stomatal limitations dominating in young leaves and nonstomatal limitations in older leaves. The transfer to high light resulted in major decrease of gm during the days following the transfer and this effect was greater as higher was the shade which leaves were developed. Acclimation to high light was linked to the leaf anatomical changes and the state of leaf development. The ratio between the gm/gsw determined the greater efficiency in water use and reduced the oxidative stress during the water stress and subsequent rehydration, suggesting the use of this ratio in breeding programs aiming to increase avoidance of water stress. Because most models to estimate gross primary production (GPP) of an ecosystem do not include gm, they are incurring an overestimation of GPP particularly under conditions of water stress because more than half of An decrease in shade-developed leaves may be due to reduction in gm. Finally, we present an analysis of the importance of how estimates of gm under water stress are affected by the refixation of CO2 that is emitted from mitochondria via photorespiration and mitochondrial respiration in light.

Reconstituição da paleoprecipitação no sul do Nordeste Brasileiro durante os dois últimos ciclos glaciais a partir da aplicação de registros isotópicos de oxigênio de estalagmites da Chapada Diamantina, Bahia.

Registros isotópicos de oxigênio obtidos em alta resolução das estalagmites CL2 e MAG das cavernas Calixto e Marota, região da Chapada Diamantina (CD) (12ºS), Estado da Bahia, sul do Nordeste brasileiro (sNEB), permitiram reconstituir as mudanças passadas da precipitação entre 165-128 e 59-39 mil anos A.P. Para a reconstituição paleoclimática considerou-se resultados de um estudo de calibração realizado em duas cavernas da CD o qual demonstrou uma relação entre composição isotópica da água meteórica e de gotejamento e sugeriu um ambiente adequado para a deposição do espeleotema em condições equilíbrio e/ou próximas com a água de gotejamento. A interpretação da paleoprecipitação através dos registros isotópicos \'\'delta\' POT.18\'O das estalagmites também foi baseada na relação entre composição isotópica da água da precipitação e a quantidade de chuva obtidos em estações da IAEA-GNIP no Brasil e de simulações das variações do \'\'delta\' POT.18\'O da chuva através do modelo climático ECHAM-4. Esses dados indicaram o efeito quantidade (amount effect) como fator preponderante de controle isotópico da água da chuva que formam os espeleotemas na CD, significando que a diminuição dos valores de \'\'delta\' POT.18\'O está associada ao aumento do volume de chuvas e vice-versa. Os registros de \'\'delta\' POT.18\'O dos espeleotemas permitiram reconstituir a variação da paleoprecipitação na escala orbital e milenar durante o penúltimo glacial bem como correlacionar mudanças na paleoprecipitação no sNEB com eventos milenares registrados na Groelândia no último glacial. Os registros da CD indicaram um aumento (diminuição) da paleoprecipitação na Bahia relacionado a diminuição (aumento) da insolação austral de verão a 10ºS durante o penúltimo glacial, similar ao observado no último ciclo precessional. Na escala orbital os registros da CD estiveram em antifase com os paleoindicadores isotópicos do Sudeste brasileiro e em fase com os valores de\'\'delta\' POT.18\'O dos espeleotemas do leste da China. Esse padrão de precipitação é similar ao observado na última glaciação e sugere que a variação na insolação de verão afetou as monções sul-americanas (MSA) promovendo mudanças na precipitação no sNEB no penúltimo glacial. Condições áridas no sNEB durante o aumento da insolação de verão estariam provavelmente associadas ao aprofundamento da subsidência de ar provocado pelo fortalecimento da circulação leste-oeste da MSA devido ao aumento das atividades convectivas na Amazônia o que teria, favorecido um posicionamento mais a sul da Zona de Convergência do Atlântico Sul (ZCAS). O oposto também ocorreria durante as fases de baixa insolação de verão quando a MSA estaria provavelmente mais desintensificada. Durante o penúltimo glacial (Terminação Glacial II) abruptas oscilações nos registros da CD para valores mais baixos de \'\'delta\' POT.18\'O indicaram um profundo aumento da precipitação coincidente com o evento Heinrich (H11). Nesse período a paleoprecipitação no sNEB esteve correlacionada negativamente com as mudanças climáticas na China e no oeste amazônico (Peru) e positivamente com o Sudeste brasileiro. Interpretou-se que as anomalias positivas da precipitação no sNEB podem ter estado relacionadas ao deslocamento para sul da Zona de Convergência Intertropical (ZCIT) bem como com a intensificação da MSA e ZCAS nesse período. Finalmente, oscilações isotópicas abruptas para valores mais altos observadas durante o estágio marinho isotópico 3 coincidentes com os eventos quentes registrados na Groelândia, denominados de eventos Dansgaard-Oeschger (DO), foram interpretados como a ocorrência de eventos muito secos no sNEB. Essas variações da precipitação na escala milenar, que estão em fase com os registros no Peru, podem ter estado relacionadas ao deslocamento para norte da ZCIT o que teria promovido uma profunda desintensificação da MSA.