uff

arXiv · 02 de setembro de 2015

Análise de Padrão de Pontos Espaciais dos Fenômenos Aéreos Não Identificados na França

Pesquisadores da Toulouse School of Economics modelaram 1.969 registros de UAP (Fenômeno Aéreo Não Identificado) observados na França entre 1951 e 2013 como padrão de pontos espaciais. O modelo estatístico identificou correlação significativa entre ocorrências não explicadas e densidade populacional, proximidade de instalações nucleares e terrenos contaminados. Algumas regiões — como a fronteira belga, a ponta da Bretanha e a Picardia — apresentam concentrações que o modelo não consegue explicar.

Visão Geral

Este artigo, publicado no repositório arXiv em setembro de 2015 (arXiv:1509.00571v1), apresenta uma análise estatística espacial dos Fenômenos Aéreos Não Identificados (UAP — Unidentified Aerial Phenomena; em francês, Phénomène Aérospatial Non identifié) registrados na França metropolitana ao longo de 60 anos. Os autores são Thibault Laurent (GREMAQ/CNRS, Toulouse School of Economics), Christine Thomas-Agnan (GREMAQ, Toulouse School of Economics) e Michaël Vaillant (Meta-Connexions). [§ p.1]

Um UAP corresponde a um fenômeno que não encontra explicação legítima no estado atual do conhecimento — na maioria das vezes por falta de informação, mas também, em casos mais raros, por limitações científicas genuínas. [§ p.1]

1. Introdução

Sistema de registro e classificação na França

Na França, quando uma pessoa testemunha um UAP, ela pode reportá-lo à Gendarmerie. O testemunho é registrado em formulário detalhado contendo: data, hora, local, duração, orientação, forma, tamanho, trajetória e distância do observador ao fenômeno. A investigação é então encaminhada ao GEIPAN (Groupe d'Études et d'Informations sur les Phénomènes Aérospatiaux Non-identifiés, http://www.geipan.fr/), unidade da agência espacial francesa CNES, cuja missão principal é validar as informações fornecidas e determinar a natureza do UAP. [§ p.2]

O GEIPAN classifica cada UAP em quatro categorias:

Até o momento do estudo, 19,5% dos UAPs permanecem indeterminados após investigação. [§ p.2]

Nova estratégia de análise

Após mais de 50 anos de progresso limitado no campo, os autores propõem uma abordagem global — baseada em análise agregada, não em testemunhos isolados — partindo da seguinte premissa: fenômenos aeroespaciais, qualquer que seja sua natureza, originam-se de duas fontes: [§ p.2]

A estratégia adotada consiste em identificar variáveis ambientais que possam se correlacionar com a presença de UAPs, de modo a reduzir o espaço de hipóteses explicativas. [§ p.2]

Seleção de variáveis

Os autores consideraram duas categorias de variáveis: [§ p.3]

Covariáveis Caracterização ambiental Riscos sistêmicos
Antrópicas Densidade populacional Número de instalações nucleares*
Número de aeroportos Número de terrenos contaminados*
Percentual de zonas úmidas
Ambientais Percentual de florestas N/A**
Insolação

*Nota dos autores: estamos na faixa de riscos convencionalmente denominados "Nuclear, Biológico e Químico" (NBC). *Na faixa de riscos naturais não foram testados: riscos sísmicos, inundações, riscos hidrogeológicos, vulcânicos e incêndios florestais. Estes poderão ser incluídos em estudos futuros. [§ p.3]

Além disso, os autores incluíram os UAP As (fenômenos identificados) como covariável, para avaliar se variações no número de UAP Ds refletem resistências ou aberturas sociopsicológicas regionais, em vez de fenômenos reais. [§ p.3–4]

Escopo geográfico e temporal

O estudo abrange a França metropolitana de 1951 a 2013, totalizando 1.969 UAPs registrados. UAPs observados na Córsega, sobre o oceano ou o Mar Mediterrâneo foram excluídos para facilitar a análise do padrão de pontos em janela conexa. UAPs com imprecisão de localização superior a 20 km (centroide do convex hull distante mais de 20 km dos vértices) foram também excluídos, representando menos de 8,1% dos casos. [§ p.4]

A análise visual dos dados por departamento revela que a intensidade dos UAPs não é homogênea no espaço: os departamentos do norte da França e aqueles com grandes metrópoles (Bouches-du-Rhône, Gironde, Isère) apresentam contagens mais elevadas. A distribuição normalizada por 100.000 habitantes confirma que a densidade populacional, isoladamente, não explica a intensidade dos UAPs. [§ p.4]

2. Estatísticas Descritivas dos Fenômenos Aéreos Não Identificados

Distribuição temporal dos UAP Ds

A proporção de UAPs classificados como categoria D variou ao longo do tempo. Nos últimos dez anos do período analisado (2000–2010 e 2010–2020), essa proporção caiu significativamente — de 0,25 (anos 1980–1990) para 0,08 (2000–2010) e 0,03 (2010–2020) —, reflexo da maior disponibilidade de informação para casos recentes. Um teste χ² indica que a distribuição dos UAPs depende do período de observação (χ² = 233; hipótese nula de independência rejeitada com p-valor < 0,05). [§ p.6]

Intensidade suavizada por kernel dos UAP Ds

Os 381 UAP Ds observados foram tratados como realização de um padrão de pontos espaciais {x₁, …, xₙ} na janela de observação A (polígono da França metropolitana; área = 540.461 km²). A intensidade média é de 0,000705 casos por km². [§ p.7]

Um teste χ² de Aleatoriedade Espacial Completa (CSR — Complete Spatial Randomness) usando contagem em grade 6×6 quadrantes fornece χ² = 198 (p < 0,05), rejeitando a hipótese de CSR. O envelope de simulação para a função K de Ripley (100 simulações) também rejeita CSR. [§ p.7]

A intensidade estimada λ*(u) em um ponto u é calculada por:

λ*(u) = Σᵢ e(xᵢ) k(xᵢ − u)

onde k é o kernel Gaussiano de suavização e e(xᵢ) é o fator de correção de borda. O parâmetro de largura de banda σ foi fixado em 20 km (bw.diggle() recomenda 18 km; bw.scott() sugere 77–100 km). [§ p.7]

A análise do mapa de intensidade revela vários agrupamentos locais com alta concentração de UAPs: [§ p.8]

3. Descrição das Covariáveis

Todas as covariáveis provêm de serviços estatísticos oficiais franceses. Os dados foram disponibilizados em diferentes resoluções espaciais (municípios, padrão de pontos, pixels) e sistemas de coordenadas distintos, exigindo processamento de harmonização. [§ p.8]

3.1 Densidade Populacional

Dados do INSEE (Censo de 1990), expressos em habitantes/km² para cada um dos 36.208 municípios da França metropolitana. O suavizador de Nadaraya-Watson é aplicado para estimar a densidade em cada pixel:

g(u) = Σᵢ k(u − xᵢ)vᵢ / Σᵢ k(u − xᵢ)

com σ = 5 km. [§ p.9]

51,5% dos UAP Ds ocorrem em zonas com densidade superior a 100 hab./km², embora essas zonas representem apenas 21,8% do território francês. Um teste de Kolmogorov-Smirnov de adequação ao modelo de Poisson (estatística D = 0,3023; p < 0,05) confirma que a densidade populacional é significativa para explicar a intensidade dos UAPs. [§ p.9]

Habitantes/km² ≤50 ]50;100] ]100;500] ]500;2000] ]2000;11224]
% de UAP Ds (de 373) 29,5 19,0 35,1 12,3 4,1
% de Pixels (de 9480) 54,2 24,0 19,1 2,3 0,4

[§ p.10]

3.2 Instalações Nucleares

A França possui 102 sítios com risco nuclear, incluindo usinas nucleares, resíduos de mineração de urânio, armazenamento e descarte de rejeitos radioativos. Para cada pixel, foi calculado o número de instalações nucleares num raio de 20 km. [§ p.11]

12,3% dos 9.480 pixels possuem ao menos uma instalação nuclear nas proximidades. Entre os pixels que contêm UAP Ds, esse percentual sobe para 21,2%, indicando correlação potencial. O teste espacial de Kolmogorov-Smirnov rejeita CSR (D = 0,1299; p < 0,05). [§ p.11]

Instalações nucleares 0 ]1;5] ]5;10] ]10;13]
% de UAP Ds (de 373) 78,8 19,0 1,6 0,6
% de pixels (de 9480) 87,7 11,9 0,3 0,1

[§ p.11]

3.3 Terrenos Contaminados

Um sítio contaminado é um local com risco permanente ou potencial à saúde humana ou ao meio ambiente, decorrente de poluição de atividade industrial anterior ou atual. Ao final de 2007, havia 3.985 sítios contaminados na janela A, para os quais o governo francês adotou medidas corretivas (dados do Ministério da Ecologia, base BASOL). [§ p.11]

49,3% dos UAP Ds ocorrem em pixels com ao menos um sítio contaminado num raio de 5 km, contra 19,7% para todos os pixels da janela A. O teste de Kolmogorov-Smirnov rejeita CSR (D = 0,3154; p < 0,05). [§ p.13]

Sítios contaminados 0 ]1;5] ]5;10] ]10;20] ]20;58]
% de UAP Ds (de 367) 50,7 31,9 7,0 6,7 3,8
% de pixels (de 9480) 80,3 16,9 1,9 0,7 0,2

[§ p.13]

3.4 Outras Covariáveis

As seguintes variáveis adicionais foram consideradas: [§ p.14]

Todos os testes de Kolmogorov-Smirnov individuais rejeitam a hipótese CSR para cada covariável (ver Tabela 6): [§ p.14–15]

Covariável Zonas Úmidas Aeroportos Floresta Sol UAP A
Estatística D 0,1902 0,2569 0,089 0,1685 0,1782
p-valor 2,155×10⁻¹² < 10⁻¹⁶ 0,004782 8,083×10⁻¹⁰ 6,198×10⁻¹¹

4. Modelo

Os autores ajustam um modelo de Poisson espacial não homogêneo usando o pacote spatstat do R, com a função ppm(). A intensidade é modelada como função log-linear das covariáveis Z: [§ p.15]

λ(u) = e^(Z(u)β)

onde β são os parâmetros [β₀ ··· β₈] associados às covariáveis:

O modelo foi ajustado pelo método da máxima pseudoverossimilhança, usando os 381 UAPs como pontos de dados e os 9.480 pixels como pontos auxiliares (dummy points). [§ p.15]

Seleção do modelo final

As variáveis foram removidas progressivamente pela menor significância estatística: primeiro uapA (p = 0,41), depois log(sun) (p = 0,35), forests (p = 0,22), aero (p = 0,11) e wetlands (p = 0,12). [§ p.15]

Modelo final retido (Tabela 8):

Coeficiente Desvio padrão p-valor
(Intercepto) −9,67973 0,21282 < 10⁻¹⁶
log(pop) 0,54115 0,04684 < 10⁻¹⁶
conta 2,39716 0,86181 0,00542
nuclear 247,07907 64,55783 0,00013

[§ p.16]

Todos os coeficientes βⱼ são positivos, indicando que maiores valores de cada covariável estão associados a maior intensidade de UAP Ds.

Exemplo numérico (pixel u*)

Para o pixel de referência u*, os valores estimados das covariáveis são: [§ p.17]

pop conta nuclear
zⱼ(u*) 5,521 0,036 0,00062
zⱼ(u*) × βⱼ 2,99 0,087 0,15
e^(zⱼ(u*)×βⱼ) 19,83 1,09 1,17

A intensidade estimada em u* é:

λ̂(u*) = e^(−9,679) × 19,83 × 1,09 × 1,17 ≈ 0,00158

Análise de Resíduos

Os resíduos brutos são definidos como: s(u) = λ*(u) − λ̂(u). Para u*: s(u*) ≈ 0,00168 − 0,00158 = 1,05 × 10⁻⁴. [§ p.17]

O mapa dos resíduos de Pearson (s(u)/λ̂(u)) revela que as áreas de 5% maior subestimação pelo modelo (em vermelho) correspondem a: [§ p.17]

As áreas de 5% maior superestimação correspondem a regiões de alta densidade populacional sem muitos UAP As próximos, como as metrópoles de Lyon, Bordeaux e Toulouse. A região de Paris é um caso atípico: há muitos UAP Ds na área, e o agrupamento é delimitado por uma zona subestimada (a cidade de Paris) e uma zona superestimada (subúrbios a oeste, com menor densidade). [§ p.17]

5. Conclusão

Quando um fenômeno é classificado como categoria D, é extremamente difícil para o GEIPAN estabelecer convincentemente sua existência objetiva, pois a opinião pública sobre o tema é frequentemente marcada por ceticismo. O objetivo deste estudo não foi propor uma nova explicação para os UAP Ds, mas determinar se existe uma ligação entre suas localizações e as covariáveis fornecidas pelos serviços estatísticos franceses, abrindo novas vias de pesquisa. [§ p.18]

Os resultados mostram que a localização dos UAP Ds pode ser parcialmente explicada por covariáveis antrópicas:

Essas correlações podem ser resultado de atividade endógena emergente ou de atividade exógena. Uma hipótese aberta é que esses sítios sensíveis possam ser locais de interesse em razão de sua conexão com questões ambientais. [§ p.18]

A análise dos resíduos revela agrupamentos de UAP Ds ainda não explicados pelo modelo. Uma nova frente de pesquisa seria incluir covariáveis que reflitam o nível de escolaridade da população, bem como dados de sociologia e psicologia sobre os testemunhos e práticas investigativas. [§ p.18]

Estatísticas Descritivas das Covariáveis

Mín. 1º Quartil Mediana Média 3º Quartil Máx.
população 2,18 24,66 45,28 102,80 89,12 10.380
terrenos contaminados 0 0 0 0,01 0,01 0,67
instalações nucleares 0 0 0 0,20 0 13
zonas úmidas 0 0,03 0,22 0,93 0,87 80,48
sol 998 1.144 1.227 1.241 1.323 1.581
aeroportos 0 0 0 0,08 0,01 14,59
floresta 0,09 11,40 24,28 31,48 45,94 99,25

[§ p.20]

Glossário

UAP
Fenômeno Aéreo Não Identificado (Unidentified Aerial Phenomena); no contexto francês, *Phénomène Aérospatial Non identifié*
GEIPAN
Groupe d'Études et d'Informations sur les Phénomènes Aérospatiaux Non-identifiés — unidade da agência espacial francesa CNES responsável pela investigação de UAPs
CSR
Complete Spatial Randomness (Aleatoriedade Espacial Completa) — hipótese de que eventos ocorrem de forma completamente aleatória no espaço, usada como referência nos testes estatísticos
Kernel suavizado
Método estatístico não paramétrico que estima a intensidade espacial de um processo de pontos usando uma função de ponderação Gaussiana (kernel), com parâmetro de largura de banda σ
Processo de Poisson não homogêneo
Modelo estatístico em que a intensidade de ocorrências varia no espaço em função de covariáveis, ao contrário do processo homogêneo onde a intensidade é constante
Máxima pseudoverossimilhança
Método de estimação de parâmetros para processos de pontos espaciais que aproxima a função de verossimilhança completa, proposto por Baddeley e Turner (2000)
Função K de Ripley
Estatística descritiva para padrões de pontos espaciais que mede a dependência espacial entre pontos a diferentes distâncias
Resíduos de Pearson
Medida de adequação do modelo calculada como s(u)/λ̂(u), onde s(u) é o resíduo bruto entre intensidade observada e estimada em cada pixel
Suavizador de Nadaraya-Watson
Estimador não paramétrico de regressão que calcula a média ponderada dos valores observados, com pesos definidos por um kernel Gaussiano
spatstat
Pacote do software R para análise de padrões de pontos espaciais, desenvolvido por Baddeley e Turner (2005), utilizado em todo o estudo
NBC
Nuclear, Biológico e Químico — categoria de riscos sistêmicos considerada no estudo ao incluir instalações nucleares e terrenos contaminados como covariáveis
Corine Land Cover
Base de dados europeia de uso e cobertura do solo, utilizada no estudo para extrair variáveis de zonas úmidas, florestas e cobertura territorial

Perguntas frequentes

Quantos UAPs foram analisados no estudo e qual o período coberto?
O estudo analisou 1.969 UAPs observados na França metropolitana entre 1951 e 2013. Destes, 381 foram classificados como UAP D (não identificados), que constituem o objeto principal da análise estatística.
Quais variáveis estatísticas o modelo identificou como significativas para explicar a distribuição dos UAPs?
O modelo final reteve três variáveis significativas: densidade populacional (p < 10⁻¹⁶), intensidade de terrenos contaminados (p = 0,00542) e intensidade de instalações nucleares (p = 0,00013). Variáveis como floresta, aeroportos, insolação, zonas úmidas e UAP As foram removidas por não atingir significância estatística.
Como o GEIPAN classifica os UAPs registrados na França?
O GEIPAN utiliza quatro categorias: A (caso explicado inequivocamente, 12%), B (caso provavelmente identificado, 29,5%), C (não identificável por falta de dados, 39%) e D (não identificável, 19,5%).
Quais regiões da França apresentam concentrações de UAPs não explicadas pelo modelo?
A análise de resíduos indica subestimação significativa nas seguintes áreas: fronteira belga, ponta da Bretanha, partes do Sudeste, regiões de Picardia e Alta Normandia, e os departamentos do Loiret e Corrèze.
O estudo conclui que os UAPs são de origem extraterrestre?
Não. O objetivo declarado do estudo não é explicar a natureza dos UAPs, mas identificar correlações entre suas localizações e variáveis ambientais e antrópicas, para abrir novas vias de pesquisa. O documento não faz afirmações sobre origem extraterrestre.
Qual a correlação identificada entre instalações nucleares e UAPs não identificados?
O modelo indica correlação positiva e estatisticamente significativa entre a proximidade de instalações nucleares e a ocorrência de UAP Ds (p-valor = 0,00013). O coeficiente estimado para essa variável é 247,08. Os autores apontam que essa relação, 'há muito tempo suspeita', é medida quantitativamente pela primeira vez neste estudo.
A proporção de UAPs não identificados mudou ao longo do tempo?
Sim. A proporção de UAP Ds em relação ao total caiu significativamente: era 0,32 nos anos 1950–1960 e chegou a 0,03 nos anos 2010–2020. Um teste χ² (valor = 233) rejeita a independência entre período e classificação (p < 0,05).

Entidades citadas

Documentos relacionados