DesG: uma nova metodologia para intervenção e coleta de dados magnetotelúricos

Esta tese apresenta duas contribuições para a interpretação geofísica, voltadas ao Método Magnetotelúrico. A primeira trata de uma nova abordagem para a interpretação MT, denominada MÉTODO DESCRITIVO-GEOLÓGICO (DesG), em alusão à incorporação explícita de informação a priori de correlação fácil com a descrição geológica tradicional. O intérprete define por meio de elementos geométricos (pontos e linhas) o arcabouço de feições geológicas bem como fornece valores de resistividade aos corpos geológicos presumidos. O método estima a distribuição de resistividade subsuperficial em termos de fontes anômalas próximas aos elementos geométricos, ajustando as respostas produzidas por estes corpos às medidas de campo. A solução obtida fornece então informações que podem auxiliar na modificação de algumas informações a priori imprecisas, permitindo que sucessivas inversões sejam realizadas até que a solução ajuste os dados e faça sentido geológico. Entre as características relevantes do método destacam-se: (i) os corpos podem apresentar resistividade maior ou menor do que a resistividade do meio encaixante, (ii) vários meios encaixantes contendo ou não corpos anômalos podem ser cortados pelo perfil e (iii) o contraste de resistividade entre corpo e encaixante pode ser abrupto ou gradativo. A aplicação do método a dados sintéticos evidencia, entre outras vantagens, a sua potencialidade para estimar o mergulho de falhas com inclinação variável, que merece especial interesse em Tectônica, e delinear soleiras de diabásio em bacias sedimentares, um sério problema para a prospecção de petróleo. O método permite ainda a interpretação conjunta das causas do efeito estático e das fontes de interesse. A aplicação a dados reais é ilustrada tomando-se como exemplo dados do COPROD2, cuja inversão produziu soluções compatíveis com o conhecimento geológico sobre a área. A segunda contribuição refere-se a desenho de experimento geofísico. Por meio de indicadores diversos, em especial a matriz densidade de informação, é mostrado que a resolução teórica dos dados pode ser estudada, que guia o planejamento da prospecção. A otimização do levantamento permite determinar os períodos e as posições das estações de medida mais adequados ao delineamento mais preciso de corpos cujas localizações são conhecidas aproximadamente.

Precipitação sazonal sobre a Amazônia oriental no período chuvoso: observações e simulações regionais com o RegCM3

O presente trabalho apresenta uma contribuição aos estudos de modelagem climática com ênfase na variabilidade pluviométrica sazonal da Amazônia oriental, durante as estações de verão e outono (DJF e MAM). Baseado nos resultados das simulações regionais do RegCM3 para um período de 26 anos (1982/83 a 2007/08) e usando domínio em alta resolução espacial (30 Km) e dois diferentes esquemas de convecção (Grell e MIT), foi investigado o desempenho do modelo em simular a distribuição regional de precipitação sazonal na Amazônia oriental, com referência a um novo conjunto de dados observacional compilado com informações de uma ampla rede integrada de estações pluviométricas. As análises quantitativas evidenciaram que o RegCM3 apresenta erros sistemáticos, sobretudo aqueles relacionados com viés seco no Amapá e norte/nordeste do Pará usando ambos os esquemas Grell e MIT, os quais apontam que o modelo não reproduz as características da ZCIT sobre o Atlântico equatorial. As simulações usando MIT, também apresentaram viés úmido no sudoeste/sul/sudeste do Pará e norte do Tocantins. Além disso, através da técnica de composições, também foi investigado o desempenho do RegCM3 em reproduzir os padrões espaciais anômalos de precipitação sazonal em associação aos episódios ENOS, e as fases do gradiente térmico sobre o Atlântico intertropical. Os resultados demonstraram que o modelo conseguiu representar realisticamente bem o padrão espacial das anomalias pluviométricas acima (abaixo) do normal em grande parte da Amazônia oriental, durante os conhecidos cenários favoráveis, i.e., condições de La Niña e gradiente de aTSM para o Atlântico sul (desfavoráveis, i.e., El Niño e gradiente de aTSM para o Atlântico norte.

Assessment of surface water in two Amazonian rivers impacted by industrial wastewater, Barcarena City, Pará State (Brazil)

No ano de 2007, águas superficiais foram coletadas a partir de 21 pontos de amostragem na cidade de Barcarena, Região Norte, Brasil: um ponto de amostragem localizado em um pequeno córrego que recebe descarga de resíduos a partir da indústria de beneficiamento do caulim e deságua no Rio Curuperê, três pontos de amostragens localizados próximos de fontes que emergem na margem esquerda e deságuam no Rio Curuperê, nove pontos de amostragem no Rio Curuperê, que deságua no Rio Dendê, e oito no Rio Dendê, afluente da margem esquerda do Rio Pará. Para todas as amostras de água foram quantificadas 14 variáveis físico-químicas e níveis de 12 metais. Os resultados nos pontos próximos das fontes do Rio Curuperê apresentaram perfil físico-químico e níveis de metais típicos para águas superficiais e esses valores foram utilizados como referência para comparar e identificar possíveis alterações nas características químicas para os demais pontos de amostragem. Quando os resultados das fontes do rio Curuperê foram comparados com os resultados do ponto próximo a descarga de resíduos industriais foram observadas fortes alterações nos valores de 6 variáveis físico-químicas (pH, condutividade elétrica (EC), total de sólidos suspensos (TDS), nitrogênio amoniacal (N-NH4), sulfato (SO4) e salinidade) e aumento em magnitude dos níveis de quatro metais (Al, Fe, Mn e Zn), caracterizando que esses resíduos eram descarregados no ambiente sem tratamentos adequados. Os resultados nos demais pontos de amostragem demonstraram que estas condições anômalas também foram encontradas ao longo dos Rios Curuperê e Dendê, principalmente durante a maré baixa. A partir desta caracterização química das águas foram identificadas condições prejudiciais aos ecossistemas aquáticos e potencial risco à saúde da população local que usa os rios para consumo, recreação e transporte.

Análise e modelagem hidrometeorológica na Bacia do Rio Tocantins em Marabá-PA

Esta pesquisa objetivou desenvolver um modelo estatístico de previsão de vazão para Marabá - PA, bem como avaliar a estrutura dinâmica atmosférica associada aos extremos do regime hidrológico da bacia do rio Tocantins. O modelo hidrológico de regressão linear múltipla utilizou as séries de observações fluviométricas e pluviométricas obtidas no banco de dados da ANA. Os testes de validação do modelo estatístico com coeficiente de Nash acima de 0,9 e erro padrão de 1,5 % e 5% nos períodos de cheia e estiagem, respectivamente, permitem que as previsões de vazão em Marabá possam ser geradas com antecedência de 2 a 4 (3 a 5) dias para o período da cheia (estiagem). Através da técnica de composições considerando todos os anos com registro de vazão acima/muito acima e abaixo/muito abaixo do normal, obtidos pela metodologia dos percentis, investigaram-se as características regionais da precipitação e a estrutura dinâmica atmosférica em cada mês (Novembro a Abril). As composições dos anos com vazão acima/muito acima mostraram que a precipitação sobre a bacia foi acima do normal em todos os meses, sendo que os padrões de grande escala indicaram a configuração associada ao fenômeno La Niña no Pacífico e condições de resfriamento no Atlântico Sul; intensificação tanto do ramo ascendente zonal da célula de Walker como do ramo ascendente meridional da célula de Hadley; intensificação da Alta da Bolívia posicionada mais a leste e anomalias negativas de ROL associadas à atuação conjunta da ZCAS e ZCIT. Inversamente, as composições dos anos com vazão abaixo/muito abaixo evidenciaram a predominância de precipitação abaixo do normal em toda bacia hidrográfica, a qual se associou com as condições de aquecimento (El Niño) sobre o Pacífico, Atlântico sul aquecido, célula de Walker e Hadley com enfraquecimento dos movimentos ascendentes, posicionamento da Alta da Bolívia mais a oeste com anomalias positivas de ROL indicando inibição da atividade convectiva tropical. Adicionalmente, uma análise quantitativa dos impactos sócio-econômicos sobre os principais núcleos da cidade de Marabá revelou que aproximadamente 10 mil pessoas (5% da população) são atingidas pela cheia do rio Tocantins com custos nas operações de enchente acima de R$ 500.000,00, considerando o caso de 2005.

Metodologia geofísica aplicada ao estudo arqueológico dos sítios Bittencourt e Jambuaçu, Estado do Pará

Medidas magnéticas, de radiação gama e eletromagnéticas com o Radar de Penetração no Solo (GPR) foram realizadas em dois sítios arqueológicos localizados nas proximidades de uma área destinada a construção de um mineroduto. O objetivo das medidas foi identificar locais para escavação arqueológica visando a descoberta de remanescentes de habitantes que povoaram a região em tempos pretéritos, auxiliando nos estudos para recompor a sua história. As medidas magnéticas evidenciaram anomalias que levaram a descoberta de uma grande quantidade de fragmentos cerâmicos e de outros artefatos. As medidas de radiação gama permitiram a discriminação entre solo arenoso e argiloso, o que representa um resultado importante, pois o solo arenoso geralmente é correlacionado a camada de ocupação arqueológica onde são encontrados os remanescentes. O Radar de Penetração no Solo confirmou as anomalias magnéticas e mostrou outras feições anômalas não relacionadas aos remanescentes arqueológicos.

Mapeamento do relevo do embasamento de bacias sedimentares através da inversão gravimétrica vinculada

Apresentamos três novos métodos estáveis de inversão gravimétrica para estimar o relevo de uma interface arbitrária separando dois meios. Para a garantia da estabilidade da solução, introduzimos informações a priori sobre a interface a ser mapeada, através da minimização de um (ou mais) funcional estabilizante. Portanto, estes três métodos se diferenciam pelos tipos de informação físico-geológica incorporados. No primeiro método, denominado suavidade global, as profundidades da interface são estimadas em pontos discretos, presumindo-se o conhecimento a priori sobre o contraste de densidade entre os meios. Para a estabilização do problema inverso introduzimos dois vínculos: (a) proximidade entre as profundidades estimadas e verdadeiras da interface em alguns pontos fornecidas por furos de sondagem; e (b) proximidade entre as profundidades estimadas em pontos adjacentes. A combinação destes dois vínculos impõe uma suavidade uniforme a toda interface estimada, minimizando, simultaneamente em alguns pontos, os desajustes entre as profundidades conhecidas pelas sondagens e as estimadas nos mesmos pontos. O segundo método, denominado suavidade ponderada, estima as profundidades da interface em pontos discretos, admitindo o conhecimento a priori do contraste de densidade. Neste método, incorpora-se a informação geológica que a interface é suave, exceto em regiões de descontinuidades produzidas por falhas, ou seja, a interface é predominantemente suave porém localmente descontínua. Para a incorporação desta informação, desenvolvemos um processo iterativo em que três tipos de vínculos são impostos aos parâmetros: (a) ponderação da proximidade entre as profundidades estimadas em pontos adjacentes; (b) limites inferior e superior para as profundidades; e (c) proximidade entre todas as profundidades estimadas e um valor numérico conhecido. Inicializando com a solução estimada pelo método da suavidade global, este segundo método, iterativamente, acentua as feições geométricas presentes na solução inicial; ou seja, regiões suaves da interface tendem a tornar-se mais suaves e regiões abruptas tendem a tornar-se mais abruptas. Para tanto, este método atribui diferentes pesos ao vínculo de proximidade entre as profundidades adjacentes. Estes pesos são automaticamente atualizados de modo a acentuar as descontinuidades sutilmente detectadas pela solução da suavidade global. Os vínculos (b) e (c) são usados para compensar a perda da estabilidade, devida à introdução de pesos próximos a zero em alguns dos vínculos de proximidade entre parâmetros adjacentes, e incorporar a informação a priori que a região mais profunda da interface apresenta-se plana e horizontal. O vínculo (b) impõe, de modo estrito, que qualquer profundidade estimada é não negativa e menor que o valor de máxima profundidade da interface conhecido a priori; o vínculo (c) impõe que todas as profundidades estimadas são próximas a um valor que deliberadamente viola a profundidade máxima da interface. O compromisso entre os vínculos conflitantes (b) e (c) resulta na tendenciosidade da solução final em acentuar descontinuidades verticais e apresentar uma estimativa suave e achatada da região mais profunda. O terceiro método, denominado mínimo momento de inércia, estima os contrastes de densidade de uma região da subsuperfície discretizada em volumes elementares prismáticos. Este método incorpora a informação geológica que a interface a ser mapeada delimita uma fonte anômala que apresenta dimensões horizontais maiores que sua maior dimensão vertical, com bordas mergulhando verticalmente ou em direção ao centro de massa e que toda a massa (ou deficiência de massa) anômala está concentrada, de modo compacto, em torno de um nível de referência. Conceitualmente, estas informações são introduzidas pela minimização do momento de inércia das fontes em relação ao nível de referência conhecido a priori. Esta minimização é efetuada em um subespaço de parâmetros consistindo de fontes compactas e apresentando bordas mergulhando verticalmente ou em direção ao centro de massa. Efetivamente, estas informações são introduzidas através de um processo iterativo inicializando com uma solução cujo momento de inércia é próximo a zero, acrescentando, em cada iteração, uma contribuição com mínimo momento de inércia em relação ao nível de referência, de modo que a nova estimativa obedeça a limites mínimo e máximo do contraste de densidade, e minimize, simultaneamente, os desajustes entre os dados gravimétricos observados e ajustados. Adicionalmente, o processo iterativo tende a "congelar" as estimativas em um dos limites (mínimo ou máximo). O resultado final é uma fonte anômala compactada em torno do nível de referência cuja distribuição de constraste de densidade tende ao limite superior (em valor absoluto) estabelecido a priori. Estes três métodos foram aplicados a dados sintéticos e reais produzidos pelo relevo do embasamento de bacias sedimentares. A suavidade global produziu uma boa reconstrução do arcabouço de bacias que violam a condição de suavidade, tanto em dados sintéticos como em dados da Bacia do Recôncavo. Este método, apresenta a menor resolução quando comparado com os outros dois métodos. A suavidade ponderada produziu uma melhoria na resolução de relevos de embasamentos que apresentam falhamentos com grandes rejeitos e altos ângulos de mergulho, indicando uma grande potencialidade na interpretação do arcabouço de bacias extensionais, como mostramos em testes com dados sintéticos e dados do Steptoe Valley, Nevada, EUA, e da Bacia do Recôncavo. No método do mínimo momento de inércia, tomou-se como nível de referência o nível médio do terreno. As aplicações a dados sintéticos e às anomalias Bouguer do Graben de San Jacinto, California, EUA, e da Bacia do Recôncavo mostraram que, em comparação com os métodos da suavidade global e ponderada, este método estima com excelente resolução falhamentos com pequenos rejeitos sem impor a restrição da interface apresentar poucas descontinuidades locais, como no método da suavidade ponderada.

Inversão de momentos de fonte em métodos potenciais

A inversão de momentos de fonte gravimétrica tridimensional é analisada em duas situações. Na primeira se admite conhecer apenas a anomalia. Na segunda se admite conhecer, além da anomalia, informação a priori sobre o corpo anômalo. Sem usar informação a priori, mostramos que é possível determinar univocamente todo momento, ou combinação linear de momentos, cujo núcleo polinomial seja função apenas das coordenadas Cartesianas que definem o plano de medida e que tenha Laplaciano nulo. Além disso, mostramos que nenhum momento cujo núcleo polinomial tenha Laplaciano não nulo pode ser determinado. Por outro lado, informação a priori é implicitamente introduzida se o método de inversão de momentos se baseia na aproximação da anomalia pela série truncada obtida de sua expansão em multipolos. Dado um centro de expansão qualquer, o truncamento da série impõe uma condição de regularização sobre as superfícies equipotenciais do corpo anômalo, que permite estimar univocamente os momentos e combinações lineares de momentos que são os coeficientes das funções-bases da expansão em multipolos. Assim, uma distribuição de massa equivalente à real é postulada, sendo o critério de equivalência especificado pela condição de ajuste entre os campos observado e calculado com a série truncada em momentos de uma ordem máxima pré-estabelecida. Os momentos da distribuição equivalente de massa foram identificados como a solução estacionária de um sistema de equações diferenciais lineares de 1a. ordem, para a qual se asseguram unicidade e estabilidade assintótica. Para a série retendo momentos até 2a. ordem, é implicitamente admitido que o corpo anômalo seja convexo e tenha volume finito, que ele esteja suficientemente distante do plano de medida e que a sua distribuição espacial de massa apresente três planos ortogonais de simetria. O método de inversão de momentos baseado na série truncada (IMT) é adaptado para o caso magnético. Para este caso, mostramos que, para assegurar unicidade e estabilidade assintótica, é suficiente pressupor, além da condição de regularização, a condição de que a magnetização total tenha direção e sentido constantes, embora desconhecidos. O método IMT baseado na série de 2a. ordem (IMT2) é aplicado a anomalias gravimétricas e magnéticas tridimensionais sintéticas. Mostramos que se a fonte satisfaz as condições exigidas, boas estimativas da sua massa ou vetor momento de dipolo anômalo total, da posição de seu centro de massa ou de momento de dipolo e das direções de seus três eixos principais são obtidas de maneira estável. O método IMT2 pode falhar parcialmente quando a fonte está próxima do plano de medida ou quando a anomalia tem efeitos localizados e fortes de um corpo pequeno e raso e se tenta estimar os parâmetros de um corpo grande e profundo. Definimos por falha parcial a situação em que algumas das estimativas obtidas podem não ser boas aproximações dos valores verdadeiros. Nas duas situações acima descritas, a profundidade do centro da fonte (maior) e as direções de seus eixos principais podem ser erroneamente estimadas, embora que a massa ou vetor momento de dipolo anômalo total e a projeção do centro desta fonte no plano de medida ainda sejam bem estimados. Se a direção de magnetização total não for constante, o método IMT2 pode fornecer estimativas erradas das direções dos eixos principais (mesmo se a fonte estiver distante do plano de medida), embora que os demais parâmetros sejam bem estimados. O método IMT2 pode falhar completamente se a fonte não tiver volume finito. Definimos por falha completa a situação em que qualquer estimativa obtida pode não ser boa aproximação do valor verdadeiro. O método IMT2 é aplicado a dados reais gravimétricos e magnéticos. No caso gravimétrico, utilizamos uma anomalia situada no estado da Bahia, que se supõe ser causada por um batólito de granito. Com base nos resultados, sugerimos que as massas graníticas geradoras desta anomalia tenham sido estiradas na direção NNW e adelgaçadas na direção vertical durante o evento compressivo que causou a orogênese do Sistema de Dobramentos do Espinhaço. Além disso, estimamos que a profundidade do centro de massa da fonte geradora é cerca de 20 km. No caso magnético, utilizamos a anomalia de um monte submarino situado no Golfo da Guiné. Com base nos resultados, estimamos que o paleopolo magnético do monte submarino tem latitude 50°48'S e longitude 74°54'E e sugerimos que não exista contraste de magnetização expressivo abaixo da base do monte submarino.

Análise comparativa de três métodos de inversão aplicados a dados gravimétricos e magnéticos em perfil

Tradicionalmente, o método dos mínimos quadrados tem sido empregado na inversão não linear de dados de campo potencial. No caso em que as observações dos campos gravimétrico ou magnético contém apenas ruído Gaussiano. O método dos mínimos quadrados não apresenta problemas. Entretanto, quando as observações são perturbadas por ruído não Gaussiano, ou mesmo por ruído não aleatório, como é o caso de muitos ruídos geológicos, o método dos mínimos quadrados torna-se bastante ineficiente, e métodos alternativos devem ser empregados a fim de produzir interpretações realísticas. Neste trabalho, uma comparação é feita entre os métodos dos mínimos quadrados, dos mínimos absolutos e do ajuste-M, aplicados à inversão não linear de dados de campo potencial. A comparação é efetuada usando-se dados teóricos, onde diversas situações geológicas são simuladas. Os resultados mostram que na presença de ruído geológico, caracterizado por pequeno corpo raso acima do corpo principal, ou por corpo grande, adjacente ao corpo principal, o ajuste-M apresenta desempenho muito superior ao dos mínimos quadrados e dos mínimos absolutos. Na presença de ruído Gaussiano, entretanto, o ajuste-M tem um desempenho inferior aos outros dois métodos. Como o ruído Gaussiano é um ruído branco, parte dele pode ser removido por um filtro passa baixa adequado, sem muita perda do sinal, o que não ocorre com o ruído geológico que contém componentes importantes de baixo número de onda. Desse modo o ajuste-M se torna uma ferramenta importante na interpretação de áreas geologicamente complexas, onde é comum a contaminação das anomalias por ruído geológico. Os três métodos em estudo são aplicados a uma anomalia magnética real causada por uma intrusão de diabásio em forma de dique, em sedimentos arenosos da formação Piauí na Bacia do Parnaíba. Os três métodos apresentaram resultados semelhantes indicando que tanto o nível de ruído Gaussiano como geológico são baixos nesta anomalia.

Equipamento eletromagnético para prospecção geofísica e modelamento reduzido

Nas últimas décadas, a prospecção por métodos eletromagnéticos vem-se constituindo numa técnica eficiente para prospecção mineral. O objetivo deste trabalho foi desenvolver um equipamento para prospecção eletromagnética quantitativa de corpos condutores, através do método dipolo-dipolo, podendo ainda ser usado em modelos reduzidos. Eletricamente, o sistema mede grandezas relacionadas ao acoplamento indutivo entre duas bobinas: transmissora e receptora. Elas são dispostas na superfície da terra, afastadas entre si, e a terra, desse modo, constitui o núcleo acoplador. Quando existem corpos condutores nas proximidades, estes são denunciados por alterações no comportamento do sinal induzido na bobina receptora. O equipamento compreende dois conjuntos: o transmissor e o receptor, além de acessórios. O transmissor gera um campo eletromagnético nas freqüências de 520 e 3.090 Hz, e um sinal de referência para o receptor, o qual é enviado através de um cabo. O receptor, inicialmente, separa o sinal induzido pelos campos secundários gerados por condutores, do campo normalmente recebido, quando a condutividade da subsuperfície é relativamente uniforme (campo primário). Em seguida, decompõe esse sinal em duas componentes ortogonais, uma em fase, e outra em quadratura com o campo primário. Através de duas escalas de precisão, as amplitudes dessas componentes são mostradas como percentagens do campo primário, com precisão de 1%. A sensibilidade do receptor é de 0,5 μV. O circuito eletrônico foi rigorosamente testado com preciso instrumental de laboratório. Em seguida, testou-se sua aplicação no Laboratório de Modelo Reduzido Eletromagnético do NCGG, refazendo-se experiências clássicas, encontradas na literatura especializada. No campo, foi experimentado próximo da cidade de Araci no Estado da Bahia, em áreas prospectadas pela "Rio Doce Geologia e Mineração S/A-DOCEGEO". Em ambos os casos, verificaram-se bons resultados.

Interpretação de anomalias magnéticas usando integrais finitas dos momentos das componentes

O objetivo do trabalho consiste em desenvolver uma metodologia para determinar os parâmetros de um corpo causador de uma anomalia magnética. Essa metodologia baseia-se na utilização de expressões analíticas, deduzidas para as integrais finitas dos momentos da anomalia de intensidade total e das componentes H_x, H_y e H_z ao longo dos três eixos de um sistema de coordenadas cartesianas. Por meio do ajuste entre essas expressões analíticas e as integrais numéricas finitas dos momentos das componentes magnéticas obtidas a partir de um levantamento geofísico, pode-se computar, através de um processo iterativo, os parâmetros de magnetização, posição, profundidade e dimensões de um corpo anômalo. No caso em que são conhecidas as medidas de somente uma componente magnética, ainda é possível aplicar o método pois as demais componentes podem ser obtidas por um esquema de filtragem matemática. A metodologia foi testada com sucesso para modelos dipolo pontual e linha de dipolos os quais são muito utilizados em Geofísica para interpretar anomalias magnéticas produzidas por corpos geológicos que possuem uma ou duas dimensões horizontais muito menores do que a profundidade.

Métodos geofísicos aplicados ao ultramafito da Serra de Quatipuru - Estado do Pará

O objetivo deste trabalho foi a aplicação de métodos geofísicos para a avaliação econômica de ocorrência de níquel, cromo e amianto, em uma área de 6,25 km² - Malha 2 - Serra de Quatipuru, localizada a cerca de 50 km a SW de Conceição do Araguaia, no sul do Estado do Pará. Inicialmente foi aplicado o método magnético em caráter de reconhecimento, usando magnetômetros de campo total e componente vertical visando localizar zonas anômalas de amianto e cromita, uma vez que estes minerais se encontram associados com a magnetita. Em seguida foi usado o método eletromagnético de "dip angle" com o objetivo de definir zonas condutoras associadas a possíveis depósitos minerais ou características estruturais. Com base na interpretação destes dois métodos, definiu-se os locais de aplicação do método de polarização elétrica induzida-resistividade, que teve caráter de detalhe visando como alvo principal depósitos de sulfetos disseminados. As possíveis ocorrências minerais na área foram correlacionadas com anomalias geofísicas e confirmadas através de testemunhos de sondagem concluindo-se não serem de importância econômica. Contudo a aplicação destes métodos geofísicos foi útil, pois possibilitou o levantamento detalhado da área, que de outra forma só poderia ser realizado por meio dos métodos mais convencionais, tais como sondagem mecânica através de perfurações, os quais são muito caros e demorados.

A capa carbonática do sudoeste do cráton amazônico, estado de Rondônia: nova ocorrência e extensão dos eventos pós-glaciação marinoana (635 Ma)

Na porção oeste da Bacia dos Parecis, Estado de Rondônia, inserida no sudoeste do Cráton Amazônico, rochas carbonáticas expostas nas bordas dos grábens Pimenta Bueno e Colorado têm sido consideradas como parte do preenchimento eopaleozoico da bacia. A avaliação das fácies/microfácies e quimioestratigrafia dessas rochas nas regiões de Chupinguaia e Pimenta Bueno, confirmou a ocorrência de dolomitos rosados que sobrepõem, em contato direto, diamictitos glaciais previamente interpretados como depósitos de leques aluviais. Trabalhos prévios reportaram excursão negativa de δ¹³C, também confirmados neste trabalho, com variações entre -4.6 e -3,8‰VPDB em Chupinguaia e média de - 3,15‰VPDB em Pimenta Bueno. Esse padrão, de sedimentação e quimioestratigráfico, ausente nas rochas paleozoicas, é comumente encontrado nos depósitos carbonáticos anômalos do Neoproterozoico. No sul do Cráton Amazônico, Estado do Mato Grosso, rochas com essas mesmas características são descritas como capas carbonáticas relacionadas à glaciação marinoana (635 Ma). Neste trabalho, consideramos que os dolomitos rosados sobre diamictitos, em Rondônia, fazem parte do mesmo contexto das capas carbonáticas encontradas no Mato Grosso. Adicionalmente, destaca-se o contato brusco e deformado do dolomito sobre o diamictito, presente em ambas as ocorrências, configurando-se uma das feições típicas das capas carbonáticas do Cráton Amazônico. Essa relação paradoxal, entre diamictito e dolomito, tem sido interpretada como produto da mudança rápida das condições atmosféricas de icehouse para greenhouse, e a deformação da base foi gerada pelo rebound isostático. A capa carbonática de Rondônia compreende duas associações de fácies (AF2 e AF3) que recobrem depósitos glacio-marinhos compostos por paraconglomerados polimíticos (Pp), e arenito seixoso laminado (Asl), da AF1. A AF2 consiste em dolomudstone/dolopackstone peloidal com laminação plana a quasi-planar e com truncamentos de baixo-ângulo (fácies Dp), megamarcas onduladas (fácies Dm) e laminações truncadas por ondas (fácies Dt), interpretada como depósitos de plataforma rasa influenciada por ondas. Esta sucessão costeira é sucedida pela AF3, que compreende as fácies: dolomudstone/dolopackstone e dolomudstone/dolograinstone com partição de folhelho (Df) e siltito laminado (Sl). A fácies Df compreende um pacote de 6 metros de dolomito com partição de folhelho, apresentando lâminas de calcita fibrosa (pseudomorfos de evaporito) e dolomitos com laminações onduladas de corrente. Sobrejacente à fácies Df, ocorre a fácies Sl, apresentando 5 metros de siltito argiloso com laminação plana. Esta associação é interpretada como depósitos de plataforma rasa influenciada por maré, sendo sobreposta discordantemente, em contato angular, por depósitos glaciais do Eopaleozoico. Os valores isotópicos de C e O são negativos e refletem o sinal primário do C. No entanto, pode-se considerar uma leve influência da diagênese meteórica no sinal. As principais quebras nos sinais negativos podem estar associadas à influência meteórica, expressa pela substituição e preenchimento de poros por calcita e pela proximidade de superfícies estratigráficas, os quais refletem alguns padrões de alteração diagenética, representados nos sinais mais negativos. Diferentemente da capa carbonática do Mato Grosso, a capa de Rondônia possui níveis de pseudomorfos de evaporito e dolomitos com partição de folhelho (ritmito), em sucessão de fácies marinha rasa, onde os dolomitos de plataforma rasa influenciada por ondas passam para ritmitos e siltitos de plataforma rasa influenciada por maré (zona de inframaré), configurando uma sucessão retrogradante. Esta nova ocorrência de capa carbonática modifica a estratigrafia da base da Bacia dos Parecis, ao passo que exclui essas rochas carbonáticas da sequência eopaleozoica. Além disso, fornece informações que permitem reconstruir melhor a paleogeografia costeira da bacia neoproterozoica que acumulou os depósitos da plataforma carbonática do Grupo Araras, bem como estende os eventos pós-marinoanos ligados à hipótese do Snowball/Slushball Earth para o sudoeste do Cráton Amazônico, exposto no Estado de Rondônia.

Contexto geológico, estudos isotópicos (C, O e Pb) e associação metálica do depósito aurífero Tocantinzinho, domínio Tapajós, Província Tapajós-Parima

O depósito Tocantinzinho, localizado em um lineamento de direção NW–SE, a SW de Itaituba (PA), é atualmente o maior depósito aurífero conhecido na Província Tapajós. Está hospedado no granito homônimo, essencialmente isótropo, no qual dominam rochas sieno e monzograníticas, que foram fraca a moderadamente alteradas por fluidos hidrotermais. Microclinização (mais precoce), cloritização, sericitização, silicificação e carbonatação (mais tardia) são os mais importantes tipos de alteração. O principal estágio de mineralização é contemporâneo à sericitização/silicificação e é representado por vênulas com sulfetos (pirita ± calcopirita ± galena ± esfalerita) e ouro associado, as quais mostram localmente trama stockwork. Além de teores expressivos de Cu, Pb e Zn, são anômalos, em algumas amostras, os de As, Bi e Mo. A relação dos teores do Au com os dos metais-base é aleatória e as razões Au/Ag variam de 0,05 a 5,0. O Au é mais enriquecido nas porções com maior abundância de sulfetos de metais-base, embora ocorra principalmente incluso na pirita. Monocristais de zircão, extraídos do granito Tocantinzinho, forneceram idade Pb-Pb média de 1982 ± 8 Ma, permitindo interpretá-lo como uma manifestação magmática precoce do arco Creporizão. Valores de δ¹³C_PDB em calcita do estágio de carbonatação, dominantemente entre -3,45 e -2,29‰, são compatíveis com fonte crustal profunda, quiçá carbonatítica, enquanto os de δ¹⁸O_SMOW (+5,97 a +14,10‰) indicam forte contribuição magmática, ainda que mascarada por influxo de águas provavelmente superficiais. Estudos de inclusões fluidas em andamento revelam a presença de fluidos aquocarbônicos, cujo CO₂ poderia ter estado dissolvido no magma granítico em vez de ser relacionado à zona de cisalhamento. Os dados até aqui disponíveis permitem classificar o depósito aurífero Tocantinzinho como do tipo relacionado à intrusão.

Nutrient and phytoplankton biomass in the Amazon River shelf waters

O estuário do rio Amazonas é notável na Plataforma Continental do Amazonas, onde a presença das águas fluviais foi detectada, mesmo durante o período da diminuição da descarga desse rio, pelos baixos valores de salinidade e altos valores de nutrientes. Contudo, a presença das águas oceânicas também foi marcante. Em relação às fases do ciclo do nitrogênio, o nitrogênio orgânico dissolvido foi a forma predominante, seguido do nitrogênio total particulado, nitrato, amônia e nitrito. Os dados de clorofila a indicaram uma área eutrófica onde há nitrogênio embora os valores da anomalia do nitrogênio inorgânico dissolvido tenham mostrado que ocorre maior remoção do que adição dessa forma nitrogenada na área em estudo. Os resultados das simulações com o modelo bidimensional MAAC-2D confirmaram que as águas ricas em nutrientes são deslocadas para noroeste (dois núcleos em 2,5°N-50°W e 4°N-51°W), pela ação de uma CNB mais forte durante o período de decaimento da vazão dos rios. No período de vazão máxima, estas lentes de águas ricas de nutrientes distribuem-se próximo de 0,5°N-48,5°W, bem como ao longo da plataforma Amazônica rasa (20m-50m profundidade, 1°N-3,5°N), como resultado do espalhamento da água doce de origem continental.