MODELO DE VELOCIDADES PARA LA REGIÓN DE AÑELO, A PARTIR DE FUNCIONES RECEPTORAS, CURVAS DE DISPERSIÓN Y LOCALIZACIÓN DE SISMOS LOCALES

CORREA OTTO SEBASTIAN; SPAGNOTTO, SILVANA; NACIF, SILVINA

Congreso; XXI Congreso Geológico Argentino; 2022

Se obtuvieron modelos de velocidad unidimensionales en la región de Añelo, ubicada dentro de la cuenca extensional Neuquina. Para ello se trabajó con datos de una red temporal de estaciones de banda ancha, desplegada por la Universidad Nacional de San Juan – Argentina tomados desde nes del mes de octubre de 2014 hasta junio de 2016. El equipamiento de las estaciones constó de Sismómetros Trillium 120P/PA de 3-componentes, very broadband y registradores Taurus. El área de estudio ocupa un sector sur de la cuenca, el Bloque Loma La Lata/Sierra Barrosa y aoramientos relacionados a la dorsal Neuquina (Silvestroy Zubiri 2008). Estructuras neotectónicas han sido documentadas en la región, indicando contracción Andina E-O principalmente (Folguera y Ramos 2011, Messager et al. 2014, Sagripanti et al. 2015, Branellec et al.2016), sin embargo, son escasos los experimentos sismológicos locales que reporten modelos de velocidad en la zona. El trabajo de Bohm et al 2002, reporta modelos localizados en Cordillera Principal y en el ante-arco chileno, y presenta escasas capas superciales, además no registra actividad sísmica en el área de estudio de este trabajo.Para la obtención de los modelos, se utilizaron dos técnicas diferentes que poseen resolución a escalas y profundidades distintas. Por un lado, se obtuvieron niveles nos superiores usando localización de sismoslocales y tiempos de arribo de las ondas P y S, y por otro, capas más gruesas y profundas, que permitieron determinar las discontinuidades de Mohorovičić y nivel de despegue utilizando inversión conjunta de fun-ciones receptoras y curvas de dispersión obtenidas con ruido. Los primeros kilómetros de profundidad del modelo se obtuvieron a partir de 115 eventos locales con un promedio de error menor a 5 km en parámetros hipocentrales. Para esto, se usó el software Velest (Kissling et al. 1995) que presenta una solución iterativa del problema acoplado de hipocentros - modelo de velocidad - correcciones de estación. Este resuelve conjuntamente la determinación de hipocentros, el cálculo del modelo de velocidad y las correcciones de estación que hagan minimizen el RMS (error medio cuadrático). Las funciones receptoras se obtuvieron utilizando la metodología de deconvolución iterativa de Ligorría y Ammon (1999) que utiliza una serie de pulsos Gaussianos en la deconvolución de las señales de entrada. Se utilizaron anchos de pulso Gaussiano de ancho 1.0 (f < 0.5 Hz) y 2.5 (f < 1.2 Hz) y se analizaron todos los telesismos de magnitudes mayores a 5.8 ocurridos en el periodo que va de Octubre de2014 a Enero de 2016. Con las funciones receptoras se hizo una estimación inicial en profundidad de las principales discontinuidades. La función del receptor no es sensible a la velocidad P y el espesor de lacorteza se estima a partir del tiempo de retardo de la fase convertida en la Moho de P a S y la relación Vp/Vs de la corteza. Para la inversión conjunta de las funciones receptoras con las curvas de dispersión y su posterior obtención de un modelo de velocidades se usó la metodología de Julia et al. (2000, 2003). Dicha metodología combina resultados independientes en una inversión linealizada que minimiza las normas LSQ ponderadas de los vectores residuales entre datos observados (función de receptor o velo-cidad de dispersión) y previsiones del modelo, con una norma de suavidad (regularización) de derivada segunda (del modelo de velocidad) nula. Las curvas de dispersión usadas fueron las curvas de dispersión obtenidas por Spagnotto et al. 2016 e integradas con curvas de dispersión regionales de largo periodo (<100 s) publicadas por Feng et al. (2004). A partir de los resultados con RF y dispersión se hizo una estimación en profundidad de las principales discontinuidades mostrada en la Fig. 1. Los modelos obtenidos arrojan claramente dos discontinuidades, unade 8 a 10 km relacionada a un posible nivel de despegue de las estructuras principales en la zona y otra de 30 a 35 km relacionada a la Moho. Los valores reportados en la zona incluyen las funciones receptoras de Yuan et al. 2006 y otros trabajos que utilizan métodos potenciales como los de Tassara y Echaurren (2012), Assumpsao et al. (2013), y Van der Meijde et al. (2013) y están en acuerdo con el valor encontrado. Por otro lado, los valores de velocidad de las ondas por debajo de la Moho son un tanto menores a los valores promedio y podrían indicar la presencia de una pluma del manto (Gianni et al. 2017) o ser consecuencia del uso de las curvas de dispersión regionales (Feng et al. 2004) que presentan estos valores anómalos en toda la región.