Trajetórias tecnológicas e sistemas agrários: qualificação territorial utilizando geoprocessamento em Paragominas e Santarém/PA

A Amazônia tem sido observada principalmente por meio do fenômeno do desmatamento, utilizando recursos tradicionais do sensoriamento remoto como a quantificação de área desflorestada e seu posterior incremento anual, que parece se constituir em uma metodologia eficaz. Ratificando este raciocínio, constatei num levantamento de 16.591 multas aplicadas pelo IBAMA/PA no período de 2000 até 2008, que mais de 85,0% das autuações estavam relacionadas apenas à componente flora;e na jurisdição da Gerência do IBAMA em Santarém, oeste do Pará, em 2008, quase 60% das multas se deu por conta de desflorestamento, identificados via sensoriamento remoto. Ressalta-se que as análises de imagens de satélites por si só não definem elementos da superfície terrestre, contribuindo pouco para o entendimento e posterior intervenção da realidade. Nesse contexto, foram investigados e vetorizados 479 estabelecimentos rurais nas regiões de Paragominas e Santarém, estado do Pará, que possuem históricos de uso e ocupação distintos, qualificando-os segundo suas trajetórias tecnológicas prevalentes, na perspectiva apresentada por Costa, concretizando um importante passo para correção das distorções no desenvolvimento econômico, agregando informação ao dado de sensoriamento remoto. Aplicaram-se recursos geotecnológicos de métricas de paisagem, construíram-se banco de dados celular integrado com estatísticas e algoritmos de otimização probabilística, associando dados de classificação não supervisionada isodata (validadas com kappa= 0,87, classificação considerada“excelente”) com os tipos de produção coletados em campo, gerando uma classificação final "integrada"(kappa= 0,78, classificação “muito boa”). Na região de Paragominas, foram qualificadas 3 tipos de trajetórias tecnológicas, a camponesa T8 (domínio de culturas temporárias), a camponesa T3 e patronal T4 (especializadas em pecuária de corte). Em Santarém, revelaram-se 2 trajetórias, a camponesa T2 (forte presença de culturas permanentes, temporárias e sistemas agro florestais) e a patronal T7 (mutação da T4, aumentando a participação das culturas temporárias). A metodologia aplicada logrou êxito, espacializando as propriedades rurais segundo seus tipos de trajetórias tecnológicas e gerando classes de uso mais representativas, como cultura temporária e pastagem, mas que na classificação via sensoriamento remoto isodata é englobada pela classe “agropecuária”, possibilitando uma visão mais realista das atividades de produção desenvolvidas na área investigada, concretizando a geração de informações espaciais integrando dados de diferentes fontes e o aumento do poder de leitura do pixel.

Dinâmica no uso da terra de 1991 a 2010 da região semiárida do estado do Piauí: cajucultura e apicultura

Pós-graduação em Agronomia (Produção Vegetal) - FCAV

Aplicação de dados multisensor (SAR e ETM+) no reconhecimento de padrões de uso e ocupação do solo em costas tropicais: Costa Amazônica, Amapá, Brasil

Os estudos ambientais necessitam de informações sobre a cobertura e o uso da terra. Este trabalho apresenta a aplicação de dados de sensores remotos orbitais (óticos e de radares) na validação de padrões de uso e cobertura do solo na planície costeira amapaense para fins de mapeamento e reconhecimento da dinâmica natural e antrópica. Esta costa é submetida a uma dinâmica intensa devido à magnitude dos processos costeiros (marés-pororoca), sua localização geográfica, influenciada pelo rio Amazonas e pela Zona de Convergência Intertropical, e processos antrópicos associados à bubalinocultura. A análise foi realizada aplicando-se dados de satélite (JERS-1, RADARSAT-1, Landsat 7 e DEM do SRTM) digitalmente processados em abordagem multisensor, multiescala e multitemporal, correlacionada com dados pretéritos e informações de campo. A análise dos produtos gerados e dados colaterais permitiu distinguir oito padrões de uso e cobertura do solo: florestas de mangue, florestas de várzeas, campos arbustivos, áreas de vegetação campestre, campo antrópico, zona de intermaré, canal estuarino e lagos, além de feições morfológicas lineares associadas a estes padrões. Estas informações são importantes para o mapeamento dos ambientes costeiros e fundamentais para o reconhecimento da dinâmica na região.

Mapeamento de ambientes costeiros tropicais (Golfão Maranhense, Brasil) utilizando imagens de sensores remotos orbitais

Este trabalho apresenta os resultados do reconhecimento e mapeamento dos ambientes costeiros da região do Golfão Maranhense, Brasil, utilizando uma abordagem metodológica que incluiu: (a) análise integrada com base no processamento digital de imagens, ópticas Landsat-4 TM e SPOT-2 HRV, de imagens SAR (Synthetic Aperture Radar) do RADARSAT-1, e dados de elevação da SRTM (Shuttle Radar Topography Mission); (b) sistema de informações geográficas; e (c) levantamentos de campo relativos à geomorfologia, topografia e sedimentologia. Os ambientes costeiros, assim mapeados foram agrupados em quatro setores: Setor 1, com pântanos salinos, pântanos de água doce, lagos intermitentes e canal estuarino; Setor 2, abrangendo tabuleiro costeiro, planície de maré lamosa, planície fluvial, planície de maré arenosa, praias de macromaré, área construída e lagos artificiais; Setor 3, com manguezal, paleodunas e planície de maré mista; e Setor 4, constituído por dunas móveis. Além disso, foram também reconhecidos lagos perenes, deltas de maré vazante e planícies de supramaré arenosas. O processamento digital e a análise visual das imagens de sensores remotos orbitais, associados ao uso de sistemas de informações geográficas, mostraram-se eficazes no mapeamento de zonas costeiras tropicais, possibilitando a geração de produtos com boa acurácia e precisão cartográfica.

Fusão de dados multisensor para a identificação e o mapeamento de ambientes flúvio-estuarinos da Amazônia

A fusão de imagens multisensor tem sido um procedimento amplamente utilizado em função da natureza complementar dos vários conjuntos de dados. Este artigo compara o desempenho de quatro métodos diferentes para fusão de imagens Landsat-7 ETM+ e RADARSAT-1 W1. A comparação foi baseada nas características espectrais das imagens, utilizando-se de análise estatística e visual dos produtos gerados. Quatro métodos foram usados para a fusão das imagens Landsat-7 ETM+ e RADARSAT-1 W1: i) fusão do SAR (radar de abertura sintética) com a tríade selecionada pelo OIF (Optimum Index Factor); ii) realce por decorrelação da tríade selecionada pelo OIF, seguida da fusão com SAR; iii) ACP (Análise por Componentes Principais) para as seis bandas ETM+ do espectro refletido (1, 2, 3, 4, 5 e 7) e posterior fusão das três primeiras componentes principais (1CP; 2CP; 3CP) com o SAR; iv) SPC-SAR (Principal Componente Seletivo - SAR). O produto SPC-SAR mostrou melhor desempenho na identificação das feições costeiras e permitiu o realce mais efetivo dos diferentes ambientes.

Reconhecimento e mapeamento dos ambientes costeiros para geração de mapas de ISA ao derramamento de óleo, Amazônia oriental

Acidentes relacionados ao transporte de petróleo constituem uma ameaça às zonas costeiras de todo o mundo, pois se trata de um evento ecologicamente desastroso podendo causar sérios danos ao meio ambiente. Na costa nordeste do Pará ocorre intenso fluxo de navios e balsas transportando petróleo e derivados, daí a necessidade de estabelecer estratégias de prevenção e gestão dos impactos ambientais ao derrame de óleo para elaboração de planos de emergência. O objetivo deste trabalho é gerar o mapa dos ambientes costeiros e o mapa de índice de sensibilidade ambiental ao derrame de óleo dos ambientes costeiros de uma porção da costa nordeste do Pará. A metodologia consistiu no processamento digital de imagens Landsat ETM+7, associadas aos dados do MDE do SRTM (shuttle radar topography mission) e aos dados de campo, em ambiente de SIG (sistema de informação geográfica). Os resultados permitiram o reconhecimento e a descrição de cinco unidades geomorfológicas e de ambientes sedimentares; identificação e hierarquização de nove índices de sensibilidade ambiental (ISA); além da geração do mapa de índice de sensibilidade ambiental ao derramamento de óleo. As técnicas de sensoriamento remoto e SIG mostraram ser importantes ferramentas no reconhecimento dos ambientes costeiros e geração dos mapas propostos.

Costa de manguezais de macromaré da Amazônia: cenários morfológicos, mapeamento e quantificação de áreas usando dados de sensores remotos

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.

Use of synthetic aperture radar for recognition of Coastal Geomorphological Features, land-use assessment and shoreline changes in Bragança coast, Pará, Northern Brazil

Imagens de radar de abertura sintética (SAR) vem sendo bem mais utilizadas do que antes nas aplicações de geociências em regiões tropicais úmidas. Nesta investigação, uma imagem RADARSAT-1, na banda C, polarização HH adquirida em 1998 foi usada para o mapeamento costeiro e avaliação da cobertura da terra na área de Bragança, norte do Brasil. Imagem do radar aerotransportado GEMS-1000, na banda X, polarização HH, adquirida em 1972 durante o projeto RADAM foi também utilizada para avaliar as variações costeiras ocorridas nas últimas três décadas. A pesquisa tem confirmado a utilidade da imagem RADARSAT-1 para o mapeamento geomorfológico e avaliação da cobertura da terra, particularmente em costas de manguezal de macromaré. Além disso, um novo método para estimar as variações da linha de costa baseado na superposição de vetores extraídos de diferentes imagens SAR, com alta acurácia geométrica, tem mostrado que a planície costeira de Bragança tem estado sujeita a severa erosão responsável pelo recuo de aproximadamente 32 km² e acreção de 20 km², resultando em uma perda de área de manguezal de aproximadamente 12 km². Como perspectiva de aplicação, dados SAR orbitais e aerotransportados provaram ser uma importante fonte de informação tanto para o mapeamento geomorfológico, quando para o monitoramento de modificações costeiras em ambientes tropicais úmidos.

Environmental sensitivity index (ESI) mapping of oil spill in the amazon coastal zone: the PIATAM Mar project

A importância do monitoramento ambiental é medida pelos vários casos de derramamentos de óleo ocorridos no mundo durante as últimas três décadas. Isto tem incentivado as empresas e órgãos do governo envolvidos na prevenção e combate a estes acidentes a aperfeiçoarem cada vez mais os métodos, tanto preventivos como corretivos, para a minimização dos danos gerados por acidentes com derramamento de óleo. Este trabalho objetiva contextualizar de forma histórica como os acidentes com derramamento de óleo propiciaram o desenvolvimento de pesquisa tecnológica a partir de parcerias entre empresas de petróleo, agências de governo, universidades e institutos de pesquisa no Brasil, em especial na zona costeira Amazônica. Como resultado, índices de sensibilidade ambiental ao derramamento de óleo (ISA) foram definidos especialmente para a Amazônia costeira, onde processos fluviais e marinhos se encontram na foz do maior rio do mundo, o rio Amazonas. Perspectivas de pesquisa e respostas de emergência a acidentes são apresentadas, a fim de se conservar a diversidade socioambiental da mais importante região tropical do planeta.

Mapas de variabilidade do potencial de produtividade obtidos por satélites como instrumentos auxiliares às decisões sobre o manejo de nematóides na cultura da soja

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

Morfodinâmica costeira e o uso da orla oceânica de Salinópolis (Nordeste do Pará)

Os aspectos morfodinâmicos relacionados à erosão ou acresção da linha de costa são alguns dos assuntos analisados na gestão das zonas costeiras que vêem sendo tratada em todo mundo no sentido de monitorar e proteger essas zonas. Esta tese objetiva analisar o comportamento da morfodinâmica costeira de Salinópolis, relacionando-o ao uso da orla oceânica. A área de estudo foi compartimentada em três setores: Oeste (praias da Corvina e do Maçarico), Central (praia do Farol Velho) e Leste (praia do Atalaia). A metodologia consistiu na: (a) aquisição e tratamento de imagens multitemporais (1988-2001-2013) do satélite Landsat 5 TM, 7 ETM e 8 OLI; (b) aplicação de entrevistas/questionários com banhistas, (c) aquisição de dados de campo durante as estações chuvosa (26, 27, 28/04/2013) e menos chuvosa (04, 05, 06/10/2013); e (d) análise laboratorial para o tratamento dos dados adquiridos em campo (topografia das praias estudadas, amostragem de sedimentos superficiais das mesmas e com o uso de armadilhas, e medições oceanográficas de ondas, marés, correntes e turbidez). Foram feitas as representações gráficas dos perfis topográficos das praias, calculados os parâmetros estatísticos granulométricos de Folk & Ward (1957), as taxas do transporte sedimentar nas praias e os parâmetros morfométricos de Short & Hesp (1982), estes últimos foram calculados com o intuito de relacioná-los aos estados morfodinâmicos de praias propostos por Wright & Short (1984) e Masselink & Short (1993). Para a classificação da costa oceânica de Salinópolis em termos de uso e ocupação foi utilizado o decreto nº 5.300 de 7 de dezembro de 2004. A partir das pesquisas sobre a urbanização na costa e das obras situadas nos ambientes costeiros foi utilizada uma matriz proposta por Farinaccio & Tessler (2010) que lista uma série de impactos ambientais, e o quadro de geoindicadores do comportamento da linha de costa proposto por Bush et al. (1999), para a identificação de locais com vulnerabilidade à erosão ou acresção. Para as condições oceanográficas em cada praia e periculosidade ao banho nas mesmas, foram integralizados os dados de ondas, de correntes, de morfodinâmica praial e questionários aplicados com banhistas. Atualmente, a orla oceânica de Salinópolis possui diferentes características quanto à utilização e conservação, abrangendo desde a tipologia de orlas naturais (Classe A) até orlas com urbanização consolidada (Classe C). A primeira ocorre nos extremos da área de estudo e, a segunda, na região da sede municipal. Quatro tipos de praias foram identificados segundo a exposição marítima e o grau das condições oceanográficas: tipo 1 (Maçarico), tipo 2 (Corvina), tipo 3 (Farol Velho) e tipo 4 (Atalaia). O trecho de costa com maiores impactos ambientais e com elevada erosão costeira localiza-se na praia do Farol Velho. O grau de periculosidade ao banho foi de 4 (praia do Maçarico) a 7 (praia do Atalaia) – médio a alto grau de risco. As praias de Salinópolis apresentam declives suaves (< 1,5°), grandes variações na linha de costa entre as estações do ano (9,6 a 88, 4 m) e volume sedimentar variável dependendo do grau de exposição das praias ao oceano aberto. Predominou o estado morfodinâmico dissipativo (Ω>5,5) para estas praias, mas com ocorrência do estado de banco e calha longitudinais (4,7<Ω<5,5) no setor oeste. As macromarés na área de estudo apresentaram altura máxima de 5,3 m (Setor Central, durante a estação menos chuvosa) e mínima de 4 m no mesmo setor, durante a estação chuvosa. As correntes longitudinais foram mais intensas no setor leste (>0,45 m/s) durante as duas estaçoes do ano. As alturas de ondas foram também maiores no setor leste (máximo de 1,05 m durante a maré enchente na estação menos chuvosa) e os períodos de ondas foram mais curtos (<4,5 s) no setor oeste. A média granulométrica obtida dos sedimentos coletados na face praial apresentou escala mais freqüente entre 2,6 a 2,8 phi, indicando a predominância de areia fina. O grau de seleção predominante dos sedimentos foi de 0,2 a 0,5 phi (muito bem selecionados e bem selecionados), e da assimetria foi de positiva (0,10 a 0,30) e de aproximadamente simétrica (-0,10 a 0,10). O grau de curtose variou desde muito platicúrtica (<0,67) a muito leptocúrtica (1,50 a 3,00). Foram observados eventos de acresção sedimentar da estação chuvosa a menos chuvosa. De 22/07/1988 a 28/08/2013 (25 anos) também houve predomínio de acresção, onde o avanço médio linear da linha de costa foi de 190,26 m. O recuo médio linear obtido para toda área de estudo foi de -42,25 m. Áreas com maior erosão são pontuais: divisas das praias da Corvina e Maçarico, e Farol Velho e Atalaia. Os traps portáteis indicaram uma maior quantidade de sedimentos transportados longitudinalmente na estação menos chuvosa (Mín. 280 g/m³: enchente, setor oeste; Máx. 1098 g/m³: vazante, setor leste). Nos traps de espraiamento, o balanço entre a quantidade de sedimentos entrando e saindo nas praias foi menor no setor central (Mín. 80 g/m³: vazante, estação menos chuvosa; Máx. 690 g/m³: enchente, estação menos chuvosa). A circulação costeira sedimentar é proveniente, principalmente, do efeito das marés, com direção governada pela enchente e vazante dos rios que atravessam a costa. Os dados indicam o transporte longitudinal de sedimentos da ilha de Atalaia e rio Sampaio para o setor oeste e as margens das faixas praiais.

Detecção de mudança e sedimentação no estuário do Rio Coreaú

O uso de novas técnicas para estudar a evolução e preenchimento de vales incisos tem fornecido, ao longo dos anos, importantes resultados para entendermos como foi a evolução costeira brasileira. Neste contexto, esta tese teve como objetivo estudar a evolução do estuário do rio Coreaú, localizado no estado do Ceará, em diferentes escalas temporais, seja “Eventual” (meses, anos), “Engenharia” anos, decádas) e Geológica” (centenas, séculos, milênios), proposta por Cowell et al. (2003), com intuíto de avaliar se as transformações/alterações ao longo dos anos foram significativas ou não. Como resultados, obteve-se no primeiro objetivo, utilizando técnicas de sensoriamento remoto, a partir de imagens dos sensores TM, ETM⁺ e OLI do satélite Landsat 5,7 e 8 e LISS-3 do satélite ResourceSat-1 de 1985 a 2013, uma alteração mínima em relação a transformações morfológicas ao longo do estuário nos últimos 28 anos (entre as escalas Eventual e de Engenharia), houve neste período um acréscimo de 0,236 km² (3%) de área, não trazendo sigificativas mudanças para o estuário. Em relação a taxa de sedimentação, correspondente ao segundo bjetivo, a partir da coleta de 9 testemunhos, de até 1 m de profundidade e utilizando o radionuclídeo ²¹⁰Pb, ao longo do estuário, obteve-se uma taxa que variou de 0,33 cm/ano a 1 cm/ano (escalas entre Engenharia e Geológica) próximo a foz do estuário, e com uma rápida sedimentação percebida na margem leste do rio, onde encontram-se sedimentos mais recentes em relação a margem oeste. Em relação ao preenchimento, terceiro e último objetivo, a partir da amostragem de testemunhos de até 18 m de profundidade, utilzando o amostrador Rammkernsonden (RKS), foram gerados perfis e seções estratigráficas que ajudaram a entender o preenchimento do vale inciso do estuário do rio Coreaú e entender que trata-se de um estuário fluvio-marinho, preenchendo os vales formados no Grupo Barreiras nos últimos 10.000 anos antes do presente. Estas análises e resultados servirão como base para comparação com outros estuários, sejam fluviais, fluvio-marinhos ou marinhos, para entendermos melhor quais os possíveis eventos que dominaram a sedimentação ao longo da costa brasileira em diferentes escalas.

Mapeamento dos parâmetros florísticos e estruturais de floresta de mangue com dados LIDAR e SRTM

Este estudo apresenta a estimativa dos parâmetros florísticos e estruturais (determinação da espécie, altura, diâmetro a altura do Peito - DAP e biomassa) do mangue a partir de informações da superfície adquiridas remotamente com os sensores Laser Detection and Range (LIDAR), Shuttle Radar Topography Mission (SRTM) e ortofotos na Ilha dos Guarás, conjunto de arquipélagos localizado a 30 km da desembocadura do rio amazonas. Para esse trabalho foram utilizadas informações do SRTM, LIDAR e fotografias aéreas processadas e ortorretificadas durante dois sobrevôos realizados entre o mês de julho e agosto de 2011. Com a ortofoto foi feito o mapa do reconhecimento de unidades geobotânicas que delimitou apenas a classe mangue. Em seguida, foi realizada a correção da altura elipsoidal para a altura ortométrica, onde a nuvem de pontos foi interpolada pelo método vizinho mais próximo, gerando Modelo Digital de Elevação (MDE) LIDAR (full points) com RMSE de 0,88 cm e por meio de uma linguagem macro foi estatisticamente separadas as informações do último pulso da superfície, conhecido também por ground points. Em seguida, os dados foram interpolados pelo método de krigeagem que gerou o valor de Modelo Digital de Superfície (MDS), o qual foi subtraído do MDE. Com base no Modelo Digital de Vegetação (MDV) foram definidos os sítios de coleta e selecionadas as árvores ascendentes, intermediárias e emergentes, porte no qual foi medido o DAP e altura. No total foram coletadas 212 amostras individuais de mangue e para assegurar o nível de acurácia do conjunto coletado, foi realizado o cálculo de RMSE entre as alturas do LIDAR e Campo, que resultou em RMSE= 1,10 m. Os modelos escolhidos para calibração LIDAR e altura de campo foi do tipo linear, com R² = 91% e RMSE= 0,98 cm e para calibração da DAP e altura de campo foi escolhido o modelo Logarítmico R² = 74,1%. Nos resultados da calibração do SRTM o modelo logarítmico também foi o mais adequado para a relação entre altura média e SRTM com R² = 91% e RMSE de 2,2 m e DAP Médio e SRTM, com R² = 88% e RMSE 2,2 cm. A partir de um inventário foi realizada a estimativa da biomassa por espécie por meio das equações alométricas de Fromard e posteriormente os resultados foram espacializados em forma de mapas com alto nível de detalhamento oriundo das informações LIDAR e SRTM corrigido e ortofotos.

Modelagem de dinâmica do uso da terra e cobertura vegetal na região de Santarém, oeste do Pará

A região de Santarém, na última década, apresentou um aumento da área de agricultura de grãos, em especial, arroz, milho e soja. Na base da estrutura fundiária da região, essa dinâmica tem concorrido para a concentração fundiária devido à substituição da pequena propriedade pela grande propriedade capitalizada. Políticas de ordenamento territorial criaram um mosaico de unidades com regras específicas de uso da terra, são unidades de conservação e diferentes modalidades de projetos de assentamentos. Este trabalho tem dois objetivos: estudar o processo de transformação da paisagem após a introdução da agricultura capitalizada de grãos e construir cenários de futuro que analisem alternativas para conter o desflorestamento e a concentração fundiária em curso. O estudo foi realizado utilizando técnicas de sensoriamento remoto e geoprocessamento, a partir de imagens Landsat 5 TM dos anos de 1999, 2004 e 2007. Técnicas de modelagem dinâmica foram empregadas para explorar cenários de futuro (2015) considerando regras de uso do território. Os resultados obtidos mostram que, até 2004, a maior parte da agricultura mecanizada foi implantada em áreas onde anteriormente eram ocupadas pela agropecuária familiar, pastagem e capoeira. Após 2004, a sua expansão se deu, principalmente, sobre áreas de floresta, em especial, dentro de projetos de assentamento. A análise das transições de uso, em diferentes modalidades de assentamento, demonstra que as regras de uso do território estabelecidas por medidas de ordenamento territorial não têm sido seguidas em muitos casos. Este trabalho apresentou como principal contribuição metodológica a incorporação de questões institucionais relacionadas à estrutura fundiária na análise de transformação da paisagem e construção de cenários. Os resultados mostram que tal abordagem é essencial para entender os processos de transformação correntes à região em questão.

Análise dos efeitos térmicos de superfície na cidade de Belém-Pará-Brasil utilizando imagens de satélite

Utilizando dados observacionais de precipitação e temperatura do ar para o período de 1967 a 2008, do qual se elaborou a climatologia destas variáveis, e analisando as condições atmosféricas influenciadas pelo fenômeno ENOS, pretende-se identificar, caracterizar e analisar os efeitos térmicos espaciais na cidade de Belém-PA, a partir de dados meteorológicos de superfície e de imagens do sensor MODIS, o qual está disposto sobre a plataforma do satélite Aqua. Observou-se que, a média da precipitação anual foi de 2978,6 mm/ano, e que apresenta tendência de aumento ao longo dos anos, comportamento semelhante observou-se para a temperatura do ar. Em geral, os resultados mostram dois núcleos de maiores intensidade de temperaturas da superfície, um na cidade de Belém e outro na cidade vizinha, Ananindeua. Estes variam espacialmente e temporalmente de intensidade. Durante eventos de La Niña, o núcleo da ilha de calor fica localizado em bairros mais próximos a baía do Guajará, enquanto que durante eventos de El Niño estes bairros apresentam temperaturas mais amenas do que os bairros mais afastados dos corpos hídricos. Observou-se ainda que, a amplitude térmica da temperatura superficial entre áreas urbana e rural variam bastante, com a maior variação de 30,8°C e a menor de 16,8°C. Neste sentido, as maiores temperaturas da superfície foram observadas nos bairros com baixo NDVI, conseqüência de uma urbanização mais densa. As superfícies urbanas e as superfícies vegetadas apresentam relações de causa e efeito muito próximas, principalmente, durante o período menos chuvoso, isto pode ser percebido pela correlação que apresenta valor acima de 50%. Este estudo apresenta resultados que auxiliam no melhor entendimento do comportamento e dos efeitos térmicos espaciais e temporais na cidade de Belém, pois o uso de imagens do satélite é de fundamental importância para a identificação e caracterização das condições ambientais climáticas e ilhas de calor urbanas.