INTRODUCCIÓN
Los arcos magmáticos en el borde occidental de Colombia han sido objeto de estudio en las últimas décadas. Fueron descritos por primera vez por Aspden et al. (1987). Esos eventos han sido atribuidos a la evolución de un sistema de subducción (Barrero, 1979; McCourt et al., 1984; Aspden et al., 1987; Grosser, 1989; Boschman et al., 2014) y a sucesivas acreciones de terrenos oceánicos al continente, bajo regímenes tectónicos compresivos y transpresivos de primer orden. (Duque-Caro, 1990; Estrada, 1995; Cediel et al., 2003; Pindell et al., 2005; Lissina, 2005; Ramos, 2009; Pindell y Kennan, 2009; Villagómez et al., 2011; Montes et al., 2012; Bayona et al., 2012; Spikings et al., 2014; Whattam y Stern, 2014; Boschman et al., 2014; Cardona et al., 2018; Montes et al., 2019).
En la Cordillera Central (CC) (Figura 1a) el magmatismo Cretácico-Paleoceno tiene, principalmente, edades de cristalización entre 99 Ma y 58 Ma (Saenz et al., 1996; Restrepo et al., 2009; Leal Mejía, 2011; Bayona et al., 2012; Villagómez y Spiking, 2013; Jaramillo et al., 2017; Cardona et al., 2020) y ha sido asociado a la subducción oblicua de la placa Farallón bajo Sudamérica y a procesos pos-colisionales generados por la subducción de la Placa Caribe (Villagómez y Spikings, 2013; Jaramillo et al., 2017; Montes et al., 2019). Sobre la Cordillera Occidental (CO) (Figura 1a) este magmatismo se ha reportado con edades desde el Coniaciano - Maastrichtiano, asociándolo a procesos extensionales en el sistema de subducción (Barrero, 1979). En trabajos recientes y a partir de análisis geoquímicos y paleomagnéticos, el magmatismo Cretácico-Paleoceno (o incluso abarcando hasta el Eoceno) es considerado intraoceánico, acrecionado al continente durante el Mioceno como parte integral del Bloque Panamá-Chocó (Duque-Caro, 1990; Cediel et al., 2003; Ramos, 2009; Bayona et al., 2012; Pindell et al., 2005; Villagómez et al., 2011; Cardona et al., 2018; Montes et al., 2019) (Figura 1b). Esta hipótesis ha sido bien acogida en la literatura, aunque no se tengan estudios detallados de los límites en el borde continental, donde se propone que se acrecionaron dichos terrenos. Esos límites tectónicos han sido atribuidos a fallas transcurrentes, como fue propuesto por Duque-Caro (1990), quien sugirió la falla de Uramita como la sutura occidental, y la zona de deformación de Istmina, ubicada entre las fallas Garrapatas e Istmina, como la sutura sur (Duque-Caro, 1990 y Cediel et al., 2003). Sin embargo, hasta la fecha no se han documentado características geológicas, geométricas y cinemáticas que permitan definirlas como zonas de sutura de la acreción del Bloque Panamá-Choco. Los estudios publicados sobre el origen y evolución de los arcos magmáticos Paleoceno-Eoceno no consideran, en el borde continental colombiano, la continuidad espacial y temporal del arco magmático controlado estructuralmente que se propone en este trabajo.

Figura 1 a) Marco tectónico regional del segmento noroccidental de Suramérica (Acosta et al., 2007). b) Esquema representativo de la sucesión de terrenos acrecionados al continente colombiano. Etapas de los procesos acrecentivos: Terreno Amaime en el Cretácico Temprano, plateau oceánico Dagua-Piñón en el Cretácico Tardío y terreno Panamá-Chocó durante el Mioceno medio (Ramos, 2009).
En los últimos cinco años, la Agencia Nacional de Hidrocarburos de Colombia (ANH), patrocinó la obtención de nuevos datos geocronológicos (Tabla 1) y geoquímicos (Tablas 2 y 3) del occidente colombiano. El objetivo de este trabajo es definir el marco tectónico en el cual ocurrió el evento magmático Paleoceno-Eoceno, basados en la relación del evento magmático con el control estructural, así como en la continuidad espacial y temporal de los cuerpos ígneos. El análisis de los datos nuevos y los recopilados (Tabla 4), incluyen dataciones U/Pb en circón y análisis de elementos mayores y tierras raras (Tablas S1 y S2, en suplemento electrónico). A partir de esa información, aquí se propone un modelo alternativo de evolución tectónica del occidente colombiano durante el lapso 60 Ma a 26 Ma.
Tabla 1 Edades isotópicas U/Pb en circones de unidades litológicas de la Cordillera Occidental. Desviación media ponderada al cuadrado (MSWD).
| Muestra | Longitud | Latitud | Unidad |
|
MSWD | Tipo de Roca |
|---|---|---|---|---|---|---|
| 1246-1 | -77.17032 | 8.32407 | Batolito de Acandí | 46.27 ± 0.28 | 1.5 | Tonalitas - granodioritas |
| 1246-4 | -77.18290 | 8.31515 | Batolito de Acandí | 46.07 ± 0.49 | 3.1 | Tonalitas - granodioritas |
| 40-006-CPC | -76.34890 | 6.34199 | Batolito de Mandé | 50.5 ± 4.1 | 2.2 | Granodiorita - tonalita |
| 90-005-CPC | -76.35251 | 6.33448 | Batolito de Mandé | 52.9 ± 1.5 | 3.8 | Tonalita horbléndica |
| 40-007-CPC | -76.33392 | 6.36684 | Mandé-Santa Cecilia la Equis | 59.2 ± 1.6 | 2.0 | Tonalita porfirítica |
Tabla 2 Resultados analíticos de óxidos mayores (% Wt) del Stock de Timbiquí. Pérdida por ignición (LOI).
| Muestra | SiO2 | TiO2 | Al2O3 | (Fe2O3)t | MnO | MgO | CaO | Na2O | K2O | P2O5 | LOI | Total |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| APO-0056-LG | 55.47 | 0.6 | 17.67 | 7.81 | 0.15 | 4.26 | 8.72 | 3.04 | 0.11 | 0.09 | -- | 97.92 |
| APO-0049-LG | 60.64 | 0.47 | 17.55 | 5.95 | 0.15 | 2.93 | 5.77 | 3.52 | 0.92 | 0.15 | -- | 98.05 |
| APO-0107-D | 56.93 | 0.52 | 17.65 | 7.79 | 0.17 | 3.39 | 8.07 | 2.96 | 0.23 | 0.13 | -- | 97.84 |
| APO-0091-P | 59.25 | 0.45 | 18.14 | 5.51 | 0.16 | 1.79 | 5.33 | 3.52 | 3.31 | 0.26 | -- | 97.72 |
| CAQ-0416-P | 42.9 | 0.84 | 18.69 | 13.04 | 0.21 | 5.96 | 12.99 | 1.97 | 0.78 | 0.56 | -- | 97.94 |
| MZD-0010-L | 53.09 | 0.58 | 18.95 | 7.93 | 0.21 | 2.94 | 8.48 | 3.46 | 2.31 | 0.43 | -- | 98.38 |
| MZD-0022-R | 64.43 | 0.28 | 16.79 | 3.64 | 0.08 | 1.05 | 4.99 | 4.12 | 1.43 | 0.13 | -- | 96.94 |
| CLM-0447-L | 52.78 | 0.7 | 18.09 | 8.72 | 0.16 | 4.98 | 9.13 | 2.98 | 0.36 | 0.12 | -- | 98.02 |
| CDG-0255-P | 54.16 | 0.85 | 16.99 | 7.98 | 0.13 | 5.07 | 10.13 | 2.51 | 0.04 | 0.12 | -- | 97.98 |
| CDG-0256-P | 52.19 | 0.63 | 16.95 | 7.8 | 0.15 | 6.87 | 12.44 | 1.57 | 0.04 | 0.06 | -- | 98.7 |
| MZD-0014-R | 41.08 | 0.98 | 18.79 | 14.08 | 0.24 | 6.13 | 14.07 | 1.38 | 0.73 | 0.89 | 1.3 | 99.67 |
| APO-0044-LG | 47.49 | 0.12 | 24.58 | 6.5 | 0.15 | 6.12 | 12.72 | 1.06 | 0.05 | 0.01 | 1 | 99.8 |
| CAQ-0346-P | 46.34 | 0.56 | 15.69 | 8.28 | 0.06 | 7.38 | 7.28 | 1.03 | 1.28 | 0.05 | 11.8 | 99.75 |
| CLM-0377-P | 53.19 | 0.66 | 17.7 | 9.03 | 0.27 | 5.14 | 7.61 | 2.28 | 0.08 | 0.09 | 3.8 | 99.85 |
| CLM-0388-P | 52.94 | 0.66 | 19.96 | 4.68 | 0.1 | 5.15 | 8.32 | 4.8 | 0.45 | 0.12 | 2.6 | 99.78 |
| CLM-0445-P | 64.47 | 0.37 | 16.68 | 4.89 | 0.13 | 2.31 | 5.35 | 3.14 | 0.28 | 0.09 | 2.2 | 99.91 |
| APO-0036-P | 58.14 | 0.35 | 18.73 | 3.33 | 0.11 | 1.54 | 6.34 | 3.62 | 1.3 | 0.19 | 6.2 | 99.85 |
| APO-0054-LG | 61.65 | 0.52 | 16.87 | 6.57 | 0.15 | 2.7 | 4.91 | 3.86 | 0.18 | 0.13 | 2.3 | 99.84 |
| APO-0065-LT | 54.39 | 0.55 | 17.44 | 7.91 | 0.19 | 3.22 | 7.25 | 2.81 | 1.77 | 0.25 | 4 | 99.78 |
| APO-0069-LT | 58.44 | 0.48 | 16.85 | 6.93 | 0.21 | 2.8 | 6.58 | 3.48 | 0.9 | 0.16 | 3 | 99.83 |
| CLM-0436-L | 65.75 | 0.34 | 16.19 | 4.4 | 0.19 | 2.05 | 3.87 | 3.98 | 0.62 | 0.12 | 2.4 | 99.91 |
| CLM-0480-P | 58.06 | 0.75 | 17.12 | 8.03 | 0.09 | 3.25 | 1.2 | 1.46 | 2.18 | 0.13 | 7.5 | 99.77 |
Tabla 3 Análisis geoquímicos de elementos traza y tierras raras en ppm, en muestras del Stock de Timbiquí.
| Muestra | Ba | Ce | Co | Cs | Cu | Dy | Er | Eu | Ga | Gd | Hf | Ho | La | Lu | Nb | Nd | Ni | Pb |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| APO-0056-LG | 71.0 | 7.6 | 19.1 | <0.1 | 52.8 | 2.6 | 1.9 | 0.6 | 13.4 | 2.1 | 1.4 | 0.7 | 3.1 | 0.3 | 0.6 | 6.0 | 8.9 | 0.2 |
| APO-0049-LG | 348.0 | 10.6 | 10.3 | 1.0 | 37.9 | 2.2 | 1.4 | 0.6 | 14.9 | 1.9 | 1.3 | 0.5 | 5.0 | 0.2 | 1.6 | 5.7 | 3.7 | 0.6 |
| APO-0107-D | 91.0 | 10.4 | 15.0 | 0.4 | 26.6 | 2.5 | 1.8 | 0.7 | 13.4 | 2.3 | 1.8 | 0.6 | 4.4 | 0.3 | 0.7 | 7.9 | 2.1 | 0.3 |
| CAQ-0416-P | 162.0 | 15.2 | 32.1 | 0.3 | 237.4 | 3.1 | 1.5 | 1.3 | 16.5 | 3.6 | 0.5 | 0.6 | 6.9 | 0.2 | 0.7 | 14.8 | 8.4 | 0.6 |
| MZD-0010-L | 510.0 | 19.4 | 14.4 | 0.2 | 27.1 | 2.7 | 1.9 | 1.0 | 16.3 | 3.4 | 1.6 | 0.6 | 9.5 | 0.3 | 1.7 | 12.2 | 0.3 | 0.2 |
| MZD-0022-R | 358.0 | 18.3 | 3.8 | <0.1 | 9.4 | 2.3 | 1.4 | 0.7 | 12.9 | 2.1 | 2.3 | 0.5 | 10.8 | 0.3 | 2.5 | 10.9 | 0.5 | 0.2 |
| MZD-0014-R | 114.0 | 15.2 | 27.0 | 0.2 | 54.2 | 4.7 | 2.2 | 1.4 | 17.3 | 4.8 | 0.7 | 0.8 | 6.4 | 0.3 | 0.7 | 14.6 | 1.0 | 0.4 |
| APO-0044-LG | 51.0 | 1.8 | 19.3 | <0.1 | 3.1 | 0.2 | 0.2 | 0.2 | 15.1 | 0.1 | 0.2 | 0.1 | 1.3 | 0.0 | <0.1 | 1.0 | 2.7 | 0.5 |
| CAQ-0346-P | 41.0 | 5.5 | 24.7 | 0.4 | 50.5 | 1.9 | 1.3 | 0.5 | 12.8 | 1.6 | 0.9 | 0.4 | 1.8 | 0.2 | 0.3 | 4.8 | 22.9 | 0.7 |
| CLM-0377-P | 25.0 | 5.1 | 22.3 | <0.1 | 137.9 | 3.0 | 2.0 | 0.6 | 14.0 | 2.3 | 1.1 | 0.7 | 1.9 | 0.3 | 0.3 | 5.2 | 5.7 | 1.0 |
| CLM-0388-P | 113.0 | 5.5 | 16.1 | 0.1 | 14.8 | 2.6 | 1.5 | 0.6 | 13.3 | 1.9 | 1.0 | 0.5 | 2.4 | 0.3 | 0.3 | 4.3 | 74.1 | 0.9 |
| CLM-0445-P | 98.0 | 7.4 | 7.5 | 0.1 | 22.4 | 2.0 | 1.3 | 0.5 | 12.9 | 1.7 | 1.5 | 0.5 | 3.5 | 0.3 | 0.6 | 5.7 | 4.7 | 0.7 |
| APO-0036-P | 334.0 | 24.7 | 4.7 | 0.5 | 17.2 | 2.1 | 1.2 | 0.9 | 17.1 | 2.4 | 1.9 | 0.4 | 13.1 | 0.2 | 3.4 | 15.7 | 0.9 | 4.5 |
| APO-0054-LG | 215.0 | 11.3 | 12.4 | <0.1 | 60.8 | 2.7 | 1.7 | 0.7 | 14.0 | 2.5 | 1.6 | 0.6 | 5.0 | 0.3 | 0.8 | 7.6 | 5.5 | 0.4 |
| APO-0065-LT | 328.0 | 23.4 | 15.8 | 0.2 | 60.9 | 3.1 | 2.0 | 1.0 | 15.7 | 3.4 | 1.6 | 0.6 | 10.6 | 0.3 | 1.5 | 15.6 | 3.5 | 1.4 |
| APO-0069-LT | 195.0 | 15.1 | 12.3 | <0.1 | 44.0 | 2.1 | 1.7 | 0.8 | 13.2 | 2.5 | 1.6 | 0.5 | 7.2 | 0.3 | 1.2 | 11.1 | 3.4 | 1.1 |
| CLM-0436-L | 160.0 | 10.0 | 4.3 | 0.6 | 4.5 | 2.0 | 1.6 | 0.6 | 13.1 | 1.9 | 1.8 | 0.5 | 4.3 | 0.3 | 0.9 | 7.8 | 2.7 | 0.8 |
| CLM-0480-P | 427.0 | 27.3 | 19.1 | 2.3 | 110.5 | 2.7 | 1.7 | 0.9 | 16.4 | 3.0 | 3.1 | 0.6 | 13.6 | 0.3 | 5.1 | 13.2 | 44.4 | 6.7 |
| CAQ-0415R | 402.0 | 20.4 | -- | 0.4 | -- | 2.7 | 1.5 | 1.4 | -- | 3.2 | -- | 0.6 | 10.4 | 0.3 | 1.5 | 13.7 | -- | 0.4 |
| CAQ-0420R | 166.0 | 15.7 | -- | 0.7 | -- | 3.2 | 1.9 | 0.9 | -- | 3.3 | -- | 0.7 | 6.6 | 0.3 | 1.9 | 11.7 | -- | 0.7 |
| CLM-0447-L | 91.0 | 9.4 | 22.5 | <0.1 | 81.0 | 2.8 | 1.8 | 0.7 | 14.0 | 2.4 | 1.6 | 0.6 | 3.9 | 0.3 | 0.8 | 7.3 | 6.3 | 0.4 |
| APO-0091-P | 656.0 | 26.3 | 9.1 | 0.4 | 52.9 | 2.8 | 1.9 | 1.1 | 16.7 | 3.1 | 2.6 | 0.6 | 14.4 | 0.3 | 3.6 | 13.9 | 1.3 | 0.2 |
| CDG-0255-P | 27.0 | 9.0 | 17.6 | 0.1 | 1.7 | 3.2 | 1.9 | 0.9 | 14.1 | 2.7 | 1.8 | 0.7 | 3.2 | 0.3 | 1.1 | 7.4 | 4.2 | 0.2 |
| CDG-0256-P | 20.0 | 5.4 | 25.8 | 0.1 | 3.4 | 2.1 | 1.5 | 0.6 | 12.8 | 1.8 | 1.0 | 0.5 | 2.2 | 0.2 | 0.7 | 5.0 | 4.2 | 0.1 |
Tabla 4 Edades isotópicas compiladas de la literatura del evento magmático Paleoceno - Eoceno, Cordilleras Occidental y Central de Colombia.
| Muestra | Longitud | Latitud | Nombre Unidad | Edad (Ma) | Método | Tipo de roca | Referencia |
|---|---|---|---|---|---|---|---|
| Cordillera Central | |||||||
| DV64 | -74.95536 | 5.98236 | Batolito Antioqueño | 58.9± 9.6 | U/Pb Cir | Granodiorita | Villagómez et al., 2011 |
| DV70_2 | -75.42586 | 6.97078 | Batolito Antioqueño | 58.1± 5.2 | U/Pb Cir | Granodiorita | Villagómez et al., 2011 |
| DV148 | -75.38492 | 6.42119 | Batolito Antioqueño | 53± 4.8 | U/Pb Cir | Granodiorita | Villagómez et al., 2011 |
| GR-II-1 | -74.91191 | 6.51065 | Batolito Antioqueño | 59.2± 1.2 | U/Pb Cir | Tonalita biotítica | Leal Mejía, 2011 |
| G-29 | -74.92009 | 6.51047 | Batolito Antioqueño | 60.7± 1 | U/Pb Cir | Tonalita biotítica | Leal Mejía, 2011 |
| R-6834 | -74.75273 | 56.49891 | Stock de Caracolí | 60.1±1.2 | U/Pb Cir | Cuarzodiorita biotítica | Leal Mejía, 2011 |
| WR-215 | -75.29833 | 5.76629 | Batolito de Sonsón | 58.6± 1 | U/Pb Cir | Granodiorita biotítica | Leal Mejía, 2011 |
| WR-216 | -75.203548 | 5.66734 | Batolito de Sonsón | 55.8± 0.9 | U/Pb Cir | Leucogranito | Leal Mejía, 2011 |
| WR-218 | -75.19388 | 5.65845 | Batolito de Sonsón | 60± 0.9 | U/Pb Cir | Leucogranito | Leal Mejía, 2011 |
| 07DV156 | -75.300009 | 5.75369 | Batolito de Sonsón | 59± 0.5 | U/Pb Cir | Granito | Cochrane, 2013 |
| BSS-5 | -75.298104 | 5.75591 | Batolito de Sonsón | 60.7± 1.4 | U/Pb Cir | Cuarzo diorita | Ordóñez-Carmona et al., 2011 |
| 69 | -75.32995 | 5.95556 | Stock de La Unión | 59± 12 | U/Pb Cir | Tonalita | Restrepo et al., 1991 |
| CH10 | -74.98982 | 5.14893 | Stock el Hatillo | 55.0±0.7 | U/Pb Cir | -- | Bustamante et al., 2016 |
| CH19 | -74.99923 | 5.27350 | Stock el Hatillo | 53.7± 2.6 | U/Pb Cir | -- | Bustamante et al., 2016 |
| CH21B | -75.00382 | 5.29017 | Stock el Hatillo | 53.6± 1.7 | U/Pb Cir | -- | Bustamante et al., 2016 |
| CH5 | -74.99378 | 5.21960 | Stock el Hatillo | 53.8± 0.3 | U/Pb Cir | -- | Bustamante et al., 2016 |
| CH8 | -74.99125 | 5.17650 | Stock el Hatillo | 53.2±1.3 | U/Pb Cir | -- | Bustamante et al., 2016 |
| GCC-11 | -74.97359 | 5.17379 | Stock el Hatillo | 54.6±0.7 | U/Pb Cir | Cuarzodiorita | Bayona et al., 2012 |
| CB1 | -75.14055 | 4.73725 | Batolito El Bosque | 55.1± 2.1 | U/Pb Cir | -- | Bustamante et al., 2016 |
| DV163 | -74.97822 | 5.17877 | Plutón de El Hatillo | 55.3±5.4 | U/Pb Cir | Granito | Villagómez et al., 2011 |
| DV04 | -74.97539 | 4.78339 | Batolito de Ibagué | 59.2± 20.2 | U/Pb Cir | Gabrodiorita | Villagómez et al., 2011 |
| DV05 | -75.26814 | 4.40769 | Batolito de Ibagué | 59.8± 16.8 | U/Pb Cir | Granodiorita | Villagómez et al., 2011 |
| ER-1 | -75.09912 | 6.48492 | Intrusivo pórfido El Rayo | 59.9± 0.9 | U/Pb Cir | Pórfido aurífero | Leal Mejía, 2011 |
| SG-002 | -75.85156 | 6.59481 | Batolito de Sabanalarga | 50.28± 16.6 | U/Pb Cir | -- | ANH-UCaldas, 2011 |
| IGM-156841 | -75.41334 | 4.18053 | Pórfido de Río Manso | 58± 10 | U/Pb Cir | Cuarzodiorita | Núñez, 1981 |
| DV30 | -76.18066 | 3.90294 | Batolito de Buga | 44.8±8.4 | U/Pb Cir | Tonalita | Villagómez et al., 2011 |
| BSB1 | -76.15173 | 3.34943 | Batolito de Santa Barbara | 58.9±0.4 | U/Pb Cir | Tonalita | Ordóñez-Carmona et al., 2011 |
| BSB2 | -76.142849 | 3.37582 | Batolito de Santa Barbara | 58.4±0.8 | U/Pb Cir | Tonalita | Ordóñez-Carmona et al., 2011 |
| 10RC73 | -76.14499 | 3.36045 | Batolito de Santa Barbara | 60.2±1.1 | U/Pb Cir | Granito | Cochrane, 2013 |
| GCC-15 | -75.40529 | 5.03409 | Granodiorita Manizales | 59.8±0.7 | U/Pb Cir | Tonalita | Bayona et al., 2012 |
| R-6834 | -74.75273 | 6.49892 | Stock de Caracolí | 60.1± 1.2 | U/Pb Cir | Cuarzodiorita biotítica | Leal Mejía, 2011 |
| Cordillera Occidental | |||||||
| CLM-0376-R | -77.41746 | 2.50920 | Diques Dioríticos | 45.6± 0.6 | U/Pb Cir | Tonalita | ANH-GRP, 2014 |
| CLM-0447-P | -77.51078 | 2.41357 | Diques Dioríticos | 44.7± 0.3 | U/Pb Cir | Diques dioríticos | ANH-GRP, 2014 |
| CDG-0255-P | -77.58847 | 2.06410 | Intrusivo Muchinca | 50.3± 0.7 | U/Pb Cir | Cuarzodiorita - cuarzomonzonita | ANH-GRP, 2014 |
| APO-0056-LG | -77.24070 | 2.74230 | Stock de Timbiquí | 46.8± 0.7 | U/Pb Cir | Cuarzodiorita | ANH-GRP, 2014 |
| WR-238 | -76.52136 | 6.726582 | Batolito de Mandé | 44.6± 0.9 | U/Pb Cir | Pórfido tonalítico | Leal Mejía, 2011 |
| WR-237 | -76.518721 | 6.723443 | Batolito de Mandé | 45.3± 1.2 | U/Pb Cir | Pórfido tonalítico | Leal Mejía, 2011 |
| XX | -77.404248 | 8.524412 | Batolito de Acandí | 49.5 | U/Pb Cir | -- | Cardona et al., 2018 |
| X | -77.592307 | 2.045487 | Stock de Timbiquí | 44.3 | U/Pb Cir | -- | Cardona et al., 2018 |
MARCO TECTÓNICO REGIONAL
La evolución tectónica del centro y sur del segmento Andino durante el Cretácico Tardío y Cenozoico, fue caracterizada por la evolución de un sistema de subducción entre las placas Farallón y Sudamérica, que dio lugar a magmatismo en la placa subducente y continuos apilamientos tectónicos bajo un régimen transpresivo-compresivo (McCourt et al., 1984; Aspden et al., 1987; Grosser, 1989; Duque-Caro, 1990; Estrada, 1995; Cediel et al., 2003; Pindell y Kennan, 2009; López, 2009; Lissina, 2005; Ramos, 2009; Whattam y Stern, 2014; Montes et al., 2012; Bayona et al., 2012; Pindell et al., 2005; Villagómez et al., 2011; Cardona et al., 2018; Montes et al., 2019; Zapata-García y Rodríguez-García, 2020) (Figura 1). Presentamos a continuación una breve perspectiva histórica de cómo han evolucionado las ideas sobre la tectónica de la margen pacífica colombiana para el Cretácico Tardío y Cenozoico.
Barrero (1979) sugiere para el Cretácico Tardío el predominio de las fases extensionales en el borde de la placa Sudamericana, bajo un ambiente de subducción. Dicho régimen tectónico dio lugar a la formación de fosas y pilares tectónicos y al ascenso fisural de magmas toleíticos bajos en potasio, actualmente catalogados como tipo meseta oceánica. Este vulcanismo basáltico es parte del Grupo Diabásico y hace parte del basamento de las cuencas “onshore” Tumaco y Chocó (Figura 1a) de la margen pacífica colombiana (Barrero, 1979).
En contraposición a las ideas de Barrero (1979), autores como McCourt et al. (1984) y Grosser (1989) proponen varias acreciones sucesivas de terrenos alóctonos con basamento oceánico y continental en el borde occidental colombiano, en un sistema de subducción. Los terrenos acrecionados al continente durante el Cretácico Tardío fueron continentales, mientras que los terrenos oceánicos fueron acrecionados durante el Mioceno medio (Etayo-Serna et al., 1983; Restrepo y Toussaint, 2020). Durante el Mioceno Tardío fue acrecionado el Bloque Panamá-Chocó, reconocido por Duque-Caro (1990) como un terreno exótico que no tiene ninguna afinidad genética con Suramérica (Figura 1b). Presenta anomalías gravimétricas positivas asociadas a cuerpos ígneos de origen oceánico y está compuesto por tres elementos litotectónicos, el Arco de Dabeiba y Baudó, la Cuenca Atrato-Chucunaque y la Zona de Deformación de Istmina (Duque-Caro, 1990). El límite sur del Bloque Panamá-Chocó es la Zona de Deformación de Istmina y al oriente podría ser la zona de falla de Uramita (Figura 1a). Con respecto a la falla de Uramita Duque-Caro (1990) comenta: “la naturaleza de esta zona de falla es aún muy difícil de precisar debido a la falta de datos estructurales más detallados. Sin embargo, la ocurrencia exótica de los foraminíferos paleocenos, originalmente provenientes del norte, junto con la naturaleza tectónica del límite oriental del Arco de Dabeiba, sugieren que el lineamiento de Uramita es el contacto entre el Arco de Dabeiba de afinidad centroamericana y el noroccidente suramericano”. Estrada (1995) se une a la interpretación de Duque-Caro (1990) de terrenos alóctonos acrecionados al NW de Colombia durante el Mioceno (terrenos Gorgona, Chocó y Cordillera Occidental), basando su interpretación de análisis paleomagnéticos realizados en basaltos almhoadillados de origen oceánico y en plutones del Terciario tardío (6-10 Ma).
Un trabajo relevante es el de Cediel et al. (2003), quienes presentan una síntesis de los trabajos que tratan sobre el margen pacífico colombiano, y que fueron realizados en las tres décadas anteriores a la publicación de su trabajo. De acuerdo con esos autores, la porción noroeste de la margen colombiana consiste en el Terreno Chocó, el cual contiene los terrenos de Cañas Gordas y Baudó, de edad cretácica, que posteriormente fueron acrecionados durante el Mioceno medio al segmento norte de la margen pacífica colombiana. Durante el Paleoceno-Eoceno en el Terreno Chocó y el terreno Cañas Gordas, se instalaron los intrusivos de Mandé y Acandí (Figura 1b) que tienen una afinidad oceánica, determinada a partir análisis geoquímicos. Estos cuerpos intrusivos no presentan relación con los batolitos emplazados en la Cordillera Central (Batolito el Bosque, Granodiorita Manizales y stock El Hatillo), debido a que dichos batolitos de la Cordillera Central tienen una afinidad autóctona. La acreción del terreno Cañas Gordas comenzó en el Mioceno como una colisión altamente destructiva con una cinemática moderadamente destral oblicua. Dos evidencias reportadas por Cediel et al. (2003) sobre la colisión destructiva son: 1. “Las estructuras asociadas con la colisión del terreno Cañas Gordas deforman el terreno Gorgona y el sistema de fallas Buenaventura, que estaban en su lugar a principios del Mioceno” 2. “La falla de Garrapatas parece controlar la sedimentación del Mioceno al Holoceno a lo largo de la margen del Pacífico”. Así entonces, con esas evidencias y basado en la sugerencia dada por Duque-Caro (1990), los autores se adhieren al modelo de acreción del Bloque Chocó en el Mioceno.
Más recientemente Montes et al. (2012) definieron varios eventos magmáticos: (1) arco magmático Cretácico Tardío hasta el Eoceno medio (58-39 Ma), asociado a subducción, (2) gap magmático desde el Eoceno tardío hasta el Oligoceno tardío (38-28 Ma), (3) reinicio del magmatismo a los 25 Ma. Estos eventos fueron parte de un proceso de colisión y acreción del borde de la Placa Caribe al continente, durante el Mioceno (Montes et al., 2019). Así mismo, Cardona et al. (2018) definieron dos arcos magmáticos Paleocénicos, no-cogenéticos, de afinidades continental y oceánica. Uno de ellos está ubicado en la parte norte de la Cordillera Occidental de Colombia y en Panamá, formado de rocas plutónicas y volcánicas con edades de 60-45 Ma, y su origen fue asociado a un ambiente tectónico tipo Aleutiano. El segundo es el arco de Timbiquí, que probablemente hizo parte de un terreno alóctono acrecionado durante el Eoceno tardío-Oligoceno (Echeverri et al., 2015; Cardona et al., 2018). Actualmente, desde el punto de vista geomorfológico, los cuerpos intrusivos están ubicados en la topografía positiva al occidente de Colombia, que es denominada la Cordillera Occidental. Desde el punto de vista genético y geológico, estos cuerpos intrusivos han sido atribuidos como parte del Bloque Panamá-Chocó. Por lo tanto, estos cuerpos son parte de la Cordillera Occidental y, según los modelos propuestos en las últimas cuatro décadas, también del Bloque Panamá-Chocó.
Zapata-García y Rodríguez-García (2020) muestran evidencias de la continuidad del arco magmático Paleoceno-Eoceno. Adicionan al Bloque Panamá-Chocó la Formación Timbiquí y la Tonalita de Napi (Figura 1b), ubicadas al sur de la Cordillera Occidental de Colombia. A partir de análisis geoquímicos, geología de campo, geocronología y petrografía, muestran que rocas volcánicas y plutónicas que constituyen el Arco Panamá-Chocó norte (Complejo Santa Cecilia-La Equis y Batolito de Acandí), central (Complejo Santa Cecilia-La Equis y Batolito de Mandé) y sur (Formación Timbiquí y Tonalita de Napi) son similares en litología, geoquímica, composición y edad (Eoceno).
El magmatismo Paleoceno temprano en la Cordillera Central ha sido atribuido a procesos relacionados con la flotabilidad del slab de subducción de la placa Caribe (Bayona et al., 2012; Cardona et al., 2018). Análisis geocronológicos de U/Pb en cuerpos intrusivos y circones detríticos, indican el aporte de sedimentos de edad Paleoceno para las cuencas intramontanas. Se ha propuesto que la parada repentina del magmatismo a principios del Eoceno medio puede estar relacionado con la dificultad de la subducción de un grueso plateau oceánico, dada por el movimiento transcurrente entre la Placa Sudamericana y Caribe (Aspden et al., 1987; Bayona et al., 2012; Pindell et al., 2005; Villagómez et al., 2011; Cardona et al., 2018).
METODOLOGÍA
Se realizaron cinco dataciones isotópicas usando el método U/Pb en circones de cuerpos ígneos (Figura 2, Tabla S3 del suplemento electrónico), en cuerpos ígneos localizados a lo largo de la CO. Se utilizó la técnica de análisis espectrométrico de relaciones isotópicas en granos individuales con ablación por punto (“spot analysis”) utilizado un sistema de laser acoplado a un espectrómetro de masas (Laser Ablation-Multicollector-Inductively Coupled Plasma-Mass Spectrometry o LA-ICP-MS) (Tabla 1, Figura 3). Se realizaron también análisis geoquímicos en roca entera, incluyendo veintidós (22) muestras para elementos mayores y veinticuatro (24) para elementos traza. Los análisis de elementos mayores y traza se llevaron a cabo mediante el método ICP-MS (Inductively Coupled Plasma-mass spectrometry) en ACMELABS (Analytical Laboratories, Vancouver, Canadá). La precisión lograda con este procedimiento es de ±2 ppm y ±5 ppm para una concentración analítica entre 50 y 5 ppm, respectivamente. Para cada muestra, fueron tomados 0.2 g y se mezclaron con 1.5 g de flujo de LiBO2. Fueron calentados a 1050 °C durante 15 min. Posteriormente se disolvieron en HNO3 al 5 %. Patrones de calibración y blancos de reactivo fueron añadidos, y las soluciones se aspiraron en un espectrómetro de emisión ICP. Los límites de detección de los elementos mayores varían entre 0.001 % y 0.004 %, y los valores de los oligoelementos son de aproximadamente 0.1 ppm. Los elementos de tierras raras tienen límites de detección entre 0.01 y 0.05 ppm. Estos resultados se presentan libres de H2O+, H2O-, CO2 y normalizados a 100 %. Las muestras con pérdida por ignición mayor a 3.5 % fueron descartadas de los análisis debido a que estas pérdidas se consideran como un índice de alteración de la roca. Los resultados de los análisis geoquímicos se manejaron y procesaron utilizando el software GCD Kit 4.1 (Janousek et al., 2006) (Tablas 2 y 3).

Figura 2 Geocronología U/Pb circón. a) y b) Batolito de Mandé, c) Mandé - Santa Cecilia la Equis, d) y e) Batolito de Acandí. [número de la muestra], desviación estándar media ponderada (MSWD).

Figura 3 Distribución, geometría y tiempo de permanencia de las rocas de edad Paleoceno y Eoceno. Los análisis reportados son realizados a partir del método U/Pb en circón. En color rojo, se muestran las edades obtenidas en este trabajo. En color azul, los análisis U/Pb recopilados de la literatura (Núñez, 1981; Restrepo et al., 1991; Leal Mejía, 2011; ANH-UCaldas, 2011; Ordóñez-Carmona et al., 2011; Bayona et al., 2012; Villagómez y Spikings, 2013; ANH-GRP, 2013; Cochrane, 2013; Bustamante et al., 2016; Cardona et al., 2018).
Se recopilaron treinta y cinco (35) análisis geocronológicos U/Pb (Tabla 4), sesenta y siete (67) análisis de elementos mayores y cuarenta y cuatro (44) de elementos traza (Tablas S1 y S2), disponibles a partir de investigaciones publicadas y trabajos e informes cartográficos del Servicio Geológico Colombiano y la ANH.
En total se realizó el análisis de cuarenta (40) datos geocronológicos U/Pb, ochenta y nueve (89) análisis de elementos mayores, y sesenta y ocho (68) de elementos traza. Todos los análisis de óxidos mayores se recalcularon al 100% teniendo en cuenta la pérdida por ignición (LOI).
RESULTADOS
Se integraron los nuevos análisis con los obtenidos de la bibliografía, para construir una base de datos geocronológicos y geoquímicos de rocas plutónicas/volcánicas del Paleoceno - Eoceno de la región estudiada. Los resultados son presentados resaltando los dominios en la CO y en la CC, enfocándose en la distribución, migración y tiempo de permanencia.
Arco magmático en la Cordillera Occidental
Los cinco (5) nuevos análisis realizados en esta investigación en conjunto con los reportados por ANH-GRP (2014), muestran que las edades de cristalización (U/Pb en circón) del arco magmático de Acandí y Mandé son muy cercanas a los intrusivos denominados Stock de Timbiquí, Diques Dioríticos, Intrusivo de Muchinca y Tonalita de Napi, ubicados en el sector sur de la CO (Figuras 3 y 4). Las composiciones de esos cuerpos intrusivos incluyen tonalita, cuarzodiorita, diorita, gabro con cuarzo y gabro horblendítico para la Tonalita de Napi. El Stock de Timbiquí está constituido por gabros, dioritas cuarzosas, andesitas, basaltos (tipo lavas almohadilladas). Los Diques Dioríticos están compuestos de diorita y cuarzodiorita. El Intrusivo de Muchinca son cuarzodioritas y microgabros. Estas rocas presentan series toleíticas a calco-alcalinas como se observa en la Figura 4.

Figura 4 a) Localización general de las unidades litológicas del evento magmático Paleoceno - Eoceno. b) Diagrama de concentración de álcalis total vs. SiO2 (TAS). c) Diagrama de clasificación de SiO2 vs. K2O. d) Diagrama AFM (A: álcalis (Na2O+K2O), F: FeO+Fe2O3, M: MgO). Los valores se presentan en % en peso. Para las rocas distribuidas en la Cordillera Occidental, los análisis son representados por triángulos de distintos colores según la unidad geológica. Los cuadrados representan los análisis realizados en intrusivos localizados a lo largo de la Cordillera Central. Son incluidos los datos reportados en la literatura (Sillitoe et al., 1982; Restrepo et al., 1991; INGEOMINAS, 1979, 1999, 2002a, 2010b; López et al., 2006; Vallejo et al., 2011; Villagómez y Spikings, 2013; Bustamante et al., 2016).
Considerando todas las edades ubicadas sobre la CO se observa que abarcan de ~59.2 Ma a 44.3 Ma. De norte (límites con Panamá) a sur (límites con Ecuador), se identifican las siguientes unidades: 1. Batolito de Acandí, con edades de 49 Ma a 46 Ma. 2. Batolito de Mandé, con edades desde 59 Ma hasta 44 Ma. 3. El Stock de Timbiquí, Diques Dioríticos y el Intrusivo de Muchinca (entre 2° a 3° de latitud N) registran edades desde 50 Ma hasta 44 Ma. (Figuras 2 y 3; Tabla 2). Composicionalmente, estos cuerpos están constituidos por variaciones de gabros hasta granitoides. En el Batolito de Mandé se observan facies más básicas hacia los bordes, diques basálticos y localmente con inclusiones máficas. (INGEOMINAS, 2002a, 2002b, 2010a, 2010b). En primera instancia, este evento magmático Paleoceno - Eoceno presenta una predominancia de magmatismo intermedio a básico, desde gabro hasta granodiorita (Figura 4b). Los diagramas de Peccerillo y Taylor (1976) y AFM muestran que el conjunto de análisis es una tendencia progresiva, de series toleíticas hacia series calcoalcalinas. (Figura 4c).
Los elementos mayores se graficaron en relación a la variación de la concentración del SiO2 (Figura 5a). Las concentraciones de SiO2 presentan valores menores a 65 wt%, una tendencia a la disminución del SiO2 a mayores concentraciones en CaO, MgO, Al2O3 y FeOt; en contraste con el aumento de SiO2 para mayores concentraciones de Na2O y K2O. En los elementos traza (Ba, Ce, La, Rb, Sr, Zr) se observa una disminución progresiva y continua al decrecimiento de SiO2. El Ba tiene valores hasta de 600 ppm, el Ce y Rb hasta 30 ppm, el La hasta 15 ppm, el Sr hasta 500 ppm y el Zr hasta 100 ppm (Figura 5b).

Figura 5 a) Diagramas de variación en %wt de elementos mayoritarios vs SiO2. Obsérvese la variación en la concentración de los diferentes óxidos mayoritarios para los cuerpos intrusivos distribuidos en la CC y CO. b) Diagramas de variación de elementos traza (ppm) vs. SiO2 (%wt). Para los análisis realizados sobre intrusivos distribuidos sobre la Cordillera Occidental se dibujan triángulos, para los de la Cordillera Central se usan cuadrados y aparecen encerrados en zonas sombreadas.
Las edades y la composición química de los cuerpos magmáticos descritas en los párrafos anteriores sugieren fuertemente que hay continuidad del conjunto magmático, desde Acandí hasta el sector sur de la CO. Esta continuidad también fue propuesta por Zapata-García y Rodríguez-García (2020) desde un punto de vista geoquímico.
Arco magmático en la Cordillera Central
En la Cordillera Central (CC) las unidades plutónicas/volcánicas están concentradas principalmente a 4.5° de latitud N y la mayor expresión está dada por el Batolito Antioqueño (Aspden et al., 1987). Las edades de cristalización (U/Pb en circón) de los cuerpos distribuidos en la CC están registradas desde el Cretácico Tardío hasta 53 Ma. Sin embargo, en este trabajo se considerarán sólo las dataciones en el rango Paleoceno - Eoceno (Figuras 3 y 4). De norte a sur, se tienen los intrusivos: 1. Batolito Antioqueño, Intrusivo pórfido el Rayo, Stock de Caracolí y Batolito de Sabanalarga, con edades que van de 60.1 Ma a 50.2 Ma. 2. Stock de la Unión con una edad de 59 Ma. 3. Batolito de Sonsón que presenta edades de 60 Ma a 55.8 Ma. 4. Stock El Hatillo varia de 53.8 Ma a 53.3 Ma. 5. Granodiorita de Manizales con edades de 60.7 Ma a 59 Ma. 6. Batolito El Bosque presenta edades al sur del cuerpo intrusivo de 55.1 Ma hasta 59.2 Ma. 7. Batolito de Buga de 44.8 Ma y otros cuerpos como el Batolito de Ibagué con edades de 59.8 a 58 8 Ma y Batolito de Santa Bárbara de 60.2 Ma a 58.4 (Restrepo et al., 1991; Ordóñez-Carmona et al., 2011; Leal Mejía, 2011; ANH-UCaldas, 2011; Bayona et al., 2012; Villagómez y Spikings, 2013; Cochrane, 2013; Bustamante et al., 2016). Esos intrusivos están compuestos principalmente por granodiorita, tonalita, cuarzodiorita, cuarzomonsodiorita, granito y cuarzomonzonita. Localmente, algunos plutones contienen xenolitos de gneis y anfibolita, así como diques aplíticos (Núñez, 1981; INGEOMINAS, 2002a, 2002c, 2002d, 2002e). La composición en general es ácida y se clasifica en los campos de las granodioritas y granitos con afinidades calcoalcalinas (Figura 4).
La concentración del SiO2 está en rangos de 75 a 65 wt%. La concentración de Na2O y K2O muestra un incremento con una tendencia positiva con respecto al SiO2. Por el contrario, el Al2O3, FeO, CaO y MgO, disminuyen su concentración con el incremento del SiO2. Los elementos traza para estos cuerpos muestran una correlación positiva con respecto al contenido de SiO2. La concentración de Ba es de 600-1000 ppm, Ce de 30-70 pmm, La 15-40 ppm, Rb 20-60 ppm, Sr 500-700 ppm y Zr desde 100 hasta 150 pmm (Figura 5).
DISCUSIÓN
La discusión que se presenta aquí se restringe al registro del Paleoceno - Eoceno, los eventos magmáticos posteriores al Eoceno quedan fuera del objetivo de este trabajo.
Tiempo de permanencia y migración del arco magmático
El evento magmático en el rango Paleoceno - Eoceno está registrado por cuerpos intrusivos distribuidos en la CC que tienen edades de 60.7 Ma a 44.8 Ma (Figura 3 y 6), y representan un gran volumen de magma. Sobre la CO, en el arco de Panamá, hasta la parte centro-norte de la cordillera, representada por el Batolito de Acandí y el Batolito de Mandé, se tienen edades desde 59.2 Ma hasta 44.6 Ma y hacia el sur, en el sector de Tumaco, las edades de los cuerpos intrusivos varían de 50.3 Ma a 44.3 Ma (Figura 3 y 5).

Figura 6 Distribución de las edades de cristalización a partir del método U/Pb en circón del occidente colombiano, en las rocas ígneas del Paleoceno - Eoceno en función de la distancia a la fosa actual. La distancia es medida desde la fosa hacia los puntos de muestreo de manera perpendicular, teniendo en cuenta la geometría de las estructuras mayores (Cordillera Central y Occidental).
Los cuerpos intrusivos con edades entre 60.7 Ma y 53 Ma se emplazaron simultáneamente en ambas cordilleras Occidental y Central (Figura 6). Sin embargo, las composiciones y geoquímica son diferentes, mientras en la CC predominan las granodioritas con series principalmente calco-alcalinas, en la CO predominan composiciones de granodiorita con variaciones composicionales hacia los bordes de los cuerpos a facies básicas, e inclusiones basálticas (Batolito de Mandé) y adicionalmente, numerosos diques andesíticos y basálticos. La geoquímica en los cuerpos intrusivos de la CO abarca desde series toleíticas hasta calcoalcalinas (Figuras 4 y 5).
A partir de los 53 Ma y hasta los 44.3 Ma, la provincia magmática Paleoceno - Eoceno se concentra en la CO y en el flanco occidental de la CC, evidenciando la migración del arco magmático hacia el occidente (Figura 6).
Edades isotópicas muestran que la actividad magmática Paleoceno-Eoceno estuvo también presente en el Arco de Panamá. La actividad magmática inició alrededor de los 70 Ma y disminuyó después de los 38 Ma.
Cinemática y control geométrico
Para el análisis cinemático de la placa Farallón con respecto a la placa Sudamericana, considerando el lapso entre 60 Ma y 44 Ma, se tomó como referencia una dirección de convergencia promedio de 65° (Figura 7). Esta dirección es equivalente a la dirección de desplazamiento (25° de azimut) de la placa Farallón con respecto al norte actual (Figura 8). Lo anterior se basa en los trabajos de Martinod et al. (2010), el cual incorpora estudios de cinemática realizados a lo largo de Suramérica, y el de Liu et al. (2008), que se refiere a la zona más cercana al límite norte de Colombia (Figura 7).

Figura 7 a) Velocidad y ángulo de convergencia de la placa Farallón - Nazca - Coiba y la placa Sudamericana para el Cenozoico. Distribución de los máximos periodos de enfriamiento según Noriega-Londoño et al. (2019), para los cuerpos plutónicos/volcánicos ubicados en la Cordillera Central y la Cordillera Occidental.

Figura 8 a) Disposición geométrica de las rocas plutónicas con su eje mayor en una dirección NE representada por las líneas continuas de color verde. En el área norte, si se hace la corrección por la rotación horaria (30°) de los cuerpos intrusivos, en su azimut regional, el eje mayor de todos los cuerpos intrusivos de este evento magmático, son paralelos a las flechas rojas. La dirección de desplazamiento de la placa Farallón es simbolizada por la flecha roja, la cual es paralela a la dirección del esfuerzo compresivo máximo σ1. b) Partición del vector de desplazamiento de la placa Farallón. La línea punteada denominada X, representa la componente ortogonal del desplazamiento y la línea punteada denominada Y, la componente transcurrente del desplazamiento. La línea azul punteada hace referencia al límite de la margen continental en el Paleoceno. c) Elipse de deformación generada para σ1 paralelo al vector de desplazamiento en dirección NE (25° de azimut) durante el Paleoceno - Eoceno. Las flechas azules indican la dirección del esfuerzo principal mínimo σ3, la línea azul continua, la dirección esperada del desarrollo de fracturas extensionales y tensionales en segundo orden de deformación.
En términos de velocidad, la placa Farallón se desplazaba en promedio a 5 cm/año entre los 65 Ma y 30 Ma (Pardo-Casas y Molnar, 1987; Daly, 1989; Somoza y Chidella, 2005). Considerando un azimut de 25° para el vector de velocidad de desplazamiento de la placa Farallón con respecto a la margen de Suramérica (Figura 8), obtenemos que la componente ortogonal a la margen es de 2.1 cm/año, lo que representa el 10 % de la velocidad total de convergencia. La componente paralela a la margen es de 4.5 cm/año, que corresponde al 90 % de la velocidad total, como se observa en la Figura 8.
El ángulo de convergencia de la placa Farallón transfiere un esfuerzo compresivo principal máximo σ1 horizontal en dirección paralela al vector de desplazamiento de la placa Farallón y un esfuerzo σ3 horizontal en dirección ortogonal a σ1, entonces se promovería el desarrollo de fracturas extensionales con dirección NE, con rotaciones progresivas en sentido horario (Figura 8c). Bajo estas condiciones se desarrollarían zonas de menor presión, o de fracturas tensionales, que pudieron controlar el emplazamiento elongado de los cuerpos plutónicos sobre la CO.
A lo largo de la CO los intrusivos son alargados, con sus ejes mayores en dirección NE, paralelos al esfuerzo principal σ1 inferido y al vector de desplazamiento de la placa Farallón, como se muestra en la Figura 8a. Por esta razón, se puede inferir que la margen continental estuvo sometida a transtensión derecha como mecanismo de primer orden y se minimizan los efectos compresionales. Una explicación para la orientación anómala de algunos batolitos que no se ajustan al esquema descrito, como es el caso del batolito de Mandé, es que fueron rotados (rotación antihoraria) durante los pulsos de la fase Andina.
En contraste, para el mismo periodo Paleoceno - Eoceno, los intrusivos emplazados en la CC son batolitos y stock con geometrías que tienden a ser circulares en su proyección de planta. Esas geometrías están relacionadas a sistemas extensionales radiales, permitiendo un proceso de emplazamiento magmático con mayor tiempo de permanencia. El emplazamiento de cuerpos intrusivos durante el Paleoceno - Eoceno en ambas cordilleras estuvo controlado por diferencia en magnitud y orientación de las zonas de extensión. Este periodo de tiempo (60.7-53 Ma) está relacionado con altos volúmenes de emplazamiento magmático y a su vez, con los menores valores de velocidad de convergencia entre las placas Farallón y Sudamericana; en consecuencia, las bajas velocidades y la alta oblicuidad de la convergencia configuran la relación entre el desarrollo de fracturas y zonas de debilidad, y facilita el emplazamiento de los magmas.
Los regímenes tectónicos transtensional-extensional pueden ser desarrollados en diferentes panoramas de interacción entre las placas tectónicas. En sistemas de subducción, es controlado principalmente por la migración negativa (alejamiento) de la trinchera con respecto a la placa subducente (rollback negativo) y por condiciones oblicuas de convergencia. Este proceso hace que la velocidad de subducción sea más rápida que el índice de convergencia, lo que implica una disminución en el acople entre las placas y en la presión de confinamiento, promoviendo el desarrollo de sistemas extensional-transtensional de primer orden en la placa subducente (Daly, 1989; Doglioni et al., 1999, Doglioni et al., 2006a; Doglioni et al., 2007; Schellart, 2008a, 2008b). Otra consecuencia del rollback negativo es el aumento en el ángulo de subducción (Doglioni et al., 2006b; Schellart, 2008b), lo que induciría a una migración de la cuña astenosférica consecuentemente con su arco magmático (subplacado), hacia el borde de la placa subducente, sin embargo, el magmatismo tendría una firma principalmente de manto hidratado (Ramos, 2009).
El régimen transtensivo-extensional es también generado en un borde continental (zona de transición). El magmatismo emplazado allí es favorecido por fallas normales, que posiblemente actúan como zonas permeables para el ascenso de basaltos syn-extensionales, con afinidades oceánicas (MORB-type), que puede estar acompañado por sucesiones de sedimentos oceánicos (Decarlis et al., 2018). En la cercanía límite océano-continente la generación de fallas normales profundas favorece el emplazamiento del subplacado, compuesto de rocas mantélicas principalmente, con posteriores procesos de serpentinización (Whitmarsh et al., 1996; Boillot y Froitzheim, 2001; Desmurs et al., 2001; Decarlis et al., 2018; Ribes et al., 2019; Chloé et al., 2017). Adicionalmente, se desarrollan estructuras extensionales simétricas y asimétricas, que localmente pueden presentar componentes de transcurrencia de segundo y tercer orden de jerarquía, y que controlan la sedimentación (Zhang et al., 2019; Osmundsen y Ebbing, 2008; Chloé et al., 2017). En el caso del borde continental colombiano durante el Paleoceno-Eoceno, o incluso desde el Cretácico (Jaramillo et al., 2017), la interacción de las placas Farallón y Sudamericana fue principalmente transcurrente, como se muestra en la Figura 7.
Un sistema de subducción oblicua con acreción de terrenos como el descrito en el marco tectónico regional habría dejado evidencias características de dicho proceso, a saber: 1. Un evento magmático Paleoceno-Eoceno en la placa subducente, con la geometría de sus ejes mayores dispuestos en una dirección perpendicular a la dirección de convergencia, que sería continuo hasta la acreción miocénica de terrenos. 2. Los límites de los terrenos acrecionados durante el Mioceno no han sido definidos y no se ha documentado la deformación generada a partir de esa acreción; así mismo, no se evidencia un prisma acrecional al oriente del bloque Panamá-Chocó. 3. No se tiene identificada la fosa del sistema de subducción que trasportó el bloque Panamá-Chocó.
Considerando los problemas expuestos en el párrafo anterior, de un modelo de acreción de terrenos, en este trabajo favorecemos un modelo alternativo retomando el postulado original de Barrero (1979). En este modelo se considera un borde de margen continental fracturado en extensión-transtensión, que facilitó y controló el emplazamiento del evento magmático desde el Cretácico Tardío como lo describió Barrero (1979), hasta el Paleoceno - Eoceno como es reportado aquí (Figura 9).

Figura 9 Representación paleogeográfica del emplazamiento de los cuerpos plutónicos y volcánicos durante el Paleoceno-Eoceno. Las rocas ubicadas al norte de la Proto-Falla de Istmina y al occidente de la falla de Uramita aparecen en su posición en el Paleoceno-Eoceno. Se reconstruyó su posición a partir de la corrección por la rotación generada en la deformación andina. Obsérvese cómo los cuerpos ígneos del Cretácico Superior, principalmente volcánicos (basaltos y dacitas), y las fallas principales Romeral, Cauca y Uramita, también presentan su eje principal paralelo al vector de convergencia. Las líneas punteadas en negro representan estructuras extensionales, interpretadas en líneas símicas y secciones estructurales (Grajales et al., 2016). El Triásico-Jurásico está conformado por plutones de sienogranitos que varían de tonalitas y cuarzomonzonitas a cuarzomonzodioritas. El Cretácico Inferior presenta plutones de gabro y peridotita serpentinizados, los volcánicos son basaltos y andesitas intercalados con lodolitas carbonosas y arenitas (Complejo Quebradagrande). El Cretácico Superior es conformado por granodioritas, tonalitas y cuarzodioritas, localmente gabros, que varía a cuarzogabros; las rocas volcánicas son compuestas por basaltos toleíticos, doleritas, picritas intercaladas con tobas básicas. El Paleoceno-Eoceno son tonalitas, granodioritas, cuarzodioritas y localmente gabros, las rocas volcánicas son basaltos intercalados con chert, lodolitas calcáreas y arenitas calcáreas, ocasionalmente calizas. Durante el Paleoceno-Eoceno la sedimentación en la plataforma está dada por secuencias continuas de shale marino con desarrollo de carbonatos. Localmente, hay abanicos deltaicos progradando hacia la plataforma en extensión.
Un esquema de la evolución y emplazamiento de los magmas en el occidente colombiano, es presentado en la Figura 10, donde se muestra el desarrollado en un margen continental fracturado, previo a la evolución del sistema de subducción. Durante el lapso 60.7-53 Ma se considera una dirección de convergencia de 25° de azimut y una velocidad de 5 cm/año, y para el lapso 53 Ma y 44 Ma una dirección de convergencia de 40° de azimut y una velocidad de 15-20 cm/año.

Figura 10 Reconstrucción esquemática de la evolución y emplazamiento de rocas plutónicas del evento magmático durante el Paleoceno-Eoceno en el occidente colombiano. El gráfico superior representa el emplazamiento durante los 65-50 Ma con una dirección de convergencia de 25° de azimut y una velocidad de 5 cm/a. Se muestra la elipse de deformación correspondiente al tensor de esfuerzos mayor, que representa la dirección de convergencia (flecha roja). Las líneas verdes dentro de la elipse simbolizan la tendencia de las fracturas tensionales dentro del sistema. El gráfico inferior representa el emplazamiento durante 50 a 33 Ma. con una dirección de convergencia de 40° de azimut y una velocidad de 15-20 cm/a. Los cuerpos del evento magmático (azules), migran hacia el occidente durante este periodo.
Implicaciones geoquímicas en la migración y el emplazamiento de intrusivos
Los parámetros físico-químicos que dominan el ascenso de los magmas son su densidad y su viscosidad. Así mismo, el mecanismo principal por el cual se genera el ascenso de magma es la propagación de fracturas, las cuales permiten una transferencia rápida de magma y el establecimiento de una red de estructuras que vuelve permeable a la corteza y que facilita el ascenso (Brown, 2004). El incremento de la densidad en los magmas es dado por la adición de FeO, MgO, TiO2 y CaO al fundido, mientras que la adición de álcalis y volátiles produce el efecto contrario (Spera, 2000; Stolper y Walker, 1980). El magmatismo que migró hacia el occidente colombiano durante el Paleoceno-Eoceno es de magmas principalmente básicos e intermedios con afinidades toleíticas, que varían progresivamente a calcoalcalinos. Comparados con los cuerpos distribuidos en la CC contienen mayores concentraciones de CaO, FeO, MgO y una menor concentración de los elementos traza (Rb, Sr, La, Zr, Ce y Ba) (Figura 4 y 5). Por lo anterior, podemos inferir que los magmas de la CO tuvieron densidad relativamente alta y viscosidad baja, por eso requirieron de fracturas tensionales que facilitaron su ascenso rápido en el borde continental. Progresivamente esos magmas fueron adquiriendo afinidades calcoalcalinas (Figura 8 y 9). El espesor de la corteza en la que se emplazan los magmas juega un papel muy importante. Chiaradia (2015), Chapman (2015) y Profeta et al. (2015) establecen campos en el gráfico Sr/Y vs. MgO que corresponden a espesores corticales promedio (Figura 11). El espesor de la corteza donde se emplazó el evento magmático Paleoceno - Eoceno es inferido por las variaciones en la razón Sr/Y y la concentración de MgO de 47 muestras tomadas en intrusivos ubicadas en la CO. Las muestras presentan valores de Sr/Y que varían a lo largo de todo el eje vertical del gráfico, aunque la mayor concentración se ubica en el rango donde Chiaradia (2015) establece un espesor de corteza entre <20-30 km (Figura 11). Las muestras que presentan mayores valores de Sr/Y, estarían relacionadas a las fases finales de emplazamiento. En la CC, por el contrario, los análisis marcan claramente un aumento del valor Sr/Y con la disminución del MgO, lo que permite interpretar espesores de la corteza continental mayores a 30 km (Figura 11). Esos espesores mayores son consistentes con las composiciones ácidas de series calco-alcalinas provenientes de magmas más evolucionados y con mayor tiempo de permanencia.

Figura 11 Gráficas Sr/Y vs. MgO de 89 análisis recopilados de la CC y CO, que muestra espesores promedio de la corteza según Chiaradia (2015), Chapman (2015) y Profeta et al. (2015), en la que se emplazaron estos cuerpos. a) Análisis realizados en cuerpos intrusivos de la CO (57 muestras). Los triángulos en celeste representan los cuerpos ubicados en el sector sur de la cordillera, como son el Stock de Timbiquí, la Tonalita de Napí, Diques Dioríticos e intrusivo de Munchica. Los triángulos rojos, los intrusivos ubicados al norte de la cordillera, Batolito de Mande y Batolito Acandí. b) Análisis de cuerpos intrusivos ubicados a 5° N sobre la CC (32 muestras).
Así entonces, el emplazamiento de los magmas de la CO está condicionado por su mayor densidad y baja viscosidad, en un medio de propagación de fracturas controladas por un régimen tectónico transtensional, desarrollado en una corteza de borde continental con un espesor < 25-30 km.
EVOLUCIÓN HACIA UN EQUILIBRIO DINÁMICO DEL MARGEN PACÍFICO COLOMBIANO
Como se ilustra en la Figura 7, entre 40 Ma y 26 Ma tuvo lugar una disminución significativa en la velocidad de convergencia de la placa Farallón, de 20 cm/año a 5 cm/año, con un ángulo promedio de convergencia de 30° de azimut, (Pardo-Casas y Molnar, 1987; Daly, 1989; Maloney et al., 2013; Sdrolias and Müller, 2006; Somoza y Chidella, 2005; Martinod et al., 2010; Liu et al., 2008). El cambio en la velocidad de la placa Farallón y la disminución de la velocidad total de la placa Sudamericana reduce la velocidad de convergencia y se espera una disminución en la eficiencia del acople entre las placas. Al reducir el acople, se disminuye la carga tectónica y, en consecuencia, el orógeno reduce su tasa de levantamiento. El proceso de equilibrio de masas da lugar a un aumento relativo en la remoción por procesos gravitacionales de la masa orogénica y de la sedimentación en los sistemas cuencales intramontanos. Este proceso de relajación del esfuerzo tectónico regional se denomina en este estudio Equilibrio Dinámico, de acuerdo al principio físico de d’Alembert.
En el ambiente intra-continental Andino colombiano, el equilibrio dinámico es manifestado por la discordancia progresiva del Eoceno medio/tardío a lo largo de las cuencas sedimentarias del Valle del Cauca Patía, Valle del Magdalena, Sinú San Jacinto entre otras, sumada a procesos masivos de peneplanización para los altos topográficos ubicados en el norte de Colombia relacionados a los procesos de la discordancia (Barrero, 1979; Cediel et al., 1998). Es también reportado en las cuencas de Yucatán (México), en Venezuela y el Caribe colombiano (James, 2009). Durante ese periodo, se reactivaron los piedemontes cordilleranos, promoviendo el depósito de secuencias de abanicos de las Formaciones Chimborazo y Peña Morada, seguida por la Formación Guachinte-Ferreira en la Cuenca Cauca Patía (Barrero et al., 2007; Ecopetrol, 2000), Formaciones San Cayetano y Maco en las cuencas Valle Inferior del Magdalena (VIM) (Pozo Betulia-1 y el Pozo ANH-Plato-1X) y Sinú San Jacinto (SSJ) (Barrero et al., 2007; ANH, 2017), y para la cuenca Valle Medio del Magdalena (VMM) en las formaciones Lizama, Esmeralda, La Paz (Gómez et al., 2005).
En los ambientes de plataforma continental, el equilibrio dinámico es evidenciado por la distribución de secuencias de calizas arrecifales y su variación de facies laterales, en toda la margen Atlántica y Pacífica colombiana. Los ambientes de plataformas de carbonatos predominaron durante el Eoceno medio-tardío, a lo largo de las márgenes Pacífica y Caribe de Colombia (pozos: ANH-P27; ANH-P17; ANH-Nueva Esperanza-1X; ANH-San Antero-1; ANH-La Estrella-1X) (ANH-GRP, 2014; ANH, 2017; Grajales et al., 2016), Venezuela y Yucatán (James, 2009). Estos sistemas están reportados igualmente en la Península Azuero y en el área Coiba al sur de Panamá, donde las secuencias abarcan edades desde el Eoceno temprano al Eoceno medio (Buchs et al., 2011). También se reportan al sur de Costa Rica, por las Formaciones Caliza Silícea de Buenavista, Calizas Punta Quepos y Caliza Punta Cuevas (Baumgartner et al., 1984). Adicionalmente, se han reportado al norte de Ecuador, en la Cuenca Esmeraldas por secuencias clásticas carbonáticas de la Formación Azúcar.
Durante el Oligoceno temprano es reportada una discordancia transgresiva en las cuencas intramontanas y de borde continental (Barrero et al., 2007). A partir de los 44 Ma y hasta los 26 Ma, no se han reportado edades de cristalización de cuerpos intrusivos/volcánicos en el occidente colombiano, determinándose así un gap magmático de 18 millones de años. Estos dos procesos podrían estar relacionados a las fases finales del equilibrio dinámico o fases iniciales de la reactivación orogénica.
La edad reportada para el siguiente evento magmático Oligoceno - Mioceno temprano es de 23 Ma y fue establecida en el Stock de Piedrancha en el sector sur de la CO, a partir de los métodos K/Ar en biotita (Álvarez et al., 1981).
CONCLUSIONES
El emplazamiento de plutones del evento magmático Paleoceno - Eoceno en la margen Pacífica colombiana es favorecido por un régimen tectónico transtensional, el cual se desarrolló sobre una margen continental fracturada, previa a la evolución del sistema de subducción. La geometría elongada de los cuerpos plutónicos presenta su eje mayor paralelo a la dirección del esfuerzo principal σ1, por lo que se infiere que las fracturas y zonas de debilidad permitieron el ascenso de cuerpos intrusivos. Las edades de cristalización U/Pb muestran que la provincia magmática se emplazó de 60 Ma a 53 Ma al norte de las cordilleras Central y Occidental. En la Cordillera Occidental se encuentra el Batolito de Mandé y en la Cordillera Central pulsos magmáticos del Batolito Antioqueño. A partir de los 53 Ma tuvo lugar una migración del arco magmático hacia el occidente, hasta los 44 Ma los cuerpos intrusivos más jóvenes son principalmente básicos e intermedios con afinidades toleíticas y se emplazaron rápidamente en un borde continental (corteza transicional) con un espesor < 25-30 km. La composición de los magmas cambió progresivamente a ser calcoalcalinas. La continuidad del cinturón de intrusivos en la Cordillera Occidental es documentada hasta el sur en el límite Colombia-Ecuador.
La disminución significativa en la velocidad de convergencia entre los 40±5 Ma y 26 Ma dio lugar a un proceso denominado Equilibrio Dinámico sobre la margen pacífica colombiana. El equilibrio dinámico es reflejado en distintos ambientes sedimentarios continentales: en ambientes de plataforma continental, por una extensa sedimentación de secuencias clásticas carbonáticas de edad Eoceno Medio-Tardío a lo largo de la margen Pacífica y Caribe de Colombia, Ecuador, Panamá y Costa Rica; en ambientes intra-continentales por la discordancia progresiva regional del Eoceno medio, y la peneplanización de altos topográficos al norte de Colombia, que favorecieron la exhumación de las cúpulas de los plutones, y refleja la reactivación de los piedemontes cordilleranos con progradaciones de abanicos.
MATERIAL SUPLEMENTARIO
Las Tablas S1, S2 y S3 se pueden descargar en el portarl web de la Revista <http://rmcg.unam.mx/>, en la versión html del resumen de de este artículo.










nueva página del texto (beta)


