Modelamiento termo-mecánico en la zona de subducción del Caribe colombiano -área Sinú San Jacinto- a partir de la interpretación de anomalías de campos potenciales

Abstract

Se propone una estructura térmica para la zona de subducción del Caribe colombiano, área Sinú-San Jacinto, en la esquina Noroccidental de Suramérica, bajo las hipótesis de subducción plana para la placa Caribe bajo la placa Suramérica y un origen Pacífico Oriental para la placa Caribe. La estructura térmica de la zona se obtiene al aplicar los códigos Temspol (cinemático) e i3Elvis (dinámico) en Matlab, que resuelven las ecuaciones de energía, masa, momento y Navier-Stoke, bajo un esquema numérico de diferencias finitas. En el modelamiento cinemático se resuelve solamente la ecuación de energía con una geometría fija y predefinida. En el modelamiento termo-cinemático-reológico, se incorporan parámetros reológicos a los nodos del modelo termo-cinemático. En el modelamiento termo-mecánico, se resuelven en forma simultánea las ecuaciones de energía, masa, momento y Navier-Stoke; aquí, la geometría no es fija, es producto de: (1) la libertad del movimiento de la placa en subducción; (2) los parámetros físicos de la placa; (3) las condiciones de frontera del modelo; y (4) las propiedades reológicas de la litología de las capas que componen el modelo. La estructura térmica propuesta, se valida mediante la comparación del flujo de calor estimado por los modelos termo-mecánico y termo-cinemático, con el flujo de calor estimado por el método de la profundidad del punto de Curie (CPD), calculado mediante la técnica espectral, en su modalidad fractal, para anomalías magnéticas. La interpretación de la estructura térmica está sujeta a restricciones impuestas por: (1) el modelamiento directo bidimensional del subsuelo a partir de señales de métodos potenciales; (2) las estimaciones de la frontera cortical frágil-dúctil, desde los modelos térmicos y los análisis isostáticos, en sus mecanismos de compensación tipo Airy y regional, en sus formas unidimensional y bidimensional, bajo las técnicas espacial y espectral; y (3) la información geológica y geofísica, principalmente de los modelos evolutivos paleo-tectónicos y palinpásticos. Se encontró que: (1) el flujo de calor estimado a partir de los modelos térmicos es coherente con el flujo de calor estimado por el método del CPD; (2) la distribución espacial del flujo de calor superficial es válida para el Caribe colombiano; (3) el flujo de calor estimado conserva la misma forma que el flujo de calor en una zona de subducción en el noroeste de Japón; y (4) la interpretación de la estructura térmica propuesta es coherente con las reconstrucciones paleo-tectónicas, donde se plantea que la evolución de la zona de convergencia entre las placas Caribe y Suramérica probablemente ha sufrido dos eventos subductivos, uno en el Cretácico-Paleógeno y otro en el Mioceno Medio. Se concluye que el régimen térmico del Caribe colombiano está caracterizado por: la subducción plana de la placa Caribe bajo la placa Suramérica, con un ángulo de subducción menor a 30° y una velocidad de convergencia promedio de 2 cm/año; una litosfera oceánica que contiene una corteza oceánica caribeña de edad Cretácica, cuya segunda fase de subducción inició hace 45 Ma, aproximadamente; y que el régimen térmico en la zona de estudio es un 20% menor que el régimen térmico calculado para una zona de subducción normal.

The objective of this work is to propose a thermal structure for the Colombian Caribbean subduction zone, Sinú-San Jacinto área, in the Northwestern corner of South America, under the hypothesis that the Caribbean plate presents flat subduction under the South American plate, with Eastern Pacific origin; associated to the Caribbean plate. The thermal structure of the study area was obtained by applying the Temspol (kinematic) and i3Elvis (dynamic) Matlab codes. Such codes solve the energy, mass, momentum, and Navier-Stokes equations, under a finite difference numerical scheme. Concerning to kinematic modeling, the geometry is fixed and predefined, involving only the energy equation; which is solved. In kinematic-rheological modeling, rheological parameters are incorporated in thermo-kinematic model’s nodes. In thermo-mechanical modeling, energy, mass, momentum, and Navier-Stokes equations are solved simultaneously; geometry is not fixed, as in kinematic modeling, it is the product of: (1) free movement of the plate in subduction; (2) physical plate parameters; (3) model boundary conditions; and (4) rheological properties based on a stratified media model. The proposed thermal structure in this work is well founded by comparing the estimated heat flux from the thermo-mechanical and thermo-kinematic models, with the heat flux estimated by the Curie point depth methodology (CPD), through the spectral technique in its fractal modality, applied to magnetic anomalies. The thermal structure interpretation is subject to restrictions such as: (1) two-dimensional direct modelling based on potential fields data (2) cortical boundary fragile-ductile based on estimations of thermal models and the isostatic analysis based on Airy and regional compensation mechanisms, in one and two-dimensional forms, using spatial and spectral techniques; and (3) geological and geophysical information, mainly from paleo-tectonic and palinspastic evolutionary models. Relevant remarks are: heat flux estimated from the thermal models is consistent with the heat flux estimated by the CPD method; spatial distribution of surface heat flux is valid for the Colombian Caribbean sea; heat flux estimated in this work retains the same shape as the heat flux in a similar subduction zone in northwestern Japan; the interpretation of the proposed thermal structure is consistent with the paleo-tectonic reconstructions as well, where the convergence zone evolution between the Caribbean and South American plates has probably undergone two subductive events, first in the Cretaceous-Paleogene, and another later in the Middle Miocene. Such evidences lead to the conclusion that the thermal regime of the Colombian Caribbean sea is characterized by the flat subduction of the Caribbean plate under the South American plate, with a subduction angle lower than 30° and an average convergence speed of 2 cm/year; the existence of an oceanic lithosphere containing a Caribbean oceanic crust of Cretaceous age, whose second subduction phase began approximately 45 Ma ago; and thermal regime in the study zone is 20% lower than the calculated thermal regime for a normal subduction zone.