T Uce 0012 58 2017

UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERÍA EN GEOLOGÍA, MINAS,

PETRÓLEOS Y AMBIENTAL
CARRERA DE INGENIERÍA EN GEOLOGÍA
“CARTOGRAFÍA GEOQUÍMICA MULTI-ELEMENTAL DE SEDIMENTOS

FLUVIALES. APLICACIÓN AL ANÁLISIS DEL POTENCIAL
METALOGÉNICO DE LA ZONA DE MACARÁ (ESCALA 1:100 000)”
Trabajo de titulación previo a la obtención del Título de Ingeniero en Geología
Autor: Marco Andrés Vásquez Guañuna
Tutor: M.Sc. Luis Felipe Pilatasig Moreno
Quito, noviembre 2017

DEDICATORIA
Este trabajo está dedicado de manera especial a mi madre, quien con tanto sacrificio y
esfuerzo me ha brindado la mejor herencia que me pudo dar, la educación.
A mis hermanas Andrea y Karina quienes siempre me han brindado su apoyo y cariño, espero
poder ser ejemplo para ellas.
A mis amig@s Mary, Lizeth, Lissette, Anita, Dyana, Jair, Frank, Alex, Andrés y Roberto con
quienes compartí tantos momentos, alegrías y tristezas; sin ellos este largo camino no
hubiese sido el mismo.
ii
AGRADECIMIENTOS
Agradezco a la Universidad Central del Ecuador, a la FIGEMPA y a cada uno de los
docentes que de manera desinteresada me impartieron sus conocimientos, experiencias y
consejos para mi formación académica.
A todos quienes conforman el INIGEMM por brindarme la oportunidad de desarrollar este

trabajo.
De manera especial agradezco al Ing. Luis Pilatasig por su tiempo, consejo y apoyo en la
realización de este proyecto.
Al Ing. Danny Burbano por brindarme su apoyo, amistad y consejo.
A la Ing. Rocío Sangucho y al Ing. Salomón Brito por su apoyo y consejo.
Por último, agradezco a todos aquellos que de una u otra manera contribuyeron a mi
formación personal, profesional y a la realización de este trabajo.
iii
DERECHOS DE AUTOR
Yo, MARCO ANDRÉS VÁSQUEZ GUAÑUNA, en calidad de autor y titular de los derechos
morales y patrimoniales del trabajo de titulación “CARTOGRAFÍA GEOQUÍMICA
MULTI-ELEMENTAL DE SEDIMENTOS FLUVIALES. APLICACIÓN AL
ANÁLISIS DEL POTENCIAL METALOGÉNICO DE LA ZONA DE MACARÁ
(ESCALA 1:100 000)”, modalidad Proyecto de Investigación, de conformidad con el Art.
114 del CÓDIGO ORGÁNICO DE LA ECONOMÍA SOCIAL DE LOS
CONOCIMIENTOS, CREATIVIDAD E INNOVACIÓN, concedo a favor de la Universidad
Central del Ecuador una licencia gratuita, intransferible y no exclusiva para el uso no
comercial de la obra, con fines estrictamente académicos. Conservo a mi favor todos los
derechos de autor sobre la obra, establecidos en la normativa citada. Así mismo, autorizo a la
Universidad Central del Ecuador para que realice la digitalización y publicación de este
trabajo de titulación en el repositorio virtual, de conformidad a lo dispuesto en el Art. 144 de
la Ley Orgánica de Educación Superior.
El autor declara que la obra objeto de la presente autorización es original en su forma de

expresión y no infringe el derecho de autor de terceros, asumiendo la responsabilidad por
cualquier reclamación que pudiera presentarse por esta causa y liberando a la Universidad de
toda responsabilidad
Marco Andrés Vásquez Guañuna
CC. 1720695467
E-mail: marcovasquezgeo@gmail.com
iv
APROBACIÓN DEL TRABAJO DE TITULACIÓN POR PARTE DEL TUTOR
Yo, Luis Felipe Pilatasig Moreno, en mi calidad de tutor del trabajo de titulación, modalidad
Proyecto de Investigación; cuyo título es: “CARTOGRAFÍA GEOQUÍMICA MULTI-
ELEMENTAL DE SEDIMENTOS FLUVIALES. APLICACIÓN AL ANÁLISIS DEL
POTENCIAL METALOGÉNICO DE LA ZONA DE MACARÁ (ESCALA 1:100 000)”,
elaborado por el señor MARCO ANDRÉS VÁSQUEZ GUAÑUNA, previo a la obtención del
Título de Ingeniero en Geología; considero que el mismo reúne los requisitos y méritos
necesarios en el campo metodológico y en el campo epistemológico, para ser sometido a la
evaluación por parte del tribunal examinador que se designe, por lo que APRUEBO, a fin de
que el trabajo investigativo sea habilitado para continuar con el proceso de titulación
determinado por la Universidad Central del Ecuador.
En la ciudad de Quito, a los 01 días del mes de noviembre del 2017
____________________________
M.Sc. Luis Felipe Pilatasig Moreno

Ingeniero en Geología
CC. 0501382493
TUTOR
v
APROBACIÓN DEL TRABAJO DE TITULACIÓN POR PARTE DEL TRIBUNAL
El tribunal constituido por: Ing. Alex Mateus, presidente del tribunal de grado oral; Ing. Elías
Ibadango e Ing. Danny Burbano como miembros, DECLARAN: Que el presente proyecto de
investigación denominado “CARTOGRAFÍA GEOQUÍMICA MULTI-ELEMENTAL
DE SEDIMENTOS FLUVIALES. APLICACIÓN AL ANÁLISIS DEL POTENCIAL
METALOGÉNICO DE LA ZONA DE MACARÁ (ESCALA 1:100 000)”, elaborado por
el señor VÁSQUEZ GUAÑUNA MARCO ANDRÉS, Egresado de la Carrera de Ingeniería
en Geología, ha sido revisado, verificado y evaluado detenida y legalmente, por lo que
califican como original y auténtico del autor.
En la ciudad de Quito, a los 20 días del mes de noviembre del 2017.
Para constancia de lo actuado firman:
_________________________
Ing. Alex Mateus
PRESIDENTE DEL TRIBUNAL
______________________ _________________________
Ing. Elías Ibadango Ing. Danny Burbano
MIEMBRO DEL TRIBUNAL MIEMBRO DEL TRIBUNAL
vi
ÍNDICE DE CONTENIDO
SIGLAS Y ABREVIATURAS ............................................................................................................ xiii

RESUMEN ........................................................................................................................................... xiv
ABSTRACT .......................................................................................................................................... xv
1. INTRODUCCIÓN .......................................................................................................................... 1
1.1. Estudios previos ...................................................................................................................... 1
1.2. Justificación ............................................................................................................................. 2
1.3. Objetivos ................................................................................................................................. 3
1.3.1. Objetivo general .............................................................................................................. 3
1.3.2. Objetivos específicos ....................................................................................................... 3
1.4. Alcance .................................................................................................................................... 4
1.5. Zona de estudio ....................................................................................................................... 4
1.5.1. Ubicación y acceso .......................................................................................................... 4
1.5.2. Clima ............................................................................................................................... 4
1.5.3. Geomorfología................................................................................................................. 5
1.5.4. Hidrografía ...................................................................................................................... 6
2. CONTEXTO GEOLÓGICO ........................................................................................................... 8
2.1. Marco Geodinámico del Ecuador y Norte de Perú .................................................................. 8
2.2. Marco Geológico Regional ..................................................................................................... 9
2.2.1. Dominio Loja-Olmos ...................................................................................................... 9
2.2.2. Cuenca Lancones-Alamor ............................................................................................... 9
2.3. Litoestratigrafía del área de estudio ...................................................................................... 10
2.3.1. Unidad Los Pilares ........................................................................................................ 10
2.3.2. Unidad Punta de Piedra ................................................................................................. 11
2.3.3. Formación Sacapalca ..................................................................................................... 11
2.3.4. Rocas Intrusivas ............................................................................................................ 12
2.4. Geología Estructural .............................................................................................................. 14
2.5. Geología Económica ............................................................................................................. 14
2.5.1. Concesiones Mineras ..................................................................................................... 14
2.5.2. Metalogenia del Sur del Ecuador .................................................................................. 15
3. MARCO TEÓRICO ...................................................................................................................... 18
3.1. Exploración Geoquímica ....................................................................................................... 18
3.2. Parámetros Geoquímicos ....................................................................................................... 18
vii
3.3. Elementos Guía ..................................................................................................................... 18
3.4. Tipos de mineralización ........................................................................................................ 19
3.4.1. Epitermales .................................................................................................................... 19
3.4.2. Pórfidos ......................................................................................................................... 21
3.4.3. Sulfuros Masivos Volcanogénicos (VMS) .................................................................... 22
3.5. Geoestadística........................................................................................................................ 23
3.5.1. Medidas de tendencia central ........................................................................................ 24
3.5.2. Medidas de tendencia no central ................................................................................... 24
3.5.3. Medidas de dispersión ................................................................................................... 24
3.5.4. Medidas de Forma ......................................................................................................... 25
3.5.5. Normalidad de la distribución de datos ......................................................................... 25
3.5.6. Variograma y Semivariograma ...................................................................................... 25
3.5.7. Modelos Teóricos de Semivariograma .......................................................................... 26
3.5.8. Interpolación.................................................................................................................. 28
3.5.9. Kriging........................................................................................................................... 28
4. MARCO METODOLÓGICO ....................................................................................................... 30
4.1. Muestreo ................................................................................................................................ 30
4.2. Análisis de laboratorio........................................................................................................... 30
4.3. Análisis Geoestadístico de Datos Geoquímicos .................................................................... 30
4.4. Análisis Geoestadístico mono-elemental .............................................................................. 31
4.4.1. Análisis Exploratorio de Datos Espaciales (AEDE)...................................................... 31
4.4.2. Tratamiento de datos ..................................................................................................... 31
4.4.3. Transformación de datos ............................................................................................... 31
4.4.4. Tratamiento de los valores atípicos (outliers)................................................................ 31
4.4.5. Análisis de tendencia de los datos ................................................................................. 32
4.4.6. Modelamiento del semivariograma ............................................................................... 32
4.4.7. Validación del método de interpolación ........................................................................ 32
4.4.8. Determinación de umbrales geoquímicos...................................................................... 32
4.4.9. Definición de anomalías de los elementos .................................................................... 33
4.5. Análisis Geoestadístico multi-elemental ............................................................................... 33
5. PRESENTACIÓN Y ANÁLISIS DE DATOS ............................................................................. 34
5.1. Geología ................................................................................................................................ 34
5.1.1. Ambientes Geológicos .................................................................................................. 34
5.1.2. Indicios de mineralización ............................................................................................. 34
5.1.3. Indicios de Oro en batea ................................................................................................ 35
5.1.4. Alteraciones ................................................................................................................... 36
viii
5.2. Geoquímica ........................................................................................................................... 37
5.3. Análisis Geoestadístico General ............................................................................................ 38
5.3.1. Análisis exploratorio de datos espaciales (AEDE) ........................................................ 38
5.3.2. Medidas estadísticas de los datos transformados .......................................................... 40
5.3.3. Tendencia de los datos .................................................................................................. 42
5.4. Análisis Geoestadístico por Ambientes Geológicos.............................................................. 42
5.4.1. Ambiente Volcánico y Magmático Cretácico (Ambiente 1) ......................................... 42
5.4.2. Ambiente Volcánico Continental Paleógeno (Ambiente 2) .......................................... 48
6. RESULTADOS ............................................................................................................................. 55
6.1. Resultados - Análisis estadístico general .............................................................................. 55
6.1.1. Distribución de los elementos ....................................................................................... 55
6.2. Resultados - Análisis estadístico por Ambientes Geológicos ............................................... 61
6.2.1. Interpolación de datos ................................................................................................... 61
6.2.2. Anomalías geoquímicas - Ambiente Volcánico y Magmático Cretácico ...................... 61
6.2.3. Anomalías geoquímicas - Ambiente Volcánico Continental Paleógeno ....................... 68
6.3. Calidad del agua del sitio de muestreo .................................................................................. 73
7. INTERPRETACIÓN Y DISCUSIÓN........................................................................................... 74
7.1. Áreas con Potencial Metalogénico ........................................................................................ 74
7.2. Discusión ............................................................................................................................... 83
8. CONCLUSIONES Y RECOMENDACIONES ............................................................................ 87
8.1. Conclusiones ......................................................................................................................... 87
8.2. Recomendaciones .................................................................................................................. 89
9. BIBLIOGRAFÍA ........................................................................................................................... 91
10. ANEXOS................................................................................................................................... 95
ix
LISTA DE FIGURAS
Figura 1: Mapa de ubicación geográfica del área de estudio. ................................................................. 5

Figura 2: Modelo digital elevaciones del área de estudio. ...................................................................... 6
Figura 3: Mapa de unidades geomorfológicas del área de estudio.. ........................................................ 6
Figura 4: Mapa de las cuencas hidrográficas .......................................................................................... 7
Figura 5. Mapa de sub-cuencas y red hidrográfica del área de estudio ................................................... 7
Figura 6: Marco Geodinámico actual del Ecuador .................................................................................. 8
Figura 7: Dominios Litotectónicos del sur de Ecuador y norte de Perú. ................................................. 9
Figura 8: Columna estratigráfica de la Cuenca Lancones-Alamor........................................................ 10
Figura 9: Mapa Geológico de Macará Escala 1:100 000. ..................................................................... 13
Figura 10: Mapa de estructuras, frontera sur del Ecuador..................................................................... 14
Figura 11: Ubicación de concesiones mineras. ..................................................................................... 15
Figura 12: Geología Regional y Ocurrencias Minerales del sur del Ecuador y norte de Perú. ............. 17
Figura 13: Modelos de mineralización de para epitermales Au-Cu de la Costa del Pacífico................ 21
Figura 14: Modelo de mineralización y alteración de un depósito tipo pórfido.................................... 22
Figura 15: Modelo de mineralización, alteración y zonación de elementos para VMS. ....................... 23
Figura 16: Comportamiento típico de un semivariograma .................................................................... 27
Figura 17: Comparación de los modelos Exponencial, Esférico y Gaussiano. ..................................... 28
Figura 18: Parámetros Geoquímicos ..................................................................................................... 33
Figura 19. Mapa de ubicación de muestras de sedimentos fluviales por Ambientes Geológicos ......... 34
Figura 20: Mapa de ubicación de indicios minerales y oro aluvial. ...................................................... 36
Figura 21: Mapa de alteraciones. .......................................................................................................... 37
Figura 22: Mapa de ubicación de muestras de sedimentos fluviales ..................................................... 37
Figura 23. Diagramas Box-Plot. Análisis General ................................................................................ 38
Figura 24: Histogramas de frecuencia de los elementos. Análisis General........................................... 40
Figura 25. Histogramas de frecuencia de los elementos transformados a Log 10. Análisis General .... 42
Figura 26. Diagramas Box-Plot de los elementos. Ambiente 1 y Ambiente 2. ..................................... 43
Figura 27: Histogramas de frecuencia. Ambiente 1 .............................................................................. 44
Figura 28. Histogramas de frecuencia de los elementos transformados a Log 10. Ambiente 1 ............ 46
Figura 29. Semivariogramas experimentales. Ambiente 1 .................................................................... 47
Figura 30: Histogramas de frecuencia. Ambiente 2 .............................................................................. 50
Figura 31. Histogramas de frecuencia de los elementos. Ambiente 2................................................... 51
Figura 32. Semivariogramas experimentales. Ambiente 2 .................................................................... 53
Figura 33: Mapa de distribución de Cobre clasificado por cuantiles. Análisis General........................ 55
x
Figura 34: Mapa de distribución de Plomo clasificado por cuantiles. Análisis General ....................... 56
Figura 35: Mapa de distribución de Zinc clasificado por cuantiles. Análisis General .......................... 56
Figura 36: Mapa de distribución de Oro clasificado por cuantiles. Análisis General ........................... 57
Figura 37: Mapa de distribución de Plata clasificado por cuantiles. Análisis General ......................... 57
Figura 38: Mapa de distribución de Mercurio clasificado por cuantiles. Análisis General .................. 58
Figura 39: Mapa de distribución de Arsénico clasificado por cuantiles. Análisis General ................... 58
Figura 40: Mapa de distribución de Antimonio clasificado por cuantiles. Análisis General ................ 59
Figura 41: Mapa de distribución de Molibdeno clasificado por cuantiles. Análisis General ................ 59
Figura 42: Mapa de distribución de Wolframio clasificado por cuantiles. Análisis General ................ 60
Figura 43: Mapa de distribución de Bario clasificado por cuantiles. Análisis General ........................ 60
Figura 44: Mapa de estimación de concentración de Cobre. Ambiente 1 ............................................. 62
Figura 45: Mapa de estimación de concentración de Plomo. Ambiente 1 ............................................ 62
Figura 46: Mapa de estimación de concentración de Zinc. Ambiente 1 ............................................... 63
Figura 47: Mapa de estimación de concentración de Oro. Ambiente 1 ................................................ 64
Figura 48: Mapa de estimación de concentración de Plata. Ambiente 1 ............................................... 64
Figura 49: Mapa de estimación de concentración de Mercurio. Ambiente 1 ........................................ 65
Figura 50: Mapa de estimación de concentración de Arsénico. Ambiente 1 ........................................ 66
Figura 51: Mapa de estimación de concentración de Antimonio. Ambiente 1 ..................................... 66
Figura 52: Mapa de estimación de concentración de Molibdeno. Ambiente 1 ..................................... 67
Figura 53: Mapa de estimación de concentración de Wolframio. Ambiente 1 ..................................... 67
Figura 54: Mapa de estimación de concentración de Bario. Ambiente 1 .............................................. 68
Figura 55: Mapa de estimación de concentración de Cobre y Plomo. Ambiente 2............................... 69
Figura 56: Mapa de estimación de concentración de Zinc y Oro. Ambiente 2 ..................................... 70
Figura 57: Mapa de estimación de concentración de Plata y Mercurio. Ambiente 2 ............................ 71
Figura 58: Mapa de estimación de concentración de Arsénico y Antimonio. Ambiente 2 ................... 71
Figura 59: Mapa de estimación de concentración de Molibdeno y Wolframio. Ambiente 2 ................ 72
Figura 60: Mapa de estimación de concentración de Bario. Ambiente 2 .............................................. 73
Figura 61: Mapa de calidad del agua del sitio de muestreo................................................................... 73
Figura 62: Mapa de anomalías geoquímicas de Macará. Ambientes 1 y 2 ........................................... 75
Figura 63. Mapa de áreas prospectivas con potencial metalogénico de Macará ................................... 75
xi
LISTA DE TABLAS
Tabla 1. Minas y prospectos metálicos del Sur del Ecuador y Norte de Perú. ...................................... 16
Tabla 2. Asociaciones de elementos y pathfinders ................................................................................ 19
Tabla 3. Principales características de los depósitos epitermales HS y LS ........................................... 20
Tabla 4. Tipos de predictores Kriging y sus propiedades...................................................................... 29
Tabla 5. Indicios minerales en la zona de Macará................................................................................. 35
Tabla 6. Resumen de unidades litoestratigráficas y número de muestras por ambientes geológicos.... 38
Tabla 7. Medidas estadísticas de los elementos. Análisis General ........................................................ 39
Tabla 8. Medidas estadísticas de los elementos transformados a Log10. Análisis General .................. 41
Tabla 9. Medidas estadísticas de los elementos. Ambiente 1................................................................ 43
Tabla 10. Medidas estadísticas de los elementos transformados a Log10. Ambiente 1 ........................ 45
Tabla 11. Parámetros de los modelos de semivariogramas elegidos. Ambiente 1 ................................ 47
Tabla 12. Parámetros de los modelos de semivariogramas elegidos Ambiente 1 ................................. 48
Tabla 13. Medidas estadísticas de los elementos. Ambiente 2 .............................................................. 49
Tabla 14. Medidas estadísticas de los elementos transformados a Log10. Ambiente 2 ........................ 51
Tabla 17. Valor de fondo y umbrales anomálicos. Análisis General y por Ambientes Geológicos...... 61
Tabla 18. Áreas prospectivas de Macará y sus asociaciones elementales. ............................................ 74
Tabla 19. Matriz de Correlación de Elementos sector Q. Mandalá ....................................................... 76
Tabla 20. Matriz de Correlación de Elementos sector Sabiango........................................................... 77
Tabla 21. Matriz de Correlación de Elementos sector San Vicente ...................................................... 78
Tabla 21. Matriz de Correlación de Elementos sector Punta de Piedra ................................................ 78
Tabla 23. Matriz de Correlación de Elementos sector Las Vigas ......................................................... 79
Tabla 24. Matriz de Correlación de Elementos sector Sauce ................................................................ 80
Tabla 25. Matriz de Correlación de Elementos sector Llanitos ............................................................ 81
Tabla 26. Matriz de Correlación de Elementos sector Nongora............................................................ 82
Tabla 27. Matriz de Correlación de Elementos sector Trapiche ........................................................... 83
Tabla 28. Anomalías geoquímicas aisladas. Ambientes 1 y 2............................................................... 83
xii
LISTA DE ANEXOS
Anexo 1. Formulas usadas en estadística .............................................................................................. 96

Anexo 2. Gráficos de análisis de tendencia de los datos. Análisis General .......................................... 97
Anexo 3. Gráficos de análisis de tendencia. Ambiente 1 ...................................................................... 98
Anexo 4. Validación cruzada de los elementos. Ambiente 1 ................................................................ 99
Anexo 5. Gráficos de análisis de tendencia. Ambiente 2 .................................................................... 100
Anexo 6. Validación cruzada de los elementos. Ambiente 2 .............................................................. 101
Anexo 7. Umbrales Anomálicos. Análisis General ............................................................................. 102
Anexo 8. Umbrales Anomálicos. Ambiente 1 ..................................................................................... 103
Anexo 9. Umbrales Anomálicos. Ambiente 2 ..................................................................................... 104
Anexo 10. Resultados de análisis químicos de indicios de mineralización en Macará ....................... 105
SIGLAS Y ABREVIATURAS
AEDE: Análisis Exploratorio de Datos Espaciales
ARCOM: Agencia de Regulación y Control Minero
BAT: Bloque Amotape Tahuín
DINAGE: Dirección Nacional de Geología de Ecuador
Fm: Formación
ICP-OES: Espectroscopía de Emisión Atómica con Plasma de Acoplamiento
INEMIN: Instituto Ecuatoriano de Minería
INGEMMET: Instituto Geológico, Minero y Metalúrgico de Perú
INIGEMM: Instituto Nacional de Investigación Geológico, Minero y Metalúrgico
Ma: Millones de años
ppb: partes por billón
ppm: partes por millón
Q(s): Quebrada(s)
U: Unidad (Litoestratigráfica)
VMS: Sulfuros Masivos Volcanogénicos (siglas en inglés)
xiii
TEMA: “Cartografía geoquímica multi-elemental de sedimentos fluviales. Aplicación al
análisis del potencial metalogénico de la zona de Macará (Escala 1:100 000)”
Tutor: Luis Felipe Pilatasig Moreno
RESUMEN
Este estudio presenta los resultados de la exploración y cartografía geoquímica multi-

elemental llevada a cabo a partir de 661 muestras de sedimentos fluviales de la zona de
Macará (sur del Ecuador), recolectados y analizados en los laboratorios del INIGEMM.
Geológicamente la zona de Macará se encuentra emplazada en el borde oriental de la Cuenca

Lancones-Alamor y está constituida por rocas volcánicas de afinidad toleítica relacionadas a
un arco de islas del cretácico (Unidad Punta de Piedra) y rocas volcánicas y volcanoclásticas
asociadas a un arco volcánico continental del paleógeno (Formación Sacapalca). Además de
rocas intrusivas de composición intermedia a ácida de afinidad calco-alcalina emplazadas en
el cretácico superior, parte del Batolito de Tangula y del Complejo Peribatolítico.
Primero se realizó un análisis geoestadístico global y luego por ambientes geológicos para los
elementos: Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba, con lo cual se determinó el umbral
anomálico y valor de fondo geoquímico para cada elemento en ambos escenarios. Las
anomalías obtenidas han sido correlacionadas con la litología, estructuras, indicios de
mineralización y alteraciones, lo que permitió establecer el contexto metalogénico de la zona
fronteriza sur del Ecuador y definir 16 áreas con potencial presencia de recursos minerales.
La zona occidental del área, litológicamente constituida por rocas predominantemente

volcánicas de afinidad toleítica, presenta tres áreas favorables para la ocurrencia de depósitos
tipo VMS relacionadas a asociaciones de los elementos Pb-Zn-Ag-Ba-Hg y Cu-Pb-Ba, entre
las cuales está Sabiango. Dentro de esta misma zona, se determinó cinco áreas prospectivas
con afinidad a mineralización epitermal por las asociaciones Cu-Zn-Sb, Cu-Ag; Au-Hg-Sb-
W, Hg-As-Sb-W, Zn-Pb-Hg y Cu-Zn-Mo-Hg-As; y tres áreas con potencial para
mineralización tipo pórfido de Cu y Cu-Mo por las asociaciones de Cu-Mo-Ag-Ba-As, Au-W
y Ag-Mo; genéticamente podrían estar relacionadas a cuerpos intrusivos menores de
composición diorítica y pórfido-andesítica. Dentro de la franja volcánica oriental se
definieron asociaciones de Cu-Zn-Sb, Au-Pb-As y Au-Zn-Ba-W que delimitan cinco áreas
prospectivas con predisposición a la ocurrencia de depósitos tipo epitermal y pórfido
posiblemente hospedados en secuencias volcánicas y volcanoclásticas.
PALABRAS CLAVE: MACARÁ / RÍO CALVAS / CARTOGRAFÍA GEOQUÍMICA /

GEOESTADÍSTICA / GEOQUÍMICA MILTI-ELEMENTAL / POTENCIAL
METALOGÉNICO.
xiv
TITLE: “Multi-elemental geochemical mapping of stream sediments. Application to
assessment of mineral resources potentiality in Macará area (Scale 1:100 000)”
Tutor: Luis Felipe Pilatasig Moreno
ABSTRACT
This study presents the main results of the exploration and multi-elemental
geochemical mapping carried out from 661 stream sediments samples over Macará area,
(south of Ecuador), collected and analyzed by INIGEMM.
Macará area is located on the eastern edge of Lancones-Alamor Basin and it is constituted by
volcanic rocks of tholeiitic affinity related to a cretaceous island arc (Punta de Piedra Unit)
and volcanic and volcanoclastic rocks related to a paleogen continental volcanic arc
(Sacapalca Formation). Besides, intermediate to acid composition intrusive rocks of the upper
cretaceous Tangula Batholith and the Peri-batholitic Complex.
A global and by geological environments geostatistical analysis were performed for the
elements Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W and Ba in order to determine the threshold
and geochemical background values for each element in both scenarios. The resulting
anomalies were correlated with lithology, structures, mineralization and hydrothermal
alterations, which allowed assessing the metallogenic framework of south border of Ecuador
and locate 16 areas with mineral resources potentiality.
The western zone, lithologically constituted by tholeiitic volcanic rocks, exhibits three areas
with potential for VMS deposits exposed by associations of Pb-Zn-Ag-Ba-Hg and Cu-Pb-Ba.
Five prospective zones for epithermal mineralization have been determined by multi-
elemental associations of Cu-Zn-Sb, Cu-Ag, Au-Hg-Sb-W, Hg-As-Sb-W, Zn-Pb-Hg and Cu-
Zn-Mo-Hg-As within this same zone; in addition, three areas for Cu and Cu-Mo porphyries
by associations of Cu-Mo-Ag-Ba-As, Au-W y Ag-Mo were defined. Both could be
genetically related to dioritic and porphyry-andesitic minor intrusions. Within the eastern
volcanic stripe, associations of Cu-Zn-Sb, Au-Pb-As and Au-Zn-Ba-W revealed five
prospective areas for porphyries and epithermal deposits probably hosted in volcanic and
volcanoclastic sequences.
KEY WORDS: MACARÁ / RÍO CALVAS / MULTI-ELEMENTAL GEOCHEMISTRY

/GEOCHEMICAL MAPPING / GEOESTATISTICS / METALLOGENIC POTENTIAL.
I CERTIFY that the above and foregoing is a true and correct translation of the original
document in Spanish.
Ing. Luis Felipe Pilatasig Moreno

Certified Translator/Tutor
ID: 0501382493
xv
1. INTRODUCCIÓN
1.1. Estudios previos

Los primeros bosquejos de la geología del sur del Ecuador fueron descritos por Wolf
(1897) y Sauer (1965). Posteriormente surgieron nuevos datos importantes, en particular de
las investigaciones geológico-mineras realizados por Villemur y la Misión Francesa en 1967.
Los trabajos de prospección geológica y geoquímica en la provincia de Loja inician en 1970

con el Programa de Exploración Sistemática de Geoquímica de la Naciones Unidas
(Operación 8). Dentro de las zonas de Macará y Rio Calvas, se analizaron 161 muestras de
sedimentos fluviales, que arrojaron concentraciones anómalas aisladas para Cu, Pb y Zn.
Análisis de sedimentos pesados determinaron la presencia de minerales de Pb, Sn y W.
Entre 1969 y 1970, como una asistencia técnica del Instituto de Ciencias Geológicas de
Londres al Servicio Nacional de Geología y Minería de Ecuador, se realizó un mapeo
fotogeológico en las provincias de Loja, El Oro, Azuay y Zamora; además, se integró una
componente de investigación minera. El análisis estadístico de 142 muestras de sedimentos
fluviales recolectadas en la zona de Macará, arrojaron valores anómalos para Cu y Zn en los
sectores del Cerro Gualanga, Algarrobillo y El Empalme. Este trabajo sirvió de base para la
elaboración de las hojas geológicas a escala 1:100 000 de Alamor, Zapotillo, Cariamanga y
Macará publicadas entre 1973 y 1974.
Entre 1987 y 1988, el Instituto Ecuatoriano de Minería (INEMIN) lleva a cabo el Sub-
proyecto Macará de Investigaciones Regionales y Estudios Específicos. Un programa de
prospección geológica y geoquímica alrededor del cantón Macará que se centró en los
sectores de Laguar, Sabiango y Angashcola, tras considerar las posibilidades que presentaban
estos sitios para la localización de depósitos primarios de oro y sulfuros masivos.
Entre los trabajos de investigación geológica más recientes se encuentra el realizado por
Jaillard et al. (1996), que describe la estratigrafía de la Cuenca Alamor-Lancones y la
correlación con su prolongación en Perú. Egüez y Poma (2001), proponen dos columnas
estratigráficas de las series occidental y oriental de la cuenca, además aportan datos
relacionados a la naturaleza de las rocas volcánicas cretácicas de la zona de Macará,
redefiniéndolas como Formación Punta de Piedra.
1
En 2005 la Dirección Nacional de Geología de Ecuador (DINAGE) y el Instituto Geológico,
Minero y Metalúrgico de Perú (INGEMMET), ejecutan el Proyecto Multinacional Andino:
Geociencias para las comunidades Andinas. Geología del Ecuador y Perú entre 3°S y 6°S. En
esta investigación se divide a la zona estudiada en Dominios Litotectónicos, según
características litológicas, estructurales y genéticas. La investigación aporta descripciones
geológicas detalladas de la zona de Macará, análisis geoquímicos para óxidos, elementos traza
y tierras raras de muestras del Batolito de Tangula y Unidad Punta de Piedra con el fin de
establecer el ambiente petrogenético de las formaciones.
En 2007 se realiza la actualización de la Hoja Geológica de Macará a escala 1:100 000

elaborada por Pilatasig L. y Gordón D. Este mapa geológico es el más actual y constituye la
principal referencia para la elaboración del modelo geológico y síntesis litoestratigráfica de la
presente investigación. En el mapa se describen las Unidades Pilares, Punta de Piedra,
Formación Sacapalca y el Batolito de Tangula.
Por el lado peruano, el INGEMMET desarrolla entre 2005 y 2012 el proyecto: Metalogenia,
Geología Económica y Potencial Minero de la Deflexión de Huancabamba, noroeste del Perú.
El área del estudio limita al norte con territorio ecuatoriano; precisamente con el área
seleccionada para la presente investigación. El estudio analiza las principales características
geológicas y metalogénicas del norte de Perú.
En 2013 el Instituto Nacional de Investigación Geológico Minero Metalúrgico (INIGEMM)

desarrolla el Proyecto: “Modelo piloto para la determinación del potencial geológico-
mineralógico de las zonas Zaruma y Cariamanga a escala 1:100 000”. El trabajo aporta
descripciones detalladas de las Unidades Punta de Piedra, Celica, Sacapalca y Batolito de
Tangula; correlacionables con la geología de la zona de Macará.
1.2. Justificación
El Ecuador cuenta con recursos minerales tanto metálicos como no metálicos los cuales
se han venido explotando a manera de minería artesanal y pequeña minería, sin aportar un
gran valor económico al país. En los últimos años se han investigado y consolidado proyectos
mineros de mediana y gran escala lo que confirma el potencial minero del país; por esta razón
es importante incrementar el conocimiento y la investigación geocientífica que permita
establecer la existencia de yacimientos minerales y justifique la explotación técnica de los
2
mismos con el objetivo de posibilitar la generación de beneficios económicos y fuentes de
trabajo para el desarrollo del país.
Desde el año 2013, el INIGEMM ejecuta el proyecto: “Investigación Geológica y

Disponibilidad de Ocurrencias de Recursos Minerales en el Territorio Ecuatoriano”, con el
objeto de incrementar el conocimiento geológico del país y promover el desarrollo sostenible
y sustentable de estos recursos. Como parte de este proyecto, durante el primer semestre de
2015, el grupo de Prospección Geoquímica del INIGEMM, llevó a cabo una campaña de
muestreo de sedimentos fluviales en el área que corresponde a las hojas topográficas de
Macará y Rio Calvas a escala 1:50 000. Las 661 muestras recolectadas fueron analizadas en el
laboratorio químico del Instituto y sus resultados se encuentran disponibles. En el año 2016,
el grupo de Ocurrencias Minerales desarrolló el inventario de mineralizaciones y alteraciones
en la zona.
Toda esta información actual, sumada a información histórica recopilada, será integrada y
correlacionada en la presente investigación, con el fin de delimitar sectores potencialmente
prospectivos y contribuir al incremento del conocimiento geológico y minero del Ecuador.
1.3. Objetivos
1.3.1. Objetivo general

Determinar anomalías y asociaciones geoquímicas correlacionables con la geología,
mineralización y alteraciones; con el fin de establecer áreas con potencial metalogénico
en la zona de Macará.
1.3.2. Objetivos específicos

• Analizar estadísticamente datos químicos de sedimentos fluviales para los elementos
Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba.
• Determinar el background, threshold y establecer zonas geoquímicamente anómalas
en el área de estudio.
• Analizar e identificar los ambientes geológicos que constituyen el área de estudio y los
tipos de depósitos minerales relacionados con cada ambiente.
• Correlacionar la geología y geoquímica del área de estudio.
• Elaborar mapas de anomalías geoquímicas mono-elementales, multi-elementales y de
áreas con potencial metalogénico de Macará.
3
1.4. Alcance
El procesamiento e interpretación geoestadística de los datos, permitió delimitar
anomalías geoquímicas correlacionables con la geología y mineralización de la zona con el fin
de estudiar el potencial metalogénico de Macará.
Los resultados obtenidos de la presente investigación contribuyeron al conocimiento del

potencial metalogénico del sur del Ecuador, permitieron delimitar áreas con potencial
económico en la zona que servirán de base para futuras campañas de exploración minera.
Los productos de la presente investigación son: mapas de anomalías geoquímicas mono-

elementales correspondientes a los elementos Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba,
mapas de anomalías geoquímicas multi-elementales y mapa de áreas con potencial
metalogénico de Macará.
1.5. Zona de estudio
1.5.1. Ubicación y acceso

El área de estudio (742.6km2) se encuentra en el extremo sur de la provincia de Loja,
abarca parte de los cantones Macará, Sozoranga, Calvas; y cubre las hojas topográficas de
Macará y Rio Calvas (Escala 1:50 000) que limitan al sur con el Perú. Figura 1
A la zona de Macará se accede por la vía Troncal de la Costa (E-25) hasta la ciudad de
Machala, Provincia de El Oro, para continuar por la vía de primer orden Machala-Alamor-
Celica-Macará. También se puede acceder por la vía Troncal de la Sierra (E-35) hasta la
ciudad de Loja, para tomar la carretera de primer orden Loja-Catacocha-Macará.
1.5.2. Clima
El clima del área está condicionado por la topografía y la altitud de las diversas zonas;
en consecuencia, se presentan dos tipos de climas. En el sector de Macará el clima es cálido
subtropical a seco, con temperaturas que varían entre 25°C y 30°C que se registran al medio
día, sobre todo en invierno. Por otro lado, en el sector de Río Calvas prevalece un clima
templado subhúmedo. Las lluvias en la zona son moderadas y varían en promedio con
precipitaciones de 1950mm/año.
El régimen de mayor pluviosidad corresponde a la estación invernal que se presenta entre

noviembre y mayo, mientras que de junio a octubre se registra la estación seca sin lluvias.
4
Esta variación en la humedad históricamente ha generado déficit hídrico en la zona
principalmente en la época seca (INAMHI, 2012).
Figura 1: a) Mapa de ubicación geográfica del área de estudio. b) Hoja de Macará escala 1:100 000, constituida
por las hojas de Macará, Río Calvas y NVII-C4 Escala 1:50 000. c) Cantones ubicados en el área de estudio.
1.5.3. Geomorfología
La geomorfología de la zona es el resultado de la larga evolución producida
principalmente por el tectonismo, plutonismo y procesos erosivos; factores que han modelado
la región hasta alcanzar el actual paisaje morfo-estructural.
El área de estudio forma parte del sistema montañoso austral interandino del Ecuador y
presenta relieves de tres categorías (UTPL, 2012), el principal es de tipo tectónico erosivo que
corresponde a levantamientos tectónicos que generan geoformas montañosas y colinadas de
diversa altura y pendiente; se destacan relieves colinados bajos hasta los de tipo montañoso
con alturas que oscilan desde 1200msnm al oeste en el sector de La Victoria, hasta 2660msnm
al este, en La Rinconada y en el Cerro Huachingue (Figura 2).
Relieves de tipo denudativo corresponden a depósitos coluviales y aluviales antiguos y

recientes (Figura 3), ubicados en el sector de Sabiango; estas geoformas han sido originadas
por la acción de la gravedad y por procesos erosivos.
Finalmente, se identifican relieves bajos de tipo deposicional, referidos a geoformas

originadas por la depositación de material acarreado por agentes erosivos como el agua y el
5
viento desde las partes altas hacia zonas bajas, cómo en la confluencia de la quebrada Guatara
y el Río Macará, cuya cota más baja está a 360 msnm, en este sector se encuentran amplias
terrazas aluviales.
Figura 2: Modelo digital elevaciones del área de estudio.
Figura 3: Mapa de unidades geomorfológicas del área de estudio. Modificado de INGEMMET (2006).
1.5.4. Hidrografía
El área de estudio se ubica en la parte norte de la Cuenca del Río Macará-Río Calvas, la
cual tiene una forma alargada con una orientación E-W y una superficie de 2833km2,
incluyendo la parte sur de la cuenca que se encuentra en territorio peruano. Esta cuenca a su
vez forma parte de la gran cuenca binacional de Catamayo-Chira. Figura 4
La Cuenca del Río Macará-Río Calvas limita al norte con la cuenca del Río Catamayo, al este
con las cuencas de los Ríos Zamora y Chinchipe (Rio Mayo), al sur y oeste en territorio
peruano con las cuencas de los Ríos Quiroz y Chira respectivamente.
6
Figura 4: Mapa de las cuencas hidrográficas que constituyen la cuenca binacional de Chira-Catamayo.
Los principales sistemas hidrográficos que atraviesan la zona de estudio están conformados
por los Ríos Pindo al este, Sabiango al oeste, Macará y Calvas al sur. Figura 5
Figura 5. Mapa de sub-cuencas y red hidrográfica del área de estudio
Los sistemas de drenajes del área están controlados tanto por la litología como por los
elementos estructurales presentes, de este modo definen un sistema denso de patrones
dendríticos a subparalelos, que generalmente corresponden a rocas ígneas y sedimentarias.
7
2. CONTEXTO GEOLÓGICO
2.1. Marco Geodinámico del Ecuador y Norte de Perú

El margen continental ecuatoriano está afectado por un límite convergente caracterizado
por la subducción de la Placa de Nazca bajo la Placa Sudamericana con dirección W-E a
razón de 70mm/año (Gutscher et al., 1999). La compresión E-W a la que está sometida la
Placa Sudamericana, ha provocado el desmembramiento de la porción noroeste de esta placa,
a la cual se ha denominado Bloque Norandino, una mini-placa independiente que se desplaza
con movimiento relativo en dirección N55ºE a razón de 10mm/año (Kellogs & Bonini, 1982)
a largo del sistema de fallas trancurrentes dextrales Guayaquil-Caracas. Figura 6
La placa de Nazca puede ser divida en dos segmentos separados por la Fractura de Grijalva
localizada entre 3ºS y 5ºS. Un segmento norte más joven de edad mioceno que se subduce con
un ángulo de 35º y un segmento sur de edad eoceno-oligoceno temprano que se subduce con
ángulos bajos (Wipf, 2006) a lo cual se le atribuye la usencia de volcanismo en esta zona.
Deflexión de
Huancabamba
Figura 6: Marco Geodinámico actual. Interacción de las placas Nazca, Cocos y Sudamericana. GG= Golfo de
Guayaquil, FZ= Zona de Fractura Grijalva, GCM= Mega Falla Guayaquil-Caracas. Gutscher et al. (1999)
La región fronteriza entre Ecuador y Perú se localiza dentro de la zona denominada

“Deflexión de Huancabamba” o “Andes de Huancabamba” que abarca el área comprendida
entre las latitudes 3°S en Ecuador hasta 8°S en Perú (Mourier et al., 1988). Esta estructura
8
representa la zona de transición entre los Andes Centrales con una orientación NNW-SSE y
los Andes del Norte de orientación NNE-SSW.
2.2. Marco Geológico Regional

La geología de la región sur del Ecuador y norte de Perú ha sido dividida en Dominios
Litotectónicos, cada uno de los cuales abarca unidades y formaciones con características
litológicas, estructurales y genéticas afines (Pilatasig et al., 2005). El área planteada para la
presente investigación se encuentra dentro de los Dominios Litotectónicos Loja-Olmos y
Cuenca Lancones-Alamor. Figura 7
Figura 7: Dominios Litotectónicos del sur de Ecuador y norte de Perú.

En rojo se enmarca el área de estudio. Tomado de Pilatasig et al. (2005)
2.2.1. Dominio Loja-Olmos

Forma parte de la zona oriental del área de estudio, está constituida de rocas y
estructuras emplazadas dentro de la Cordillera Real y su continuación en territorio peruano
denominado Complejo Metamórfico de Olmos. Comprende mayormente franjas sub-paralelas
de rocas metamórficas de bajo a alto grado con un rango de edad que va de paleozóico a
jurásico, cortados por cuerpos intrusivos tipo S del triásico y tipo I del jurásico. En el sur, las
rocas metamórficas se encuentran cubiertas por depósitos cenozoicos de origen volcánico
asociados a un arco volcánico continental y depósitos de cuencas continentales.
2.2.2. Cuenca Lancones-Alamor

Forma parte de la zona occidental del área de estudio. También denominada Cuenca
Celica-Lancones o Alamor-Punta de Piedra (INIGEMM, 2013), se ubica en el extremo
suroeste de Ecuador, entre el Macizo de Amotape-Tahuín y el Terreno Loja-Olmos.
9
La cuenca se habría formado durante el cretácico inferior, producto de un adelgazamiento de
la corteza continental y presenta dos regiones litológicamente bien representadas. La Región
Oriental de origen volcánico con escasas intercalaciones sedimentarias y la Occidental de
origen netamente sedimentario (Figura 8); el contacto entre las dos regiones se interpreta
como transicional o interdigitado. Hacia el oeste las secuencias sedimentarias de la cuenca
yacen discordantes sobre las rocas metamórficas del Macizo Amotape-Tahuín y sobre rocas
volcánicas hacia la frontera con Perú (Egüez y Poma, 2001). En el este, el contacto entre las
secuencias de la cuenca y el Terreno Loja-Olmos es tectónico (Pilatasig et al., 2005).
Figura 8: Columna estratigráfica de la Cuenca Lancones-Alamor. Incluye las Unidades litológicas presentes en
las Regiones Oriental y Occidental de la Cuenca. Tomado de Pilatasig et al. (2005)
2.3. Litoestratigrafía del área de estudio

La zona de Macará está constituida por rocas volcánicas y volcanosedimentarias de las
Unidades Pilares y Punta de Piedra, relacionadas a un arco volcánico oceánico del cretácico; y
rocas volcánicas y volcanoclásticas asociadas a un arco volcánico continental del paleógeno.
Además de rocas intrusivas de composición intermedia a ácida emplazadas en el cretácico
superior representadas por el Batolito de Tangula y stocks subvolcánicos recientes. Figura 9
2.3.1. Unidad Los Pilares

Esta unidad es una extensión del miembro superior de la Formación San Pedro del Perú,
aflora al sureste de Macará, sector el Tabacal, donde forma un anticlinal con eje NE/SW y
flancos con buzamientos fuertes. Constituye una secuencia sedimentaria de ambiente marino,
10
compuesta por intercalaciones de areniscas de grano fino, lutitas negras y niveles delgados de
cherts con microplegamientos y bandeamientos (Pilatasig et al., 2005). Un fósil de
Inoceramus sp. de la Formación San Pedro en la Quebrada Pilares (frontera Ecuador-Perú,
fuera del área) indica una edad cretácica inferior para esta unidad (INGEMMET, 1987).
2.3.2. Unidad Punta de Piedra

Kennerley (1973) define inicialmente a las rocas volcánicas del sur del país como
Formación Piñón; posteriormente son denominadas Formación Celica por Lebrat et al.
(1979). Egüez y Poma (2001) redefinen a las rocas volcánicas aflorantes al este y noreste de
Macará como Formación Punta de Piedra. Esta Unidad predomina en la parte central,
cubriendo aproximadamente un 50% del mapa geológico de Macará. Observaciones de
campo, delimitan una potencia entre 1500 y 2000m (INIGEMM, 2016a). La litología típica
incluye gabros, basaltos, andesitas, andesitas basálticas y pilow lavas asociadas a un arco
volcánico toleítico formado sobre corteza continental adelgazada durante el cretácico inferior
tardío, Aptiense-Albiense (Pilatasig et al., 2005). INIGEMM (2013) interpreta a las rocas de
esta unidad como procedentes de un arco de islas.
Los gabros ocurren como pequeños cuerpos aflorantes al este de Macará constituidos
principalmente por piroxeno (80%); hornblenda, clorita y epidota se presentan como
minerales de alteración propilítica incipiente (Pilatasig et al., 2005). Los basaltos ocurren
principalmente en el sector del Río Sabiango y en el Cerro Punta de Piedra. En los bordes del
Batolito, las rocas basálticas presentan biotita secundaria y cuarzo producto de metamorfismo
de contacto.
Metavolcánicos débilmente foliados afloran hacía el este de la zona, mismos que

posiblemente sean ventanas de la Unidad Alao-Paute de edad Jurásica, indicando un posible
basamento metamórfico en esta parte de la cuenca (Pilatasig et al., 2005). La unidad también
presenta una componente volcanoclástica formada por hialoclastitas, tobas y brechas; y
volcanosedimentos representados por areniscas y limolitas.
2.3.3. Formación Sacapalca

Esta formación definida inicialmente por Kennerley (1973), agrupa las secuencias
volcánicas y volcanosedimentarias de afinidad calco-alcalina que conforman una franja
alargada ubicada al este del área, entre la Cuenca Lancones-Alamor y los metamórficos del
Terreno Loja-Olmos. A esta formación se le atribuye una edad paleoceno-eoceno temprano
(66.9±5.8Ma Hungerbuhler, 1997; 44-55Ma Jaillard et al., 1996).
11
Pilatasig et al. (2007) discriminan facies de tobas riolíticas, tobas purpuras y lavas andesito-
basálticas. INIGEMM (2013), separa a la Formación Sacapalca en las Unidades Changaimina,
que representa una facie lávica andesítica; Unidad La Cruz constituida por tobas andesíticas,
aglomerados y brechas; Unidad Río Guayucu representado por tobas purpuras y
ocasionalmente lentes de arenisca; y la Unidad Ahuaca formada por tobas ácidas y stocks
dacíticos.
2.3.4. Rocas Intrusivas
2.3.4.1. Batolito de Tangula

Complejo plutónico de gran dimensión que intruye a las secuencias aflorantes en la región
oriental de la Cuenca Lancones-Alamor. Regionalmente tiene forma de herradura abierta
hacia el sur con prolongaciones irregulares hacia territorio peruano. En el este se encuentra
sobreyacido por las rocas de la Fm. Sacapalca, sin embargo, en el sector de Bella María
(Fuera del área de estudio) el contacto es tectónico (INIGEMM, 2013).
El batolito presenta una afinidad calco-alcalina emplazado en un arco volcánico continental

(INIGEMM, 2013) y está conformado por al menos 3 fases de intrusión caracterizadas por el
enriquecimiento de feldespato potásico, hornblenda y cuarzo (Pilatasig et al., 2005). La
litología típica incluye dioritas, granodioritas ricas en hornblenda y granitos. Se reportan dos
edades para el batolito de 111±3Ma (Kennerley et al., 1973) y 93±1.1Ma (Ullrich, 2004).
2.3.4.2.Complejo Peribatolítico
Un Complejo Peribatolítico, constituido por stocks pequeños de diorita y filones pórfido-
andesíticos rodea la parte norte del batolito (Egüez y Poma, 2001); posiblemente representa a
la fase final de emplazamiento del batolito y estaría relacionado con mineralización en la zona
(INIGEMM, 2013). De manera similar a lo reportado en Perú donde un borde gabro-diorítico
rodea a cuerpos intrusivos ácidos e intermedios emplazados en la parte central del Complejo
Plutónico de las Lomas correlacionable con el Batolito de Tangula (INGEMMET, 2012).
2.3.4.3. Pórfidos
Varios cuerpos porfiríticos de composición andesítica, dacítica y riolítica están presentes
en el área. Subvolcánicos andesíticos ocurren como diques y cuerpos circulares que cortan las
secuencias volcanoclásticas cenozoicas hacia el este (INIGEMM, 2016a).
12
Figura 9: Mapa Geológico de Macará Escala 1:100 000. Modificado de Pilatasig et al. (2007)
13
2.4. Geología Estructural
Las fallas principales presentan una dirección preferencial NE-SW; mientras que el
fracturamiento presenta dos direcciones, una paralela a las estructuras mayores y otra
perpendicular. El fallamiento de las unidades de la cuenca ha favorecido el emplazamiento de
diques, vetas y vetillas según las direcciones NNE-SSO, NNW-SSE y E-O (INIGEMM,
2016a).
En el sector Bella María (fuera del área), una estructura de tipo inversa con buzamiento hacia
el oeste cabalga las rocas intrusivas de Tangula sobre las secuencias de Sacapalca. Pliegues
anticlinales con eje NNE-SSW y fuertes buzamientos afectan a la U. Los Pilares y Fm.
Sacapalca lo que representaría el ambiente compresivo de la cuenca (Pilatasig et al., 2007).
La tectónica transpresiva, se refleja en lineamientos y estructuras regionales de orientación E-

W que ocurren hacia el norte de la cuenca, en el BAT y hacia el sur de la cuenca en territorio
peruano; cinemática que se manifiesta en el área de estudio en estructuras antitéticas y
ortogonales entre sí (INIGEMM, 2016a). Figura 10
Figura 10: Mapa de estructuras, frontera sur del Ecuador. Modificado de Pilatasig et al. (2007)
2.5. Geología Económica

2.5.1. Concesiones Mineras
El catastro minero de la ARCOM indica que en la zona estudio existen concesiones
inscritas y en trámite para minería metálica, principalmente para exploración de yacimientos
14
tipo primario de Au, Cu y Ag; y otras para oro aluvial en especial en la cuenca del Rio
Sabiango; lo que muestra el interés que existe por parte de las compañías de realizar
exploración en esta zona. Figura 11
Figura 11: Ubicación de concesiones mineras. Catastro minero actualizado a agosto de 2017 (ARCOM)
2.5.2. Metalogenia del Sur del Ecuador

Paladines (1980) incluye las mineralizaciones presentes en los sectores de Los Linderos,
Rio Playas, El Huato, y Macará (U) dentro de la denominada Sub-Provincia Polimetálica del
Altiplano, que incluye principalmente depósitos tipo pórfido cuprífero que se habrían
desarrollado del cretácico superior al mioceno. INIGEMM (2013) relaciona la mineralización
presente en los sectores de Rio Playas y El Huato a la última fase de intrusión del Batolito de
Tangula caracterizado por un complejo intrusivo-filoniano diorítico y pórfido-andesítico que
rodea la parte norte del Batolito.
En el Mapa Tectonometalogénico del Ecuador del año 1993, se reconocen 8 zonas

mineralizadas alrededor de Macará, principalmente de tipo Pórfido (Cu±Mo±Au),
Epitermales (Au+Ag+Zn+Pb±Sb±Mn), Vetas (Au) y stockworks posiblemente relacionados
al magmatismo del Batolito de Tangula (Tabla 1, figura 12). En el sector también se
reconocen trabajos de oro aluvial en el Rio Sabiango y Quebrada Laguar (Merlyn y
Mosquera, 1988). PRODEMINCA (2000) relaciona las mineralizaciones del sur del país con
el Evento Geotectónico Celica del cretácico superior (115-85Ma) caracterizado por
volcanismo y plutonismo de cuenca pull-apart; y plutonismo y sedimentación continental.
Por su parte, el carácter metalogénico del norte de Perú ha sido bien definido. Allí se
identifican diversos tipos de depósitos metálicos que incluyen VMS de Pb-Zn-Cu hospedados
15
en secuencias volcánicas y volcanosedimentarias del albiense, pórfidos Cu-Mo hospedados en
granitoides del cretácico superior, vetas epitermales de Au-Ag hospedadas en secuencias
volcánicas del albiense-cenomaniense y skarns de magnetita (Rodríguez et al., 2008).
Nombre X Y País Asociación Edad de Tipo de depósito

mineralización
Au-
Angascola 631443 9511494 Ecuador Cretácico Relacionado a intrusivo
Polimetálico
Zapotillo 582171 9514737 Ecuador Au Cretácico Relacionado a intrusivo
Río Playas 633461 9546752 Ecuador Cu-Mo Cretácico Relacionado a intrusivo
El Huato 645896 9543047 Ecuador Cu-Mo Cretácico Relacionado a intrusivo
Los Linderos 599131 9527303 Ecuador Cu-Mo Cretácico Relacionado a intrusivo
Changaimina 666745 9531340 Ecuador Polimetálico Cenozoico Epitermal
Laguar 608689 9522309 Ecuador Au Cretácico Relacionado a intrusivo
Cariamanga 657273 9520068 Ecuador Polimetálico Cenozoico Epitermal
Sabiango 631873 9515760 Ecuador Au-Cu Cretácico Relacionado a intrusivo
s/n 629339 9513924 Ecuador Au - Oro Aluvial
Los Colorados 667857 9514353 Perú Polimetálico Cenozoico Epitermal
Pampas
612969 9514927 Perú Au Cretácico sup-Paleoceno Epitermal
Redondas
El Suyo 607214 9506725 Perú Au Cretácico sup-Paleoceno Epitermal
Servilleta 603775 9507585 Perú Au Cretácico sup-Paleoceno Epitermal
Rinconada 633315 9483230 Perú Cu-Mo Cretácico superior Pórfido Cu-Mo
Lanchipamba 659880 9489792 Perú Cu Cretácico superior Veta
Lagunas 631060 9471588 Perú Au Cretácico superior Epitermal
Bolsa del
593019 9505407 Perú Au Cretácico sup-Paleoceno Epitermal
Diablo
El Noque 581510 9494956 Perú - Cretácico sup-Paleoceno Skarn
Alumbre 575358 9499850 Perú Au Cretácico sup-Paleoceno Epitermal
Petrobayo 575028 9493765 Perú Cu Cretácico sup-Paleoceno VMS Pb-Zn-Cu
Pilares 567685 9490325 Perú Au Cretácico sup-Paleoceno Epitermal
La Bocana 573043 9487680 Perú Zn Cretácico superior Pórfido Cu-Mo
Pampa
571323 9486952 Perú - Cretácico superior Skarn
Quemadas
La Saucha 599303 9485629 Perú Pb Cretácico sup-Paleoceno VMS Pb-Zn-Cu
Lagartos 566561 9479411 Perú Cu-Mo Cretácico superior Pórfido Cu-Mo
Chancadora 556374 9473987 Perú Cu-Mo Cretácico superior Pórfido Cu-Mo
Somate 562724 9470614 Perú Cu Cretácico sup-Paleoceno VMS Pb-Zn-Cu
Tambogrande 572646 9454276 Perú Cu-Zn Cretácico inferior VMS Cu-Zn-Au
Cerro
590770 9460957 Perú Cu Cretácico sup-Paleoceno VMS Pb-Zn-Cu
Colorado
El Papayo 593019 9467836 Perú Pb Cretácico sup-Paleoceno VMS Pb-Zn-Cu
Orquetas 578953 9485804 Perú Cu-Mo Cretácico superior Pórfido Cu-Mo
La Copa 604499 9491688 Perú Pb Cretácico sup-Paleoceno VMS Pb-Zn-Cu
Tabla 1. Minas y prospectos metálicos del sur del Ecuador y norte de Perú.
Fuente: CODIGEM (1993), Pilatasig et al. (2007), INGEMMET (2012)
16
GEOLOGÍA REGIONAL Y OCURRENCIAS MINERALES
DEL SUR DEL ECUADOR Y NORTE DE PERÚ
Figura 12: Geología Regional y Ocurrencias Minerales del sur del Ecuador y norte de Perú. Compilado de: CODIGEM (1993), Egüez y Poma (2001), Pilatasig et al. (2007), INGEMMET (2012), INIGEMM (2013)
17
3. MARCO TEÓRICO
3.1. Exploración Geoquímica

Método de investigación minero basado en la distribución espacial, abundancia,
asociaciones y propiedades químicas de los elementos en la naturaleza. Tiene como objetivo
la localización de anomalías geoquímicas posiblemente asociadas a cuerpos mineralizados.
Las campañas geoquímicas se clasifican como de reconocimiento y de detalle, en función del

propósito y la escala; también se pueden clasificar de acuerdo al material muestreado (suelo,
vegetación, sedimentos, agua, roca y gas). El planeamiento óptimo de una campaña
geoquímica debe basarse en el conocimiento de la distribución, abundancia y principios que
gobiernan la migración de los elementos (Rose, Hawkes & Webb, 1979).
3.2. Parámetros Geoquímicos

El background o valor de fondo geoquímico se refiere a la concentración natural de un
elemento. Puede ser determinado para cada elemento, para cada área y para cada tipo de
material muestreado.
Una anomalía geoquímica representa la abundancia por encima de lo normal de un elemento

en una asociación o ambiente que no le corresponde; geológicamente aparece cuando se
adiciona algún componente a las rocas a partir de una fuente externa (García y Moreno,
2006), esto ocasiona enriquecimiento de elementos y contrastes con las concentraciones
originales de la litología. Estadísticamente constituye los valores iguales o más altos al
threshold o umbral anomálico, mismo que está en función de la concentración normal más
dos veces la desviación estándar (Lepeltier, 1969).
3.3. Elementos Guía

En exploración geoquímica, un elemento medido con el objetivo de detectar un cuerpo
mineral es nombrado elemento indicador. En la mayoría de situaciones, el elemento indicador
es un componente económicamente valioso del mineral buscado (target element), por
ejemplo, Cu para una mena de Cu. Sin embargo, si el componente valioso es difícil de
analizar, es inmóvil o arroja resultados difíciles de interpretar; otros elementos asociados al
mineral o a un tipo de mineralización pueden ser útiles; tales elementos se denominan
18
pathfinder o elemento explorador, que pueden ser detectados más fácilmente, como es el caso
del As como pathfinder para depósitos de Au. Tabla 2
Tipo de depósito mineral Componentes mayores Elementos asociados

Depósitos Magmáticos
Mena de cromita Cr Ni, Fe, Mg
Capas de magnetita Fe V, Ti, P
Sulfuros de Cu-Ni inmiscibles Cu, Ni, S Pt, Co, As, Au
Intrusión en capas de Pt-Ni-Cu Pt, Ni, Cu Cr, Co, S
Óxidos de Fe-Ti inmiscibles Fe, Ti P
Carbonatita de Nb-Ta Nb, Ta Na, Zr, P
Pegmatitas de metales raros Be, Li, Cs, Rb B, U, Th, Tierras raras
Depósitos Hidrotermales
Pórfido de cobre Cu, S Mo, Au, Ag, Re, As, Pb, Zn, K
Pórfido de molibdeno Mo, S W, Sn F, Cu
Skarn de magnetita Fe Cu, Co, S
Skarn de Cu Cu, Fe, S Au, Ag
Skarn de Pb-Zn Pb, Zn, S Cu, Co
Skarn de W-Mo-Sn W, Mo, Sn F, S, Cu, Be, Bi
Vetas de metales base Pb, Zn, Cu, S Ag, Au, As, Sb, Mn
Greisen de Sn-W Sn, W Cu, Mo, Bi, Li, Rb, Si, Cs, Re, F, B
Vetas de sulfuro de estaño Sn, S Cu, Pb, Zn, Ag, Sb
Vetas de Co-Ni-Ag Co, Ni, Ag, S As, Sb, Bi, U
Epitermal de metales preciosos Au, Ag Sb, As, Hg, Te, Se, S, U
Depósito de mercurio Hg, S Sb, As
Veta de uranio U Mo, Pb, F
Cobre en basaltos Cu Ag, As, S
VMS de Cu Cu, S Zn, Au
VMS Zn-Cu-Pb Zn, Pb, Cu, S Ag, Ba, Au, As
Formaciones de Fe ricas en Au-As Au, As, S Sb
Mississippi Valley Pb-Zn Zn, Pb, s Ba, F, Cd, Cu, Ni, Co, Hg
Mississippi Valley de fluorita F Ba, Pb, Zn
Arenisca urinífera U Se, Mo, V, Cu, Pb
Capas rojas de Cu Cu, S Ag, Pb
Calcrete de U U V
Tipo Sedimentario
Lutita de Cu Cu, S Ag, Zn, Pb, Co, Ni, Cd, Hg
Arenisca cuprífera Cu, S Ag, Co, Ni
Tabla 2. Asociaciones de elementos y pathfinders útiles en la identificación de yacimientos minerales.
Traducido de Rose, Hawkes & Webb, (1979)
3.4. Tipos de mineralización

Basado en la configuración geológica y tectónica del área de estudio, se puede
establecer que presenta afinidad para la ocurrencia de depósitos epitermales, pórfidos y VMS.
3.4.1. Epitermales
Los depósitos epitermales se forman a profundidades menores a 1km (Corbett, 2013)
por circulación de fluidos hidrotermales con temperaturas entre 100 y 300ºC que depositan
minerales producto de complejas reacciones químicas. Por lo general ocurren en arcos
volcánicos relacionados a márgenes de subducción (arcos continentales y arcos de islas)
19
aunque también ocurren asociados a cuencas pull-apart. Pueden ser clasificados en Alta
Sulfuración (HS) y Baja Sulfuración (LS) en base a la naturaleza y distribución de la
alteración hidrotermal, pH y temperatura del fluido hidrotermal, forma del depósito y texturas
de minerales de mena y ganga (PRODEMINCA, 2000a). Los depósitos epitermales de alta
sulfuración se derivan de fluidos ácidos ricos en sulfuro generados en ambientes
hidrotermales volcánicos; a diferencia de los depósitos epitermales de baja sulfuración
caracterizados por fluidos reducidos casi neutros empobrecidos en sulfuro producto de una
mayor interacción con agua de origen meteórico. En la tabla 3, se resume las principales
características de estos tipos de depósito.
Baja sulfuración Alta Sulfuración

Margen continental volcano-plutónico y
Ambiente Margen continental volcano-plutónico y arcos
arcos oceánicos y tras-arcos. Campos
tectónico oceánicos y tras-arcos.
volcánicos continentales
Ambiente Sistemas de fallas regionales, grabens, Sistemas de fallas regionales, grabens,
geológico / calderas silícicas, estratovolcanes calderas silícicas, estratovolcanes andesíticos.
estructural andesíticos Pueden sobreyacer a sistemas porfídicos
Forma del Remplazamiento masivo finogranular y

Vetas y stockworks dominantes.
depósito diseminados, lentes y masas irregulares.
(proximal) silicificación, argilítica avanzada,
Zonación de (proximal) silicificación, sericítica
argilítica,
alteración argilítica, propilítica (distal).
propilítica (distal).
Texturas de Cuarzo coloforme, crustiforme, texturas
Reemplazamientos de sílice masiva y residual,
acrestadas, bladed, reemplazo de carbonato,
mena sílice geodal, buggy
cuarzo drusi.
Pirita, enargita-luzonita, bornita, covellina,
Pirita, oro electrum, oro, plata, argentita. oro, oro electrum.
Minerales de Subordinados: calcopirita, galena, esfalerita,
Subordinados: esfalerita, galena, calcopirita,
mena tetraedrita. tennantita, marcasita, arsenopirita,
tetraedrita, piragirita, selenuros.
sulfosales de plata, telururos.
Minerales de Cuarzo, calcedonia, amatista, carbonato
Cuarzo, barita, yeso, anhidrita
ganga (Mn), adularía, sericita, barita, fluorita.
Minerales de Alunita, caolinita, dickita, pirofilita, jarosita,
Sericita, illita, adularía, clorita, esmectita.
alteración esmectita.
Au, Ag. Locales o menores Zn, Pb, Cu, Cu, Au, As. Locales o menores Ag, Pb, Bi,
Metales
Mo, Sb, As, Te, Se, Hg, Ba, F, Mn. Sb, Mo, Sn, Zn, Te, Hg, W, B.
Zonación de Superior: As, Sb, Hg, B, TI, Au, Ag. Superior: As, Sb, Hg, B, TI, Au, Ag.
metales Inferior: Cu, Pb, Zn, Bi, Te, Se, Co. Inferior: Cu, Pb, Zn, Bi, Te, Se. Co.
Tabla 3. Principales características de los depósitos epitermales HS y LS. Tomado de PRODEMINCA., (2000a)
Estos depósitos en general presentan una zonación hacia la parte superficial de los elementos
As, Sb, Hg, B, TI, Au, Ag mientras que Cu, Pb, Zn, Bi, Te, Se, Co están mayormente
presentes hacia la parte interna del sistema. Figura 13
20
Figura 13: Modelos de mineralización de para epitermales Au-Cu de la costa del Pacífico.
Tomado de Cobertt & Leach, (1997)
3.4.2. Pórfidos
Los depósitos porfídicos son sistemas de gran tonelaje, generalmente de mineralización
dispersa y baja ley; se asocian a rocas intrusivas intermedias a acidas comúnmente de textura
porfídica. Se caracterizan por eventos múltiples de intrusión, enjambres de diques y brechas
intrusivas (PRODEMINCA, 2000b).
En cuanto a la tectónica están asociados a márgenes de subducción de volcanismo activo. Los

pórfidos se clasifican según el elemento principal que contengan, es así que se tiene: Pórfidos
de Cu-Au asociados a arcos de islas, pórfidos de Cu-Mo asociados a márgenes continentales y
series magmáticas calco-alcalinas diferenciadas. Figura 14
La mineralización y alteración forma amplias zonas con cambios laterales y una zonación
bien definida iniciando con alteración potásica hacia el núcleo de mayor temperatura del
sistema (biotita, feldespato potásico), alteración clorítica-sericítica (clorita, sericita, illita,
hematita) y fílica hacia la parte superior (cuarzo, sericita); argílica avanzada hacia la parte
externa (cuarzo, alunita, pirofilita, dickita, caolinita) y propilítica (clorita, epidota, calcita) en
las partes más distales y de menor temperatura (Sillitoe, 2010).
21
Figura 14: Modelo de mineralización y alteración de un depósito tipo pórfido. Tomado de McQueen, (s/f)
3.4.3. Sulfuros Masivos Volcanogénicos (VMS)

Los Sulfuros masivos Volcanogénicos (VMS) son depósitos minerales constituidos por
masas de diversos tamaños de sulfuros polimetálicos (60% de la masa mineral) dispuestos de
manera estratiforme o lenticular entre niveles de rocas volcánicas y volcanosediementarias
depositadas originalmente en fondo oceánico; también pueden presentar mineralización en
stockwork a lo largo del sistema de alimentación del VMS (PRODEMINCA, 2000c). Estos
depósitos tienen un origen volcánico exhalativo y son característicos de zonas de tectónica
extensional el cual permite la emisión de magma desde el manto y el ascenso de fluidos
hidrotermales asociados: dorsales oceánicas, zonas extensionales intracontinentales, arcos de
islas especialmente en cuencas tras-arco y rifts continentales.
Los VMS son las mayores fuentes de Zn, Cu, Pb, Ag, Au y fuentes significativas de Co, Sn,
Se, Mn, Cd, In, Bi, Te, Ga y Ge (Galley, 2007). Existen varias clasificaciones para estos
depósitos, la más aceptada plantea tres tipos principales: Tipo Kuroko (Cu-Zn±Au±Ag)
asociado a volcanismo bimodal de lavas toleíticas y secuencias félsicas y calco-alcalinas
situadas en zonas oceánicas de convergencia de placas, Tipo Chipre (Cu±Zn±Au) asociado a
secuencias basálticas de complejos ofiolíticos generados en zonas de creación de corteza,
Tipo Besshi (Cu-Zn±Au±Ag) relacionado a sedimentos clásticos y rocas volcánicas básicas
de intraplaca (Sawkins, 1976 en PRODEMINCA, 2000c).
22
El modelo ideal de mineralización de un VMS presenta una zonación bien definida
constituida por un núcleo central enriquecido en pirita (Fe) y calcopirita (Cu) rodeado por una
zona rica en galena (Pb) y esfalerita (Zn). En la zona periférica se encuentran zonas ricas en
esfalerita, galena y baritina (Ba). Figura 15
Figura 15: Modelo de mineralización, alteración y zonación de elementos para VMS.

Tomado de McQueen, (s/f)
3.5. Geoestadística
Conjunto de técnicas estadísticas utilizadas para el análisis de datos o fenómenos que
varían espacialmente. Por lo general, la variación espacial de estos fenómenos es compleja
para ser descrita por un modelo matemático sencillo y aunque parezcan completamente
aleatorios, tienen cierto grado de correlación espacial.
La primera etapa que se debe cumplir en Geoestadística es el Análisis Exploratorio de Datos

Espaciales (AEDE) para posteriormente identificar de manera experimental la estructura de
autocorrelación espacial y direccional del conjunto de datos mediante el semivariograma y su
uso en la predicción a través de técnicas de interpolación.
El AEDE se fundamenta en una estadística descriptiva apoyado del uso de gráficas de

frecuencia y análisis de parámetros estadísticos como: medidas de tendencia central y no
central, medidas de dispersión y medidas de forma.
23
3.5.1. Medidas de tendencia central
Intentan identificar el dato más representativo de la distribución del conjunto y
corresponden a la media, mediana y moda.
• Media (μ o X). Es el resultado de la suma de todos los valores divididos para el
número de datos. Es influenciada por todos los valores que asume la variable y es
sensible a los valores extremos. (Fórmula 1.1; Ver Anexo 1)
• Mediana (Me). Es el valor que divide al conjunto de datos en dos partes iguales. No
es sensible a los valores extremos debido a que depende únicamente del orden de los
datos. (Fórmula 1.2)
• Moda (Mo). Es el valor que presenta mayor frecuencia absoluta, es decir que se repite
más veces en la distribución de datos y puede no ser única. (Fórmula 1.3)
3.5.2. Medidas de tendencia no central

• Cuartiles. Son tres valores que distribuyen la serie de datos en cuatro tramos iguales
correspondientes al 25% de los resultados.
• Deciles. Son nueve valores que distribuyen la serie de datos en diez tramos iguales
correspondientes al 10% de los resultados.
• Percentiles. Son noventa y nueve valores que distribuyen la serie de datos en cien
tramos iguales correspondientes al 1% de los resultados.
3.5.3. Medidas de dispersión

Indican la mayor o menor concentración de los datos con respecto a las medidas de
tendencia central. Dan una idea sobre la homogeneidad o cuan agrupados están los datos.
• Desviación estándar (S). Mide la dispersión de los valores con respecto a la media.
Es fuertemente influenciada por los valores extremos. Una desviación estándar alta
indica que los datos están lejos de la media, y una desviación baja indica que los datos
están agrupados cerca de la media. (Fórmula 1.4)
• Varianza (S2). Describe la variabilidad de la distribución. Mide la distancia entre los
valores de la serie y la media. Mientras más se aproxima a cero, más concentrados
están los datos entorno a la media; y mientras mayor sea la varianza, los datos se
encuentran más dispersos. (Fórmula 1.5)
24
• Coeficiente de variación (CV). Mide la representatividad de la media. Valores
extremos concluyen que la media no es representativa, es decir, existen valores entre
las observaciones que se separan de las demás. (Fórmula 1.6)
3.5.4. Medidas de Forma

Evalúan la forma de la distribución de los datos respecto a una curva modelo de
distribución normal.
• Coeficiente de Curtosis (K). Mide cuán apuntada o achatada es la curva con respecto
a una distribución normal. Si los datos están muy concentrados hacia la media es
leptocúrtica (curtosis >0). Si los datos están muy dispersos, la distribución es
platicúrtica (curtosis <0). La distribución normal tiene una distribución mesocúrtica y
su curtosis es igual a 0. (Fórmula 1.7)
• Coeficiente de Sesgo o Asimetría (CS). Es la forma que toma la curva en el eje de las
abscisas, evalúa el grado de inclinación de la curva. (Fórmula 1.8)
3.5.5. Normalidad de la distribución de datos

Un análisis geoestadístico será óptimo, únicamente si la variable en estudio sigue una
distribución normal de datos (Campana de Gauss), respecto a los siguientes parámetros:
• Respecto a la media, mediana y moda. Para que una distribución sea normal; su
media, mediana y moda deben ser valores iguales o similares.
• Respecto al coeficiente de sesgo. Si el CS es menor a 0.5 se acepta como una
distribución normal y se puede aplicar el método geoestadístico. Si el CS está entre 0.5
y 1, es necesaria una transformación de datos de tipo raíz cuadrada. Si Cs es mayor a
1, es necesaria una transformación de datos de tipo logarítmico. Realizado la
transformación, es necesario calcular nuevamente los parámetros estadísticos y
realizar la prueba de normalidad.
• Respecto al coeficiente de variación. Si CV es menor a 100, no hay problema de
datos extremos. Si CV esta entre 10 y 200, los efectos causados por los valores
extremos son tolerables. Si CV es superior a 200, los datos extremos influyen
fuertemente.
3.5.6. Variograma y Semivariograma

Es una función que mide el grado de diferencia de las variables en relación a la distancia
de separación de las mismas, es decir, la autocorrelación espacial de los datos (Moreno,
25
2008). Su función se denota como 2γ(h) en la que la varianza de los incrementos de la
variable regionalizada es infinita. Utilizando la definición teórica de la varianza en términos
del valor esperado de una variable aleatoria, se tiene:
Fórmula 1.9
La mitad del variograma γ(h), se conoce como semivarianza y caracteriza las propiedades de
dependencia espacial. Se estima por el método de momentos, a través del semivariograma
experimental, que se calcula mediante:
Fórmula 1.10
Donde Z(x) es el valor de la variable en un lugar x; Z(x+h) es otro valor de muestra separado
del anterior una distancia h y n es el número de parejas que se encuentran separadas por dicha
distancia.
3.5.7. Modelos Teóricos de Semivariograma

Es necesario ajustar el modelo de la semivarianza con el fin de que generalice lo
observado en el semivariograma experimental a cualquier distancia. Los modelos teóricos de
semivarianza se dividen en no acotados (lineal, logarítmico, potencial) y acotados (esférico,
exponencial, gaussiano) (Giraldo, 2008). Este segundo grupo garantiza que la covarianza de
los incrementos es finita, por lo cual son ampliamente usados cuando hay evidencia de que
presentan buen ajuste. Los modelos tienen tres parámetros comunes:
• Efecto Pepita (Nugget). Se denota C0 y representa una discontinuidad puntual del

semivariograma en el origen (Figura 16). Puede producirse por errores de medición en
la variable o a la escala de la misma. En ocasiones puede indicar que parte de la
estructura espacial se concentra en distancias inferiores a las observadas.
• Meseta (Sill). Se denota C1 ó (C0 + C1) cuando la pepita es diferente de cero. Es la
cota superior del semivariograma o el límite del semivariograma cuando la distancia h
tiende a infinito; la meseta puede ser o no finita.
• Rango. Se denota como a. Corresponde a la distancia a partir de la cual dos
observaciones son independientes. El rango se interpreta como la zona de influencia. Existen
26
modelos de semivariograma en los que no existe una distancia finita en la que dos
observaciones sean independientes; por ello se llama rango efectivo a la distancia para la cual
el semivariograma alcanza el 95% de la meseta. Entre más pequeño sea el rango, más cerca se
está del modelo de independencia espacial.
Figura 16: Comportamiento típico de un semivariograma y la representación de sus parámetros básicos.

Tomado de Giraldo, (2002)
3.5.7.1. Modelo Esférico

Tiene un crecimiento rápido cerca al origen, pero los incrementos marginales van
decreciendo para distancias grandes, hasta que para distancias superiores al rango los
incrementos son nulos. Su expresión matemática es la siguiente:
Fórmula 1.11
3.5.7.2. Modelo Exponencial

Este modelo se aplica cuando la dependencia espacial tiene un crecimiento exponencial
respecto a la distancia entre las observaciones. El valor del rango es igual a la distancia para la
cual el semivariograma toma un valor igual al 95% de la meseta. Este modelo es ampliamente
usado. Su expresión matemática es la siguiente:
Fórmula 1.12
27
3.5.7.3. Modelo Gaussiano
Al igual que en el modelo exponencial, la dependencia espacial se desvanece solo en
una distancia que tiende a infinito. El principal distintivo de este modelo es su forma
parabólica cerca al origen (Figura 17). Su expresión matemática es:
Fórmula 1.13
Figura 17: Comparación de los modelos Exponencial, Esférico y Gaussiano. Tomado de Giraldo, (2002)
3.5.8. Interpolación
Se define como un procedimiento que permite calcular el valor de una variable en una
posición del espacio (punto no muestral) conociendo los valores de esa variable en otras
posiciones del espacio (puntos muestrales) (Bosque, 1992 en Moreno, 2008). Los métodos de
interpolación se pueden clasificar en directos (deterministas) y analíticos (geoestadístico) en
función de mecanismo matemático de interpolación.
Los métodos analíticos tienen como objetivo generar superficies que incorporan las
propiedades estadísticas de los datos muestrales y que proporcionan una medida de error de
las mismas, siendo este último un indicador de una buena o mala predicción (Moreno, 2008).
Todos estos métodos geoestadísticos forman parte de la familia del Kriging.
3.5.9. Kriging
Encierra un conjunto de métodos de predicción espacial fundamentados en la
minimización del error cuadrático medio de predicción. En la tabla 4 se mencionan los tipos
de Kriging y algunas de sus propiedades.
28
Tipo de Nombre Propiedades
predictor
• Simple • Son óptimos si hay normalidad multivariada.
Lineal • Ordinario • Independiente de la distribución son los mejores
• Universal predictores linealmente insesgados.
• Indicador • Son predictores óptimos.

• Probabilístico
No lineal • Log Normal,
Trans-Gaussiano
• Disyuntivo
Tabla 4. Tipos de Predictores Kriging y sus propiedades. Tomado de Giraldo, (2002)
El Kriging asume que el fenómeno que se está estudiando toma valores en todas partes del
área de estudio, es decir, tiene continuidad espacial. El Kriging Ordinario o Puntual, asume
que la variable es estacionaria y que no tiene tendencia, se denota por la ecuación: 𝑍(𝑠) =
𝜇(𝑠) + 𝜀(𝑠); donde, Z(s) es la variable de interés; 𝜇(𝑠) una constante (media) desconocida,
𝜀(𝑠); los errores aleatorios; s simplemente indican el lugar de unas determinadas coordenadas
espaciales X y Y.
29
4. MARCO METODOLÓGICO
4.1. Muestreo
Durante el primer semestre del año 2015, el grupo de Geoquímica del INIGEMM
recolectó un total de 661 muestras de sedimentos fluviales a partir de drenajes de primer y
segundo orden que cubren las hojas topográficas 1:50 000 de Macará y Río Calvas, con una
densidad de muestreo de 0.9 muestras por km2.
Los equipos recorrieron las quebradas tomando muestras de sedimentos, en puntos sin
contaminación visible a una distancia mínima de 100 a 150m de las uniones entre drenajes. Se
recolectó aproximadamente 300g de sedimento tamizado en húmedo en el sitio con una malla
Nº80 (177µm) que fue guardado y codificado en bolsas de papel kraft no ácido, para ser
entregadas al laboratorio del Instituto.
4.2. Análisis de laboratorio

En el Laboratorio Químico del INIGEMM, las muestras fueron secadas a 60ºC,
pulverizadas y disueltas en agua regia para luego ser sometidas a análisis por ICP-OES
(Espectroscopia de Emisión Atómica con Plasma de Acoplamiento) para 34 elementos (Hg-
Ag-Al-As-B-Ba-Be-Ca-Cd-Ce-Co-Cr-Cu-Fe-K-Li-Mg-Mn-Mo-Na-Nb-Ni-Pb-Sb-Sc-Se-Sn-
Sr-Ti-Tl-V-W-Zn-Zr). La concentración de Au es determinada mediante Ensayo al Fuego en
el orden de ppb.
4.3. Análisis Geoestadístico de Datos Geoquímicos

Comprende el procesamiento de la información y tratamiento geoestadístico de los datos
geoquímicos. Además de las correlaciones y el análisis multi-elemental mediante coeficientes
de correlación de Pearson.
Se realizó un primer análisis geoestadístico general con el total de muestras que cubre el área
de estudio. En un segundo análisis, se dividió el área de estudio según los Ambientes
Geológicos. Para ello se seleccionan muestras que espacialmente se encuentren ubicadas
dentro de cada dominio y se realiza el análisis individual para cada Ambiente.
30
4.4. Análisis Geoestadístico mono-elemental
4.4.1. Análisis Exploratorio de Datos Espaciales (AEDE)
La primera etapa del estudio de variables georreferenciadas corresponde al AEDE, el
cual permite identificar: valores extremos, distribución, variabilidad, localización y
correlación de los datos con el fin de establecer la validez de aplicación de la teoría
geoestadística y el procedimiento de predicción más conveniente.
Para el AEDE se realiza una estadística descriptiva y gráficos exploratorios para el mejor
entendimiento de los datos. Tablas y diagramas de frecuencia definen medidas estadísticas
como: tendencia central, medidas de posición no central, medidas de dispersión y medidas de
forma. Diagramas Box-Plot (cajas y bigotes) permiten identificar observaciones extremas o
valores atípicos. Gráficos de dispersión son útiles para la detección de relaciones entre las
variables y para la identificación de tendencias en el valor promedio de la variable en la
región.
4.4.2. Tratamiento de datos

Un análisis geoestadístico es óptimo, únicamente si la variable en estudio sigue una
distribución normal (Campana de Gauss). De ser necesario, la muestra debe ser tratada.
Para comprobar si la variable posee una distribución normal, se aplican pruebas de

normalidad con respecto a varios factores (media, mediana, moda, coeficiente de variación y
coeficiente de sesgo). Si la curva no cumple con estas propiedades se realizan
transformaciones a las variables para que la variación de datos no sea representativa.
4.4.3. Transformación de datos

Algunos métodos de interpolación como el Kriging Ordinario, Simple y Universal,
requieren que los datos se ajusten a una distribución normal. Las transformaciones más
comunes son: logaritmo neperiano, raíz cuadrada (box-cox) y arcoseno; su elección depende
del grado de asimetría y ubicación de los valores extremos.
4.4.4. Tratamiento de los valores atípicos (outliers)

Los datos atípicos distorsionan el tratamiento estadístico, debido a que pueden aportar
ruido a la muestra. Su identificación se lo realiza seleccionando aquellos con mayor número
de veces su desviación estándar con respecto a la media o gráficamente mediante diagramas
31
box-plot. Una vez identificados, deben ser eliminados o substituidos por los valores máximos
o mínimos con el objetivo de acotar la muestra.
4.4.5. Análisis de tendencia de los datos

El análisis de tendencia proporciona gráficos en tres dimensiones que ayudan a
identificar tendencias globales de los datos en función de la longitud y latitud, para plantear
un modelo de comportamiento de las variables. Este análisis ayuda a establecer la mejor
dirección de correlación.
4.4.6. Modelamiento del semivariograma

Permite analizar la autocorrelación espacial de los datos. Se elige el modelo teórico que
mejor se ajusta a la distribución de los datos para predecir valores y calcular errores.
4.4.7. Validación del método de interpolación

Cuando se predicen valores en lugares no muestrales, se dispone de dos herramientas
para comprobar la validez del modelo: validación cruzada y validación. La Validación
cruzada compara una curva de regresión teórica con la curva de regresión de los valores
predecidos. La interpolación se considera como válida, cuando los resultados del error medio
estandarizado y el error cuadrático medio son cercanos a cero y el error cuadrático medio
estandarizado es cercano a uno.
4.4.8. Determinación de umbrales geoquímicos

Para la determinación de anomalías, umbral anomálico y valor de fondo se utilizan
métodos estadísticos analíticos y gráficos.
Lepeltier (1969), define al fondo geoquímico (Background) como los valores inferiores a la
media geométrica y al umbral anomálico (Threshold) como la media geométrica más dos
veces la desviación estándar, mismo que corresponde aproximadamente al 2.5% más alto de
la frecuencia acumulada. Este criterio se basa en una distribución normal de la variable y en
que los valores mayores al umbral son controlados por procesos diferentes a los que originan
los valores cercanos a la media (Figura 18). También plantea el uso de hojas de probabilidad
log-normalizadas y que el umbral anomálico corresponde al punto de inflexión en la curva de
concentración versus probabilidad acumulada.
32
Figura 18: Parámetros Geoquímicos
4.4.9. Definición de anomalías de los elementos

En este apartado se analiza el tratamiento estadístico y la representación de las muestras
anómalas de los elementos de interés económico y elementos indicadores de mineralización.
Se generan mapas de anomalías geoquímicas mono-elementales para elementos mayores e

indicadores: Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba. Los rangos de anomalías, valor del
background y threshold, obtenidos del paso anterior, son plasmados en el mapa mediante la
reclasificación de rangos.
4.5. Análisis Geoestadístico multi-elemental

El objetivo del análisis multi-elemental es reducir la dimensión del gran número de variables
(elementos) a un número más reducido (grupos de elementos). El análisis realizado se basa en
la superposición espacial de anomalías geoquímicas, mismas que permiten plantear conjuntos
de elementos genéticamente relacionados. Las asociaciones obtenidas pueden traducirse en la
existencia de asociaciones mineralógicas, correlacionables a un estilo de mineralización o a
un tipo yacimiento mineral en específico.
Para establecer asociaciones multi-elementales a partir de resultados de química de roca de

muestras de la zona, se utilizó el método de Factores de correlación de Pearson, con este
método se puede establecer conjuntos de elementos que se correlacionan y que son
dependientes entre sí.
Los resultados de las correlaciones se presentarán a manera de tablas de correlación, y en

mapas de distribución de asociaciones multi-elementales.
33
5. PRESENTACIÓN Y ANÁLISIS DE DATOS
5.1. Geología
5.1.1. Ambientes Geológicos
Las unidades litoestratigráficas presentes en el área de estudio han sido agrupadas en 2
Ambientes Geológicos considerando el criterio litológico, tectónico y geocronológico.
El Ambiente Volcánico y Magmático Cretácico constituye el 75% de la zona occidental y

abarca a las rocas volcánicas y volcanosedimentarias de las Unidades Pilares y Punta de
Piedra, relacionadas a un arco volcánico oceánico del cretácico y a las rocas intrusivas de
composición intermedia emplazadas en el cretácico superior asociadas al Batolito de Tangula.
Las unidades volcánicas toleíticas y los intrusivos continentales no han sido separados en dos
ambientes distintos debido a su edad similar, extensión irregular, y a que los drenajes cortan
indistintamente a ambas litologías, haciendo difícil definir la procedencia del sedimento al
dividir las muestras en poblaciones.
El Ambiente Volcánico Continental Paleógeno constituye el 25% de la zona oriental del área
y abarca las rocas volcánicas y volcanoclásticas de afinidad calco-alcalina asociadas al
volcanismo de Sacapalca; además de stocks subvolcánicos recientes. Figura 19
Figura 19. Mapa de ubicación de muestras de sedimentos fluviales por Ambientes Geológicos
5.1.2. Indicios de mineralización

INIGEMM (2016b) reconoce veinte indicios de mineralización en el área de Macará,
principalmente de tipo epitermal, pórfido y VMS. Tabla 5 y figura 20
34
Tipo de Posible sistema de
Sector X Y Asociación
ocurrencia mineralización
El Papayal 633837 9503798 Indicio Au-Cu-Ag Epitermal
Las Vigas 631274 9507127 Indicio Au-Ag-Cu (As-Sb-Ba); Cu-Au Epitermal
Porotillo 637315 9511595 Indicio Au-Cu (Pb-Zn) Epitermal
Paltos 630875 9509602 Indicio Au, Cu Epitermal
El Coco 626117 9512110 Indicio Au, Cu Epitermal
El Carmen 621950 9515372 Indicio Au Epitermal
Mandalá 620603 9520731 Indicio As (Ag-Pb-Au) Epitermal
Sauce 637659 9507136 Indicio Ba (Pb-Cu-Zn) VMS
Delicia Sur 636916 9513803 Indicio Ag (Cu-Mo-Sb) Pórfido
Llanitos 640699 9513316 Indicio Au-Ag; Zn (Ag-Cu-As) Epitermal
Punta De Piedra 645363 9516535 Indicio Ag-Au-Cu (Pb-Zn-As-Sb) Epitermal
San Vicente 645149 9517067 Indicio Mo (Ag-Pb-As); Zn-Pb Pórfido/Epitermal
Trapiche 656230 9506852 Indicio Au-Cu-Pb; Ag-Au(As-Pb-Zn) Epitermal
Nongora 658552 9512019 Indicio Ag-Au-As(Pb-As) Epitermal
Naypungo 665223 9515284 Indicio Au Mineralización
Artón 650408 9515792 Indicio Au Mineralización
Taparuca 656560 9515498 Indicio Au Mineralización
Q. Los Paltones 662707 9505314 Indicio Au-Ag-Mo-Sb Mineralización
Curichanga 626393 9521124 Indicio Au Pórfido
Sabiango 632733 9520332 Indicio Ba-As (Cu-Pb-Zn) VMS
Tabla 5. Indicios minerales en la zona de Macará. Modificado de INIGEMM (2016b)
La mineralización epitermal se caracteriza por la presencia de vetas y vetillas de 0,10m hasta

1.5m de ancho, con dirección preferencial NE-SW y NW-SE y que cortan mayormente a las
rocas volcánicas de la U. Punta de Piedra. Están constituidas por cuarzo, calcita, óxidos de
hierro y sulfuros como pirita y calcopirita principalmente, a veces bornita y calcosina.
Malaquita se presenta a veces como mineral secundario. La principal alteración relacionada es
la propilítica y en menor proporción argílica (INIGEMM, 2016b).
Pórfidos dioríticos y andesíticos cortados por vetillas mineralizadas con pirita, calcopirita,
magnetita; además de diseminados se encuentran principalmente en la cercanía del Batolito de
Tangula y presentan comúnmente alteración propilítica de bajo grado y argílica (INIGEMM,
2016b).
Depósitos tipo VMS podrían estar restringidos a la U. Punta de Piedra asociada a volcanismo
oceánico. La ocurrencia de gossans en la parte norte del área, sector Sabiango (INIGEMM,
2016b) podrían estar relacionados a este tipo de depósitos.
5.1.3. Indicios de Oro en batea

Los trabajos de bateo en drenajes se realizan para verificar la existencia de oro aluvial.
En la Figura 20 se muestran los resultados de partículas de oro en batea recopilados de los
35
principales trabajos de exploración realizados en la zona, lo cual indicaría la presencia de
depósitos de oro primario dentro del área de estudio.
Las mayores concentraciones de partículas de oro se encuentran difundidas en drenajes de la

parte central y oriental. Resaltan las concentraciones en los afluentes de la Q. Paltos al oeste
de La Victoria, Q. San Fernando al norte de Sabiango, Qs. San Lapo y Pedregal en la parte
central; y al este, en las Qs. La Paccha, San Vicente y Quisanga, afluentes del Rio Pindo.
SABIANGO
LA VICTORIA
PERÚ
Figura 20: Mapa de ubicación de indicios minerales y oro aluvial. INEMIN (1988) e INIGEMM (2016b)
5.1.4. Alteraciones
La influencia de los cuerpos intrusivos y subvolcánicos emplazados en la U. Punta de
Piedra, define en la roca de caja, alteraciones hidrotermales de baja temperatura. La alteración
de mayor ocurrencia en esta unidad es la propilítica de incipiente a pervasiva que se desarrolla
hacia el centro del área. Figura 21
La alteración argílica eventualmente está acompañada de la alteración propilítica y es muy

localizada (sector El Papayal, sur del área), donde se presentan cantidades menores de arcilla
(caolinita y esmectita) reemplazando a las plagioclasas y acompañada por cuarzo secundario y
pirita fina. En los bordes del intrusivo, es común encontrar silicificación, donde los minerales
primarios que conforman las andesitas basálticas se encuentran parcial o totalmente
remplazados por sílice (INIGEMM, 2016a).
En ciertos sectores hacia el este, las tobas acidas contienen pirita diseminada e incluso se
identifica la presencia de las arcillas alunita y dickita confirmadas mediante análisis espectral
36
(INIGEMM, 2016a), las cuales son características de zonas de alta temperatura, que al estar
acompañadas de sílice configuran una alteración argílica avanzada.
Figura 21: Mapa de alteraciones. INIGEMM (2016a)
5.2. Geoquímica
Los datos químicos analizados corresponden a 661muestras de sedimentos fluviales con
valores de concentraciones para los elementos Cu, Pb, Zn, Au, Ag, Hg, Sb, As, Mo, W y Ba.
Las muestras poseen una buena distribución espacial, cubriendo toda el área de estudio con
una densidad promedio de 0.9 muestras por km2. Figura 22
Figura 22: Mapa de ubicación de muestras de sedimentos fluviales
Se realizó un primer análisis estadístico general para toda el área con el total de todas las
muestras y un segundo análisis por ambientes geológicos con el objetivo de contrastar los
resultados. Tabla 6
37
Área Nº de Densidad de
Ambientes Geológicos Unidad/Formación
(km2) muestras muestreo
• Unidad Los Pilares
Ambiente Volcánico y
504.6 • Unidad Punta de Piedra 472 0.94
Magmático Cretácico
• Batolito de Tangula
Ambiente Volcánico
238 • Formación Sacapalca 189 0.79
Continental Paleógeno
Total 742.6 661 0.90
Tabla 6. Resumen de unidades litoestratigráficas y número de muestras por ambientes geológicos
5.3. Análisis Geoestadístico General

5.3.1. Análisis exploratorio de datos espaciales (AEDE)
Se ha realizado una estadística descriptiva con el objetivo de visualizar y entender la
variabilidad, distribución y correlación espacial de los datos químicos en toda el área de
estudio, para ello se han elaborado: histogramas de frecuencia y gráficos de tendencia de
datos.
La determinación de los outliers (valores atípicos) para cada elemento, se realizó gráficamente
mediante diagramas box-plot o cajas y bigotes (Figura 23). Se han suprimido valores atípicos
extremos para el caso de los elementos: Cu, Pb, Zn, Hg, Au, Sb, Mo y Ba, al considerar que
pueden distorsionar fuertemente la estadística.
Figura 23. Diagramas Box-Plot de los elementos Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba. Análisis General
38
En la tabla 7 se presentan los principales parámetros estadísticos que describen la distribución
y variabilidad de los 11 elementos estudiados.
Parámetros Cu Pb Zn Au Ag Hg As Sb Mo W Ba
Estadísticos (ppm) (ppm) (ppm) (ppb) (ppm) (ppm) (ppm) (ppm) (ppm) (ppm) (ppm)
Datos Analizados 659 659 660 660 661 660 661 660 660 661 659
Mínimo 3,900 5,500 5,600 1,000 1,600 0,100 2,800 3,047 0,900 1,000 22,610
Máximo 227,100 115,900 355,000 927,000 13,970 0,639 63,110 29,450 10,960 17,560 462,100
Media 58,193 16,696 61,891 35,026 1,957 0,109 6,720 10,958 2,302 6,321 120,554
Mediana 51,070 13,210 57,655 12,000 1,600 0,100 2,800 10,695 1,947 6,173 106,700
Moda 40,240 5,500 25,330 1,000 1,600 0,100 2,800 2,800 0,900 1,000 106,700
Varianza 906,52 148,45 1335,67 9057,73 1,537 0,00128 66,343 13,533 2,735 6,949 3786,223
Desviación Est. 30,109 12,184 36,547 95,172 1,240 0,036 8,145 3,679 1,654 2,636 61,532
Q1 38,590 9,386 38,540 3,000 1,600 0,100 2,800 8,305 1,223 4,561 81,463
Q2 71,540 20,170 75,363 22,000 1,710 0,100 7,123 13,138 2,616 7,946 145,500
Coef. de Sesgo 1,750 3,076 2,927 6,001 6,593 7,399 3,315 0,668 2,544 0,409 1,642
Curtosis 5,046 14,647 15,653 41,978 50,928 81,700 13,476 1,125 7,974 0,674 3,966
Coef. de Variacion
51,739 72,976 59,050 271,720 63,346 32,963 121,204 33,571 71,840 41,703 51,041
(%)
Tabla 7. Medidas estadísticas de los elementos. Análisis General
Con las medidas estadísticas obtenidas y los histogramas de frecuencia elaborados (Figura
24), se verifica la normalidad de la distribución de los datos. Se concluye qué, la media,
mediana y moda en los elementos Cu, Pb, Zn, Au, As, Sb, W y Ba, no presentan valores
similares; a diferencia de los elementos Ag, Hg, Mo, cuyos valores son similares o varían en
una unidad.
El coeficiente de sesgo es mayor a cero en todos los elementos, por lo que sus histogramas de
frecuencia presentan asimetría positiva; a excepción del W, cuyo histograma es casi simétrico.
El coeficiente de variación de los datos oscila entre 32% y 73%; para los elementos Cu, Pb,
Zn, Ag, Hg, Sb, Mo, W y Ba, lo cual indica que los efectos causados por los valores extremos
son aceptables. Para el caso del As, el CV es superior al 100% lo cual se considera aun
tolerable. El Au presenta un CV de 245%, lo que indica que los datos extremos tienen una
fuerte influencia, por lo que los datos deberían ser transformados.
39
Cu Pb Zn
C u
C u
PB PB
Au Ag Hg
C u C u
C u
PB PB PB
As Sb Mo
C u
C u C u
PB PB PB
W Ba
C u
C u
PB PB
Figura 24: Histogramas de frecuencia de los elementos. Análisis General
Por lo expuesto, es necesario realizar una transformación de tipo Logaritmo de base 10 para
los 11 elementos, con el objetivo de que la distribución de los mismos sea normal, es decir
adquiera forma de Campana de Gauss.
5.3.2. Medidas estadísticas de los datos transformados

Con las medidas estadísticas de los datos transformados (Tabla 8), se verifica la
normalidad de los mismos; y se concluye qué, la media, mediana y moda de todos los
elementos tienen valores similares con variaciones menores a 0.5. El CS mejora ampliamente
respecto a los datos originales, con excepción del Hg, As y Ag, cuyos histogramas presentan
un fuerte sesgo a la izquierda, debido al gran número de muestras con concentraciones bajas.
El CV varía entre 8% y 70% para todos los elementos, lo que indica que los valores extremos
no presentan mayor influencia.
40
Log Log Log Log Log Log Log Log Log Log Log
Parámetros Estadísticos
Cu Pb Zn Au Ag Hg As Sb Mo W Ba
Datos Analizados 659 659 660 660 661 660 661 660 660 661 659
Mínimo 0,591 0,740 0,748 0,000 0,204 -1,000 0,447 0,484 -0,046 0,000 1,354
Máximo 2,356 2,064 2,550 2,967 1,145 -0,194 1,800 1,469 1,040 1,245 2,665
Media 1,713 1,145 1,731 0,976 0,260 -0,976 0,666 1,014 0,285 0,753 2,031
Mediana 1,708 1,121 1,761 1,079 0,204 -1,000 0,447 1,029 0,289 0,790 2,028
Moda 1,605 1.107 1,834 1.049 0,361 -0.869 0.956 1.049 0.310 0.82 2.047
Varianza 0,046 0,061 0,055 0,472 0,018 0,007 0,105 0,023 0,061 0,050 0,044
Desviación Est. 0,214 0,246 0,234 0,687 0,135 0,085 0,324 0,153 0,247 0,224 0,209
Q1 1,586 0,972 1,586 0,477 0,204 -1,000 0,447 0,919 0,088 0,659 1,911
Q2 1,855 1,305 1,877 1,342 0,233 -1,000 0,853 1,119 0,418 0,900 2,163
Coef. de Sesgo -0,221 0,555 -0,297 0,162 3,605 4,336 1,386 -0,492 0,572 -1,291 -0,090
Curtosis 1,094 0,253 1,069 -0,286 15,809 21,615 0,890 0,522 0,164 2,251 0,451
Coef. de Variación (%) 12,488 21,514 13,496 70,392 51,844 8,761 48,627 15,049 86,713 29,707 10,305
Tabla 8. Medidas estadísticas de los elementos transformados a Log10. Análisis General
Los histogramas de frecuencia de los datos transformados se aproximan en lo posible, a la

distribución normal. Las curvas para los elementos Cu, Zn, Sb, y Ba son leptocúrticas y
presentan simetría, mientras que el W presenta asimetría negativa. Las curvas de Au, Pb, Hg,
Mo, Ag y As presentan curvas leptocúrticas con asimetría positiva; para el caso del Hg y Ag
esto se debe al gran número de datos con concentraciones bajas que presentan esto elementos.
Figura 25
Cu Pb Zn
C u C u
PB PB
Au Ag Hg
C u
C u C u
PB PB PB
As Sb Mo
C u C u C u
PB PB PB
41
W Ba
C u C u
PB PB
Figura 25. Histogramas de frecuencia de los elementos transformados a Logaritmo de base 10. Análisis General
5.3.3. Tendencia de los datos

Útil para definir tendencias globales de los datos. Las curvas de tendencia para Pb, Zn y
Mo, Ba y Au indican una dirección de los datos hacia el sureste; el W presenta tendencia
hacia el norte y el As hacia el noreste; por lo que estos datos presentan una anisotropía en esa
dirección. Por otro lado, las curvas de tendencia de los elementos Cu, Sb, Ag, Hg, Au no
poseen una tendencia bien definida, con crecimientos de concentraciones hacia la parte
central, lo que se interpreta como un comportamiento isotrópico. Anexo 2
5.4. Análisis Geoestadístico por Ambientes Geológicos
5.4.1. Ambiente Volcánico y Magmático Cretácico (Ambiente 1)

5.4.1.1. Análisis exploratorio de datos espaciales (AEDE)
Se ha determinado y suprimido de la estadística al menos un valor atípico extremo para
los elementos Cu, Pb, Zn, Ba, As, al considerar que pueden distorsionar fuertemente la
estadística (Figura 26). En la tabla 9 se presentan los principales parámetros estadísticos que
describen la distribución geoquímica de los elementos estudiados y localizados en el
Ambiente 1.
42
Figura 26. Diagramas Box-Plot de los elementos Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba.
Ambiente 1 (Z1) y Ambiente 2 (Z2)
Con las medidas estadísticas obtenidas y los histogramas de frecuencia, se verifica la

normalidad de la distribución de los datos, y se concluye qué, la media, mediana y moda en
los elementos Cu, Pb, Zn, Mo, W, Ba, Au no son similares; a diferencia de los elementos Sb,
Ag, As y Hg cuyos valores varían en una unidad. El coeficiente de sesgo es mayor a 0 en
todos los elementos, por lo que sus histogramas de frecuencia presentan asimetría positiva.
Figura 27
Parámetros Cu Pb Zn Mo W Sb Ba Au Ag As Hg
estadísticos (ppm) (ppm) (ppm) (ppm) (ppm) (ppm) (ppm) (ppb) (ppm) (ppm) (ppm)
Datos analizados 469 470 470 472 472 472 470 471 472 469 472
Mínimo 3,9 5,5 5,6 0,9 1 3,307 22,61 1 1,6 2,8 0,1
Máximo 206,5 94,45 280,2 10,92 17,56 29,45 409,2 927 13,92 26,61 0,304
Media 63,200 14,195 54,916 2,170 6,431 11,129 120,438 34,955 1,865 3,641 0,104
Mediana 59,13 11,43 48,1 1,8785 6,286 10,85638 105,1849 12 1,6 2,8 0,1
Moda 40.24 5.5 25.33 0.9 1 11.11 106.7 1 1.6 2.8 0.1
Varianza 1046,144 111,820 1181,549 2,405 7,699 14,685 4093,613 9135.013 1,488 5,959 0,000
Desviación est. 32,344 10,574 34,374 1,551 2,775 3,832 63,981 95,577 1,220 2,441 0,020
Q1 39,985 8,099 32,853 1,098 4,483 8,373 75,503 2,000 1,600 2,800 0,100
Q2 79,265 16,335 67,140 2,507 8,250 13,308 149,500 22,000 1,600 2,837 0,100
Coef. de sesgo 1,307 3,734 2,692 2,735 0,269 0,719 1,375 5,970 7,058 4,828 6,345
Curtosis 2,841 20,122 11,826 9,671 0,178 1,158 2,341 41,160 54,872 30,921 46,195
Coef. de
51,178 74,495 62,593 71,468 43,145 34,435 53,124 273,426 65,416 67,044 19,206
variación (%)
Tabla 9. Medidas estadísticas de los elementos. Ambiente 1
El coeficiente de variación de los datos oscila entre 19% y 75%; lo cual se considera
aceptable e indica que los efectos causados por los valores extremos son tolerables. En el caso
del Au, el CV es de 288% lo que indica que los datos extremos tienen una fuerte influencia.
43
Por lo expuesto, es necesario realizar una transformación de tipo Logaritmo de base 10 para
los 11 elementos, con el objetivo de que la distribución de los mismos sea normal.
Cu Pb Zn
Au Ag Hg
As Sb Mo
W Ba
Figura 27: Histogramas de frecuencia de Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba. Ambiente 1
5.4.1.2. Medidas estadísticas de los datos transformados

Con las medidas estadísticas de los datos transformados (Tabla 10), se concluye que la
media, mediana y moda de todos los elementos tienen valores similares con variaciones
menores a 0.3.
El CS mejora ampliamente respecto a los datos originales, con excepción del Ag, Hg y As,
cuyos histogramas presentan un fuerte sesgo a la izquierda, debido al gran número de
muestras con concentraciones bajas de estos elementos. El CV oscila entre 5% y 92% para
todos los elementos, lo que indica que los datos extremos no presentan una fuerte influencia.
44
Parámetros Log Log Log Log Log Log Log Log Log Log Log
estadísticos Cu Pb Zn Mo W Sb Ba Au Ag As Hg
Mínimo 0,591 0,740 0,748 -0,046 0,000 0,519 1,354 0,000 0,204 0,447 -1,000
Máximo 2,315 1,975 2,447 1,038 1,245 1,469 2,612 2,967 1,144 1,425 -0,517
Media 1,745 1,081 1,673 0,261 0,755 1,020 2,024 0,969 0,241 0,518 -0,987
Mediana 1,772 1,058 1,682 0,274 0,798 1,036 2,022 1,079 0,204 0,447 -1,000
Moda 1.790 1.032 1.682 0.320 0.847 1.047 2.050 1.072 0.371 0.614 -0.917
Varianza 0,052 0,053 0,057 0,059 0,057 0,024 0,050 0,476 0,015 0,026 0,003
Desviación est. 0,229 0,231 0,240 0,242 0,239 0,155 0,225 0,690 0,124 0,162 0,057
Q1 1,602 0,908 1,517 0,041 0,652 0,923 1,878 0,301 0,204 0,447 -1,000
Q2 1,899 1,213 1,827 0,399 0,916 1,124 2,175 1,342 0,204 0,453 -1,000
Coef. de sesgo -0,481 0,803 -0,076 0,608 -1,301 -0,430 -0,138 0,174 4,826 2,712 5,323
Curtosis 0,987 0,931 0,798 0,237 1,925 0,372 0,081 -0,312 25,952 7,440 30,451
Coef. de variación
13,106 21,370 14,316 92,570 31,648 15,231 11,097 71,21 51,451 31,341 5,763
(%)
Tabla 10. Medidas estadísticas de los elementos transformados a Log10. Ambiente 1
Los histogramas de frecuencia de los datos transformados se aproximan a una distribución

normal. Las curvas para los elementos Cu, Zn, Ba son leptocúrticas y presentan simetría; las
curvas de Pb, Mo, Au presentan curvas leptocúrticas con asimetría positiva; los elementos W,
Sb presentan asimetría negativa. Las curvas para los elementos Ag, Hg y As poseen asimetría
positiva y un pico en su distribución en la primera familia, debido a la gran cantidad de datos
con baja concentración. Figura 28
Cu Pb Zn
Au Ag Hg
Sb Mo
As
45
W Ba
Figura 28. Histogramas de frecuencia de los elementos Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba
transformados a Logaritmo. Ambiente 1
5.4.1.3. Tendencia de los datos

Las curvas de tendencia para Au, Ag, Cu, Pb, Zn, As, Sb, Ba indican una dirección de los
datos hacia el sureste; W presenta tendencia hacia el norte y Mo hacia el sur; lo que indica
que estos datos tienen una anisotropía en esa dirección. Por otro lado, las curvas de tendencia
de Hg no poseen una tendencia bien definida de concentración de los datos, con crecimientos
hacia la parte central, lo que se interpreta como un comportamiento isotrópico. Anexo 3
5.4.1.4. Modelo geoestadístico

Para el caso de los elementos Pb, Zn, Ba y Au, se seleccionó el modelo Gaussiano, para el
Cu, Mo, Sb y W se eligió el modelo Exponencial y para los elementos Hg, As, Ag el modelo
Esférico (Figura 29). En la Tabla 11 se resumen los parámetros de los modelos de
semivariogramas escogidos.
Cu Pb
Zn Mo
46
W Ba
Sb Au
Ag Hg
As
Figura 29. Semivariogramas experimentales (cruces azules) ajustados al semivariograma teórico (línea azul) de
los elementos Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba. Ambiente 1
Elemento Modelo Partial Sill Nugget Sill Lag Zise Ánisotropía (º) Rango (m)
Cu Exponencial 0.032 0.0186 0.0506 500 20º 4000
Pb Gaussiano 0.028 0.025 0.53 500 190º 4500
Zn Gaussiano 0.033 0.028 0.061 500 150º 5000
Mo Exponencial 0.028 0.025 0.053 500 133º 2500
Sb Exponencial 0.013 0.010 0.023 500 25º 2000
W Exponencial 0.023 0.018 0.041 500 310º 4000
Ba Gaussiano 0.023 0.021 0.04 500 120º 4000
Ag Esférico 0.010 0.003 0.013 500 95º 2000
As Esférico 0.018 0.015 0.033 500 112º 3000
Au Gaussiano 0.065 0.038 0.103 500 64º 3000
Hg Esférico 0.002 0.001 0.003 500 - 2000
Tabla 11. Parámetros de los modelos de semivariogramas elegidos. Ambiente 1
47
Los variogramas presentan un rango promedio de 3273m, distancia a la cual se estabiliza la
curva. La pepita no supera el 50% de la meseta, lo que indica que no existen mayores errores
en la toma de los datos, procesamiento o análisis químico.
5.4.1.5. Validación de la interpolación

La validación de la interpolación se realizó mediante Validación Cruzada. Como se
observa en el Anexo 4, ambas curvas presentan variación en la pendiente debido a los errores
en la predicción. En el caso de los elementos Ag, As y Hg, la variación es más alta, debido a
la alta concentración de valores bajos, quedando reducido el número de valores medios y altos
para ser comparados. En todos los casos, el error medio cuadrático y error medio estándar son
lo más bajos posible y el error medio cuadrático estándar cercano a 1 (Tabla 12); por lo que la
interpolación es considerada válida.
Ecuación de Error medio Error cuadrático Error cuadrático medio

Elemento
regresión estandarizado medio estandarizado
Cu 0.4720X+0.9249 0.0020 0.1818 1.0177
Pb 0.3651X+0.6680 -0.0099 0.2084 1.1443
Zn 0.4406X+0.9323 0.0005 0.1846 1.0114
Mo 0.3348X+0.1597 -0.0053 0.2104 1.0594
Sb 0.2443X+0.7737 -0.0015 0.1406 1.0444
W 0.4020X+0.4695 -0.0005 0.1926 1.1140
Ba 0.5014X-1.0133 -0.0034 0.1713 1.0709
Ag 0.2723X+0.1632 -0.0007 0.1105 1.3022
As 0.1053X+0.4541 -0.0024 0.1883 1.2575
Au 0.1189X-0.8784 -0.0071 0.6621 1.0084

Cu 0.0127X-0.9808 -0.0052 0.0590 1.3360
5.4.2. Ambiente Volcánico Continental Paleógeno (Ambiente 2)

5.4.2.1. Análisis exploratorio de datos espaciales (AEDE)
Se ha determinado y suprimido de la estadística al menos un valor atípico extremo a cada
elemento al considerar que pueden distorsionar fuertemente la estadística (Figura 26; pág. 43).
En la tabla 13 se presentan los principales parámetros estadísticos que describen la

distribución y variabilidad de los elementos estudiados y ubicados en el Ambiente 2.
48
Parámetros Cu Pb Zn Au Ag Hg As Sb Mo W Ba
estadísticos (ppm) (ppm) (ppm) (ppb) (ppm) (ppm) (ppm) (ppm) (ppm) (ppm) (ppm)
Mínimo 12,650 8,176 27,330 1 1,600 0,100 2,800 3,047 0,900 1 37,980
Máximo 77,950 78,530 168,600 302 5,248 0,639 53,410 20,380 10,960 13,490 370,600
Media 44,457 22,242 75,337 25,176 2,084 0,120 13,379 10,495 2,637 5,985 119,410
Mediana 43,840 18,190 73,520 12 1,611 0,100 10,425 10,485 2,096 5,974 108,300
Moda 38,590 12,440 66,310 1 1,600 0,100 2,800 13,750 1,072 4,193 108,300
Varianza 158,855 136,577 452,903 2302,566 0,505 0,003 105,786 10,304 3,488 4,364 2465,048
Desviación est. 12,604 11,687 21,282 47,985 0,711 0,058 10,285 3,210 1,868 2,089 49,649
Q1 36,373 13,715 62,990 3,000 1,600 0,100 6,141 8,118 1,505 4,646 88,663
Q2 51,008 27,32 82,01 21 2,419 0,100 17,425 12,448 2,976 7,229 137,150
Coef. de sesgo 0,203 1,745 1,348 3,845 1,675 5,020 1,592 0,281 2,186 0,327 2,043
Curtosis 0,114 4,042 3,601 15,999 2,649 35,507 2,711 0,176 5,334 0,624 6,537
Coef. de variación
28,351 52,543 28,249 190,600 34,107 48,422 76,877 30,585 70,812 34,902 41,579
(%)
Tabla 13. Medidas estadísticas de los elementos. Ambiente 2
La media, mediana y moda en los elementos Cu, Pb, Zn, Au, As, Sb y Ba, no presentan
valores similares; a diferencia de los elementos Ag, Hg, Mo y W, cuyos valores son similares
o varían en una unidad. El coeficiente de sesgo es mayor a 0 en todos los elementos, por lo
que sus histogramas de frecuencia presentan asimetría positiva, con excepción del Cu, Sb y
W, cuyos histogramas son casi simétricos. Figura 30
El coeficiente de variación de los datos oscila entre 28% y 77%; lo cual se considera
aceptable e indica que los efectos causados por los valores extremos son tolerables. En el caso
del Au, el CV es de 190,6% lo que indica que los datos extremos tienen influencia, pero aun
es tolerable al ser menor al 200%.
Para el caso de los elementos Cu, Sb, W sus CS se aproximan a cero, el CV es menor a 35%,
su media y mediana son similares, por lo que su distribución se acepta como normal, a pesar
de que la moda sea diferente en al menos 1 a 5 unidades. Para el resto de elementos es
necesario realizar una transformación de tipo Logaritmo de base 10, con el objetivo de que la
distribución de los mismos sea normal.
Cu Pb Zn
49
Au Ag Hg
As Sb Mo
Ba
W
Figura 30: Histogramas de frecuencia de Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba. Ambiente 2
5.4.2.2. Medidas estadísticas de los datos transformados

Con las medidas estadísticas de los datos transformados (Tabla 14), se verifica la
normalidad de los mismos; y se concluye que la media, mediana y moda de todos los
elementos tiene valores similares con variaciones menores a 0.3. El CS mejora ampliamente
respecto a los datos originales, con excepción del Hg y Ag, cuyos histogramas presentan un
fuerte sesgo a la izquierda, esto es debido al gran número de muestras con concentraciones
bajas de estos elementos. El CV oscila entre 6% y 73% para todos los elementos, lo que
indica que los datos extremos no presentan una fuerte influencia.
Parámetros Log
Log Pb Log Zn Log Au Log Ag Log Hg Log As Log Ba
Estadísticos Mo
Datos analizados 184 183 182 183 184 184 184 184
Mínimo 0,913 1,437 0,000 0,204 -1,000 0,447 -0,046 1,580
Máximo 1,895 2,227 2,480 0,720 -0,194 1,728 1,040 2,569
Media 1,299 1,861 0,954 0,299 -0,945 1,009 0,342 2,046
Mediana 1,260 1,866 1,079 0,207 -1,000 1,018 0,321 2,035
Moda 1,180 1,857 1.189 0.310 -0.818 0.992 0.301 2.008
Varianza 0,040 0,014 0,419 0,016 0,017 0,106 0,063 0,025
Desviación est. 0,199 0,118 0,647 0,126 0,129 0,325 0,251 0,160
Q1 1,137 1,799 0,477 0,204 -1,000 0,788 0,177 1,948
Q2 1,436 1,914 1,322 0,384 -1,000 1,241 0,474 2,137
Coef. de sesgo 0,512 -0,074 -0,043 1,154 2,792 -0,019 0,510 0,306
50
Curtosis -0,286 1,792 -0,541 0,271 8,426 -0,697 0,073 1,056
Coef. de variación 15,321 6,330 67,861 42,101 13,667 32,235 73,522 7,796
(%)
Tabla 14. Medidas estadísticas de los elementos transformados a Log10. Ambiente 2
Los histogramas de frecuencia de los datos transformados se aproximan a una distribución

normal en lo posible. Las curvas para los elementos Mo, As, Au, Ag y Hg presentan asimetría
positiva y un pico en su distribución en la primera familia con una gran frecuencia, debido a
la gran cantidad de datos de baja concentración. Las curvas del Zn, Ba son leptocúrticas y
presentan simetría. Figura 31
Zn Pb Mo
Au Ag Hg
As Ba
Figura 31. Histogramas de frecuencia de los elementos

Pb, Zn, Au, Ag, Hg, As, Mo y Ba transformados. Ambiente 2
5.4.2.3. Tendencia de los datos

Las curvas de tendencia para Zn, Ag y Ba indican una dirección de los datos hacia el
suroeste; Cu, Pb, W, Au presenta tendencia hacia el noroeste y As hacia el norte; lo que indica
que estos datos tienen una anisotropía en esa dirección. Por otro lado, las curvas de tendencia
de los elementos Mo, Hg y Sb no poseen una tendencia bien definida, con crecimientos de
concentraciones hacia la parte central, lo que se interpreta como un comportamiento
isotrópico. Ver Anexo 5
51
5.4.2.4. Modelo geoestadístico
Para el caso de los elementos Cu, Sb, Ag y As, se seleccionó el modelo Gaussiano, para
el Pb, Zn, Mo, W, Ba y Au se eligió el modelo Exponencial, y para el Hg el modelo Esférico
(Figura 31). En la Tabla 15 se resumen los parámetros de los modelos de semivariogramas
escogidos.
Cu Pb
Zn Mo
W Ba
Sb Au
52
Ag Hg
As
Figura 32. Semivariogramas experimentales (cruces azules) ajustados al semivariograma teórico (línea
azul) de los elementos Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba. Ambiente 2
Los semivariogramas presentan un rango promedio de 3040m, distancia a la cual se estabiliza

la curva. El nugget no supera el 50% de la meseta, lo que indica que no existen mayores
errores en la toma de los datos, procesamiento o análisis químico.
Elemento Modelo Partial sill Nugget Sill Lag zise Anisotropía (º) Rango (m)
Cu Gaussiano 104.37 90 193.37 500 345º 4250
Pb Exponencial 0.0243 0.0170 0.0413 500 345º 3000
Zn Exponencial 0.0071 0.0100 0.0171 500 240º 4000
Mo Exponencial 0.0402 0.0219 0.0621 500 - 1900
Sb Gaussiano 0.0072 0.0170 0.0242 450 - 2800
W Exponencial 3.6381 1.000 4.6381 500 80º 4000
Ba Exponencial 0.0112 0.0143 0.0255 500 230º 4000
Ag Gaussiano 0.0127 0.0090 0.0217 500 215º 2000
As Gaussiano 0.0444 0.0350 0.0794 500 360º 3000
Au Exponencial 0.3849 0.1000 0.4849 500 290º 2500
Hg Esférico 0.0122 0.0080 0.0202 500 - 2000
5.4.2.5. Validación de la interpolación

Como se observa en el Anexo 6 las curvas presentan variación en la pendiente debido a
los errores en la predicción. En el caso de los elementos Ag y Hg, la variación es más alta,
53
debido a la alta concentración de valores bajos, quedando reducido el número de valores
medios y altos para ser comparados.
En todos los casos, el error medio cuadrático y error medio estándar son lo más bajo posible y
el error medio cuadrático estándar cercano a 1 (Tabla 16); por lo que la interpolación se
considera válida.
Error medio Error cuadrático Error cuadrático medio

Elemento Ecuación de regresión
estandarizado medio estandarizado
Cu 0.2183X+35.4173 0.0020 0.1818 1.0177
Pb 0.3024X+0.9007 -0.0075 0.1859 1.0741
Zn 0.1594X+1.5704 -0.0036 0.1316 1.0970
Mo 0.2610X+0.2436 0.0085 0.2212 0.9976
Sb 0.1083X+0.8950 -0.0035 0.1472 1.0208
W 0.3677X+3.7542 1.9277 -0.0269 1.2224
Ba 0.2719X+1.4885 0.0029 0.1477 1.0194
Ag 0.1855X+0.2335 -0.0006 0.1514 1.2731
As 0.5608X+0.4485 -0.0014 0.2531 1.1271
Au 0.0891X+0.8887 -0.0026 0.7241 1.3627
Hg 0.1257X-0.8362 0.0026 0.1493 1.2531
Los umbrales anomálicos (Threshold) se determinaron por el método gráfico de Lepeltier,

(1969), el cual plantea el uso de hojas de probabilidad log-normalizadas y que el umbral
anomálico corresponde al punto de inflexión en la curva de concentración versus probabilidad
acumulada.
Los umbrales se determinaron para cada elemento, para toda el área de estudio para cada
Ambiente Geológico según consta en los Anexos 7, 8 y 9.
54
6. RESULTADOS
6.1. Resultados - Análisis estadístico general

6.1.1. Distribución de los elementos
El análisis general se realizó con el objetivo de conocer el comportamiento y la
variabilidad de los datos en toda el área de estudio. Para la representación gráfica se
elaboraron mapas de concentraciones según los tramos definidos por los percentiles 99%,
95%, 90%, 75% y 50%.
Cobre. Las mayores concentraciones de cobre se ubican en la parte occidental y central; a lo

largo de la Q. Mandalá donde se presentan concentraciones de hasta 206.5ppm
(618350;9519600), y en un tramo de 2km de la Q. Jorupe, donde se presentan concentraciones
de 198.8ppm (622446;9517868) hasta 302.5ppm (622597;9518515). Un valor aislado de
199.1ppm está presente en un afluente de la Q. Curichanga (627018;9517365). La zona
central se caracteriza por un total de 20 muestras que superan las 150ppm ubicadas en su
mayoría en las Qs. Tumbunuma, Zambi, Palacras, Hueco de la Rueda y Gualguana; además
de dos muestras de 227.1ppm (645268;9516734) y 338ppm (645288;9516993) en la Q.
Yanga. La zona oriental se caracteriza por un contenido bajo de Cu, con concentraciones que
no superan las 70ppm. Figura 33
Figura 33: Mapa de distribución de Cobre clasificado por cuantiles. Análisis General
Plomo. Las mayores concentraciones de Pb se ubican en la zona central, en la Q. Zambi

donde se presentan valores de hasta 127.8ppm (643568; 9512858) y 131ppm (644123;
9512436), además de concentraciones aisladas de 94.45ppm (643846; 9520150) en el norte y
55
88.28ppm (636082; 9507228) en el sur. En la parte oriental, las concentraciones oscilan entre
42ppm y 115.9ppm (659780; 9512813) en los afluentes de las Qs. Nongora, San Juan y Santa
Lucia. Las concentraciones más bajas se presentan en el oeste, donde las muestras no superan
las 25ppm, particularmente en las cuencas del Rio Sabiango y Qs. Jorupe, Palto, Vadeal y
Guatara. Figura 34
Figura 34: Mapa de distribución de Plomo clasificado por cuantiles. Análisis General
Zinc. Las mayores concentraciones de Zn se encuentran hacia el sur de la Victoria, en el

límite fronterizo con Perú, a lo largo de quebradas que confluyen en el Rio Macará con
valores de hasta 207ppm (629842; 9505544). En la parte central, en la Q. Zambi se registran
tres valores de 255ppm (643568; 9512858), 258.1ppm (644105; 9512463) y 280.2ppm
(644904; 9512220); y concentraciones aisladas de 321.6ppm (642917; 9518827) y 355ppm
(643335; 9520209) en el norte y 497.3ppm (656015; 9506112) en el sur. Las concentraciones
más bajas se encuentran en la parte oeste alrededor de Macará y en la cuenca del R. Sabiango.
Figura 35
Figura 35: Mapa de distribución de Zinc clasificado por cuantiles. Análisis General
56
Oro. Existen 508 muestras que superan el límite detección, considerándose al resto como
parte del fondo geoquímico. Las mayores concentraciones de oro se encuentran de manera
dispersa hacia el centro y oriente. Un valor aislado de 922ppb se registra en la Q. Curichanga
(625673;9515189) afluente del Rio Sabiango. Concentraciones de 841ppb (641093;9519492)
y 1091ppb (639465;9520675) están presentes en un drenaje en el extremo norte del área y un
valor aislado de 894ppm (656072;9506694) en un afluente de la Q. Purón afluente del Rio
Calvas. Las concentraciones más bajas se ubican en los alrededores de la ciudad de Macará
donde los valores varían entre 1 y 80ppb. Figura 36
Figura 36: Mapa de distribución de Oro clasificado por cuantiles. Análisis General
Plata. Existen 181 muestras que superan el límite detección, considerándose al resto como
parte del fondo geoquímico. Las mayores concentraciones se encuentran hacia la parte este.
Un valor aislado de 13.97ppm se ubica en la Q. Pindo (663934;9503991) en el límite oriental
del área; concentraciones de hasta 13.92 se encuentran en los afluentes de la Q. Samanamaca
(651681;9510217) cercano a su confluencia con el Rio Calvas. Las menores concentraciones
dominan la parte oeste con valores que oscilan entre 1.6ppm y 7ppm. Figura 37
Figura 37: Mapa de distribución de Plata clasificado por cuantiles. Análisis General
57
Mercurio. En general las concentraciones de Hg son bajas debido a la alta volatilidad y
movilidad de este elemento, existiendo un gran número de muestras con concentración
inferior a las 0.1ppm que constituyen parte del fondo geoquímico y que dominan la parte
occidental y central. Existe una muestra aislada con una concentración alta de 1.05ppm
(664803;9516219) en el extremo oriental de la zona en un afluente de la Q. El Paccha.
Concentraciones menores de 0.304ppm (645268; 9516734) y 0.639ppm (660794; 9503687) se
ubican en las Qs. Yanga y La Chirimoya respectivamente. Figura 38
Figura 38: Mapa de distribución de Mercurio clasificado por cuantiles. Análisis General
Arsénico. Existen 368 muestras que superan el límite detección, considerándose al resto
como parte del fondo geoquímico. Las mayores concentraciones de As se encuentran hacia la
parte oriental. En el extremo NE en la Q. Mariposa, se presentan valores de hasta 60.74ppm
(663587;9520016) y concentraciones que oscilan entre 37 y 53ppm en los afluentes de la Q.
Santa Lucia. Un único valor alto de 63ppm se ubica en la Q. El Batán a 8km al NE de Macará
(625488;9518036). Las menores concentraciones se encuentran en la parte central y
occidental. Figura 39
Figura 39: Mapa de distribución de Arsénico clasificado por cuantiles. Análisis General
58
Antimonio. Las mayores concentraciones de Antimonio se encuentran de manera dispersa
hacia la parte norte, en la cuenca alta de las Qs. Mandala, Jorupe, Curichanga y Chingulle con
concentraciones de 23.68ppm (620300; 9520925), 21.87ppm (622489; 9518748), 22.68ppm
(627916; 9518191), 29.45ppm (652750; 9519102) respectivamente. Además de una muestra
aislada con valor de 47.96ppm (663587; 9520016) en la Q. Mariposa en el extremo noreste de
la zona. Las menores concentraciones en encuentran hacia la zona sur. Figura 40
Figura 40: Mapa de distribución de Antimonio clasificado por cuantiles. Análisis General
Molibdeno. Las mayores concentraciones de Mo se encuentran hacia el este, en los afluentes

de la Q. Samanamaca con valores de 10.64ppm (651505;9511337) y 10.92ppm
(651104;9512088). Valores aislados de 10.96ppm (655567;9516481), 12.94ppm
(662850;9515050) y 10.12 ppm (664800;9508050) se ubican en las Qs. San Juan, San Vicente
y La Cucaracha respectivamente. Las menores concentraciones se ubican en la parte central y
occidental, alrededor de las poblaciones de Macará, Sabiango y La Victoria, con valores que
no superan las 4.5ppm. Figura 41
Figura 41: Mapa de distribución de Molibdeno clasificado por cuantiles. Análisis General
59
Wolframio. Las mayores concentraciones de W se ubican en la parte central, con valores de
hasta 13.7ppm (636350; 9511800) en un afluente de la Q. Guando al note del poblado de La
Victoria; y en la parte norte, en la Q. Chingulle donde se presentan valores de hasta 15.81ppm
(652509; 9518890) y 17.56ppm (652750; 9519102). Además de una concentración aislada de
13.88ppm (640962; 9520241). Las concentraciones más bajas se ubican tanto al este como al
oeste de la zona con valores que oscilan entre 1 y 10ppm. Figura 42
Figura 42: Mapa de distribución de Wolframio clasificado por cuantiles. Análisis General
Bario. Las mayores concentraciones de Ba se encuentran hacia el centro de la zona; un

conjunto de valores altos está presente en los afluentes de la Q. Tolungue con valores que
oscilan entre 250 y 761ppm (639202; 9506438) al suroeste de La Victoria; y valores aislados
en las Qs. Yanga y Artón, con concentraciones de hasta 1145ppm (645288; 9516993). Las
menores concentraciones se encuentran hacia la zona oeste en los alrededores de Macará y la
cuenca baja del Rio Sabiango con valores menores a las 150ppm. Figura 43
Figura 43: Mapa de distribución de Bario clasificado por cuantiles. Análisis General
60
6.2. Resultados - Análisis estadístico por Ambientes Geológicos
6.2.1. Interpolación de datos
Se ha realizado una estadística descriptiva univariante de los 11 elementos para cada
Ambiente Geológico por separado, con el objetivo de visualizar y entender la variabilidad,
distribución y correlación espacial del conjunto de datos geoquímicos en las dos poblaciones.
Previo a la interpolación mediante Kriging Ordinario se realizó el análisis estructural de los
datos y se ajustó el modelo de semivariograma para elegir el que genere el menor error en la
predicción de los datos
Para la representación gráfica se elaboraron mapas contornos basados en playas de colores

según los tramos definidos por los umbrales relativos a los percentiles 95%, 90%, 75%, 50%
y 25%. El umbral anomálico para cada elemento se lo ha definido según el método gráfico de
Lepeltier, (1969). El valor de fondo corresponde a la media o al 50% de la frecuencia
acumulada. En la Tabla 17 se comparan los resultados obtenidos en los dos análisis.
General Ambiente 1 Ambiente 2

Background Threshold Background Threshold Background Threshold
Elemento
(ppm) (ppm) (ppm) (ppm) (ppm) (ppm)
Cu 51.1 120,00 59 110,00 43.7 69,50
Pb 13.8 40,00 11.4 27,90 18.2 40,00
Zn 57.5 130,00 42.5 90,00 73.5 115,00
Au (ppb) 15.0 182,30 12 125,60 12 104,10
Ag 2.2 3,20 1.6 3,00 1.61 3,20
Hg 0.13 0,30 0.1 0,18 0.1 0,27
As 4.0 21,20 2.8 10,10 10.3 30,10
Sb 10.7 16,90 10.8 19,00 10.5 14,00
Mo 2.1 4,60 1.9 5,00 2.1 4,80
W 6.2 11,90 6.3 9,80 5.9 9,00
Ba 105.2 250,00 105.2 317.00 108.1 321.00
Tabla 17. Valor de fondo y umbrales anomálicos. Análisis General y por Ambientes Geológicos
6.2.2. Anomalías geoquímicas - Ambiente Volcánico y Magmático Cretácico

Constituye el 75% de la zona occidental y abarca a las rocas volcánicas y
volcanosedimentarias de las Unidades Pilares, Punta de Piedra y las rocas intrusivas del
Batolito de Tangula. Contiene 472 muestras de sedimentos fluviales con una densidad de
0.94muestras/km2.
61
Cobre. Se definieron 7 zonas anomálicas para Cu que superan el umbral de 110ppm presentes
en la parte central y sur; ubicadas dentro de la U. Punta de Piedra misma que posee una
concentración normal 78 a 95ppm (INIGEMM, 2013). Existen 2 zonas de mayor
concentración, la primera se ubica 5km al NE de Macará, está cortada por las Qs. Jorupe y
Mandalá (622446;9517868) y presenta un dominio de concentraciones entre 100 y 302.5ppm.
La segunda se ubica hacia el E de la zona, a 14km de Sabiango, está cortada por la Q. Yanga
y Tumbunuma (645268;9516734) y presenta valores de hasta 338ppm. El valor de fondo
corresponde a 59ppm distribuido al norte de Macará, Cuenca del Rio Sabiango y en el límite
oriental. Figura 44
Figura 44: Mapa de estimación de concentración de Cobre. Ambiente 1
Plomo. Se establecieron 3 zonas anomálicas para Pb (Figura 45) que superan el umbral de
28ppm ubicadas en los bordes del Batolito de Tangula, mismo que presenta una concentración
promedio de 6.4ppm (INIGEMM, 2013).
Figura 45: Mapa de estimación de concentración de Plomo. Ambiente 1
62
La primera anomalía se localiza en el límite norte del área (644479;9520451)con
concentraciones entre 30 y 94ppm. La segunda anomalía se ubica hacia el E, atravezado por la
Q. Zambi (643561;9512949), con concentraciones entre 50 y 131ppm. La tercera se ubica a
2.5km al SE de La Victoria (636737;9507282) y presenta concentraciones hasta 88.23ppm. El
valor de fondo corresponde a 11.4ppm difundido al oeste de la zona.
Zinc. Se definieron 5 zonas anomálicas para Zn que superan el umbral de 90ppm ubicadas
dentro de la U. Punta de Piedra, misma que presenta una concentración normal entre 82 y
108ppm (INIGEMM, 2013).
Las dos primeras se ubican al NE de Macará (622457;9519321) y N de Sabiando

(632284;9520845) respectivamente con concentraciones que oscilan entre 90 y 170 ppm. La
anomalía 3 se ubica en el extremo N (643412;9520512) y presenta concentraciones hasta
355ppm. La anomalía 4 se halla cortada por las Qs. Zambi y Pedregal (641968;9513289), y
presenta concentraciones hasta 260ppm. La quinta anomalía se ubica al sur de La Victoria, se
encuentra cortada por las Qs. Simoras, Hueco de la Rueda y Palacras y presenta
(630553;9505034) concentracines hasta 207ppm. El valor de fondo corresponde a 42.5ppm
difundido alrededor de Macará, cuenca del rio Sabiango y en el extremo este de zona. Figura
46
Figura 46: Mapa de estimación de concentración de Zinc. Ambiente 1
Oro. Se definieron 3 zonas anomálicas para Au que superan el umbral de 125.6ppb ubicadas
dentro de la U. Punta de Piedra. La primera anomalía se ubica junto a la confluencia del Rio
Sabiango y Macará (626378;9512463), presenta concentraciones hasta 198ppb. La segunda
anomalía se ubica 7km al NE de Sabiango (639809; 9520226) y presenta concentraciones
63
sobre las 800ppb. La tercera nomalía se ubica el NE de La Victoria, está atravesada por la Q.
Limón (638173;9511924) y presenta concentraciones superiores a 450ppb. El valor de fondo
corresponde a 12ppb difundido tanto al NE, NW como S del área. Figura 47
Figura 47: Mapa de estimación de concentración de Oro. Ambiente 1
Plata. Se definieron 5 zonas anomálicas para Ag que sobrepasan el umbral de 3.4ppm

ubicadas hacia la parte oriental dentro del Batolito de Tangula. Dos anomalías se ubican en el
límite N, a 7Km al NE de Sabiando con concentraciones entre 5 y 12ppm. Una anomalía
ubicada en la parte central (643365;9513114) presenta concentraciones entre 3.5 y 10.3ppm.
Dos anomalías ubicadas al SE, cerca a la confluencia de la Q. Samanamaca y el R. Calvas,
presenta concentraciones entre 4.8 y 13.9ppm. El valor de fondo corresponde 1.6ppm que
cubre mayormente la zona oeste. Figura 48
Figura 48: Mapa de estimación de concentración de Plata. Ambiente 1
64
Mercurio. En general las concentraciones de Hg son bajas, debido a la alta mobilidad y
volatilidad del elemento. Se reconocen 6 zonas anomálicas para Hg que superan el umbral de
1.8ppm ubicadas en su mayoría hacia el NE. La primera se ubica cerca a la frontera con Perú,
al sur de La Victoria, atravesada por la Q. Palacras (631983;9504027). Cuatro anomalías se
ubican al E de Sabiango, y presentan concentraciones superiores a las 0.2ppm. La sexta
anomalía se ubica en el extremo NE cortado por la Q. Chingulle (651152;9520392) y
presenta concentraciones superiores a 0.2ppm. El valor de fondo corresponde a 0.10ppm
difundido en la mayor parte del área. Figura 49
Figura 49: Mapa de estimación de concentración de Mercurio. Ambiente 1
Arsénico. Se definió 5 zonas anomálicas para As que superan el umbral de 10.1ppm ubicadas
tanto al sur como al norte. Las anomalías 1 y 2 se ubican hacia el S (630521;9505582 y
644411;9509402) y presentan concentraciones hasta 23ppm. Las anomalías 3, 4 y 5 se ubican
hacia la parte N y corresponden a los valores atípicos extremos que varían entre 40 y 63ppm.
El valor de fondo corresponde a 2.8ppm distribuido en la zona central. Figura 50
Antimonio. Existen cuatro zonas anómalas de Sb que superan el umbral de 18.9ppm. La

primera se ubica 5km al NE de Macará (623807;9520101), está atravezada por la Q. Jorupe y
presenta concentraciones sobre las 20ppm. Dos anomalías se ubicas a 8km al E de Sabiango
(639682;9519936 y 639136;9517605) y presenta concentraciones entre 21y 23ppm. La cuarta
anomalía se ubica en el extremo NE de la zona (663842;9520052), está cortada por la Q.
Chingulle y posee concentraciones hasta 29ppm. El valor de fondo corresponde 10.8ppm
difundido en gran parte del área. Figura 51
65
Figura 50: Mapa de estimación de concentración de Arsénico. Ambiente 1
Figura 51: Mapa de estimación de concentración de Antimonio. Ambiente 1
Molibdeno. Se identifican 4 zonas anómalas para Mo que superan el umbral de 5ppm

ubicadas sobre el Batolito de Tangula, mismo que presenta una concentración promedio
inferior a 2ppm (Pilatasig et al., 2005). La primera anomalía se ubica hacia el N, está cortada
por al Q. Urama (645327;9517727) y presenta concentraciones hasta 8.7ppm. La segunda
anomalía se ubica hacia el SW de La Victoria en la zona fronteriza con Perú
(628659;9506953), está cortada por las Qs. Simoras y Hueco de la Rueda y presenta
concentraciones hasta 7ppm. La anomalía más importante se ubica hacia el este, está cortada
por las Qs. Anton y San Juan, cerca a su desembocadura en el R. Calvas (651478;9511029),
se caracteriza por concentraciones superiores a 10ppm. El valor de fondo corresponde 1.9ppm
distribuido en la parte central y norte del área. Figura 52
66
Figura 52: Mapa de estimación de concentración de Molibdeno. Ambiente 1
Wolframio. Se establecieron 5 zonas anomálicas para W que sobrepasan el umbral de 9.8ppm

y que en su mayoría se ubican hacia la zona N, mayormente dentro del Batolito de Tangula
que presenta una concentración promedio de 2ppm (Pilatasig et al., 2005). La primera
anomalía se ubica 5 km al N de Macará (618610;9520544), se encuentra atravesada por la Q.
Paccha y agrupa concentraciones entre 11 y 12.2ppm. La segunda anomalía se ubica a 10km
al NE de Sabiango (639727;9520097), y presenta concentraciones entre 12 a 13.88ppm. La
anomalía 3 se ubica en el centro del área, a 2.5km al NE de La Victoria (636370;9511945),
está atravesada por la Q. Limón, afluente del R. Sabiango y presenta concentraciones entre 11
y 13.7ppm. Las anomalías 4 y 5 se ubican en el extremo NE del área, respectivamente con
concentraciones que oscilan entre 11 y 17.56pm. El valor de fondo corresponde a 6.3ppm
distribuido en la zona sur. Figura 53
Figura 53: Mapa de estimación de concentración de Wolframio. Ambiente 1
67
Bario. Se identifican 3 zonas anomálicas que superan el ubral de 317ppm ubicadas en la
parte este dentro de la U. Punta de Piedra, misma que presenta concentraciones normales
entre 200 y 300ppm. Concentraciones altas que caen dentro del Batolito de Tangula no han
sido consideradas como anómalas debido a que el mismo presenta concentraciones promedio
entre 680 y 780ppm (Pilatasig et al., 2005). La primera anomalía se ubica hacia el N, a 10km
hacia al E de Sabiango (643398;9520461), presenta concentraciones hasta 410ppm. La
anomalía 2 se encuentra atravezada por las Qs. Pugllo y Sauce, y presenta concentraciones
entre 250 y 760ppm. La tercera anomalía se encuentra cortada por la Q. Anton
(638766;9507882) y presenta concentraciones que varían entre 230 y 370ppm. El valor de
fondo corresponde a 105.2ppm difundido hacia la parte oeste y norte. Figura 54
Figura 54: Mapa de estimación de concentración de Bario. Ambiente 1
6.2.3. Anomalías geoquímicas - Ambiente Volcánico Continental Paleógeno

Constituye el 25% de la zona oriental del área y abarca las rocas volcánicas y
volcanoclásticas de afinidad calco-alcalina asociadas al volcanismo de Sacapalca. Contiene
189 muestras de sedimentos fluviales con una densidad de 0.79 muestras/km2.
Cobre. Se identifiacan 2 anomalías para Cu que superan el umbral de 69.5ppm. La primera se

ubica en el W (655657;9516993), cortada por las Qs. Pinine, Santa Lucia y San Juan; y
presenta concentraciones hasta 105ppm . La segunda anomalía se presenta en el extremo NE
del área, a 2km de la ciudad de Cariamanaga (664547;9519893) con concentraciones hasta
78ppm. El valor de fondo corresponde a 43.7ppm presente mayormente hacia el sur. Fig 55a
Plomo. Se reconocen 3 zonas anomálicas para Pb que superan el umbral de 40ppm. La

primera presenta concentraciones hasta 116ppm y está atravezada por las Qs. Piedra Liza y
68
Cofradía (656800;9512442). La segunda anomalía presenta concentraciones entre 43 y 68ppm
y se encuentra atravezada por las Qs. Santa Lucia y Nongora (660166;9512336). La tercera se
ubica en el extremo NE de la zona a 3km de Cariamanga (664378;9520253), presenta
concentraciones hasta 78ppm y está atravezada por la Q. Mariposa. El valor de fondo
corresponde a 18.2ppm distribuido tanto al sur como al norte. Figura 55b
a b
Figura 55: Mapa de estimación de concentración de Cobre y Plomo. Ambiente 2
Zinc. Cuatro zonas anomálicas para Zn superan el umbal de 115ppm. La primera se ubica en
la parte central, está atravezada por las Qs. San Juan y Potrerillos y presenta concentraciones
hasta 160ppm. La segunda anomalía se ubica cerca a la frontera con Perú (656266;9505878),
presenta concentraciones que oscillan entre 118 y 497ppm y es cortada por las Qs. Puron y La
Chirimoya, afluentes del R. Calvas. Dos anomalías ubicadas en el extremo este se encuentran
cortadas por las Qs. Cucaracha y Sanguillin respectivamente con valores entre 100 y 120ppm.
El valor de fondo corresponde a 73.5ppm distribuido en la zona central. Figura 56a
Oro. Se estableciron 4 zonas anomálicas para Au que superan el umbral de 104ppb. La

primera anomalia se ubica al E de Cariamanga (664666;9519966), esta atravezada por la Q.
Mariposa y presenta valores hasta 237ppb. La anomalia 2 se ubica en la zona central
(659031;9512142) y presenta concentraciones entre 300 y 600ppb. Una tercera anomalia se
ubica en la zona central-sur, cortada por la Q. Cachaco (658559;9509399), presenta
concentraciones sobre las 200ppm. La cuarta anomalía se ubica en el límite fonterizo con
Perú (655824;9506692), con concentraciones hasta 800ppb. El valor de fondo corresponde a
12ppm distribuido en gran parte del área. Figura 56b
69
a b
Figura 56: Mapa de estimación de concentración de Zinc y Oro. Ambiente 2
Plata. Se reconocieron 3 anomalías que sobrepasan el umbral de 3.2ppm. La anomalía 1 está

atravesada por la Q. Cachaco (656237;9508470) y agrupa concentraciones entre 3.3 y
5.2ppm. La anomalía 2 está cruzada por la Q. San Vicente (663354;9514818) y presenta
concentraciones hasta 10.3ppm. La tercera anomalía se ubica al SE, (664013;9503832) está
atravesada por la Q. Calvas con concetraciones hasta 13.97ppm. El valor de fondo
corresponde a 1.61ppm difundido en la parte norte. Figura 57a
Mercurio. Se identifican 3 zonas anomalicas para Hg que superan el umbral de 0.27ppm. Las
anomalias 1 y 2 se ubican en los alrededores de Cariamanga con concentraciones superiores a
0.29ppm. La anomalía 3, está cortada por la Q. Paccha (664680;9515778) y presenta
concentraciones hasta 1.05ppm. El valor de fondo corresponde a 0.10ppm difundido en la
mayro parte del área. Figura 57b
Arsénico. Las concentraciones de As en el Ambiente 2 son mas altas que en el Ambiente 1,

presentandose valores hasta 60.7ppm. Existen 3 zonas anomálicas para As que superan el
umbral de 30ppm. La anomalía 1 ubicada en el extremo NE de la zona, al E de Cariamanga
(664584;9520120) está atravezada por la Q. Mariposa con concentraciones superiores a
30ppm. La anomalía 2 ubicada en la naciente de la Q. La Cofradía (660806;9515191) con
concentraciones hasta 50ppm. La tercera anomalía cortada por la Q. Santa Lucia
(656861;9512651) y presenta concentraciones sobre las 40ppm. El valor de fondo
corresponde a 10.3ppm difundido en la parte sur. Figura 58a
70
Antimonio. Cuatro zonas anomálicas sobrepasan el umbral de 14ppm. La primera anomalía
se ubica a 2km al E de Cariamanga, está atravezada por la Q. Mariposa y presenta
concentraciones hasta 48ppm. La segunda anomalía se ubica en el limite E de la zona,
atravezada por la Q. Camayos (656456;9520231) y posee concentraciones entre 16 y 20 ppm.
Dos anomalía se ubican en el extremo E del área en las Qs. San Joaquín y Camayos
respectivamente. El valor de fondo corresponde a 10.5ppm el cual domina la parte central de
la zona. Figura 58b
a b
Figura 57: Mapa de estimación de concentración de Plata y Mercurio. Ambiente 2
a b
Figura 58: Mapa de estimación de concentración de Arsénico y Antimonio. Ambiente 2
71
Molibdeno. Se definieron 3 zonas anomálicas para Mo que superan el umbral de 4.8ppm. La
primera anomalía está cortada por la Q. San Vicente (662957;9515175) y presenta
concentraciones por encima de las 9ppm. Dos anomalías se ubican en el extremo este del área
y presenta concentraciones entre 7 y 10ppm. El valor de fondo corresponde a 2.1ppm
difundido en la parte sur y norte. Figura 59a
Wolframio. Existen cuatro zonas anomálicas para W que sobrepasan el umbral de 9ppm. En
su mayorá se ubican en la parte y agrupan concentraciones sentre 10 y 11.55ppm. Una cuarta
anomalía se ubica en el sur, en la zona fornteriza con Perú, y agrupa concentraciones hasta
17.32ppm. El valor de fondo corresponde a 5.9ppm difundido al sur del área. Figura 59b
a b
Figura 59: Mapa de estimación de concentración de Molibdeno y Wolframio. Ambiente 2
Bario. Se identifican tres zonas anomálicas para Ba que superan el umbral de 321ppm. La
primera cubre está atravezada por la Q. San José (655275;9510682) y presenta
concentraciones dispersas hasta 370ppm. La segunda anomalía se ubica en la zona fronteriza
con Perú (656433;9505618) y presenta concentraciones entre 180 y 462ppm.
El valor de fondo corresponde a 108.1ppmdifundido en la parte norte. Figura 60
72
Figura 60: Mapa de estimación de concentración de Bario. Ambiente 2
6.3. Calidad del agua del sitio de muestreo

El análisis de la calidad del agua de los sitios de muestreo de sedimentos fluviales de Macará,
indica la existencia de 401 muestras tomadas en lugares con agua cristalina que representa el
61% del total de muestras y 260 muestras de lugares con agua turbia que representa el 39% de
las muestras. El mayor número de muestras con agua turbia se ubica hacia el este en el
Ambiente 2, esto puede influir en la calidad de la muestra y en los resultados del análisis
químico puesto que el agua turbia indica exceso de caudal de agua y lavado de material lo que
puede incrementar la concentración de los elementos en el sedimento. Considerando esto, se
eliminó en el análisis estadístico los valores atípicos extremos con el objetivo de evitar errores
y sesgos en los resultados. Figura 61
Figura 61: Mapa de calidad del agua del sitio de muestreo
73
7. INTERPRETACIÓN Y DISCUSIÓN
En base al análisis estadístico realizado según ambientes geológicos, se han definido

anomalías geoquímicas para cada elemento, mismas que han sido agrupadas por
superposición de las mismas y asociaciones que existen entre estos elementos. El mapa de
anomalías geoquímicas se presenta en la Figura 62.
7.1. Áreas con Potencial Metalogénico

Con base en los resultados del análisis de la geoquímica, geología, alteraciones, indicios
de mineralización y por medio de la superposición de mapas se definieron dieciséis áreas
prospectivas con potencial metalogénico en el área de estudio, mismas que se detallan en la
Tabla 18.
Ambiente Áreas Superficie Posible mineralización

Elementos
Geológico prospectivas (km2) asociada
Área 1 9.96 Cu, Zn, Sb Epitermal
Área 2 11.36 Zn VMS
Área 3 9.50 Cu, Ag, Au, Hg Sb, W Epitermal
Área 4 9.50 Pb, Zn, Ag, Ba, Hg VMS/Epitermal
Ambiente Área 5 11.50 Cu, Mo, Ag, Ba, As Epitermal/Pórfido
Volcánico y Área 6 6.82 Hg, As, Sb, W Epitermal
Magmático Área 7 16.00 Cu, Zn, Mo, Hg, As Epitermal
Cretácico Área 8 9.00 Au, W Pórfido
Área 9 10.80 Cu, Pb, Ba VMS
Área 10 10.10 Zn, Pb, Hg Epitermal
Área 11 6.54 Ag, Mo Pórfido
Área 12 10.30 Cu, Zn, Sb Epitermal
Ambiente
Área 13 8.82 Au, Cu, Pb, As, Hg, Sb Epitermal
Volcánico
Área 14 12.42 Au, Pb, As Epitermal
Continental
Área 15 6.16 Au, Zn, Ba, W Epitermal
Paleógeno
Área 16 11.16 Zn, Sb, Mo Epitermal/Pórfido
Tabla 18. Áreas prospectivas de la zona de Macará y sus asociaciones elementales.
Se establecieron once zonas prospectivas dentro del Ambiente Volcánico y Magmático

Cretácico, mismo que según las asociaciones elementales definidas, presenta aptitud para la
ocurrencia de depósitos tipo VMS hospedados en rocas de la U. Punta de Piedra de afinidad
toleítica y depósitos tipo pórfido y epitermales posiblemente asociados a los distintos pulsos
intrusivos del Batolito de Tangula, esta mineralización posiblemente se encuentre hospedada
en las rocas volcánicas de la U. Punta de Piedra y en el mismo Batolito. Figura 63
74
ANOMALÍAS GEOQUÍMICAS
AMBIENTE 2
AMBIENTE 1
PERÚ
Figura 62: Mapa de anomalías geoquímicas de Macará. Ambientes 1 y 2
Au, Cu, P, AS, Hg, Sb, W
Figura 63. Mapa de áreas prospectivas con potencial metalogénico de Macará
75
Por otra parte, dentro del Ambiente Volcánico Continental Paleógeno, se definieron cinco
áreas prospectivas con predisposición a la ocurrencia de depósitos de tipo pórfido y epitermal
posiblemente hospedados en las secuencias volcánicas y volcanoclásticas relacionadas al
volcanismo de Sacapalca.
A continuación, se describe cada área prospectiva:
Área 1. Se ubica al NE de la ciudad de Macará, abarca mayormente la microcuenca de la Q.

Jorupe. Presenta una asociación de Cu, Zn y Sb posiblemente relacionada a mineralización de
tipo epitermal. Litológicamente se caracteriza por la presencia de tobas y lavas andesito-
basálticas correspondientes a la U. Punta de Piedra.
A 1.4km hacia el E, en la Q. Mandalá, INIGEMM (2016b) reporta la presencia oro aluvial

además de rodados de riolita con alteración argílica y presencia de óxidos de hierro; y rodados
de andesita con alteración propilítica además de pirita y calcopirita diseminados. Análisis
químicos a muestras de roca del sector indican concentraciones de Cu hasta 357ppm, Au
520ppb, Ag 10.44ppm, Pb 122ppm, Zn 1143ppm, Sb 49.66ppm y As 202ppm.
El análisis multi-elemental por factores de correlación de Pearson (Tabla 19) a los resultados
químicos de 6 muestras de roca del sector de Mandalá, indica una fuerte correlación entre: As
(Ag, Pb, Mo, Au) y Sb (Cu-Ba); asociaciones posiblemente relacionadas a depósitos tipo
epitermal lo que fundamenta los resultados obtenidos a través del análisis de sedimentos
fluviales.
Ag As Au Ba Cu Mo Pb Sb Zn
Ag 1
As 0,8917 1
Au 0,7984 0,5044 1
Ba -0,2669 0,0119 -0,3538 1
Cu -0,2755 -0,0073 -0,3792 0,9937 1
Mo 0,9683 0,8465 0,8796 -0,2802 -0,3100 1
Pb 0,5881 0,8791 0,1055 0,3892 0,3631 0,5252 1
Sb -0,4824 -0,2266 -0,4557 0,9181 0,8871 -0,4534 0,1951 1
Zn -0,4406 -0,4483 -0,2983 -0,3158 -0,2777 -0,3930 -0,5230 -0,3568 1
Tabla 19. Matriz de Correlación de Elementos sector Q. Mandalá
Área 2. Se ubica N del poblado de Sabiango, abarca la parte superior de las microcuencas de
las Qs. San Fernando, Bellavista y Sabiango. Constituye una anomalía de Zn, además de la
presencia de oro aluvial en las Qs. Bellavista y San Fernando.
76
En el sector, INIGEMM (2016b) reporta la presencia de tobas con alteración cuarzo-sericítica
y argílica intercaladas con niveles de óxidos de hierro (hematita y gohetita), aparentemente un
gossan relacionado a la oxidación de la parte externa de un VMS. Análisis químicos a
muestras del sector indican valores de Cu hasta 2562ppm, Au 390ppb, Ag 7.22ppm, Pb
388.9ppm y Zn 2310ppm.
químicos de 7 muestras de roca del sector de Sabiango, indica una fuerte correlación entre:
Ba-As (Cu-Pb-Zn-Mo); asociación posiblemente relacionada a depósitos tipo VMS lo que
fundamenta los resultados obtenidos a través del análisis de sedimentos fluviales.
Ag 1
As -0,3217 1
Au -0,7188 0,5767 1
Ba -0,2139 0,9810 0,4854 1
Cu -0,4730 0,9710 0,7156 0,9244 1
Mo -0,3862 0,9592 0,5626 0,9678 0,9208 1
Pb -0,2549 0,9922 0,5123 0,9894 0,9385 0,9712 1
Sb -0,5375 0,2225 0,8842 0,1646 0,3542 0,2526 0,1764 1
Zn -0,3546 0,8981 0,4996 0,8236 0,9136 0,7650 0,8522 0,0648 1
Tabla 20. Matriz de Correlación de Elementos sector Sabiango
Área 3. Se ubica al este de Sabiango, sector Villa Nueva, abarca las microcuencas de las Qs.
Surumbe y Gualguana. Presenta una asociación de Au, Ag, Cu, Sb, Hg, W posiblemente
asociada a un depósito tipo epitermal. Litológicamente corresponde al Batolito de Tangula y a
la U. Punta de Piedra la cual está afectada por alteración propilítica. En los drenajes que
cruzan el área se reporta oro aluvial.
Área 4. Se ubica al este de Sabiango. Presenta una asociación de Pb, Zn, Ag, Ba, Hg
posiblemente relacionada con una mineralización tipo VMS. Litológicamente el área está
dominada por lavas andesíticas basálticas de la U. Punta de Piedra. Esta asociación también
puede relacionarse a un sistema epitermal.
Área 5. Ubicado hacia el centro de la zona, abarca parte de las microcuencas de las Qs.
Urama, Yanga y Calguamine. Presenta una asociación de Cu, Mo, Ag, Ba, As correlacionable
con mineralización tipo epitermal y/o pórfido. INIGEMM (2016b) reporta en el sector, lavas
77
andesíticas basálticas con alteración propilítica y un pórfido diorítico cortadas y vetillas
mineralizadas de pirita, calcopirita, magnetita y óxidos de hierro. Aledaño al sector, en las Qs.
Yanga y Seca hay presencia de oro aluvial. Análisis químicos a muestras del sector indican
valores de Cu hasta 2367ppm; Au 6100ppb y Ag 416,4ppm.
químicos de 6 muestras de roca del sector de San Vicente, indica una fuerte correlación entre:
Mo (Ag-As-Pb) y Zn-Pb; asociación posiblemente relacionada a depósitos tipo epitermal y/o
pórfido.
Ag 1
As 0,6379 1
Au -0,8099 -0,0979 1
Ba -0,7482 -0,7241 0,4752 1
Cu 0,1576 -0,3986 -0,4571 -0,2928 1
Mo 0,7024 0,8872 -0,3390 -0,6593 -0,4600 1
Pb 0,3318 0,2863 -0,3868 -0,1440 -0,4809 0,6875 1
Sb 0,2773 0,7293 0,2338 -0,6812 0,0233 0,3672 -0,3744 1
Zn 0,2393 0,0439 -0,4399 0,0284 -0,3752 0,4949 0,9663 -0,5952 1
Tabla 21. Matriz de Correlación de Elementos sector San Vicente
químicos de 11 muestras de roca del sector de Punta de Piedra, indica una fuerte correlación
entre: Ag-Au-Cu(Pb-Zn-As-Sb); asociación posiblemente relacionada a depósitos tipo
epitermal lo que fundamenta los resultados obtenidos a través del análisis de sedimentos
fluviales
Ag 1
As 0,8346 1,0
Au 0,9408 0,9240 1,0
Ba -0,3521 -0,4179 -0,3744 1,0
Cu 0,9937 0,8071 0,9177 -0,3428 1,0
Mo 0,3317 0,4891 0,2894 -0,2314 0,3887 1,0
Pb 0,9982 0,8325 0,9447 -0,3426 0,9954 0,3485 1,0
Sb 0,9950 0,8456 0,9281 -0,3748 0,9843 0,3377 0,9879 1,0
Zn 0,8083 0,6195 0,7880 -0,3605 0,8427 0,4349 0,8321 0,7632 1
Tabla 22. Matriz de Correlación de Elementos sector Punta de Piedra
78
Área 6. Se ubica hacia el NE del área, abarca la microcuenca de la Q. Chingulle. Presenta una
asociación de Hg, Sb, As y W, típicos indicadores de mineralización tipo epitermal.
Litológicamente corresponde al batolito de Tangula.
Área 7. Ubicado en el límite S, sector El Papayal, presenta una asociación de Cu, Zn, Mo,
Hg, As posiblemente asociado a una mineralización tipo epitermal y pórfido. Las rocas que
afloran en la zona son de tipo volcánicas, con presencia de brechas volcánicas y tobas
andesíticas con alteración propilítica de la U. Punta de Piedra y volcanosedimentos de la U.
los Pilares. INIGEMM (2016b) reporta en el sector mineralización asociada a vetillas de
cuarzo, pirita, calcopirita, calcita y gohetita con una dirección preferencial NE-SW además de
pirita y calcopirita diseminada en la roca encajante. Análisis químicos a muestras del sector
indican valores de Cu hasta 158,1ppm, Au 2500ppb, Pb 833ppm y Zn 805ppm.
químicos de 3 muestras de roca del sector de Las Vigas, indica una fuerte correlación entre:
Au-Ag-Cu(As-Sb-Ba), Cu-Au y Cu(Pb-Zn); asociación posiblemente relacionada a depósitos
tipo epitermal lo que fundamenta los resultados obtenidos a través del análisis de sedimentos
fluviales
Ag 1
As 0.7051 1
Au 0.8212 0.9837 1
Ba 0.8776 0.9588 0.9942 1
Cu 0.8264 0.9820 1.0000 0.9952 1
Mo -0.9935 -0.6196 -0.7507 -0.8172 -0.7568 1
Pb 0.7808 0.9936 0.9977 0.9848 0.9971 -0.7044 1
Sb 0.7499 0.9979 0.9933 0.9753 0.9922 -0.6695 0.9988 1
Zn 0.7741 0.9947 0.9969 0.9828 0.9962 -0.6968 0.9999 0.9993 1
Tabla 23. Matriz de Correlación de Elementos sector Las Vigas
Área 8. Se ubica en la zona central y presenta anomalías de Au y W posiblemente asociado a

mineralización tipo pórfido, además de la presencia de oro aluvial en los afluentes de la Q.
Limón que cruzan la zona. Litológicamente el área está caracterizada por andesitas con débil
alteración propilítica de la U. Punta de Piedra intruidas por el Batolito de Tangula.
INIGEMM (2016b) reporta en el sector mineralización de vetas y vetillas de cuarzo,
malaquita, óxidos de manganeso y hierro con una dirección preferencial E-W; NE-SW que
79
cortan tanto a los volcánicos como al intrusivo. Análisis químicos a muestras del sector
indican valores de Cu hasta 250ppm y Au 4300ppb.
Área 9. Se ubica hacia el SE de La Victoria, presenta una asociación Cu, Pb, Ba posiblemente
asociada a una mineralización tipo VMS, litológicamente el área de encuentra dominada por
tobas y lavas andesito-basálticas de la U. Punta de Piedra. Dentro del área, INIGEMM
(2016b) reporta mineralización asociada a vetas de cuarzo y vetillas de carbonatos y
oxidación de gohetita y jarosita, con una dirección preferencial NW-SE que se encuentran
cortando las andesitas. Análisis químicos a muestras del sector indican valores de 1309ppm
de Cu, 631ppm de Zn en una veta de cuarzo y sulfuros.
químicos de 3 muestras de roca del sector de Sauce, indica una fuerte correlación entre:
Ba(As-Au), Cu(Zn-Sb), Pb(Au-As-Ba), si bien las asociaciones no están bien definidas,
posiblemente se relacionen a depósitos tipo VMS lo que fundamenta los resultados obtenidos
a través del análisis de sedimentos fluviales
Ag 1
As -0,442 1
Au -0,842 0,856 1
Ba -0,849 0,849 1,000 1
Cu -0,674 -0,365 0,169 0,181 1
Mo 0,873 0,051 -0,473 -0,484 -0,948 1
Pb -0,580 0,987 0,928 0,923 -0,211 -0,109 1
Sb -0,631 -0,417 0,113 0,125 0,998 -0,929 -0,266 1
Zn -0,806 -0,175 0,360 0,371 0,980 -0,992 -0,015 0,968 1
Tabla 24. Matriz de Correlación de Elementos sector Sauce
Área 10. Se ubica en la parte central al NE de la Victoria, presenta una asociación de Zn, Pb,
Hg posiblemente asociado a una mineralización tipo epitermal. Litológicamente el área
presenta un dominio de lavas andesito-basálticas con alteración propilítica de grado medio de
la U. Punta de Piedra cortadas por vetas de cuarzo con pirita, calcopirita y esfalerita, con
direcciones preferenciales NE-SW y E-W; además de stockworks con presencia de limonita y
hematita, intruidas por el Batolito de Tangula (INIGEMM 2016b). En la quebrada s/n del
sector El Pedregal que atraviesa el área se presenta oro aluvial. Análisis químicos a muestras
del sector indican valores de Cu hasta 250ppm; Au 19.23ppm y Ag 53.76ppm.
80
químicos de 13 muestras de roca del sector de Llanitos, indica una fuerte correlación entre:
Au-Ag, Zn(Ag-Cu-As) y As-Pb, posiblemente se relacionen a depósitos tipo epitermal lo que
fundamenta los resultados obtenidos a través del análisis de sedimentos fluviales
Ag 1
As 0,3224 1
Au 0,8430 0,3490 1
Ba -0,4153 0,0786 -0,2619 1
Cu 0,2319 0,2004 -0,1525 -0,1194 1
Mo -0,1082 -0,4396 -0,2498 -0,3667 0,1339 1
Pb -0,2999 0,2952 -0,1939 0,2491 0,3537 -0,1574 1
Sb -0,1398 0,1137 -0,0092 0,1033 0,3051 0,2041 0,4864 1
Zn 0,5319 0,5091 0,3360 -0,3313 0,5347 -0,2846 0,1362 0,1046 1
Tabla 25. Matriz de Correlación de Elementos sector Llanitos
Área 11. Se ubica hacia el E de la zona, abarca la cuenca baja de la Q. Samanamaca, presenta
una asociación de Ag y Mo posiblemente asociada a un depósito tipo pórfido.
Litológicamente el área está dominada por granodioritas asociadas al Batolito de Tangula. Se
registra presencia de oro aluvial en los drenajes que atraviesan la zona.
Área 12. Se ubica al SW de Cariamanga, abarca las microcuencas de las Qs. San Juan y s/n.
Presenta una asociación de Cu, Zn, Sb, W posiblemente asociada a una mineralización tipo
epitermal. Litológicamente el área se encuentra dominada por tobas con alteración argílica y
andesitas con alteración subpropilítica de la Fm. Sacapalca.
Área 13. Se ubica en el extremo NE de la zona, abarca las microcuencas de las Qs. Mariposa
y Arrayán. Presenta una asociación de Au, Cu, Pb, As, Hg, Sb posiblemente asociada a
mineralización tipo epitermal. Oro aluvial se presenta en la Q. Arrayán al sur del área.
Litológicamente el área está dominada por lavas andesíticas con alteración argílica de la Fm.
Sacapalca.
Área 14. Se ubica en la parte central-este abarcando la cuenca media de la Q. Nongora donde
se presenta oro aluvial. Presenta una asociación de Au, Pb, As posiblemente relacionado a
mineralización tipo epitermal. Litológicamente el área está dominada por lavas andesíticas
con alteración argílica y tobas riolíticas de la Fm. Sacapalca cortadas por una falla de
dirección SW-NE. INIGEMM (2016b) reporta en el sector rocas volcánicas alteradas con
81
presencia de hematita, limonita, gohetita y cristales de yeso. Análisis químicos a muestras del
sector indican valores de Cu hasta 242,6ppm, Au 192ppb, As 170.2ppm, Sb 41.54ppm, estos
dos últimos valores respaldan la idea de mineralización tipo epitermal en la zona.
químicos de 4 muestras de roca del sector de Nongora, indica una fuerte correlación entre:
Ag-Au-As(Pb-Sb); posiblemente se relacionen a depósitos tipo epitermal lo que fundamenta
los resultados obtenidos a través del análisis de sedimentos fluviales.
Ag 1
As 0,9607 1
Au 0,6194 0,4100 1
Ba 0,7406 0,8827 0,1591 1
Cu 0,5566 0,6554 0,3493 0,8694 1
Mo -0,7121 -0,7097 -0,6690 -0,7702 -0,9291 1
Pb 0,4861 0,2247 0,9053 -0,1627 -0,0690 -0,3049 1
Sb 0,5041 0,2562 0,9704 -0,0614 0,1129 -0,4696 0,9762 1
Zn 0,5877 0,7566 -0,2712 0,7495 0,3306 -0,1900 -0,3411 -0,3840 1
Tabla 26. Matriz de Correlación de Elementos sector Nongora
Área 15. Se ubica hacia el sur, en la zona fronteriza con Perú, sector Loma Trapiche. Presenta
una asociación de Au, Zn, Ba, W. Litológicamente el área está dominada por tobas y lavas
andesíticas de la Fm. Sacapalca. INIGEMM (2016b) reporta en el sector la presencia de tobas
oxidadas con hematita, limonita, gohetita, óxidos de manganeso y yeso. Análisis químicos a
muestras del sector indican valores de Cu hasta 250ppm; Au 520ppb y Ag 10,44 ppm.
químicos de 4 muestras de roca del sector de Trapiche, indica una fuerte correlación entre:
Au-Cu (Pb), Ag-Au (As-Pb-Zn); posiblemente se relacionen a mineralización tipo epitermal.
Ag 1
As 0,9981 1
Au 1,0000 0,9986 1
Ba -0,1463 -0,1703 -0,1507 1
Cu 1,0000 0,9981 1,0000 -0,1463 1
Mo 1,0000 0,9981 1,0000 -0,1463 1,0000 1
Pb 0,9572 0,9386 0,9548 -0,1265 0,9572 0,9572 1
Sb -1,0000 -0,9981 -1,0000 0,1463 -1,0000 -1,0000 -0,9572 1
82
Zn -0,1776 -0,2356 -0,1859 0,1351 -0,1776 -0,1776 0,1143 0,1776 1
Tabla 27. Matriz de Correlación de Elementos sector Trapiche
Área 16. Se ubica en el extremo E de la zona, abarca las microcuencas de las Qs. Camayos,
Cucaracha y Los Linderos. Presenta una asociación de Zn, Sb, Mo posiblemente relacionada a
mineralización tipo epitermal y/o pórfido. Litológicamente el área está dominada lavas
andesíticas de la Fm. Sacapalca.
Existen varias anomalías aisladas en ambos Ambiente Geológicos que no están relacionadas
con muestras adyacentes o entre elementos, por lo que no fueron incluidas dentro de las áreas
prospectivas; estas anomalías aisladas se detallan en la Tabla 28.
Ambiente Elemento Concentración (ppm) Ubicación X Y

W 10.70-13.50 Q. Paccha 618614 9520573
As 60.74-63.11 Q. El Batán 625449 9518133
Ambiente 1 Au 150-193 (ppb) Q. s/n; E de Macará 624667 9513095
Mo 6.63-10 Q. s/n; Frontera con Perú 648922 9509618
Cu-As 114.6-198.50; 23.89-35.21 Q. s/n; E de La Victoria 644741 9509168
Hg 0.291-0.302 Q. s/n; Sur de Cariamanga 659563 9519119
As 36.87-53.40 Q. La Cofradía 660904 9515565
Ag-Mo 9.46-13.91; 10.12-12.94 Q. San Vicente 663025 9515139
Hg 1.00-1.05 Q. La Paccha 664555 9515901
Ambiente 2
Ag 4.0-7.57 Q. Cachaco 655255 9505930
Au 205-735 (ppb) Q. Cachaco 658326 9509304
Ag 13.92-13.97 Q. Calvas 664091 9503957
Zn 121.3-253.7 Q. Sanguillín 666047 9504338
Tabla 28. Anomalías geoquímicas aisladas. Ambientes 1 y 2
7.2. Discusión
En base a la configuración geológica, tectónica, anomalías y asociaciones geoquímicas,
indicios de mineralización y alteraciones, se establece que el área de estudio presenta afinidad
para la ocurrencia de depósitos epitermales, pórfidos y VMS; lo que concuerda con lo
expuesto por Gossens, (1969); Paladines, (1980); Merlyn y Mosquera, (1988); CODIGEM
(1993) y Pillajo, (2008), quienes plantean la existencia de depósitos tipo pórfido cuprífero y
mineralización de oro en los sectores de Los Linderos, Río Playas, El Huato y Laguar; los
cuales se habrían desarrollado en el cretácico; además de mineralización polimetálica
desarrollada en el paleógeno y neógeno en los sectores de Changaimina, Cariamanga y Los
83
Colorados. Por otra parte, no se ha reportado la presencia de depósitos VMS en los
alrededores de Macará, a diferencia de lo que ocurre en el norte peruano donde existen varios
de estos depósitos (Rodríguez et al., 2008).
En base al criterio litológico, tectónico y geocronológico, la zona de estudio se divide en dos

ambientes geológicos, el Ambiente Volcánico y Magmático Cretácico ubicado al oeste y el
Ambiente Volcánico Continental Paleógeno ubicado al este. Si bien dentro del primer
ambiente coexisten unidades volcánicas toleíticas asociadas a un arco de islas e intrusivos
calco-alcalinos asociados a procesos de subducción (INIGEMM, 2013), estos no fueron
separados en dos ambientes distintos debido a su similitud en cuanto a la edad, extensión
irregular y a que los drenajes cortan indistintamente a ambas litologías, haciendo difícil
definir la procedencia del sedimento fluvial al dividir las muestras en poblaciones.
Dentro de la zona occidental, se definieron asociaciones geoquímicas de Pb-Zn-Ag-Ba-Hg

(Área 2), Zn (Área 4) y Cu-Pb-Ba (Área 9) presentes en los sectores de Sabiango y Sauce,
relacionadas a rocas basálticas de la Unidad Punta de Piedra y que son correlacionables con
mineralización de tipo VMS, como ocurre en el norte de Perú, donde se reporta la ocurrencia
de depósitos VMS de Cu-Zn-Au hospedados en rocas volcánicas submarinas de la Formación
Ereo correlacionable con la U. Punta de Piedra, siendo el más importante el depósito de
Tambogrande particularmente caracterizado por una cobertura extensa (550x200m) de gossan
(INGEMMET, 2012). De acuerdo a INIGEMM (2016b), en el Ecuador, niveles de gossan han
sido reportados en el sector de Sabiango. Las asociaciones geoquímicas de Zn-Pb-Ba
definidas a partir de sedimentos fluviales y rocas, concuerdan con la zona externa de los VMS
caracterizada por enriquecimiento de esfalerita, galena, pirita y baritina (McQueen, s/f) lo cual
indica que, de existir algún depósito, el sistema se encontraría conservado. Sin embargo, las
áreas propuestas dentro del área de estudio se ubican entre los 1500 a 2600m de altitud,
mientras que los principales depósitos VMS del norte de Perú se ubican entre los 100 y
700msnm, esto como consecuencia de un levantamiento del segmento ecuatoriano de la
Cuenca Lancones-Alamor lo que implica que los posibles depósitos estén erosionados, por lo
tanto, no están conservados o son sistemas incompletos ubicados dentro de techos colgantes.
Dentro de esta misma zona, se definieron tres áreas favorables para la ocurrencia de depósitos
tipo pórfido de Cu y Cu-Mo por las asociaciones geoquímicas de Cu-Mo-Ag-Ba-As (Área 5),
W-Au (Área 8) y Mo-Ag (Área 11) ubicados en los límites del Batolito, lo que concuerda con
lo establecido por INIGEMM (2016b), que plantea la presencia de pórfidos dioríticos y
84
andesíticos cortados por vetillas mineralizadas con pirita, calcopirita, magnetita; además de
diseminados, comúnmente con alteración propilítica de bajo grado y argílica. Las alteraciones
presentes, características de la zona externa de un sistema de pórfido (Sillitoe, 2010) indican
que los posibles depósitos se encuentran conservados. Sin embargo, a consecuencia del
levantamiento del segmento ecuatoriano de la Cuenca Lancones-Alamor es posible que estos
sistemas estén erosionados y no se hayan conservado.
Asimismo, se definieron cinco áreas prospectivas para mineralización epitermal

caracterizadas por los típicos indicadores As-Sb-Hg acompañados de Au, Ag, Cu, Pb y Zn
(Áreas 1, 3, 6, 7 y 10) dispuestos de forma dispersa, lo que concuerda con lo establecido por
INIGEMM (2016b), que plantea la presencia de vetas y vetillas con cuarzo, óxidos de hierro,
sulfuros de hierro y cobre, cortando mayormente a las rocas de la U. Punta de Piedra, siendo
la principal alteración relacionada la propilítica y en menor proporción argílica. El contexto
geológico indica que depósitos tipo pórfido y epitermal en el área, podrían asociarse a la
actividad magmática del Batolito de Tangula, sin embargo, aparentemente este intrusivo sería
estéril y no generaría mineralización (Pilatasig y Gordón, 2007). Más bien la mineralización
estaría relacionada al emplazamiento de cuerpos menores y filones dioríticos y pórfido-
andesíticos (INIGEMM, 2013) como ocurre en el sector de Cariamanga (fuera del área) donde
un complejo peribatolítico de composición similar rodea la parte norte y oeste del Batolito
(Egüez y Poma, 2001), el mismo que sería el responsable de la mineralización en las zonas de
Río Playas y El Huato, de manera similar a lo que ocurre en Perú, donde depósitos tipo
pórfido de Cu-Mo y epitermales de Au-Ag están relacionados a múltiples intrusiones del
Complejo Plutónico Las Lomas del cretácico, también caracterizado por un borde gabro-
diorítico que rodea a cuerpos intrusivos ácidos e intermedios (INGEMMET., 2012).
La zona oriental está constituida por rocas volcánicas y volcanoclásticas asociadas al arco
continental Sacapalca del paleógeno. De acuerdo a los modelos planteados por Corbett &
Leach, (1997); Corbett, (2013) y Sillitoe, (2010), este arco representa un ambiente propicio
para la generación de mineralización epitermal y pórfido, lo que concuerda con la
asociaciones de Cu-Zn-Sb, Au-Cu-Pb-As-Hg-Sb, Au-Pb-As, Au-Zn-Ba-W definidas, además
de indicios minerales reconocidos por INIGEMM, (2016b), con lo que se establece cinco
áreas prospectivas con predisposición a la ocurrencia de depósitos de tipo epitermal (Áreas
12, 13, 14 y 15) y una para pórfido (Área 16). De acuerdo a CODIGEM (1993), dentro de esta
franja se registra mineralización polimetálica desarrollada en el paleógeno y neógeno en los
85
sectores de Changaimina, Cariamanga y Los Colorados; similar a lo que se encuentra en Perú
donde se reporta un sub-cinturón mioceno medio-plioceno inferior caracterizado por
mineralización de pórfidos de Cu-Mo (Ag-Pb), skarns de Pb-Zn-Cu-Ag y depósitos
polimetálicos relacionados a intrusivos miocénicos (INGEMMET, 2012).
Basado en los resultados obtenidos, estudios previos y por correlación con lo reportado en el
norte de Perú, se establece al menos tres épocas metalogénicas para la zona sur del país. La
primera del cretácico inferior-tardío caracterizado por un volcanismo toleítico y que habría
generado mineralización tipo VMS de Cu-Pb-Zn y Cu-Zn-Au hospedado en la U. Punta de
Piedra. La segunda con un rango de diferencia de edad corta, del cretácico superior-temprano,
relacionada a la última fase de la actividad magmática del Batolito de Tangula que habría
generado mineralización tipo epitermal de metales preciosos Au, Ag y metales base Cu, Zn,
Pb, Mo; y mineralización tipo pórfido de Cu y Cu-Mo(Au). La tercera época del paleoceno al
mioceno relacionada al volcanismo de Sacapalca, habría originado mineralización tipo
epitermal de Cu-Zn, Au principalmente.
86
8. CONCLUSIONES Y RECOMENDACIONES
8.1. Conclusiones
• El análisis estadístico por ambientes geológicos indicó variaciones representativas en

el valor de fondo y umbral anomálico para los elementos Cu, Pb, Zn, Hg, As, Sb, y
Ba, esto como resultado de las diferentes concentraciones normales de los tipos
litológicos presentes en cada ambiente. Estas variaciones revelaron zonas anómalas
que no son apreciadas al considerar un umbral anomálico único para toda el área en el
análisis global.
• La zona de estudio está constituida por dos ambientes geológicos, el Ambiente

Volcánico y Magmático Cretácico ubicado al oeste, dentro del cual coexisten rocas
volcánicas toleíticas (Unidad Punta de Piedra) asociadas a un arco de islas e intrusivos
tipo I relacionados a procesos de subducción (Batolito de Tangula), y el Ambiente
Volcánico Continental Paleógeno ubicado al este que abarca a las rocas volcánicas y
volcanoclásticas de afinidad calco-alcalina de la Formación Sacapalca.
• Tres épocas metalogénicas se establecen para la zona sur del país. La primera del
cretácico inferior-tardío caracterizado por un volcanismo toleítico y que habría
generado mineralización tipo VMS de Cu-Pb-Zn y Cu-Zn-Au genéticamente
relacionada a la U. Punta de Piedra. La segunda con un rango de diferencia de edad
corta, del cretácico superior-temprano, relacionada a la última fase de la actividad
magmática de Tangula que habría generado mineralización tipo epitermal de metales
preciosos Au, Ag y metales base Cu, Zn, Pb y Mo; y mineralización tipo pórfido de
Cu y Cu-Mo(Au). La tercera época del paleoceno al mioceno relacionada al
volcanismo de Sacapalca y que habría originado mineralización tipo epitermal de Cu-
Zn, Au principalmente.
• Posibles depósitos VMS formados a nivel del fondo marino y que en la actualidad se
ubican entre los 1500 y 2600m de altitud dentro del área de estudio, probablemente no
se encuentren conservados o sean sistemas incompletos ubicados dentro de techos
colgantes a consecuencia del levantamiento, intrusión y procesos erosivos que han
87
afectado al segmento norte de la Cuenca Lancones-Alamor al menos en los últimos
100Ma, lo que le resta interés económico a la zona de Macará, sin embargo, son
necesarios mayores estudios.
• Dentro del Ambiente Volcánico y Magmático Cretácico, se definieron tres áreas

favorables para la ocurrencia de depósitos tipo VMS por las asociaciones de Pb-Zn-
Ag-Ba-Hg y Cu-Pb-Ba presentes en los sectores de Sabiango al norte y Sauce en el
centro, genéticamente relacionados al volcanismo basáltico toleítico. En estos sectores
no se ha reportado afloramientos de sulfuros masivos, pero sí la presencia de gossans
que se relacionan a la parte más superficial de este tipo de mineralización, como
ocurre en territorio peruano (VMS de Tambogrande).
• Dentro del Ambiente Volcánico y Magmático Cretácico se definieron cinco zonas

prospectivas con afinidad a mineralización epitermal por las asociaciones Cu-Zn-Sb,
Cu-Ag; Au-Hg-Sb-W, Hg-As-Sb-W, Zn-Pb-Hg y Cu-Zn-Mo-Hg-As correlacionables
con indicios de mineralización en vetas y vetillas con cuarzo, óxidos de hierro,
sulfuros de hierro y cobre que cortan a las rocas volcánicas de la U. Punta de Piedra.
Además, tres áreas con potencial para mineralización tipo pórfido de Cu y Cu-Mo
definidas por las asociaciones de Cu-Mo-Ag-Ba-As, Au-W y Ag-Mo.
• Mineralización tipo pórfido y epitermal presentes en la zona oeste se relaciona al

emplazamiento de cuerpos menores y filones dioríticos y pórfido-andesíticos presentes
en el área como ocurre en el sector de Cariamanga con el Complejo Peribatolítico de
composición diorítica y las mineralizaciones de Río Playas y El Huato.
• Dentro de la franja volcánica del paleógeno ubicada al este, se establecieron

asociaciones de Cu-Zn-Sb, Au-Pb-As y Au-Zn-Ba-W que definen cinco áreas
prospectivas con predisposición a la ocurrencia de depósitos de tipo epitermal y
pórfido. El volcanismo y magmatismo relacionado a este arco representa un ambiente
propicio para la generación de este tipo de depósitos.
88
8.2. Recomendaciones
• Los resultados obtenidos mediante el análisis estadístico por ambientes geológicos han
sido favorables. Se recomienda la utilización de esta metodología de tratamiento de los
datos por poblaciones en futuras campañas de exploración geoquímica. En campañas
de exploración geoquímica a detalle es factible incluso el tratamiento de los datos
según unidades litoestratigráficas.
• Realizar un estudio tectónico para establecer las tazas de levantamiento, exhumación y

erosión que han afectado al segmento norte de la Cuenca Lancones-Alamor y como
han influenciado en la conservación de posibles depósitos minerales en el sur del país.
• Realizar trabajos de geología en la zona de Macará para establecer la continuidad del

Complejo Peribatolítico hacia el sur y su relación con mineralización tipo pórfido y
epitermal presentes en la zona occidental del área de estudio.
• Considerando las asociaciones geoquímicas e indicios minerales presentes en las

zonas de Sabiango y Sauce (Áreas 2, 4 y 9) que sugieren la posible presencia de
mineralización tipo VMS, se recomienda realizar una campaña de follow-up de las
anomalías geoquímicas aumentando la densidad de muestreo de sedimentos fluviales
sobre una base topográfica escala 1:25 000. Además de una campaña de geofísica con
el método cargabilidad-resistividad y gravimetría, mismos que presentan una buena
respuesta a este tipo de depósitos. Estos métodos han sido empleados con éxito en la
exploración de VMS en el norte peruano. En estas áreas es poco recomendable el
método electromagnético puesto que el carácter basáltico de la Unidad Punta de Piedra
podría enmascarar la presencia de minerales magnéticos en el subsuelo.

zonas de Q. Jorupe, Villa Nueva, Q. Chingulle, El Papayal y La Victoria al oeste y Q.
San Juan, Q. Mariposa, Q. Nongora y Loma Trapiche al este (Áreas 1, 3, 6, 7, 10, 12,
13, 14 y 15) que sugieren la posible presencia de sistemas epitermales, se recomienda
realizar una campaña de follow-up de las anomalías geoquímicas aumentando la
densidad de muestreo de sedimentos fluviales sobre una base topográfica de mayor
escala. El estudio debe ser acompañado con muestreo y análisis químico de rocas para
la determinación de elementos pathfinders (Sb, Hg, As) y análisis espectrales para la
89
determinación de minerales arcillosos indicadores de alteración hidrotermal.
Asimismo, es recomendable realizar estudios geofísicos terrestres de polarización
inducida y resistividad, métodos que presentan buenos resultados en la exploración de
mineralización en vetas.

zonas de San Vicente, Q. Limón y Q. Samanamaca ubicados al oeste y Llano Grande
al este (Áreas 5, 8, 11 y 16) que sugieren la posible presencia de mineralización tipo
pórfido, se recomienda realizar una campaña de muestreo y análisis químico de rocas
para la determinación de elementos pathfinders (Cu, Mo, W, Sn) y análisis espectrales
para la determinación de minerales arcillosos asociados a los distintos halos de
alteración característicos de estos depósitos. También es recomendable realizar una
campaña de gravimetría y magnetometría terrestre en mallas de 200 o 100 m sobre los
targets. Al aplicar métodos magnéticos en estos depósitos, es recomendable prestar
atención a anomalías negativas que pueden estar asociadas a la alteración intensa de la
roca caja y de los minerales magnéticos, especialmente en niveles inferiores.
• Se recomienda realizar una campaña de gravimetría aeroportada a lo largo de la franja

volcánica ubicada al este, lo que permitirá definir posibles intrusiones y cuerpos
mineralizados que se ubican bajo la cobertura volcánica y volcanosedimentaria de la
Formación Sacapalca y que no se evidencian en superficie.
90
9. BIBLIOGRAFÍA
Carranco, J. Ruiz, G. (2015). Prospección Geológica y Geoquímica de la Hoja Topográfica

Zamora (E 1:50.000). Tesis de pregrado, Universidad Central del Ecuador.
CODIGEM-BGS. (1993). Mapa Tectonometalogénico del Ecuador escala 1:1000000.
Corbett, G. (2013). World Gold Pacific Rim Epithermal Au-Ag. Australasian Institute of
Mining and Metallurgy.
Corbett, G., & Leach, T. (1997). Short course manual: Southwest Pacific rim gold-copper
systems: Structure, alteration and mineralization. Special Publications of the Society of
Economic Geologists.
Egüez A. y Poma O. (2001). La Cuenca Alamor-Lancones en el contexto geodinámico de los

Andes de Huancabamba, SW Ecuador. Cuartas Jornadas en Ciencias de la Tierra,
Memorias, Escuela Politécnica Nacional. pp 19-22.
Galley, A., Hannington, M., & Jonasson, I. (2007). Volcanogenic massive sulphide deposits.
Mineral Deposits of Canada: A Synthesis of Major Deposit-Types, District
Metallogeny, the Evolution of Geological Provinces, and Exploration Methods:
Geological Association of Canada, Mineral Deposits Division, Special Publication No.
5, 141-161.
García O, Moreno G. (2006). Herramientas de Análisis Estadístico y espacial aplicadas en la

separación Línea Base-Anomalía Geoquímica. Implicaciones en Exploración de
Recursos Naturales. Parte. Parte 1: Revisión Metodológica. Departamento de
Geociencias. Universidad Nacional de Colombia. Bogotá-Colombia
Giraldo, R. (2002). Introducción a la Geoestadística: Teoría y Aplicación. Departamento de

Estadística, Universidad Nacional de Colombia, 97 p. Bogotá-Colombia
Gossens, P. J. (1969). Mapa índice mineralógico de la República del Ecuador. Quito.
Gutscher M, Malavielle J, Lallemand S, Collot J. (1999). Tectonic segmentation of the North

Andean margin: impact of the Carnegie Ridge collision. Earth and Planetary Science
Letters, 168, France. pp.255-270
Hungerbuhler, D. (1997). Neogene basin in the Andes of southern Ecuador: Evolution,

deformation and regional tectonic implications. Diss. ETHZ No. 12371. Swiss Federal
Institute of Technology Zurich.
INAMHI. (2012). Anuario Meteorológico. Instituto Nacional de Meteorología e Hidrología.

Quito-Ecuador
91
INEMIN. (1988). Sub-proyecto Macará. Informe inédito.
INGEMMET. (1987). Reyes L, Caldas J. Geología de los cuadrángulos de Las Playas, La

Tina, Las Lomas, Ayabaca, San Antonio, Morropón, Huancabamba, Olmos y
Pomahuaca. Boletín, Serie A: Carta Geológica Nacional, N° 39. Lima, Perú
INGEMMET. (2006). Vílchez, M., Núñez, S., Valenzuela, G. Estudio Geoambiental de la

Cuenca del Rio Chira-Catamayo. Boletín Nº 31 Serie C: Geodinámica e Ingeniería
Geológica. Lima-Perú
INGEMMET. (2012). Rodríguez, I., Villareal, E., Valencia, M., y Sánchez, V. Metalogenia,
Geología Económica y Potencial Minero de la Deflexión de Huancabamba: Noroeste
del Perú. Boletín Nº 29-Serie B: Geología Económica. Lima- Perú
INIGEMM. (2013). Ibadango, E., Torres, C., Chiriboga, V., Vélez, T. (2013). Memoria
técnica hoja geológica de Cariamanga. Proyecto “Modelo piloto para la determinación
del potencial geológico-mineralógico de las zonas Zaruma y Cariamanga a escala
1:100.000. Quito-Ecuador.
INIGEMM. (2016a). Burbano, D., Carranco, A. Memoria Técnica Preliminar Hoja Geológica
de Macará Escala 1:100 000. Dirección de Geología-Geología Regional. Quito-Ecuador.
Informe no publicado.
INIGEMM. (2016b). Lomas, W., Pinzón, F., Quiñónez, K. Memoria técnica de ocurrencias de
recursos minerales Bloque 6 (4°S-5°S). Hojas topográficas escala 1:100 000.
(Cariamanga, Gonzanamá, Alamor, Zapotillo, Macará, Las Aradas y Zumba). Informe
no publicado.
Jaillard E, Ordóñez M, Berrones G, Bengstone P, Bonhomme M, Jiménez N, Zambrano I.

(1996). Sedimentary and tectonic evolution of the arc zone of Southwestern Ecuador
during Late Cretaceous and Early Tertiary times. Journal of South American Earth
Sciences 9, pp 131-140
Kellogs J, Bonini W. (1982). Subduction of the Caribbean plate and basement uplifts in the
overriding South American plate: Tectonics, v. 1, pp. 251-276
Kennerley J. (1973). Geology of Loja Province Southern Ecuador. Institute of Geological

Sciences-Overseas Division-Photogeological Unit. London- United Kingdom.
Lebrat M., Mégard F., Dupuy c., Dostal J. (1979). Geochemistry and tectonic setting of pre-
collision Cretaceous and Paleogene volcanic rocks of Ecuador. Geol. Soc. Amer.
Bull.99, pp 569-578.
Lepeltier, C. (1969). A Simplified Statistical Treatment of Geochemical Data by Graphical

Representation. Economic Geology. Vol. 64. pp.538-550.
92
Mateus, A. (2011). Análisis Geoestadístico por el Método de Kriging Ordinario Aplicado a
Muestras de Sedimentos Fluviales en el Distrito Minero Azuay. Tesis de maestría,
Universidad San Francisco de Quito, Quito, Ecuador.
McQueen, K. G. (s/f). Ore deposit types and their primary expressions. Regolith Expression
of Australian Ore Systems, 1-14. Recuperado de
http://www.munkhiinvolcan.mn/Images/1_OreDepositTypes.pdf
Merlyn M, Mosquera L. (1988). Geología y Mineralización del indicio Aurífero de Laguar,

Cantón Macará, Provincia de Loja. Quinto Congreso Ecuatoriano de Geología, Minas,
Petróleos y ciencias Afines. Quito-Ecuador
Moreno A. (2008). Sistemas y Análisis de Información geográfica. Manual de

autoaprendizaje con ArcGIS. Segunda Edición. Alfaomega. México DF
Mourier, T., Laj, C., Roperch, P., Mitourd, P., Farfán Medrano, A., (1988). An accreted
continental terrene in northwestern Peru. Earth and Planetary Science Letters,88, pp.
182-192., Amsterdam, Netherlands.
Paladines A. (1980). Mapa metalogénico del Ecuador. Escala 1.1000 000. Ministerio de
Recursos Naturales y Energéticos. Dirección General de Geología y Minas. Quito-
Ecuador.
Paladines A. (1989). Zonificación Geotectónica y Metalogénica del Ecuador. Quito-Ecuador
Pilatasig L, Gordón D, Palacios O y Sánchez J. (2005). Geología del Ecuador y Perú entre
3°S y 6°S. Proyecto Multinacional Andino: Geociencias para las comunidades Andinas.
Ecuador-Canadá- Perú. DINAGE-INGEMMET
Pilatasig. L, Gordón. D, Palacios. O y Sánchez. J. (2007). Mapa geológico binacional región

sur del Ecuador y norte del Perú. Proyecto Multinacional Andino: Geociencias para las
comunidades Andinas. Ecuador-Canadá- Perú. DINAGE-INGEMMET
Pilatasig, L., Gordón, P. (2007). Mapa Geológico del Ecuador, Escala 1:100 000, Hoja 40
Macará. Dirección Nacional de Geología, Ministerio de Minas y Petróleos, Ecuador
Pillajo, E. (2008). Mapa del Potencia Aurífero Primario de la República del Ecuador. Quito.
PRODEMINCA. (2000a). Evaluación de Distritos Mineros del Ecuador. Depósitos

epitermales en la Cordillera Andina. DINAGE, BGS
PRODEMINCA. (2000b). Evaluación de Distritos Mineros del Ecuador. Depósitos porfídicos

y epi-mesotermales relacionados con intrusiones de las Cordilleras Occidental y Real.
DINAGE, BGS
PRODEMINCA. (2000c). Evaluación de Distritos Mineros del Ecuador. Sulfuros masivos

alijados en volcanitas. DINAGE, BGS
93
PRODEMINCA. (2002). Manual of Exploration of Metalliferous deposits in Ecuador.
Assessment of Ore Districts from Ecuador. DINAGE, BGS
Rodríguez, I., Quispe, J., Sánchez, V., Villarreal, E., Jaimes, F. (2008). Metalogenia de la
cuenca Lancones: Noroeste del Perú-Sur de Ecuador. INGEMMET
Rose A, Hawkes H, Webb J. (1979). Geochemistry in Mineral Exploration. Second Edition.

Academic Press. London
Sillitoe R. (2010). Porphyry Copper Systems. Society of economic Geologist. Economic

Geology. V. 105. Pp, 3-41
Sinclair A. (1974). Selection of threshold values in Geochemical data using probability

graphs. Department of Geological Sciences. University of British Columbia.
Vancouver-Canadá.
Ullrich T. D. (2004). Summary Report on 40ar/39ar Dating for MAP: GAC Activity Ec-50
Ecuador. Samples Dg-56 And Ga-019. Argon Geochronology Laboratory. Pacific
Centre for Isotopic and Geochemical Research Earth & Ocean Sciences, University of
British Columbia.
United Nations. (1971). Survey of Metallic and Non-Metallic Minerals. Systematic

Geochemical Exploration of Operation Nº 8, Quito, Ecuador.
UTPL. (2012). Proyecto: “Generación de Geoinformación para la Gestión de Territorio a

Nivel Nacional Escala 1:25 000”. Geopedología y Amenazas Geológicas. Memoria
Técnica-Subcomponente 2.1: “Geomorfología”.
Villemur J. (1967). Estudio de reconocimiento geológico-mineralógico de la provincial de

Loja. Cooperación Técnica Francesa. Servicio Nacional de Geología y Minas. Quito-
Ecuador
Wipf M. (2006). Evolution of the Western Cordillera and Coastal Margin of Peru: Evidence
from low-temperature Thermochronology and Geomorphology. Swiss Federal Institute of
Technology Zurich
94
10. ANEXOS
MEDIDAS DE TENDENCIA CENTRAL

PARÁMETRO DATOS AGRUPADOS DATOS NO AGRUPADOS
MEDIA
(Fórmula 1.1)
MEDIANA
(Fórmula 1.2)
MODA
(Fórmula 1.3)
95
MEDIDAS DE DIRPERSIÓN
PARÁMETRO DATOS AGRUPADOS DATOS NO AGRUPADOS
DESVIACÓN
ESTANDAR
(Fórmula 1.4)
VARIANZA
(Fórmula 1.5)
COEFICIENTE DE
VARIACIÓN
(Fórmula 1.6)
MEDIDAS DE FORMA
CURTOSIS
(Fórmula 1.7)
COEFICIENTE DE
ASIMETRÍA
(Fórmula 1.8)
Anexo 1. Formulas usadas en estadística
96
Cu Pb
Zn Mo
W Ba
Ag As
Hg Au
Sb
Anexo 2. Gráficos de análisis de tendencia de los datos. Análisis General
97
Au
Ag
Cu Pb
Zn Mo
W As
Sb Hg
Ba
Anexo 3. Gráficos de análisis de tendencia de Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba. Ambiente 1
98
Cu Pb
Zn Mo
W Ba
Sb Au
Ag Hg
As
Anexo 4. Validación cruzada de los elementos Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba. Ambiente 1
99
Cu Au Ag
Pb Zn Mo
W As Sb
Hg Ba
Anexo 5. Gráficos de análisis de tendencia de Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba. Ambiente 2
100
Cu Pb
Zn Mo
W Ba
Sb
Au
Ag Hg
As
Anexo 6. Validación cruzada de los elementos Cu, Pb, Zn, Au, Ag, Hg, As, Sb, Mo, W y Ba. Ambiente 2
101
Anexo 7. Umbrales Anomálicos. Análisis General
102
Anexo 8. Umbrales Anomálicos. Ambiente 1
103
Anexo 9. Umbrales Anomálicos. Ambiente 2
104
Anexo 10. Resultados de análisis químicos de indicios de mineralización en Macará
Sabiango
COORDENADAS Ag As Au Cu Mo Ba Pb Zn Sb Cr Ni
MUESTRAS
X Y ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm
A-BR-056 632693 9520520 1,78 182,9 0,390 860 23,94 623,1 59,77 552 552 93,8 2,1
A-BR-057 632733 9520332 1,60 249,4 0,245 1091 12,56 348,0 55,83 1459 1459 84,6 2,1
A-BR-058 632733 9520332 2,92 1130 0,293 2562 75,18 4723,0 388,9 2310 2310 45,9 2,1
A-BR-043 632490 9521887 5,65 30,24 0,040 197,1 6,71 498,8 12,04 74,17 34,41 307,7 70,16
A-BR-044 632226 9521917 7,22 18,82 0,030 102,6 3,37 40,7 26,57 65,80 32,41 102,9 35,61
A-BR-045 632627 9521402 3,12 5,41 0,026 20,49 15,45 1,6 24,88 33,60 10,38 186,3 25,12
A-BR-046 632504 9521312 3,51 5,37 0,076 379,9 16,44 206,4 12,00 46,12 26,43 200,5 69,33
Las Vigas
COORDENADAS Ag As Au Ba Cr Cu Mo Ni Pb Sb Zn
MUESTRA
A-JS-019 631462 9505821 3.72 16.89 0.345 275.7 65.16 28.55 5.60 9.42 33.15 4.67 15.48
A-JS-020 631274 9507127 4.87 411.40 2.536 631.5 67.65 158.10 4.07 11.53 833.40 24.08 805.40
A-KQ-016 631827 9504872 1.66 62.16 0.14 190.50 125.80 15.13 10.04 6.44 25.46 5.54 17.70
Mandalá
MUESTRA
A-BR-099 620394 9520904 1,6 2,8 0,159 137,4 16,02 14,4 0,9 2,1 11,5 19,60 186,4
A-BR-100 620555 9520890 2,48 97,41 0,083 1644,0 9,82 2193 0,9 2,1 84,2 49,66 61,3
A-BR-101 620603 9520731 2,50 56,83 0,138 179,5 30,87 204,4 4,08 5,87 25,2 3,15 1143,0
A-BR-102 620603 9520731 4,62 49,77 0,193 192,3 34,11 357,0 4,79 6,21 15,5 2,8 362,8
A-KQ-020 619714 9520467 8,78 250 0,19 243,70 128,00 250 24,98 2,1 122,00 4,31 47,40
A-KQ-021 619714 9520467 10,44 202,00 0,52 274,80 70,51 250 39,83 2,1 66,69 4,54 48,02
El Sauce
MUESTRA
A-JS-008 637695 9507634 1,6 2,8 0,057 668,7 34,05 1309,00 0,9 11,36 8,51 20,65 34,89
A-JS-012 636583 9506929 2,05 3,88 0,099 746,0 122,4 257,30 7,60 6,29 23,09 8,56 17,26
A-JS-013 636583 9506929 7,75 2,8 0,012 582,2 183,4 56,29 13,21 3,68 5,5 7,18 7,92
105
Llanitos
MUESTRAS
X ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm ppm
A-KQ-001 641282 9514207 1,63 184,20 0,22 70,14 116,30 103,10 8,74 2,1 57,92 10,02 94,60
A-KQ-002 640629 9513578 3,15 73,50 2,09 131,40 95,81 111,50 5,01 17,55 72,79 3,96 58,66
A-KQ-003 640665 9513300 2,55 130,70 4,38 500 250 66,70 0,9 19,52 31,13 8,81 154,10
A-KQ-004 640667 9513304 1,6 163,00 0,16 500 228,30 54,19 3,70 10,36 11,08 6,39 90,79
A-KQ-005 640669 9513310 3,76 250 1,34 500 149,00 55,99 0,9 6,32 44,68 4,50 51,53
A-KQ-006 640699 9513316 53,76 250 19,23 50,26 80,85 71,70 3,99 29,42 11,76 7,39 173,80
A-KQ-007 640700 9513318 22,72 248,70 0,61 49,62 180,00 250 8,40 18,35 23,83 3,97 250
A-KQ-008 640704 9513319 27,90 113,40 4,75 58,59 54,58 148,50 3,59 13,63 33,81 4,08 250
A-KQ-009 640243 9513087 5,13 180,80 2,04 42,16 250 118,50 2,63 54,54 77,07 12,09 198,50
A-KQ-010 640243 9513087 9,43 250 2,47 500 250 250 8,79 31,07 144,50 14,56 185,40
A-KQ-011 640243 9513087 3,41 107,50 0,53 500 33,66 49,81 3,15 8,20 71,44 4,71 63,57
A-KQ-018 639073 9513005 16,56 34,86 0,49 359,80 107,30 221,10 8,14 2,67 19,95 9,89 71,12
A-KQ-026 639344 9514327 2,23 2,8 0,01 55,21 186,80 50,13 22,47 15,57 5,60 8,02 24,29
Punta de Piedra
MUESTRAS
A-KQ-012 641231 9518892 4,34 27,26 0,02 466,50 130,70 250 13,38 29,51 250 2,3 250
A-KQ-013 641221 9518891 5,15 22,33 0,02 485,70 139,60 250 15,90 25,55 250 3,95 250
A-KQ-014 645363 9516535 7,49 93,26 1,97 328,20 82,67 62,79 3,30 17,10 250 2,3 250
A-WL-004 643099 9518189 5,41 22,38 0,004 127,3 97,06 79,79 7,557 17,23 66,88 11,43 102,8
A-WL-005 641405 9519055 5,41 188,6 1,109 173,3 126,7 79,79 17,61 1,721 24,22 16,9 26,2
A-WL-006 641259 9519092 5,41 22,38 0,03 910,2 19,4 79,79 5,54 5,365 7,41 11,43 24,48
A-WL-007 641257 9519097 5,41 22,38 0,044 450,4 58,44 79,79 5,54 3,003 7,41 11,43 29,11
A-WL-009 644659 9512022 5,41 22,38 0,01 129,7 34,59 205,8 11,86 8,522 33,09 11,43 311,2
A-WL-010 645020 9509453 5,41 22,38 0,018 283,7 29,25 178,7 5,54 8,644 16,53 11,43 39,33
A-WL-011 643017 9513222 416,4 305,7 6,109 53,6 264,1 2367 14,84 6,678 5853 163,5 634,5
A-WL-012 642992 9513231 5,41 22,38 0,011 103,4 51,88 79,79 5,54 12 11,65 11,43 68,56
106
San Vicente
MUESTRAS
A-JS-041 645209 9517020 1,6 6,84 0,042 1689,0 35,28 607,00 0,9 3,54 18,99 19,81 102,60
A-KQ-015 645415 9517294 5,41 37,39 0,011 271,6 34,28 2183 5,54 2,202 10,14 26,21 52,89
A-KQ-016 645149 9517067 5,41 43,72 0,018 1415 23,49 824,4 5,54 1,45 9,028 19,62 74,71
A-KQ-017 645204 9517019 5,41 65,72 0,012 975,9 134,1 220,1 17,54 3,66 175,1 13,55 735,2
A-WL-013 653932 9511650 5,41 157,9 0,035 254,3 108 112 19,39 3,29 47,36 37,96 48,11
A-WL-014 655674 9511718 1,6 6,84 0,042 1689,0 35,28 607,00 0,9 3,54 18,99 19,81 102,60
Trapiche
MUESTRAS
A-KQ-021 656230 9506852 10,44 202,00 0,52 274,80 70,51 250 39,83 2,1 66,69 4,54 48,02
A-KQ-022 656230 9506852 5,41 26,79 0,004 52,89 26,61 79,79 5,54 7,604 26,99 11,43 85,4
A-KQ-023 656221 9506855 5,41 38,67 0,009 63,05 42,16 79,79 5,54 8,334 9,103 11,43 24,48
A-KQ-024 656126 9506731 5,41 28,19 0,004 1176 19,42 79,79 5,54 1,45 18,9 11,43 61,31
Nongora
MUESTRAS
A-JS-030 662976 9515168 6,11 149,00 0,197 113,4 103,10 87,85 11,08 2,1 52,99 41,54 27,10
A-JS-031 663169 9515214 1,6 2,8 0,025 57,9 288,30 17,45 20,58 3,07 13,25 14,80 5,6
A-KQ-025 657515 9511157 5,41 168,5 0,11 191,9 11,63 242,6 5,54 1,45 15,6 21,11 52,9
A-KQ-027 658552 9512019 5,41 170,2 0,021 156,8 22,03 79,79 16,69 3,791 11,56 11,43 93,97
107

T Uce 0012 58 2017

Cargado por

Copyright:

Formatos disponibles

T Uce 0012 58 2017

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

T Uce 0012 58 2017

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD CENTRAL DEL ECUADOR

FACULTAD DE INGENIERÍA EN GEOLOGÍA, MINAS,

CARRERA DE INGENIERÍA EN GEOLOGÍA

“CARTOGRAFÍA GEOQUÍMICA MULTI-ELEMENTAL DE SEDIMENTOS

Trabajo de titulación previo a la obtención del Título de Ingeniero en Geología

Autor: Marco Andrés Vásquez Guañuna

Tutor: M.Sc. Luis Felipe Pilatasig Moreno

Quito, noviembre 2017

esfuerzo me ha brindado la mejor herencia que me pudo dar, la educación.

poder ser ejemplo para ellas.

hubiese sido el mismo.

Agradezco a la Universidad Central del Ecuador, a la FIGEMPA y a cada uno de los

docentes que de manera desinteresada me impartieron sus conocimientos, experiencias y

consejos para mi formación académica.

A todos quienes conforman el INIGEMM por brindarme la oportunidad de desarrollar este

realización de este proyecto.

Al Ing. Danny Burbano por brindarme su apoyo, amistad y consejo.

A la Ing. Rocío Sangucho y al Ing. Salomón Brito por su apoyo y consejo.

formación personal, profesional y a la realización de este trabajo.

El autor declara que la obra objeto de la presente autorización es original en su forma de

Marco Andrés Vásquez Guañuna

CARRERA DE INGENIERÍA EN GEOLOGÍA

APROBACIÓN DEL TRABAJO DE TITULACIÓN POR PARTE DEL TUTOR

En la ciudad de Quito, a los 01 días del mes de noviembre del 2017

M.Sc. Luis Felipe Pilatasig Moreno

CARRERA DE INGENIERÍA EN GEOLOGÍA

APROBACIÓN DEL TRABAJO DE TITULACIÓN POR PARTE DEL TRIBUNAL

En la ciudad de Quito, a los 20 días del mes de noviembre del 2017.

Para constancia de lo actuado firman:

Ing. Alex Mateus

PRESIDENTE DEL TRIBUNAL

Ing. Elías Ibadango Ing. Danny Burbano

MIEMBRO DEL TRIBUNAL MIEMBRO DEL TRIBUNAL

SIGLAS Y ABREVIATURAS ............................................................................................................ xiii

Figura 1: Mapa de ubicación geográfica del área de estudio. ................................................................. 5

Anexo 1. Formulas usadas en estadística .............................................................................................. 96

AEDE: Análisis Exploratorio de Datos Espaciales

ARCOM: Agencia de Regulación y Control Minero

BAT: Bloque Amotape Tahuín

DINAGE: Dirección Nacional de Geología de Ecuador

ICP-OES: Espectroscopía de Emisión Atómica con Plasma de Acoplamiento

INEMIN: Instituto Ecuatoriano de Minería

INGEMMET: Instituto Geológico, Minero y Metalúrgico de Perú

INIGEMM: Instituto Nacional de Investigación Geológico, Minero y Metalúrgico

Ma: Millones de años

ppb: partes por billón

ppm: partes por millón

VMS: Sulfuros Masivos Volcanogénicos (siglas en inglés)

Autor: Marco Andrés Vásquez Guañuna

Tutor: Luis Felipe Pilatasig Moreno

Este estudio presenta los resultados de la exploración y cartografía geoquímica multi-

Geológicamente la zona de Macará se encuentra emplazada en el borde oriental de la Cuenca

La zona occidental del área, litológicamente constituida por rocas predominantemente

PALABRAS CLAVE: MACARÁ / RÍO CALVAS / CARTOGRAFÍA GEOQUÍMICA /

Autor: Marco Andrés Vásquez Guañuna

Tutor: Luis Felipe Pilatasig Moreno

KEY WORDS: MACARÁ / RÍO CALVAS / MULTI-ELEMENTAL GEOCHEMISTRY

Ing. Luis Felipe Pilatasig Moreno