Como entender e modelar estruturas internas da terra na prática
A maioria dos manuais apresenta a Terra como uma sequência ordenada de camadas com limites bem definidos. Na prática, esses limites são zonas de transição com espessura de dezenas a centenas de quilômetros, e a resolução dos modelos varia drasticamente dependendo da distribuição dos sismógrafos e do tipo de dado usado.
O que compõe as estruturas internas da terra
A crosta varia entre 5 e 70 km de espessura, sendo mais fina no oceânico e mais grossa nos continentes. O manto se estende até cerca de 2900 km de profundidade e abriga descontinuidades importantes: a de 410 km, associada à transformação de olivina em wadsleyíta, e a de 660 km, onde os minerais se separam em perovskita de silicato de magnésio e Bridgmanita. A descontinuidade deovich marca a base da crosta. O núcleo externo é líquido, composto principalmente de ferro e níquel, e é responsável pelo campo magnético terrestre por meio do efeito dínamo. O núcleo interno é sólido, com raio de aproximadamente 1220 km, e exibe anisotropia sísmica que varia conforme a direção de propagação das ondas.
Uma coisa que poucos mencionam: a velocidade das ondas P e S aumenta com a pressão, mas diminui com a temperatura. Isso significa que uma região de alta velocidade pode ser tanto fria quanto quimicamente distinta. Separar esses dois efeitos sem dados complementares de física mineral é praticamente impossível.
Método de imageamento sísmico
O imageamento por tomografia sísmica funciona de forma semelhante a uma TC médica, mas com uma diferença crítica: os raios sísmicos vêm de fontes naturais distribuídas de maneira irregular ao redor do planeta. Isso gera zonas de sombra enormes, especialmente no manto inferior abaixo de regiões sem estações sísmicas densas. O processo começa com a coleta de tempos de chegada de fases sísmicas, como P, S, PKP e SKS. Cada fase carrega informação diferente sobre o trajeto que percorreu. Fases que atravessam o núcleo, como PKP, sensibilizam o núcleo externo, enquanto fases que refletem na interface núcleo-manto, como SKS, não penetram no núcleo líquido.
A inversão propriamente dita minimiza a diferença entre os tempos observados e os tempos calculados por um modelo inicial. O modelo é atualizado iterativamente usando a matriz jacobiana, que relaciona variações de velocidade com variações nos tempos de chegada. A regularização é essencial aqui — sem ela, o modelo tende a produzir anomalias espúrias de pequena escala que não têm correspondência física. Uma limitação importante: a resolução típica em modelos globais como o SEED ou o LLNL-2D é da ordem de 100 a 200 km no manto superior e piora para 300 a 500 km no manto inferior. Se você precisa resolver estruturas mais finas, como placas subduzidas com menos de 50 km de espessura, modelos globais não são suficientes.
👉 Clique no botão abaixo para saber mais sobre o assunto!
Um problema real que encontrei
Trabalhei com dados de estações broadband da rede brasileira dos anos 1990 para mapear a litosfera sul-americana. Os tempos de chegada de P pareciam normais à primeira vista, mas as anomalias de velocidade no manto superior apresentavam padrões que não faziam sentido geológico. A causa era um erro sistemático de temporização de aproximadamente 300 ms em algumas estações, comum em sistemas de aquisição analógica da época. A solução foi usar estações vizinhas que observavam os mesmos eventos teleseísmicos. Ao calcular os tempos relativos entre pares de estações e comparar com as previsões do modelo IASP91, identifiquei quais estações tinham desvio consistente. Apliquei correções de atraso individual e o modelo resultante mostrou estruturas compatíveis com a evolução tectônica conhecida da plataforma sul-americana.
Duas coisas que aprendi na prática
A primeira é sobre a descontinuidade de 660 km. Ela não é uma superfície plana — sua profundidade varia conforme o gradiente térmico regional. Em zonas de subducção, onde o material é frio, a descontinuidade afunda até 700 km ou mais. Em pontos quentes, como no centro de bacias oceânicas jovens, ela sobe para cerca de 600 km. Tratar essa descontinuidade como referência fixa em estudos regionais gera distorções significativas. A segunda é sobre a distinção entre perturbações térmicas e composicionais. Uma anomalia de alta velocidade no manto pode ser interpretada como uma placa fria subduzida ou como um domo mantélico composiçãolemente distinto. Sem dados de ondas de corpo combinados com modos normais e informações geoquímicas, fica difícil decidir. Modelos que reportam apenas imagens de velocidade sem discutir a ambiguidade térmica-composicional estão incompletos.
A fronteira entre núcleo e manto, conhecida como descontinuidade de Gutenberg, também não é lisa. A zona D'' na base do manto, com espessura variável entre 200 e 300 km, contém estruturas complexas como prováveis vestígios de placas subduzidas acumuladas ou anomalias térmicas associadas a grandes províncias de baixo cisalhamento (LLSVPs). Esses corpos têm implicações diretas na convecção do manto e no vulcanismo de pontos quentes, mas sua natureza exata permanece debatida.
Onde encontrar dados
Modelos sísmicos globais prontos para uso estão disponíveis no IRIS DMC, que hospeda versões do ak135, PREM e modelos tomográficos mais recentes. Para dados brutos de sismogramas, o serviço ADS (Array Data Search) do IRIS permite baixar formatos SAC, SEG-Y e miniSEED de milhares de estações. O banco de dados ISC fornece tempos de chegada publicados para eventos globais, útil para cross-checking. Software como Spececho, ObsPy e SeisFlow facilita o processamento e a análise. O ObsPy é gratuito e em código aberto, com boa documentação para leitura, filtragem eicks de sismogramas. O SeisFlow é mais voltado para imageamento de estrutura crostica usando métodos de múltipla correlação.
Erros comuns ao trabalhar com modelos de camadas
Interpolar modelos globais em escalas regionais sem considerar a diferença de resolução é um erro frequente. Um modelo com resolução de 2 graus pode parecer detalhado, mas em escala local perde completamente a capacidade de resolver falhas e descontinuidades menores que 200 km. Sempre verifique a resolução espacial do modelo antes de usá-lo para interpretar estruturas crosticas. Outro erro é confundir velocidades sísmicas com composição rochosa diretamente. A relação entre velocidade e composição depende de pressão, temperatura e estado de fusão. A mesma rocha pode produzir velocidades diferentes em profundidades distintas. Usar curvas de velocidade de laboratório sem corrigir para condições in-situ leva a interpretações erradas de litologia.
O manto terrestre não é um meio isotrópico simples. A anisotropia sísmica, causada pelo alinhamento preferencial de cristais de olivina durante o fluxo do manto, afeta tanto ondas S quanto P. Ignorar esse efeito em estudos de estrutura profunda pode introduzir vieses sistemáticos na estimação de velocidades.
