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

Cuando la Tierra se sacude, incluso un pequeño temblor transporta una historia oculta sobre cómo se rompió el suelo. Los sismólogos han sido capaces durante mucho tiempo de localizar dónde comenzó un terremoto y cuánta energía liberó, pero comprender exactamente cómo se liberó esa energía a lo largo del tiempo ha seguido siendo un rompecabezas difícil. Esta historia es contada por algo llamado la función de tiempo de la fuente, un registro del latido de la ruptura desde su primera grieta hasta su parada final. Conocer este ritmo es crucial porque revela la velocidad de la ruptura, la dirección en la que viajó y los cambios de estrés que ocurrieron en las profundidades del subsuelo. Sin embargo, extraer este ritmo de las ondas desordenadas registradas en la superficie es notoriamente difícil. Los métodos tradicionales dependen de encontrar un terremoto más pequeño y similar para usarlo como referencia —un proceso que puede ser subjetivo e impreciso— o intentan ejecutar simulaciones masivas, pesadas para la computadora, que son demasiado lentas para un uso rutinario.

En un nuevo estudio, los investigadores Wojciech Dębski y Kamil Waśkiewicz, del Instituto de Geofísica de la Academia Polaca de Ciencias, han propuesto una forma diferente de escuchar estos ritmos ocultos. Desarrollaron una técnica basada en la simple idea física de que el sonido y las ondas sísmicas pueden viajar hacia atrás en el tiempo de la misma manera que viajan hacia adelante. Al tomar los registros de la superficie, invertirlos en el tiempo y enviarlos de vuelta al interior de la tierra usando un modelo computacional, las ondas se enfocan naturalmente en el lugar donde comenzó el terremoto. A medida que estas ondas invertidas convergen, reconstruyen el patrón original de liberación de energía. El equipo probó esta idea utilizando simulaciones por computadora de pequeños terremotos causados por el hombre, tales como los provocados por la minería, utilizando datos sintéticos libres de ruido. Sus resultados muestran que este método puede recuperar con éxito el cronometraje detallado de la ruptura sin necesidad de inversiones matemáticas complejas o de encontrar un terremoto de referencia. Ofrece una forma más rápida y robusta de entender la mecánica de los pequeños eventos sísmicos, convirtiendo un problema matemático difícil en un proceso físico directo.

El núcleo de este nuevo enfoque se basa en una propiedad fundamental de la física de ondas: la reversibilidad temporal. En un medio perfectamente elástico, donde la energía no se pierde en forma de calor o fricción, las ecuaciones que describen cómo se mueven las ondas funcionan de la misma manera ya sea que el tiempo avance hacia adelante o hacia atrás. Imagine un video de una piedra cayendo en un estanque; si reproduce el video en reversa, las ondas viajan de regreso hacia el centro y la piedra salta fuera del agua. Aunque usted no puede revertir realmente el tiempo en el mundo real, las matemáticas permiten simular este proceso. Los investigadores utilizaron este concepto para crear un algoritmo que llaman TRSTF. El proceso comienza tomando las ondas sísmicas registradas en varias estaciones en la superficie. Estas ondas son luego invertidas en el tiempo y alimentadas de nuevo en un modelo computacional de la corteza terrestre. A medida que las ondas invertidas viajan a través del modelo, rebotan en los límites e interactúan con la compleja geología, convergiendo finalmente de regreso en la ubicación original de la fuente.

Cuando estas ondas invertidas se encuentran en la fuente, interfieren entre sí. En el lugar exacto del terremoto, las ondas se alinean perfectamente, creando una señal fuerte y enfocada. En cualquier otro lugar, las ondas llegan desfasadas y se cancelan entre sí, dejando solo un ruido tenue. Este efecto de enfoque es la clave del método. La señal que emerge en el punto focal es una reconstrucción directa de la función de tiempo de la fuente, el mismísimo ritmo de la liberación de energía del terremoto. Los investigadores encontraron que esta técnica funciona notablemente bien en sus experimentos sintéticos libres de ruido, incluso cuando el modelo del interior de la Tierra no se conoce perfectamente. En sus experimentos por computadora, simularon tres escenarios diferentes para imitar configuraciones de monitoreo del mundo real: uno donde los sensores se colocaron cerca del nivel de la actividad minera, uno donde los sensores estaban en la superficie sobre el evento, y un tercero con una disposición de superficie diferente. En cada caso, el método recuperó con éxito la forma del patrón de ruptura original.

El estudio también exploró cómo los errores en el modelo de la velocidad del sonido de la Tierra afectan los resultados. En el mundo real, los geólogos no conocen la velocidad exacta de las ondas sísmicas en cada punto del subsuelo; dependen de estimaciones. Los investigadores probaron su método introduciendo errores aleatorios y variaciones en bloque en sus modelos computacionales, simulando una situación en la que el modelo de velocidad era solo un 90 por ciento preciso. Los resultados indicaron que el uso de modelos de velocidad inexactos introduce errores sistemáticos en la función de tiempo de la fuente reconstruida. El método resultó ser más sensible a los errores organizados y en bloque que a los errores aleatorios y dispersos. Sin embargo, a pesar de estos errores, la señal recuperada seguía siendo reconocible y proporcionaba información útil, demostrando que el método mantiene su utilidad incluso con modelos imperfectos, aunque la precisión no es inmune a las inexactitudes del modelo.

Uno de los hallazgos más sorprendentes fue cómo la disposición de los sensores influyó en la calidad del resultado. Los investigadores descubrieron que tener sensores más cerca de los límites del área que se estudia realmente ayudaba a la reconstrucción. Esto se debe a que las reflexiones de estos límites actúan como sensores virtuales adicionales, proporcionando más información al sistema. En uno de sus casos de prueba, donde los sensores se colocaron cerca del borde superior del modelo, la señal reconstruida fue más nítida y precisa que en las otras configuraciones, a pesar de que los datos brutos parecían más distorsionados debido a la interferencia. Este resultado contraintuitivo sugiere que el método se beneficia de la complejidad del campo de ondas, utilizando las reflexiones a su favor en lugar de tratarlas como ruido.

Los investigadores también midieron características específicas de las señales reconstruidas, como la rapidez con la que comenzó la ruptura y cuánto duró. Encontraron que, si bien las señales reconstruidas eran ligeramente más anchas que las originales, probablemente debido al número limitado de sensores utilizados en la simulación, la forma general y el cronometraje se preservaron. La capacidad de estimar la duración de la ruptura es vital porque ayuda a calcular la caída de estrés dinámico, una medida de cuánta tensión se liberó durante el evento. Esta información es esencial para comprender la mecánica de la sismicidad inducida, como los terremotos causados por la minería o la inyección de fluidos. El estudio mostró que el nuevo método podía estimar estos parámetros con una precisión comparable, o mejor que las técnicas existentes, pero con mucho menos esfuerzo computacional.

A diferencia de la inversión de forma de onda completa, que requiere ejecutar miles de simulaciones y complejos bucles de optimización para encontrar el mejor ajuste, este nuevo método requiere solo una única pasada de las ondas invertidas a través del modelo. Esto lo hace computacionalmente eficiente, capaz de ejecutarse en computadoras de escritorio estándar para modelos pequeños o en supercomputadoras para escalas regionales. También evita la necesidad de seleccionar un terremoto de referencia, un paso que puede introducir sesgo humano e incertidumbre en los métodos tradicionales. Los investigadores enfatizan que su enfoque es particularmente adecuado para eventos sísmicos pequeños, donde la fuente es lo suficientemente pequeña como para ser tratada como un único punto. Reconocen que el método asume que el suelo se comporta como un material elástico perfecto, lo que significa que no absorbe energía, y que la fuente es pequeña en comparación con la longitud de onda de las ondas. Estas son suposiciones razonables para muchos eventos relacionados con la minería, pero los investigadores señan que se necesita más estudio para comprender cómo se desempeña el método en entornos más complejos del mundo real con una pérdida de energía significativa.

El artículo concluye resaltando el potencial de esta técnica para convertirse en una herramienta estándar en sismología. Al proporcionar una forma rápida, confiable y físicamente intuitiva de recuperar la función de tiempo de la fuente, abre la puerta al análisis rutinario de pequeños eventos sísmicos que anteriormente eran demasiado difíciles de caracterizar en detalle. Los investigadores sugieren que este método podría integrarse en los sistemas de monitoreo existentes para proporcionar información inmediata sobre la naturaleza de la sismicidad inducida. Aunque el estudio actual se limitó a simulaciones por computadora, la base teórica es sólida y los resultados son prometedores. El siguiente paso, según los autores, será probar el método con datos sísmicos reales y refinar la comprensión de sus limitaciones e incertidumbres. Por ahora, el estudio es una prueba de concepto de que la reversión del tiempo puede ser una lente poderosa para mirar hacia el pasado de la Tierra, revelando la historia oculta de cómo se rompe el suelo.