DOI: https://doi.org/10.18273/revbol.v41n2-2019003

 

 Artículos científicos

 

Aportes al conocimiento del Magmatismo de la Cordillera Central de Colombia en su Flanco Oriental: Área geotérmica de San Diego, Samaná, Caldas

 

Contributions to the Magmatism knowledge of the Central Cordillera of Colombia in its Eastern Flank: Geothermal Area of San Diego, Samaná, Caldas

 

Jesús Bernardo Rueda-Gutiérrez1, * - ORCID: 0000-0002-5575-4702

 

1Grupo de Exploración e Investigación en Recursos Geotérmicos, Servicio Geológico Colombiano, Bogotá, Colombia. (*) jbrueda@sgc.gov.co; jbrueda86@gmail.com

 

 

Forma de citar: Rueda-Gutiérrez, J.B. (2019). Aportes al conocimiento del Magmatismo de la Cordillera Central de Colombia en su Flanco Oriental: Área geotérmica de San Diego, Samaná, Caldas. Boletín de Geología, 41(2), 45-70. DOI: 10.18273/ revbol.v41n2-2019003

 

Creative Commons: https://creativecommons.org/licenses/by/4.0/

 

Trabajo recibido: mayo 16 de 2018

Trabajo aceptado: abril 30 de 2019

 

RESUMEN

Fueron caracterizados en el área geotérmica de San Diego, cinco cuerpos plutónicos, cuyas características químicas posicionan a cuatro unidades (Stock de Florencia, Batolito de Sonsón, Complejo de Samaná y Alaskita de Samaná) en el campo de los granitos tipo I, y a la unidad Intrusivo Gnéisico como granito tipo S. Todos presentan firmas químicas de Rb, Y y Nb típicas de granitoides de subducción, de serie calcoalcalina y tipo magnésicos. Dataciones U-Pb mediante LA-ICP-MS en circones arrojaron edades de 132,8±2,6 Ma para la Alaskita de Samaná, 129,8±9,7 Ma para el Complejo Ígneo de Samaná, 54,6±4,4 Ma para el Stock de Florencia y 47±11 Ma para el Batolito de Sonsón.

Nueve cuerpos porfiríticos andesiticos (Pl+Anf) de los cuales cinco fueron reconocidos en las labores de campo de este trabajo, se expresan algunos de ellos con evidentes formas dómicas y matrices vítreas. Estos cuerpos porfiríticos presentan mayor grado de alteración hidrotermal (fílica y argílica) en especial en el sector de Puente Linda, controlados estructuralmente, con presencia de brechas orientadas este-oeste a lo largo del rio Samaná. Todos estos cuerpos con firmas adakiticas (Sr entre 335,31 y 1160,25 ppm) presentan anomalía positiva de EuN/√(SmN*GdN) desde 1,34 a 2,07. El pórfido de Puente Linda (Pl + Qz + Anf ± Cb ± K ± Dck) tiene sobreimpuesta la mayor alteración en extensión y grado; datación U-Pb arrojó una edad de 2,1±1,9 Ma. Por otro lado, cinco cuerpos aflorantes de naturaleza dacítica (Qz + Pl + Anf) fueron reconocidos y caracterizados, cuya geoquímica contrastante con las demás intrusiones, presentando anomalías negativas de EuN/√(SmN*GdN) desde 0,18 a 0,64. Finalmente, dos estructuras volcánicas, domo de San Diego (Qz + Pl + Bt) y Domos del volcán El Escondido (Pl + Anf) se reconocen como evidencias de actividad volcánicas reciente; dataciones Ar-Ar registraron una edad Plateau de 89±4,4 Ka para San Diego y de 153,7±38,2 Ka para El Escondido.

Palabras clave: San Diego; geotermia; dataciones U-Pb; dataciones Ar-Ar; geoquímica.

 

ABSTRACT

Five plutonic bodies were characterized in the San Diego geothermal area, whose chemical characteristics position four units (Florencia Stock, Sonsón Batholith, Samaná Complex and Samaná Alaskite) in the field of granites type I, and the Gneissic Intrusive as granite type S. All have chemical signatures of Rb, Y and Nb typical of subduction environment granitoids, calcoalcaline series and magnesium type. The U-Pb data by LA-ICP-MS in zircons showed ages of 132.8±2.6 my for the Samaná Alaskite, 129.8±9.7 my for the Samaná Igneous Complex, 54.6±4.4 my for the Florencia Stock and 47±11 my for the Sonsón Batholith.

Nine andesitic porphyritic bodies (Pl+Anf) of which five were recognized in this work, are expressed some of them with evident domic forms and vitreous matrices. These porphyritic bodies present a greater degree of hydrothermal alteration (phyllic and argillic), especially in the Puente Linda sector, structurally controlled, with the presence of oriented breccias (EW) along the Samaná River. All of these bodies with adakitic signatures (Sr between 335.31 and 1160.25 ppm) have positive anomaly EuN/√(SmN*GdN) from 1.34 to 2.07. The Puente Linda porphyry (Pl + Qz + Anf ± Cb ± K ± Dck) has superimposed the greater alteration in extension and degree; U-Pb dating in zircons, showed an age of 2.1±1.9 Ma. On the other hand, five dacite bodies (Qz + Pl + Anf) were equally recognized and characterized, whose geochemistry contrasted with that of other intrusions, with negative anomalies EuN/√(SmN*GdN) from 0.18 to 0.64. Finally, two volcanic structures, San Diego dome (Qz + Pl + Bt) and El Escondido volcano domes (Pl + Anf) are recognized as evidence of recent volcanic activity; Ar-Ar dating recorded a Plateau age of 89±4.4 Ka for San Diego and 153.7±38.2 Ka for El Escondido.

Keywords: San Diego; geothermal; U-Pb dating; Ar-Ar dating; geochemistry.

 

 

INTRODUCCIÓN

 

 

En el marco de las labores de campo en el sector de San Diego, Municipio de Samaná, adelantados por el Servicio Geológico Colombiano (SGC), relacionadas a la exploración e investigación de los recursos geotérmicos en el territorio nacional, fueron caracterizadas las unidades aflorantes de este sector, permitiendo adelantar la caracterización geoquímica, petrográfica (mineralogía y caracterización de alteraciones) y geocronológica (dataciones U-Pb y Ar-Ar) de los cuerpos intrusivos plutónicos, subvolcánicos y volcánicos, dentro del área geotérmica de San Diego.

Trabajos anteriores como el de (Feininger et al. 1972) son la base de la cartografía del área y del reconociendo de basamento metamórfico (Complejo Cajamarca), los cuerpos intrusivos y la única unidad sedimentaria que aflora en la zona de estudio (FIGURA 1). Otros trabajos tuvieron como objetivo el reconocimiento y la caracterización de los volcanes de la zona y sus depósitos asociados. Trabajos como el de (Toro 1989), (Toro 1991) hace la primera caracterización del Volcán de (San Diego y Florez 1987) quien caracteriza uno de los depósitos de esta zona. Más adelante (Monsalve y Ortíz 2013), (Monsalve et al. 2014), (Monsalve et al.2017) y (Monsalve 2015), caracterizan las unidades eruptivas del mismo volcán de San Diego y del recién descubierto volcán El Escondido.

Los primeros indicios de actividad geotérmica en San Diego son reflejados por las manifestaciones hidrotermales registradas en superficie y que fueron en principio referenciadas en el desarrollo del Proyecto de Investigación Geotérmica en la región del Macizo Volcánico del Ruiz (CHEC, 1968). Nuevamente (CHEC, 1981) hace referencia a una anomalía térmica relacionada con actividad volcánica con productos que llaman pórfidos en el área de San Diego.

La actividad termal de la zona ha sido en principio referenciada por parte de CORPOCALDAS, pero ya hoy en día, se tiene reporte de 15 manifestaciones termales, nueve en cercanías a focos volcánicos reconocidos y seis en el costado occidental de la zona de estudio, sector Puente Linda (SGC, 2013).

La presencia de volcanismo reciente (fuente de calor), cuerpos porfiríticos que por sus características sugieren un emplazamiento muy cercano a la superficie con edades pleistocenas, además de la alteración hidrotermal al occidente del área y cuerpos controlados por fallas y fracturamiento (debilidad del subsuelo), hacen necesaria la caracterización a detalle de estas unidades que junto a la evidencia que estas características proporcionan y el análisis de las fuentes termales, posicionan la zona de trabajo como de interés para la exploración de recursos geotérmicos.

 

 

LOCALIZACIÓN Y CONTEXTO GEOLÓGICO

 

 

El área geotérmica de San Diego abarca un área de 460 km2, localizada en el municipio de Samaná, Caldas, compartiendo el extremo occidental con el Corregimiento de Puente Linda, Departamento de Antioquia y al este con el Municipio de Norcasia; geográficamente está localizado en el flanco oriental de la Cordillera Central de Colombia entre las latitudes 5°44’ y 5°29’ N y longitudes 75°09’ y 74°53’, y hace parte del extremo norte de la franja de volcanes de esta cordillera. El área cubre las planchas cartográficas en escalas 1:25.000 168IIIB (IGAC, 1979), 188IA, 188IB, 188IC, 188ID (IGAC, 1961) y 188IIIA, 188-IIIB (FIGURA 1).

La geología está dominada por el Complejo Cajamarca con predominio de esquistos alumínicos, cuarcitas, esquistos verdes y en menor proporción anfibolitas (Feininger et al., 1972). La única unidad sedimentaria corresponde al horizonte fosilífero de Berlín de edad cretácea (Albiano) Etayo, en (Rueda y Rodríguez, 2016) compuesta por shales, areniscas, conglomerados y limolitas calcáreas. Esta unidad ha sido objeto de exploración de uranio por parte de la compañía U3O8 Corp. y casos de estudios más detallados en (Pimiento 2011) quien caracteriza la mineralogía del uranio y algunos elementos menores económicamente explotables, (Ríos-Guerrero 2012) quien analiza la estratigrafía y (Caceres 2012) quien describe a detalle el origen de las concentraciones anómalas de uranio, así como vanadio, molibdeno, fosforo y zinc.

Cuatro grandes cuerpos ígneos intruyen el basamento metamórfico del Complejo Cajamarca: el Stock de Florencia (Tcdf) de edad paleógena, el Batolito de Sonsón (KTcds) de edad Cretácico superior descrito en (González 2002), la Alaskita de Samaná (Kas) y el Complejo Ígneo de Samaná (Kds) los dos de edad Cretácea. Trabajos como el de (Loaiza 2012) caracterizan la totalidad de esta última unidad. Además, aflora un cuerpo intrusivo foliado denominado “Intrusivo Gnéisico de Naranjales” (TRin) caracterizados como cuerpos intrusivos néisicos sintectónicos (Feininger et al., 1972); (González, 1990); (Vinasco et al., 2006); (López et al., 2007), de edad triásica (Vinasco et al., 2006), pero cuyo conocimiento geocronológico es limitado.

La actividad magmática pleistocena está representada por pequeños y numerosas intrusiones hipoabisales de naturaleza andesítica-dacítica y dacítica (Feininger et al., 1972). El vulcanismo reciente está asociado a los volcanes de San Diego (Florez, 1987); (Toro, 1989), (Toro, 1991); (Monsalve y Ortíz, 2013) y El Escondido (Monsalve y Ortíz, 2013); (Monsalve et al., 2017), constituyéndose en el volcanismo más reciente del extremo norte de los Andes Colombianos.

 

Estructuralmente la zona está dominada por la falla Palestina (Feininger, 1970), que atraviesa el centro de la zona en dirección NNE-SSW, junto a fallas paralelas como la falla Patio Bonito y la falla San Diego (Page, 1986). La falla Rio Dulce propuesta por (González, 1980) para la plancha 187, está dada por el trazo del rio Samaná a la altura del sector de Puente Linda.

 

 

FIGURA 1. Ubicación, dataciones y geología del área de estudio. Fallas mostradas solo dentro del área de estudio (Feininger, 1970); (Page, 1986). Geología editada de (Barrero y Vesga 1976), (González et al. 1980a), (González et al. 1980b), (Feininger et al. 1972) y (Rueda y Rodríguez 2016). La línea amarilla representa un límite inferido separando el Batolito de Sonsón en eventos Permo-Triásico al sur y Paleógeno al norte en (Leal-Mejía 2011). La línea azul casi vertical representa el límite  nferido en el mapa geológico escala 1:1.000.000 de (Gómez et al. 2015b) donde al este se  epresenta las edades más antiguas y al  este y arriba de la línea azul se representan las edades Paleógenas. Dataciones fuera del área de estudio tomadas de (Gómez et al. 2015a).

 

 

METODOLOGÍA

Los recorridos para el mapeo del área se realizaron partiendo de la información cartográfica de las memorias explicativas, realizadas por (Feininger et al. 1972) para la plancha 188, La Dorada y 168, Argelia y, de la información de los mapas geológicos de las mismas planchas (Barrero y Vesga 1976); (Feininger et al. 1970). Para esto, tanto las estaciones como las muestras colectadas en ese trabajo fueron georreferenciadas planeando así nuevos recorridos y mayor detalle en la cartografía. Aunque cuerpos nuevos fueron reconocidos, los contactos de las unidades ya descritas en su momento no sufrieron cambios considerables, especialmente para el complejo de rocas metamórficas.

 

 

Análisis químicos (ICP-MS, FRX y DRX) y dataciones U-Pb (LA-ICP-MS) y Ar-Ar

 

 

Un total de 56 muestras de roca de origen intrusivo y volcánico fueron enviadas a los laboratorios del SGC. Cada muestra fue en principio molida mediante una trituradora de mandíbula y posteriormente pulverizada con un molino de tungsteno hasta llegar a un tamaño tamiz malla 200. Un peso aproximado de 100 gr del producto de esta preparación fue rotulado y entregados a los laboratorios químicos del SGC para su posterior análisis.

La caracterización de los elementos mayores se hizo mediante fluorescencia de rayos X (FRX) a 38 muestras con un equipo PANalytical AXIOS Minerals, para esto fueron pesados 10g de muestra a la que posteriormente se le agregó un aglomerante, se homogenizó y se llevó al equipo prensador para formar la pastilla que posteriormente es analizada. Los resultados reportados se expresan en %wt para cada uno de los elementos mayores (SiO2, Al2O3, Na2O, K2O, P2O5, CaO, MgO, TiO2, Fe2O3, y FeO) incluyendo el LOI, exceptuando el manganeso que es presentado en ppm.

Para conocer la química de los elementos, menores y traza se empleó la técnica ICP-MS (Inductively Coupled Plasma-MassSpectrometry) la cual se hizo con un Espectrómetro de Masas con Plasma de Acoplamiento Inductivo con celda de colisión KED, Marca Perkin Elmer REF. NexION 300D. Para el análisis fueron separados 0,1g de muestra haciendo la digestión con ácido clorhídrico, nítrico y perclórico para posteriormente llevar a volumen de 10 ml con solución de ácido nítrico diluido al 3% v/v. Los resultados son presentados en ppm.

Para la caracterización de las fases cristalinas de 35 muestras se utilizó la técnica de difracción de rayos X (DRX) con un equipo X-PERT PRO sin monocromador con fuente de rayos X de Cu. Los resultados recibidos son cualitativos y semicuantitativos y no reportan minerales amorfos, por lo que el 100% de la muestra es presentada como cristalina. Por otro lado, la técnica se hizo en polvo desorientado, haciendo difícil la identificación de minerales arcillosos.

Los circones de las seis muestras que se dataron mediante U-Pb fueron separados teniendo en cuenta el procedimiento de separación del laboratorio de huellas de fisión del SGC (Peña-Ureña et al., 2018). El proceso de datación fue hecho en las instalaciones del SGC sede CAN a cargo del Grupo de Tecnologías Nucleares mediante un espectrómetro de masas (ICP-MS) Thermo Element 2 monocolector acoplado a un sistema Laser Photon Machine de 193 nm. Mediante imágenes de catodoluminiscencia, se ubicaron los puntos para su posterior ablación de acuerdo a las estructuras de cada circón. Las condiciones analíticas fueron: una energía al 55%, un diámetro de ablación de 20 micras, tiempo de ablación de 1:05 minutos. El software para el análisis fue Isoplot 3.5.

Tres muestras fueron datadas mediante Ar-Ar en los laboratorios de la universidad del Estado de Oregon, USA. La técnica se hizo con un espectrómetro de masas ARGUS-VI usando el método de calentamiento gradual. Para el caso de las muestras de este trabajo, donde las rocas tenían una edad tentativa menor a 2 Ma, se hizo una irradiación de 0,5 horas en el reactor TRIGA CLICIT de la Universidad del Estado de Oregon.

RESULTADOS

 

Geología local

 

Las rocas más antiguas que afloran en el área corresponden al basamento metamórfico de esa parte de la Cordillera Central agrupadas en el denominado Complejo Cajamarca (Feininger et al., 1972). Estas rocas se constituyen, en la zona de estudio, de esquistos alumínicos, esquistos verdes dispuestos como lentes de metros hasta kilómetros a lo largo de la foliación dominante (N-S) y esquistos grafitosos localizados en el costado occidental, filitas, cuarcitas aflorantes en el costado central y sur del área y anfibolitas distribuidas al norte del área.

La unidad TRin corresponde a un cuerpo de composición granítica con posterior metamorfismo conocido como Intrusivo Gnéisico por (Feininger et al., 1972) o como intrusivos miloníticos, gneises sintectónicos entre otros términos. Esta unidad corta el costado oriental del área de estudio en dirección casi N-S; en numerosos afloramientos luce como una roca intrusiva de textura masiva de tamaño de grano medio a fino (FIGURA 2A y 2C), aunque al norte del área, en cercanías al rio Samaná la textura metamórfica se hace más evidente con una clara foliación y porfidoblastos de plagioclasas (FIGURA 2B), a simple vista luce cuarzo, plagioclasa, biotitas y granate.

Cuatro cuerpos plutónicos hacen parte de la geología del área de estudio, entre los que están el Complejo Ígneo de Samaná (Kds) y la Alaskita de Samaná (Kas) que aunque son mapeadas por separado, es común encontrar afloramientos leucocráticos dentro del Complejo Ígneo de Samaná; el complejo ígneo presenta diferentes litologías, predominan granodioritas, monzogranitos y en menor proporción cuerpos básicos a manera de pequeños diques (FIGURA 2C, 2D, 2E y 2F). Esta litología es vista solo en el área de trabajo, pero registros como en (Loaiza 2012) referencia que al menos el 60% del cuerpo es de litología diorítica-gabroica.

La unida Alaskita (Kas) se muestra siempre muy meteorizada, en muestra de mano son visibles cristales de feldespatos, cuarzos y algunas biotitas oxidadas (FIGURA 2G). Está en contacto fallado con la secuencia metamórfica al occidente y la unidad sedimentaria al oriente. Afloramientos muestran el carácter intrusivo de cuerpos ácidos en la unida sedimentaria Kib, pero estos se tratarían de las intrusiones dacíticas (Td) y no del leucogranito, igualmente ningún evento térmico se observa en la unidad sedimentaria que suprayace la intrusión.

Afloran también el Stock de Florencia (Tcdf) y parte del Batolito de Sonsón (KTcds) los cuales presentan características muy similares, aunque el primero presenta un tamaño de grano más grande (FIGURA 2I y 2H). El Batolito (KTcds) está en contacto con el pórfido de Puente Linda (Tadp) y muestra una extensa e intensa alteración que dificulta su diferenciación. Las dos unidades atraviesan el basamento metamórfico, tanto a cuarcitas como esquistos alumínicos, pero solo se evidenció un ligero metamorfismo de contacto entre el Stock y la unidad metamórfica de cuarcitas; esta misma unidad Tcdf, sirve igualmente de basamento a las intrusiones dómicas del complejo de domos del volcán El Escondido, y claramente al propio cráter del volcán, evidenciado en numerosos fragmentos angulares de granodioritas dentro de la mayoría de los depósitos volcánicos.

 

 

 

FIGURA 2. Muestras de mano de las unidades intrusivas aflorantes. A. Intrusivo Gnéisico (TRin) con textura ígnea. B. Intrusivo Gnéisico de Samaná con textura gnéisica y augen de plagioclasa. C. Diques máficos en la unidad TRin. D. E. y F. Complejo Ígneo de Samaná (Kds). G. Alaskita o leucogranito de Samaná (Kas). H. Batolito de Sonsón (KTcds). I. Granodiorita del Stock de Florencia (Tcdf).

Al menos 4 cuerpos porfiríticos dacíticos (Td) se encuentran en el área de estudio, dos de ellos reconocidos en (Feininger et al. 1972). Se presentan comúnmente alterados (matriz y fenocristales poco reconocibles), color blanco con fenocristales de cuarzo y asociados geográficamente a la falla Palestina en el contacto entre el complejo metamórfico y la unidad sedimentaria “Horizonte Fosilífero de Berlín (Kib)”.

Un total nueve cuerpos porfiríticos de tamaños no mayores a 3 km2 de los cuales cuatro fueron reconocidos en (Feininger et al. 1972) y llamados pórfidos andesíticos-hornblendíticos (Tad) se encuentran en el área de estudio. En muestra de mano presentan grandes fenocristales blancos de plagioclasa euhedrales de hasta 2 cm y algunos más pequeños de hornblenda, la matriz es de color gris y varia la relación matriz - fenocristales. La mayor diferencia entre ellos radica en el tamaño de los fenocristales de plagioclasa (de 2 mm a casi 1 cm) y de los anfíboles (menos de 1 mm hasta casi 6 mm) (FIGURA 3B, 3C y 3E). Están comúnmente alterados pasando a ser de color amarillento y blanco, y en ocasiones con una matriz endurecida (FIGURA 3E y 3F) o muy fácil de disgregar (FIGURA 3H).

Dentro de este grupo Tad es catalogado según (Feininger et al. 1972) el cuerpo dómico del Maar de San Diego el cual no es claramente hornblendítico sino más bien biotítico con un carácter más dacítico (FIGURA 3J y 3K). Igualmente el complejo de domos del volcán El Escondido de Florencia presenta un aspecto muy similar a los demás cuerpo porfiríticos pero estos son claramente más recientes y correspondería a la actividad intracaldérica del volcán (Monsalve, 2015); (Monsalve et al., 2017) (FIGURA 3L y 3M).

Excluyendo los cuerpos dómicos relacionados a las dos estructuras volcánicas, algunas de las demás litologías porfiríticas del grupo Tad presentan relieves elevados con geoformas dómicas (FIGURA 3A, 3G y 3I) sugiriendo una intrusión que llegó a la superficie, lo cual se ve reflejado también en la petrografía, donde se observan las matrices muy vítreas.

En cercanías al volcán El Escondido en dirección sur, una estructura cilíndrica resalta de la topografía cubierta por vegetación alta del Parque Natural Selva de Florencia (PNNSF). En la FIGURA 3N se observa una pared inclinada casi 90°, justo en ese punto la roca que se observó es una brecha con líticos porfiríticos y algunos otros metamórficos (FIGURA 3O). Antes de llegar a ese punto, solo se controlan rodados de roca porfirítica (FIGURA 3P) y casi ningún afloramiento es visible, pero es claro que se trata de un cuerpo de naturaleza intrusiva. Por su forma, textura y debido la actividad volcánica de la zona, pero sin tener un estricto control de campo más al sur de dicha geoforma y dados los escasos afloramientos, se plantea que se trate de una diatrema.

 

Petrografía

 

La unidad Intrusivo Gnéisico muestra variación composicional entre granitos (monzogranitos), granodioritas y tonalitas, pueden tener estructura foliada con porfidoblastos de plagioclasa alterada (FIGURA 4A). Predomina el cuarzo con 24 a 45%, seguido de plagioclasa entre 16 y 47%, y cantidades variables de feldespato potásico entre 5 y 38%, biotita 2 y 15% con trazas de granate, piroxenos y minerales de alteración como sericita y clorita (FIGURA 4B).

La unidad Complejo de Samaná claramente muestra variación composicional y mineralógica con litologías félsicas de granodioritas-monzogranitos con cuarzo entre 35 y 50%, plagioclasa entre 17 y 35% (generalmente tipo albita) y feldespato potásico entre 15 y 17% (tipo sanidina). Las litologías gabrodioritas presentan cuarzo en un 2,3%, plagioclasa entre 48 y 59% y anfíbol tipo hornblenda entre 22 y 36% con cantidades variables de sericita <20% y minerales opacos (FIGURA 4C y 4D). En la unidad Alaskita de Samaná predominan los granitos leucocráticos con feldespato potásico con 34,3%, plagioclasa 26% y cuarzo 24,3%, comúnmente alteradas a sericita y cantidades variables de biotita (FIGURA 4E).

 

En el Stock de Florencia predominan las granodioritas, cuarzodioritas y tonalitas de grano medio a grueso con plagioclasa tipo albita principalmente entre 45 y 62%, cuarzo entre 16 y 37%, feldespato potásico entre 1,7 y 13%, biotita entre 4 y 13% y trazas de anfíboles tipo hornblenda, sericita y clorita (FIGURA 4F). Para el Batolito de Sonsón en la zona de estudio, se observó predominancia de granodioritas de grano fino con plagioclasa en un 35,6%, cuarzo 34%, feldespato potásico con 16% y biotita 9,6% con variación en el contenido de anfíbol, clorita y sericita (FIGURA 4G).

Los cuerpos subvolcánicos son de composición andesítica y dacítica. Las andesitas tienen fenocristales (>2 mm) y microfenocristales (2-0,5 mm) de plagioclasas (predominio de albita y labradorita mediante análisis por DRX) microfenocristales de anfíboles hexagonales y cuarzo bipiramidal; la matriz varía de microcristalina a criptocristalina, con microcristales de anfíbol, con porcentajes entre 50 y 90%, en algunas rocas la matriz vítrea predomina sobre la cristalina (FIGURA 4J y 4K). Las dacitas tienen microfenocristales de cuarzo bipiramidal (FIGURA 4H) entre 6 y 10%, en menor proporción de plagioclasa entre 4 y 7% y hornblenda puede estar presente. La matriz de estos cuerpos varía de 51% a 80% y comúnmente se muestra alterada.

El pórfido de Puente Linda tiene matriz microlítica con fenocristales y microfenocristales de plagioclasa entre 5,8-13,8% y presenta relleno de carbonatos en fracturas, cuarzo puede estar presente en (0,5-5,6%), hornblenda entre 0,31-3,2% parcialmente absorbida por la matriz o alterada a carbonatos (FIGURA 5C).

 

 

FIGURA 3. Aspecto de las unidades subvolcánicas (Tad) y volcánicas. A. y B. Cuerpo del sector noroccidental del área de estudio. C. D. E. F. Muestra de mano y afloramientos de los cuerpos porfiríticos alterados en el sector de Puente Linda a lo largo del Rio Samaná. G. H. I. Afloramiento y muestra de mano de los cuerpos porfiríticos en sector de la Quiebra. J. K. Panorámica del domo de San Diego y detalle de la roca. L. M. Panorámica del complejo de domo del volcán Escondido (Florencia) y detalle de la roca. N. y O. Panorámica de la diatrema y detalle de muestra de mano. P. Detalle de la roca alrededor de la diatrema.

 

 

 

El domo de San Diego es dacítico con presencia de fenocristales de plagioclasa y en menor proporción microfenocristales y microcristales de cuarzo anhedral entre 6-8%, biotita tabular rojiza y parda entre 15 y 25% y plagioclasa subhedral que varía entre 6 y 10% (FIGURA 4N) dentro de una matriz criptocristalina félsica. Según DRX las fases de feldespato reportadas para este cuerpo son sanidina y anortita.

El complejo de domos del volcán El Escondido es andesítico con fenocristales y microfenocristales de hornblenda euhedrales hexagonales (FIGURA 4M) (magnesohormblenda identificada por DRX) y plagioclasa subhedral según DRX tipo albita y en menor proporción andesina, con microcristales de piroxenos siendo rodeados por microcristales de hornblendas (FIGURA 4L). La matriz es microcristalina y vítrea para algunos de los líticos dentro de los depósitos de origen volcánico.

La geoforma cilíndrica al sur (¿Diatrema?) microscópicamente muestra líticos entre 0,1 y mayores a 1,5 mm de roca porfirítica andesítica con cristales angulares de plagioclasa y anfíboles, el otro grupo de líticos son metamórficos con composiciones cuarzosas y algunos otros esquistosos con evidente foliación superando algunos los 3 cm (FIGURA 4O). También son visibles cristales libres comúnmente de plagioclasas con formas bien definida y algunos anfíboles, todo esto dentro de una matriz criptocristalina de color marrón que no supera el 40% de la roca.

 

 

FIGURA 4. Microfotografías cuerpos plutónicos, hipoabisales y volcánicos. A. y B. Unidad TRin con clorita (Chl), granate (Gr) y sericita (Ser). C. Unidad Complejo de Samaná, facie granítica. D. Unidad Complejo de Samaná, facie diorítica-gabroica. E. Alaskita de Samaná. F. Unidad Stock de Florencia. G. Batolito de Sonsón. H. Pórfidos dacíticos con cuarzo bipiramidal. I. Cuerpos andesiticos. J. K. Cuerpo andesítico con matriz vítrea. L. M. Domo de El Escondido con plagioclasa (Pl) y Anfíbol (Anf) tipo hornblenda principalmente y piroxeno (Px). N. Domo volcán San Diego con Qz, Pl y Bt en matriz microcristalina de Qz. O. Brecha al sur del área de estudio (LtV: Lítico volcánico. LtM: Lítico metamórfico).

 

Alteraciones

 

Gran parte del área de estudio esta demarcada por una intensa meteorización de las unidades aflorantes, predominando la oxidación de la unidad ígnea (TRin) y esquistos sericíticos (TRes) y unidades saprolizadas en el caso del Stock de Florencia (Tcdf) y Batolito de Sonsón (KTcds). Los esquistos en todas sus variedades presentan pequeños sectores con alteraciones que hacen cambiar el color de la roca, y cuyas extensiones varían desde pocos metros a cientos. Los cuerpos de composición acida presentan una coloración blanca de textura arenosa y los intermedios-ácidos (Tad) se muestra mayormente alterados con textura arcillosa de color naranja. En el sector norte se muestra una meteorización fuerte en la unidad gnéisica en la mayoría de su extensión y alteraciones con presencia de clorita y sericita. En el sector donde aflora en mayor extensión el Complejo Ígneo de Samaná (Kds), el leucogranito (Kas) y esquistos sericíticos con lentes de esquistos verdes (TRev) presenta una apariencia alterada superficialmente, especialmente para el cuerpo alaskítico que por su composición ácida y por el intemperismo se muestra siempre fácil de disgregar. Las manifestaciones de alteración hipógena más notorias se localizan en el sector de Puente Linda. Esta alteración se ve igualmente en el Batolito de Sonsón en el contacto con el cuerpo de Puente Linda donde en muestra de mano la roca luce muy arcillosa, pero con evidencias del carácter granítico de la roca. Dentro del propio cuerpo porfirítico Tadp, se ven sectores de alteración argílica y fílica con una marcada mineralización de sulfuros por zonas. Este cuerpo presenta sectores controlados estructuralmente con intensa alteración y brechas presentando una alteración sericita+cuarzo+pirita. Los tipos de alteración más común se resumen en la TABLA 1 y se ilustran en la FIGURA 5

 

 

TABLA 1. Resumen de las alteraciones y resultados de difracción de rayos X para algunas de las muestras del área de estudio, especialmente en el sector de Puente Linda (Qz: Cuarzo. Ser: Sericita. Dck: Dickita. Ill: Illita. Cb: Carbonato. K: Caolinita. Ch: Clorita).

 

 

.

Cloritización y Seritización (débil): Reconocida por una evidente cloritización de micas y algunos anfíboles, mientras que la sericita ataca a los feldespatos sódicocálcicos y potásicos. Se presenta en los grandes intrusivos, especialmente en la unidad gnéisica. Los resultados de DRX reporta la facie moscovita que haría referencia tanto a la moscovita primaria que es evidente en sección delgada como a sericita, que es igualmente reconocible bajo el microscopio. Igualmente, la fase clinocloro es reconocida por DRX, el cual se trata de un mineral del grupo de las cloritas. Fílica (Sericítica): Afecta ligeramente a los grandes intrusivos (TRin, KTcds, Kds y Kas) reemplazando la plagioclasa por sericita y micas por cloritas, en ocasiones en sus núcleos o remplazándola totalmente como es el caso del intrusivo gnéisico. Se da mayormente en los cuerpos andesíticos en el sector de Puente Linda y en el propio pórfido de Puente Linda, donde altera los fenocristales y la matriz. Argílica: Representada por caolinita (K) y cuarzo (Qz) en los cuerpos porfiríticos y acompañada de sericita y dickita (Dck) reemplazando mayormente plagioclasa, también se presentan carbonatos (Cb) remplazando mayormente anfíboles y plagioclasa. Se expone en la mayoría de los cuerpos andesíticos y en el de Puente Linda, exceptuando el cuerpo al sur del área de estudio, el cual se muestra mucho menos alterado que los demás.

 

 

FIGURA 5. Microfotografía de algunas alteraciones en los cuerpos porfiríticos. A. y B. Pórfidos dacíticos con cuarzo y algunas plagioclasa y anfíboles, comúnmente alterados a caolinita (K). C. Cuerpos andesiticos (Tadp) con carbonatos (Cb) reemplazando anfíboles (Anf). D. E. F. Cuerpo andesítico con matriz vítrea, cristales oxidados y reemplazamiento de plagioclasa por sericita.

 

 

Geoquímica

 

Los resultados de geoquímica de roca total (TABLA 2 y 3) de los cuerpos plutónicos representados en 25 muestras dejan hacer las siguientes observaciones.

El cuerpo Intrusivo Gnéisico muestra una concentración de SiO2 que va desde 65,21 hasta 76,4 %wt, los valores de Al2O3 van de aproximadamente de 12 a 13,5%, el K2O va desde 2,2 a 2,79%. Las variaciones del CaO son evidentes en estos cuerpos, tanto para lo muestreados en el área de estudio como los reportados en otros lugares, los valores tienen picos de 8 %wt y bajos de menos de 1 %wt. El diagrama TAS clasifica las muestras como granitos y granodioritas (FIGURA 6A).

En la unidad Complejo Ígneo de Samaná los valores de SiO2 van de 68,15 a 74,18 %wt para las muestras más acidas y el Al2O3 va de 14,02 a 18,14. Los contenidos de Na2O y K2O están en el rango de 3,38 y 4,16 y entre 3,4 y 0,75 respectivamente y el Ti2O tiene cantidades de 0,34 a 1,98. Los valores de las muestras básicas muestran una variación importante en SiO2 con 53,13 %wt y en Na2O y K2O, con 4,1 y 0,74 %wt respectivamente. El Ti2O se presenta con 2,94 %wt muy por encima que los valores de la roca caja (rocas acidas) al igual que el Fe2O3 con 13,721 %wt. La clasificación TAS ubica las rocas del complejo ígneo en el campo de los gabros hasta los granitos, pasando por dioritas y cuarzodiorita (FIGURA 6A y 6B).

La unidad Alaskita de Samaná presenta valores de SiO2 entre 75,56 y 79,26 %wt, Al2O3 de 12 a 18 %wt. El Na2O y K2O van de 2,56 a 3,37 y 1,56 a 5,08 respectivamente. Para el CaO se dan valores de 2,5 hasta por debajo del límite de detección, así como el MgO con valores menores a 0,2 %wt.

La unidad Stock de Florencia muestra valores de SiO2 entre 61,83 y 67,31 %wt y 72,42 para una muestra, el Al2O3 está entre el 16 y 20 %wt, para el Na2O y el K2O no sobrepasan el 3,7 y 1,88 %wt respectivamente, con un valor alto para la misma muestra (GR-246-1) de 0,2 de Na2O y 3,2 %wt de K2O. Según clasificación TAS las rocas se plotean en el campo de las granodioritas-diorita.

Las rocas del Batolito de Sonsón presentan contenidos de SiO2 que varían de 68,68 y 74,81 %wt de la roca total, el Al2O3 tiene una variación de más de 5 %wt entre la muestra con mayor y menor contenido de Al2O3. Para el Na2O los valores van de 3,09 a 3,62 y el K2O de 2,4 y 4,3 %wt. El CaO se mantiene casi constante con valores de 1,2 y 1,1 %wt y se clasifican como rocas graníticas y cuarzodioríticas según diagrama TAS.

 

 

 

TABLA 2. Geoquímica de elementos mayores.

 

 

En general todas las muestras, según el contenido de K2O vs SiO2 se ubican en el campo de la serie calcoalcalina (FIGURA 6C) en el diagrama de (Peccerillo y Tylor 1976). La clasificación según Índice de Saturación de Alúmina (ASI) posiciona la mayoría de las muestras en el campo de rocas peraluminosas (FIGURA 6D) lo cual aportaría evidencia en principio de una fuente de fundido de un protolito de naturaleza pelítica o que el magma es muy diferenciado, pero como lo explica (Miller 1985), el carácter peralumínico de este tipo de rocas podría estar ligado a fundido félsico metalumínico. Las muestras con un comportamiento fuertemente peraluminoso explicaría la concentración de alúmina presente en fases como moscovita, cordierita, granate, por otro lado, los comportamientos débilmente peralumínicos estarían relacionados a crecimiento de biotita alumínica. Lo anterior explicaría y reafirmaría el componente pelítico de la unidad Intrusivo Gnéisico (TRin) la cual muestra una mayor dispersión de rocas hacia el campo peraluminoso y cuya mineralogía está representada por moscovita, biotita, feldespato potásico, granate entre otros, y donde además, para el intrusivo gnéisico de Abejorral es reportada la paragénesis de cordierita+moscovita y para el de Rio Verde son reportados porfidoblástos de cordierita y feldespato potásico (Vinasco et al., 2006). Las demás unidades muestran un carácter metalumínico o ligeramente peralumínico. Los valores de ASI menores a 1 reflejarían un exceso de calcio ajustado en fases como anfíboles tipo hornblenda (Frost et al., 2001) vista en cierta abundancia en la mineralogía de la mayoría de las unidades. El diagrama de tipos de granitos de (Chappell y White 1992) de la FIGURA 6E revela igualmente el posible carácter pelítico del fundido (granitoides tipo S) de la unidad gnéisica (TRin) en relación con los demás cuerpos graníticos del área de estudio (Granitoides tipo I) y se confirma con la clasificación de tipo de protolito de (Norman et al. 1992) para algunas de las muestras, mostrado en la FIGURA 6F. Según clasificación de (Frost et al. 2001) los granitos de este trabajo se ubican en el campo de los tipos magnésicos y calco-alcalino y/o cálcicos (FIGURA 6G y 6H) característico de granitos cordilleranos o granitos de arco volcánico para (Pearce et al. 1984). En la FIGURA 6G, se ubican en el área del promedio para rocas Mesozoicas graníticas tipo cordilleranas de América del Norte (Área morada). En esta misma figura son promediados leucogranitos peralumínicos o granitos tipo S producto de fundido metasedimentario. En este caso las muestras de la unidad gnéisica (símbolos verdes) tiene una mayor dispersión cayendo solo unas pocas en el promedio de rocas en (Frost et al. 2001), aunque se presentan tanto como granitos magnésicos y férricos, mostrando como característica especial, el alto contenido de SiO2. Tanto la clasificación de Thiéblemont (1999) con elementos inmóviles como el Zr y Nb (FIGURA 6I) como la clasificación mediante Y, Nb y Rb de (Pearce et al. 1984) (FIGURA 6J y 6K) ubican las muestras en el campo rocas derivadas de magmatismo en zonas de subducción especialmente para las rocas más evolucionadas de carácter acido.

 

 

TABLA 3. Geoquímica de elementos menores. Ver continuación tabla 3

 

 

 ContinuaciónTABLA 3. Geoquímica de elementos menores.

 

 

FIGURA 6. Clasificaciones químicas para los cuerpos plutónicos. A. Clasificación TAS para rocas plutónicas (Cox et al., 1979). B. Clasificación para rocas plutónicas saturadas en sílice (Middlemost, 1994). C. y D. Saturación de alúmina de (Maniar y Piccoli 1989). E. Clasificación de tipo de granitos según valores Na2O y K2O de granito (Chappell y White, 1992). F. Tipo de protolito de (Norman et al. 1992). G. y H. Clasificaciones de (Frost et al. 2001) teniendo presente cantidad de FeO-MgO y Na2O, K2O y CaO respectivamente. I. Discriminación tectónica diagrama logarítmico de (Thiéblemont 1999). J. y K. Discriminación tectónica (Pearce et al., 1984). GAV: Granitos de Arco Volcánico.

 

 

Respecto a los cuerpos porfiríticos, volcánicos y subvolcánicos, se colectaron 11 muestras de los pórfidos andesíticos (Tad), cuatro de los dacíticos (Td), seis del pórfido de Puente Linda (Tadp) y cuatro de los domos volcánicos (TABLAS 2 y 3). Se realizó el filtro de alteración de (Wilt 1995) donde la mayoría de muestras caen en el campo de rocas frescas (FIGURA 7A y 7B), la excepción son las rocas dacíticas (Td) cuya mineralogía es claramente alterada, vista tanto en muestra de mano como bajo el microscopio.

Los cuerpos andesíticos muestran valores de SiO2 entre 65,78 y 75,53 %wt. Al2O3 entre 15,51 y 18,14. El CaO está entre 0,92 y 3,87, Na2O entre 2,24 y 4,77 y K2O entre 4,79 y 1,85. El contenido de SiO2 para las muestras JR-271-1 y JR-286-1 dan por debajo del promedio, con un rango de 61 a 58 %wt respectivamente, igualmente los porcentajes de Al2O3 muestra valores elevados de 20,9 y 23,59 por encima de los otros cuerpos con excepción de la muestra IA-002-1 con 20,41 %wt en Al2O3 y 61,68 en SiO2. Estas mismas muestras registras valores altos en MgO de 2,37 y 4,28 por encima del promedio de 1,28 para los demás cuerpos porfiríticos. La clasificación TAS posiciona estos cuerpos desde las andesitas hasta las dacítas y riolitas, aunque petrográficamente claramente caen en el campo de las andesitas con hornblenda, químicamente el contenido de SiO2 varia ampliamente seguramente por la variación de cuarzo modal, que en algunos casos, es más alto que los valores de plagioclasa (FIGURA 7C).

Para la unidad pórfido de Puente Linda los valores de SiO2 van de 59,07 a 68,50 %wt. Los porcentajes de Al2O3 muestran valores de 15,64 a 19,92 %wt aumentando a medida que aumenta la sílice, el CaO presenta variaciones importantes entre las muestras, por un lado hay valores mayores 5 %wt y otro grupo que va de 1,29 a 2,95 %wt. El Na2O y K2O muestran rangos de entre 2,50 a 3,98 y 2,14 a 3 %wt respectivamente. Las muestras según clasificación TAS caen en el campo de las andesitas y dacitas (FIGURA 7C).

 

 

FIGURA 7. Clasificaciones químicas para los cuerpos porfiríticos. A. y B. Filtros de alteración (Wilt, 1995). C. Clasificación TAS. D. Clasificación de rocas porfiríticas (Winchester y Floyd, 1977). E. Clasificación de rocas subalcalinas (Peccerillo y Tylor, 1976). F. y G. Discriminación tectónica (Pearce et al., 1984). GAV: Granitos de Arco Volcánico. H. Clasificación de cuerpos porfiríticos en diagrama de Defant y Drummond (1990) para adakitas. I. Clasificación de cuerpos porfiríticos en diagrama de (Martin 1999) para rocas adakiticas. En azul y rojo se muestran las rocas del Nevado de Ruiz y del este de Manizales

 

 

Los valores de la química de los cuerpos dacíticos (Td) como se esperaría para cuerpos tan ácidos, muestra concentraciones de SiO2 que sobrepasan el 75 %wt para tres de las muestras (79,82 a 84,28 %wt) y solo una registra un valor menor de 71,32 % wt. El contenido de Al2O3 está entre el 10 y 17 % wt, para Na2O y el K2O los valores no son mayores al 4 y 3,7 % wt respectivamente. El CaO en todas las muestras da por debajo del límite de detección. Según clasificación TAS las muestras de los cuerpos porfiríticos dacíticos por su alto contenido de SiO2 se salen del área del diagrama.

Se utilizó el diagrama de (Winchester y Floyd 1977), que permite hacer la clasificación con elementos inmóviles. La FIGURA 7D muestra que la mayoría de muestras grafican dentro del campo de las andesitas, para los cuerpos dacíticos se ve una dispersión entre riolitas y traqui-andesitas.

Las muestras del complejo de domos del volcán El Escondido caen en el campo de las andesitas, mientras que para el domo de San Diego se ubican en límites de las traqui-andesitas y basaltos alcalinos (FIGURA 7C).

Según el diagrama de (Peccerillo y Tylor 1976) la mayoría de las muestras caen en el campo de rocas de serie subalcalina (FIGURA 7E), calcoalcalina con cuatro muestras en el campo de la serie calcoalcalina alta en potasio (FIGURA 7F). Por último, todas las muestras caen en el campo de los granitoides de arco volcánico según el comportamiento de Rb vs Y+Nb y Nb vs Y (FIGURA 7G y 7H). Para las rocas porfiríticas volcánicas y subvolcánicas se graficó el diagrama de clasificación de rocas adakiticas (FIGURA 7I), donde se evidencia que todos los cuerpos dómicos (domos de los volcanes San Diego y El Escondido), los cuerpos denotados como porfiríticos andesíticos y el pórfido de Puente Linda se ubican en el campo de las adakitas.

El patrón de REE para los cuerpos volcánicos y subvolcánicos muestra lo siguiente. El patrón de REE en la unidad andesitas (Tad) (FIGURA 8A) de manera general presenta enriquecimiento de las LREE respecto a los las HREE, anomalía positiva en Eu con valores de EuN/√(SmN*GdN) desde 1,34 a 2,07. Para la unidad dacitica (Td) (FIGURA 8A) se nota un patrón con una mayor abundancia y dispersión de las HREE y con anomalía negativa en Eu con un valor de EuN/√(SmN*GdN) desde 0,18 a 0,64. Para la unidad pórfido de Puente Linda (FIGURA 8B) la dispersión de las concentraciones es menor que en los cuerpos anteriores mostrando homogeneidad en las LREE y dos patrones en la HREE. Una de estas rocas que muestra un patrón de HREE diferente, hace referencia a una roca colectada del contacto con esquistos verdes y cuya alteración era evidente.

Las cuatro muestras de los cuerpos volcánicos (FIGURA 8C) reflejan el mismo patrón de REE con ligeras variaciones en las cantidades muy seguramente debido al grado de diferenciación. Una de las muestras refleja un comportamiento anómalo especialmente en las HREE aunque muestra la misma tendencia. Todas muestran una anomalía ligeramente positiva de EuN/√(SmN*GdN) 1,45 a 4,76 para San Diego y de 1,65 a 1,11 para El Escondido.

La mayoría de los cuerpos plutónicos muestran anomalías negativas en Eu exceptuando algunas de la unidad TRin. Esta última unidad muestra variación en la HREE respecto a otras muestras de este cuerpo, pero en lugares diferentes, y aunque están mapeados como TRin, es claro que la mineralogía varía entre cuerpos; que este cambio sea de las HREE reflejaría variación en la presencia de granate en la fuente.

A diferencia de lo anterior, al menos dos muestras de la unidad Complejo de Samaná presentan un comportamiento casi planar y leve anomalía en Eu con valores de EuN/√(SmN*GdN) de 1,06 y ligeramente más enriquecidas en HREE, mientras que la misma unidad para las muestras graníticas registran valores de EuN/√(SmN*GdN) de 0,42 a 0,54 (FIGURA 8E).

Para el Stock de Florencia se observa un patrón de enriquecimiento de LREE respecto a las HREE es notoria una leve anomalía tanto positiva como negativa en Eu donde los valores de EuN/√(SmN*GdN) van de 0,43 a 0,8 para las anomalías negativas y de 1,19 y 1,42 para las positivas (FIGURA 8D). Los resultados del Batolito de Sonsón de las muestras analizadas aquí presentan un enriquecimiento de LREE respecto a las HREE con valores de EuN/ (SmN*GdN)1/2 de 1,26 y 0,78 mostrando una leve anomalía positiva y negativa.

 

FIGURA 8. Abundancia de tierras raras (REE) vs Condrita de (Nakamura 1974) para los cuerpos: A. y B. Subvolcánicos. C. Volcánicos. D. E. y F. Plutónicos.

Geocronología

circones para cinco muestras y resultados Ar-Ar para A continuación se presentan los resultados de la geocronología mediante U-Pb por LA-ICP-MS en     dos muestras de los cuerpos dómicos, cuyos resultados se resumen en la TABLA 4 y FIGURAS 9 y 10.

 

 

TABLA 4. Resultados de edades U-Pb para cinco cuerpos plutónicos y Ar-Ar para tres volcánicos  

 

 

El evento más antiguo datado corresponde a la unidad Intrusivo Gnéisico a la cual se le hizo ablación a 63 puntos en un total de 46 cristales de los cuales 15 presentaba una zonación oscilatoria visible en las imágenes de catodoluminiscencia característica de circones ígneos. Al menos 21 ablaciones reportaron valores con porcentajes de discordias por encima del 100% con edades entre 200 mil y 2,5 Ma, las cuales fueron descartadas como la edad del evento. Dos edades recientes de 86,16±5,30 y 78,43±5,51 Ma, cuyas relaciones Th/U presentan valores de 0,52 y 0,57, además su zonado concéntrico, sugerirían un origen ígneo. Los valores en los extremos de los cristales (último zonado o cierre del cristal) pueden ser interpretados como la fecha del evento metamórfico que afecto a este cuerpo, cuya edad dio 243±19 Ma. En (Barrero y Vesga 1976) reportan una edad de 207±7 Ma para este mismo cuerpo muy cerca al punto de muestreo de este trabajo mientras que (Vinasco et al. 2006) realizaron dataciones U-Pb y Ar-Ar para varios de los cuerpos denotados como granitos gnéisicos con edades de 230; 274±9,6; 302±4; 290±10 Ma. La variación de edades puede ser producto de la herencia de pulsos más recientes o la variación del metamorfismo en distintos lugares del mismo cuerpo.

El leucogranito o Alaskita de Samaná registra una edad de 132,8±2,6 Ma (FIGURA 9B) como resultado de 40 ablaciones en igual número de circones (FIGURA 9A) donde al menos cinco de estos cristales presentaban una zonación oscilatoria con relaciones Th/U de entre 0,5 y 0,9 (FIGURA 9C). Para el complejo de Samaná se realizaron 75 ablaciones en un total de 72 circones (FIGURA 9D[UdW1] ) de los cuales a tres cristales se les realizó ablación en el núcleo y las zonaciones que mostraban con relaciones Th/U de entre 0,3 y 0,7 (FIGURA 9F). Al menos tres cristales registraron edades de mayores a 295 Ma consideradas como heredadas. Haciendo los recálculos la edad reportada es de 129,8±9,7 Ma (FIGURA 9E). Se tiene el reporte de edades menores p. ej. 119±10 en facie gabroica (Vesga y Barrero, 1978) mediante K-Ar (edad de enfriamiento), pero esta facie al menos en el área de norte de este trabajo se presenta como diques que cortan la facie tonalítica más general y cuya datación reporta la edad más alta.

Todas estas unidades cortan la secuencia metamórfica (complejo Cajamarca) cuya edad de metamorfismo sería anterior al Triásico (Vinasco et al., 2006), por lo que algunos de los circones heredados en las muestras datadas podrían provenir de dicha unidad y/o de núcleos más antiguos.

El Stock de Florencia se dató en circones con zonaciones sectorizadas y concéntricas y con relaciones Th/U típicas de circones de origen ígneo entre 0,3 y 0,8 (FIGURA 9I). Se reportó una edad de 54,6±4,4 Ma (FIGURA 9H), como resultado de la ablación de 64 puntos en igual número de cristales (FIGURA 9G), así como edades heredadas desde 150 hasta 500 Ma. Para el Batolito de Sonsón se realizaron 45 ablaciones en igual número de circones (FIGURA 9J), de ese total nueve mostraron edades de entre 115 y 290 Ma. La mayoría de los circones mostró una zonación sectorizada y unos pocos una zonación concéntrica o muy débilmente zonada con morfología de circones alargada y cuya relación de Th/U fue bastante dispersa (FIGURA 9L). Se reportó una edad de 47±11 Ma como edad mínima de cristalización (FIGURA 9K).

Para el pórfido de Puente Linda se realizó ablación a un total de 88 puntos (FIGURA 9M) para el mismo número de circones con zonaciones poco evidentes, al menos 6 de estos circones mostraban un núcleo aparente donde se procedió a hacer la ablación. Cuatro edades registradas sobrepasan los 40 Ma, por lo que son consideradas como heredadas. Los valores de la relación Th/U se enmarcan entre 0,1 y 0,5 (FIGURA 9O) con edades entre 1,7 y 2,3 Ma, y se consideran evidencia de datos que reflejan el evento ígneo. La edad reportada es de 2,1±1,9 Ma (FIGURA 9N).

El domo de San Diego, uno de los domos de Florencia y un cuerpo porfirítico en el mismo sector de Florencia pero al oeste del complejo de domos del volcán El Escondido se dataron mediante Ar-Ar (FIGURA 10).

Para el domo de San Diego se registró una edad de 89±4,4 Ka (FIGURA 10A). Teniendo presente lo propuesto en (Monsalve 2015) y (Monsalve et al. 2014), los depósitos piroclásticos de este volcán estarían asociados a una actividad posterior a la intrusión del cuerpo dómico y por lo cual deberían presentar edades mucho más jóvenes que esta estructura. Por otro lado, es posible pensar en que el domo tuviera su formación antes de la actividad freática del maar y que los depósitos que mantelan desde la superficie al domo, no provengan de este, sino más bien fueron dirigidos desde el cráter hacia el NE, mantelando así la estructura dómica. Trabajos como el de (Woodwar-Clyde Consultants 1980) registran este domo con una edad menos a 50.000 años y (Toro 1991) reporta una edad de 18.700 años con un error considerable de ±10.000.

 

 

FIGURA 9. Histograma de edades, diagrama de concordia y relación Th/U (sin edades heredadas) respectivamente de izquierda a derecha. A, B y C. Unidad Alaskita de Samaná. D, E y F. Unidad Complejo Ígneo de Samaná. G, H y I. Unidad Stock de Florencia. J, K y L. Unidad Batolito de Sonsón. M, N y O. Unidad Pórfido de Puente Linda.

 

 

FIGURA 10. Diagrama de edad Plateau para los cuerpos dómicos. A. Domo San Diego. B. Domo del Florencia (Volcán El Escondido). C. Estructura al oeste del Domo de Florencia.

 

 

La datación de uno de los domos de Florencia (Volcán El Escondido) con una edad de 153,7±38,5 Ka (FIGURA 10B) representa un problema puesto que dicho complejo está asociado a la última actividad intracaldérica del volcán El Escondido cuyos depósitos volcánicos están datados en alrededor de los 33 Ka (Monsalve et al., 2014; Monsalve 2015). La datación del domo supone que la actividad efusiva fue mucho antes que la actividad que dio lugar a los depósitos de los alrededores del volcán. Por otro lado, la estructura al W del domo refleja una edad de 55±23,6 Ka (FIGURA 10C), cuyo error daría un margen a favor para correlacionar dicha estructura, con la actividad explosiva del volcán.

 

 

DISCUSIÓN

 

 

Las intrusiones plutónicas, sobvolcánicas y volcánicas de firma calcoalcalina y afinidad adakitica para los cuerpos recientes (unidades porfiríticas menores a 2 Ma) registran una intensa evolución magmática cuyos pulsos se emplazaron en varios niveles de la corteza. La unidad Intrusivo Gnéisico refleja un evento sobreimpuesto de deformación (sintectónica), planteando la discusión de la forma como se le debe llamar a este cuerpo. Las intrusiones cretácicas (complejo de Samaná y leucogranito o Alaskita de Samaná) muestran edades muy similares y se pueden interpretar como una evolución de un mismo pulso dando origen a la variedad de litologías dentro del complejo, como a la fase más evolucionada representada por granitos leucocraticos. Los cuerpos paleógenos (Batolito de Sonsón y Stock de Florencia) muestran similitudes en petrografía (relación porcentual de cuarzo, plagioclasa y feldespatos, microclina y albita para el Stock y para el Batolito de Sonsón ortoclasa, albita, labradorita) y geoquímica (REE y elementos traza, con anomalías negativas en Nb y Ti) sugiriendo un mismo pulso con diferencias geocronológicas muy mínimas en su cristalización.

De gran interés resulta la edad registrada para el Batolito de Sonsón, pues los antecedentes reportan dos edades tanto Paleógenos como Triásicas. Esta edad abre la discusión respecto a las otras dataciones registradas en otras partes del Batolito y por diferentes autores. Por un lado (González 1980), (González 2002) mediante K-Ar en roca total registró valores de 170±10 Ma y 160±4 Ma. (Pérez 1967) mediante K/Ar reporta una edad de 69,3±3 Ma, (Ordoñez et al. 2001) reporta edad de cristalización de 60,7±1,4 Ma mediante SHRIMP U-Pb en zircones y (Leal-Mejía 2011) mediante el método U-Pb en circones con LA-MC-ICP-MS reporta una edad de 58,6±0,9 Ma para una granodiorita, facie principal del Batolito y 55,8±0,9 Ma y 60,0±0,09 Ma para dos leucogranitos, todas estas dataciones realizadas al norte del Batolito. Tanto (Ordoñez et al. 2001) como (Leal-Mejía 2011), soportan la idea de una edad Paleógena para la parte norte del batolito que difiere química e isotópicamente de la parte sur de esta misma intrusión línea amarilla y morada en (FIGURA 1). Esta nueva fecha pondría al Batolito de Sonsón al menos es su parte central y más oriental en el Paleoceno y no en el Jurásico, lo que reflejaría claramente que se trata de un cuerpo diferente al denotado como Batolito de Sonsón.

Las edades de los cuerpos dómicos relacionados a la actividad reciente se relacionan con dataciones de algunos cuerpos porfiríticos de la zona con edades entre 400.000 años y poco más de 1 Ma (Leal-Mejía, 2011) lo que demuestra la continua actividad del arco magmático de la cordillera central y por lo tanto el continuo aporte de calor a la zona de estudio.

Aunque las fuentes termales se relacionan, por ubicación geográfica y química, con la actividad volcánica del maar de San Diego y El Escondido, en el sector de Puente Linda, donde también hay manantiales termales y, teniendo en cuenta las características de los cuerpos subvolcánicos, no se descarta la presencia de una anomalía térmica (intrusión) sin una expresión superficial.

 

 

CONCLUSIONES

 

 

La actividad ígnea del costado oriental de la Cordillera Central refleja intrusiones de batolitos y stocks graníticos y granodioríticos relacionados a diferentes arcos magmáticos y series calcoalcalina formando magmas tipo I e Intrusivo Gnéisico Triásico tipo S. En general la química de las REE para los cuerpos graníticos muestra enriquecimiento en LREE y empobrecimiento en HREE, se diferencia para rocas gabróicas-dioríticas de la unida Kds un comportamiento más planar (menos diferenciado) en las REE. Anomalías negativas o ligeramente positivas en Eu en la mayoría de muestras son visibles. Estos pulsos registran edades de cristalización desde el Cretáceo inferior para las unidades Kds y Kas con 129,8±9,7 Ma y 132,8±2,6 Ma respectivamente, hasta el Paleógeno con la unidad Tcdf que reporta 54,6±4,4 Ma y el Batolito de Sonsón con 47±11 Ma.

Actividad más reciente representada por los cuerpos porfiríticos de carácter subvolcánico y dómico (volcánicos), de composición andesítica con hornblenda, dacíticos con hornblenda y biotita y el carácter vítreo de algunas muestras, son clara muestra de la continuidad de la actividad magmática cerca de la superficie. Químicamente estas rocas se dan en el campo de las andesitas hasta el límite dacíticoriolítico. La geoquímica posiciona estos cuerpos como vulcanismo de subducción de serie calcoalcalina de carácter adakítico. El comportamiento de las REE es similar en la mayoría con una diferencia notable en la anomalía negativa de Eu para los cuerpos dacíticos y positiva para los demás cuerpos andesíticos (Tad) incluyendo los domos volcánicos, reflejando esto la ausencia de plagioclasa para los primeros y presencia para los segundos, siendo este un mineral que acomoda más Eu en su estructura, y aclarando que la composición química reportada aquí es de roca total que incluye la matriz y los fenocristales. Las dataciones Ar-Ar posicionan los cuerpos dómicos en los 89±4,4 Ka para el domo de San Diego y 153,7±38,5 Ka para uno de los domos del volcán El Escondido, lo que evidencia la continua actividad efusiva en la zona y cuyos resultados preliminarmente suponen que dicha actividad podría ser anterior a los eventos explosivos que produjeron la pila de depósitos volcánicos de la zona.

Los cuerpos de carácter subvolcánico de composiciones acidas e intermedias con intensa alteración, dan evidencian de interacción fluido-roca que reafirman la debilidad estructural y apertura de la roca que proporcionan el paso de fluidos calientes (Falla Rio Dulce) para el sector occidental. De igual manera estructuras como el volcán San Diego podrían ser el resultado de desarrollo de fallas (Palestina, San Diego y Patio Bonito) y la debilidad en los límites entre los cuerpos intrusivos y un basamento pelítico fuertemente foliado. Por otro lado, la presencia de termales en límites de los cuerpos intrusivos o cercanos a ellos, reflejan anomalías térmicas cuyas fuentes podrían ser relacionadas a la evidente actividad volcánica de la zona o a pulsos magmáticos sin expresión superficial.

AGRADECIMIENTOS

Se agradece al Servicio Geológico Colombiano el cual cumpliendo con los objetivos misionales, permitió la realización de este trabajo dentro del marco de los estudios realizados por el Grupo de Investigación y Exploración de Recursos Geotérmicos; se agradece a la geóloga Gina Z. Rodríguez quien hizo parte de las labores de campo del área geotérmica de San Diego cuyos resultados se plasman en el documento de (Rueda y Rodríguez 2016); al geólogo Gabriel Rodríguez del SGC sede Medellín por la asesoría suministrada. A todos los profesionales y técnicos del laboratorio químico del SGC así como a los miembros del grupo de Geotermia que hicieron parte de una u otra forma de este trabajo (Miguel Beltrán, Claudia Alfaro, Iván Ortiz, Maria L. Monsalve, Iván Adarme entre muchos más).

 

 

REFERENCIAS

 

Barrero, D., y Vesga, C. (1976). Geología de la Plancha 188, La Dorada. Mapa Geológico. INGEOMINAS.

Caceres, A. (2012). Genesis of the sediment-hosted uraniferous phosphate deposit in the Berlin Project, Central Cordillera, Colombia and its implications for exploration. M.Sc. Thesis, Queen´s University, Ontario, Canada.

Chappell, B.W., and White, A.J.R. (1992). I- and S-type granites in the Lachlan Fold Belt. Earth and Environmental Science Transactions of the Royal Society of Edinburgh, 83(1-2), 1-26.

 

CHEC. (1968). Proyecto de investigación geotérmica en la Región del Macizo Volcánico del Ruiz. Central Hidroeléctrica de Caldas S.A., Manizales, Colombia.

CHEC. (1981). Investigación geotérmica en la Región del Macizo Volcánico del Ruiz (Colombia). Central Hidroeléctrica de Caldas S.A. Fase I. Informe Final.

Cox, K.G., Bell J.D.,,  and Pankhurst R.J. (1979). The interpretation of igneous rock. London: Allen and Unwin.

Barrero, D. y Vesga J. (1976). Compiladores Geología de la Plancha 188, La Dorada. Mapa Geológico. INGEOMINAS.

Feininger, M.J. and Drummond, M.S. (1990). Derivation of some modern arc magmas by melting of young subducted lithosphere. Nature, 347, 662-665. doi: 10.1038/347662a0.

Feininger, T. (1970). The Palestina fault, Colombia. Geological Society of America Bulletin, 81, 12011216.

Feininger, T., Barrero, D., y Castro, N. (1972). Geología de parte de los departamentos de Antioquia y Caldas (Subzona II-B). Boletín Geológico, 20(2), 1-173.

Feininger, T., Barrero, D., Castro, N., Ramírez, O., Lozano, H., y Vesga, J. (1970). Geología de la Plancha 168, Argelia. Mapa Geológico. INGEOMINAS.

Frost, R.B., Barnes, C.G., Collins, W.J., Arculus, R.J., Ellis, D.J., and Frost, C.D. (2001). A geochemical classification for granitic rocks. Journal of Petrology, 42(11), 2033-2048. doi: 10.1093/petrology/42.11.2033.

Gómez, J., Montes, N.E., Alcárcel, F.F., y Ceballos, J.A. (2015a). Catálogo de dataciones radiométricas de Colombia en ArcGIS y Google Earth. En: J. Gómez, M.F. Almanza. (eds). Compilando la geología de Colombia: Una visión 2015 (pp. 63-419). Servicio Geológico Colombiano, Publicaciones Geológicas Especiales 33. Bogotá, Colombia.

Gómez, J., Montes, N.E., Nivia, Á., y Diederix, H. Compiladores. (2015b). Mapa Geológico de Colombia 2015. Escala 1:1.000.000. Servicio Geológico Colombiano, Bogotá, Colombia.

González, H. (1980). Geología de las planchas 167, Sonsón y 187, Salamina, del mapa geológico de Colombia (Antiguo Cuadrángulo J-8). Boletín Geológico, 23(1), 1-174.

González, H. (1990). Mapa geológico del departamento de Caldas. Escala 1:250.000. Memoria explicativa. INGEOMINAS.

González, H. (2002). Catálogo de unidades litoestratigráficas de Colombia. Batolito de Sonsón (Jts). Instituto de Investigación e Información Geocientífica, Minero-Ambiental y Nuclear.

González, H., Agudelo, S., y Calle, Z. (1980a). Geología de la Plancha 187, Salamina. Mapa Geológico. INGEOMINAS.

González, H., Agudelo, S., y Calle, Z. (1980b). Geología de la Plancha 167, Sonsón. Mapa Geológico. INGEOMINAS.

IGAC. (1961). Carta general del departamento de Caldas. Planchas 188IA, 188IB, 188IC, 188ID, escala 1:25.000. Instituto Geográfico Agustín Codazzi, Bogotá.

IGAC. (1979). Carta general del departamento de Antioquia. Planchas. 168IIIB, escala 1:25.000.  Instituto Geográfico Agustín Codazzi, Bogotá.

Leal-Mejía, H. (2011). Phanerozoic gold metallogeny in the Colombia Andes: A tectono-magmatic approach. Tesis Doctoral, Universitat de Barcelona, España.

Loaiza, D.G. (2012). Caracterización geológica y geoquímica con fines de exploración, análisis petrográfico y metalográfico e identificación de zonas anómalas de mineralizaciones polimetálicas de plomo, zinc, cobre (Pb, Zn, Cu) y metales preciosos: oro, plata (Au-Ag) en el pórfido dacítico (Adakitico) “El Morro” y el “Complejo Ígneo de Samaná” (Granodiorita, Alaskita, Agmatitas) al NE del Municipio de Samaná, Departamento de Caldas. Tesis, Universidad de Caldas, Manizales, Colombia.

López, J., Cuéllar, MA., Aguirre, R., Valencia, M., y Sánchez, A. (2007). Evidencias petrográficas y de campo de una intrusión sintectónica en la Cordillera Central de Colombia: El Caso de la Milonita Granítica del Guacaica. IX Congreso Colombiano de Geología, Bucaramanga, Colombia.

Maniar, P.D., and Piccoli, P.M. (1989). Tectonic discrimination of granitoids. Geological Society of America Bulletin, 101(5), 635-643. doi: 10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2.

Maritza, L., Valencia, J., Ceballos, J., Narváez, L., Pulgarín, B., Correa, A., Murcia, H., Zuluaga, I., Rueda, J., y Pardo, N. (2015). Geología y estratigráfica del Complejo Volcánico Nevado del Ruiz. Informe Técnico. Servicio Geológico Colombiano.

Martin, H. (1999). The adakitic magmas: modern analogues of Archean granitoids. Lithos, 46(3), 411-429. doi: 10.1016/S0024-4937(98)00076-0.

Middlemost, E.A.K. (1994). Naming materials in the magma/igneous rock system. Earth-Science Reviews, 37(3-4), 215-224. doi: 10.1016/00128252(94)90029-9.

Miller, C.F. (1985). Are strongly peraluminous magmas derived from pelitic sedimentary sources?. Journal of Geology, 93(6), 673-689.

Monsalve, M.L. (2015). Vulcanismo del área geotérmica de San Diego. Informe de Avance. Servicio Geológico Colombiano.

Monsalve, M.L., y Ortíz, I. (2013). Reconocimiento vulcanológico área geotérmica de San Diego. Servicio Geológico Colombiano.

Monsalve, M.L., Ortiz, I., and Norini, G. (2014). Deposits associated to San Diego Maar (Colombia). IAVCEI-5th International Maar Conference. Querétaro.

Monsalve, M.L., Ortiz, I., and Norini, G. (2017). El Escondido, a newly identified silicic Quaternary volcano in the NE region of the northern volcanic segment (Central Cordillera of Colombia). Journal of Volcanology and Geothermal Research. doi: 10.1016/j.jvolgeores.2017.12.010.

Nakamura, N. (1974). Determination of REE, Ba, Mg, Na, and K in carbonaceos and ordinary chondrites. Geochemical et Cosmochimical Acta, 38(5), 757775. doi: 10.1016/0016-7037(74)90149-5.

Norman, M., Leeman, W., and Mertzman, S.A. (1992). Granites and rhyolites from the northwestern USA: temporal variation in magmatic processes and relations to tectonic setting. Earth and Environmental Science Transactions of the Royal Society of Edinburgh, 83(1-2), 71-81. doi: 10.1017/S0263593300007768.

Ordoñez, O. (2001). Caracterização isotópica Rb-Sr e Sm-Nd dos principais eventos magmáticos nos Andes Colombianos. Ph.D. thesis, Universidade de Brasilia, Brasilia, Brasil.

Page, W.D. (1986). Seismic Geology and Seismicity of Northwestern Colombia. Compilación ISA. Medellín.

Pearce, J.A., Harris, B.W., and Tindle, A.G. (1984). Trace elements discrimination diagrams for the tectonic interpretation of granitic rocks. Journal of Petrology, 25(4), 956-983. doi: 10.1093/petrology/25.4.956.

Peccerillo, A., and Taylor, S.R. (1976). Geochemistry of Eocene calc-alkaline volcanic rocks from the Kastamonu area, Northern Turkey. Contributions to Mineralogy and Petrology, 58(1), 63-81. doi: 10.1007/BF00384745.

Peña-Ureña, M.L., Muñoz-Rocha, J.A., y Urueña, C.L. (2018). Laboratorio de geocronología en el Servicio Geológico Colombiano: Avances sobre datación U-Pb en circones mediante la técnica LA-ICP-MS. Boletín Geológico, 44, 39-56.

Pérez, G. (1967). Determinación de la edad absoluta de algunas rocas de Antioquia por métodos radioactivos. Dyna, 84, 27-31.

Pimiento, R. (2011). Mineralogía y petrografía de la mineralización de uranio en fosforitas del cretácico inferior, sinclinal de Berlín, Cordillera Central, (departamento de Caldas, Colombia). Tesis, Universidad Industrial de Santander, Bucaramanga, Colombia.

Rueda, J.B. y Rodríguez, G.Z. (2016). Geología del área geotérmica de San Diego. Servicio Geológico Colombiano. Informe Técnico, Bogotá.

 

Ríos-Guerrero, L.J. (2012). Estratigrafía y ambiente de depósito del Sinclinal de Berlín, Municipio de Caldas (Colombia). Tesis de Grado, Universidad Industrial de Santander, Bucaramanga, Colombia.

SGC. (2013). Inventario Nacional de Manantiales. Servicio Geológico Colombiano. Consultado el 14 de diciembre de 2017. http://hidrotermales.sgc.gov.co/

Thiéblemont, D. (1999). Discrimination entre magmatismes calco-alcalins mantellique et crustal, l´exemple des Andes. Comptes Rendus de l'Academie des Sciences- Series IIA - Earth and Planetary Science, 329(4), 243-250. doi: 10.1016/S1251-8050(99)80242-0.

Toro, G. (1989). Caracterización del volcanismo de San Diego y estudios de los depósitos de San Diego (Caldas) y Nariño (Antioquia), Colombia. VI Congreso Colombiano de Geología. Bucaramanga, Colombia.

Toro, G. (1991). Geología de la zona del volcán de San Diego, departamento de Caldas Colombia. I Simposio Magmatismo Andino. Manizales, Colombia.

Vesga, C.J., y Barrero, D. (1978). Edades K/Ar en rocas ígneas y metamórficas de la Cordillera Central de Colombia y su implicación geológica. II Congreso Colombiano de Geología, Bogotá, Colombia.

Vinasco, C.J., Cordani, U.G., González, H., Weber, M., and Pelaez, C. (2006). Geochronological, isotopic, and geochemical data from Permo-Triassic granitic gneisses and granitoids of the Colombian Central Andes. Journal of South American Earth Sciences, 21(4), 355-371. doi: 10.1016/j.jsames.2006.07.007.

Wilt, J.C. (1995). Correspondence of alkalinity and ferric/ferrous ratios of igneous rocks associated with various types of porphyry copper deposits. In: F.W. Pierce, J.G. Bolm (eds.). Porphyry copper deposits of the American Cordillera (pp. 180-200). Arizona Geological Society Digest 20.

Winchester, J.A., and Floyd, P.A. (1977). Geochemical discrimination of different magma series and their differentiation products using immobile elements. Chemical Geology, 20, 325-343. doi:10.1016/0009-2541(77)90057-2.

Woodward-Clyde Consultants. (1980). Prelyminary Seismic Hazard Study San Carlos Proyect Colombia. Report to Interconexion Electrica S.A. and Integral. LTDA. Medellín, Colombia. Stage A, phase I.

 

 [UdW1]