RCAAP Repository

A obtenção da temperatura da superfície da terra (Ts) tem sido objeto de estudo de vários pesquisadores. Neste trabalho, a principal referência para estimar Ts é a temperatura de brilho (Tb) obtida a partir da conversão do número digital (ND) em radiância monocromática da banda termal do TM - Landsat 5. Para a obtenção da Ts deve-se conhecer, além da temperatura brilho, a emissividade de cada pixel, a transmitância e a temperatura média do ar (Ta). Neste sentido, apresenta-se uma correção atmosférica da temperatura de brilho. Utilizou-se uma imagem TM - Landsat 5 do dia 04 de dezembro de 2000, com cobertura sobre áreas irrigadas de Petrolina, PE (09º24'S; 40º30'W; 376 m) e dados de uma radiossondagem sobre a referida área. De acordo com os resultados, a diferença entre a Ts obtida com correção apenas da emissividade e com correções da emissividade e da atmosfera, resulta em diferenças inferiores a 1ºC. No entanto, a diferença entre a Ts corrigida e a Tb, é de aproximadamente 5ºC.

A migração pré-empilhamento de seção de tiro comum baseada na equação da onda tem produzido imagens sísmicas de excelente qualidade. Entretanto, este processo exige uma alta demanda computacional, principalmente em aplicações sísmicas 3-D, e seu alto custo impede o seu uso no dia-a-dia do processamento sísmico. Assim, métodos que reduzem o número de migrações requeridas para um experimento sísmico completo, como a migração de onda planas, são bastante atrativos no processamento de dados sísmicos. Os métodos de migração de ondas planas podem produzir imagens de alta qualidade, comparável às fornecidas pela migração de tiro comum com um menor custo computacional. A migração 2-D pré-empilhamento implementada neste trabalho é realizada sobre seções de parâmetro de raio constante, obtidas através da decomposição dos dados ordenados em famílias de receptor comum, aplicando-se a transformada tau-p. Em seguida, essas seções de ondas planas são extrapoladas em profundidade usando-se operadores de extrapolação Split-Step. No processo de migração, as imagens sísmicas em profundidade são obtidas somando-se as freqüências de interesse de todas as seções de parâmetro de raio migradas. 0 método de ondas planas foi testado com sucesso em dois conjuntos de dados sísmicos sintéticos. Os resultados alcançados com este novo método de migração de ondas planas são comparáveis aos fornecidos pela migração de tiro comum e possibilitando assim a obtenção de imagens de subsuperfície de boa qualidade a um baixo custo computacional.

Medições gamaespectrométricas, de resistividade elétrica e de susceptibilidade magnética foram tomadas em solos agrícolas situados próximos a cidade de Maringá, noroeste do estado do Paraná, sul do Brasil, objetivando caracterizar a distribuição espacial das concentrações de radionuclídeos (K, eU e eTh), da resistividade aparente e da susceptibilidade magnética. Na área agrícola ocorrem três tipos de solos: Nitossolo Vermelho Eutrófico, Latossolo Vermelho Eutroférrico textura argilosa, ambos derivados da Formação Serra Geral (basaltos do Cretáceo Inferior) e Latossolo Vermelho Distrófico textura média, oriundo do retrabalhamento das formações Serra Geral e Goio-Erê, esta última composta por arenitos do Grupo Caiuá (Cretáceo Superior). Constatou-se que as concentrações de radionuclídeos e os valores de susceptibilidade magnética nos solos mais argilosos foram maiores do que nos solos mais arenosos, em função, principalmente, da maior adsorção dos primeiros e da maior disponibilidade de minerais magnéticos. Os teores médios em ppm e Bq Kg-1 de K, eU e eTh na área com atividade antrópica foram respectivamente de 1766-54,75; 0,83-10,22 e 1,78-7,27, significativamente maiores do que em áreas virgens sem aplicação de fertilizantes (1101-34,15 de K; 0,14-1,69 de eU e 1,31-5,36 de eTh, respectivamente em ppm e Bq Kg-1). Houve correlações entre os teores de urânio e argila, urânio e susceptibilidade magnética, urânio e matéria orgânica e entre resistividade elétrica e argila. Foram também observadas variações nas concentrações de radionuclídeos em diferentes formulações de adubos utilizados nas culturas de soja e trigo. A resistividade elétrica aparente, cuja magnitude variou de 25 a 647 Ohm.m, e a susceptibilidade magnética, com oscilações entre 0,28 e 1,10 × 10-3 SI, decorrentes dos teores de argila e de minerais magnéticos, respectivamente, se mostraram importantes na discriminação dos solos na área de estudo, podendo ser incorporadas como ferramentas rápidas e de baixos custos no mapeamento sistemático de solos.

The Princeton Ocean Model (POM) was adapted to São Sebastião Channel (SSC) by nesting three numerical grids so as to study the seasonal changes of its circulation and thermohaline structure. Through monthly average conditions, the numerical model simulated reasonably the typical conditions of spring, summer, autumn and winter, filling the bottom channel with South Atlantic Central Water (SACW) in spring and in summer. In autumn and in winter this mass of water is not found, nevertheless its weaker signs appear in autumn, season in which the stronger signs of Tropical Water (never more than 50%) are found. Northeasterly winds on the Southeast Continental Shelf (SECS) grid are essential so that the SACW enters SSC. The bottom circulation obtained in SSC is essentially to northeast and associated to the intrusion of the SACW forced first by the northeasterly winds in SECS and second by the pressure force always bigger at the south entrance than in the north entrance and always bigger in summer than in winter. The superficial currents are southwestward, weakening in autumn and intensifying towards summer, with its maximum in this season.

De acordo com o Atlas de Manguezais, este ecossistema representa 8% de toda a linha de costa do planeta e um quarto da linha de costa da zona tropical, perfazendo um total de 181.077 km². Este trabalho objetiva quantificar a extensão dos manguezais de macromaré da costa nordeste do Pará e noroeste do Maranhão, aqui denominada de Costa de Manguezais de Macromaré da Amazônia (CMMA). O reconhecimento das áreas de manguezais e sua quantificação foram realizados a partir da utilização de imagens Landsat-7 ETM+, adquiridas em 1999 e 2000 e de um sistema de informações geográficas (SIG). A integração dos dados de sensores remotos, dados geológicos e oceanográficos permitiu o reconhecimento de cinco setores geomorfológicos, que abrangem uma superfície total de 7.591 km² de manguezais. Esta área representa a maior faixa de manguezais contínuos do planeta e corresponde a 56,6% dos manguezais do Brasil. Medidas prioritárias de conservação dos manguezais da Amazônia devem ser tomadas e pesquisas que busquem a melhor compreensão deste complexo e importante ecossistema devem ser financiadas e desenvolvidas.

Este trabalho apresenta um estudo de detalhe no aqüífero fraturado na região de Lindóia, Estado de São Paulo. A partir do estudo regional de favorabilidade de ocorrência de água subterrânea definiu-se a área de estudo a qual foi mapeada como excelente favorabilidade. O estudo de detalhe consistiu no mapeamento de fraturas em fotografias aéreas (escala 1:25.000) e realização de levantamento geofísico de eletrorresistividade (caminhamento elétrico e sondagem elétrica vertical). Realizou-se a inversão dos dados geofísicos e a análise integrada desses dados com os dados geológicos de campo, dados de poços, fontes e fraturas mapeadas em fotografias aéreas. Os resultados mostram que foi possível otimizar os trabalhos de campo com a obtenção de resultados satisfatórios em relação à definição de locais mais favoráveis para a locação de poços.

O embasamento adjacente às bacias sedimentares de Camamu e Almada é caracterizado, principalmente, por rochas das facies granulito e anfibolito, com idades arqueanas e paleoproterozoicas, pertencentes ao orógeno Itabuna-Salvador-Curaçá. As unidades em maior proporção neste contexto são os granulitos metatonalíticos associados aos granulitos básicos e metamonzoníticos pertencentes ao cinturão Itabuna. Em menor área, tem-se o granito de Teolândia e o granodiorito de Moenda associado aos anfibolitos da banda de Ipiaú, os charnockitos e charnoenderbitos do bloco Jequié, os sienitos neoproterozóicos e os diques máficos. Os teores de K, U e Th nessas rochas variam de 0,02 a 6,33% para o K, < 0,2 a 9,10 ppm para o U e < 0,4 a 64,38 ppm para o Th. Esses teores são mais altos nos charnockitos, no granodiorito de Moenda, no granito de Teolândia e nos sienitos, têm valores intermediários nos metatonalitos e metamonzonitos e são mais baixos nos granulitos básicos. As taxas de produção de calor nessas rochas são mais altas onde os teores de K, U e Th são maiores, variando de 0,58 a 5,57 µW m-3. As áreas de abrangência dessas litologias são, contudo, pequena quando comparada com aquelas dos metatonalitos, metamonzonitos e sienitos onde as taxas variam de 0,10 a 1,44 µW m-3, 0,23 a 5,55 µW m-3 e 0,60 a 2,24 µW m-3, respectivamente. Os granulitos básicos possuem as menores taxas, entre 0,06 a 0,36 µW m-3. A distribuição das litologias nas bordas das duas bacias sugere que, no embasamento sob seus sedimentos mais recentes, deve predominar os metatonalitos, seguidos pelos metamonzonitos e com uma participação significante de sienitos no caso da bacia de Almada. Para estas litologias, as taxas volumétricas de produção de calor, com o intervalo de um desvio padrão, são 0,41 ± 0,30 µW m-3 para os metatonalitos, 0,71 ± 0,57 µW m-3 para os metamonzonitos e 1,20 ± 0,51 µW m-3 para os sienitos.

Esse estudo interpreta os perfis acústicos de subfundo e modelos 3D submarinos de alguns fiordes adjacentes ao Campo de Gelo Patagônico Sul, Chile. As bacias aprofundadas mostram a morfologia irregular e distintas formas associadas ao sistema sedimentar de zonas de grounding line, de línguas de gelos flutuantes. Refletores acústicos exibem duas fácies principais: caóticas e estratificadas, segundo sua configuração interna e geometria externa. Os sedimentos provavelmente foram depositados durante a deglaciação durante o Holoceno. A geometria dos depocentros e as características dos refletores acústicos mostram evidências da influência da batimetria e da topografia pré-existentes. Devido às grandes profundidades das bacias, as margens terminais podem ter sido um misto de condições flutuantes e aterradas no fundo marinho, embora possa predominar o regime glacial temperado.

Foi desenvolvida neste trabalho uma nova técnica de modelagem analógica para IP-resistividade em laboratório. Para tanto, foram utilizados modelos geométricos cilíndricos e tabulares. Para a construção destes modelos usou-se cimento, areia quartzítica e grafita. Como a grafita apresenta um intenso fenômeno IP, variou-se o teor de grafita contido nos modelos. Assim, para diferentes formas e medidas dos modelos, experimentaram-se diferentes conteúdos de partículas polarizáveis. As medidas foram feitas no sentido de verificar o comportamento da resposta IP e resistividade dos modelos quanto à forma, conteúdo de grafita, profundidade e orientação do modelo estudado.

Este trabalho avaliou a atividade dos radionuclídeos em rochas sedimentares e ígneas da bacia do Rio Corumbataí, estado de São Paulo, e, conseqüentemente, a taxa de exposição. Esta bacia possui uma grande área com cultura de cana-de-açúcar, onde há intensa aplicação de fertilizantes fosfatados e também abriga o mais importante pólo cerâmico da América do Sul (Pólo Cerâmico de Santa Gertrudes). A caracterização química das amostras mostrou uma concentração de Ca > Mg > K > Na, sendo essa relação devido à composição mineralógica de cada formação. A média ponderada da taxa de exposição (35 nGyh-1) para as formações estudadas representa apenas 64% da média mundial devido à radiação gama terrestre, sendo que a porção sudeste da área exibe a mais alta taxa de exposição. Os resultados obtidos podem ser utilizados em estudos ambientais envolvendo a bacia do Rio Corumbataí, inclusive aqueles sobre o uso do solo.

A interpretação de dados de investigações geofísicas rasas pode ser complexa devido às incertezas a respeito dos alvos em subsuperfície. Para reduzir as ambigüidades e auxiliar a interpretação dos dados geofísicos de superfície foi instalado no campus da USP o Sítio Controlado de Geofísica Rasa do IAG-USP. Diferentes alvos foram instalados em várias profundidades ao longo de 7 linhas. Os alvos instalados no Sítio Controlado simulam alvos encontrados em estudos geotécnicos, ambientais e arqueológicos. Neste artigo são apresentadas as principais etapas envolvidas na construção do Sítio Controlado de Geofísica Rasa do IAG/USP e alguns resultados GPR ao longo de uma linha caracterizada por tambores metálicos vazios. A caracterização geofísica, sob condições controladas de campo, permitiu uma calibração das respostas geofísicas para os tambores metálicos instalados no SCGR-IAG. O perfil GPR de 250 MHz apresentou uma excelente concordância com a modelagem GPR 2D através do método das diferenças finitas no domínio do tempo. A migração-fk melhorou a resolução e o posicionamento espacial dos tambores metálicos. Os resultados GPR 3D na forma de cubo 3D, diagrama de cercas e depth slices, mostraram-se muito importantes para melhorar a caracterização dos alvos em subsuperfície. O SCGR-IAG constitui-se numa ferramenta importante para o Ensino, Pesquisa e Extensão Universitária em Geofísica, e será de grande importância para a nossa comunidade. Com a instalação deste Sítio Controlado, foi dado um importante passo em busca do "Estado da Arte em Geofísica Rasa" no Brasil, onde o desenvolvimento de novos equipamentos, software de modelagem e de inversão poderão ser testados e suas respostas calibradas.

A partir da análise do posicionamento por GPS de 38 Referências de Nível, num perfil localizado no Estado de São Paulo, são sugeridos procedimentos de uso do GPS em levantamentos geofísicos terrestres. Como referência, foi realizado o posicionamento relativo estático, com receptores de dupla freqüência, coletando dados durante sessões de 60 minutos. As coordenadas geodésicas foram determinadas por meio de diversas estratégias de processamento, variando o tipo de receptores utilizados (simples e dupla freqüência), o número de bases (1 e 3), o comprimento das bases (2 a 250 km) e a duração da sessão (10 a 60 minutos). As altitudes ortométricas foram determinadas no modo relativo, usando as 5 estações da Rede GPS do Estado de São Paulo (RGSP) mais próximas. A análise estatística dos resíduos das coordenadas e altitudes obtidas, em relação aos valores de referência, gerou informações sobre a acurácia dos posicionamentos geográfico e altimétrico por GPS em função dos parâmetros analisados. Com base nos resultados obtidos, foram definidos os métodos de posicionamento GPS para levantamentos geofísicos locais e regionais, conforme suas características. Assim, para se obter acurácia de 10 cm, com bases maiores de que 30 km, deve-se usar receptores de dupla freqüência e sessões de no mínimo 30 minutos. Com bases menores do que 30 km podem ser usados receptores de dupla ou simples freqüência e sessões de 30 minutos.

Pseudo walkaway noise tests were accomplished in study area in São Paulo city, Brazil, in order to evaluate the potentiality of seismic reflection technique in shallow geological/geotechnical investigation in urban environment, where there is usually the presence of different types as well as intensities of noises and the paved surface. Data acquisition was conducted with a 24 channel seismograph with 24 bit A/D converter, 28 Hz and 100 Hz geophones, as well as sledge hammer and seismic rifle as sources. The best results were obtained with: (i) 100 Hz geophones with 0.18 m spike and hammer source with steel plate; (ii) 100 Hz geophones coupled through clay on asphaltic coverage with hammer source applied on asphalt. Laboratory tests were also accomplished in shake table in order to evaluate geophone responses when they were coupled with different types of clay. Tests indicated clay employment is a good alternative to couple geophones, notwithstanding kaolinitic clays are less indicated.

Neste trabalho apresentamos a aplicação do operador de extrapolação em profundidade em duas etapas com interpolação (PSPI-SS), operador que faz uso da correção Split-Step em cada passo de extrapolação, conjuntamente com a condição de imagem de tempo de excitação, para realizar a migração pré-empilhamento em profundidade de dados sísmicos ordenados em famílias de tiro comum. Os métodos de migração pré-empilhamento que realizam a extrapolação dos dados ao longo dos receptores e a continuação para baixo das fontes sísmicas com operadores diferentes, são chamados de métodos híbridos de migração. Nesses métodos o campo de onda de cada fonte sísmica é determinado através do cálculo dos tempos de chegada da fonte a partir da aproximação de alta freqüência. Aqui utilizamos apenas os tempos das primeiras chegadas, que são calculados a partir da solução da equação iconal por diferença finitas (TGRID), por um método puramente geométrico (FAST) e também pelo método de traçamento de raios paraxiais (PARAX). Na implementação proposta, o campo de velocidades usado para o cálculo dos tempos de chegada precisa ser suavizado e o seu grau de suavização é determinado a partir dos resultados apresentados pela migração convencional Kirchhoff. A extrapolação dos dados dos receptores é feita usando-se o campo de velocidades original sem nenhuma suavização, através do operador PSPI-SS. Para testar os novos métodos de migração usamos os dados sintéticos Marmousi e também um conjunto de dados reais do Golfo do México. Os resultados das migrações com os métodos híbridos apresentaram imagens de boa qualidade, mostrando-se superiores a obtida pela migração Kirchhoff e comparáveis à fornecida pelo método PSPI-SS completo (operador PSPI-SS tanto para os tiros como para os receptores), com um desempenho computacional superior ao obtido com o método PSPI-SS completo.

O algoritmo desenvolvido neste trabalho é baseado no método das diferenças finitas aplicado à equação acústica da onda, assumindo que a Terra se comporta como um meio acústico, que permite implementar numericamente uma modelagem sísmica empregando malhas regulares em modelos que representam meios geológicos bidimensionais. As derivadas segundas presentes na equação da onda foram obtidas por expansões da série de Taylor de quarta ordem, para o caso espacial, e de segunda, para o temporal. Os modelos foram representados por dois tipos distintos de parametrizações: por blocos e polinomial trigonométrica. São apresentadas simulações do fenômeno de propagação de ondas em vários modelos, representados pelas duas parametrizações mencionadas, possibilitando, assim, a geração de sismogramas sintéticos para comparação. Os sismogramas obtidos, com a parametrização polinomial, mostram alguns comportamentos não desejados como: produção de reflexões artificiais, atenuação de reflexões, supressão de difrações, e pequenas alterações nos tempos de trânsito e amplitudes calculados. Por outro lado, como vantagens desta parametrização proposta, temos: economia de espaço de memória (pois um modelo complicado de velocidade pode ser representado por uma pequena quantidade de coeficientes, isto é: o modelo não é armazenado em um arquivo, mas é comprimido, ou contido, numa fórmula matemática ou numa representação analítica); produção de modelos de velocidades suavizadas (úteis na migração sísmica); sintetização do modelo num pequeno conjunto de coeficientes que se tornam parâmetros do modelo a serem estimados em possíveis procedimentos de inversão (independentemente do grau de complicação da geometria do modelo e da variação da velocidade sísmica); melhor representação de modelos realísticos por polinômios trigonométricos em comparação com a parametrização por blocos; e redução das ambigüidades normalmente presente nas técnicas de inversão, se uma quantidade suficiente de dados for disponível. Cabe acrescentar que: o efeito de Gibbs, presente em representações polinomiais, é contornado, no método das diferenças finitas, através da escolha conveniente dos nós da malha ou pelo aumento do número de coeficientes do polinômio.

A integração de dados geológicos e geofísicos permitiu um melhor entendimento da inter-relação entre estruturas oceânicas e continentais. Fenômenos tectônicos e magmáticos do Meso-Cenozóico do Sudeste Brasileiro podem ter sido condicionados pelo prolongamento das zonas de fratura oceânicas em direção à margem continental. O foco deste trabalho é a evolução da Zona de Fratura de Vitória-Trindade (ZFVT). A Zona de Fratura de Vitória-Trindade (ZFVT), localizada na latitude de 18º40' na Cordilheira Mesoatlântica (CMA), foi mapeada utilizando-se sísmica de reflexão, gravimetria e isócronas magnéticas. A ZFVT na CMA corresponde a um deslocamento axial em torno de 6 km, configurando uma descontinuidade de segunda ordem. O segmento assísmico da ZFVT é definido por uma depressão gravimétrica de direção E-W. Sua extensão para oeste corresponde à Cadeia Vitória-Trindade, representada por altos topográficos individualizados, que limitam desníveis crustais de mais de 400 m e províncias sedimentares de características e espessuras distintas. Esta região cortada pela intrusão de diques (alguns aflorando na superfície) e sills, indica instabilidade tectônica com a movimentação de falhas e ocorrência de estruturas transpressivas (estruturas em flor), elevando toda a coluna sedimentar incluindo o assoalho oceânico. Estas observações sugerem reativações recentes da crosta oceânica com esforços alternantes, ora compressivos ora distensivos, associados a vulcanismos e falhamentos normais. Evidências geológicas e geofísicas, além do mapeamento do alinhamento de epicentros de terremotos, sugerem que o prolongamento da ZFVT está relacionado ao Alto de Vitória e ao tectonismo e magmatismo na porção emersa, representado pelo Lineamento Sismo-Magmático Alcalino Trindade-Paxoréu / Alto do Paranaíba. Propõe-se que a ZFVT atuou e continua atuando como conduto para o magmatismo da Pluma Mantélica de Trindade.

Um dos problemas existentes em investigações ambientais é a determinação da extensão do contaminante no solo. Este trabalho mostra alguns resultados obtidos com os métodos eletromagnético indutivo, eletrorresistividade e georadar (GPR - Ground Penetrating Radar) aplicados para determinar e localizar áreas potencialmente contaminadas com resíduos industriais de um curtume. A área está localizada na cidade de Várzea Grande - MT, e atualmente está interditada e desativada para a deposição de resíduos. A aquisição de dados consistiu em perfis de caminhamento elétrico, caminhamento eletromagnético indutivo e perfis common offset de GPR. No caminhamento elétrico foi utilizado o arranjo dipolo-dipolo, com espaçamento entre os dipolos de 4 metros e 5 níveis de investigação. Os dados de condutividade elétrica aparente foram adquiridos com o equipamento EM-34, nos modos dipolo magnético horizontal e vertical a cada 2 metros, com espaçamento entre as bobinas de 10 metros. Nas aquisições de GPR foram utilizadas antenas não-blindadas de 200 MHz, empilhamento de 128 sinais e traços coletados a cada 0,25 metros. Após as coletas dos dados, foram realizadas as classificações e delimitações das áreas de acordo com o padrão das anomalias de cada método. As áreas foram classificadas em Áreas Contaminadas - AC (sigmaa > 40 mS/m, ro < 30 ohm.m e alta atenuação do sinal do georadar), Áreas Suspeitas de Contaminação - AS (sigmaa de 20 a 40 mS/m, ro de 30 a 50 ohm.m) e Áreas Não-Contaminadas - AN (sigmaa < 20 mS/m, ro > 50 ohm.m e baixa atenuação do sinal do GPR). Os resultados integrados possibilitaram a elaboração de mapas de contaminação, indicando uma área potencialmente contaminada com um sentido preferencial NW-SE. As informações geofísicas aliadas à informação geológica de superfície e de poço mostraram uma análise preliminar eficiente na identificação de áreas potencialmente contaminadas.

Neste trabalho são apresentados os resultados obtidos a partir da utilização do método GPR - Ground Penetrating Radar na caracterização de tambores plásticos enterrados no Sítio Controlado de Geofísica Rasa do IAG instalados dentro do campus da Universidade de São Paulo - USP. Os tambores plásticos simulam estudos de contaminação ambiental em áreas urbanas. Os tambores enterrados vazios, com água doce e com água salgada, representam respectivamente contaminantes resistivos, água proveniente de reatores nucleares e contaminantes condutivos. Na primeira etapa, os trabalhos de campo foram realizados antes da disposição dos alvos em subsuperfície, visando obter dados referentes ao "background " do terreno, ou seja, dos horizontes geológicos anteriores às perturbações. Foram construídos modelos sintéticos bidimensionais para simular numericamente as respostas geofísicas, visando auxiliar na interpretação dos dados e diminuir as incertezas nos resultados. Na segunda etapa, foram obtidos dados de GPR após a instalação dos alvos. A análise dos resultados foi realizada através da integração dos dados reais de GPR, das modelagens numéricas GPR 2D e das profundidades conhecidas dos alvos, os quais apresentaram uma excelente concordância, permitindo identificar com clareza todos os alvos e distingui-los quanto à assinatura geofísica. Os resultados qualitativos e quantitativos sobre o Sitio Controlado do IAG poderão ser utilizados em locais onde há pouca ou nenhuma informação sobre a subsuperfície.

O principal objetivo deste estudo foi detectar mudanças no microclima da cidade de Manaus (AM), Brasil e procurar indicar as possíveis causas destas flutuações e tendências observadas. Foram utilizados dados meteorológicos e da geotermia rasa. Através de dados meteorológicos determinamos um incremento nas temperaturas médias do ar de 0,27 ± 0,04ºC, durante os últimos 80 anos, baseado na normal climatológica de 1950 a 1979. Este incremento não tem sido contínuo, mas modulado pelas freqüências de meses com El Niño, La Niña e Neutros, detectados no Oceano Pacifico tropical, Região Niño 3.4. Através de dados geotérmicos determinamos um incremento de 3,17 ± 0,53ºC, que se estima esteja relacionado com o uso do solo, desde a segunda metade do século XIX. Verificamos que as perturbações transientes do subsolo de locais com cobertura vegetal são bons indicadores das variabilidades climáticas anuais recentes e que as análises dos perfis sob locais sem cobertura vegetal indicam mudanças no microclima atribuídas às atividades antrópicas, tais como, o desflorestamento e urbanização.

Este trabalho apresenta resultados de estudos visando à estimativa dos teores de umidade de solos superficiais com GPR em uma área teste na cidade de São Paulo/SP, Brasil, e à análise da acurácia dos processos envolvidos. Foram empregadas duas metodologias diferentes: i) ajuste linear dos tempos de chegada da onda direta no solo para diferentes afastamentos; ii) diferenças dos tempos de chegada da onda aérea e da onda direta no solo para um único afastamento. Utilizamos o equipamento da RAMAC Mala Geosciences com antenas de 50 MHz, 100 MHz e 200 MHz. A aquisição constituiu-se de diversos perfis multi-offset (múltiplos afastamentos) com o emprego da geometria WARR (Wide Angle Reflection and Refraction ) ao longo de um perfil em três períodos climáticos distintos (seco, intermediário e úmido). Sendo o incremento entre os traços dos perfis multi-offset 0,2 m para a antena de 50 MHz e 0,1 m para as antenas de 100 MHz e 200 MHz. Concomitantemente ao levantamento GPR foram coletadas amostras de solo a intervalos de 0,5 m desde a superfície do terreno até a profundidade aproximada de 4,0 m para a determinação de seus teores de umidade em laboratório. Os resultados do GPR apresentaram boa correlação com os dados obtidos através das análises de laboratório, sendo que os obtidos com as antenas de freqüência central de 200 MHz propiciaram os melhores resultados. Os dois métodos empregados apresentaram diferentes acurácias. O método da diferença dos tempos de chegada da onda aérea e da onda direta no solo foi o menos preciso devido à ocorrência de um deslocamento de tempo da reta de ajuste dos tempos de chegada da onda direta no solo no afastamento zero. Esses resultados mostraram que, nesses casos, é importante aplicar um fator de correção para melhorar a acurácia dos resultados finais.