Univ Federal Ceara
Programa de Pós-graduação em Estatística e Experimentação Agropecuária; UFLA; brasil; Departamento de Ciências Exatas
DIAGNOSTIC TECHNIQUES OF LOCAL INFLUENCE IN SPATIAL ANALYSIS OF SOYBEAN YIELD
#geoestatistic#maximum likelihood#restrict maximum likelihood#COVARIANCE#MODELS#Agricultural Engineering
Mapeamento do índice de vegetação da diferença normalizada em lavoura de algodão; Mapping of the normalized difference vegetation index in cotton field
Motomiya, Anamari Viegas de Araujo; Molin, José Paulo; Motomiya, Wagner Rogerio; Rojo Baio, Fábio Henrique
#Gossypium hirsutum L.#geoestatística#sensoriamento remoto#variabilidade espacial#sensor óptico ativo#Gossypium hirsutum L.#geoestatistics#remote sensing#spatial variability#active optical sensor
Os dados de sensoriamento remoto em campo podem fornecer informações detalhadas sobre a variabilidade de parâmetros biofísicos ligados à produtividade em grandes áreas e apresentam potencial para o monitoramento destes parâmetros, ao longo de todo o ciclo de desenvolvimento da cultura. Este trabalho objetivou mapear a variabilidade espacial do índice de vegetação da diferença normalizada (NDVI) e seus componentes, em duas lavouras comerciais de algodão (Gossipium hirsutum L.), utilizando sensor óptico ativo, em nível terrestre. Os dados foram coletados utilizando-se sensor instalado em um pulverizador autopropelido agrícola. Um receptor GPS foi acoplado ao sensor, para a obtenção das coordenadas dos pontos de amostragem. As leituras foram realizadas em faixas espaçadas em 21,0 m, aproveitando-se as passadas do veículo no momento da pulverização de agroquímicos, e os dados submetidos à análise estatística clássica e geoestatística. Mapas de distribuição espacial das variáveis foram elaborados pela interpolação por krigagem. Observou-se maior variabilidade espacial do NDVI e da reflectância espectral da vegetação na região do infravermelho próximo (IVP) (880 nm) e do visível (590 nm) na lavoura com maior estresse fisiológico...
Uso de geoestatística para aprimorar o controle de pragas na citricultura.; Use of geoestatistical to improve the pests control in citrus plants.
Okumura, Miriam Harumi
Publicado em 24/10/2007 Português
#Citricultura#Citrus plants#Controle de pragas#Geoestatística#Geoestatistics#Mineração#Pest control#Tecnology transfer
A preocupação com o desenvolvimento de novas tecnologias que priorizem o menor impacto ambiental, é cada vez maior. Esta tendência também ocorre na agroindústria, onde a preocupação com a redução de uso de defensivos agrícolas já é grande. Os problemas fitossanitários sempre foram barreiras que impediram o desenvolvimento da produtividade da citricultura, apenas ácaro da leprose (Brevipalpus phoenicis) representa cerca de 80 % dos custos destinados ao controle de pragas. O desenvolvimento de novas tecnologias, que colaborem com o controle de pragas e doenças já existentes, e outras que podem surgir, são uma necessidade deste setor. Nesta pesquisa foram utilizadas ferramentas informatizadas, utilizadas com sucesso no setor mineiro, no estudo dos ácaro da falsa ferrugem (Phyllocoptruta oleiva), ácaro branco (Plyphagotarsonemus latus) e ácaro da leprose (Brevipalpus phoenicis) com o objetivo de aprimorar o planejamento e gerenciamento da aplicação localizada de defensivos agrícolas, permitindo a verificação das variações espaciais e temporais dos fatores relacionados ao controle dessas pragas que prejudicam a produção citrícola. Estas ferramentas possuem análises geoestatísticas, que são primordiais para a modelagem espacial de fenômenos naturais...
Estimativa da vulnerabilidade de aqüíferos utilizando sistemas de informação geográfica e geoestatística-UGRHI-PCJ; Aquifer vulnerabilty assessment using geographic information system and geoestatistics - UGRHI-PCJ
Russo, Alexandre de Silvio
#Aquifer#Aquíferos#GIS#SIG#Vulnerabilidade#Vulnerability
A estimativa da vulnerabilidade natural de aqüíferos é uma das primeiras tarefas a ser realizada para o planejamento, controle e prevenção da qualidade das águas subterrâneas e criação de perímetros de proteção para as mesmas. Nesse estudo o modelo GOD (G - groundwater confinement, O - overlay strata e D - depth to groundwater) (FOSTER et al. 1998) foi adotado para a estimativa da vulnerabilidade natural dos aqüíferos pertencentes à Unidade de Gerenciamento dos Recursos Hídricos do Estado de São Paulo das Bacias Hidrográficas dos Rios Piracicaba, Capivari e Jundiaí, sendo realizada uma extensa pesquisa sobre as características fisiográficas e ambientais da área de estudo, e critérios hidrogeológicos adotados para a estimativa de vulnerabilidade de aqüíferos. O modelo GOD utiliza como parâmetros o tipo de confinamento da água subterrânea, a litologia da zona vadosa e a profundidade do nível da água. Para a elaboração do respectivo mapa, foram adotados como base de dados o Mapa Geológico do Brasil ao Milionésimo (CPRM, 2004) e o Mapa de Águas Subterrâneas do Estado de São Paulo (DAEE, IG, IPT e CPRM, 2005). Apesar da notoriedade do modelo GOD e de sua extensa aplicação no Brasil e América Latina...
Krigagem com deriva externa aplicada à avaliação de recursos minerais de calcário e de minério laterítico; Kriging witch external drift applied to evaluation of mineral resources of limestone and lateritic ore
Fernandes, Jorge Augusto Basilio
#Applied geoestatistics#Avaliação de recursos minerais#Evaluation of mineral resources#Geoestatística aplicada#Geoestatística multivariada#Krigagem com deriva externa#Kriging with external drift#Multivariate geostatistics#Residuals variograms#Variograma de resíduos
Essa dissertação é um trabalho na área de geoestatística aplicada à avaliação de recursos minerais. Esta aplicação pode ser considerada o passo anterior à avaliação das reservas minerais. Os conceitos básicos são fundamentados na geoestatística atual, com amplo uso de recursos de informática. O objetivo principal foi de testar a técnica da krigagem com deriva externa na avaliação de recursos minerais, principalmente no que tange o ganho de resultado com a utilização de mais de uma variável, ou seja, verificar se o uso de variáveis auxiliares pode ser feito de maneira amigável com bom custo/benefício. Foram escolhidos dois depósitos, com geologias distintas entre si, um depósito de metal base com minério saprolítico e outro depósito de calcário. Em cada um foram selecionadas duas variáveis, uma variável principal e outra secundária, no depósito laterítico as variáveis foram 'SiO IND.2' e Fe e no depósito de calcário as variáveis foram CaO e 'SiO IND.2'. Foi elaborada a modelagem tridimensional dos depósitos no programa Datamine. Os modelos foram transferidos para o programa Isatis e serviram como base para a estimativa das variáveis. Essa modelagem foi necessária para que se pudesse dar aos blocos estimados a geometria do corpo do minério. A estimativa por krigagem com deriva externa foi comparada à krigagem ordinária...
Epidemiologia e controle de fuligem e sujeira de mosca em macieiras (Malus x domestica Borkh) no Sul do Brasil; Epidemiology and control of sooth blotch and flyspeck on apples (Malus x domestica Borkh.) in southern Brazil
#Maçã#Mosca#Epidemiologia#Doença de planta#Controle químico
A cultura da macieira no Brasil ocupa uma área de 35 mil hectares, concentrada nos estados do Sul do país. Apesar do alto patamar técnico atingido, a cultura sofre o ataque e perdas por inúmeras doenças, dentre as quais, incluem-se a Fuligem e a Sujeira de Mosca (F&SM). Classificadas como doenças secundárias, reduzem o valor de comercialização dos frutos in natura, pela formação de manchas escuras na epiderme dos frutos como resultado da colonização epífita por fungos na superfície das maçãs. No Brasil, dados referentes à biologia de F&SM restringem-se a testes de eficiência de fungicidas ou provêm de estudos com outros patossistemas. O objetivo deste trabalho foi o de elucidar aspectos epidemiológicos destas doenças, fornecendo as bases para o manejo integrado. Foram realizados estudos de disponibilidade, dinâmica e controle de inóculo, de análise de distribuição espacial, de modelagem do progresso temporal e uma avaliação de um modelo de previsão de F&SM na região Nordeste do Rio Grande do Sul. As pesquisas foram realizadas nos ciclos 2006/07 e 2007/08 em Vacaria-RS. Foi verificado que o inóculo de F&SM esteve disponível durante todo o período de formação dos frutos, sugerindo que ciclos secundários são responsáveis pelo incremento da doença. A incidência de F&SM apresentou relação positiva com a intensidade das chuvas (mm/chuva e mm/h)...
Dalchiavon, Flávio Carlos; Passos e Carvalho, Morel de
Tipo: Artigo de Revista Científica Formato: 541-552
#Precision agriculture#geoestatistics#Glycine max#spatial variability
The soybean is the crop most cultivated in Brazil, with great socioeconomic importance. In the agriculture year 2008/09 in Selviria County, Mato Grosso do Sul State, in the Brazilian Savannah, was analyzed the production components and the soybean yield cultivated in a Typic Acrustox on no-tillage. The main purpose objective was select among the production components number of pods per plant, number of grains per pod, number of grains per plant, mass of a thousand grains, mass of grains per plant and population of plants, which of the best linear and spatial correlation aiming explain the soybean yield variability. The irregular geostatistical grid was installed to collect of data, with 120 sampling points, in an area of 8.34 ha. The values of spatial dependence range to be utilized should be among 38.1 and 114.7 meters. The model of the adjusted semivariograma was predominantly the spherical. of the lineal and spatial point of view, the number of pods per plant and the mass of grains per plant they were correlated in a direct way with the soybean yield, demonstrating be the best components to esteem her.
de Assis Silvan, Samuel; de Souza Lima, Juliao Soares; Batista de Oliveira, Rone; Soares de Souza, Gustavo; de Assis e Silva, Michel
Fonte: Univ Federal Ceara Publicador: Univ Federal Ceara
#spatial variability#geoestatistics#soil conservation#Coffea canephora Pierre x froenher
This work aimed to study the space behavior of the water erosion in a red-yellow latosol. Then a study was developed in an area with colinon coffee cultivation in an Experimental Farm of Bananal do Norte of INCAPER in Cachoeiro de Itapemirim - ES. Soil samples were obtained from 0,0 to 0,20 m depth in an irregular grid with 109 samples. The analyzed variables were granulometric fractions, erodibility (K), natural erosion potential (PNE), soil loss (A) and erosion risk (RE). All the variables showed space dependency with moderate index of space dependency and similar standard of space distribution. The soil loss is related with the space distribution of the granulometric fractions.
Modelagem da variação horária da temperatura do ar em Petrolina, PE, e Botucatu, SP
Ramos, Clovis M. C.; da Silva, Alessandra F.; Sartori, Anderson A. da C.; Zimback, Celia R. L.; Bassoi, Luis H.
Tipo: Artigo de Revista Científica Formato: 959-965
#geoestatistics#temporal dependency#climatic element
The air temperature is a climatic element of interest for several agricultural related studies. Hence, daily air temperature from two locations in Brazil (Botucatu and Petrolina, with different climate classification), hourly recorded from January 2004 to December 2005, were analysed to verify if 24 h long air temperature variation can be predicted with 6 h interval records. Three data set were selected, with different first time record of air temperature, and analyses performed using geoestatistics to confirm the temporal dependency of records. Then, ordinary krigging was used to model air temperature variation within 1 h interval, and spherical model presented the best approach, with range of 11.6 h for Botucatu and 12.7 for Petrolina. For both locations, the time of air temperature recording was the most importnat factor on accuracy of estimation of air temperature variation. Records at 2, 8, 14 and 20 h and at 4, 10, 16 and 22 h presented the best approaches, but for records at 9, 15 and 21 h, the approach was worst.; A temperatura do ar é um elemento climático de interesse para diversos estudos relacionados à produção agrícola. A partir de dados de temperatura do ar de Botucatu, SP, e Petrolina, PE (diferentes classificações climáticas)...
De Mello, Gláucia; Bueno, Célia R. P.; Pereira, Gener T.
Tipo: Artigo de Revista Científica Formato: 315-322
#Geoestatistics#Soil conservation#Soil map#USLE (Universal Soil Loss Equation)
The objective of this work was to identify the spatial variability of the natural erosion potential, soil loss and erosion risk in two intensely cultivated areas, in order to assess the erosion occurrence patterns. The soil of the area located at Monte Alto, São Paulo state, was classified as Paleudalf (PVA) with moderately slope, with different managements. The soil of the area located at Jaboticabal, São Paulo state, was classified as Haplortox(LV) with gentle slope and cultivated with sugarcane. A irregular grid was imposed on the experimental areas. Soil samples were obtained from 0-0.2 m depth at each grid point: 88 samples in Monte Alto area (1465 ha) and 128 samples at Jaboticabal area (2597 ha). In order to obtain the values of the studied variables USLE was applied at each grid point. Descriptive statistics were calculated, and geoestatistical analyses were performed for defining semivariograms. Kriging techniques to develop map showing spatial patterns in variability of selected soil attributes were used. All variables showed spatial dependence. The PVA soil showed higher erosion risk due to the slope and atual management compared to the soil LV.
Silva, Samuel de Assis; Lima, Julião Soares de Souza; Souza, Gustavo Soares de; Oliveira, Rone Batista de
Tipo: Artigo de Revista Científica Formato: 75-81
#Geoestatistics#Handling areas#Interpolation methods
The obtaining of the correct space distribution for attributes of the soil is relevant in the agricultural planning, in what concerns to the installation and maintenance of the cultures. The objective of that work was to compare statistical interpolation methods (ordinary krigagem) and deterministic methods (inverse square distance) in the estimate of CTC and V% in a distrophic yellow-red Latossolo. The study was accomplished in the State of Experimental Hands on of the Instituto Capixaba de Pesquisa, Assistência Técnica e Extensão Rural (INCAPER), in an irregular grading with 109 points. The data were collected in the layer of 0 - 0,20 m in the projection of the cup of the plants, in the superior part of the slope. The performance of the interpolators was obtained and compared using the criterion of the medium mistake. The observations are dependent in space until a maximum reach of 14,1 m, considering the isotropy. IDW presented larger mistake in the estimate of the data; however its difference in relation to KRIG was small for both variables.
Tipo: Dissertação de Mestrado Formato: 127 f. : il.
#Cartografia#Imagens multiespectrais#Redes neurais (Computação)#Multispectral image#High spatial resolution#Artificial neural networks classification#Geoestatistics analysis
Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES); Pós-graduação em Ciências Cartográficas - FCT; Imagens multiespectrais de alta resolução espacial podem se constituir em uma fonte de dados adequada para o mapeamento de processos associados ao desenvolvimento de culturas agrícolas, como a detecção de plantas daninhas. A aerofotogrametria convencional e imagens de satélite de alta resolução espacial são alternativas para aquisição dessas imagens. Entretanto, devido ao custo elevado da aquisição destas imagens, tem sido desenvolvido, pelo Departamento de Cartografia da Faculdade de Ciências e Tecnologia da UNESP de Presidente Prudente, um Sistema de Sensoriamento Remoto Aerotransportado (SRA), capaz de oferecer resolução espacial sub-métrica. Este trabalho considerou a hipótese de que a partir de imagens adquiridas pelo Sistema é possível discriminar graus de infestação de plantas daninhas em culturas de café. Para investigar esta hipótese, foi realizado o mapeamento de plantas daninhas utilizando dois diferentes métodos: classificação de imagens multiespectrais (classificação por redes neurais artificiais - RNA) e análise geoestatística (krigagem por indicação com dados indiretos). Os mapas temáticos foram submetidos à análise da qualidade temática. A krigagem por indicação apresentou resultados suavizados e confusos...
Tipo: Dissertação de Mestrado Formato: xii, 67 f. : il.
#Cana-de-açúcar#Pragas - Controle biologico#Geologia - Métodos estatísticos#Distribuição de probabilidade binomial negativa#Liberação inundativa#Parasitóide#Biological contraol#Negative binomial distribution#Inundative liberation#Geostatistic#Parasitoid
Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq); Pós-graduação em Agronomia (Entomologia Agrícola) - FCAV; distribuição espacial do parasitismo, número de machos, fêmeas e total de adultos emergidos da geração F2 do parasitóide larval de Diatraea saccharalis (Fabricius), Cotesia flavipes (Cameron), foi estudada em um talhão comercial de 100x105 m de cana-de-açúcar da variedade CTC 3, com 7 meses, dividido em 100 parcelas de 105 m2, usando a metodologia de liberação de 4 pontos por hectare, com 1500 parasitóides espaçados 50 m entre eles. Em cada parcela foram colocados 10 colmos de cana com uma larva de D. saccharalis em seu interior fixados no solo para avaliação do parasitismo e número de insetos emergidos. Os índices de dispersão utilizados foram: razão variância/média (I), índice de Morisita (Iδ), coeficiente de Green (Cx) e expoente k da distribuição binomial negativa. A distribuição binomial negativa foi o modelo mais adequado para representar a distribuição de freqüência do parasitóide, já que a variância foi superior à média para todas as variáveis avaliadas. A análise dos índices de agregação mostrou que esse inseto se distribui de maneira agregada no campo. O uso da Geoestatística permitiu quantificar e mapear a distribuição de C. flavipes...
Variabilidade espacial de atributos do solo e da produtividade de soja
Beraldo, José Marcos Garrido
Tipo: Dissertação de Mestrado Formato: iii, 62 f. : il.
#Solos - Manejo#Produtividade agrícola#Soja#Geostatistics#Precision farming
Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES); Pós-graduação em Agronomia (Produção Vegetal) - FCAV; O presente trabalho teve como objetivo avaliar a distribuição espacial de atributos do solo antes e após a aplicação em doses variadas de calcário e fertilizante fosfatado e estudar as correlações entre os atributos químicos e físicos do solo com a produtividade da soja. Foram coletadas amostras de solo geo-referenciadas na intensidade de uma amostra por hectare na profundidade de 0-0,20 m em três talhões. Foram determinados: pH, teores de matéria orgânica, P, K, Ca, Mg, H+Al, textura do solo e calculado a CTC e o V%. O solo foi corrigido com calcário e fertilizado com adubo fosfatado, utilizando uma distribuidora dotada de um sistema para aplicação em doses variadas, procedendo-se a semeadura da soja. A colheita foi monitorada utilizando uma colhedora equipada com um monitor de produtividade. Após a colheita coletou-se amostras de solo nos mesmos pontos de amostragem. Os dados foram analisados com análise estatística descritiva e geoestatística. O manejo do solo através da aplicação variada de calcário e fósforo diminuiu a variabilidade espacial dos teores de P e os valores de V%. O P foi o atributo do solo que mais se correlacionou com a produtividade da soja. O mapa de produtividade possibilitou identificar onde as baixas e altas produtividades estão localizadas na área.; The objective of this work was to evaluate soil properties maps developed by different methods...
Metodos auxiliares para diagnostico da necessidade de subsolagem de solos agricolas
#Solos#Penetrometro#Solos#Engenharia
A subsolagem dos solos agrícolas é uma prática utilizada para se minimizar os efeitos causados pela sua compactação. Vem sendo geralmente realizada em todo o campo antes de ser implantada uma nova cultura. Por ser uma operação de preparo de solo que consome uma quantidade elevada de energia deve ser otimizada. Uma forma de otimizar essa operação é realiza-la apenas na área onde ocorre resistência mecânica igual ou superior à considerada crítica em função do tipo de cultura, solo, umidade, etc. Para avaliar os níveis de resistência mecânica à penetração e auxiliar na tomada de decisão quanto a necessidade de se proceder a operação de subsolagem foram propostas três metodologias. As metodologias propostas baseiam-se no volume de solo a ser mobilizado em função da resistência mecânica à penetração considerada crítica para o caso. O primeiro método, denominado probabilístico, utiliza-se de técnicas estatísticas para a definição do volume de solo a ser mobilizado. O segundo, denominado determinístico, utiliza-se do cálculo da porcentagem do volume de solo a ser mobilizado e sua distribuição espacial. Finalmente, o terceiro, denominado geoestatístico, utiliza-se de técnicas geoestatísticas para a determinação do volume de solo a ser mobilizado e da variabilidade espacial da resistência mecânica à penetração. Por meio da comparação entre as curvas de volume de solo a ser mobilizado...
Determinação de esquemas de amostragem para avaliação de propriedades fisicas do solo sob diferentes sistemas de preparo
#Fisica do solo#Agricultura
A metodologia atual de avaliação dos efeitos do preparo do solo (com ferramentas larga ou estreita), tem sido limitada e não tem conseguido capaz de detectar diferenças estatísticas significativas e nem tampouco orientar a coleta de dados quanto à forma e o número de amostras para avaliar propriedades físicas do solo. A utilização de uma nova proposta metodológica de avaliação, tendo a geoestatística mostrou que os diferentes tipos de implementos causaram modificações específicas nos principais parâmetros físicos do solo analisados. Todos os parâmetros analisados quando comparados por grupos de cinco e quarenta e cinco amostras registraram diferenças estatísticas ao nível de 5% de probabilidade. Quando relacionados ao esquema de amostragem, também apresentaram diferenças estatísticas significativas entre o modelo convencional de amostragem e o modelo proposto neste trabalho. A maioria dos parâmetros físicos do solo em suas respectivas profundidades, demonstraram não possuir dependência espacial e apresentaram efeito pepita puro. O tratamento com escarificador apresentou o maior número de modelos geoestatísticos ajustados. A profundidade de 0 a 200 mm apresentou o maior percentual de modelos geoestatísticos onde os semivariogramas mostraram dependência espacial dos dados.; The current methodology of evaluation the soil tillage effects by using wide and narrow tines has been limited and it has not been able to detect significant statistical differences and has not be able to guide the sample size requirements for assessing soil physical properties. The use of a new approach based on geoestatistics showed that different type of agricultural implements caused specific modifications in the main physical parameters. All studied parameters when compared by groups of five and forty five samples...
Molinero, Giulliano Renato
#Mecânica dos solos#Engenharia civil#Solos - análise#Solos - erosão
Tese (doutorado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Civil e Ambiental, 2007.; Este trabalho apresenta uma proposta metodológica para o desenvolvimento de modelagem para previsão de evolução de processos erosivos ao longo de rodovias. No desenvolvimento da pesquisa foram estudadas as rodovias BR-060, BR-080, DF-220, e DF-290. Para tanto, inicialmente foi desenvolvida uma modelagem determinística para a previsão pontual dos processos erosivos localizados às margens das rodovias, com base nas informações obtidas em ensaios com amostras de solos coletados em campo e cartas temáticas obtidas pelo projeto Radam Brasil. Na seqüência essas informações foram aprimoradas a partir da fotointerpretação de imagens de satélite de alta resolução - 0,60 m – Quick Bird aliadas a pesquisas de campo. Assim, as informações contidas nas cartas temáticas do projeto Radam Brasil foram melhor compartimentadas. A partir da modelagem determinística foi desenvolvida uma proposta de espacialização do modelo para trabalhos em escala regional. Para isso as variáveis pontuais (umidade higroscópica e massa específica aparente seca dos solos) foram espacializadas a partir de análises geoestatísticas...
Geoestatística aplicada às normais climatológicas de temperaturas médias compensadas no Brasil
Garcia, Thiago Junior Furtado
Fonte: Programa de Pós-graduação em Estatística e Experimentação Agropecuária; UFLA; brasil; Departamento de Ciências Exatas Publicador: Programa de Pós-graduação em Estatística e Experimentação Agropecuária; UFLA; brasil; Departamento de Ciências Exatas
#Temperatura#Temperature#Semivariograma#Semivariogram#Krigagem#Kriging#Mudanças climáticas#Climatic changes#Geoestatistics#climatological standard normals#Estatística
The study and the concern over the climate have been becoming more and more frequent. That because life in general is dependent upon favorable conditions of climate on the planet. That work studied the Average Temperature Normals in the period of 1931-1960 to 1961-1990 by utilizing of the geostatistical methodology (a branch of spatial statistics). It was aimed to study and understand then spatial distribution of the Normals in the Brazilian territory for each of those periods aiming to evaluate possible change of the means of temperature in the country. From the use of Geoestatistics, it was possible to know deterministic characteristics (trends) and spatial distribution characteristics such as: semivariances and semivariogram, fit model of semivariograms, anisotropy, model validation and study of the best spatial interpolation model, Kriging. Through Kriging, two interpolation maps for the average temperatures respective to the two periods were generated. On a third map, generated with the differences between the two first ones, it was found that the regions affected by changes of means of temperature. We conclude the increase of the average temperatures around 0.2 to 1.33oC. In the studied period (1931-1990), the world means increased about 0.55oC. Thus...