Retrieving the Seismic Moment Release Rate of Small Seismic Events Using the Adjoint Method - Theory and Numerical Validation

Quando a Terra treme, até mesmo um tremor minúsculo carrega uma história oculta sobre como o solo se rompeu. Os sismólogos há muito conseguem determinar onde um terremoto começou e quanta energia ele liberou, mas compreender exatamente como essa energia foi liberada ao longo do tempo permaneceu um quebra-cabeça difícil. Essa história é contada por algo chamado função do tempo da fonte, um registro do batimento cardíaco da ruptura, desde sua primeira fissura até sua parada final. Conhecer esse ritmo é crucial porque revela a velocidade da ruptura, a direção em que ela viajou e as mudanças de tensão que ocorreram nas profundezas do subsolo. No entanto, extrair esse ritmo das ondas desordenadas registradas na superfície é notoriamente difícil. Os métodos tradicionais ou dependem de encontrar um terremoto menor e semelhante para usá-lo como referência — um processo que pode ser subjetivo e impreciso — ou tentam executar simulações massivas, pesadas para o computador, que são lentas demais para uso rotineiro.

Em um novo estudo, os pesquisadores Wojciech Dębski e Kamil Waśkiewicz, do Instituto de Geofísica da Academia Polonesa de Ciências, propuseram uma maneira diferente de ouvir esses ritmos ocultos. Eles desenvolveram uma técnica baseada na ideia física simples de que o som e as ondas sísmicas podem viajar para trás no tempo da mesma forma que viajam para frente. Ao pegar os registros da superfície, invertê-los no tempo e enviá-los de volta para o interior da terra usando um modelo computacional, as ondas naturalmente se refocam no ponto onde o terremoto começou. À medida que essas ondas revertidas convergem, elas reconstroem o padrão original de liberação de energia. A equipe testou essa ideia usando simulações computacionais de pequenos terremotos causados pelo homem, como aqueles desencadeados pela mineração, utilizando dados sintéticos livres de ruído. Seus resultados mostram que este método pode recuperar com sucesso o tempo detalhado da ruptura sem a necessidade de inversões matemáticas complexas ou de encontrar um terremoto de referência. Ele oferece uma maneira mais rápida e robusta de compreender a mecânica de pequenos eventos sísmicos, transformando um problema matemático difícil em um processo físico direto.

O cerne desta nova abordagem reside em uma propriedade fundamental da física de ondas: a reversibilidade do tempo. Em um meio perfeitamente elástico, onde a energia não é perdida para o calor ou fricção, as equações que descrevem como as ondas se movem funcionam da mesma forma, quer o tempo esteja movendo-se para frente ou para trás. Imagine um vídeo de uma pedra caindo em um lago; se você reproduzir o vídeo em reverso, as ondulações viajam de volta para dentro e a pedra salta para fora da água. Embora você não possa realmente reverter o tempo no mundo real, a matemática permite que os cientistas simulem esse processo. Os pesquisadores usaram esse conceito para criar um algoritmo que chamam de TRSTF. O processo começa pegando as ondas sísmicas registradas em várias estações na superfície. Essas ondas são então revertidas no tempo e alimentadas de volta em um modelo computacional da crosta terrestre. À medida que as ondas revertidas viajam através do modelo, elas ricocheteiam em limites e interagem com a geologia complexa, eventualmente convergindo de volta ao local original da fonte.

Quando essas ondas revertidas se encontam na fonte, elas interferem umas com as outras. No local exato do terremoto, as ondas se alinham perfeitamente, criando um sinal focado e forte. Em qualquer outro local, as ondas chegam fora de sincronia e se cancelam, deixando apenas um ruído fraco. Esse efeito de foco é a chave do método. O sinal que emerge no ponto focal é uma reconstrução direta da função do tempo da fonte, o próprio ritmo da liberação de energia do terremoto. Os pesquisadores descobriram que esta técnica funciona de forma notável em seus experimentos sintéticos livres de ruído, mesmo quando o modelo do interior da Terra não é perfeitamente conhecido. Em seus experimentos computacionais, eles simularam três cenários diferentes para mimetizar configurações de monitoramento do mundo real: um onde os sensores foram colocados próximos ao nível da atividade de mineração, um onde os sensores estavam na superfície acima do evento, e um terceiro com uma configuração de superfície diferente. Em todos os casos, o método recuperou com sucesso a forma do padrão de ruptura original.

O estudo também explorou como erros no modelo da velocidade do som da Terra afetam os resultados. No mundo real, os geólogos não conhecem a velocidade exata das ondas sísmicas em cada ponto do subsolo; eles dependem de estimativas. Os pesquisadores testaram seu método introduzindo erros aleatórios e variações em blocos em seus modelos computacionais, simulando uma situação onde o modelo de velocidade era apenas 90 por cento preciso. Os resultados indicaram que o uso de modelos de velocidade imprecisos introduz erros sistemáticos na função do tempo da fonte reconstruída. O método provou ser mais sensível a erros organizados, em blocos, do que a erros aleatórios e dispersos. No entanto, apesar desses erros, o sinal recuperado permaneceu reconhecível e forneceu informações úteis, demonstrando que o método mantém a utilidade mesmo com modelos imperfeitos, embora a precisão não seja imune a imprecisões do modelo.

Uma das descobertas mais marcantes foi como o arranjo dos sensores influenciou a qualidade do resultado. Os pesquisadores descobriram que ter sensores mais próximos dos limites da área sendo estudada realmente ajudava a reconstrução. Isso ocorre porque as reflexões desses limites atuam como sensores virtuais adicionais, fornecendo mais informações ao sistema. Em um de seus casos de teste, onde os sensores foram colocados perto da borda superior do modelo, o sinal reconstruído foi mais nítido e preciso do que nas outras configurações, mesmo que os dados brutos parecessem mais distorcidos devido à interferência. Este resultado contraintuitivo sugere que o método se beneficia da complexidade do campo de ondas, usando as reflexões a seu favor, em vez de tratá-las como ruído.

Os pesquisadores também mediram características específicas dos sinais reconstruídos, como a rapidez com que a ruptura começou e quanto tempo durou. Eles descobriram que, embora os sinais reconstruídos fossem ligeiramente mais largos que os originais, provavelmente devido ao número limitado de sensores usados na simulação, a forma geral e o tempo foram preservados. A capacidade de estimar a duração da ruptura é vital porque ajuda os cientistas a calcular a queda de tensão dinâmica, uma medida de quanta tensão foi liberada durante o evento. Essa informação é essencial para compreender a mecânica da sismicidade induzida, como terremotos causados por mineração ou injeção de fluidos. O estudo mostrou que o novo método pode estimar esses parâmetros com uma precisão comparável ou melhor que as técnicas existentes, mas com muito menos esforço computacional.

Ao contrário da inversão de forma de onda completa, que requer a execução de milhares de simulações e loops de otimização complexos para encontrar o melhor ajuste, este novo método requer apenas uma única passagem das ondas revertidas pelo modelo. Isso o torna computacionalmente eficiente, capaz de rodar em computadores desktop padrão para modelos pequenos ou supercomputadores para escalas regionais. Ele também evita a necessidade de selecionar um terremoto de referência, um passo que pode introduzir viés humano e incerteza nos métodos tradicionais. Os pesquisadores enfatizam que sua abordagem é particularmente bem adequada para pequenos eventos sísmicos, onde a fonte é pequena o suficiente para ser tratada como um único ponto. Eles reconhecem que o método assume que o solo se comporta como um material elástico perfeito, o que significa que não absorve energia, e que a fonte é pequena em comparação com o comprimento de onda das ondas. Estas são suposições razoáveis para muitos eventos relacionados à mineração, mas os pesquisadores observam que mais estudos são necessários para entender como o método se comporta em ambientes do mundo real mais complexos, com perda significativa de energia.

O artigo conclui destacando o potencial desta técnica de se tornar uma ferramenta padrão na sismologia. Ao fornecer uma maneira rápida, confiável e fisicamente intuitiva de recuperar a função do tempo da fonte, ela abre as portas para a análise rotineira de pequenos eventos sísmicos que anteriormente eram difíceis de caracterizar em detalhes. Os pesquisadores sugerem que este método pode ser integrado aos sistemas de monitoramento existentes para fornecer insights imediatos sobre a natureza da sismicidade induzida. Embora o estudo atual tenha sido limitado a simulações computacionais, a base teórica é sólida e os resultados são promissores. O próximo passo, segundo os autores, será testar o método em dados sísmicos reais e refinar a compreensão de suas limitações e incertezas. Por enquanto, o estudo serve como uma prova de conceito de que a reversão do tempo pode ser uma lente poderosa para olhar para o passado da Terra, revelando a história oculta de como o solo se rompe.