In the past, REE and high-field strength elements have been accepted as rather immobile elements, unaffected by hydrothermal alteration process and metamorphism. Recently it has been realized that REE could be mobilized under certain conditions. The detailed descriptions of REE mobility in the hydrothermal fluids and hydrothermally-altered rocks of some mineralized systems have been discussed elsewhere ( - ). The mobility of REE in the hydrothermal fluids is significantly controlled by low pH, high water/rock ratios and the ability of complexing ions such as , F−, Cl−, and (  ). In particular, they are preferentially complexed by Cl− and in acid conditions  . HREE are favorably concentrated in later products of hydrothermal systems as they form more stable complexes with some ligands and exist in solution longer than LREE ( ).
We examine the behavior of the trace and REE in the hydrothermal alteration of Cijulang high sulfidation epithermal gold deposit, West Java, Indonesia. The study on REE behavior was carried out in the alteration facies of silicic, advanced argillic, and argillic. The main objective of this study is to understand the physico-chemical conditions of hydrothermal fluids responsible for hypogene wall-rock alteration during the evolution of high-sulfidation epithermal system in the Cijulang area.
2. Geological Setting
The Islands of Java and Sumatra combined to form as a part of the Sunda-Banda magmatic arc which stretches from Aceh in the western part to Banda in the eastern part  (Figure 1). This magmatic arc, the longest of all magmatic arcs in Indonesia, has been created during Cenozoic due to the subduction of north-moving Indo-Australian Plate underneath the Southeast Asian Plate and tectonically active since Eocene ( ). The western segment which is located to the west of the central Java was developed on the continental crust whereas
Figure 1. Map showing the distribution of magmatic arcs in Indonesia .
the eastern segment which is situated to the east of the central Java was built upon the thinner island arc crust .
The Sunda-Banda magmatic arc constitutes volcanics and volcaniclastic rocks of Tertiary to Quaternary age. They are sometimes intruded by smaller plutonic masses with similar composition as volcanics. Numerous and still-active calc-alkaline magmatism occur throughout the magmatic arc . The volcanism of the entire Java arc has shifted from south to north over time ( ) which is evident by the E-W alignment of deeply dissected Miocene-Pliocene volcanic centers along the southern coast of Java and a parallel E-W alignment of juvenile and still active Quaternary volcanoes further north along the central Java. However, the magmatic affinities of the volcanics in the magmatic arc remained calc-alkaline with time until recent time ( ).
The Sunda-Banda Magmatic Arc served as a host for several hydrothermal deposits such as porphyry copper-gold, epithermal gold-silver, skarngold-copper, and volcanogenic massive sulfides. West Java is a well-known gold district in this magmatic arc and one of the most gold-producing areas in Indonesia. Present day operating gold mines are located in the three distinct regions of West Java: Bayah Dome Complex, Honje Igneous Complex and Southern Mountains. Recently several occurrences of epithermal gold and base metals mineralization have also been reported. Hence, much attention has been paid to these regions as they still stand as districts of future potential mineral resources. Cijulang is one of the interesting areas located in the Southern Mountains ( ) of West Java. Southern Mountains is regarded as the present-day forearc between the Quaternary volcanic chain and the Java Trench. It is comprised of volcanic and volcaniclastic rocks of Tertiary to Quaternary age.
Cijulang gold deposit was discovered by PT Aneka Tambang Tbk in 1996. Previous works ( ) indicated that gold mineralization in the Cijulang area shows typical features of high-sulfidation epithermal system by the presence of enargite-gold mineralization and accompanying acid sulfate alteration. Even though some researches and exploration works have been carried out in the Cijulang area, the behavior of REE and trace elements in different hydrothermal alteration facies and physico-chemical conditions of hydrothermal fluids during the evolution of high-sulfidation epithermal system have yet to be studied. Thus, the main objective of this study is to investigate the behavior of the REE and trace elements in the hydrothermal alteration facies of Cijulang area, West Java, Indonesia.
The local geology of the Cijulang area consists of andesite lava, lapilli tuff with minor crystal-lithic tuff, andesite and andesitic breccias (Figure 2). Andesite lava is regarded as a part of the Upper Oligocene to Middle Miocene Jampang Formation whereas lapilli tuff is considered to be a member rock of the Upper Miocene Koleberes Formation. Andesite and andesitic breccias belong to Mt. Kendeng Lava (Pleistocene). Geological structure in the study area is dominated by NE-SW and N-S trending strike-slip faults. The Cikahuripan Fault, the N-S trending dextral fault, played a significant role in the development of high sulfidation epithermal system in the research area (Figure 2).
Figure 2. Map showing the geology and hydrothermal alteration of the Cijulang area, West Java, Indonesia.
3. Hydrothermal Alteration
In Cijulang area, hydrothermal alteration is associated with pyrite-enargite-gold mineralization and covers over an area of approximately 2 km long and 1 km wide north-south trending elongate zone (Figure 2). Hypogene alteration primarily occurs in lapilli tuff and andesite lava. Based on field investigation, petrological, and mineralogical studies, four different hydrothermal alteration facies are identified: Silicic, advanced argillic, argillic, and propylitic. A characteristic alteration sequence and zonation of advanced argillic, argillic and propylitic outward from silicic core has resulted from the progressive cooling and neutralization of hot acidic magmatic fluids .
Silicic alteration occurs in the lapilli tuff (Figure 3(a)) and consists of both vuggy and massive quartz. This alteration is well-recognized by weather-resistant bodies of quartz. Vuggy zone comprises residual silica remaining after removal of its reactive components such as primary minerals and fragments
Figure 3. Various hydrothermal alteration types in Cijulang area. (a) vuggy silica alteration in lapilli tuff, (b) advanced argillic alteration in magmatic hydrothermal breccias, (c) argillic alteration in lapilli tuff, and (d) propylitic alteration in andesite lava.
of the host rock by leaching of highly acidic magmatic vapor . The vuggy silica core is generally formed by initial leaching by hot, extremely acidic magmatic fluids (temperature ~300˚C and pH < 2 ) with subsequent silica deposition and recrystallization in the vugs, open spaces and the matrix resulting massive quartz. Vugs are sometimes filled with clays (kaolinite, dickite, and pyrophyllite) (Figure 3(b)) and ores (pyrite and enargite). Silica zone often attains yellow or rusty color in places where oxidation of sulfides has taken place.
Advanced argillic alteration is the most widespread alteration and occurs in lapilli tuff and hydrothermal breccias (Figure 3(c)). In an outcrop, this alteration zone is strongly bleached and contains a variable amount of kaolinite, dickite, pyrophyllite, diaspore, illite and pyrite. Pyrophyllite, kaolinite and dickite are the dominant minerals in the outcrops whereas alunite is major constituents in the drill core sample (Figure 3(d)). Advanced argillic alteration occurs as massive replacement, fractures and vugs fillings (Figure 3(d)) in the host rock and as matrix in the magmatic hydrothermal breccia. The dominant mineral assemblages are quartz-pyrophyllite-kaolinite, quartz-dickite-pyrophyllite, quartz-kaolinite, quartz-alunite-dickite-pyrite, and quartz-alunite-kaolinite-pyrite.
Argillic alteration occurs as intermediate and discontinuous zones between advanced argillic and propylitic zones. It occurs in lapilli tuff and andesite lava. This zone is enriched in quartz and clay minerals such as kaolinite, illite, illite-smectite, smectite, muscovite, and chlorite. This alteration is quite pervasive and generally white and soft (Figure 3(e)). The common mineral assemblages are quartz-illite-kaolinite-pyrite and quartz-illite-montmorillonite-pyrite.
Propylitic alteration is the marginal facies of the hydrothermal alteration system and hosted by lapilli tuff and andesite lava (Figure 3(f)). It is characterized by quartz, chlorite, epidote, illite, smectite, pyrite, hematite, goethite, carbonate, and magnetite. In contrast to other alteration types, primary texture of the host rock is generally preserved. The common assemblage is quartz-chlorite-illite- epidote-calcite.
4. Analytical Techniques
A total of 75 samples, involving a spectrum of the least- to the most-altered rocks and fresh/un-altered rock from both surface and subsurface (drill cores section), were collected from several localities. Characteristics of the rocks such as lithology, mineralogy, and alteration of samples were examined in the field.
Mineralogical studies were carried out by optical microscopy in combination with X-ray (powder) diffraction (XRD) analysis. Bulk-rock and clay fraction samples (<2 µm) were determined using a Rigaku RINT-2100 Diffractometer. The analysis was performed at a scan speed of 2˚/min with CuKα radiation (λ = 1.542Å) from 2˚ to 65˚/2θ. Count data for the random powder mounts and clay fractions were collected from 2˚ to 65˚/2θ and 2˚ to 30˚/2θ respectively. Major and elemental compositions were determined by X-Ray Fluorescence (XRF) Spectrometer under pressed pellets using a RIGAKU RIX-3100 (Series VR 25006). The X-ray machine is operated at 50 kV and 50 mA, Scanning speed: automatic and 4˚/min for the determination of the major element composition. For quality control, the standard sample JA-3 was used as standard sample.
The mineral composition of alteration minerals was determined using SHIMADZU SS-550 Scanning Electron Microscope equipped with a Genesis-2000 EDS Spectrometer. Analytical conditions were 15 kV accelerating voltage, 10 mA beam current, and 3 μm beam diameter. XRD and XRF analyses were performed at the laboratory of Economic Geology, Department of Earth Resources Engineering, Kyushu University, Japan whereas SEM-EDS and EPMA analyses were carried out at The Center for Advanced Instrumental Analysis, Kyushu University, Japan. Analyses for trace elements and REE were carried out for 12 samples taken from three different alteration zones and un-altered precursor rock by instrumental neutron activation analysis (INAA), Inductively Coupled Plasma (ICP) and Inductively Coupled Plasma Mass Spectrometry (ICP-MS) methods at ACME Laboratories, Canada. Trace elements and REE elements were reported on a 0.1 g sample and concentrations were determined following a lithium metaborate/tetraborate fusion and dilute nitric acid digestion.
Detection limits for all major oxides are 0.01 wt%, except MnO, TiO2 and S which are 0.001 wt%. Detection limits for trace elements and REE are: Sc, Be, Co, Ni, Zn, Cu, Ga, Ge, Nb, Sn, W, 1 ppm; Rb, Sr, Y, Mo, 2 ppm; Ba, Se, 3 ppm; Zr, 4 ppm; V, Pb, As, Cr, Ir, 5 ppm; Cr, 20 ppm; La, Ce, Nd, Sm, Gd, Tb, Dy, Ho, Er, Yb, Ta, Tl, Th, U, Sc, 0.1 ppm; In, Hf, Sb, 0.2 ppm; Ag, 0.3 ppm; Bi, 0.4 ppm; Cd, Sb, Cs, As, Br, 0.5 ppm; Lu, 0.04 ppm; Pr, Eu, Tm, 0.05 ppm; Au, 2 ppb. Loss on ignition (LOI) was determined for all samples by measuring the weight difference after ignition at 1000˚C. The results of chemical analyses were presented in Table 1.
Table 1. Results of ICPMS analyses for different hydrothermally-altered and un-altered rocks from Cijulang area.
Note: L.O.I.: Loss on ignition; bdl: Below detection limit.
5.1. The Behavior of the REE and Other Trace Elements
Abundances of REE in studied samples from three different hydrothermal alteration facies (silicic, advanced argillic, and argillic) and un-altered rocks were normalized to chondrite values  and the resulting REE patterns of silicic (CSI-1 and CSI-2), advanced argillic (CAD-1, CAD-2, CAD-3 and CAD-4) and argillic facies rocks (CA-1, CA-2, CA-3 and CA-4) were compared with those of un-altered rocks (CUA-1 and CUA-2) which are considered to be the precursor rocks of hydrothermal alteration (Figures 4(a)-(c)). X-ray diffraction patterns of three different alteration facies are also shown (Figures 5(a)-(f)). The chondrite-normalized REE patterns of the un-altered volcaniclastic rocks indicate variable enrichment of LREEs with respect to HREEs and exhibit negative Eu anomalies which are typical characteristics of the average upper continental crust .
The ratios of (La/Sm)cn, (La/Yb)cn, (Tb/Yb)cn and (Eu/Eu*)cn were calculated for the different alteration facies and the precursor rocks (Table 1). Comparisons were made between precursor rocks and each alteration facies in order to understand the possible fractionation processes where REE could have been encountered during hydrothermal alteration (Figures 6(a)-(d)). (La/Yb)cn ratio was calculated in order to observe the enrichment of LREE over HREE. Fractionation among LREE was indicated by (La/Sm)cn and those between MREE and HREE were expressed by (Tb/Yb)cn. (Eu/Eu*)cn ratio was also calculated by using the equation (Eu/Eu* = Eu ). The geochemical behavior of REE and their ratios will be discussed in the following sections.
5.1.1. Silicic Alteration
As a result of strong acid leaching and mobilization of almost all oxides and elements, silicic alteration facies rocks are characterized by quartz (Figure 5(a)). REEs in this facies are strongly depleted if compared to the un-altered precursor rocks (Figure 4(a)). A strong depletion is more pronounced in MREE and HREE (Figure 4(a)). The degree of impoverishment of HREE is related to the degree of silicification as the most intensely-silicified rocks seem to be characterized by higher depletion of HREE.
Silicic facies is characterized by high ratio of (La/Sm)cn, moderately high ratio of (Tb/Yb)cn and low ratios of (La/Yb)cn and (Eu/Eu*)cn (Figure 6). A strong depletion of all REEs in this facies could be explained by the very low pH condition  and the abundances of complexing ions (Cl−, F−, and ) in the hydrothermal solutions (   ). The depletion of REEs might also be favored by the absence of secondary minerals that can host the REE in their structures.
5.1.2. Advanced Argillic Alteration
Based on the result of XRD analysis, advanced argillic facies rocks are characterized by alunite, pyrophyllite, kaolinite and dickite (Figures 5(b)-(d)). On the
Figure 4. Chondrite-normalized REE patterns the un-altered rocks (CUA-1 and CUA-2) and hydrothermally-altered rocks. (a) silicic alteration (CSI-1 and CSI-2), (b) advanced argillic alteration (CAD-1, CAD-2, CAD-3 and CAD-4), (c) argillic alteration (CA-1, CA-2, CA-3 and CA-4). Data in the figures are from Table 1; normalization values are taken from Sun and McDonough .
Figure 5. The result of XRD analyses for different alteration facies. (a) silicic alteration; ((b), (c) and (d)) advanced argillic alteration; (e) and (f) argillic alteration. Mineral abbreviation: Qz-quartz; Alu-alunite; Py-pyrite; Prl-pyrophyllite; Dck-dickite; Kln-kaolinite; Ilt-illite.
chondrite-normalized REE diagram, this facies rocks show LREE enrichment and its enrichment factor reaches about 110 times of chondritic level in the sample CAA1 (Figure 4(b)).
In contrast, both MREE and HREE concentrations are strongly depleted in this facies rocks. The (La/Yb)cn ratios are quite high, whereas (La/Sm)cn, (Tb/Yb)cn and (Eu/Eu*)cn ratios are moderately high (Figure 6). Fractionation
Figure 6. Variations in (La/Sm)cn, (La/Yb)cn, (Tb/Yb)cn ratios displaying REE fractionation in the hydrothermally-altered rock compared to un-altered equivalent original rocks. The (Eu/Eu*)cn ratio shows the variation of Eu anomaly. (Eu/Eu* = Eu × ).
of REE is quite obvious in this facies. Hydrothermal fluids responsible for the development of the advanced argillic alteration were derived from the partial neutralization and dilution of fluids which involved in the silicic alteration. The presence of alunite indicates a high activity of the sulfate ion in the hydrothermal solutions .
5.1.3. Argillic Alteration
The rocks of argillic alteration facies are composed of quartz; illite and kaolinite (Figure 5(e) and Figure 5(f)). REE patterns of this facies rocks do not show significant variation with respect to un-altered precursor rock (Figure 4(c)). This is also evident by the higher ratios of (La/Sm)cn in the studied samples (Figure 6).
5.1.4. Isocon Method
Isocon method of Grant  was applied in order to determine the mobility of elements in different alteration facies with respect to precursor rock, as the chemical analyses of an original and an altered rock cannot be directly compared to estimate gains and losses until a reference frame or assumption has been made . In this study, some components are assumed to be remained as immobile during alteration. Selected elements were plotted against TiO2 (Figure not shown) and those with a correlation coefficient r ≥ 0.6 were considered as immobile. La, Pr, Nd, Ho, Tm, Lu, and P2O5 were found to be immobile and are used as mass changes indicators. The geometry of isocon diagram is such that the elements below the straight isocon line are depleted whereas those above are enriched by the hydrothermal alteration (Figures 7(a)-(c)).
Figure 7. Isocon diagrams of three different alteration facies and histograms showing the percentage of gains and losses of selected elements during hydrothermal alteration as calculated from isocon method . (a) argillic, (b) an advanced argillic, and (c) silicic alteration. Values of +4 to +3 = a gain of 75% to 100%, +3 to +2 = 50% to 75%, +2 to +1 = 25% to 50%, and +1 to 0 = 0% to 25% respectively, and 0 = immobile elements. Values of 0 to −1 = losses of 0% to 25%, −1 to −2 = 25% to 50%, −2 to −3 = 50% to 75%, and −3 to −4 = 75% to 100%.
In the argillic rocks, oxides and element such as SiO2, K2O, Sn, Tl, Cs, Y, Pb, Zr, As, and S are enriched, however, the remaining components are depleted (Figure 7(a)). The advanced argillic rocks show strong enrichments of As, Sr, Ba, Pb, and SiO2 (Figure 7(b)). The concentration of other oxides and elements such as Sn, Th, P2O5, S and U are also slightly enriched whereas most major oxides (K2O, Al2O3, CaO, MnO, Na2O, MgO, and Fe2O3) and elements (Nd, Zr, Nb, Co, Cr, Y, Zn, and Sc) are depleted in the advanced argillic facies rocks. In the silicic alteration facies (Figure 7(c)), except for SiO2, As, Pb, Sn, U, Th and S, other oxides and elements are strongly removed from the rocks. The anomalous concentration of arsenic is observed in this facies rocks.
Based on the currently available data, the following summary and conclusions were made to interpret the development of hydrothermal system, and the REE behavior in different alteration facies of the Cijulang area:
A characteristic alteration sequence and zonation of an advanced argillic, argillic and propylitic alteration outward from the silicic core has developed from the progressive cooling and neutralization of hot acidic magmatic fluid with the host rocks. Silicic alteration facies is identified by the presence of residual and amorphous silica. All REE are strongly depleted in the silicic facies and a remarkable depletion of REE in this facies is probably caused by high acidity and abundance of ion and Cl− and F− complexing ions in the hydrothermal solutions. Advanced argillic alteration facies is mineralogically composed of variable amounts of quartz, pyrophyllite, alunite, dickite, kaolinite and diaspore. The enrichment of LREE along with depletion of HREE is observed in this facies. The immobile nature of LREE might be caused by the bounding of these elements into alunite structure. This confirms an important role of the occurrences of secondary minerals which control the behavior of REE in the hydrothermal paragenesis. The major constituents of argillic facies are quartz, illite, smectite, and pyrite. The relative immobility of REE in this facies is possibly related to the higher pH and lower water/rock ratio. The presence of phyllosilicates such as kaolinite, smectite and illite/sericite could also be a possible factor to fix the REE into their structures. Mass balance calculation indicates that REE are relatively mobile in the silicic facies whereas they are immobile in both advanced argillic and argillic facies.
The result of the study shows how REE behavior in the hydrothermal alteration paragenesis gives important information concerning to the physico-chemical conditions during the leaching processes. A strong variability of LREE, MREE and HREE concentrations of different alteration facies indicates variable degree of alteration and different fluid/rock ratios. It shows the selective nature of REE and trace elements even under strongly acidic conditions. The results also represent a useful tool to study HS epithermal ore deposits because it helps understand the mechanisms during the development of ‘main stage hypogene wall-rock alteration.
This research was financially supported by AUN/SEED-Net and JICA. The authors acknowledge the members of GeoMin Unit, PT Aneka Tambang Tbk for their support during field investigation in the Cijulang area.
 Lewis, A.J., Palmer, M.R., Sturchio, N.C. and Kemp, A.J. (1997) The Rare Earth Element Geochemistry of Acid-Sulphate and Acid-Sulphate-Chloride Geothermal Systems from Yellowstone National Park, Wyoming, USA. Geochimicaet Cosmochimica Acta, 61, 695-706.
 Alderton, D.H.M., Pearce, J.A. and Potts, P.J. (1980) Rare Earth Element Mobility during Granite Alteration: Evidence from Southeast England. Earth and Planetary Science Letters, 49, 149-165.
 Taylor, R.P. and Fryer, B.J. (1982) Rare Earth Element Geochemistry as an Aid to Interpreting Hydrothermal Ore Deposits. In: Evans, A.M., Ed., Metallization Associated with Acid Magmatism, John Wiley, New York, 57-365.
 Palacios, C.M., Hein, U.F. and Dulski, P. (1986) Behaviour of Rare Earth Elements during Hydrothermal Alteration at the Buena Esperanza Copper-Silver Deposit, Northern Chile. Earth and Planetary Science Letters, 80, 208-216.
 Lottermoser, B.G. (1990) Rare Earth Element and Heavy Metal Behaviour Associated with the Epithermal Gold Deposit on Lihir Island, Papua New Guinea. Journal of Volcanology and Geothermal Research, 40, 269-289.
 Arribas, A.J., Cunningham, C.G., Rytuba, J.J., Rye, R.O., Kelly, W.C., Podwisocki, M.H., McKee, E.H. and Tosdal, R.M. (1995) Geology, Geochronology, Fluid Inclusions, and Isotope Geochemistry of the Rodalquilar Gold Alunite Deposit, Spain. Economic Geology, 90, 795-822.
 Fulignati, P., Gioncada, A. and Sbrana, A. (1999) Rare Earth Element (REE) Behaviour in Alteration Facies of the Active Magmatic-Hydrothermal System of Volcano (Aeolin Island, Italy). Journal of Volcanology and Geothermal Research, 88, 325-342.
 Salvi, S., Fontan, F., Monchoux, P., Williams-Jones, A.E. and Moine, B. (2000) Hydrothermal Mobilization of High Field Strength Elements in Alkaline Igneous Systems: Evidence from the Tamazeght Complex (Morocco). Economic Geology, 95, 559-576.
 Jiang, N., Sun, S., Chu, X., Mizuta, T. and Ishiyama, D. (2003) Mobilization and Enrichment of High-Field Strength Elements during Late- and Post-Magmatic Processes in the Shuiquangousyenitic Complex, Northern China. Chemical Geology, 200, 117-128.
 Karakaya, N. (2009) REE and HFS Element Behavior in the Alteration Facies of the ErenlerDaği Volcanics (Konya, Turkey) and Kaolinite Occurrence. Journal of Geochemical Exploration, 101, 185-208. https://doi.org/10.1016/j.gexplo.2008.07.001
 Parsapoor, A., Khalili, M. and Mackizadeh, M.A. (2009) The Behaviour of Trace and Rare Earth Elements (REE) during Hydrothermal Alteration in the Rangan Area (Central Iran). Journal of Asian Earth Sciences, 34, 123-134. https://doi.org/10.1016/j.jseaes.2008.04.005
 Küpeli, ü. (2010) Trace and Rare-Earth Element Behavior during Alteration and Mineralization in the Attepe Iron Deposits (Feke-Adana, Southern Turkey). Journal of Geochemical Exploration, 105, 51-74. https://doi.org/10.1016/j.gexplo.2010.04.001
 Wood, S.A. (1990) The Aqueous Geochemistry of the Rare Earth Elements and Yttrium: Theoretical Prediction in Hydrothermal Solutions to 350 ºC at Saturation of Water Vapour Pressure. Chemical Geology, 88, 99-125.
 Haas, J.R., Shock, E.L. and Sassani, D.C. (1995) Rare Earth Elements in Hydrothermal Systems: Estimates of Standard Partial Modal Thermodynamic Properties of Aqueous Complexes of the Rare Earth Elements at High Pressures and Temperatures. Geochimicaet Cosmochimica Acta, 59, 4329-4350. https://doi.org/10.1016/0016-7037(95)00314-P
 Takahashi, Y., Tada, A. and Shimizu, H. (2004) Distribution of Pattern of Rare Earth Ions between Water and Montmorillonite and Its Relation to the Sorbed Species of the Ions. Analytical Sciences, 20, 1301-1306. https://doi.org/10.2116/analsci.20.1301
 Carlile, J.C. and Mitchell, A.H.G. (1994) Magmatic Arcs and Associated Gold Copper Mineralization in Indonesia, Journal of Geochemical Exploration, 50, 91-142.
 Rangin, C., Jolivet, L. and Pubellier, M. (1990) A Simple Model for the Tectonic Evolution of Southeast Asia and the Indonesian Region for the Past 43 m.y. Bulletin de la Société Géologique de France, 8, 889-905. https://doi.org/10.2113/gssgfbull.VI.6.889
 Whitford, D.J., Nicholls, I.A. and Taylor, S.R. (1979) Spatial Variations in the Geochemistry of Quaternary Lavas across the Sunda Arc in Java and Bali. Contributions to Mineralogy and Petrology, 70, 341-356. https://doi.org/10.1007/BF00375361
 Nicholls, I.A., Whitford, D.J., Harris, K.L. and Taylor, S.R. (1980) Variation in the Geochemistry of Mantle Source for Tholeitic and Calc-Alkaline Mafic Magma, Western Sunda Volcanic Arc, Indonesia. Chemical Geology, 30, 177-199.
 Soeria-Atmadja, R., Maury, R.C., Bellon, H., Pringgoprawiro, H., Polve, M. and Priadi, B. (1994) Tertiary Magmatic Belts in Java. Journal of Southeast Asian Earth Sciences, 9, 13-27.
 Verdiansyah, O., Bangun, P. and Rahmat, I. (2012) High-Sulfidation Epithermal Gold Occurrences in Cijulang Area, Garut, West Java. 41th IAGI Annual Convention and Exhibition, Yogyakarta, 17-20 September 2012.
 Tun, M.M., Warmada, I.W., Idrus, A., Harijoko, A., Al-Furqan, R. and Watanabe, K. (2015) Characteristics of Hydrothermal Alteration in Cijulang Area, West Java, Indonesia. Journal of Applied Geology, 7, 1-9.
 Sun, S.S. and McDonough, W.S. (1989) Chemical and Isotopic Systematic of Oceanic Basalts: Implications for Mantle Composition and Processes. In: Saunders, A.D. and Norry, M.J., Eds., Magmatism in the Ocean Basins, The Geological Society of London, London, Special Publication, 313-345. https://doi.org/10.1144/GSL.SP.1989.042.01.19