Um tremor de magnitude minúscula sacode o solo a dez metros de profundidade, perto do polo sul da Lua. A onda se espalha em todas as direções pelo regolito solto, viaja a pouco mais de cem metros por segundo e, ao encontrar o piso de uma cratera de cem metros de diâmetro que jamais viu o Sol, faz algo que nunca faria em terreno seco: ela volta. Mais de noventa por cento da energia é rebatida para trás, como se a onda tivesse batido em uma parede de concreto enterrada. Do outro lado da cratera, um sismômetro plantado na escuridão perpétua registra quase nada. Um vazio no sinal. E esse vazio é exatamente a assinatura do tesouro que meio século de exploração lunar tenta localizar com precisão.

O cenário acima não é ficção nem exercício de imaginação livre. Ele foi calculado, célula por célula, em um volume simulado de oito quilômetros por oito quilômetros por um quilômetro de profundidade, e publicado em 31 de julho de 2026 na revista Science Advances por Harrison P. Lisabeth, do Laboratório Nacional Lawrence Berkeley, Nicholas Schmerr, da Universidade de Maryland, e Matthew Siegler, da Universidade do Havaí. O trabalho responde a uma pergunta que a indústria petrolífera terrestre resolveu há décadas, mas que ninguém havia formulado seriamente para outro mundo: como seria o som do gelo lunar? Não o som audível, que o vácuo proíbe, e sim o modo como as ondas elásticas atravessam, refletem, contornam e se despolarizam ao encontrar água congelada enterrada no solo da Lua.

A resposta importa porque a água mudou de status na exploração espacial. Ela deixou de ser curiosidade científica para virar infraestrutura. Fora da Terra, os recursos energéticos que sustentam uma base permanente não são hidrocarbonetos, mas hidrogênio, oxigênio e outros voláteis — e todos eles, no polo sul lunar, estão presos na mesma molécula. Um quilo de gelo lunar extraído e eletrolisado vira propelente para retornar à órbita, oxigênio para respirar, água para beber e blindagem contra radiação. É por isso que tanto o programa Artemis quanto os planos chineses de base lunar convergiram para a mesma faixa de terreno acidentado nos oitenta e poucos graus de latitude sul, onde o Sol nunca sobe muito acima do horizonte e o fundo de certas crateras permanece congelado desde antes de existirem oceanos na Terra.

O problema é que saber que existe água não é o mesmo que saber onde cavar. As evidências acumuladas desde 2008 formam uma cadeia sólida, porém indireta. Naquele ano, análises de vidros vulcânicos trazidos pelas missões Apollo revelaram água aprisionada no interior lunar. Em 2009, o impacto controlado da sonda LCROSS contra o piso da cratera Cabeus levantou uma pluma de detritos onde havia gelo em quantidade mensurável. O mapeamento de hidrogênio pelas regiões polares mostrou concentrações compatíveis com depósitos de gelo. Em 2021, o observatório aerotransportado SOFIA detectou água molecular na porção iluminada da superfície. Dados espectroscópicos recentes indicam algo entre 100 e 400 miligramas de água por grama de regolito, valor muito acima do que se supunha uma geração atrás.

Toda essa constelação de medições compartilha uma limitação incômoda: sensoriamento remoto enxerga a superfície e os primeiros centímetros logo abaixo dela. Nêutrons, radar e espectroscopia infravermelha estabelecem a presença de hidrogênio e de água, mas restringem muito mal a abundância em profundidade. Para uma companhia de mineração terrestre, esse nível de informação equivale a saber que existe petróleo em algum lugar de uma bacia sedimentar de mil quilômetros quadrados. É o suficiente para justificar investimento em prospecção, e insuficiente para escolher onde posicionar a broca. O desafio prático da próxima década será identificar depósitos que sejam razoavelmente exploráveis em escala industrial, com volume, continuidade e acessibilidade conhecidos antes que o primeiro equipamento pesado saia da Terra.

A geofísica terrestre resolve esse tipo de pergunta com sísmica. Dispara-se uma fonte controlada, registra-se o campo de ondas em uma rede de sensores e reconstrói-se a estrutura do subsolo a partir dos tempos de chegada, das amplitudes e das conversões entre tipos de onda. A técnica localizou reservatórios de hidrocarbonetos, mapeou aquíferos, monitorou permafrost em degelo e identificou câmaras magmáticas. Sua transposição para a Lua parecia natural — e permaneceu, por décadas, sem base quantitativa, porque ninguém havia modelado explicitamente as consequências sísmicas de distribuições macroscópicas de gelo dentro do subsolo raso lunar. É essa lacuna que o novo trabalho fecha, costurando quatro peças que costumavam viver em comunidades científicas separadas: microtomografia de raios X em síncrotron, teoria de meios efetivos, modelagem térmica de terreno real e simulação tridimensional de propagação de ondas.

Há um motivo adicional para desconfiar de extrapolações ingênuas da experiência terrestre. O regolito lunar se comporta de forma acústica peculiar, sem paralelo em materiais comuns do nosso planeta. As medições feitas em amostras trazidas pela Apollo revelaram atenuação sísmica anomalamente baixa: o fator de qualidade Q, que mede o quanto um material preserva a energia de uma onda em vez de convertê-la em calor, chegou a valores iguais ou superiores a 3.100 na rocha lunar 70215. Comparado a qualquer sedimento seco terrestre, isso é extraordinário. A explicação aceita envolve adesão eletrostática notável entre partículas no ambiente seco da Lua, mecanismo que desaparece assim que voláteis adsorvidos recobrem as superfícies dos grãos e dissipam essas forças. O solo lunar também é altamente inelástico e fortemente dependente da pressão de confinamento.

A consequência conceitual é elegante: se voláteis controlam o comportamento acústico do regolito, então o regolito é, por natureza, um detector de voláteis. Basta aprender a ler o que ele diz. E ler significa, primeiro, escrever a relação matemática entre quanto gelo existe dentro dos poros e quão rápido uma onda elástica atravessa aquele volume. Essa relação não é única nem óbvia, porque depende da geometria com que o gelo se distribui entre os grãos. Gelo que apenas cimenta os pontos de contato entre partículas rígidas produz um efeito mecânico completamente diferente de gelo que preenche cavidades isoladas ou de gelo que envolve grãos flutuantes em uma matriz congelada contínua.

Para tratar esse problema, a física de rochas usa modelos de meio efetivo, que estimam as propriedades elásticas de um agregado complexo a partir das propriedades de seus componentes e das frações de volume de cada um. Dois limites teóricos delimitam o espaço de possibilidades. O modelo de Voigt corresponde à condição de deformação uniforme, na qual todas as fases sofrem a mesma deformação, e fornece o limite superior de rigidez — aproximação razoável quando o gelo é escasso e o empacotamento dos grãos sustenta a carga. O modelo de Reuss corresponde à condição de tensão uniforme, com todas as fases submetidas à mesma tensão, e fornece o limite inferior — situação plausível quando o gelo é abundante ou se concentra em lentes nas quais os grãos ficam suspensos em uma matriz congelada. A média entre os dois recebe o nome de Voigt-Reuss-Hill.

Uma analogia doméstica ajuda a fixar a diferença. Imagine dois blocos de materiais distintos comprimidos por uma prensa. Se eles estão lado a lado, ambos encolhem igualmente e o conjunto é tão rígido quanto a média das rigidezes: é o caso Voigt. Se estão empilhados um sobre o outro, a mesma força atravessa os dois em sequência, o mais macio cede primeiro e o conjunto é tão frouxo quanto a média das flexibilidades: é o caso Reuss. Materiais granulares naturais, como o regolito, situam-se em algum ponto intermediário, porque contêm arranjos dos dois tipos misturados em escalas microscópicas. Encontrar esse ponto intermediário com base física, e não por chute, é a diferença entre um modelo confiável e uma curva bonita sem significado.

Um refinamento adicional vem dos limites de Hashin-Shtrikman, formulados em 1963. Sob certas hipóteses de isotropia, eles definem as fronteiras teóricas mais estreitas possíveis para os módulos elásticos efetivos de um compósito de duas fases, dependendo apenas dos módulos dos constituintes e de suas frações volumétricas. O valor real do compósito precisa cair dentro dessa faixa, e a média entre os limites superior e inferior costuma fornecer boa estimativa das variações verdadeiras. É um instrumento consagrado na modelagem de sedimentos complexos e, particularmente relevante para o caso lunar, já foi aplicado com sucesso a permafrost não consolidado e a permafrost salino saturado, materiais em que gelo e grãos convivem em geometrias variadas.

Escolhido o arcabouço teórico, restava alimentá-lo com números reais. Os módulos elásticos de um material podem ser calculados a partir de sua densidade e das velocidades das ondas acústicas que o atravessam. O módulo de cisalhamento, que mede a resistência à deformação por tensões tangenciais, é o produto da densidade pelo quadrado da velocidade da onda S. O módulo de compressibilidade, que mede a resistência a mudanças de volume sob pressão hidrostática, sai da combinação entre as velocidades das ondas P e S. Medir as duas velocidades em laboratório, portanto, entrega diretamente os dois números que o modelo precisa.

O experimento foi montado em uma coluna cilíndrica de acrílico com 107,5 milímetros de diâmetro interno e 150 milímetros de altura, equipada com sensores ultrassônicos e de temperatura. Os transdutores foram posicionados em lados opostos, de modo que o sinal atravessasse o diâmetro da amostra. Um gerador de pulsos de onda quadrada forneceu excitação de 100 a 400 volts com frequência central de 1 megahertz, e os sinais transmitidos foram capturados em osciloscópio digital. Para cada amostra, o material foi pesado, despejado no cilindro, agitado por batidas para assentar os grãos, compactado e medido em altura. Na amostra saturada, água destilada foi despejada até o topo do empacotamento e o conjunto foi novamente agitado para expulsar bolhas de ar.

Dois experimentos paralelos foram conduzidos: um com a amostra totalmente saturada de água, levada a −10 °C, e outro com uma versão miniaturizada da célula contendo material nominalmente seco, levado a −40 °C. Essa assimetria de temperaturas não é descuido, e sim limite físico dos freezers disponíveis — o maior deles alcança temperaturas menos baixas, mas oferece estabilidade térmica superior. A amostra seca não foi levada ao forno e nenhuma tentativa foi feita para expulsar a água adsorvida, de modo que as superfícies dos grãos permaneceram hidratadas pela umidade atmosférica ambiente. As medições ultrassônicas foram tomadas durante o aquecimento lento das amostras, com pelo menos 24 horas de estabilização em cada patamar de temperatura.

Da amostra congelada saiu uma densidade de 2,103 gramas por centímetro cúbico e, das velocidades medidas nessa condição, um módulo de cisalhamento de 7,16 gigapascals e um módulo de compressibilidade de 22,19 gigapascals. Para o extremo oposto, o regolito sem gelo, os autores recorreram a velocidades ultrassônicas do conjunto de dados de pós de rocha da Apollo 17, obtidas sob condição hidrostática de 1 bar, que aproximam bem o ambiente de laboratório nominalmente seco. Essas amostras têm composição semelhante à do simulante usado no estudo, com vidro basáltico e plagioclásio. Combinadas a uma densidade de regolito seco de 1,510 grama por centímetro cúbico, elas devolvem módulo de cisalhamento de 0,84 gigapascal e módulo de compressibilidade de 4,41 gigapascals.

Os dois pares de números contam a história inteira em miniatura. Encher os poros de gelo multiplica a rigidez ao cisalhamento por um fator próximo de oito e meio, e a resistência à compressão por cerca de cinco. Um material que se comportava como areia fofa passa a se comportar como rocha. Nenhuma outra transformação disponível no ambiente lunar produz salto comparável nas propriedades mecânicas do solo, e é justamente essa magnitude que torna a sísmica promissora: contrastes assim não passam despercebidos por instrumento nenhum. Quando os autores aplicam os diferentes modelos de meio efetivo entre esses dois extremos, a dispersão entre as curvas previstas é ampla, o que reforça a exigência de escolher o modelo com base na microestrutura real, e não por conveniência matemática.

Determinar a microestrutura real exigiu ir ao Advanced Light Source, o síncrotron do Laboratório de Berkeley, e fazer microtomografia de raios X in situ do material experimental na linha de luz 8.3.2. As amostras foram varridas com luz monocromática de 32 quiloeletronvolts, através de uma cadeia óptica com cintilador de LuAG dopado com cério, objetivas Mitutoyo de 2 e 10 aumentos com longa distância de trabalho e abertura numérica de 0,055, e detector sCMOS de 2.560 por 2.160 pixels. O arranjo entregou tamanhos de pixel de 3,5 e 0,65 mícron, com campos de visão laterais de 5 e 0,9 milímetro. Na prática, isso significa enxergar cada grão individual e cada filete de gelo entre grãos em três dimensões, sem precisar cortar a amostra.

O material analisado foi o simulante JSC-1A, extraído de um cone de cinzas no Arizona e composto por cinza basáltica grossa, predominantemente vítrea, com alguma fração cristalina. Preparar uma amostra gelada realista é mais delicado do que parece: não basta molhar e congelar. Os pesquisadores usaram um protocolo de sinterização térmica que reproduz melhor a microestrutura esperada, submetendo o material a cinco ciclos térmicos até atingir equilíbrio textural, com as imagens capturadas a −40 °C. Sem esse cuidado, o gelo se organizaria de maneira artificial, e todo o edifício de modelagem construído sobre a microestrutura ficaria comprometido na base.

As imagens mostraram empacotamentos com porosidade total de aproximadamente 37%, valor próximo ao medido para o regolito lunar abaixo de cinco metros de profundidade pelo radar de penetração da sonda chinesa Chang’e-3. A maior parte do espaço vazio corresponde a poros grandes entre grãos, mas os próprios grãos também são porosos, com cavidades da ordem de alguns mícrons, desconectadas entre si. Essas cavidades internas provavelmente não funcionam como canais de fluxo, embora possam servir de esconderijo para voláteis aprisionados. O espaço poroso principal aparece parcialmente preenchido por gelo, que ocupa tanto os contornos entre grãos quanto o interior dos poros maiores.

A distribuição do gelo dentro da amostra é heterogênea, com bolsões completamente secos convivendo com regiões em que a fração local de gelo chega a 90%. Essa granulosidade importa: um depósito lunar não será uniforme, e a sísmica registrará a média ponderada de um mosaico. A análise estatística das imagens segmentadas indicou que grãos, gelo e poros compartilham distribuições de razão de aspecto amplas e largamente sobrepostas, com a maioria dos objetos agrupada em valores moderados, entre 0,45 e 0,55, sem evidência de população alongada ou direcionalmente enviesada. Os histogramas de orientação são aproximadamente planos ao longo de todo o intervalo angular, sem alinhamento preferencial persistente.

Traduzindo: o material não tem trama. Não existe orientação privilegiada, camadas achatadas nem segregação espacial entre as fases. A microestrutura é isotrópica e relativamente homogênea entre grãos, gelo e espaço vazio. Essa conclusão, aparentemente técnica, é a chave que autoriza o passo seguinte, porque os limites de Hashin-Shtrikman pressupõem justamente isotropia. A geometria observada — contatos multifásicos entre grãos e geometrias de contato diversas — corresponde ao que a média de Hashin-Shtrikman descreve melhor, e foi essa a escolha adotada para o modelo de velocidades. Persiste incerteza considerável sobre a textura verdadeira do gelo no regolito lunar, mas as imagens fornecem uma estimativa de primeira ordem sobre como o gelo pode existir dentro de empacotamentos de grãos.

A segunda perna do trabalho trata de onde o gelo pode estar, e essa pergunta é térmica antes de ser sísmica. A equipe adaptou uma família de modelos de estabilidade de gelo desenvolvida originalmente para interpretar dados do radiômetro Diviner, depois aplicada a Mercúrio, à Lua e a possíveis sítios de pouso lunares. O método usa traçado de raios para monitorar as trocas de radiação entre facetas triangulares construídas sobre a topografia real do terreno, o que permite calcular quanto de luz solar direta, de luz refletida pelas paredes vizinhas e de radiação infravermelha reemitida chega a cada pedaço de chão. Cada faceta é sustentada por um modelo térmico de elementos finitos com cem camadas, representando o primeiro metro e meio da superfície.

Abaixo desse metro e meio, as variações diurnas de temperatura praticamente desaparecem, e o modelo trata a condição inferior como fluxo de calor constante, adotado em 15 miliwatts por metro quadrado. A condução lateral de calor é desprezível na camada superficial, mas passa a importar em profundidade, à medida que a densidade cresce e a condutividade térmica sobe de valores da ordem de um milésimo de watt por metro por kelvin para 0,1 watt por metro por kelvin ou mais. Para estender o modelo além da camada superficial, os autores aplicaram a temperatura calculada em 1,5 metro como condição de contorno de uma malha tridimensional de elementos finitos com resolução aproximada de um metro, prolongando o cálculo até 500 metros de profundidade.

O alvo escolhido foi a região de Mons Mouton, um platô elevado próximo ao polo sul lunar, com o modelo centrado nas coordenadas 85,4096 graus sul e 31,1630 graus leste em projeção estereográfica polar. A escolha não é aleatória: Mons Mouton é o sítio de pouso previsto para o veículo VIPER, sigla em inglês para Rover de Exploração Polar Investigador de Voláteis, projetado especificamente para procurar gelo com uma broca de um metro. A região abriga abundância de regiões permanentemente sombreadas, ou permanently shadowed regions, áreas cujo fundo jamais recebe luz solar direta por causa da inclinação axial mínima da Lua combinada com a topografia local. São candidatas naturais a reter gelo por tempo geológico.

Os resultados térmicos justificam o entusiasmo dos planejadores de missão. A condutividade térmica extremamente baixa do regolito cria temperaturas notavelmente estáveis a 1,5 metro de profundidade em toda a região, imunes às oscilações diurnas da superfície. Áreas extensas onde o gelo permaneceria estável contra sublimação por períodos superiores a um ou dois bilhões de anos aparecem em múltiplas crateras com cerca de cem metros de diâmetro. As taxas de perda caem para aproximadamente um metro por bilhão de anos em temperaturas abaixo de cerca de 100 kelvins na superfície, ou em até 145 kelvins quando o depósito está enterrado sob vários metros de regolito. Enterrar preserva, e a Lua enterra sozinha, por impacto e por lenta agitação do solo.

Com a curva de módulos elásticos em função da fração de gelo de um lado e o mapa térmico tridimensional do outro, o passo seguinte foi converter temperatura em conteúdo de gelo, e conteúdo de gelo em velocidade sísmica. O modelo de velocidades resultante mostra comportamento fortemente não linear. Frações modestas de gelo alteram pouco a velocidade, porque o gelo ainda não conecta mecanicamente os grãos; à medida que o preenchimento se aproxima da saturação, a velocidade da onda P dispara de algumas centenas de metros por segundo para valores acima de três mil metros por segundo, e a onda S acompanha em escala menor. A razão entre as duas velocidades permanece quase constante em torno de 1,78 na maior parte do intervalo e sobe para aproximadamente 2,0 apenas nas frações mais altas.

Essa curvatura tem consequência prática direta. Ela significa que a sísmica não é um instrumento linear de medição de umidade, mas um detector particularmente sensível justamente na faixa que interessa à mineração: depósitos com poros bem preenchidos, exatamente aqueles que valeria a pena explorar. Umidade dispersa produz sinal fraco; blocos congelados produzem sinal forte. Como ferramenta de prospecção, esse viés é uma vantagem, porque o instrumento reage com intensidade máxima ao alvo economicamente relevante e ignora o ruído de fundo do regolito levemente hidratado que cobre boa parte da superfície lunar.

De posse dessa relação, a equipe construiu quatro realizações tridimensionais distintas do subsolo de Mons Mouton, todas partindo do mesmo modelo de fundo de densidade e velocidade derivado de perfis de profundidade publicados, e diferindo apenas nas hipóteses sobre como o gelo se distribui. A origem, a longevidade, o transporte e a distribuição de voláteis no regolito continuam sendo áreas ativas de pesquisa, e criar múltiplos cenários é a maneira honesta de lidar com essa incerteza. O limiar de 140 kelvins foi adotado como guia, com base em modelagem de longevidade de gelo em regolito coberto por poeira, e o limiar de 80 kelvins como condição mais restritiva.

Na primeira realização, todo o espaço poroso abaixo de 140 kelvins foi preenchido com gelo, criando crateras com pisos congelados e um conteúdo total de 6,2% de gelo em peso. Na segunda, o preenchimento total ocorreu abaixo de 80 kelvins, com gelo distribuído estocasticamente na faixa entre 80 e 140 kelvins, em quantidade inversamente proporcional à temperatura, produzindo pisos congelados cercados por halos de gelo e um total de 7,6% em peso. Na terceira, todo o espaço poroso abaixo de 140 kelvins recebeu conteúdo de gelo estocasticamente variável, também inversamente proporcional à temperatura, gerando halos centrados nas crateras e 6,0% em peso. Na quarta, o mesmo tratamento foi aplicado somente abaixo de 80 kelvins, resultando em gelo disperso restrito aos pisos de cratera e apenas 1,3% em peso.

Um detalhe merece atenção porque ancora todo o exercício na realidade observacional. Embora existam localmente regiões com alto conteúdo de gelo nos modelos, o percentual total de água em peso varia entre 2% e 8%, faixa compatível com o que a missão LCROSS mediu ao levantar a pluma de detritos da cratera Cabeus. Os cenários simulados não são otimismo especulativo: eles cabem dentro do envelope daquilo que já foi medido diretamente na Lua, o que dá às previsões sísmicas subsequentes um grau de credibilidade que modelos puramente teóricos não teriam.

A etapa final foi soltar ondas dentro desses quatro mundos digitais. As simulações usaram o Wave Propagation Program, código de diferenças finitas de segunda ordem em malha cartesiana e curvilínea, projetado para computação paralela de alto desempenho. O domínio simulado mede 8.000 por 8.000 metros e alcança 1.000 metros de profundidade, com espaçamento de grade de 4 metros. Para as menores velocidades presentes no modelo, de 116 metros por segundo, esse espaçamento fornece cerca de trinta pontos por comprimento de onda e garante campo de ondas preciso até 3,5 hertz, faixa que se sobrepõe à sensibilidade dos experimentos sísmicos da era Apollo.

A fonte é um tremor lunar de magnitude Mw igual a −4,5, representado por um tensor de momento de falha transcorrente posicionado no centro do domínio, a dez metros de profundidade, com direção de zero grau, inclinação do plano de falha de 35 graus, rake de 90 graus e frequência de canto de 5 hertz. Magnitude negativa não é erro de digitação: trata-se de um evento minúsculo, comparável a uma pancada forte, o tipo de fonte que uma missão robótica poderia gerar com um pequeno explosivo, com um impactador ou até com o próprio movimento de um veículo. A simulação corre por cinquenta segundos, ponto em que as ondas começam a interagir com as bordas da caixa e a criar reflexões espúrias, e registra instantâneos do campo de ondas a cada segundo.

Sensores virtuais de três componentes foram distribuídos a cada cem metros ao longo de uma linha oeste-leste de oito quilômetros de extensão, reproduzindo o arranjo de um levantamento sísmico terrestre convencional. A visualização adotada pelos autores é engenhosa: o campo de ondas de três componentes é traduzido em cores, com vermelho representando movimento vertical, verde o movimento radial e azul o movimento transversal, e a intensidade luminosa proporcional à amplitude. Regiões que aparecem esbranquiçadas nos gráficos indicam onde o espalhamento despolarizou as ondas, misturando as três direções de movimento até que nenhuma domine.

Os resultados são drásticos e dependem fortemente da geometria do gelo. Nos dois primeiros modelos, em que existem regiões de gelo contíguo, as velocidades sísmicas elevadas do regolito congelado produzem uma fronteira quase totalmente refletora, com coeficientes de reflexão superiores a 0,9 tanto para ondas P quanto para ondas S. Mais de 90% da amplitude sísmica é espalhada de volta a partir da região gelada, na forma de reflexão retroespalhada que aparece nos registros como uma onda propagando-se para trás e chegando de cinco a dez segundos depois na estação. Um sismômetro colocado dentro de qualquer das regiões permanentemente sombreadas de baixa temperatura se encontra, então, dentro de uma zona de sombra sísmica.

Dentro dessa zona, as amplitudes caem para algo entre 10% e 20% dos valores registrados nas áreas livres de gelo ao redor, dependendo do modelo. A energia que ainda é detectada do outro lado da região congelada não atravessou o obstáculo: ela o contornou por difração ao redor das bordas e do fundo, reaparecendo além das crateras com amplitude reduzida. Essas zonas de sombra se estendem por qualquer região com gelo difuso, e um levantamento tridimensional que mapeie a geometria das sombras funciona como diagnóstico direto da presença de gelo. É o mesmo princípio que permite a um médico enxergar um osso em uma radiografia: o que interessa não é o que chega ao detector, e sim a forma exata daquilo que foi bloqueado.

Quando um halo aleatório de gelo é introduzido no subsolo, como nos modelos dois e três, o campo de ondas se torna progressivamente mais espalhado, e as amplitudes embranquecem nos gráficos de polarização onde a energia das ondas superficiais é dispersada. Esses dois cenários produzem a maior atenuação do campo de ondas, com uma coda de energia despolarizada e zonas de sombra que cobrem 75% do perfil sísmico. Apenas a porção norte do domínio, com menos gelo, preserva uma frente de onda reconhecível. Mesmo o modelo com a menor quantidade de gelo distribuído ainda exibe zona de sombra atenuada na região oeste, velocidades elevadas na área gelada e espalhamento mínimo a leste, onde as anomalias de gelo praticamente inexistem.

A sensibilidade do método também aparece nas medidas mais convencionais de velocidade. Em todos os casos simulados, qualquer onda que trafegue pelas regiões impregnadas de gelo apresenta velocidade superior a mil metros por segundo, valor de duas a três vezes maior que o dos modelos de regolito circundante. Quanto maior a abundância de gelo, mais rápida a onda. Medir a velocidade do regolito entre duas estações posicionadas dentro de uma região permanentemente sombreada, e comparar com o resultado obtido fora dela, indica se há gelo presente; a magnitude da diferença indica quanto. Nos cenários mais ricos, a velocidade sobe para a faixa de 3.500 a 4.000 metros por segundo, contra cerca de mil metros por segundo no cenário mais pobre. Regiões sem gelo apresentam variações pequenas de velocidade entre estações vizinhas.

O conjunto desses achados compõe algo que se aproxima de uma gramática de leitura sísmica do gelo lunar, com regras distinguíveis umas das outras. Detecção de uma zona de sombra associada a uma região sombreada, sem perda de polarização da fonte, aponta para gelo fortemente ligado, com transições nítidas entre regiões congeladas e secas — um bloco compacto com fronteiras bem definidas. Despolarização do campo de ondas somada à perda de amplitude na estação interna à região sombreada aponta para gelo amplamente distribuído e difuso, contrastando com geometria mais esparsa no regolito circundante. Cada assinatura corresponde a um tipo diferente de depósito, e a diferença entre os dois tipos é exatamente a diferença entre um alvo de mineração e uma curiosidade científica.

A proposta operacional derivada do trabalho é direta o suficiente para caber em uma missão robótica de porte modesto. Um conjunto de sismômetros posicionados dentro e fora de uma região permanentemente sombreada, acoplado a uma fonte sísmica controlada disparada dentro ou fora dessa mesma região, permitiria distinguir entre os vários cenários de gelo no regolito e caracterizar a distribuição do depósito. Fontes passivas próximas, como impactos naturais e tremores lunares rasos, forneceriam comparação adicional dos efeitos de propagação. O que a rede mede, no fim, é a velocidade sísmica do regolito dentro e fora da sombra, e velocidades elevadas indicam a presença de gelo com sensibilidade à sua abundância.

Os autores reconhecem que outros mecanismos podem elevar velocidades sísmicas, como lâminas de fundido de impacto congeladas ou antigos derrames de lava enterrados. A discriminação vem da correlação espacial: se a anomalia de velocidade coincide com a área onde o modelo térmico prevê temperaturas baixas o bastante para preservar gelo, a coincidência entre geofísica e termodinâmica constitui evidência forte de recurso enterrado. Efeitos de espalhamento são mais circunstanciais, porque heterogeneidades não relacionadas a gelo dentro do regolito e do megaregolito também produzem espalhamento; se o espalhamento estiver localizado sobre uma região sombreada específica, o padrão pode ser usado para isolar a presença de gelo.

A honestidade metodológica do trabalho aparece nas limitações que ele próprio enumera. Os modelos de meio efetivo assumem grau de homogeneidade microestrutural que pode simplificar demais a complexidade real do regolito lunar, onde anisotropia e conectividade heterogênea entre fases provavelmente desempenham papéis importantes. As simulações omitem topografia realista e espalhamento em múltiplas escalas, fatores conhecidos por influenciar a propagação de ondas na Lua. Os experimentos de laboratório não reproduzem o vácuo extremo, a ciclagem térmica e o intemperismo espacial que atuam sobre o material lunar real ao longo de bilhões de anos.

Essas ausências, contudo, não invalidam as conclusões centrais, e os autores explicam por quê. Qualquer hipótese de heterogeneidade aleatória adicionada ao modelo produziria uma coda sísmica que seria, ela mesma, modulada pela heterogeneidade imposta pelo gelo. Os efeitos primários — zonas de sombra e aumento de velocidade — permaneceriam em dados simulados que incorporassem espalhamento realista. As medições da era Apollo já haviam demonstrado que o campo de ondas lunar é altamente espalhado por heterogeneidade induzida por impactos, ao mesmo tempo em que os tempos de chegada permanecem, em geral, preservados. É sobre os tempos de chegada e sobre a geometria das sombras que a proposta se apoia.

A dimensão histórica dessa pesquisa remonta a julho de 1969. As missões Apollo instalaram na Lua tanto experimentos sísmicos passivos quanto experimentos de fonte controlada, incluindo o experimento sísmico ativo da Apollo 14 e o perfilamento sísmico lunar da Apollo 17. Os registros revelaram um mundo que reverbera: um único impacto produzia oscilações que duravam dezenas de minutos, comportamento sem equivalente terrestre, atribuído à combinação de baixa atenuação com espalhamento intenso no regolito fraturado. Aquela rede de sismômetros operou até 1977 e definiu quase tudo o que sabemos sobre a estrutura interna da Lua. Meio século depois, o mesmo tipo de instrumento reaparece com outro objetivo: não mais entender o interior profundo, e sim prospectar os primeiros centenas de metros em busca de água.

O calendário de missões torna essa transição iminente. O sismômetro embarcado na missão indiana Chandrayaan-3 já registrou sinais no sítio de pouso lunar, demonstrando que instrumentos compactos e de baixo custo funcionam na superfície. A Estação Sísmica do Lado Oculto levará monitoramento de longa duração à face que nunca vemos. A Estação de Monitoramento Ambiental Lunar, prevista para ser instalada por astronautas durante a Artemis III nas proximidades do polo sul, oferecerá a primeira oportunidade de amostrar sismicamente o subsolo polar meridional com equipamento moderno. Instrumentos inerciais a bordo de veículos de superfície também estão sendo avaliados como sensores sísmicos secundários.

O caso do VIPER ilustra a turbulência dessa nova fase. Projetado para percorrer a região de Mons Mouton com uma broca de um metro chamada TRIDENT, o veículo foi cancelado em julho de 2024 por questões orçamentárias, com o hardware praticamente montado, e ressuscitado em setembro de 2025, quando a NASA contratou a Blue Origin para entregá-lo à superfície lunar com o módulo Blue Moon MK1, lançado por um foguete New Glenn, em missão prevista para o final de 2027. Enquanto isso, o experimento PRIME-1, que levava uma versão do mesmo sistema de perfuração a bordo do módulo Athena, chegou a Mons Mouton em março de 2025 e tombou de lado dentro de uma cratera, impedindo que a broca penetrasse o solo.

Esse acidente resume o argumento econômico a favor da sísmica. Perfurar é caro, lento, arriscado e pontual: cada furo amostra alguns centímetros quadrados de terreno e depende de o veículo estar exatamente no lugar certo, na orientação certa, com energia suficiente. Um levantamento sísmico cobre quilômetros quadrados com um punhado de sensores e uma fonte, e produz uma imagem tridimensional do volume inteiro. Na Terra, nenhuma empresa perfura antes de fazer sísmica, porque a sísmica é a etapa que transforma uma área de interesse em um alvo específico. Transportar essa lógica para a Lua reduz o custo esperado de cada quilo de água extraída, e o custo por quilo é a variável que decide se uma base lunar permanente é viável ou não.

As implicações se estendem além da prospecção de recursos. As propriedades sísmicas alteradas pela presença de gelo influenciam o comportamento mecânico geral do regolito, informação necessária para projetar fundações, plataformas de pouso, taludes escavados e qualquer estrutura que precise se apoiar naquele solo. Um piso de cratera com poros congelados suporta carga de maneira completamente diferente de um piso de regolito solto, e engenheiros que ignorarem essa diferença projetarão errado. O mesmo arcabouço serve ainda para estudos comparativos de planetologia, aplicando-se à estabilidade de gelo em outros corpos sem atmosfera, como as regiões polares permanentemente sombreadas de Mercúrio, já mapeadas com técnicas de modelagem térmica aparentadas.

Há uma qualidade rara nesse trabalho que merece registro. Ele não descobriu gelo, não observou a Lua, não trouxe amostras. O que fez foi construir uma ponte quantitativa entre quatro disciplinas — física de rochas experimental, imageamento por síncrotron, modelagem térmica planetária e sismologia computacional — para responder a uma pergunta que só faz sentido no contexto de uma economia lunar que ainda não existe. É engenharia de conhecimento antecipada, do tipo que precisa estar pronta antes que a missão decole, porque instrumentos são projetados anos antes do lançamento e a escolha de frequências, espaçamentos e geometrias de arranjo depende exatamente de saber o que se espera medir.

O artigo, submetido em junho de 2025 e aceito em junho de 2026 após revisão longa, encerra com uma constatação sóbria: esforços experimentais e de modelagem contínuos serão necessários para traduzir plenamente esses achados em estratégias confiáveis de exploração. Não há promessa de solução completa, apenas a demonstração de que o caminho é percorrível. Os dados brutos das varreduras tomográficas ficaram disponíveis publicamente em repositório aberto, e o trabalho inteiro foi publicado sob licença Creative Commons, permitindo que qualquer grupo verifique, reproduza e estenda os resultados.

Existe algo apropriado na ideia de que a água lunar será encontrada pelo silêncio que ela produz. Durante bilhões de anos, esses depósitos ficaram guardados nos únicos lugares do sistema solar interno onde o Sol nunca chega, protegidos por uma geometria orbital que mantém certos fundos de cratera abaixo de cem kelvins desde antes que a vida existisse na Terra. Encontrá-los não exigirá luz — nada ali reflete luz. Exigirá bater no chão e escutar com atenção suficiente para perceber onde o eco não volta. A próxima geração de exploradores lunares carregará brocas, escavadeiras e eletrolisadores, mas o instrumento que dirá a eles onde parar o veículo talvez seja o mais antigo e mais simples de todos: um sensor de vibração, apoiado na poeira, esperando um tremor.