Introduction
The Sierra Nevada de Santa Marta Massif (SNSM) (Figure 1) contains significant deposits of natural resources, including coal, marble, limestone, barite, and iron. Mining and quarrying activities in the region constitute an important source of raw materials and production, primarily through the extraction of these commodities. Among them, iron (Fe) represents a particularly promising target for exploration. Several projects have been initiated to evaluate the geological potential of iron ores, aiming to contribute to the socioeconomic development of the La Guajira department. These initiatives have included ore evaluation studies, extraction processes, and their potential applications.

Source: Author
Figure 1 Left. Satellite Google Earth image of the Sierra Nevada de Santa Marta Massif (SNSM), showing the study area. Right. Satellite Google Earth image showing the distribution of sampling localities for iron ore occurrences hosted in Los Tábanos Rhyodacite of San Juan del Cesar (La Guajira).
However, comprehensive research on the petrology and ore deposit geology of iron occurrences in this region remains limited. In particular, the genetic types of these deposits are still a matter of debate. In this context, mineralogical characterization using integrated analytical techniques is therefore crucial for understanding the genesis of these ores. Iron deposits can undergo multiple transformations under varying environmental conditions, resulting in a wide range of mineralogical and microstructural varieties (Cornell & Schwertmann, 2003). Due to the irregular distribution and uncommon concentration of these deposits, no universal classification exists for the geochemical anomalies they present. A useful approach to classify metallic deposits, based on geological environment, age, morphology, and mineralogy, distinguishes among iron formations, ironstones, and bog iron deposits (Stanton & Stanton, 1972).
The main iron-bearing minerals in these deposits include hematite (Fe₂O₃), magnetite (Fe₃O₄), limonite (2Fe₂O₃·3H₂O), siderite (FeCO₃), and pyrite (FeS₂) (Ussher et al., 2004). Geological evidence also indicates volcanic activity toward the southern sector of La Guajira, extending from the Triassic to the early Cretaceous. This magmatic event is represented by Los Tábanos Rhyodacite (Jt) (García-González et al., 2009).
The purpose of this study is to present and discuss data concerning the iron ores hosted in Los Tábanos Rhyodacite of San Juan del Cesar (La Guajira). Specifically, the study aims to characterize the occurrence of these ores and assess their petrological significance. We suggest that this mineralization may serve as an important case study, providing new insights into volcanic-hosted Fe mineral occurrences in Colombia.
Regional Geology
The SNSM is located in the northwestern part of Colombia, covering an area of 13.700 km2, with a maximum altitude of 5.900 m (Cardona et al., 2010). The SNSM forms a triangular massif of several lithologies. Within the framework of the Colombian Andes, it has been the subject of numerous studies (Campbell, 1965; Cardona et al., 2010a; Cardona et al., 2010; Cordani et al., 2005; Doolan & MacDonald, 1970; Herrera et al., 2008; Irving, 1975; Kellogg, 1984; MacDonald et al., 1971; Molina et al., 2006; Montes et al., 2010; Ordóñez et al., 2002; Radelli, 1962; Restrepo-Pace et al., 1997; C. Tschanz et al., 1974; C. M. Tschanz et al., 1969). Its geological evolution is the result of complex tectonic interactions among the Caribbean, South American, and Pacific plates during the Late Mesozoic to Cenozoic, which has triggered block rotation, basin formation and tectonic uplift in northern Colombia (Colletta et al., 1997; Kellogg, 1984; Montes et al., 2005, 2010). The massif is bounded by major faults and surrounded by Tertiary basins and the structural depression of the Cesar and Ranchería rivers. Its geological complexity reflects the effects of three major orogenic episodes and at least three periods of faulting. Two fault systems of particular importance in the tectonic evolution of northwestern South America affect the massif: the sinistral Santa Marta-Bucaramanga Fault and the dextral Oca Fault. Other major bounding structures include the he Oca fault (Irving, 1975), the Santa Marta-Bucaramanga fault (Campbell, 1965), the Cerrejón thrust sheet (Kellogg, 1984) and the Romeral suture (Montes et al., 2010). The massif is mainly composed of crystalline rocks and comprises three geotectonic provinces (Sierra Nevada, Sevilla, and Santa Marta), having a defined outboard younging pattern from east to west. The Geotectonic Province of Sierra Nevada (GPSN) includes ca. 1.0-1.2 Ga high-grade metamorphic rocks (granulites, gneisses, and amphibolites) that were affected during the Grenvillian orogenic event (Herrera et al., 2008; Restrepo-Pace et al., 1997). Jurassic plutons and vulcanites intrude and cover these metamorphic rocks. Minor Carboniferous and Late Mesozoic sedimentary sequences rest in unconformity towards the southeast (Tschanz et al., 1974; Tschanz et al., 1969). The Geotectonic Province of Sevilla (GPSE) corresponds to a polymetamorphic complex that includes Paleozoic gneisses and schists with Permian millonitized peraluminous granitoids (Cardona et al., 2010b; Herrera et al., 2008). The Geotectonic Province of Santa Marta (GPSM) comprises an inner sub-belt of Cretaceous imbricated metamorphic rocks (greenschists and amphibolites) and an outer sub-belt of Mesozoic amphibolites, greenschists and phyllites separated by the lower to middle Cenozoic Santa Marta Batholith (Doolan, 1970; Doolan & MacDonald, 1970).
Many tectonic models for the southern Caribbean postulate an accretion during the Late Cretaceous and/or Tertiary rocks of oceanic island arc of Mesozoic age. However, there are still numerous geological controversies yet to be resolved in relation to the studies performed by Radelli (1962) and Tschanz et al. (1969). In the study area, lithologies range from Precambrian to Paleocene ages. These rocks include the oldest geological unit corresponding to Los Mangos Granulite (PEm); the volcanic rocks of the Golero Rhyolites (Jg), the volcano-sedimentary rocks of Los Tábanos Rhyodacite and Los Clavos Ignimbrite (Jt and Jlc); the plutonic igneous rocks of the Pueblo Bello and Patillal Batholith in its three phases, two granitic phases and one quartzmonzonitic phase (Jpbp); the sedimentary rocks represented by Cretaceous limestones and shales, the Guatapurí and Corual Formations (Kcsi, Tjg, PTc). Figure 2 shows a generalized geological map of the study area. Los Tábanos Rhyodacite are the rocks of interest in the present study. They include rhyolites of alkaline and calc-alkaline affinity with high potassium, as well as alkaline trachytes and some dacites.

Source: Adapted and modified after (Colmenares et al., 2007)
Figure 2 Generalized geological map of the study area
Methodology
A total of four samples were collected from several outcrops representing different lithological units (to be specified), which were subsequently prepared for laboratory analysis. A preliminary macroscopic description was followed by a visual inspection using a Zeiss Stemi DV4 binocular stereomicroscope, in order to identify the textural and structural features of the materials. For detailed characterization, the following analytical techniques were employed: X-ray powder diffraction (XRPD), field emission gun-environmental scanning electron microscopy coupled with energy dispersive X-ray spectroscopy (FEG-ESEM/EDS), Fourier transform infrared spectroscopy with attenuated total reflection (FTIR-ATR), and X-ray fluorescence (XRF). The application of these methods has been previously described in detail in several works by Ríos and collaborators (Bonilla-Jaimes et al., 2016; Castellanos-Alarcón et al., 2016; Lizcano-Cabeza et al., 2015b, 2015a; Ríos et al., 2014) and is summarized below.
The bulk mineralogical composition was determined by XRPD using a BRUKER D8 ADVANCE diffractometer operating in Da Vinci geometry. The instrument was equipped with a Cu-Kα1 X-ray tube (λ = 1.5406 Å), a one-dimensional LynxEye detector (aperture angle of 2.93°), a 0.6 mm divergent slit, two soller slits (primary and secondary, 2.5° each), and a nickel filter. Samples were ground in an agate mortar to a particle size of <50 μm and mounted on polymethylmethacrylate (PMMA) holders using the front-filling technique. Data were collected under the following conditions: 40 kV and 30 mA, over the 2θ range of 3.5-70°, with a step size of 0.01526° (2θ) and a counting time of 1 s/step. Phase identification was performed using the PDF-2 crystallographic database (ICDD) and the Search-Match program. Unit-cell parameters, atomic positions, peak broadening factors, and phase concentrations were refined by the Rietveld method using the MDI RIQAS software.
The microtextural features and chemical composition of the samples were examined using back-scattered electron (BSE) imaging and EDS analysis on a FEI QUANTA 650 FEG-ESEM. Analytical conditions included a magnification of 100-800×, working distance (WD) of 10.2-10.8 mm, accelerating voltage of 20 kV, BSED-ETD detector, combined BSE-SE mode, horizontal field width (HFW) ranging from 2.98 mm to 373 μm, and an EDAX APOLO X EDS detector with a resolution of 126.1 eV (Mn Kα).
The structural characterization of functional groups was carried out by FTIR-ATR spectroscopy using a THERMO SCIENTIFIC IS50 spectrometer equipped with a diamond crystal, covering the spectral range of 400-4000 cm⁻¹.
Finally, the elemental composition of the samples was determined by wavelength-dispersive X-ray fluorescence (WDXRF) using a BRUKER S8 TIGER spectrometer (4 kW, Rh tube). The instrument was equipped with a flow sensor, scintillation crystals, interchangeable collimators, and specialized analyzers to minimize spectral overlap and improve detection sensitivity. For this analysis, samples were ground in an agate mortar, sieved to <38 μm, and then calcined. The resulting powders were remilled and placed in plastic sample holders (34 mm diameter). Spectral data and elemental concentrations were processed using the QUANT-EXPRESS software, based on the fundamental parameter method, achieving detection limits in the mg/kg range for mineral samples.
Results
Fieldwork occurrence
The iron ores of interest in this study occur within Los Tábanos Rhyodacite (Jt), represented by a volcanic-sedimentary sequence of rocks with a porphyritic texture and acid-to-intermediate composition, including andesitic, dacitic and trachytic tuffs, as well as lytic tuffs with blocks of volcano-sedimentary fragments. These rocks can be classified as highly fractured rhyolites and trachytes (Colmenares et al., 2007). Figure 3C illustrates locally flow structures, each of different color and texture, dominated by rocks of pink to very light purple color, with small phenocrysts of potassium feldspar). Figure 3A shows a pink to orange-pink color, with sporadic pink feldspar phenocrysts, as well as typical light pink, pink orange, very light purple or reddish gray aphanitic felsic rocks with several flow structures, often containing small feldspar phenocrysts. According to (Colmenares et al., 2007), based on a K-Ar dating in a crystal of sanidine of an obsidian 142 ± 6 Ma assigned an Upper Jurassic age for Los Tábanos Rhyodacite. Figure 3 illustrates the occurrence of outcrops of Los Tábanos Rhyodacite hosting iron ores. Field observations of the outcrops reveal rocks ranging from fresh to a moderately weathered. These rocks present their joint sets and other randomly orientated joints, which display similar characteristics with little variations. The openings of these discontinuity structures generally are classified as moderately spaced. The changes in the joint sets´ properties induced by weathering in Los Tábanos Rhyodacite are mainly of mineralogical type, developing late-stage oxidation halos (~5 mm) slightly stained with magnetite and hematite, which progressively decrease inward fresh rock.
Petrography
The host rock (rhyodacite) exhibits a hollocrystalline texture with an inequigranular distribution of grains, and color index of 35%. It is composed of plagioclase 50% quartz 25%, and feldspar 25%, with secondary minerals as hematite and magnetite (/25%). Hand specimens of Los Tábanos Rhyodacite show a freshly broken surface characterized by differential chemical weathering—probably mostly oxidation—progressing inward (Figure 4). This case accurately illustrates chemical weathering by oxidation, a process in which water and oxygen infiltrate the rock, reacting with ferrous minerals and generating oxides that modify its internal structure, color, and resistance. At the macroscopic scale, the analyzed samples of iron ores reveal a massive appearance of a thin skin mainly composed of magnetite-hematite of brown reddish color. The oxidation crust around the outside of the rock displays an inner white layer where the rhyodacite's minerals are beginning to break down and an outer red layer where magnetite and hematite are formed.
Scanning electron microscopy
Figure 5 presents SEM photomicrographs of the analyzed iron ores. Energy dispersive X-ray spectrometry (EDS) on scanning electron microscopy offer advanced analysis options for quantitative mineralogy and ore characterization (Salge et al., 2001). In this study, metal oxides were analyzed by measuring both the metal and oxygen peaks. However, distinguishing between hematite and magnetite requires careful evaluation of analytical consistency and precision. In addition, it is necessary to properly clean and polish samples and account for sample charging, sample homogeneity and other systematic sources of error (Konopka, 2012).
The EDS analyses proved the presence of iron oxide minerals. Figure 6 shows a typical EDS spectrum of the analyzed iron ores showed the presence of high concentration of Fe (64.35-70.80 wt%) and O (21.84-29.19 wt%). The peak of C is attributed to the carbon coating. The peaks of Si, Mg, Al, Ca and Mn (results not shown) represent the contribution of quartz and actinolite of the rhyodacite hosting rock. According to Paciornik and Gomes (2009), the discrimination of hematite and magnetite in backscattered electron images requires an extreme contrast in the image, which precludes the segmentation of other phases. On the other hand, even automated SEM systems with X-ray energy dispersion microanalysis (EDS) cannot differentiate between hematite and magnetite. Thus, the quantitative microstructural characterization of iron ores remains a challenge.
X-RAY diffraction
The analyzed samples contain 18.9-21.2% of magnetite, 21.6-24.2% of hematite, 9.4-8.7% of quartz and 6.4-8.6% of actinolite as crystalline mineral phases, and 27.3-41.9% of amorphous and other mineral phases. The X-ray diffraction pattern obtained by quantitative analysis is presented in Figure 7. The iron oxides (magnetite and hematite) have different crystal structures and, therefore, they have different patterns of X-ray diffraction. XRD pattern reveals that analyzed sample corresponds to a mixture of the two oxides.
Raman spectroscopy
Figure 8 illustrates the Raman spectra of the iron ores hosted in Los Tábanos Rhyodacite. The assignment of vibration band positions and the identification of the mineral phases were based on combined Raman data for the main vibration bands of iron oxides (e.g., De Faria et al., 1997; Shebanova & Lazor, 2003), which for Fe3O4 are 193 (week), 306 (week), 538 (week), 668 (strong); for α-Fe2O3 are 225 (strong), 247 (week), 299 (strong), 412 (strong), 497 (week), 613 (medium); for γ-Fe2O3 are 350 (strong), 500 (strong), 700 (strong). The characteristic peak positions of magnetite (Fe3O4), hematite (α-Fe2O3) and maghemite (γ-Fe2O3) were determined in the Raman region of interest in this investigation: 50-1500 cm-1. In the Raman spectra, the characteristic peaks of magnetite (Fe3O4) at 200, 463, 579, 661 and 669 cm-1 and hematite (α-Fe2O3) at 222, 242, 288, 405, 496 and 1300 cm-1 and maghemite (γ-Fe2O3) at 391 cm-1 are observed.
Chemical analysis
Iron ores are disseminated Los Tábanos Rhyodacite. The concentration of Fe2O3 ranges from 64.21 to 78.35% with sample YJ-1 and YJ-3 revealing the highest (78.35%) and lowest (64.21%) concentration of Fe2O3. This variation can be attributed to the dissemination of the iron ores in the rock. Al2O3 show a very low concentration (0.56-0.72%). The concentrations of CaO and MgO are 1.64-1.96% and 2.17-2.54%, respectively. K2O and TiO2 show a constant concentration (0.02%). MnO concentration ranges from 0.28 to 0.33. The concentration of SiO2 (15.40-30.17%) can be mainly related to the presence of quartz in the host rock. The major oxide results generally reflect the mineralogical composition of the iron ores disseminated in the rhyodacite (Table 1).
Origin of the iron ores
Weathering is the physical disintegration and chemical decomposition of igneous rocks into soil, dissolved chemical components, and solid chemical residues (Formoso, 2006). These rocks are most stable under the pressure, temperature, and environmental conditions present during their crystallization. When exposed to near-surface conditions, however, they are no longer in equilibrium and become progressively affected by weathering.
Los Tábanos Rhyodacite’s rocks crystallized under pressure and temperature conditions markedly different from those at the surface. Consequently, its primary minerals are unstable and undergo chemical reactions that alter their original composition (Figure 9). Physical weathering causes the disintegration of igneous rocks without chemical change. It can occur due to temperature, pressure, frost, etc. Los Tábanos Rhyodacite rocks were affected by physical weathering resulting in changing solid, large rock into smaller, movable unconsolidated debris, which was probably driven by changes in pressure and temperature conditions and organic activity. Chemical weathering is a gradual and ongoing process as the mineralogy of igneous rocks adjusts to the near surface environment, which occurs mainly at the surface of these rocks with the formation of secondary minerals develop from the original rock-forming minerals.
Los Tábanos Rhyodacite’s rocks were affected by oxidation, one of the more common and visible forms of chemical weathering. This reaction involves oxygen combining with minerals, particularly those containing iron to form iron oxide (rust). Therefore, the weathered region consists of rusty spots weaken these rocks and they break apart, which can be explained taking into account that iron oxide is much more fragile than iron, so. As a consequence of this process, the mineral compositions and the color of these rocks changed in the weathered region. Chemical weathering becomes more effective as the surface area of the rock increases (Sharma, 2014). In the case of rocks belonging to the Rhyodacite’s Los Tábanos Unit, it is likely that the alterations observed are mainly due to chemical reactions occurring on their exposed surfaces. However, it should be noted that the smaller the particle size, the greater the surface-to-volume ratio, which increases the specific area available for chemical interaction. According to Lapakko et al. (2006), the kinetics of these reactions are usually directly proportional to the surface area of the mineral exposed to the reactive medium.
In this study, a clear oxidation process is recognized in Los Tábanos Rhyodacite, driven by the interaction of iron and oxygen in the presence of water and the main products of oxidation are magnetite and hematite. It strongly affected these rocks, which supposes the presence of abundant water. Oxygen molecules (O2) collide with iron atoms of rock forming-minerals on the surface of the rock, which react to form iron oxide (rust) (Figure 10). Oxygen is a highly reactive gas that combines readily with other elements, such as iron. Rock-forming minerals of Los Tábanos Rhyodacite reacted directly with oxygen in the air to form iron oxides. However, the formation of rust after reaction of oxygen molecules and iron atoms was probably carried out in the presence of water, which is referred in this study as oxidation, a term for reactions in which iron atoms lose electrons and become more positively charged.
Evidently, Los Tábanos Rhyodacite´s rocks contained iron-rich minerals that were breakdown when exposed to the atmosphere, making them susceptible to chemical weathering. The main chemical reactions associated with chemical weathering are dissolution, hydrolysis, and oxidation (Ahluwalia, 2017). Oxidation occurred when oxygen reacted with iron-rich minerals to form iron oxide minerals. The absorption is usually from O2 dissolved in soil water and that present in the atmosphere and the oxidation is more effective in the presence of moisture and results in hydrated oxides (Rathinasamy & Saliha, 2017). According to these authors, during the oxidation processes, the following main chemical changes occur: (1) oxidation, solution and removal of the material, and (2) in situ transformation of rock-forming minerals into oxidized compounds.
4FeO (Ferrous oxide) + O2 → 2Fe2O3 (Ferric oxide)
4Fe3O4 (Magnetite) + O2 → 6Fe2O3 (Hematite)
2Fe2O3 (Hematite) + 3H2O → 2 Fe2O3.3H2O (Limonite)
As noted by Arikan and Aydin (2012), an increase in Fe2O3 content with increasing weathering indicates that oxidation plays a central role not only in the weathering of rocks of interest in the present study but also in the chemical transformation of iron minerals.
Rocks from Los Tábanos Rhyodacite are characterized by late-stage oxidation halos (~5 mm), which are common along joints. The exposure of these rocks to the weathering and its surface area affected its rate of weathering. The chemical weathering of these rocks occurs by the solution of Fe-rich mineral phases and tends to be concentrated at the rock surface or along joints, and a key element is the presence of water. Dissolution obviously occur, also acts as a medium for transporting acids, etc. Rainwater is acidic by nature due to atmospheric CO2, which reacts with water to form carbonic acid. However, a hydration process appear to have had little effect on the rock-forming minerals, as indicated by the absence of limonite; thus, there is no evidence of the conversion of iron oxides to iron hydroxides. The occurrence of iron oxides magnetite (Fe3O4) and hematite (Fe2O3) hosted in Los Tábanos Rhyodacite of San Juan del Cesar (La Guajira), SNSM, consist of more than 60% magnetite-hematite but those present a structural pattern similar to that of the lava mafic flows. Iron minerals, such as those reported in this study, can represent a rare case of an effusive Fe oxide liquid (Dare et al., 2015), however, fieldwork and petrographic evidences suggest that these probably are the result of the replacement of rhyodacite flows, with the development of magnetite-hematite pseudomorphs after rock-forming minerals.
Conclusion
The occurrence of iron minerals hosted in Los Tábanos Rhyodacite of San Juan del Cesar (La Guajira), Sierra Nevada of Santa Marta Massif, reveals that iron minerals occur in a freshly broken surface as a result of a chemical weathering process—probably mostly oxidation. The main products of oxidation are magnetite and hematite (α-Fe2O3) with traces of maghemite (γ-Fe2O3). The combination of analytical techniques, such as scanning electron microscopy, X-ray diffraction, and Raman spectroscopy is very useful for rapid, semiquantitative characterization of iron ore mineralogy. These techniques were well concurring with each other. Los Tábanos Rhyodacite´s rocks were affected by the oxidation as a result of the reaction of oxygen with minerals, particularly those containing iron to form iron oxide (rust).























