Estudo dos potenciais impactos das mudanças climáticas e de alterações na cobertura vegetal nos recursos hídricos na Região Central da Amazônia

Este trabalho busca avaliar e quantificar os impactos na disponibilidade hídrica local na região da Amazônia Central, decorrentes de possíveis mudanças climáticas e modificações na cobertura vegetal, por meio de um experimento de simulação numérica com o modelo de biosfera Common Land Model (CLM), no modo “off line”. Os resultados de 9 modelos climáticos acoplados Oceano-Atmosfera para três cenários de alterações climáticas do Painel Intergovernamental sobre Mudanças Climáticas (IPCC-AR4) foram utilizados para compor a base das forçantes climáticas do modelo de biosfera CLM. Em relação ao uso da terra, foram utilizados cenários de dinâmica de desmatamento para a região no caso “business as usual” previstos para cada ano do período de 2001 a 2050. A área de estudo compreende o domínio da drenagem da Bacia do Rio Cuieiras na Amazônia Central. A partir dos resultados dos modelos analisaram-se, para cada conjunto de simulações, as incertezas das projeções em relação à precipitação e a temperatura e o impacto no ciclo hidrológico terrestre, considerando a variabilidade entre os modelos e os cenários de emissões de CO2; bem como as alterações nas componentes do balanço de água e energia a superfície associada a variações progressivas na cobertura de floresta e sua substituição por pastagem. Os resultados indicam que diante de um cenário de mudanças climáticas que resultem em uma diminuição (aumento) persistente na precipitação média anual, tanto o escoamento quanto o armazenamento de água no solo serão diretamente afetados. Em relação às alterações na cobertura vegetal tanto as componentes dos balanços de água e energia foram significativamente afetados pela substituição da floresta por pastagem, apresentando redução na evapotranspiração e aumentos no armazenamento de água no solo e do escoamento total.

Caracterização de regimes de umidade em regiões tropicais: comparação entre floresta e savana

Este trabalho tem como objetivo investigar um método de classificação de regimes de umidade, baseado na caracterização de diferentes “estados” da Camada Limite Atmosférica Tropical (CLAT), tanto acima de uma área de floresta quanto acima de uma área de savana, de acordo com a metodologia proposta por Mahrt (1991). A partir dessa classificação é elaborado um aperfeiçoamento da mesma que incorpora tanto uma análise da estabilidade termodinâmica da CLAT para a área de floresta, quanto uma investigação sobre a contribuição da Energia Potencial Disponível para Convecção (CAPE) para a classificação dos regimes de umidade. Para essas análises foram utilizados dados de radiossondagens e de torres micrometeorológicas, coletados durante o período menos chuvoso de cada região. Esses dados foram obtidos durante experimentos de campo realizados nas áreas de estudo, sendo que para a área de floresta (Caxiuanã) se utilizaram os dados do experimento COBRA-PARÁ (realizado no período de 06 a 13/11 de 2006), enquanto que para a área de savana (Daly Waters) se utilizaram os dados do experimento “KOORIN” (realizado no período de 15/07 a 13/08 de 1974, na Austrália). A comparação entre os regimes de umidade de cada região sugere que, caso houvesse a savanização da Amazônia, a mesma apresentaria uma atmosfera seca, com a maior parte da energia utilizada para o aquecimento dessa atmosfera, com baixos índices de evapotranspiração, menores taxas de precipitação e inexistência de CAPE. Por outro lado, a análise da estabilidade da atmosfera para Caxiuanã mostrou que, contrariamente ao observado em experimentos na Amazônia Ocidental, na Amazônia Oriental, durante o experimento COBRA-PARÁ, os maiores valores de CAPE ocorreram às 18:00 HL, possivelmente, em decorrência da convergência de umidade que provém da baía de Caxiuanã através da circulação de brisa terrestre. Isso indica que nessa região os máximos de CAPE estiveram associados predominantemente aos campos de umidade e não aos de temperatura. Para essas condições de CLAT “perturbada” o espaço de fase proposto por Mahrt (1991) não caracteriza muito bem regimes de umidade associados a grandes valores da CAPE.

Dinâmica do carbono em pequenas bacias de drenagem sob uso de agricultura familiar na Amazônia Oriental

O presente estudo teve como objetivo avaliar a dinâmica do carbono em uma região da Amazônia Oriental, cujo uso da terra predominante é a agricultura familiar; a unidade espacial adotada neste estudo foram três pequenas bacias de drenagem. A dinâmica do carbono foi avaliada a partir de medidas hidrológicas e biogeoquímicas em águas dos igarapés Cumaru, Pachibá e São João entre junho de 2006 a maio de 2007. O ambiente aquático e a hidrogeoquímica fluvial foram caracterizados a partir de medidas in situ da condutividade elétrica, temperatura, pH e concentração de oxigênio dissolvido. Amostras de água foram coletadas e analisadas para determinação do carbono orgânico dissolvido (COD) e pressão parcial do dióxido de carbono (pCO 2 ). A partir dos valores de pCO 2 , foram calculadas as concentrações de carbono inorgânico dissolvido (CID). Já os fluxos de C02 foram medidos in situ e também calculados a partir do pC0 2 . Utilizando-se medidas de vazão instantânea a cada campanha mensal de campo, calcularam-se fluxos anuais de COD. A caracterização dos solos e do uso da terra nas porções estudadas das bacias, assim como os índices pluviométricos e fluviométricos, foram considerados na interpretação dos resultados. Podem-se enumerar como principais resultados o seguinte: 1) As características físico-químicas das águas fluviais das bacias estudadas retrataram seus solos ácidos, a vegetação ripária, e processos hidrológicos biogeoquímicos no ambiente aquático e terrestre, e com certa variabilidade sazonal; 2)0 pH e o oxigênio dissolvido se correlacionaram positivamente com o carbono dissolvido na coluna d'água; 3) O transporte de COD por unidade de área foi elevado quando comparado com outras bacias amazônicas, e mais intenso em períodos chuvosos; 4) O transporte de COD e a evasão de C0 2 pareceram responder positivamente à presença de vegetação secundária e floresta densa, e negativamente às atividades agropecuárias; e 5) As taxas de evasão de C0 2 foram elevadas comparando-as a outros rios amazônicos, e corroboram a hipótese de que pequenas bacias são importante fontes de C0 2 para atmosfera na região.

Matemática e Cartografia: como a cartografia pode contribuir no processo de ensino-aprendizagem da matemática?

O Ensino da Matemática vem sofrendo transformações, ao longo dos últimos anos, com o intuito de tornar-se mais prazeroso e menos desgastante para o aluno, contribuindo dessa forma para a melhoria do seu aprendizado. A busca do aperfeiçoamento do ensino da Matemática, através de metodologias, alternativas é a tônica no momento. A modelagem matemática representa uma dessas metodologias. Utilizando seus preceitos e orientações, o presente trabalho objetiva associar a cartografia como modelo matemático, as ferramentas contemporâneas que alicerçam o processo de ensino-aprendizagem da Matemática. A cartografia, entendida como um conjunto de tecnologias de interpretação e de representação gráfica da superfície terrestre, cujos fundamentos são essencialmente matemáticos, tem-se constituído uma grande ferramenta de representação de dados e informações da superfície por diversos profissionais. Neste trabalho, advoga-se a favor da sua aplicação através da modelagem matemática para o processo de fortalecimento do ensino-aprendizagem da Matemática.

Projeto e construção de um robô móvel AGV/ROV não-holonômico com habilidade para navegação autônoma do tipo wall-following

Esta dissertação de mestrado apresenta o projeto e a construção de um robô móvel terrestre denominado LOGBOT, com tração de movimento do tipo diferencial – com duas rodas motoras e uma roda livre para manter a estabilidade de sua estrutura em relação à superfície. O controle do robô dispõe dos modos de telemetria e autônomo. No modo de controle por telemetria (ROV), a comunicação do robô com a estação de controle é feita por radiofreqüência a uma distância de até um quilometro em ambientes externos, e até cem metros em ambientes internos. No modo de controle autônomo (AGV), o robô tem habilidade para navegar em ambientes internos e desconhecidos usando sempre a parede à sua esquerda como referência para a trajetória de seu movimento. A seqüência de movimentos para execução da trajetória é enviada para a estação de controle que realiza análises de desempenho do robô. Para executar suas tarefas no modo autônomo, a programação do robô conta com um agente inteligente reativo, que detecta características do ambiente (obstáculos, final de paredes, etc.) e decide sobre qual atitude deve ser executada pelo robô, com objetivo de contornar os obstáculos e controlar a velocidade de suas rodas. Os problemas de erro odométrico e suas correções com base no uso de informações sensoriais externas são devidamente tratados. Técnicas de controle hierárquico do robô como um todo e controle em malha fechada da velocidade das rodas do robô são usadas. Os resultados mostraram que o robô móvel LOGBOT é capaz de navegar, com estabilidade e precisão, em ambientes internos no formato de um corredor (wall following).

Padrões na estruturação de comunidades de mamíferos terrestres na América do Sul

Na análise das causas que levam ao aparecimento de padrões das estruturas em comunidades de mamíferos terrestres da América do Sul, sempre foi um desafio conseguir um perfil exato dos processos evolutivos observáveis em larga escala. Neste trabalho procura-se analisar os padrões de distribuição da mastofauna terrestre na América do Sul, procurando correlacionar as diferentes faunas com diferentes gráficos utilizando um refinamento do método de cenogramas. Levaram-se em conta variáveis ecológicas que teriam influenciado a distribuição das espécies e o tamanho que apresentam (alimentação, locomoção e porte). Verificou-se também a possibilidade de aplicação do método para a análise da fauna de mamíferos terrestres extintos na América do Sul. Os resultados apontam para a consistência do uso de tal metodologia nos conjuntos atuais e uma boa possibilidade de sua aplicação em comunidades extintas. Para tanto, considerações devem ser feitas acerca de lacunas de conhecimento sobre os padrões de distribuição das mastofaunas atuais, como as relações ecológicas e evolutivas entre os mamíferos e o continente sul-americano.

Colonização por anfíbios e lagartos de áreas reflorestadas no Platô Saracá, região de Porto Trombetas-Pará

Poucos estudos têm considerado a colonização pela fauna em áreas reflorestadas após mineração. Para determinar os padrões de colonização por anfíbios e lagartos de áreas de reflorestamento em Porto Trombetas, Pará, foram examinadas a composição, riqueza e abundância de espécies, e as características biológicas de anfíbios e lagartos que ocupam áreas reflorestadas. Também foi avaliado o efeito da estrutura da vegetação e da distância da floresta nativa sobre a comunidade de anfíbios dos reflorestamentos. Anfíbios e lagartos foram amostrados ao longo de oito campanhas em oito áreas de reflorestamento e quatro áreas de floresta nativa através de procura ativa e com a utilização de poças artificiais para a reprodução de anfíbios. Foram registradas 20 espécies de anfíbios e 20 espécies de lagartos, sendo 14 espécies de anfíbios e 11 de lagartos em reflorestamentos e 19 espécies de anfíbios e 16 de lagartos em floresta nativa. Entre os anfíbios, Leptodactylus sp., Osteocephalus oophagus e Allobates femoralis foram as espécies mais abundantes nos dois ambientes e entre os lagartos, Gonatodes humeralis e Leposoma guianense foram as espécies mais abundantes em reflorestamentos e floresta nativa, respectivamente. Espécies de anfíbios de reprodução terrestre ou que utilizam pequenos corpos d’água temporários para a desova e lagartos arborícolas foram os grupos mais abundantes nos reflorestamentos. Espécies fossoriais e semifossoriais de anfíbios e lagartos de liteira foram os principais grupos ausentes nos reflorestamentos, sugerindo que o atual estágio da sucessão da vegetação ainda não oferece microhábitats apropriados para algumas espécies. A riqueza de espécies de anfíbios foi maior em áreas com maior cobertura do dossel. Áreas com maior cobertura de dossel tiveram maior abundância de Leptodactylus sp., Osteocephalus oophagus e Allobates femoralis. Apenas quatro espécies de anfíbios utilizaram as poças artificiais para desova e não houve relação significativa do número de espécies que utilizaram estas poças com a distância para a floresta nativa ou com a cobertura do dossel. Osteocephalus oophagus desovou em poças a maiores distâncias e A. femoralis em poças mais próximas em relação à floresta nativa. Os resultados evidenciam que a fauna de anfíbios e lagartos nas diferentes áreas de reflorestamento é um subconjunto da fauna da floresta nativa, e compreendem espécies florestais que indicam a importância relativa destas áreas para a conservação da fauna local.

Composição, riqueza e abundância de espécies de anfíbios na região do médio rio Xingú

Este estudo teve como objetivo inventariar a fauna de anfíbios nos diferentes hábitats da região do Médio Rio Xingu. Para tanto, foram selecionadas seis localidades, sendo as amostragens realizadas em três incursões, entre novembro de 2007 e março de 2008, totalizando 48 dias. Foram identificados seis tipos de hábitats na região amostrada: floresta de terra-firme, floresta sazonalmente alagável (várzea), margem do rio, “pedral”, lagoa e área alterada. Para a coleta de dados foram utilizadas duas metodologias: armadilhas de interceptação e queda e procura ativa. As armadilhas foram instaladas apenas nas áreas de floresta de terra-firme, enquanto a procura ativa foi empregada em todos os tipos de hábitats identificados. As amostragens resultaram no registro de 56 espécies de anfíbios e outras oito espécies foram identificadas e registradas em estudos anteriores na área. A diversidade e riqueza de espécies foram maiores na localidade Caracol, onde predominou a floresta de terra-firme, e menor na localidade Ilha Grande, dominada por floresta sazonalmente alagável. Foram registrados nove modos reprodutivos na área de estudo, sendo todos observados na floresta de terra-firme e apenas três nos “pedrais”, o que parece refletir a baixa heterogeneidade ambiental dessa área. Para comparar as localidades estudadas foram realizadas duas análises de similaridade, uma para cada metodologia de coleta utilizada. A análise de similaridade dos dados de procura ativa apontou maior semelhança na composição das espécies entre as áreas de floresta de terra-firme do que entre os outros hábitats. A análise de agrupamento realizada entre a composição de espécies obtida nesse estudo e outros levantamentos realizados na Amazônia agrupou esta área com outra também localizada no Médio Xingu.

Influência de estruturas geológicas bidimensionais no campo geoeletromagnético na presença do eletrojato equatorial

A Terra atua como um grande magneto esférico, cujo campo assemelha-se àquele gerado por um dipolo magnético. Este campo apresenta mudanças de intensidade que variam com a localização e a hora local. A parte principal do campo geomagnético se origina no interior da Terra através de processos eletromagnéticos. Extensivos estudos mostraram ainda que existem contribuições de origem externa ao planeta, principalmente de origem solar. Dentre estas fontes há anomalias do campo magnético que surgem a partir de um aumento diurno da corrente elétrica em uma estreita faixa da ionosfera, de direção leste-oeste, centrada no equador magnético e denominada Eletrojato Equatorial (EEJ). Ocasionalmente estas correntes podem apresentar reversões de fluxo, sendo denominadas Contra-Eletrojato (CEJ). Vários autores têm estudado os efeitos do EEJ e CEJ sobre as observações geoeletromagnéticas. Eles estão interessados no efeito combinado do EEJ e estruturas geológicas condutivas 1-D e 2-D. Nestes trabalhos a estrutura 2-D sempre se apresentava paralela ao eletrojato, o que é uma hipótese bastante restritiva ao se modelar ambientes geológicos mais realistas, em que corpos bidimensionais podem ter qualquer strike em relação ao EEJ. Neste trabalho apresentamos a solução deste problema sem esta restrição. Assim, mostramos os campos geoeletromagnéticos devidos a estruturas bidimensionais que possuam strike oblíquo em relação ao EEJ, através de perfis dos campos elétrico e magnéticos calculados na superfície e formando direção arbitrária à heterogeneidade condutiva 2-D. Com esta resposta avaliamos ainda qual a influência que estruturas bidimensionais exercem sobre a resposta magnetotelúrica, sob influência do Eletrojato Equatorial. Durante o desenvolvimento deste trabalho, utilizamos o método de elementos finitos, tendo por fonte eletromagnética o EEJ e o CEJ, que por sua vez foram representados por uma combinação de distribuições gaussianas de densidade de corrente. Estas fontes foram decompostas nas direções paralela e perpendicular à estrutura 2-D, resultando nos modos de propagação TE₁ e TE₂ e TM acoplados, respectivamente. Resolvemos o modo acoplado aplicando uma Transformada de Fourier nas equações de Maxwell e uma Transformada Inversa de Fourier na solução encontrada. De acordo com os experimentos numéricos realizados em um modelo interpretativo da Anomalia Condutiva da Bacia do Parnaíba, formado por uma enorme estrutura de 3000 ohm-m dentro de um corpo externo condutivo (1 ohm-m), concluímos que a presença do CEJ causa uma inversão na anomalia, se compararmos com o resultado do EEJ. Concluímos também que para as frequências mais altas as componentes do campo elétrico apresentam menor influência da parte interna do corpo 2-D do que da parte externa. Já para frequências mais baixas este comportamento se observa com as componentes do campo magnético. Com relação à frequência, vimos os efeitos do “skin-depth”, principalmente nas respostas magnéticas. Além disso, quando a estrutura 2-D está paralela ao eletrojato, o campo elétrico é insensível à estrutura interna do modelo para todos os valores de frequência utilizados. Com respeito ao ângulo θ_h entre a heterogeneidade e a fonte, vimos que o modo TM se manifesta naturalmente quando θ_h é diferente de 0°. Neste caso, o modo TE é composto por uma parte devido à componente da fonte paralela à heterogeneidade e a outra devido à componente da fonte perpendicular, que é acoplada ao modo TM. Assim, os campos calculados têm relação direta com o valor de θ_h. Analisando a influência do ângulo entre a direção do perfil dos campos e o strike da heterogeneidade verificamos que, à medida que θ_h se aproxima de 90°, os campos primários tornam-se variáveis para valores de θ_p diferentes de 90°. Estas variações causam uma assimetria na anomalia e dão uma idéia da inclinação da direção do perfil em relação aos corpos. Finalmente, concluímos que uma das influências que a distância entre o centro do EEJ e o centro da estrutura 2-D, causa sobre as componentes dos campos está relacionado às correntes reversas do EEJ e CEJ, pois a 500 km do centro da fonte estas correntes têm máxima intensidade. No entanto, com o aumento da distância, as anomalias diminuem de intensidade. Nas sondagens MT, nós também usamos o EEJ e o CEJ como fonte primária e comparamos nossos resultados com a resposta da onda plana. Deste modo observamos que as componentes do campo geoeletromagnético, usadas para calcular a impedância, têm influência do fator de acoplamento entre os modos TE₂ e TM. Além disso, esta influência se torna maior em meios resistivos e nas frequências mais baixas. No entanto, o fator de acoplamento não afeta os dados magnetotelúricos em frequências maiores de 10^-2 Hz. Para frequências da ordem de 10^-4 Hz os dados MT apresentam duas fontes de perturbação: a primeira e mais evidente é devido à presença fonte 2-D (EEJ e CEJ), que viola a hipótese da onda plana no método MT; e a segunda é causada pelo acoplamento entre os modos TE₂ e TM, pois quando a estrutura bidimensional está obliqua à fonte 2-D temos correntes elétricas adicionais ao longo da heterogeneidade. Concluimos assim, que o strike de uma grande estrutura condutiva bidimensional relativamente à direção do EEJ ou CEJ tem de fato influência sobre o campo geomagnético. Por outro lado, para estudos magnetotelúricos rasos (frequências maiores de 10^-3 Hz) o efeito do ângulo entre a estrutura geológica 2-D e a direção do EEJ não é tão importante. Contudo, em estudos de litosfera frequências menores de 10^-3 Hz) o acoplamento entre os modos TE₂ e TM não pode ser ignorado.

Uma nova abordagem para interpretação de anomalias gravimétricas regionais e residuais aplicada ao estudo da organização crustal: exemplo da Região Norte do Piauí e Noroeste do Ceará

A interpretação de anomalias gravimétricas é de grande importância no estudo de feições geológicas que ocorrem na crosta terrestre. Esta interpretação é, no entanto, dificultada pelo fato das anomalias gravimétricas serem resultantes da soma total dos efeitos produzidos por todos os contrastes de densidades de subsuperfície. Desse modo, com o intuito de separar efeitos de feições mais profundas de efeitos de feições mais rasas, bem como a caracterização da geometria desses dois conjuntos de feições, apresentamos um método de separação das componentes regional e residual do campo e a subsequente interpretação de cada componente. A separação regional-residual de dados gravimétricos é efetuada através da aproximação do campo regional por um polinômio ajustado ao campo observado por um método robusto. Este método é iterativo e usa como aproximação inicial a solução obtida através do ajuste polinomial pelo método dos mínimos quadrados. O método empregado minimiza a influência de observações contendo forte contribuição do campo residual no ajuste do campo regional. A componente regional obtida a partir da separação regional-residual é transformada em um mapa de distâncias verticais em relação a um nível de referência. Esta transformação compreende duas etapas. A primeira consiste na obtenção da continuação para baixo da componente regional, que é pressuposta ser causada por uma interface suave separando dois meios homogêneos, representando a interface crosta-manto, cujo contraste de densidade é supostamente conhecido. A segunda consiste na transformação do mapa de continuação para baixo em um mapa de distâncias verticais entre o nível de continuação (tomado como nível de referência) e a interface. Este método apresenta duas dificuldades. A primeira está ligada à instabilidade, havendo portanto a necessidade do emprego de um estabilizador o que acarreta a perda de resolução das feições que se desejam mapear. A segunda, inerente ao método gravimétrico, consiste na impossibilidade da determinação das profundidades absolutas da interface em cada ponto, bastando entretanto o conhecimento da profundidade absoluta em um ponto, através de informação independente, para que todas as outras profundidades absolutas sejam conhecidas. A componente residual obtida a partir da separação regional-residual é transformada em um mapa de contrastes de densidade aparente. Esta transformação consiste no cálculo do contraste de densidade de várias fontes prismáticas através de uma inversão linear pressupondo que as fontes reais estejam das a uma placa horizontal, com contrastes de densidade variando apenas nas direções horizontais. O desempenho do método de separação regional-residual apresentado foi avaliado, através de testes empregando dados sintéticos, fornecendo resultados superiores em relação aos métodos dos mínimos quadrados e da análise espectral. O método de interpretação da componente regional teve seu desempenho avaliado em testes com dados sintéticos onde foram produzidos mapeamentos de interfaces bem próximas das estruturas reais. O limite de resolução das feições que se desejam mapear depende não só do grau do polinômio ajustante, como também da própria limitação inerente ao método gravimétrico. Na interpretação da componente residual é necessário que se postule ou tenha informação a priori sobre a profundidade do topo e espessura da placa onde as fontes estão supostamente confinadas. No entanto, a aplicação do método em dados sintéticos, produziu estimativas razoáveis para os limites laterais das fontes, mesmo na presença de fontes interferentes, e pressupondo-se valores para profundidade do topo e espessura da placa, diferentes dos valores verdadeiros. A ambiguidade envolvendo profundidade do topo, espessura e densidade pode ser visualizada através de gráficos de valores de densidade aparente contra profundidade do topo presumida para a placa para vários valores postulados para a espessura da placa. Estes mesmos gráficos permitem, pelo aspecto das curvas, a elaboração de uma interpretação semi-quantitativa das profundidades das fontes reais. A seqüência dos três métodos desenvolvidos neste trabalho foi aplicada a dados gravimétricos da região norte do Piauí e noroeste do Ceará levando a um modelo de organização crustal que compreende espessamentos e adelgaçamentos crustais associados a um evento compressivo que possibilitou a colocação de rochas densas da base da crosta a profundidades rasas. Este modelo ê compatível com os dados geológicos de superfície. É ainda sugerida a continuidade, por mais 200 km em direção a sudoeste, do Cinturão de Cisalhamento Noroeste do Ceará por sob os sedimentos da Bacia do Parnaíba, com base nas evidências fornecidas pela interpretação da anomalia residual. Embora esta seqüência de métodos tenha sido desenvolvida com vistas ao estudo de feições crustais de porte continental, ela também pode ser aplicada ao estudo de feições mais localizadas como por exemplo no mapeamento do relevo do embasamento de/bacias sedimentares onde os sedimentos são cortados por rochas intrusivas mais densas.

Paleomagnetismo de rochas vulvânicas do Nordeste do Brasil e a época da abertura do Oceano Atlântico Sul

Na primeira parte deste trabalho foram desenvolvidos estudos de magnetismo de rochas e paleomagnetismo em amostras de rochas vulcânicas do Nordeste brasileiro. As idades das amostras compreende os períodos Jurássico e Cretáceo. Com este objetivo foram amostradas quatro áreas tendo sido estudado um total de 496 amostras em 55 sítios. Para a coleta foi utilizada uma perfuradora portátil que extrai amostras de 2.5 cm de diâmetro. A orientação das amostras foi feita por meio de uma bússola magnética e de um clinômetro. Os espécimes foram submetidos a desmagnetizações por campo magnético alternado e em alguns poucos casos foi empregada a desmagnetização térmica. Atribuindo-se peso unitário a cada sítio foi determinada a direção média da magnetização remanescente característica de cada uma das áreas estudadas. As rochas vulcânicas do período Jurássico, localizadas na borda oeste da Bacia do Maranhão (Porto Franco-Estreito), apresentaram uma direção media em que D= 3.9°, I= -17.9° com α₉₅= 9.3°, k= 17.9, N= 15 e todos os sítios apresentaram polaridade normal. Para esta área foi determinado o polo paleomagnético de coordenadas 85.3°N, 82.5°E (A₉₅= 6.9º) que se localiza próximo a outros polos paleomagnéticos conhecidos para esse período. As rochas da borda leste da Bacia do Maranhão (Teresina-Picos-Floriano) de idade cretácica inferior apresentaram uma direção média de magnetização remanescente característica tal que D= 174.7°, I= +6.0º com α₉₅= 2.8º, k= 122, N= 21 e todos os sítios apresentaram polaridade reversa. O polo paleomagnético associado a elas apresentou por coordenadas 83.6°N, 261.0°E (A₉₅=1.9°) e mostrou concordância com outros polos sul americanos de mesma idade. No Rio Grande do Norte foi estudado um enxame de diques toleíticos também de idade cretácica inferior, cuja direção média da magnetização remanescente característica encontrada foi D= 186.6º, I= +20.6º com α₉₅= 14.0° e k= 12.9, N= 10. Os sítios desta área apresentaram magnetizações com polaridades normal a reversa. O polo paleomagnético obtido se localiza em 80.6°N e 94.8°E com A₉₅= 9.5°. O estudo das rochas vulcânicas da província magnética do Cabo de Santo Agostinho indicou para a região um valor de D= 0.4º, I= -20.6º com α₉₅= 4.8° e k= 114, N= 9 para a magnetização remanescente característica. Todos os sítios apresentaram polaridade normal e o polo paleomagnético determinado apresentou as seguintes coordenadas: 87.6ºN, 135ºE com A₉₅= 4.5º. Foi discutida a eliminação da variação secular das direções obtidas, de forma que cada polo apresentado nesta dissertação é verdadeiramente um polo paleomagnético. A análise dos minerais magnéticos portadores da remanência, efetuada por curvas termomagnéticas ou por difração de Raio-X, indicou na maior parte das ocorrências, a presença de titanomagnetita pobre em titânio. A presença de maguemita e algumas vezes hematita, na maior parte das vezes resultado de intemperismo, não anulou a magnetização termoremanente associada à época de formação da rocha, que foi determinada após a aplicação de técnicas de desmagnetização aos espécimes. Pelas curvas termomagnéticas obteve-se, para a maioria das amostras, uma temperatura de Curie entre 500 e 600ºC. Os casos mais freqüentes indicaram a ocorrência de titanomagnetita exsolvida, em que foram observadas a presença de uma fase próxima à magnetita e outra fase rica em titânio, próxima à ilmenita, resultado de oxidação de alta temperatura. A segunda parte do trabalho diz respeito à determinação da época de abertura do oceano Atlântico Sul por meio de dados paleomagnéticos. Entretanto ao invés de se utilizar o procedimento comumente encontrado na literatura, e que se baseia nas curvas de deriva polar aparente de cada continente, foi aplicado um teste estatístico que avalia a probabilidade de determinada posição relativa entre os continentes ser válida ou não, para determinado período em estudo. Assim foi aplicado um teste F a polos paleomagnéticos da África e da América do Sul, dos períodos Triássico, Jurássico, Cretáceo Inferior e Cretáceo Médio-Superior, tendo sido estudadas situações que reconstituem a posição pré-deriva dos continentes e configurações que simulem um afastamento entre eles. Os resultados dos testes estatísticos indicaram, dentro de uma probabilidade de erro de menos de 5%, que a configuração pré-deriva de Martin et al (1981) é compatível com os dados paleomagnéticos do Triássico, mas apresenta uma diferença significativa para os paleopolos de Jurássico, Cretáceo Inferior, Cretáceo Médio-Superior. Outras reconstruções pré-deriva testadas apresentaram o mesmo resultado. A comparação entre os polos paleomagnéticos da América do Sul e da África, segundo uma reconstrução que admite uma pequena abertura entre os continentes, como a proposta por Sclater et al (1977) para 110 m.a. atrás, indicou que os dados do Triássico não são compatíveis com este afastamento. Por outro lado os paleopolos do Jurássico e do Cretáceo Inferior, embora mais antigos que a data sugerida pela reconstrução, são consistentes com esta separação dentro de uma probabilidade de erro de menos de 5%. Os dados do Cretáceo Médio-Superior se mostraram consistentes com a reconstrução sugerida para 80 m.a. atrás por Francheteau (1973) e que propõe uma separação maior entre os continentes. Com base na premissa de deslocamentos de blocos continentais rígidos a análise dos resultados obtidos indicou que América do Sul e África estavam unidas por suas margens continentais opostas no período Triássico e que uma pequena separação entre estes continentes, provavelmente devida a uma rutura inicial, ocorreu no Jurássico e se manteve, então, aproximadamente estacionária até o início do Cretáceo Inferior. Esta conclusão difere da maior parte dos trabalhos que discutem a abertura do oceano Atlântico Sul. Os dados do Cretáceo Médio-Superior são compatíveis com um afastamento rápido e significativo entre os continentes naquele período.

Região do espaço que mais influencia em medidas eletromagnéticas no domínio da frequência: caso de uma linha de corrente sobre um semi-espaço condutor

Localizar em subsuperfície a região que mais influencia nas medidas obtidas na superfície da Terra é um problema de grande relevância em qualquer área da Geofísica. Neste trabalho, é feito um estudo sobre a localização dessa região, denominada aqui zona principal, para métodos eletromagnéticos no domínio da freqüência, utilizando-se como fonte uma linha de corrente na superfície de um semi-espaço condutor. No modelo estudado, tem-se, no interior desse semi-espaço, uma heterogeneidade na forma de camada infinita, ou de prisma com seção reta quadrada e comprimento infinito, na direção da linha de corrente. A diferença entre a medida obtida sobre o semi-espaço contendo a heterogeneidade e aquela obtida sobre o semi-espaço homogêneo, depende, entre outros parâmetros, da localização da heterogeneidade em relação ao sistema transmissor-receptor. Portanto, mantidos constantes os demais parâmetros, existirá uma posição da heterogeneidade em que sua influência é máxima nas medidas obtidas. Como esta posição é dependente do contraste de condutividade, das dimensões da heterogeneidade e da freqüência da corrente no transmissor, fica caracterizada uma região e não apenas uma única posição em que a heterogeneidade produzirá a máxima influência nas medidas. Esta região foi denominada zona principal. Identificada a zona principal, torna-se possível localizar com precisão os corpos que, em subsuperfície, provocam as anomalias observadas. Trata-se geralmente de corpos condutores de interesse para algum fim determinado. A localização desses corpos na prospecção, além de facilitar a exploração, reduz os custos de produção. Para localizar a zona principal, foi definida uma função Detetabilidade (∆), capaz de medir a influência da heterogeneidade nas medidas. A função ∆ foi calculada para amplitude e fase das componentes tangencial (H_x) e normal (H_z) à superfície terrestre do campo magnético medido no receptor. Estudando os extremos da função ∆ sob variações de condutividade, tamanho e profundidade da heterogeneidade, em modelos unidimensionais e bidimensionais, foram obtidas as dimensões da zona principal, tanto lateralmente como em profundidade. Os campos eletromagnéticos em modelos unidimensionais foram obtidos de uma forma híbrida, resolvendo numericamente as integrais obtidas da formulação analítica. Para modelos bidimensionais, a solução foi obtida através da técnica de elementos finitos. Os valores máximos da função ∆, calculada para amplitude de H_x, mostraram-se os mais indicados para localizar a zona principal. A localização feita através desta grandeza apresentou-se mais estável do que através das demais, sob variação das propriedades físicas e dimensões geométricas, tanto dos modelos unidimensionais como dos bidimensionais. No caso da heterogeneidade condutora ser uma camada horizontal infinita (caso 1D), a profundidade do plano central dessa camada vem dada pela relação p_o = 0,17 δ_o, onde p_o é essa profundidade e δ_o o "skin depth" da onda plana (em um meio homogêneo de condutividade igual à do meio encaixante (σ₁) e a freqüência dada pelo valor de w em que ocorre o máximo de ∆ calculada para a amplitude de H_x). No caso de uma heterogeneidade bidimensional (caso 2D), as coordenadas do eixo central da zona principal vem dadas por d_o = 0,77 r₀ (sendo d_o a distância horizontal do eixo à fonte transmissora) e p_o = 0,36 δ_o (sendo p_o a profundidade do eixo central da zona principal), onde r₀ é a distância transmissor-receptor e δ_o o "skin depth" da onda plana, nas mesmas condições já estipuladas no caso 1D. Conhecendo-se os valores de r₀ e δ_o para os quais ocorre o máximo de ∆, calculado para a amplitude de H_x, pode-se determinar (d_o, p_o). Para localizar a zona principal (ou, equivalentemente, uma zona condutora anômala em subsuperfície), sugere-se um método que consiste em associar cada valor da função ∆ da amplitude de H_x a um ponto (d, p), gerado através das relações d = 0,77 r e p = 0,36 δ, para cada w, em todo o espectro de freqüências das medidas, em um dado conjunto de configurações transmissor-receptor. São, então, traçadas curvas de contorno com os isovalores de ∆ que vão convergir, na medida em que o valor de ∆ se aproxima do máximo, sobre a localização e as dimensões geométricas aproximadas da heterogeneidade (zona principal).

Aspectos fundamentais das medidas e interpretação de registros paleomagnéticos em rochas sedimentares da Formação Longá-Bacia do Parnaíba

Considerando que o número de dados paleomagnéticos para o Paleozóico Superior, principalmente Devoniano da América do Sul, é ainda insuficiente não permitindo que uma Curva de Deriva Polar confiável possa ser construída, foram amostrados 43 níveis estratigráficos da Formação Longá, os quais foram estudados paleomagneticamente com o propósito de contribuir para que uma Curva mais confiável possa ser construída. A amostragem foi feita segundo o método dos blocos ao longo dos perfis Teresina, Barras, Batalha, na rodovia PI-13 e Floriano, Nazaré do Piauí, Oeiras nas rodovias PI-24 e BR-230, no estado do Piauí. Os tratamentos foram iniciados no laboratório de Paleomagnetismo do IAG-USP e complementados no do NCGG-UFPa. Utilizou-se as técnicas de desmagnetização progressiva por campos alternados até 700 e/ou temperaturas até 670-700ºC. A interpretação dos resultados foi feita por meio dos diagramas vetoriais de Zijderveld, pelas curvas J-T/C de variação da intensidade magnética com os campos ou temperaturas, e pelos gráficos de variação na direção do vetor magnetização. Os cálculos de direção média e polos foram feitos segundo a Estatística de Fisher (1953). Foram identificadas 4 direções de magnetização remanente: 1. Uma secundária de origem química (CRM) e polaridade reversa, cujo mineral responsável é a hematita produzida provavelmente por alteração deutérica a partir da magnetita. Esta magnetização (identificada pela letra B) quando datada paleomagneticamente (coordenadas do polo 80ºS, 3°E, A₉₅ = 13.6°) indicou idade correspondente ao intervalo Carbonífero-Permiano. 2. Uma componente isotérmica (IRM) dura de espectro totalmente superposto à magnetização inicial que não foi afetada por nenhum dos tratamentos. Esta magnetização foi denominada D e apresentou direção muito estável em torno do ponto de declinação = 234.23º e inclinação = 41.94º. 3. Um grupo de direções de magnetização de origem viscosa (VRM) moles, identificados pela letra C, cuja direção média é dada por: declinação = 15º e inclinação = -20º. Foram removidas a temperatura entre 300 – 600ºC. 4. Finalmente uma magnetização principal, de polaridade normal, denominada A, provavelmente de origem detrítica (DRM) cujo correspondente polo paleomagnético (de coordenadas: 48ºS, 331.7ºE; A₉₅ = 9.9º) mostrou-se compatível com a idade da Formação (Devoniano Superior). Esta magnetização foi considerada inicial. Os polos paleomagnéticos correspondentes às magnetizações A e B, juntamente com outros da América do Sul, foram rotacionados para a África segundo a configuração pré-deriva de Smith e Hallam (1970) e comparados a polos Africanos e Australianos de mesma idade, mostrando-se coerentes. Suas polaridades também estão em acordo com as escalas magnetoestratigráficas publicadas por Irving e Pullaiah (1976) e Khramov e Rodionov (1981).

Avaliação da influência do eletrojato equatorial na interpretação de dados magnetotelúrico bidimensionais

O campo magnetotelúrico em regiões equatoriais viola a condição de ondas planas por causa de uma fonte fortemente concentrada na direção E-W na ionosfera, denominada eletrojato equatorial. No presente trabalho, procurou-se analisar a resposta magnetotelúrica de fontes que simulam o efeito do eletrojato equatorial. Foram considerados dois tipos de fontes para simular o eletrojato: uma linha infinita de corrente e uma distribuição gaussiana de densidade de corrente em relação a uma das coordenadas horizontais. A resistividade aparente foi obtida através da relação de Cagniard e comparada com os resultados de ondas planas. É mostrada também a comparação entre a fase da impedância na superfície, para os três tipos de fontes (ondas planas, eletrojato gaussiano e linha de corrente). O problema de meios com heterogeneidades laterais foi resolvido em termos de campos secundários, sendo as equações diferenciais solucionadas através da técnica de elementos finitos bidimensionais. Os resultados mostram que o eletrojato tem pouca influência nas respostas (resistividade aparente e fase) de estruturas geológicas rasas. Entretanto, a influência pode ser considerável nas estruturas profundas (maior que 5000 m), principalmente se suas resistividades são altas (maior que 100 Ω.m). Portanto, a influência do eletrojato equatorial deve ser considerada na interpretação de dados magnetotelúricos de bacias sedimentares profundas ou no estudo da crosta terrestre.

Desenvolvimento de um magnetometro a precessão nuclear

O objetivo desta tese foi desenvolver um magnetômetro à precessão nuclear para prospecção geofísica e estações-base magnéticas. O magnetômetro à precessão nuclear mede a intensidade total do campo magnético. Seu funcionamento é baseado na ressonância magnética nuclear. A medida de campo é feita pela de terminação da freqüência de precessão de núcleos de hidrogênio – prótons - de líquidos não viscosos no campo magnético terrestre. O magnetômetro é constituído de duas partes: o sensor e o instrumento de medida. O sensor é uma bobina solenoidal, cujo núcleo é preenchido com o líquido. Três líquidos diferentes foram testados; água, propanol e um querosene sintético. Optou-se pelo uso do querosene porque oferece maior amplitude no sinal de precessão, dando, conseqüentemente, maior relação sinal/ ruído. O sistema de medida contém os circuitos de sintonia e amplificação do sinal e, os circuitos lógicos para a programação da operação e contagem da freqüência de precessão. Cada ciclo de medida tem duração de 3 segundos, sendo 2,3s para a polarização e 0,7s para a recepção do sinal. São possíveis dois modos de operação: manual, reciclando automaticamente e por controle remoto. O sinal de precessão é amplificado seletivamente em uma das 14 faixas de sintonia, que cobrem medidas entre 22000 e 95000 gammas. A freqüência de precessão é multiplicada por um fator de 64 e contada durante um tempo igual a 0,36699s, determinado com base na razão giromagnética do próton. O número de pulsos contados é numericamente igual ao valor do campo magnético em gammas. A resposta pode ser lida em mostradores digitais ou na saída BCD paralela quando operando por controle remoto. A precisão da medida é de 1 gamma. O instrumento foi testado no campo para avaliar a relação sinal/ruído, gradiente suportável e consumo de potência. Nos testes de aplicação do protótipo, foram obtidos dados de variação diurna e realizaram-se levantamentos magnético de reconhecimento e detalhe em um sítio arqueológico na Ilha de Marajó, Pará. As respostas dos testes foram comparados com dois magnetômetros comerciais - o GP-70, McPhar e o G-816, Geometrics e, ainda, com dados do Observatório Magnético de Tatuoca-Pa. Em todos os casos, a comparação dos dados mostrou bom desempenho do magnetômetro em teste.