Abstract

Tinakula is the first seafloor massive sulfide deposit described in the Jean Charcot troughs and is the first such deposit described in the Solomon Islands—on land or the seabed. The deposit is hosted by mafic (basaltic-andesitic) volcaniclastic rocks within a series of cinder cones along a single eruptive fissure. Extensive mapping and sampling by remotely operated vehicle, together with shallow drilling, provide insights into deposit geology and especially hydrothermal processes operating in the shallow subsurface. On the seafloor, mostly inactive chimneys and mounds cover an area of ~77,000 m2 and are partially buried by volcaniclastic sand. Mineralization is characterized by abundant barite- and sulfide-rich chimneys that formed by low-temperature (<250°C) venting over ~5,600 years. Barite-rich samples have high SiO2, Pb, and Hg contents; the sulfide chimneys are dominated by low-Fe sphalerite and are high in Cd, Ge, Sb, and Ag. Few high-temperature chimneys, including zoned chalcopyrite-sphalerite samples and rare massive chalcopyrite, are rich in As, Mo, In, and Au (up to 9.26 ppm), locally as visible gold. Below the seafloor, the mineralization includes buried intervals of sulfide-rich talus with disseminated sulfides in volcaniclastic rocks consisting mainly of lapillistone with minor tuffaceous beds and autobreccias. The volcaniclastic rocks are intensely altered and variably cemented by anhydrite with crosscutting sulfate (± minor sulfide) veins. Fluid inclusions in anhydrite and sphalerite from the footwall (to 19.3 m below seafloor; m b.s.f.) have trapping temperatures of up to 298°C with salinities close to, but slightly higher than, that of seawater (2.8–4.5 wt % NaCl equiv). These temperatures are 10° to 20°C lower than the minimum temperature of boiling at this depth (1,070–1,204 m below sea level; m b.s.l.), suggesting that the highest-temperature fluids boiled below the seafloor. The alteration is distributed in broadly conformable zones, expressed in order of increasing depth and temperature as (1) montmorillonite/nontronite, (2) nontronite + corrensite, (3) illite/smectite + pyrite, (4) illite/smectite + chamosite, and (5) chamosite + corrensite. Zones of argillic alteration are distinguished from chloritic alteration by large positive mass changes in K2O (enriched in illite/smectite), MgO (enriched in chamosite and corrensite), and Fe2O3 (enriched in pyrite associated with illite/smectite alteration). The δ18O and δD values of clay minerals confirm increasing temperature with depth, from 124° to 256°C, and interaction with seawater-dominated hydrothermal fluids at high water/rock ratios. Leaching of the volcanic host rocks and thermochemical reduction of seawater sulfate are the primary sources of sulfur, with δ34S values of sulfides, from –0.8 to 3.4‰, and those of sulfate minerals close to seawater sulfate, from 19.3 to 22.5‰.

The mineralization and alteration at Tinakula are typical of a class of ancient massive sulfide deposits hosted mainly by permeable volcaniclastic rocks with broad, semiconformable alteration zones. Processes by which these deposits form have never been documented in modern seafloor massive sulfide systems, because they mostly develop below the seafloor. Our study shows how hydrothermal fluids can become focused within permeable rocks by progressive, low-temperature fluid circulation, leading to a large area (>150,000 m2) of alteration with reduced permeability close to the seafloor. In our model, overpressuring and fracturing of the sulfate- and clay-cemented volcaniclastic rocks produced the pathways for higher-temperature fluids to reach the seafloor, present now as sulfate-sulfide veins within the footwall. In the geologic record, the sulfate (anhydrite) is not preserved, leaving a broad zone of intense alteration with disseminated and stringer sulfides typical of this class of deposits.

Introduction

Volcanogenic massive sulfide (VMS) deposits have been classified according to host-rock lithology, bulk composition, and tectonic setting (e.g., Barrie and Hannington, 1999; Franklin et al., 2005). In extensional volcanic arc and back-arc environments, VMS deposits are typically hosted by volcanic successions including both clastic and coherent rock types (cf. Franklin et al., 2005; Hannington et al., 2005; Galley et al., 2007). Volcanic facies play a first-order control on the sizes and shapes of the deposits, the flow of hydrothermal fluids, and the nature of the footwall alteration zones. Coherent lavas that formed by effusive eruptions are relatively impermeable, and hydrothermal fluids tend to be focused along faults with minimal influx of cold seawater and form discordant pipe-like alteration zones (e.g., Ansel and Amulet deposits, Noranda district, Québec: Kerr and Gibson, 1993; Gibson and Galley, 2007). In contrast, volcaniclastic-dominated strata have high permeability and porosity, promoting lateral fluid flow, widespread mixing with seawater, and broadly stratiform alteration zones. There are numerous examples of ancient volcanicla-stic-hosted VMS deposits, including the Mattabi deposit in the Sturgeon Lake Caldera Complex, Canada (cf. Franklin et al., 1975; Morton and Franklin, 1987; Holk et al., 2008), the Horne and Bouchard-Hebert deposits in the Noranda mining district, Canada (cf. Gibson et al., 1993; Galley et al., 1995), the Boundary deposit in the Newfoundland Appalachians, Canada (cf. Piercey et al., 2014), and the Rosebery deposit in the Mount Read Volcanics belt of western Tasmania (cf. Large et al., 2001). In these deposits, the massive sulfide ore-bodies occur as stratiform lenses that are interpreted to form in whole or in part below the seafloor. Underlying the massive sulfide lenses, stringer mineralization is vertically restricted, occurring as semiconformable zones of disseminated sulfide mineralization. In the immediate footwall of the VMS deposits, zones of diffuse pipe-like alteration occur along synvolcanic faults, consisting of extreme leaching, silica mobility, and Fe metasomatism (e.g., Morton and Franklin, 1987). Typical alteration assemblages are sericite and chlorite, occurring as matrix cement and pervasive replacement of glass-rich clasts (e.g., Gifkins and Allen, 2001). These diffuse pipe-like zones merge outward into more widespread semiconformable alteration zones (e.g., Franklin et al., 1975). The geometry of these zones can be distinguished using oxygen isotope mapping, where increasing δ18O values of altered rocks are associated with decreasing temperatures of hydrothermal activity (e.g., King et al., 2000; Holk et al., 2008). Importantly, silicified horizons tend to form near the seafloor, acting as impermeable caps that isolate high-temperature circulation from shallow, low-temperature hydrothermal regimes (e.g., Franklin et al., 1981; Galley, 1993). This highlights the importance of early hydrothermal processes, such as silicification, in controlling the pathways of subsequent fluid flow (cf. Gibson et al., 1983; Skirrow and Franklin, 1994). While hydrothermal fluid flow regimes and associated alteration processes are well studied in the ancient rock record, these processes are largely unknown in modern settings due a lack of seafloor drilling to investigate the third dimension of these systems.

The Tinakula seafloor massive sulfide (SMS) deposit is the first identified modern analog of a mafic volcaniclastic-dominated deposit of the type described by Morton and Franklin (1987). It is located in the Jean Charcot troughs in the northern New Hebrides back arc, centered at 167°36′50″E/12°12′10″S (Fig. 1). This part of the Jean Charcot troughs is characterized by enhanced magmatism where regional structures, including the EW-trending Tikopia fracture zone, extend into the back arc (Fig. 1A). The transition from the arc to the back arc is marked by an abrupt change from tectonic extension associated with normal faulting to magmatic-tectonic extension where dike-fed lavas have erupted along linear fissures (Fig. 1B). Unlike hydrothermal systems in many deep back-arc basins, which are localized along mid-ocean ridge-type spreading centers (e.g., Hannington et al., 2005), hydrothermal activity at Tinakula is focused along a row of young mafic cinder cones on the flank of the Starfish stratovolcano (Figs. 1B, 2). Cinder cones are common volcanic landforms but have not been widely described on the modern seafloor. These features form from the buildup of volcanic tephra associated with Strombolian-style eruptions (i.e., fire-fountaining; cf. Cashman, 2004). Submarine pyroclastic deposits associated with fire-fountaining eruptions have been described in the ancient rock record (e.g., Mueller and White, 1992). This style of eruption may be more common in arc-proximal settings where mafic magmas have higher water contents compared to those in other tectonic settings (e.g., Gill et al., 1990).

Fig. 1.

Seafloor morphology of the Jean Charcot troughs. (A) Global marine altimetry from Sandwell et al. (2014), with back-arc troughs deeper than 2,800 m b.s.l. indicated by gray hatched pattern. (B) RESON SeaBat 8160 50-kHz multibeam bathymetry of the Starfish area (50 m resolution) collected by Neptune Minerals, Inc., showing transition from tectonic extension in the arc (forming NNW-trending normal faults, black lines) to tectonic-magmatic extension in the back arc (forming N-trending volcanic ridges, dotted black lines). Location of Tinakula deposit (Fig. 2) is shown, along with other sites of hydrothermal venting indicated by yellow stars (after Beaulieu and Szafranski, 2018). Globe inset in lower right corner shows location of research area in the South Pacific Ocean.

Fig. 1.

Seafloor morphology of the Jean Charcot troughs. (A) Global marine altimetry from Sandwell et al. (2014), with back-arc troughs deeper than 2,800 m b.s.l. indicated by gray hatched pattern. (B) RESON SeaBat 8160 50-kHz multibeam bathymetry of the Starfish area (50 m resolution) collected by Neptune Minerals, Inc., showing transition from tectonic extension in the arc (forming NNW-trending normal faults, black lines) to tectonic-magmatic extension in the back arc (forming N-trending volcanic ridges, dotted black lines). Location of Tinakula deposit (Fig. 2) is shown, along with other sites of hydrothermal venting indicated by yellow stars (after Beaulieu and Szafranski, 2018). Globe inset in lower right corner shows location of research area in the South Pacific Ocean.

Fig. 2.

Geology of the Tinakula deposit area interpreted from ROV dive videos and sampling overlain on hill- and slope-shaded RESON SeaBat 7125-ROV2 400-kHz multibeam bathymetry (inset), with locations of study sites and drill holes indicated.

Fig. 2.

Geology of the Tinakula deposit area interpreted from ROV dive videos and sampling overlain on hill- and slope-shaded RESON SeaBat 7125-ROV2 400-kHz multibeam bathymetry (inset), with locations of study sites and drill holes indicated.

In this paper, we present a comprehensive three-dimensional study of an actively forming massive sulfide deposit, with seafloor and subseafloor components, in a nascent backarc setting dominated by bimodal mafic volcaniclastic rocks. We investigate how features of the volcaniclastic substrate have influenced the mineralization and alteration at Tinakula, including the possible roles of boiling, fluid mixing, duration of hydrothermal activity, host-rock compositions, and magmatic contributions to the hydrothermal system.

Regional Geology and Exploration History

The Jean Charcot troughs in the northern New Hebrides (Vanuatu-Solomon Islands) constitute one of the youngest back arcs in the world, forming only ~4 m.y. ago in response to eastward subduction and rollback of the Indo-Australian plate beneath the Pacific plate (Fig. 1; Monjaret et al., 1991; Pelletier et al., 1998). Unlike the simple graben morphologies of other young back-arc rifts, such as the Coriolis troughs to the south (cf. Anderson et al., 2016), back-arc extension in the Jean Charcot troughs is characterized by a complex series of horst-and-graben structures that are widely spaced up to ~120 km from the arc front and have variable orientations (Fig. 1A; Pelletier et al., 1993; Maillet et al., 1995). Seafloor depths range from ~800 to 3,600 m below sea level (m b.s.l.), and individual grabens range from 20 to 65 km long and 5 to 15 km wide. The Jean Charcot troughs are bound to the south by the Hazel Holme fracture zone, a seismically active ridge that is considered to be an extensional zone related to opening of the North Fiji basin (Maillet et al., 1995). A second EW-trending extensional zone (Tikopia fracture zone; Fig. 1A) occurs north of the Hazel Holme fracture zone. Like the Hazel Holme fracture zone, the Tikopia fracture zone displays a slow-spreading trough morphology to depths of ~4,200 m b.s.l. east of Tikopia Island. To the west of this island, the presence of volcanic ridges and aligned seamounts suggests that the basement structure that controls the location of volcanic activity extends at least as far as the volcanic arc. Where the Tikopia fracture zone crosses the back-arc region, there is a notable bathymetric high, interpreted to be recent volcanic complexes that obscure the horst-and-graben morphology of the back arc (Fig. 1A; Récy et al., 1986). The Jean Charcot troughs are separated from the North Fiji basin by the West Tikopia Ridge and the Duff Ridge, which are interpreted to be part of the prerift fossil volcanic arc (Auzende et al., 1995). There is no consensus on the northern extent of the Jean Charcot troughs, although here we consider these to include the NW-SE–trending trough between the Duff Ridge and the Reef Islands. The current rate of extension of the troughs increases northward from 36 to 55 mm yr1 to 70 to 81 mm yr1 (Bergeot et al., 2009); however, this extension may be distributed over the entire width of the back arc, in which case individual troughs could be opening much more slowly.

Research cruises in the late 1980s to early 1990s provided most of the geochemical data that we now have for the Jean Charcot troughs (summarized by Maillet et al., 1995), but high-resolution (50 m or better) regional bathymetric data are lacking, and most of the region is unexplored, particularly in the north. The first indication of hydrothermal venting in the Jean Charcot troughs was a methane anomaly detected in the water column near the Starfish seamount during the KAIYO 89 cruise (Urabe et al., 1990). This anomaly was investigated on the SOLAVENTS cruise, which recovered fresh and altered rhyodacite samples by dredging and recorded a number of temperature anomalies in the water column near the seafloor (McConachy et al., 2002). Additional cruises in the area included the SEAPSO Leg 2 cruise (Récy et al., 1986), a GLORIA side-scan sonar survey (Johnson et al., 1992), and the SOPACMAPS cruise (Kroenke, 1995). More recent work has been aimed at exploration for SMS deposits. The Tinakula deposit was discovered by Neptune Minerals, Inc., in September 2011 during the ASHR Leg 2 cruise aboard the MV Dorado Discovery. It is located in the Temotu Province within the exclusive economic zone of the Solomon Islands, approximately 890 km east-southeast of Guadalcanal. The discovery was initially made by locating the source of a seafloor-hugging hydrothermal plume, followed by dives with a remotely operated vehicle (ROV), which located sulfide chimneys. Additional mapping, ROV grab sampling, and shallow drilling occurred during four subsequent (2011–2012) cruises.

Host Rocks, Volcanology, and Structure

The Tinakula SMS deposit is located at depths of 1,070 to 1,204 m b.s.l. on the flank of the Starfish stratovolcano (Fig. 1B). The geology of the Tinakula area has been interpreted from seafloor observations and sampling during ROV dives, drill core logging, and analysis of newly acquired high-resolution (20 cm) bathymetry. The bathymetry details reveal a strong structural control on volcanic and hydrothermal features (Fig. 2).

Volcanoes in the Tinakula area are composed of basaltic andesite (53.4–57.3 wt % SiO2) and dacite (63.0–67.9 wt % SiO2) (Anderson, 2018) and include composite cones, cinder cones, and the submarine equivalent of maars (craters without cones) (see App. Fig. A1). The morphologies of these edifices are broadly similar to those of subaerial volcanoes, but the edifices are consistently smaller than their subaerial equivalents. In the central part of the map area, the volcanic cones are aligned along a large (~1.6 km long) NS-trending fissure (Fig. 2, inset). The southern part of the fissure is dominated by a large composite cone with bimodal lava flows and blocky volcanic talus. This is followed northward by a 900-m-long row of cinder cones, on which the deposit is centered, and two small maar-like craters at the northern end of the fissure. The volcanic material making up the cinder cones consists of variably altered basaltic andesite with vertical and lateral facies variations over several meters, including changes in clast size, shape, composition, and degree of sorting (Figs. 3, 4; App. Figs. A2-A18). The breccia facies include the following: (1)

Fig. 3.

Volcanic and hydrothermal features of the footwall of the Tinakula deposit. (A) Variably clay altered, monomict, matrix-supported hyaloclastite breccia with flow foliation and local jigsaw-fit, sulfate-cemented with minor iron oxides (dark red) and chlorite alteration (green) (sample P003). (B) Intensely clay altered, monomict, clast-supported volcanic breccia (hyaloclastite), weakly sulfate-cemented (sample MA13-074). (C) Variably clay altered, polymict, clast-supported lapillistone, sulfate-cemented (sample P002). (D) Intensely clay altered, monomict, clast-supported lapillistone, sulfate-cemented and weakly graded (sample MA13-092). (E) Intensely clay altered, monomict, matrix-supported lapillistone with fine-grained beds and larger blocks, sulfate-cemented with late sulfate + sulfide veins (sample MA13-087). (F) Variably altered, polymict, matrix-supported proximal reworked volcanic breccia, silicified, with preserved primary porosity (sample MA13-023).

Fig. 3.

Volcanic and hydrothermal features of the footwall of the Tinakula deposit. (A) Variably clay altered, monomict, matrix-supported hyaloclastite breccia with flow foliation and local jigsaw-fit, sulfate-cemented with minor iron oxides (dark red) and chlorite alteration (green) (sample P003). (B) Intensely clay altered, monomict, clast-supported volcanic breccia (hyaloclastite), weakly sulfate-cemented (sample MA13-074). (C) Variably clay altered, polymict, clast-supported lapillistone, sulfate-cemented (sample P002). (D) Intensely clay altered, monomict, clast-supported lapillistone, sulfate-cemented and weakly graded (sample MA13-092). (E) Intensely clay altered, monomict, matrix-supported lapillistone with fine-grained beds and larger blocks, sulfate-cemented with late sulfate + sulfide veins (sample MA13-087). (F) Variably altered, polymict, matrix-supported proximal reworked volcanic breccia, silicified, with preserved primary porosity (sample MA13-023).

Fig. 4.

Microscopic volcanic and hydrothermal features of the footwall of the Tinakula deposit, shown by transmitted-light photomicrograph in plane-polarized light (TL-PPL) or cross-polarized light (TL-XPL) or reflected-light photomicrograph in plane-polarized light (RL-PPL). (A) TL-PPL photomicrograph of polymict clast types, including glass-rich volcanic clasts and tube pumice, cemented by amorphous silica (sample MA13-023). (B) TL-XPL photomicrograph of clay-altered scoriaceous clasts with intense clay alteration of clasts and matrix (sample MA13-068). (C) TL-PPL photomicrograph of volcanic clasts and glass shards surrounded by gypsum after anhydrite (sample MA13-040). (D) TL-PPL photomicrograph of tube pumice with late chlorite infill (sample MA13-053). (E) TL-PPL photomicrograph of intensely clay altered clast with multiple generations of clay infilling vesicles (sample MA13-044). (F) TL-XPL photomicrograph of clay-altered clast, with clay lining of vesicles and late sphalerite (red) (sample MA13-072). (G) RL-PPL photomicrograph of fine-grained pyrite disseminated throughout altered clast and late sulfides (pyrite, sphalerite, chalcopyrite) surrounded by anhydrite (blades) (sample MA13-045). (H) TL-XPL photomicrograph of anhydrite altered to gypsum (sample MA13-040). Abbreviations: Am-Si = amorphous silica, Anh = anhydrite, Ccp = chalcopyrite, Chl = chlorite, Gp = gypsum, GS = glass shard, Py = pyrite, Sp = sphalerite.

Fig. 4.

Microscopic volcanic and hydrothermal features of the footwall of the Tinakula deposit, shown by transmitted-light photomicrograph in plane-polarized light (TL-PPL) or cross-polarized light (TL-XPL) or reflected-light photomicrograph in plane-polarized light (RL-PPL). (A) TL-PPL photomicrograph of polymict clast types, including glass-rich volcanic clasts and tube pumice, cemented by amorphous silica (sample MA13-023). (B) TL-XPL photomicrograph of clay-altered scoriaceous clasts with intense clay alteration of clasts and matrix (sample MA13-068). (C) TL-PPL photomicrograph of volcanic clasts and glass shards surrounded by gypsum after anhydrite (sample MA13-040). (D) TL-PPL photomicrograph of tube pumice with late chlorite infill (sample MA13-053). (E) TL-PPL photomicrograph of intensely clay altered clast with multiple generations of clay infilling vesicles (sample MA13-044). (F) TL-XPL photomicrograph of clay-altered clast, with clay lining of vesicles and late sphalerite (red) (sample MA13-072). (G) RL-PPL photomicrograph of fine-grained pyrite disseminated throughout altered clast and late sulfides (pyrite, sphalerite, chalcopyrite) surrounded by anhydrite (blades) (sample MA13-045). (H) TL-XPL photomicrograph of anhydrite altered to gypsum (sample MA13-040). Abbreviations: Am-Si = amorphous silica, Anh = anhydrite, Ccp = chalcopyrite, Chl = chlorite, Gp = gypsum, GS = glass shard, Py = pyrite, Sp = sphalerite.

monomict breccia characterized by clasts with glassy curvi-planar margins, flow-banded phenocrysts, and local jigsaw-fit textures (Figs. 3A, B, 4B); (2) monomict breccia composed of vesicular volcanic clasts associated with variable amounts of juvenile clasts, glass shards, and crystal fragments and rare blocks or bombs (Figs. 3C-E, 4C); and (3) polymict breccia associated with mass-flow bedforms with rounded to subangular clasts, ranging from relatively unaltered to completely clay altered clasts (Fig. 3F). The predominant facies is massive to crudely bedded monomict breccia (lapillistone) with clasts containing up to 20% vesicles (uniformly distributed) and no evidence of welding. Basaltic-andesitic lava flows were the last products to be erupted along the row of cinder cones at Tinakula. These facies are consistent with cinder cones on land, which form via Strombolian-type eruptions, progressing from cone-building activity to lava flows during the waning stages as the magma loses most of its gas (Vergnoille and Mangan, 2000). Cinder cones of the type at Tinakula are generally considered to be monogenic, having formed during one eruptive cycle over a period of days to months (e.g., Wood, 1980).

Dark-colored mafic(?) volcaniclastic sediment covers most of the map area around the Tinakula deposit, interpreted to be fallout from nearby eruptions. The material ranges from coarse sand- to pebble-sized, with a few larger volcanic fragments or bombs (with elongate vesicles and quenched margins) and no fine ash-sized fraction. The sediment ranges in thickness from 1.5 to >8 m and is massive to crudely bedded with weak grading. Abundant ripple marks are observed in the sediment, with large-scale dunes (up to 7 m wide and 124 m long) visible in the bathymetric data.

Two large (55–87 m tall and 327–512 m wide) composite cones with dacite lava flows in the northwest corner of the mapped area (Fig. 2) appear to be unrelated to the central fissure. Large dacitic lava flows in the southwest corner of the map have a distinctly different morphology, with sinuous ridges that are characteristic of channelized flow. The relative ages of the flows are interpreted from ROV dive observations of sediment thickness, with higher degrees of sediment associated with older relative ages. The ages of these flows relative to other eruptions in the area are unknown.

Deposit Geology

The Tinakula SMS deposit consists of a series of over 900 chimneys and mound complexes distributed over a strike length of ~1,050 m along the row of cinder cones and composite cones (Fig. 2). The chimneys and mounds are located mainly at the tops and along the flanks of the cinder cones, as well as inside the summit craters (Fig. 5). They also occur in the center of each of the maar-like craters in the north. West of the row of cinder cones, over a strike length of ~400 m, the chimneys and mounds are not spatially associated with observable volcanic features at the seafloor but instead appear to occur along two partially buried northeast-southwest and northwest-southeast structures (between sites F and D and between sites F and E; Fig. 2).

Fig. 5.

3-D close-up of the bathymetry of study sites in the Tinakula area (RESON SeaBat 7125-ROV2 400-kHz multibeam bathymetry; 0.2-m resolution), with chimney complexes and mounds (purple) and locations of drill holes (red circles) examined in this study. (A) Site A, looking northwest. (B) Site B, looking northwest. (C) Sites C and D, looking northeast. (D) Sites C and D, looking northwest. (E) Site E, looking northwest. (F) Site F, looking northeast. Locations of cross sections (Figs. 10, 11) are shown by thick black lines.

Fig. 5.

3-D close-up of the bathymetry of study sites in the Tinakula area (RESON SeaBat 7125-ROV2 400-kHz multibeam bathymetry; 0.2-m resolution), with chimney complexes and mounds (purple) and locations of drill holes (red circles) examined in this study. (A) Site A, looking northwest. (B) Site B, looking northwest. (C) Sites C and D, looking northeast. (D) Sites C and D, looking northwest. (E) Site E, looking northwest. (F) Site F, looking northeast. Locations of cross sections (Figs. 10, 11) are shown by thick black lines.

The chimneys at Tinakula are <1 to 15 m high (typically 2–8 m) and include typical multispired sulfide structures with central vent orifices (Fig. 6A, B) and, more commonly, chimneys with irregular morphologies (bulbous, drooping, with smooth or popcorn-textured outer surfaces; Fig. 6C, D). The bulbous chimneys lack central orifices and have a morphology similar to those of beehives or flange-type chimneys (e.g., Koski et al., 1994; Zhu et al., 2007). The chimneys are commonly surrounded by sulfide talus and grow on top of low-relief hydrothermal mounds that typically rise ~2 to 3 m from the seafloor. The extent of these mounds is largely unknown due to burial by volcaniclastic sediments, but they appear to make up the bulk of the mineralization at the seafloor, accounting for about two-thirds of the total amount of identified massive sulfide. Individual chimney mound complexes can be traced for up to 70 m on the seafloor. Where exposed, the mound material is distinctly laminated with layers of barite and sulfide minerals up to 0.5 m thick, and they have unexpectedly smooth surfaces (Fig. 6C). The layered barite-sulfide mounds are most similar to features described by Hein et al. (2014) from the East Diamante caldera in the Mariana arc. The sulfide deposits are mostly inactive, and the surrounding volcaniclastic sediment is widely indurated and locally mineralized (Fig. 6D). Active venting is limited to a few sites with shimmering water at the base of chimneys or along fractures in the indurated sediments, associated with bacterial mats (Fig. 6A). Venting of white to gray smoke from chimney spires is rare (Fig. 6B). Fe-Mn oxides are widespread along fractures in the seafloor, extending up to 100 m away from the chimney complexes and mounds (Fig. 2). Abundant <1-m-tall Fe and Mn oxyhydroxide chimneys and silica-rich chimneys are associated with bacterial mats (Fig. 6E) and occur along the fractures (Fig. 6F).

Fig. 6.

Chimney morphologies encountered during ROV dive operations. (A) Partly oxidized inactive spired chimneys with bacterial mat, field of view ~15 m. (B) Weakly venting spire with minor oxidation, field of view ~3 m. (C) Inactive bulbous chimney on top of layered barite-sulfide mound, field of view ~5 m. (D) Inactive bulbous chimney partially buried by volcaniclastic sediment (indurated), field of view ~3 m. (E) Siliceous iron oxide chimney with small worm casts, field of view ~1 m. (F) Iron and manganese oxyhydroxides forming dendritic oxide chimneys in fractures within indurated volcaniclastic sediment, field of view ~5 m.

Fig. 6.

Chimney morphologies encountered during ROV dive operations. (A) Partly oxidized inactive spired chimneys with bacterial mat, field of view ~15 m. (B) Weakly venting spire with minor oxidation, field of view ~3 m. (C) Inactive bulbous chimney on top of layered barite-sulfide mound, field of view ~5 m. (D) Inactive bulbous chimney partially buried by volcaniclastic sediment (indurated), field of view ~3 m. (E) Siliceous iron oxide chimney with small worm casts, field of view ~1 m. (F) Iron and manganese oxyhydroxides forming dendritic oxide chimneys in fractures within indurated volcaniclastic sediment, field of view ~5 m.

The surface expression of hydrothermal precipitates covers ~77,000 m2, including ~13,800 m2 of active chimneys and mounds, 51,800 m2 of inactive chimneys and mounds, and 9,300 m2 of partially buried chimneys and mounds (Fig. 2). Areas of low-temperature venting, dominated by Fe and Mn oxides, occur over an additional 130,000 m2 (91,000 m2 of Fe oxides and 38,000 m2 of Mn oxides), although these areas commonly overlap those of the chimneys and mounds. The area of subseafloor sulfide mineralization is not completely known, so the full size of the deposit could be much larger. The areal extent of this deposit is significant compared to others discovered along mid-ocean ridges and along arcs and back-arc basins; fewer than 10% have dimensions larger than 5,000 m2 (Hannington et al., 2010). Modern sulfide deposits with comparable sizes (by area) include the Solwara 1 deposit in the Eastern Manus basin (90,000 m2) and the Mir zone, Ashadze-1, and Ashadze-2 deposits on the Mid-Atlantic Ridge (50,000 m2 each) (Hannington et al., 2010, 2011). It should be emphasized, however, that the large areal footprint of the Tinakula deposit may not necessarily correlate with a large tonnage of massive sulfide.

Methods

Ship-based multibeam bathymetry operating at 50 kHz (RESON SeaBat 8160) was collected over 685.6 km2 in the Starfish-Tinakula area, and the data were gridded with a cell size of 50 m. Higher-resolution bathymetry was collected over ~2.2 km2 that includes the Tinakula deposit, using an ROV-mounted multibeam system (RESON SeaBat 7125-ROV2 at 400 kHz with inertial navigation) and gridded with a cell size of 0.2 m. Bottom observations and grab samples were collected during 13 ROV dives. Shallow seafloor drilling, including push-coring and rotary diamond drilling, was carried out using the Benthic Geotech portable remotely operated drill PROD-3. Approximately 270 m of drill core was logged and sampled to establish between-hole correlations along the length of the NS-trending eruptive fissure (sites A-D; Fig. 2; App. Figs. A2-A18). Mineralogical and chemical data are reported here on 41 ROV grab samples and 21 drill core samples of chimney and mound material and on 61 drill core samples of clay-altered volcaniclastic rocks. Polished thin sections of the sulfide-sulfate chimney samples and epoxy-impregnated clay-altered volcaniclastic samples were prepared by Dettmar Dissection Technology (Bochum, Germany). Clay-altered samples were cut using an industry-standard water-free technique that employs kerosene in order to prevent clay swelling, slaking, and dissolution of anhydrite. Thin sections were analyzed by transmitted- and reflected-light microscopy.

Whole-rock geochemistry was determined on 45 bulk samples of sulfide-sulfate material and 47 samples of fine (<63 μm) clay-altered volcaniclastic material at Activation Laboratories (Ontario, Canada). For the sulfide samples, a sodium peroxide fusion was used; major elements (Al, Ca, K, Mg, P, Si, Ti) were analyzed by inductively coupled plasma-optical emission spectrometry (ICP-OES) and trace elements (B, Ba, Be, Bi, Cd, Co, Cs, Cu, Ga, Ge, In, Li, Mn, Mo, Nb, Ni, Pb, Rb, Sn, Sr, Ta, Th, U, V, W, Y, Zn, and rare earth elements [REEs]) by inductively coupled plasma-mass spectrometry (ICP-MS). For the clay-altered samples, a lithium metaborate/tetraborate fusion was used; the major elements (SiO2, Al2O3, Fe2O3(T), MnO, MgO, CaO, Na2O, K2O, TiO2, P2O5) were analyzed by ICP-OES, and trace elements (Ba, Be, Bi, Co, Cs, Ga, Ge, Hf, In, Mo, Nb, Rb, Sr, Ta, Th, Tl, U, V, W, Y, Zr, REEs) by ICP-MS, with additional elements (Ag, Cd, Cu, Ni, Pb, Zn) analyzed by ICP-OES following a total acid digestion technique. In both sample types, cold-vapor atomic absorption spectrometry (AAS) was used to determine Hg; infrared (IR) detection was used to determine total S. Instrumental neutron activation analysis (INAA) was used to determine Au, Ag, As, Ba, Br, Cr, Fe, Hf, Ir, Lu, Na, Sb, Sc, and Zn in the sulfides and Au, As, Br, Cr, Ir, Sc, Se, and Sb in the clay-altered samples. Duplicates and standard reference materials were analyzed after every 10 to 15 samples and indicated a precision of better than 10% and accuracy of better than 5% for most elements in samples where the analytical values are well above the limit of detection, with values degrading rapidly as the limits of detection are approached.

Sulfide mineral compositions were determined in 16 sulfide-sulfate samples using a JEOL-JXA 8200 electron probe microanalyzer (EPMA) equipped with five wavelength dispersive spectrometers and one energy dispersive detector at the GEOMAR Helmholtz Centre for Ocean Research in Kiel. A 15-keV accelerating potential and 50-nA probe (cup) current were used for the analyses. Standards used were chalco-pyrite (FeKα, SKα, CuKα), rhodonite (MnKα) gallium arsenide (GaLα, AsKβ), sphalerite (ZnKα), silver telluride (AgLα), indium antimonide (InLα, SbLα), cadmium sulfide (CdLα), gold (AuLα), cinnabar (HgMα), and galena (PbMα).

Alteration mineralogy was determined by X-ray diffraction (XRD) analysis of bulk samples of the altered footwall rocks. Fine (<63 μm) and coarse (>63 μm) fractions were separated by wet sieving, followed by separation into silt (2–63 μm) and clay (<2 μm) by settling in cylinders of standing water, following the technique of Moore and Reynolds (1997). The mineralogy of the <2-μm clay fraction from 58 samples was determined using a Philips PW 1710 X-ray diffractometer, equipped with a Co tube and an automatic divergence slit and monochromator, located at GEOMAR. Operating conditions were 40 kV and 35 mA. Oriented mounts were prepared by wet suspension followed by air drying and were measured with a 2θ scanning angle of 3° to 40° at a scan rate of 1 s per 0.01° step. Samples were then saturated with ethylene glycol following a standard vaporization technique and reanalyzed. MacDiff v.4.2.6 software was used for mineral identification and display of the XRD data.

Short-wave infrared (SWIR) reflectance spectroscopy was also used to identify alteration assemblages and compositional variations in the altered footwall rocks. Spectral measurements of fine (<63 μm) fractions were conducted using a Fieldspec® Pro spectroradiometer at the University of Ottawa, with a spectral bandwidth of 3 to 10 nm and a scan time of 100 ms. Duplicate measurements were made for every sample; white reference measurements were made every fifth sample. The data were normalized, and mineral phases were automatically identified using The Spectral Geologist (TSG) 7 software and a built-in reference library (e.g., Jinfeng et al., 2014). The spectra were then manually assessed to confirm mineral identifications based on the shape and wavelength minima of the absorption features.

Five samples consisting of hydrothermally altered and sulfate-cemented volcaniclastic material were selected from drill hole TKAD001 for fluid inclusion studies. Doubly polished thick sections (~130 μm) were prepared at Dettmar Dissection Technology using a water-free technique. Fluid inclusion petrography in barite, anhydrite, and sphalerite is described in detail in Anders (2016). Microthermometric measurements were made on a Linkam™ THMS 600 heating and freezing stage at GEOMAR. Salinities were calculated from freezing-point depression values following Bodnar and Vityk (1994). Homogenization temperatures were corrected for pressure at the sample depth (1,107.61 m b.s.l.; 108.65 bar) following Zhang and Frantz (1987). Accuracy was determined by measuring phase transitions in synthetic fluid inclusions (Syn Flinc®); replicate measurements of homogenization and freezing-point depression temperatures were performed every five to 10 analyses, resulting in estimated errors of less than ±0.2°C for salinity measurements and ±4°C for homogenization temperature measurements. The critical points, bulk molar volumes, densities, and isochores were calculated using the program FLINCOR (Brown, 1989).

Oxygen and hydrogen isotope analyses were performed on 14 samples of monomineralic clays (<2 μm), free from Fe and Mn oxides as determined by XRD. The analyses were carried out by the Queen’s Facility for Isotope Research (Queen’s University, Ontario, Canada). Oxygen was extracted from 5-mg samples at 550° to 600°C using the BrF5 technique of Clayton and Mayeda (1963). Samples were analyzed via dual inlet on a Thermo-Finnigan DeltaPlus XP isotope ratio mass spectrometer (IRMS). For hydrogen isotope measurements, samples were weighed into silver capsules, degassed for 1 h at 100°C, and then crushed and loaded into a zero blank autosampler. Hydrogen isotope compositions were measured using a Finnigan thermos-combustion elemental analyzer (EA) coupled to a Thermo-Finnigan DeltaPlus XP continuous-flow IRMS. All values are expressed in the conventional δ notation as per mil (‰) deviation from V-SMOW (Vienna-Standard Mean Ocean Water), with analytical precision and accuracy for δ18O values of ±0.1‰ (1σ) and for δD of ±3‰ (1σ) and ± 1.5‰ (1σ), respectively.

Sulfur isotope measurements were conducted on 57 handpicked mineral separates, including anhydrite, gypsum, barite, chalcopyrite, and sphalerite, at the Westfälische Wilhelms-Universität Münster, Germany (described in detail by Wienkenjohann, 2015). For 35 pure monomineralic samples (as defined optically), the sulfur isotope composition was determined directly without prior chemical treatment, using a Carlo Erba EA interfaced to a ThermoFinnigan Delta Plus mass spectrometer (EA-IRMS). Repeat analyses showed a precision of better than ±0.3‰ (1σ), and standard measurements (IAEA-S1, S2, S3, NBS127, Ag2S-lab) gave an accuracy of better than ±0.2‰ (1σ). For 22 samples of intergrown sulfide-sulfate mineral assemblages, a mineral-specific, wet-chemical pretreatment protocol was applied, yielding four different sulfur phases. Briefly, anhydrite and gypsum were dissolved using a 10% NaCl-solution for 24 h, and dissolved sulfate was precipitated as BaSO4 using standard procedures (e.g., Dogramaci et al., 2001). BaSO4 precipitates were then converted into silver sulfide using Thode solution (a mixture of HI, HCl, and H3PO2; Thode and Monster, 1961). Sphalerite sulfur was liberated as hydrogen sulfide (acid volatile sulfide [AVS]) with HCl (25%) for 1 h at 85°C. Subsequently, disulfides such as chalcopyrite and pyrite were extracted using chromous-chloride solution (chromium-reducible sulfide [CRS] fraction, cf. Canfield et al., 1986). Each time, resulting hydrogen sulfide was precipitated first as zinc sulfide (using zinc acetate) and ultimately as silver sulfide (using silver nitrate) as described in Canfield et al. (1986) and Rice et al. (1993). Following the quantitative removal of all sulfide minerals using acidic (HCl) CrCl2 solution (cf. Canfield et al., 1986), acid-insoluble sulfate (barite) was liberated as hydrogen sulfide, following Thode and Monster (1961), and ultimately precipitated as silver sulfide. All silver sulfide precipitates, representing water-soluble sulfates (WSS) such as anhydrite and gypsum (Ag2SWSS), acid-insoluble sulfates (AISs) such as barite (Ag2SAIS), AVSs (i.e., sphalerite), and CRSs (i.e., chalcopyrite and pyrite), were subjected to fluorination IRMS, following the procedure outlined by Ono et al. (2006), using a ThermoScientific Mat 253 mass spectrometer. Results are reported in the standard delta notation as per mil difference relative to Vienna-Canyon Diablo Troilite (V-CDT). Precision and accuracy were better than ±0.03‰ (1σ).

The 226Ra/Ba dating of hydrothermal barite from 23 sulfidesulfate samples was carried out following the methods of Jamieson et al. (2013), modified after Ditchburn et al. (2004). Samples were crushed and powdered, with one split used to determine Ba concentrations by INAA (described above) and another split used to determine 226Ra activity via gamma spectrometry at the Elemental Tracer Laboratory at the University of Ottawa. The daughter product of 226Ra is 222Rn, with a halflife of 3.8 d. A daughter isotope from the decay of 222Rn, 214Pb, was measured as a proxy to avoid peak overlap with 235U. Prior to measurements, the sample tubes were sealed for 21 d to prevent the escape of 222Rn gas and obtain secular equilibrium between 226Ra and 214Pb. Measurements were made using an XLB-GWL-SV digital gamma spectrometer (DSPec, Ortec®) containing a high-purity germanium crystal, 54.7 × 66.2 mm, with an active well depth of 40 mm. Counting times ranged from 2 to 12 h, depending on the activity of 226Ra in each sample. Spectral signatures were calibrated using IAEA-312 reference material and analyzed using a custom program developed and provided by Peter Appleby (University of Liverpool) to correct for background, detector efficiency, sample volumes, and matrix effects. Ages were determined relative to a zero-age end member selected from one of the actively venting chimneys and assumed to represent the initial 226Ra/Ba ratio (1,022 ± 68 Bq/kg·wt %; sample GMA13-024). All ages are considered minimum values, as incorporation of young barite in the samples cannot be ruled out (e.g., Ditchburn and de Ronde, 2017).

Drill Sites

Shallow seafloor drilling was focused along the main NS-trending fissure (sites A-E), targeting areas with accumulations of chimneys and mounds and spatially associated with buried circular structures that may be manifestations of the cinder cones, with an additional site to the west of this structure (site F) (Fig. 2). Three subseafloor rock types were encountered in most drill holes, in order of increasing depth:

  • (1) unconsolidated volcaniclastic sediment (cover sequence),

  • (2) mineralized intervals consisting of mixed chimney talus, and (3) intensely clay altered volcaniclastic material (footwall sequence). In addition, several drill holes intersected intervals of oxidized ferruginous mud associated with rounded to subrounded chimney fragments, interpreted to be paleoseafloor horizons.

Site A

Site A is a local topographic high at the northern end of the row of craters. It rises ~10 m above background sediments over an area of ~100 × 80 m (Figs. 2, 5A). Chimneys and mounds are distributed around what may be a buried crater, partially covered by volcaniclastic material. The volcaniclastic cover sequence (up to 13.8 m thick) is underlain by 0.7 to 2.6 m of chimney talus, followed by intensely clay altered monomictic lapillistone that is crudely bedded with both matrix-supported and clast-supported intervals, extending to the bottom of the drill holes. The northernmost drill hole, TKAD001, is located next to an active chimney complex and reaches a depth of 22.3 m below seafloor (b.s.f.; App. Fig. A2). The upper ~14 m consists of variably altered volcaniclastic cover material with three discrete oxidized intervals, followed by altered lapillistone. The upper part of the clay-altered lapillistone sequence is unconsolidated to poorly consolidated, with increasing cementation by anhydrite downhole, and is cut by a few anhydrite veins (2–5 mm wide). Disseminated sulfides (pyrite, sphalerite, and minor chalcopyrite, galena, covellite, and bornite) occur with anhydrite that has been partly hydrated to gypsum, both in the cement and in the vein fill. To the south of this drill hole, TKAD015 and TKAD027 were drilled in close proximity (~3 m apart) to a maximum depth of 7.3 m b.s.f., next to a partially buried chimney complex. The volcaniclastic cover is <1 m thick and contains fragments of chimney material, likely sourced by mass wasting from nearby (unknown) chimneys (App. Figs. A3, A4). The mineralized interval is a ~1.5- to 2.6-m-thick layer of chimney talus. As in TKAD001, intensely clay altered and poorly sorted lapillistone occurs below the chimney talus in both drill holes. TKAD028, located ~20 m east of TKAD015/027, intersected 0.5 m of volcaniclastic cover material, followed at depth by a narrow (0.7 m) sulfide interval dominated by barite and sphalerite at the top and becoming pyrite rich downhole (App. Fig. A5). This mineralized interval appears to correlate between the drill holes. The clay-altered lapillistone at the bottom of the hole is similar but better sorted than that found in TKAD015/027. Approximately 30 m to the east, TKAD029 encountered only the variably altered volcaniclastic cover material to a depth of 7.8 m b.s.f. (App. Fig. A6). This drill hole is located at the edge of the low-relief cinder cone where the volcaniclastic cover is thickest.

Site B

Site B is another crater-like feature (70 × 80 m) on a cone at the north end of the row of cinder cones (Figs. 2, 5B). Chimneys and mounds occur along the crater walls and on the outer flanks of the cone and are variably covered by volcaniclastic material. Here, massive sulfide chimneys occur on top of chimney talus. At the western edge of the crater, TKAD016 intersected thin (0.2 m) volcaniclastic cover, followed by a layer of pyritic chimney talus ~0.4 m thick and then poorly sorted and crudely bedded monomict, matrix-supported, and intensely clay altered lapillistone (App. Fig. A7). The upper part of the lapillistone is cemented mainly by gypsum with a few large crystals up to 7 mm long; the lower part is mostly unconsolidated with patchy anhydrite cement that has been partially hydrated to gypsum, a few crosscutting veinlets of anhydrite, weak silicification, and disseminated sulfide minerals (pyrite with minor sphalerite near the top and minor to trace chalcopyrite near the bottom). At the southern wall of the crater, TKAD026 intersected massive sphalerite-barite chimney material to a depth of 1.4 m b.s.f. (App. Fig. A8). The substrate of lapillistone contains minor patchy barite cement in the upper part and anhydrite cement in the lower part, with gray-white clay alteration that generally decreases downhole. At the center of the crater, the mineralized interval of chimney talus is thicker (up to ~3.9 m) in TKAD013, cemented by anhydrite-rich and clay-altered mud and finegrained sulfide minerals (App. Fig. A9). This interval was not recovered in TKAD025 (located ~3.6 m away) owing to poor drilling conditions (App. Fig. A10). In the center of the crater, intensely clay altered lapilli with minor disseminated sulfides occur to ~4 m b.s.f., below which is a chaotic mixture of matrix-supported volcanic and chimney fragments. The volcanic clasts are polymict (altered and unaltered), subangular to angular, and variably infilled by late-stage silica. Several short intervals within the lapilli units contain gypsum crystals up to 3 cm long. On the upper rim of the volcanic cone in the east, a thin (0.6 m) cover of unconsolidated volcaniclastic material is directly underlain by altered monomict lapillistone, with increasing sulfate cement and decreasing alteration intensity downhole to a depth of 7.4 m b.s.f. (see TKAD012, App. Fig. A11). At this site, there are fewer disseminated sulfides compared to the other drill holes and decreasing sphalerite abundance downhole.

Site C

Site C is a large volcanic cinder cone, rising to ~70 m above the seafloor and with a diameter of ~120 × 90 m (Figs. 2, 5C, D). Chimney complexes and mounds occur on the flank of this cone in the vicinity of a circular slump feature. TKAD032 intercepted unconsolidated volcaniclastic cover material to 1.8 m, followed by intensely altered volcaniclastic breccia to 6.9 m b.s.f. (App. Fig. A12). The altered sequence consists of monomict lapillistone with subrounded to subangular vesicular clasts but variable clay alteration. The uppermost lapillistone is matrix supported with variable sulfate cement and veins and disseminated pyrite, sphalerite, and chalcopyrite (decreasing downhole). Several horizons are clast supported (up to 60% clasts), with angular monomict clasts, and appear to be autobreccias. One gypsum-rich interval occurs at 6.4 to 6.6 m b.s.f. TKAD033 was drilled at the edge of a small crater or collapse feature near the summit of the cone. The recovered rock types are similar to those in TKAD032, with a transition between the volcaniclastic cover and the clay-altered footwall occurring at 1.1 m b.s.f. (App. Fig. A13). Variable alteration intensity gives this rock a polymict appearance.

Site D

Site D is an irregularly shaped topographic high, ~70 × 80 m in diameter, in the saddle between two cinder cones (Figs. 2, 5C, D). Chimneys occur on the north flank, where there is local venting of milky fluids from fractures near the base of the chimneys. On the northern side of the chimney cluster, TKAD017 intersected unconsolidated volcaniclastic cover material to ~1.4 m, separated from an interval of chimney talus by a layer of oxidized mud and chimney fragments (App. Fig. A14). The chimney talus extends to a depth of 5.2 m b.s.f. but is interrupted by a second layer of oxidized mud at depths of 2.6 to 3.3 m b.s.f. The altered footwall comprises subrounded, matrix-supported, monomict lapilli-sized clasts with minor disseminated pyrite and patchy anhydrite + gypsum cement that increases downhole to a depth of 7.4 m b.s.f. A few subrounded fragments of chimneys are also found in the footwall. The upper oxidized layer between the cover sequence and the mineralized interval can be traced ~20 m eastward to TKAD030 (App. Fig. A15). The lower oxidized layer was not recovered at TKAD030; instead, a buried sphalerite-rich chimney was intersected. The drilled sphalerite-rich chimney is underlain by chimney talus, which appears to be continuous with material in TKAD017. Altered lapillistone to 5.3 m b.s.f. consists of several beds that are distinguished by changes in grain size and proportion of clasts to matrix and by a higher proportion of sulfate cement. The two oxidized intervals encountered in TKAD017 were also traced 25 m to the south in TKAD031 (App. Fig. A16). The volcaniclastic cover sequence is absent in this drill hole, but a large chalcopyrite- and sphalerite-rich chimney block that likely came from a chimney complex located upslope to the south was drilled. The uppermost layer of chimney talus was similarly encountered here, extending to a similar depth (~5 m b.s.f.). The cover sequence is also absent in TKAD018 to the east. This hole started in oxidized chimney talus that is continuous with the uppermost oxidized interval in the other drill holes (App. Fig. A17). The chimney talus here is only ~0.6 m thick, grading to the clay-altered footwall at 1.2 m b.s.f. The clay-altered lapillistone at this location contains a higher proportion of fine-grained material (tuff). In the southeast, TKAD019 encountered only unconsolidated volcaniclastic cover material (to 6.1 m b.s.f.) that accumulated within a topographic low at the periphery of a chimney cluster (App. Fig. A18). Within the thick accumulation of cover material, there are several distinct intervals of intensely oxidized sediment, which may be continuations with the oxidized layers in the other holes.

Site E

Site E is on the southernmost cinder cone near several small craters (Figs. 2, 5E). TKAD005 is near an inactive chimney complex on the eastern margin of the northernmost crater; TKAD003 is on a small ridge south of TKAD005, near a scarp that may be a crater wall or slump and adjacent to a site of active venting. Both drill holes recovered moderately sorted, clast-supported, monomict lapillistone with weak grading and intense clay alteration. The matrix consists of anhydrite cement with minor disseminated sphalerite and pyrite that decrease in abundance downhole and are locally oxidized to hematite.

Site F

Site F is located west of the main NS-trending fissure, where chimneys and mounds are aligned in a circular pattern, but constructional volcanic features are absent (Figs. 2, 5F). TKAD007 was drilled next to a site of active venting in a sequence of clast-supported, clay-altered monomict lapillistone characterized by moderately sorted subangular clasts. The lapillistone is cemented by anhydrite, which increases in abundance downhole and is locally fracture controlled. The breccia also contains minor disseminated pyrite and fracture-controlled patchy iron oxide. A few centimeter-wide pyrite-chalcopyrite veinlets associated with the most intense clay alteration occur at the top of the recovered core. Inactive chimney complexes and mounds surround the area of active venting and appear to be aligned along structural corridors (Fig. 2). TKAD008 is located along one of these NW-SE–trending structures, 102 m south of TKAD007, on a low-relief ridge with numerous inactive chimney complexes and mounds. The upper part is matrix-supported lapillistone that grades downhole into a clast-supported lapillistone similar to that in TKAD008. Increasing clast size and angularity and decreasing alteration intensity downhole suggest that this drill hole may be at the periphery of the main hydrothermal upflow zone.

Seafloor and Subseafloor Mineralization

Isolated chimney complexes throughout the Tinakula area are most common on steep slopes that have not been completely buried by the volcaniclastic cover. Four main chimney types are recognized: (1) chalcopyrite rich, (2) zoned (with variable proportions of chalcopyrite and sphalerite), (3) sphalerite-rich chimneys, and (4) barite rich and/or silica altered (details in App. Table A1; App. Fig. A19). Unlike the subseafloor breccia facies, the seafloor chimneys contain only minor anhydrite and gypsum. Different chimney types commonly occur together; however, higher-temperature Cu-rich chimneys are more abundant near sites A and D, and lower-temperature barite-rich chimneys are common at site B (App. Fig. A20). Mineralized volcanic breccia was recovered from the base of several chimney complexes, where the chimneys have overgrown volcaniclastic material.

The Cu-rich chimneys comprise abundant, massive, coarsegrained chalcopyrite with minor pyrite (Fig. 7A-C). Locally, minor Fe-poor sphalerite and barite line the chimney conduits. In one sample, trace amounts of microscopic gold are present as inclusions in the chalcopyrite (sample GMA13-048; Fig. 7C). Distinctly zoned chimneys have layers that are chalcopyrite, pyrite (± marcasite), and sphalerite rich (Fig. 7D-F). Both Fe-rich and Fe-poor sphalerite occur in different layers. The low-Fe sphalerite commonly has inclusions of fine-grained Pb-As-Sb sulfosalts (jordanite) and minor galena. Late-stage, coarse tabular barite, amorphous silica, and iron oxyhydroxides occupy open spaces.

Fig. 7.

Features of chimney types at the Tinakula deposit with high porosities reflecting fluid channels. (A) Hand sample photo of chalcopyrite-rich chimney (sample GMA13-048). (B) Reflected-light in plane-polarized light (RL-PPL) photomicrograph of chalcopyrite-rich chimney with pyrite lining former vent conduit (sample GMA13-038). (C) RL-PPL photomicrograph of chalcopyrite-rich chimney with fine-grained gold inclusions within chalcopyrite (sample GMA13-048). (D) Hand sample photo of zoned chimney (sample GMA13-031). (E) RL-PPL photomicrograph of zoned chimney with chalcopyrite, marcasite (filling former vent conduit), and sphalerite (sample GMA13-050). (F) RL-PPL photomicrograph of zoned chimney with banded chalcopyrite and sphalerite (sample GMA13-031). (G) Hand sample photo of sphalerite-rich chimney (sample GMA13-012). (H) RL-PPL photomicrograph of dendritic to acicular barite, followed by sphalerite with minor late barite (sample GMA13-022). (I) RL-PPL photomicrograph of sphalerite-rich chimney, with early barite + spherulitic pyrite, followed by coarse-grained barite, and marcasite replaced by sphalerite (zone refining), with a late generation of globular sphalerite + galena (not shown) (sample GMA13-058). (J) Hand sample photo of barite-rich chimney (sample GMA13-045). (K) Transmitted-light in plane-polarized light (TL-PPL) photomicrograph of banded dendritic barite + pyrite, with minor sphalerite, and late tabular barite (sample GMA13-018B). (L) TL-PPL photomicrograph of dendritic to radial barite, tabular barite, and late amorphous silica with very fine grained inclusions of iron and manganese oxyhydroxides (sample GMA13-009B). Abbreviations: Am-Si = amorphous silica, Au = gold, Brt = barite, Ccp = chalcopyrite, Fe-MnOx = iron and manganese oxyhydroxides, Mrc = marcasite, Py = pyrite, Sp = sphalerite.

Fig. 7.

Features of chimney types at the Tinakula deposit with high porosities reflecting fluid channels. (A) Hand sample photo of chalcopyrite-rich chimney (sample GMA13-048). (B) Reflected-light in plane-polarized light (RL-PPL) photomicrograph of chalcopyrite-rich chimney with pyrite lining former vent conduit (sample GMA13-038). (C) RL-PPL photomicrograph of chalcopyrite-rich chimney with fine-grained gold inclusions within chalcopyrite (sample GMA13-048). (D) Hand sample photo of zoned chimney (sample GMA13-031). (E) RL-PPL photomicrograph of zoned chimney with chalcopyrite, marcasite (filling former vent conduit), and sphalerite (sample GMA13-050). (F) RL-PPL photomicrograph of zoned chimney with banded chalcopyrite and sphalerite (sample GMA13-031). (G) Hand sample photo of sphalerite-rich chimney (sample GMA13-012). (H) RL-PPL photomicrograph of dendritic to acicular barite, followed by sphalerite with minor late barite (sample GMA13-022). (I) RL-PPL photomicrograph of sphalerite-rich chimney, with early barite + spherulitic pyrite, followed by coarse-grained barite, and marcasite replaced by sphalerite (zone refining), with a late generation of globular sphalerite + galena (not shown) (sample GMA13-058). (J) Hand sample photo of barite-rich chimney (sample GMA13-045). (K) Transmitted-light in plane-polarized light (TL-PPL) photomicrograph of banded dendritic barite + pyrite, with minor sphalerite, and late tabular barite (sample GMA13-018B). (L) TL-PPL photomicrograph of dendritic to radial barite, tabular barite, and late amorphous silica with very fine grained inclusions of iron and manganese oxyhydroxides (sample GMA13-009B). Abbreviations: Am-Si = amorphous silica, Au = gold, Brt = barite, Ccp = chalcopyrite, Fe-MnOx = iron and manganese oxyhydroxides, Mrc = marcasite, Py = pyrite, Sp = sphalerite.

The Zn-rich chimneys are most abundant and consist of dendritic sphalerite and barite with pyrite and marcasite (Fig. 7G-I). These chimneys lack clearly defined vent conduits but have alternating layers of low- and high-Fe sphalerite, similar to beehive chimneys that form by diffuse fluid venting (e.g., Koski et al., 1994). The low-Fe sphalerite commonly contains inclusions of jordanite and galena; fine-grained inclusions of chalcopyrite are present at grain boundaries but are rare. Late-stage tabular barite is followed paragenetically by minor amounts of amorphous silica and trace iron and manganese oxyhydroxides. The Ba-rich chimneys also form bulbous and beehive-like structures with variable amounts of amorphous silica filling or lining open spaces (Fig. 7J-L).

Chalcopyrite and pyrite show little compositional variation, although As is relatively enriched in both pyrite (up to 1.42 wt %) and marcasite (up to 2.71 wt %) (App. Table A2). Mercury is also present above detection limits in pyrite, together with Pb and Sb. FeS contents of sphalerite range from 0.02 to 13.98 mol %, with more than 50% of the analyses reporting <2 mol % FeS (Fig. 8). The most Fe-rich sphalerite is in chalcopyrite-rich assemblages, and the widest range of FeS values occur in the dendritic barite-rich chimneys. Cadmium (up to 2.04 wt %), In (up to 1.31 wt %), Ga (up to 1.03 wt %), and Hg (up to 1.28 wt %) are likely present in solid solution in the sphalerite (e.g., Cook et al., 2009). Cadmium concentrations in sphalerite decrease from an average of 0.33 wt % in the chalcopyrite-rich chimneys (n = 16) to 0.03 wt % in the dendritic barite-rich chimneys (n = 79), suggesting a strong temperature dependence, which is also reflected in the Fe content. Galena contains up to 1.23 wt % Ag, mainly in the barite-rich chimneys. The Pb-As-Sb sulfosalts have up to 11.91 wt % As and 0.28 wt % Sb, which is close to those of the As end member of the jordanite-geocronite series (cf. Biagioni et al., 2016).

Fig. 8.

Frequency histograms of FeS (mol %) concentrations in sphalerite (uncorrected for chalcopyrite inclusions) from different chimney types at the Tinakula deposit (n = 598), calculated following Scott and Barnes (1971). Abbreviations as in Figure 7.

Fig. 8.

Frequency histograms of FeS (mol %) concentrations in sphalerite (uncorrected for chalcopyrite inclusions) from different chimney types at the Tinakula deposit (n = 598), calculated following Scott and Barnes (1971). Abbreviations as in Figure 7.

Bulk compositions of the seafloor chimneys and mineralized volcanic breccias are listed in Table 1. The seafloor chimneys contain abundant Zn up to 52.60 wt % (avg 19.20 wt %; n = 36) and Cu up to 30.00 wt % (avg 5.82 wt %; n = 36) but generally low Pb (<0.01–2.93 wt % Pb; avg 0.35 wt %; n = 36; Fig. 9) compared to other modern back-arc basin sulfide deposits (cf. Hannington et al., 2005; Monecke et al., 2016). Calcium concentrations are low in the chimneys (avg 0.06 wt %; n = 36), but Ba (avg 6.49 wt %; n = 36) and Sr (avg 2,157 ppm; n = 36; up to 0.85 wt % in barite-rich samples) are high. The SiO2 contents average 11.67 wt % (n = 36), but one chimney contains 75.09 wt % SiO2 together with 1.75 wt % Ba. The Cu-rich chimneys have high Mo (avg 343 ppm; n = 2) and In (avg 262 ppm; n = 2), with high Au (9,260 ppb) present in one sample. The Zn-rich samples contain very high Ga (up to 1,540 ppm), Ge (up to 99 ppm), Hg (up to 95 ppm), and Se (up to 20 ppm). Concentrations of Ag, Sb, As, and Au are also higher than in the Cu-rich chimneys, although high As is found in both Zn- and Cu-rich samples. Mercury concentrations are especially high in samples containing marcasite and late sphalerite + galena. The four sampled barite-rich chimneys also have high Pb (avg 7,073 ppm), As (avg 299 ppm), Sb (avg 82 ppm), Au (avg 2,038 ppb), and Ag (avg 121 ppm). Statistical correlations between the different trace elements clearly reflect the higher (>250°C) and lower temperatures (<250°C) of the Cu-rich and Zn-rich assemblages, respectively (App. Figs. A21, A22), as well as various lattice substitutions in the corresponding minerals (e.g., Cd in sphalerite, Sr in barite, As and Sb in Pb-rich sulfosalts, etc.; App. Table A2; cf. Hannington et al., 2005; Monecke et al., 2016). Gallium is incorporated mainly in galena in low-temperature assemblages, whereas in higher-temperature assemblages it mainly resides in sphalerite with moderate to high Fe contents. Arsenic is high in low-temperature assemblages, reflecting the presence of As in low-temperature minerals such as galena and sulfosalts. In a single sample containing secondary bornite (chalcopyrite-rich chimney GMA13-010), Ag concentrations are 1,253 ppm. Contents of Co, Ni, Sn, W, Bi, and Te are generally below detection limits in all samples, unlike other arc-related hydrothermal systems that commonly contain discrete Bi and Te minerals, e.g., Monecke et al. (2016).

Fig. 9.

Ternary diagram of Cu, Pb, and Zn concentrations in bulk whole-rock chimney samples from the Tinakula deposit, classified according to sample type (see text for classification). Abbreviations: Anh = anhydrite, Gp = gypsum, Volc. = volcanic; others as in Figure 7.

Fig. 9.

Ternary diagram of Cu, Pb, and Zn concentrations in bulk whole-rock chimney samples from the Tinakula deposit, classified according to sample type (see text for classification). Abbreviations: Anh = anhydrite, Gp = gypsum, Volc. = volcanic; others as in Figure 7.

Table 1.

Average Whole-Rock Geochemistry of Sulfide-Sulfate Samples by Chimney Type at Tinakula Deposit

  Chalcopyrite-richZoned chalcopyrite-richDendritic sphalerite-richBarite-rich
  Ccp n = 2Ccp + Py ± Mrc n = 3Ccp + Mrc + Sp n = 5Ccp + Sp n = 3Sp + Brt n = 13Sp + Brt + Py n = 1Sp + Brt + Mrc n = 2Sp + Mrc/Py n = 3Brt + Mrc ± Py n = 2
STotalwt %33.335.434.427.623.019.821.930.135.8
Auppb4,6921,0142,6703,5171,1334,2705,4553,5103,905
Fewt %25.6527.9721.7812.052.882.106.6014.8022.80
Cuwt %30.0017.0310.1410.520.680.810.381.520.51
Moppm34385536614110466711
Inppm2625957125142755<0.2
Nippm<10<10<10<10<101020<10<10
Bippm<2<2<2<2<2<2<2<2<2
Coppm2<0.2<0.2<0.2<0.2<0.2<0.221
Znwt %1.991.541.5421.5030.9227.1015.6125.023.39
Cdppm1446115.129631,0081,120601761155
Gappm2321682913165182613
Geppm92472746951305021
Seppm985411131211137
Tlppm1261232178384711
Pbppm444281,4011,9274,4415,3006,6002,44412,600
Hgppb53511,08030,12210,69719,1456,48019,05035,10027,300
Asppm545681,091579641425567625488
Sbppm9.355158193143179364299119
Agppm424266138139124272112210
Teppm<6<6<6<6<6<6<6<6<6
Snppm146.3195941<0.5<0.5
Wppm714231<0.7<0.7<0.7
Bawt %0.050.651.104.069.499.1513.492.338.58
Cawt %0.020.030.050.040.070.060.110.080.08
Srppm32094341,4452,7812,7306,2351,5394,135
Mgwt %0.030.02<0.01<0.01<0.01<0.01<0.01<0.01<0.01
Mnppm1020541189362273112,338578
Kwt %<0.1<0.1<0.1<0.1<0.1<0.1<0.1<0.1<0.1
Nawt %0.090.150.110.190.070.210.100.160.04
Siwt %0.63.33.25.14.39.15.06.72.3
Alwt %0.220.210.180.100.030.170.170.090.02
Tiwt %<0.01<0.01<0.01<0.01<0.01<0.01<0.01<0.01<0.01
Vppm374327351959203924
Pwt %0.0530.0330.0170.022<0.0050.012<0.0050.0530.007
  Chalcopyrite-richZoned chalcopyrite-richDendritic sphalerite-richBarite-rich
  Ccp n = 2Ccp + Py ± Mrc n = 3Ccp + Mrc + Sp n = 5Ccp + Sp n = 3Sp + Brt n = 13Sp + Brt + Py n = 1Sp + Brt + Mrc n = 2Sp + Mrc/Py n = 3Brt + Mrc ± Py n = 2
STotalwt %33.335.434.427.623.019.821.930.135.8
Auppb4,6921,0142,6703,5171,1334,2705,4553,5103,905
Fewt %25.6527.9721.7812.052.882.106.6014.8022.80
Cuwt %30.0017.0310.1410.520.680.810.381.520.51
Moppm34385536614110466711
Inppm2625957125142755<0.2
Nippm<10<10<10<10<101020<10<10
Bippm<2<2<2<2<2<2<2<2<2
Coppm2<0.2<0.2<0.2<0.2<0.2<0.221
Znwt %1.991.541.5421.5030.9227.1015.6125.023.39
Cdppm1446115.129631,0081,120601761155
Gappm2321682913165182613
Geppm92472746951305021
Seppm985411131211137
Tlppm1261232178384711
Pbppm444281,4011,9274,4415,3006,6002,44412,600
Hgppb53511,08030,12210,69719,1456,48019,05035,10027,300
Asppm545681,091579641425567625488
Sbppm9.355158193143179364299119
Agppm424266138139124272112210
Teppm<6<6<6<6<6<6<6<6<6
Snppm146.3195941<0.5<0.5
Wppm714231<0.7<0.7<0.7
Bawt %0.050.651.104.069.499.1513.492.338.58
Cawt %0.020.030.050.040.070.060.110.080.08
Srppm32094341,4452,7812,7306,2351,5394,135
Mgwt %0.030.02<0.01<0.01<0.01<0.01<0.01<0.01<0.01
Mnppm1020541189362273112,338578
Kwt %<0.1<0.1<0.1<0.1<0.1<0.1<0.1<0.1<0.1
Nawt %0.090.150.110.190.070.210.100.160.04
Siwt %0.63.33.25.14.39.15.06.72.3
Alwt %0.220.210.180.100.030.170.170.090.02
Tiwt %<0.01<0.01<0.01<0.01<0.01<0.01<0.01<0.01<0.01
Vppm374327351959203924
Pwt %0.0530.0330.0170.022<0.0050.012<0.0050.0530.007
  Silicified barite-richMineralized volcanic brecciaMineralized volcanic breccia-stringerMixed chimney talus
  Am-Si + Brt + Mrc/Py n = 2Sp + Brt + Volc. Clasts n = 1Sp + Brt + Volc. clasts n = 3Py + Brt + Volc. clasts n = 1Ccp + Brt + Volc. clasts n = 2Sp + Py + Brt + Gp n = 1Ccp + Sp + Py + Brt/Anh n = 1
STotalwt %5.110.217.413.724.922.231.7
Auppb1712212,3971,220255356552
Fewt %2.351.395.0911.7018.950.9122.20
Cuwt %0.020.050.775.4010.180.867.86
Moppm910291716322143
Inppm<0.2294844<0.230
Nippm<10<10<10<10<10<10<10
Bippm<2<2<2<2<2<2<2
Coppm262711<0.22
Znwt %0.233.3916.742.092.4535.605.16
Cdppm320562771091,610263
Gappm49203319324
Geppm2227911196632
Seppm15108101212
Tlppm12627315<0.149
Pbppm1,5453,8906,9189779539217
Hgppb35,88025,20023,1031,77090513,7001,070
Asppm112221830441255659306
Sbppm4647243191353014
Agppm3256110<5<5401bdl
Teppm<6<6<6<6<6<6<6
Snppm<0.5<0.5<0.511311
Wppm110112558
Bawt %15.1815.406.020.090.719.420.97
Cawt %0.110.640.130.090.090.040.02
Srppm3,1375,0602,03861641,740303
Mgwt %0.030.270.010.080.04<0.01<0.01
Mnppm6327258123611192
Kwt %<0.1<0.1<0.10.80.3<0.1<0.1
Nawt %0.230.430.190.260.140.100.04
Siwt %27.012.215.625.514.51.710.9
Alwt %0.141.790.552.951.300.090.13
Tiwt %<0.010.110.050.330.10<0.01<0.01
Vppm2710944160941866
Pwt %<0.0050.019<0.0050.0090.020<0.0050.012
  Silicified barite-richMineralized volcanic brecciaMineralized volcanic breccia-stringerMixed chimney talus
  Am-Si + Brt + Mrc/Py n = 2Sp + Brt + Volc. Clasts n = 1Sp + Brt + Volc. clasts n = 3Py + Brt + Volc. clasts n = 1Ccp + Brt + Volc. clasts n = 2Sp + Py + Brt + Gp n = 1Ccp + Sp + Py + Brt/Anh n = 1
STotalwt %5.110.217.413.724.922.231.7
Auppb1712212,3971,220255356552
Fewt %2.351.395.0911.7018.950.9122.20
Cuwt %0.020.050.775.4010.180.867.86
Moppm910291716322143
Inppm<0.2294844<0.230
Nippm<10<10<10<10<10<10<10
Bippm<2<2<2<2<2<2<2
Coppm262711<0.22
Znwt %0.233.3916.742.092.4535.605.16
Cdppm320562771091,610263
Gappm49203319324
Geppm2227911196632
Seppm15108101212
Tlppm12627315<0.149
Pbppm1,5453,8906,9189779539217
Hgppb35,88025,20023,1031,77090513,7001,070
Asppm112221830441255659306
Sbppm4647243191353014
Agppm3256110<5<5401bdl
Teppm<6<6<6<6<6<6<6
Snppm<0.5<0.5<0.511311
Wppm110112558
Bawt %15.1815.406.020.090.719.420.97
Cawt %0.110.640.130.090.090.040.02
Srppm3,1375,0602,03861641,740303
Mgwt %0.030.270.010.080.04<0.01<0.01
Mnppm6327258123611192
Kwt %<0.1<0.1<0.10.80.3<0.1<0.1
Nawt %0.230.430.190.260.140.100.04
Siwt %27.012.215.625.514.51.710.9
Alwt %0.141.790.552.951.300.090.13
Tiwt %<0.010.110.050.330.10<0.01<0.01
Vppm2710944160941866
Pwt %<0.0050.019<0.0050.0090.020<0.0050.012

Abbreviations: Am-Si = amorphous silica, Anh = anhydrite, bdl = below detection limit, Brt = barite, Ccp = chalcopyrite, Gp = gypsum, Mrc = marcasite, Py = pyrite, Sp = sphalerite, Volc. = volcanic

Subseafloor mineralization consists of (1) unconsolidated and locally hydrothermally cemented chimney talus (Fig. 10A, B), (2) chimney talus, mixed with altered volcanic clasts, that is generally hydrothermally cemented (Fig. 10C), (3) chimney talus associated with oxidized intervals (Fig. 10D), (4) disseminated to patchy sulfides (chalcopyrite, pyrite, sphalerite) within altered volcaniclastic rocks (Fig. 10E, F), and (5) sulfide minerals associated with late sulfate veins that extend to maximum depths of ~21 m b.s.l. (e.g., Figs. 3E, 4G, H; App. Fig. A2). Drilling shows that the majority of the mineralization at Tinakula comprises chimney talus that occurs directly beneath the unconsolidated volcaniclastic material (generally between 0.6 and 5.3 m b.s.f.). The oxidized intervals contain variable amounts of sulfide material (mostly chimney talus). Disseminated and vein-controlled sulfide mineralization (pyrite, chalcopyrite, and sphalerite) occur in the upper ~0.5 to 2.6 m of most drill holes (e.g., Fig. 4F, G). Seven samples of mineralized volcanic breccia down to 5.8 m b.s.f. contain mostly pyrite (9.7–30.1 wt % sulfide S) with minor anhydrite accounting for ~0.1 to 0.5 wt % sulfate S. Taking into account the dilution from wall rock, mineralized portions of the samples contain slightly higher Co and Se and lower Cd, Hg, Sb, and Ag than the seafloor chimneys, consistent with higher temperatures of mineralization (e.g., Monecke et al., 2016).

Fig. 10.

Subseafloor mineralization at the Tinakula deposit. (A) Fragment of barite- and sphalerite-rich chimney with a 3-cm-wide barite rim preserved in drill core (sample MA13-011). (B) Chimney talus with clasts of various sizes, shapes, and mineralogy, cemented by fine-grained sulfides and minor sulfates (sample MA13-030). (C) Clastic massive sulfides with sphalerite, chalcopyrite, and pyrite, and altered volcanic clasts, hydrothermally cemented by barite and fine-grained amorphous silica (sample MA13-032). (D) Oxidized unit characterized by abundant orange iron oxides and chimney fragments (bottom right: fragment of barite- and silica-rich chimney), with variable rounding from transport (sample MA13-022). (E) Barite-rich silicified volcaniclastic breccia with abundant chalcopyrite + pyrite (sample P011). (F) Lapilli-rich volcaniclastic tuff intensely altered to white clay, with disseminated pyrite, chalcopyrite, and minor sphalerite (sample MA13-013). Abbreviations: Si = amorphous silicon; others as in Figure 7.

Fig. 10.

Subseafloor mineralization at the Tinakula deposit. (A) Fragment of barite- and sphalerite-rich chimney with a 3-cm-wide barite rim preserved in drill core (sample MA13-011). (B) Chimney talus with clasts of various sizes, shapes, and mineralogy, cemented by fine-grained sulfides and minor sulfates (sample MA13-030). (C) Clastic massive sulfides with sphalerite, chalcopyrite, and pyrite, and altered volcanic clasts, hydrothermally cemented by barite and fine-grained amorphous silica (sample MA13-032). (D) Oxidized unit characterized by abundant orange iron oxides and chimney fragments (bottom right: fragment of barite- and silica-rich chimney), with variable rounding from transport (sample MA13-022). (E) Barite-rich silicified volcaniclastic breccia with abundant chalcopyrite + pyrite (sample P011). (F) Lapilli-rich volcaniclastic tuff intensely altered to white clay, with disseminated pyrite, chalcopyrite, and minor sphalerite (sample MA13-013). Abbreviations: Si = amorphous silicon; others as in Figure 7.

The chimney talus and oxidized intervals appear continuous between drill holes at sites A, B, and D, but mineralization was not encountered in the two drill holes at site C (Figs. 11, 12). Since the drill platform required relatively flat slopes, drilling at this site was limited to the areas beside the known chimney mounds. Therefore, the occurrence of chimney talus below the mounds was not confirmed at the location. The buried sulfide mineralization is estimated to have an area at least three times larger than that of the seafloor chimneys and mounds, although more drilling is required to establish the full extent and continuity.

Fig. 11.

Cross sections of the geology and alteration types interpreted from drill core (collars shown by circles), seafloor morphologies, and seafloor observations from (A) site A (looking north) and (B) site B (looking north-northwest). Cross section locations are shown in Figure 5A, B. Temperatures are calculated from δ18O data. Abbreviations: Cham = chamosite, Corr = corrensite, I/S = illite/smectite, Non = nontronite, Py = pyrite. Datum: WGS84 Zone 58S.

Fig. 11.

Cross sections of the geology and alteration types interpreted from drill core (collars shown by circles), seafloor morphologies, and seafloor observations from (A) site A (looking north) and (B) site B (looking north-northwest). Cross section locations are shown in Figure 5A, B. Temperatures are calculated from δ18O data. Abbreviations: Cham = chamosite, Corr = corrensite, I/S = illite/smectite, Non = nontronite, Py = pyrite. Datum: WGS84 Zone 58S.

Fig. 12.

Cross sections of the geology and alteration interpreted from drill core (collars shown by circles), seafloor morphologies, and seafloor observations from (A) site C (looking north) and (B) site B (looking north). Cross section locations are shown in Figure 5C, D. Temperatures are calculated from δ18O data. Abbreviations as in Figure 11. Datum: WGS84 Zone 58S.

Fig. 12.

Cross sections of the geology and alteration interpreted from drill core (collars shown by circles), seafloor morphologies, and seafloor observations from (A) site C (looking north) and (B) site B (looking north). Cross section locations are shown in Figure 5C, D. Temperatures are calculated from δ18O data. Abbreviations as in Figure 11. Datum: WGS84 Zone 58S.

Subseafloor Alteration

Below the mineralized interval, the volcaniclastic breccias have been intensely hydrothermally altered, with clay minerals occurring as replacements of volcanic material and infills of pore spaces (e.g., Fig. 4D-F). The predominant lithology of the footwall is monomict to polymict lapillistone that is crudely bedded with both matrix- and clast-supported horizons. Alteration and late infill by patchy sulfate cement obscures many of the primary volcanic textures, resulting in a high degree of uncertainty regarding facies interpretations. The altered volcaniclastic material extends to maximum drilled depths of 22.32 m b.s.f.; however, the full depth and areal extent of the alteration is unknown, with most of the drilling having been focused along the central row of cinder cones (sites A-E; Fig. 2). Alteration minerals in the footwall breccias identified by XRD include dioctahedral smectite (montmorillonite and nontronite), chlorite (mainly chamosite, with minor clinochlore), minor chloritoid, illite, and mixed-layer clays: illite/smectite, muscovite/illite (rectorite), and chlorite/smectite (corrensite) (Table 2; App. Fig. A23). These minerals constitute five alteration facies that are distinguished in drill core by different colors (Table 2): (1) montmorillonite/nontronite, with yellow-orange, yellow-gray, and red-brown hues; (2) nontronite + corrensite, with greenish-yellow to greenish-gray hues; (3) illite/smectite + pyrite, with yellowish-gray to olive-gray hues; (4) illite/smectite + chamosite, with bluish-white, olive, olive-gray, and greenish-yellow hues; and (5) chamosite + corrensite, with greenish-yellow to pale olive, with minor brown to reddish-brown hues.

Table 2.

Clay Minerals Identified by XRD from Samples of Altered Volcaniclastic Footwall from Tinakula Deposit, with Drill Holes Arranged According to Site from North to South (locations shown on Fig. 2)

Location Drill holeSampleDepth (m b.s.f.)Munsell clay colorSmectiteChloriteMixed-layer claysIlliteFeldsparNonclay mineralsAlteration assemblage
MntNonCldClcChamCorrAlie(?)RectI/SIllPgPlQzPyMrcHemAta(?)AnilBrtCel
Site A (North)
TKAD001MA13-0784.54–4.68Moderate yellowish brown+ +          + +        Mnt/Non + Pl
 MA13-0827.00–7.10Light yellowish                     
   orange+ +          + +        Mnt/Non + Pl
 P00210.39–10.50Dusky yellow + +     (?)   +        Mnt/Non + Pl
 MA13-08512.90–13.00Dark yellowish                     
   orange + +         + +        Mnt/Non + Pl
 MA13-08614.00–14.12Moderate reddish(?)+ +            Tr     Mnt/Non + Pl
   brown                     
 P00314.81–14.95Greenish gray + +Tr  +              Non + Corr
 P00415.65–15.80Light greenish gray  Tr  + + + + Tr          Non + Corr
 MA13-08716.05–16.22Pale greenish yellow(?)+ +   Tr              Non + Corr
 MA13-08817.40–17.50Light greenish gray    + ++ +  + +     Tr     I/S + Cham
 MA13-08917.90–18.00Light gray    + +  + +    Tr       I/S + Cham
 MA13-09018.40–18.57Light greenish gray    + +  + +TrTrTr         I/S + Cham
 P00519.11–19.30Ver)/light gray    +  + +TrTrTr  (?)      I/S + Cham
 MA13-09120.77–21.00Pale greenish yellowTr   +   + +    (?)      I/S + Cham
 MA13-09221.30–21.45Light greenish gray    + +Tr  Tr +         I/S + Cham
TKAD015MA13-0503.49–3.63Pale oliveTr   Tr   + + Tr   (?)     I/S + Cham
 MA13-0514.34–4.43Pale greenish yellow    + +   +           Cham + Corr
 MA13-0526.42–6.52Grayish yellow green    + +   Tr           Cham + Corr
 MA13-0537.23–7.34Grayish yellow green    + +   Tr           Cham + Corr
TKAD027MA13-0562.50–2.60Light olive    + +  + +      (?) (?)    (?)I/S + Cham
 MA13-0574.60–4.70Light greenish gray    + +   Tr           I/S + Cham
 MA13-0586.50–6.59Grayish yellow green    + +   +           I/S + Cham
TKAD028MA13-0994.50–4.60Pale greenish yellow    + +  (?)+    Tr      I/S + Cham
TKAD029MA13-0953.60–3.70Light yellowish                     
   orange + +         + +        Non + Pl
Site B
TKAD016MA13-0342.02–2.13Greenish gray      (?)+ +      Tr     I/S + Py
 MA13-0355.65–5.71Pale greenish yellow    +  + +(?)           I/S + Cham
 MA13-0367.00–7.05Yellowish gray    + +(?)+             Cham + Corr
TKAD025MA13-0395.50–5.60Pale olive    + +   Tr           I/S + Cham
 MA13-0406.91–7.00Pale olive(?)   + +Tr              Cham + Corr
TKAD026MA13-0672.90–3.00Light bluish gray    Tr  + +     Tr      I/S + Cham
 MA13-0683.60–3.68Light greenish grayTr   + +  (?)+ +           I/S + Cham
 MA13-0695.64–5.75Light greenish gray    + +  (?)+           I/S + Cham
TKAD012MA13-0622.81–2.90Pale olive    Tr+ +              Cham + Corr
 MA13-0634.25–4.36Light greenish yellow   Tr + +    Tr         Cham + Corr
 MA13-0646.49–6.60Greenish gray     + +  (?)     Tr     Cham + Corr
Site C
TKAD032MA13-0722.05–2.17Yellowish gray + +            Tr     Mnt/Non + Pl
 MA13-0733.90–4.00Greenish gray + +   Tr        Tr Tr     Non + Corr
 MA13-0745.31–5.40Grayish yellow + +   Tr        Tr     Non + Corr
 MA13-0756.42–6.54Greenish gray + +   Tr       Tr(?)     Non + Corr
TKAD033MA13-0253.00–3.10Dusky yellow + +         + + Tr   (?)  Mnt/Non + Pl
 MA13-0264.76–4.88Moderate brown + + (?)       + +   Tr    Mnt/Non + Pl
 MA13-0277.30–7.40Moderate brown   Tr + +     Tr   (?)    Cham + Corr
Site D
TKAD017MA13-0206.90–7.00Yellowish grayTr      Tr+ ++   +   +  I/S + Py
TKAD030MA13-0055.06–5.21Grayish yellow + +     + Tr   Tr   Tr  I/S + Py
 MA13-0066.45–6.60Yellowish gray       + ++Tr     + +  I/S + Py
TKAD018MA13-0143.00–3.10Light grayish                     
   yellow+       + +Tr   Tr   Tr  I/S + Py
 MA13-0155.00–5.10Yellowish grayTr      Tr+ +Tr     Tr Tr  I/S + Py
 MA13-0166.60–6.71Light grayish                     
   yellow+   Tr   + ++     Tr Tr  I/S + Cham
TKAD031MA13-0445.64–5.69Grayish YellowTr       + +Tr   (?)   Tr  I/S + Py
 MA13-0456.60–6.67Grayish YellowTr       + +Tr   (?)   Tr  I/S + Py
Site E (South)
TKAD005P02020.55–20.67Moderate greenish gray(?)   + +   + +(?)          I/S + Cham
TKAD003P01211.11–11.20Light greenish                     
   gray(?)   + +   +           I/S + Cham
 P01313.40–13.50Bluish white    + +   + +           I/S + Cham
 P01416.41–16.50Bluish white    + +   +           I/S + Cham
Site F (West)
TKAD007P0213.50–3.60Light olive gray     + +  (?)    Tr      I/S + Py
 P02211.23–11.34Bluish white    + +  (?)Tr           I/S + Cham
 P02322.35–22.45Pale greenish                     
   yellow    + +   Tr(?)          I/S + Cham
TKAD008P0257.50–7.56Moderate reddish brownTr   ++ +         (?)    Cham + Corr
 P02418.13–18.20Moderate reddish brown + +             (?)    Mnt/Non + Pl
Location Drill holeSampleDepth (m b.s.f.)Munsell clay colorSmectiteChloriteMixed-layer claysIlliteFeldsparNonclay mineralsAlteration assemblage
MntNonCldClcChamCorrAlie(?)RectI/SIllPgPlQzPyMrcHemAta(?)AnilBrtCel
Site A (North)
TKAD001MA13-0784.54–4.68Moderate yellowish brown+ +          + +        Mnt/Non + Pl
 MA13-0827.00–7.10Light yellowish                     
   orange+ +          + +        Mnt/Non + Pl
 P00210.39–10.50Dusky yellow + +     (?)   +        Mnt/Non + Pl
 MA13-08512.90–13.00Dark yellowish                     
   orange + +         + +        Mnt/Non + Pl
 MA13-08614.00–14.12Moderate reddish(?)+ +            Tr     Mnt/Non + Pl
   brown                     
 P00314.81–14.95Greenish gray + +Tr  +              Non + Corr
 P00415.65–15.80Light greenish gray  Tr  + + + + Tr          Non + Corr
 MA13-08716.05–16.22Pale greenish yellow(?)+ +   Tr              Non + Corr
 MA13-08817.40–17.50Light greenish gray    + ++ +  + +     Tr     I/S + Cham
 MA13-08917.90–18.00Light gray    + +  + +    Tr       I/S + Cham
 MA13-09018.40–18.57Light greenish gray    + +  + +TrTrTr         I/S + Cham
 P00519.11–19.30Ver)/light gray    +  + +TrTrTr  (?)      I/S + Cham
 MA13-09120.77–21.00Pale greenish yellowTr   +   + +    (?)      I/S + Cham
 MA13-09221.30–21.45Light greenish gray    + +Tr  Tr +         I/S + Cham
TKAD015MA13-0503.49–3.63Pale oliveTr   Tr   + + Tr   (?)     I/S + Cham
 MA13-0514.34–4.43Pale greenish yellow    + +   +           Cham + Corr
 MA13-0526.42–6.52Grayish yellow green    + +   Tr           Cham + Corr
 MA13-0537.23–7.34Grayish yellow green    + +   Tr           Cham + Corr
TKAD027MA13-0562.50–2.60Light olive    + +  + +      (?) (?)    (?)I/S + Cham
 MA13-0574.60–4.70Light greenish gray    + +   Tr           I/S + Cham
 MA13-0586.50–6.59Grayish yellow green    + +   +           I/S + Cham
TKAD028MA13-0994.50–4.60Pale greenish yellow    + +  (?)+    Tr      I/S + Cham
TKAD029MA13-0953.60–3.70Light yellowish                     
   orange + +         + +        Non + Pl
Site B
TKAD016MA13-0342.02–2.13Greenish gray      (?)+ +      Tr     I/S + Py
 MA13-0355.65–5.71Pale greenish yellow    +  + +(?)           I/S + Cham
 MA13-0367.00–7.05Yellowish gray    + +(?)+             Cham + Corr
TKAD025MA13-0395.50–5.60Pale olive    + +   Tr           I/S + Cham
 MA13-0406.91–7.00Pale olive(?)   + +Tr              Cham + Corr
TKAD026MA13-0672.90–3.00Light bluish gray    Tr  + +     Tr      I/S + Cham
 MA13-0683.60–3.68Light greenish grayTr   + +  (?)+ +           I/S + Cham
 MA13-0695.64–5.75Light greenish gray    + +  (?)+           I/S + Cham
TKAD012MA13-0622.81–2.90Pale olive    Tr+ +              Cham + Corr
 MA13-0634.25–4.36Light greenish yellow   Tr + +    Tr         Cham + Corr
 MA13-0646.49–6.60Greenish gray     + +  (?)     Tr     Cham + Corr
Site C
TKAD032MA13-0722.05–2.17Yellowish gray + +            Tr     Mnt/Non + Pl
 MA13-0733.90–4.00Greenish gray + +   Tr        Tr Tr     Non + Corr
 MA13-0745.31–5.40Grayish yellow + +   Tr        Tr     Non + Corr
 MA13-0756.42–6.54Greenish gray + +   Tr       Tr(?)     Non + Corr
TKAD033MA13-0253.00–3.10Dusky yellow + +         + + Tr   (?)  Mnt/Non + Pl
 MA13-0264.76–4.88Moderate brown + + (?)       + +   Tr    Mnt/Non + Pl
 MA13-0277.30–7.40Moderate brown   Tr + +     Tr   (?)    Cham + Corr
Site D
TKAD017MA13-0206.90–7.00Yellowish grayTr      Tr+ ++   +   +  I/S + Py
TKAD030MA13-0055.06–5.21Grayish yellow + +     + Tr   Tr   Tr  I/S + Py
 MA13-0066.45–6.60Yellowish gray       + ++Tr     + +  I/S + Py
TKAD018MA13-0143.00–3.10Light grayish                     
   yellow+       + +Tr   Tr   Tr  I/S + Py
 MA13-0155.00–5.10Yellowish grayTr      Tr+ +Tr     Tr Tr  I/S + Py
 MA13-0166.60–6.71Light grayish                     
   yellow+   Tr   + ++     Tr Tr  I/S + Cham
TKAD031MA13-0445.64–5.69Grayish YellowTr       + +Tr   (?)   Tr  I/S + Py
 MA13-0456.60–6.67Grayish YellowTr       + +Tr   (?)   Tr  I/S + Py
Site E (South)
TKAD005P02020.55–20.67Moderate greenish gray(?)   + +   + +(?)          I/S + Cham
TKAD003P01211.11–11.20Light greenish                     
   gray(?)   + +   +           I/S + Cham
 P01313.40–13.50Bluish white    + +   + +           I/S + Cham
 P01416.41–16.50Bluish white    + +   +           I/S + Cham
Site F (West)
TKAD007P0213.50–3.60Light olive gray     + +  (?)    Tr      I/S + Py
 P02211.23–11.34Bluish white    + +  (?)Tr           I/S + Cham
 P02322.35–22.45Pale greenish                     
   yellow    + +   Tr(?)          I/S + Cham
TKAD008P0257.50–7.56Moderate reddish brownTr   ++ +         (?)    Cham + Corr
 P02418.13–18.20Moderate reddish brown + +             (?)    Mnt/Non + Pl

Amounts indicated are major (++), minor (+), and trace (Tr); uncertain identifications are indicated by (?)

Mineral abbreviations: Alie = aliettite, Anh = anhydrite, Ata = paratacamite(?), Brt = barite, Cel = celestine, Cham = chamosite, Clc = clinochlore, Cld = chloritoid, Corr = corrensite, Hem = hematite, Ill = illite, I/S = illite/smectite, Mnt = montmorillonite, Mrc = marcasite, Non = nontronite, Pg = paragonite, Ph = phillipsite, Pl = plagioclase, Py = pyrite, Qz = quartz, Rect = rectorite

The alteration assemblages vary with location and depth, generally progressing downward from the seafloor in the order montmorillonite/nontronite to nontronite + corrensite, to illite/smectite + pyrite, to illite/smectite + chamosite, to chamosite + corrensite. Footwall breccias from site A are dominated by chamosite + corrensite alteration, with localized illite/smectite + chamosite alteration that varies in abundance over short distances and is interpreted to reflect proximity to inferred synvolcanic faults (Fig. 11A). Alteration at site B is semiconformable, varying downhole from illite/smectite + pyrite, to illite/smectite + chamosite, to chamosite + corrensite (Fig. 11B). Site C is characterized by nontronite near the surface, followed by nontronite + corrensite and chamosite + corrensite at depth (Fig. 12A). Alteration at site D varies from illite/smectite + pyrite near the surface, to illite/smectite + chamosite at depth (Fig. 12B). End-member illite and rectorite is present in all samples from site D, reflecting high temperatures (>220°C) and a progression from argillic to phyllic alteration at depth (e.g., Inoue et al., 1990). The variations in alteration assemblages identified at different sites along the NS-trending fissure appear to reflect variable thermal gradients and permeability structure.

Smectite is ubiquitous throughout the footwall in every alteration assemblage as discrete montmorillonite/nontronite replacing volcanic glass and/or mixed-layer clays. Mixed-layer illite/smectite clays range in composition from 50% illite to ~91% illite, with an average of 74% illite (n = 30) (App. Table A3). The proportion of illite increases progressively downhole, indicating a process of continuous transformation of smectite to illite (e.g., Inoue et al., 1990). Chlorite may form directly from smectite or through an intermediate corrensite step (e.g., Shau et al., 1990). At Tinakula, chamosite appears before corrensite in the downhole progression of alteration assemblages, suggesting that corrensite is not an intermediate phase. Corrensite appears after illite/smectite, following a transition from illite/smectite + chamosite to chamosite + corrensite. This suggests that corrensite forms as an alteration product of the mixed-layer illite/smectite clays.

Alteration minerals were also characterized by SWIR spectroscopy (Table 3; App. Fig. A24). The white micas are identified by the distinct Al-OH absorption feature between 2,180 and 2,228 nm, with wavelengths of 2,189 to 2,195 nm for paragonite and 2,200 to 2,213 nm for muscovite and illite (this study; Thompson et al., 1999). In addition, small absorption features are present at 2,350 and 2,450 nm for the white micas. Montmorillonite is distinguished by a narrow and symmetrical 2,200-nm feature with a small shoulder on the H2O + OH absorption feature at 1,450 to 1,460 nm. Chlorite is identified by the strong absorption feature at 2,340 nm, the wavelength of which increases from 2,324 nm for Mg-rich compositions to >2,350 nm for Fe-rich compositions. Gypsum is identified by the OH + H2O feature at 1,400 nm, together with the OH feature at 1,800 nm. The alteration mineral assemblages observed in the SWIR analysis progress downhole from smectite (montmorillonite/nontronite), to white mica (paragonite → paragonitic illite → muscovitic illite → illite), to chlorite (e.g., App. Fig. A25), consistent with the XRD results. However, SWIR analyses with both montmorillonite and white mica reflect the presence of mixed-layer illite/smectite clays rather than discrete end members, with progressive increases downhole in aluminum content that shift the Al-OH absorption feature to lower wavelengths (e.g., Velde, 1980).

Table 3.

Clay Minerals Identified by SWIR Spectroscopy from Samples of Altered Volcaniclastic Footwall from Tinakula Deposit, with Drill Holes Arranged from North to South (see Fig. 2)

Location Drill holeSampleDepth (m b.s.f.)SmectiteChloriteWhite micaSulfateAlO(OH)
Mnt (%)Non (%)Chl (%)Pg (%)Pgl (%)Msl (%)Ms (%)Gp (%)Dsp (%)
Site A (North)
TKAD001MA13-086-C14.00–14.12+++      + 
 MA13-087-C16.05–16.22+++   ++    
 MA13-088-C17.40–17.50    (?) +++  
 P005-C19.11–19.30(?)   +++    
 MA13-091-CM20.77–21.00  ++   ++  
TKAD015MA13-052-C6.42–6.52  (?)+++   + 
 MA13-053-C7.23–7.34  (?)+++   ++ 
TKAD028MA13-099-C4.50–4.60       +++ + + *
Site B
TKAD016MA13-034-C2.02–2.13       +++ + + *
 MA13-036-C7.00–7.05(?) +++++     
TKAD025MA13-039-C5.50–5.60(?) (?)    +++ 
 MA13-040-C6.91–7.00  +++    + 
TKAD026MA13-067-C2.90–3.00++  +++     
 MA13-068-C3.60–3.68(?) (?)+++     
 MA13-069-C5.64–5.75   (?)*   (?)* 
TKAD012MA13-062-C2.81–2.90++ (?) ++    
Site C
TKAD032MA13-072-C2.05–2.17+++        
 MA13-073-C3.90–4.00+++        
 MA13-074-C5.31–5.40+++ (?) ++    
 MA13-075-C6.42–6.54+++ (?) ++    
TKAD033MA13-026-C4.76–4.88+++        
 MA13-026-M4.76–4.88+++        
 MA13-027-C7.30–7.40+++        
 MA13-027-M7.30–7.40        + + + *
Site D
TKAD030MA13-005-C5.06–5.21(?)  +++     
 MA13-006-C6.45–6.60(?)  +++     
TKAD018MA13-015-C5.00–5.10   ++   +++ 
TKAD031MA13-044-C5.64–5.69 ++  +++    
Site E (South)
TKAD003P013-C13.40–13.50++ (?) +++    
 P014-C16.41–16.50(?) ++ ++    
Site F (West)
TKAD007P022-C11.23–11.34(?) ++ +++    
Location Drill holeSampleDepth (m b.s.f.)SmectiteChloriteWhite micaSulfateAlO(OH)
Mnt (%)Non (%)Chl (%)Pg (%)Pgl (%)Msl (%)Ms (%)Gp (%)Dsp (%)
Site A (North)
TKAD001MA13-086-C14.00–14.12+++      + 
 MA13-087-C16.05–16.22+++   ++    
 MA13-088-C17.40–17.50    (?) +++  
 P005-C19.11–19.30(?)   +++    
 MA13-091-CM20.77–21.00  ++   ++  
TKAD015MA13-052-C6.42–6.52  (?)+++   + 
 MA13-053-C7.23–7.34  (?)+++   ++ 
TKAD028MA13-099-C4.50–4.60       +++ + + *
Site B
TKAD016MA13-034-C2.02–2.13       +++ + + *
 MA13-036-C7.00–7.05(?) +++++     
TKAD025MA13-039-C5.50–5.60(?) (?)    +++ 
 MA13-040-C6.91–7.00  +++    + 
TKAD026MA13-067-C2.90–3.00++  +++     
 MA13-068-C3.60–3.68(?) (?)+++     
 MA13-069-C5.64–5.75   (?)*   (?)* 
TKAD012MA13-062-C2.81–2.90++ (?) ++    
Site C
TKAD032MA13-072-C2.05–2.17+++        
 MA13-073-C3.90–4.00+++        
 MA13-074-C5.31–5.40+++ (?) ++    
 MA13-075-C6.42–6.54+++ (?) ++    
TKAD033MA13-026-C4.76–4.88+++        
 MA13-026-M4.76–4.88+++        
 MA13-027-C7.30–7.40+++        
 MA13-027-M7.30–7.40        + + + *
Site D
TKAD030MA13-005-C5.06–5.21(?)  +++     
 MA13-006-C6.45–6.60(?)  +++     
TKAD018MA13-015-C5.00–5.10   ++   +++ 
TKAD031MA13-044-C5.64–5.69 ++  +++    
Site E (South)
TKAD003P013-C13.40–13.50++ (?) +++    
 P014-C16.41–16.50(?) ++ ++    
Site F (West)
TKAD007P022-C11.23–11.34(?) ++ +++    

Amounts indicated are >50% (+++), 25–50% (++), <25% (+); * = unreliable results

(?) indicates manual search of the spectral library; all others by Spectral Assist automated identification using TSG 7 software

Abbreviations: C = clast separate, Chl = chlorite, CM = clast + matrix, Dsp = diaspore, Gp = gypsum, M = matrix separate, Mnt = montmorillonite, Ms = muscovite, Msl = muscovitic illite, Pg = paragonite, Pgl = paragonitic illite

Alteration Geochemistry

Least altered volcanic rocks from the Tinakula deposit range from basaltic andesite and andesite to dacite, with basaltic andesite most abundant within and near the volcanic cones. Representative whole-rock geochemical analyses of altered volcaniclastic rocks from drill core are given in Table 4, with the associated alteration assemblages indicated. The least altered host rocks, described by Anderson (2018), were identified petrographically by a lack of alteration minerals, loss on ignition (LOI) <1.2 wt %, low base metal values (<183 ppm Cu + Zn + Pb), and no obvious Na2O depletion (>4.0 wt %). Data for the three least altered basaltic andesite samples lie on linear trends, with those for the altered samples in plots of TiO2 versus Zr and TiO2 versus Al2O3 (Fig. 13). Excursions from this linear trend (e.g., Zr >160 ppm at ~1.5 wt % TiO2; Fig. 13A) may reflect variability in the parental rock compositions, as volcanic rocks in this area show progressive fractionation from basaltic andesite to andesite (e.g., Anderson, 2018). The concentration ranges of immobile and mobile elements for each alteration assemblage are shown in Figure 14. The narrow range of concentrations of Zr, Y, Sc, Nb, Ti, and Al reflect the fact that these elements are least mobile. Concentrations of Ba, Si, Al, Fe, S, Mg, Ca, Na, and K show significant variation, consistent with their inferred mobility during hydrothermal alteration. High Ba, Ca, and S(T) reflect the abundance of barite and anhydrite/gypsum in the samples. High Fe2O3(T) and S(T) are related to pyrite occurrence. Illite-smectite alteration is characterized by relative enrichments in K2O compared to the least altered basaltic andesite, whereas chamosite + corrensite alteration is dominated by MgO addition. The bulk compositions are transitional between illite/smectite + pyrite to illite/smectite + chamosite to chamosite + corrensite end members (Fig. 15).

Fig. 13.

Immobile element ratios of altered volcanic rocks and least altered precursor (average of three samples, from Anderson, 2018) showing linear trends of relative mass gains and losses, with increasing alteration intensity reflected by increasing distance from average least altered precursor composition: (A) Zr versus TiO2 and (B) Al2O3 versus TiO2.

Fig. 13.

Immobile element ratios of altered volcanic rocks and least altered precursor (average of three samples, from Anderson, 2018) showing linear trends of relative mass gains and losses, with increasing alteration intensity reflected by increasing distance from average least altered precursor composition: (A) Zr versus TiO2 and (B) Al2O3 versus TiO2.

Fig. 14.

Box plots of selected whole-rock immobile and mobile elements associated with each alteration assemblage (full data in Anderson, 2018). Top of each box represents the 75th percentile, bottom is the 25th percentile, central line represents median, bars show range of data, and outliers are indicated by stars (*). Abbreviation: Mnt = montmorillonite; others as in Figure 11.

Fig. 14.

Box plots of selected whole-rock immobile and mobile elements associated with each alteration assemblage (full data in Anderson, 2018). Top of each box represents the 75th percentile, bottom is the 25th percentile, central line represents median, bars show range of data, and outliers are indicated by stars (*). Abbreviation: Mnt = montmorillonite; others as in Figure 11.

Fig. 15.

Ternary plots of whole-rock geochemistry distinguishing alteration types. (A) MgO-Na2O-K2O distinguishing chlorite + smectite from illite alteration. (B) Al2O3-MgO-Na2O+K2O distinguishing chlorite from illite alteration (modified from MacDonald et al., 1996). General trends associated with increasing depth are shown. Abbreviations: B-A = basaltic andesite, Mnt = montmorillonite; others as in Figure 11.

Fig. 15.

Ternary plots of whole-rock geochemistry distinguishing alteration types. (A) MgO-Na2O-K2O distinguishing chlorite + smectite from illite alteration. (B) Al2O3-MgO-Na2O+K2O distinguishing chlorite from illite alteration (modified from MacDonald et al., 1996). General trends associated with increasing depth are shown. Abbreviations: B-A = basaltic andesite, Mnt = montmorillonite; others as in Figure 11.

Table 4.

Whole-Rock Geochemical Analyses of Select Altered Volcaniclastic Rocks from Tinakula Deposit with Associated Alteration Assemblages, Normalized to 100% LOI-Free, and of Average Least Altered Basaltic-Andesite Composition (average of three samples: GMA13-057, P0016, and P0040; Anderson, 2018)

ElementDetection limitMethodMA13-086-C TKAD001 Mnt/NonMA13-073-C TKAD032 Non + CorrMA13-006-C TKAD030 I/S + PyMA13-088-C TKAD001 I/S + ChamMA13-036-C TKAD016 Cham + CorrAverage least altered n = 3 B-A
SiO2 (%)0.01FUS-ICP47.0351.7241.6552.3435.2754.39
TiO20.001FUS-ICP1.5861.7182.4671.9831.8041.30
Al2O30.01FUS-ICP19.2821.9327.0023.8022.8517.27
Fe2O3(T)0.01FUS-ICP11.8913.6920.248.9923.189.98
MnO0.001FUS-ICP0.0780.0450.0100.2530.0790.189
MgO0.01FUS-ICP13.387.650.636.429.113.57
CaO0.01FUS-ICP4.200.650.691.163.438.50
Na2O0.01FUS-ICP1.582.341.501.051.334.09
K2O0.01FUS-ICP0.850.245.723.942.850.55
P2O50.01FUS-ICP0.140.040.070.060.110.14
S(T)0.01CS1.347.0013.302.7312.700.01
Cu (ppm)1TD-ICP1112371281412,25073
Zn1TD-ICP16251,08013592812098
Ag0.3TD-ICP<0.33.11.30.60.6<0.3
As0.5INAA8.856.3226.024.528.41.75
Bi0.1FUS-MS<0.1<0.1<0.1<0.1<0.1<0.1
Cd0.5TD-ICP<0.51.4<0.52.1<0.5<0.5
Ga1FUS-MS171728232217
Ge0.5FUS-MS1.70.9<0.51.0<0.51.3
Au (ppb)2INAA<242948251<2
Hg51G5902142980.69
In (ppm)0.1FUS-MS<0.10.30.30.3<0.20.1
Sb0.2INAA11.310.43.41.3<30.5
Se3INAA<3<3<3<3<3<3
Sn1FUS-MS<1<1<1117
Tl0.05FUS-MS<0.051.140.400.760.340.07
Ir5INAA<5<5<5<5<5<5
Mo2FUS-MS<2128<27<2
Co1FUS-MS92711171723
Cr5INAA<526<513<57
Ni1TD-ICP5245768
Sc0.1INAA21.624.832.628.821.224.8
V5FUS-ICP90316350360283351
Ba2FUS-ICP1307,1942,2281,559395108
Pb3TD-ICP675224030<3<3
Rb1FUS-MS1446456286
Cs0.1FUS-MS0.90.20.30.30.20.2
Sr2FUS-ICP3542022348858278
Zr1FUS-ICP1019116315613989
Hf0.1FUS-MS2.42.13.63.93.32.0
Nb0.2FUS-MS4.63.86.66.75.82.9
Ta0.01FUS-MS0.220.180.430.360.220.18
Th0.05FUS-MS0.840.671.251.351.030.68
U0.01FUS-MS6.2913.2016.612.908.150.38
Be1FUS-ICP22<11<1<1
W0.5FUS-MS1.00.6178.02.8124.00.6
Y0.5FUS-MS27.820.838.433.337.129.2
La0.05FUS-MS6.164.9875.5014.9015.206.89
Ce0.05FUS-MS14.9012.10157.0031.9033.9016.10
Pr0.01FUS-MS2.191.7618.804.374.452.42
Nd0.05FUS-MS11.008.8674.0019.2020.1011.0
Sm0.01FUS-MS3.042.8913.604.895.233.32
Eu0.005FUS-MS1.321.417.582.532.671.27
Gd0.01FUS-MS3.903.258.245.585.904.15
Tb0.01FUS-MS0.690.61.130.951.010.73
Dy0.01FUS-MS4.713.917.486.126.614.80
Ho0.01FUS-MS1.040.771.451.311.371.00
Er0.01FUS-MS2.982.254.503.774.112.85
Tm0.005FUS-MS0.4680.3330.6560.5450.6080.427
Yb0.01FUS-MS3.102.274.393.394.123.72
Lu0.002FUS-MS0.4850.3650.7260.5210.6550.596
ElementDetection limitMethodMA13-086-C TKAD001 Mnt/NonMA13-073-C TKAD032 Non + CorrMA13-006-C TKAD030 I/S + PyMA13-088-C TKAD001 I/S + ChamMA13-036-C TKAD016 Cham + CorrAverage least altered n = 3 B-A
SiO2 (%)0.01FUS-ICP47.0351.7241.6552.3435.2754.39
TiO20.001FUS-ICP1.5861.7182.4671.9831.8041.30
Al2O30.01FUS-ICP19.2821.9327.0023.8022.8517.27
Fe2O3(T)0.01FUS-ICP11.8913.6920.248.9923.189.98
MnO0.001FUS-ICP0.0780.0450.0100.2530.0790.189
MgO0.01FUS-ICP13.387.650.636.429.113.57
CaO0.01FUS-ICP4.200.650.691.163.438.50
Na2O0.01FUS-ICP1.582.341.501.051.334.09
K2O0.01FUS-ICP0.850.245.723.942.850.55
P2O50.01FUS-ICP0.140.040.070.060.110.14
S(T)0.01CS1.347.0013.302.7312.700.01
Cu (ppm)1TD-ICP1112371281412,25073
Zn1TD-ICP16251,08013592812098
Ag0.3TD-ICP<0.33.11.30.60.6<0.3
As0.5INAA8.856.3226.024.528.41.75
Bi0.1FUS-MS<0.1<0.1<0.1<0.1<0.1<0.1
Cd0.5TD-ICP<0.51.4<0.52.1<0.5<0.5
Ga1FUS-MS171728232217
Ge0.5FUS-MS1.70.9<0.51.0<0.51.3
Au (ppb)2INAA<242948251<2
Hg51G5902142980.69
In (ppm)0.1FUS-MS<0.10.30.30.3<0.20.1
Sb0.2INAA11.310.43.41.3<30.5
Se3INAA<3<3<3<3<3<3
Sn1FUS-MS<1<1<1117
Tl0.05FUS-MS<0.051.140.400.760.340.07
Ir5INAA<5<5<5<5<5<5
Mo2FUS-MS<2128<27<2
Co1FUS-MS92711171723
Cr5INAA<526<513<57
Ni1TD-ICP5245768
Sc0.1INAA21.624.832.628.821.224.8
V5FUS-ICP90316350360283351
Ba2FUS-ICP1307,1942,2281,559395108
Pb3TD-ICP675224030<3<3
Rb1FUS-MS1446456286
Cs0.1FUS-MS0.90.20.30.30.20.2
Sr2FUS-ICP3542022348858278
Zr1FUS-ICP1019116315613989
Hf0.1FUS-MS2.42.13.63.93.32.0
Nb0.2FUS-MS4.63.86.66.75.82.9
Ta0.01FUS-MS0.220.180.430.360.220.18
Th0.05FUS-MS0.840.671.251.351.030.68
U0.01FUS-MS6.2913.2016.612.908.150.38
Be1FUS-ICP22<11<1<1
W0.5FUS-MS1.00.6178.02.8124.00.6
Y0.5FUS-MS27.820.838.433.337.129.2
La0.05FUS-MS6.164.9875.5014.9015.206.89
Ce0.05FUS-MS14.9012.10157.0031.9033.9016.10
Pr0.01FUS-MS2.191.7618.804.374.452.42
Nd0.05FUS-MS11.008.8674.0019.2020.1011.0
Sm0.01FUS-MS3.042.8913.604.895.233.32
Eu0.005FUS-MS1.321.417.582.532.671.27
Gd0.01FUS-MS3.903.258.245.585.904.15
Tb0.01FUS-MS0.690.61.130.951.010.73
Dy0.01FUS-MS4.713.917.486.126.614.80
Ho0.01FUS-MS1.040.771.451.311.371.00
Er0.01FUS-MS2.982.254.503.774.112.85
Tm0.005FUS-MS0.4680.3330.6560.5450.6080.427
Yb0.01FUS-MS3.102.274.393.394.123.72
Lu0.002FUS-MS0.4850.3650.7260.5210.6550.596

Method abbreviations refer to laboratory methods used by ActLabs: 1G = aqua regia Hg cold vapor flow injection technique, CS = total sulfur by IR, FUS = fusion, TD = total digestion; ICP = inductively coupled plasma, INAA = instrumental neutron activation analysis, MS = mass spectrometry Mineral abbreviations: B-A = basaltic-andesite, Cham = chamosite, Corr = corrensite, I/S = illite/smectite, Mnt = montmorillonite, Non = nontronite, Pl = plagiolcase, Py = pyrite

1

In the original data, several analyses of Cu and Zn that were over limits were reanalyzed with a different method; averages shown here include both the TD-ICP and 4 acid ICP-OES data

The chondrite-normalized REE profiles for the altered volcaniclastic rocks and the least altered basaltic andesite samples (Fig. 16) show pronounced positive Eu anomalies, as Eu becomes concentrated in hydrothermal fluids after it is released from plagioclase during water-rock interactions. Smectite-rich assemblages (nontronite/montmorillonite and nontronite + corrensite) show slight REE depletion and have flat profiles. The illite/smectite + pyrite assemblage displays significant light REE (LREE) enrichment and steeper profiles, whereas the chamosite-rich assemblages (illite/smectite + chamosite and chamosite + corrensite) show moderate LREE enrichment and flat profiles. The change in slope of the LREEs during the main stage of alteration (illite/smectite + pyrite) relative to the least altered rock indicates strong mobility. In contrast, the heavy REEs (HREEs) have slopes similar to those of the least altered rocks, which indicates that they are immobile during hydrothermal alteration. The HREEs are simply moved up or down relative to the least altered values due to the total mass gain or loss of the rock during alteration.

Fig. 16.

Chondrite-normalized REE plots of different alteration assemblages compared to least altered basalts/andesites (Anderson, 2018). (A) Least altered basalt/andesite samples. (B) Nontronite/montmorillonite-altered samples with REE depletion and flat profiles. (C) Nontronite + corrensite-altered samples with REE depletion and flat profiles. (D) Illite/smectite + pyrite-altered samples with LREE enrichment in some samples and downward-sloping REE profiles. (E) Illite/smectite + chamosite-altered samples with moderate LREE enrichment in some samples and flat to downward-sloping profiles. (F) Chamosite + corrensite-altered samples with slight LREE enrichment in most samples and slightly downward sloping profiles. All alteration types are associated with pronounced positive Eu anomalies. Chondrite data from Boynton (1984).

Fig. 16.

Chondrite-normalized REE plots of different alteration assemblages compared to least altered basalts/andesites (Anderson, 2018). (A) Least altered basalt/andesite samples. (B) Nontronite/montmorillonite-altered samples with REE depletion and flat profiles. (C) Nontronite + corrensite-altered samples with REE depletion and flat profiles. (D) Illite/smectite + pyrite-altered samples with LREE enrichment in some samples and downward-sloping REE profiles. (E) Illite/smectite + chamosite-altered samples with moderate LREE enrichment in some samples and flat to downward-sloping profiles. (F) Chamosite + corrensite-altered samples with slight LREE enrichment in most samples and slightly downward sloping profiles. All alteration types are associated with pronounced positive Eu anomalies. Chondrite data from Boynton (1984).

Major element mass changes were calculated for all alteration assemblages following the single precursor method of MacLean and Barrett (Maclean, 1990; Maclean and Barrett, 1993; Barrett and MacLean, 1999) using the average concentration of the three least altered samples as the precursor. In this method, the ratio of the concentration of an immobile element in an altered rock to the concentration of that element in the least altered precursor is proportional to the total mass change in the rock due to hydrothermal alteration. We used the proportion of the least mobile element, Zr, in the precursor sample relative to the altered sample to calculate the enrichment factor (EF). The reconstructed composition (RC) of the altered sample was calculated for each element by multiplying its concentration in the altered sample (normalized to 100% volatile free) by the EF. Finally, the mass change for each element in the sample was calculated by subtracting the precursor concentration of that element in the least altered sample from the RC. The largest mass changes are associated with moderate to intense alteration, generally increasing with depth (App. Figs. A26, A27). Mass changes for each assemblage are outlined below (Fig. 17; Table 5):

  1. Montmorillonite/nontronite alteration is associated with Si and Na losses in smectite-rich samples, Ca loss in all samples, no gains in Al or Fe (either conserved or lost), and no change in K.

  2. Nontronite + corrensite alteration is associated with Si, Ca, and Na losses, Mg gains, no losses of Al or Fe (either conserved or gained), and no change in K.

  3. Illite/smectite + pyrite alteration is associated with losses in Si, Mg, Na, and Ca and gains in both Fe (due to pyrite) and K.

  4. Illite/smectite + chamosite alteration is associated with Si and Na losses, K gains, and variable gains and losses in Fe and Mg (Mg losses where illite/smectite is more abundant than chamosite; Fe losses where illite/smectite and chamosite are equally abundant).

  5. Chamosite + corrensite alteration is associated with Si and Na losses, Mg gains, Fe gains when chamosite is more abundant than corrensite and Fe losses when corrensite is more abundant, K gains when minor illite/smectite is present, and no changes in K where no illite/smectite is present.

Fig. 17.

Average mass changes for each alteration mineral assemblage. Abbreviation: Mnt = montmorillonite; others as in Figure 11.

Fig. 17.

Average mass changes for each alteration mineral assemblage. Abbreviation: Mnt = montmorillonite; others as in Figure 11.

Table 5.

Mass Changes Associated with Alteration Mineral Assemblages in Altered Volcaniclastic Rocks Relative to at Least Altered Basaltic Andesite Precursor

Alteration SampleSiO2TiO2Al2O3Fe2O3MnOMgOCaONa2OK2OComment
Mnt/Non ± Pl
 MA13-086-CL0000GLL0Non-rich
 MA13-025-CG00000L00Pl-rich
 MA13-026-CL0LL0LLL0Non = Pl
 MA13-072-CG0LL0LLL0Non-rich
Non + Corr
 MA13-073-CL0GG0GLL0Corr > Non
 MA13-074-CL0GG0GLL0Corr > Non
 MA13-075-CL0000GLL0Corr > Non
I/S + Py
 MA13-087-CL0LL0LLLGMinor Corr
 MA13-034-CL00G0LGLGRect
 MA13-019-CML0GG0GLLG 
 MA13-020-CML0GG0LLLGI/S, Mnt, Rect
 MA13-044-CL0LG0LLLGI/S
MA14-045-CL0GG0LLLGI/S, Mnt, Pl?
 MA13-005-ML0LG0LLLGNon, Rect
 MA13-006-CL0LG0LLLGRect, I/S
 MA13-014-CL0LG0LLLGMont, I/S
 MA13-015-CL0LG0LLLGI/S, Mont, Rect
I/S + Cham
 MA13-088-CL0LL00LLGI/S = Cham
 MA13-090-CL0LG00LLGCham > I/S
 MA13-091-CMG00G00LLGCham > I/S
 MA13-050-CG0GG0LLLGI/S > Cham
 MA13-056-CML0GG0LLLGI/S > Cham
 MA13-057-CL0GG0GGLGCham > I/S
 MA13-058-CL00G0GLLGCham > I/S
 MA13-035-CL0LL00LLGRect = Cham
 MA13-039-CL0LG0GGLGCham > I/S
 MA13-067-CL0GG0LLLGI/S > Cham
 MA13-068-CL0LL0GLLGI/S = Cham
 MA13-069-CL0LL0GLLGI/S = Cham
 MA13-016-CL0GG0LGLGI/S > Cham
Cham + Corr
 MA13-051-CML0LG0GLLGMinor I/S
 MA13-052-CL0GG0GGLGMinor I/S
 MA13-053-CL0GG0GGLGMinor I/S
 MA13-036-CL0LG0GLLGMinor I/S(?)
 MA13-040-CL0LG0GLL0Cham > Corr
 MA13-062-CL0LL0GLL0Corr > Cham
 MA13-063-CL0LL0GLL0Corr > Cham
 MA13-064-CL0GG0GGL0Cham > Corr
 MA13-027-CML0GL00LL0Corr > Cham
Alteration SampleSiO2TiO2Al2O3Fe2O3MnOMgOCaONa2OK2OComment
Mnt/Non ± Pl
 MA13-086-CL0000GLL0Non-rich
 MA13-025-CG00000L00Pl-rich
 MA13-026-CL0LL0LLL0Non = Pl
 MA13-072-CG0LL0LLL0Non-rich
Non + Corr
 MA13-073-CL0GG0GLL0Corr > Non
 MA13-074-CL0GG0GLL0Corr > Non
 MA13-075-CL0000GLL0Corr > Non
I/S + Py
 MA13-087-CL0LL0LLLGMinor Corr
 MA13-034-CL00G0LGLGRect
 MA13-019-CML0GG0GLLG 
 MA13-020-CML0GG0LLLGI/S, Mnt, Rect
 MA13-044-CL0LG0LLLGI/S
MA14-045-CL0GG0LLLGI/S, Mnt, Pl?
 MA13-005-ML0LG0LLLGNon, Rect
 MA13-006-CL0LG0LLLGRect, I/S
 MA13-014-CL0LG0LLLGMont, I/S
 MA13-015-CL0LG0LLLGI/S, Mont, Rect
I/S + Cham
 MA13-088-CL0LL00LLGI/S = Cham
 MA13-090-CL0LG00LLGCham > I/S
 MA13-091-CMG00G00LLGCham > I/S
 MA13-050-CG0GG0LLLGI/S > Cham
 MA13-056-CML0GG0LLLGI/S > Cham
 MA13-057-CL0GG0GGLGCham > I/S
 MA13-058-CL00G0GLLGCham > I/S
 MA13-035-CL0LL00LLGRect = Cham
 MA13-039-CL0LG0GGLGCham > I/S
 MA13-067-CL0GG0LLLGI/S > Cham
 MA13-068-CL0LL0GLLGI/S = Cham
 MA13-069-CL0LL0GLLGI/S = Cham
 MA13-016-CL0GG0LGLGI/S > Cham
Cham + Corr
 MA13-051-CML0LG0GLLGMinor I/S
 MA13-052-CL0GG0GGLGMinor I/S
 MA13-053-CL0GG0GGLGMinor I/S
 MA13-036-CL0LG0GLLGMinor I/S(?)
 MA13-040-CL0LG0GLL0Cham > Corr
 MA13-062-CL0LL0GLL0Corr > Cham
 MA13-063-CL0LL0GLL0Corr > Cham
 MA13-064-CL0GG0GGL0Cham > Corr
 MA13-027-CML0GL00LL0Corr > Cham

Notes: G = gain, L = loss, 0 = no change

Mineral abbreviations: Cham = chamosite, Corr = corrensite, I/S = illite/smectite, Mnt = montmorillonite, Non = nontronite, Pl = plagioclase, Py = pyrite, Rect = rectorite

Fractionation factors after (1) Sheppard and Gilg (1996), (2) Savin and Lee (1988), and (3) Buatier et al. (1995); initial δ18Owater and δDwater calculated for equilibrium temperature of 240°C; hydrogen fractionation factor for chlorite-water (α = 0.954) after Marumo et al. (1980)

Mass dilution effects caused by anhydrite and pyrite precipitation were partially mitigated through handpicking of individual altered clasts, avoiding veins and cement/infill. Despite this procedure, not all of these materials could be omitted, thus contributing CaO, Fe2O3, and S and accounting for some apparent loss of SiO2. The main difference between the low-temperature alteration stage (dominated by illite/smectite + pyrite) and the higher-temperature alteration stage (dominated by chamosite + corrensite) is the mass change of MgO (higher in chamosite) and K2O (higher in illite/smectite), typical of VMS-style alteration in permeable volcaniclastic rocks (e.g., Shikazono et al., 1995; Franklin et al., 2005).

Conditions of Alteration and Mineralization

The high permeability of the volcanic substrate at Tinakula has had a large impact on the composition, size, and morphology of the seafloor mineralization and subseafloor alteration. The conditions of alteration and mineralization have been constrained in this study by fluid inclusion microthermometry, stable isotope data, and dating of the hydrothermal barite.

Fluid inclusion microthermometry

Fluid inclusions were investigated in five samples of hydrothermally altered volcaniclastic rocks recovered in drill hole TKAD001 between 14.8 and 19.3 m b.s.f. Details of the petrography and microthermometry are provided in Appendix Table A4. Fluid inclusions occur in three distinct paragenetic stages: early anhydrite in veins and rimming altered volcanic clasts, sphalerite at the margins of the early anhydrite veins in the deepest samples (MA13-090 and P005), and later infills of anhydrite and barite.

Primary and secondary inclusions in large (up to 1.4 mm), tabular anhydrite crystals are generally 10 to 20 μm, reaching maximum diameters of 60 μm (Fig. 18A). All inclusions are liquid dominated with 1 to 10 vol % vapor (Fig. 18B). Large crystals of barite (up to 1 mm), after anhydrite, contain mostly secondary inclusions, between 4 and 40 μm in size and with a maximum of 10 vol % vapor (Fig. 18C). The inclusions in sphalerite (after anhydrite) include both primary liquid-vapor inclusions up to 50 μm in diameter with 5 to 30 vol % vapor (Fig. 18E) and groups of smaller secondary inclusions (<10 μm) distributed along healed cracks (Fig. 18F). Two inclusions in sphalerite contained a dark solid phase (3–5 vol %) with a tabular habit, likely a sulfide mineral.

Fig. 18.

Transmitted-light in plane-polarized light (TL-PPL) photomicrographs of fluid inclusions from samples from drill hole TKAD001. (A) Single primary fluid inclusion (red circle) in tabular anhydrite next to secondary inclusions formed as a result of crack healing (sample P004_E_g03). (B) Two primary fluid inclusions in anhydrite with a vapor phase (red circle; sample P003_A_g01). (C) A large primary inclusion in barite (red circle; sample P004_R_g01). (D) Large primary inclusion within sphalerite (sample P005_P_g03). (E) Row of secondary fluid inclusions in sphalerite with negative crystal shapes (sample MA13-090_I_g01). (F) Sphalerite crystal containing numerous secondary fluid inclusions, including large one within red circle (sample P005_G_g01). Fluid inclusion groups are outlined in Appendix Table A4.

Fig. 18.

Transmitted-light in plane-polarized light (TL-PPL) photomicrographs of fluid inclusions from samples from drill hole TKAD001. (A) Single primary fluid inclusion (red circle) in tabular anhydrite next to secondary inclusions formed as a result of crack healing (sample P004_E_g03). (B) Two primary fluid inclusions in anhydrite with a vapor phase (red circle; sample P003_A_g01). (C) A large primary inclusion in barite (red circle; sample P004_R_g01). (D) Large primary inclusion within sphalerite (sample P005_P_g03). (E) Row of secondary fluid inclusions in sphalerite with negative crystal shapes (sample MA13-090_I_g01). (F) Sphalerite crystal containing numerous secondary fluid inclusions, including large one within red circle (sample P005_G_g01). Fluid inclusion groups are outlined in Appendix Table A4.

Estimated trapping temperatures of primary fluid inclusions in anhydrite reach 295°C with salinities close to that of seawater (2.8–3.8 wt % NaCl equiv). Inclusions in barite have slightly lower trapping temperatures (median 218°C) with a bimodal distribution from 202° to 220°C and 256° to 266°C (Fig. 19A) and salinities from 2.8 to 3.4 wt % NaCl equiv (Fig. 19B). However, it is known that barite can stretch during heating (e.g., Ulrich and Bodnar, 1988), so these homogenization temperatures may be biased toward higher values. Trapping temperatures of primary fluid inclusions in sphalerite are up to 298°C (Fig. 19C) with salinities of 2.6 to 4.5 wt % NaCl equiv (Fig. 19D).

Fig. 19.

Homogenization temperatures and salinities of fluid inclusions (A-D) from TKAD001 and oxygen isotopes and calculated temperatures (E, F) in clay alteration assemblages. (A) Entrapment temperatures (Te) of inclusions in barite and anhydrite from sample P004 (15.7 m b.s.f.). (B) Salinities of inclusions in barite and anhydrite from sample P004 (15.7 m b.s.f.). (C) Entrapment temperatures (Te) of inclusions in sphalerite from samples MA13-090 (18.5 m b.s.f.) and P005 (19.3 m b.s.f.). (D) Salinities of inclusions in sphalerite from samples MA13-090 (18.5 m b.s.f.) and P005 (19.3 m b.s.f.). (E) Frequency distribution histograms for δ18O values in clay mineral assemblages. (F) Calculated temperatures from oxygen isotopes for the different clay mineral phases. Abbreviations: Anh = anhydrite, Brt = barite, Sp = sphalerite.

Fig. 19.

Homogenization temperatures and salinities of fluid inclusions (A-D) from TKAD001 and oxygen isotopes and calculated temperatures (E, F) in clay alteration assemblages. (A) Entrapment temperatures (Te) of inclusions in barite and anhydrite from sample P004 (15.7 m b.s.f.). (B) Salinities of inclusions in barite and anhydrite from sample P004 (15.7 m b.s.f.). (C) Entrapment temperatures (Te) of inclusions in sphalerite from samples MA13-090 (18.5 m b.s.f.) and P005 (19.3 m b.s.f.). (D) Salinities of inclusions in sphalerite from samples MA13-090 (18.5 m b.s.f.) and P005 (19.3 m b.s.f.). (E) Frequency distribution histograms for δ18O values in clay mineral assemblages. (F) Calculated temperatures from oxygen isotopes for the different clay mineral phases. Abbreviations: Anh = anhydrite, Brt = barite, Sp = sphalerite.

The highest inclusion trapping temperatures are only ~10° to 20°C lower than the minimum temperature of boiling for Tinakula fluids (308°–311°C at 1,070–1,204 m water depth; Fig. 20), and the highest salinities clearly exceed that of seawater. Although coexisting vapor-rich and brine fluid inclusions associated with boiling were not encountered in the investigated minerals, elevated salinities in several inclusions likely reflect mixing between a brine (derived from boiling) and seawater. This strongly suggests that the highest-temperature fluids have boiled at depth and cooled by mixing on ascent to the seafloor. Cooling of deep fluids from a higher temperature at depth (e.g., 350°C for typical end-member hydrothermal fluids) would have resulted in the deposition of most of the Cu-rich sulfides below the seafloor, which could account for the predominance of Zn-rich sulfide in the chimneys (e.g., Hannington et al., 2005).

Fig. 20.

Seawater boiling curve showing maximum temperatures of submarine hydrothermal vents, with examples from different tectonic settings (Hannington et al., 2005) with data from Stoffers et al. (2006). Depth range of Tinakula deposit (gray) corresponds to maximum temperature range of 308° to 311°C. Depths of mid-ocean ridge (MOR) upflow and reaction zones are approximate; C.P. = critical point of seawater.

Fig. 20.

Seawater boiling curve showing maximum temperatures of submarine hydrothermal vents, with examples from different tectonic settings (Hannington et al., 2005) with data from Stoffers et al. (2006). Depth range of Tinakula deposit (gray) corresponds to maximum temperature range of 308° to 311°C. Depths of mid-ocean ridge (MOR) upflow and reaction zones are approximate; C.P. = critical point of seawater.

Oxygen and hydrogen isotopes in clay minerals

Clean separates of monomineralic nontronite (n = 2), chamosite (n = 9), corrensite (n = 2), and rectorite (n = 1) were analyzed for oxygen and hydrogen isotopes, spanning sample depths of 2.8 to 22.5 m b.s.f. (Table 6; Fig. 19E, F). Values of δ18OV-SMOW range from 12.3 to 12.5‰ for nontronite, 2.2 to 8.3‰ for chamosite, 9.6 to 11.1‰ for corrensite, and 5.9‰ for one sample of rectorite. The large range and high δ18O values reflect substantial shifts relative to those of hydrothermally modified seawater (e.g., Huston, 1999). Temperatures were calculated following Savin and Lee (1988) assuming δ18O of the hydrothermal fluid is close to that of typical mid-ocean ridge hydrothermal vent fluids (δ18O = 1.3‰; Shanks et al., 1995; Jean-Baptiste et al., 1997). The calculated temperatures range from 124° to 251°C, with an error of ±10% (Fig. 19F), showing a systematic increase from nontronite (136°–138°C), to chamosite (124°–205°C), corrensite (207°–227°C), and rectorite (256°C), with the highest temperatures close to trapping temperatures of fluid inclusions in the anhydrite and sphalerite. The δ18O value of the hydrothermal fluid in equilibrium with the analyzed clay minerals could have been closer to 5.0‰, if we assume a temperature of formation of 240°C (Table 6), following Sheppard and Gilg (1996). The large range of calculated δ18OH2O values at the inferred temperatures of alteration is consistent with boiling of the hydrothermal fluids (e.g., Huston, 1999; see Fig. 21).

Fig. 21.

δ18O-δD diagram showing isotopic values for the hydrothermal fluid at Tinakula equilibrated at a temperature of 240°C (see Table 6), compared to data from Lau basin (Mottl et al., 2011), PACMANUS deposit in Manus basin (Reeves et al., 2011), and Kuroko ore-forming fluids (Ohmoto, 1996). Fields for seawater and magmatic water (Sheppard, 1986; Taylor, 1992) and global meteoric water line (Craig, 1961) are also shown. Data from Tinakula closely follow adiabatic boiling path from Huston (1999), associated with increasing δ18O values and slightly decreasing δD values with progressive boiling.

Fig. 21.

δ18O-δD diagram showing isotopic values for the hydrothermal fluid at Tinakula equilibrated at a temperature of 240°C (see Table 6), compared to data from Lau basin (Mottl et al., 2011), PACMANUS deposit in Manus basin (Reeves et al., 2011), and Kuroko ore-forming fluids (Ohmoto, 1996). Fields for seawater and magmatic water (Sheppard, 1986; Taylor, 1992) and global meteoric water line (Craig, 1961) are also shown. Data from Tinakula closely follow adiabatic boiling path from Huston (1999), associated with increasing δ18O values and slightly decreasing δD values with progressive boiling.

Table 6.

Oxygen and Hydrogen Isotope Measurements and Calculated Temperatures of Formation of Clay Minerals from Footwall of Tinakula Deposit and Initial Oxygen and Hydrogen Isotope Values for Hydrothermal Fluid Calculated at 240°C

SiteDrill holeSampleDepth (m b.s.f.)MineralδDrock(‰ V-SMOW)δ18Orock(‰ V-SMOW)Fractionation factorsTemperature (°C)δ18Owater, initial (% V-SMOW)δDwater, initial (% V-SMOW)
ABReference
ATKAD001MA13_08716.1–16.2Nontronite–8612.52.554.0511366.9 
ATKAD027MA13_0574.6–4.7Chamosite–534.93.1012.6331644.8–11.3
ATKAD015MA13_0537.2–7.3Chamosite–534.13.1012.6331755.0–10.5
BTKAD025MA13_0395.5–5.6Chamosite–658.33.1012.6331249.2–23.1
BTKAD025MA13_0406.9–7.0Chamosite–522.23.1012.6332053.1–9.8
BTKAD026MA13_0672.9–3Rectorite–695.92.4584.1922560.8 
BTKAD026MA13_0695.6–5.8Chamosite–615.13.1012.6331616.0–18.7
BTKAD012MA13_0622.8–2.9Corrensite–859.64.288.7832272.1 
BTKAD012MA13_0634.3–4.4Corrensite–8011.14.288.7832073.6 
CTKAD032MA13_0756.4–6.5Nontronite–9212.32.554.0511386.7 
ETKAD003POO1211.1–11.2Chamosite–564.63.1012.6331695.4–14.1
ETKAD003POO1416.4–16.5Chamosite–574.03.1012.6331774.9–15.1
FTKAD007POO2211.2–11.3Chamosite–545.83.1012.6331526.7–11.7
FTKAD007POO2322.4–22.5Chamosite–533.33.1012.6331874.2–10.6
SiteDrill holeSampleDepth (m b.s.f.)MineralδDrock(‰ V-SMOW)δ18Orock(‰ V-SMOW)Fractionation factorsTemperature (°C)δ18Owater, initial (% V-SMOW)δDwater, initial (% V-SMOW)
ABReference
ATKAD001MA13_08716.1–16.2Nontronite–8612.52.554.0511366.9 
ATKAD027MA13_0574.6–4.7Chamosite–534.93.1012.6331644.8–11.3
ATKAD015MA13_0537.2–7.3Chamosite–534.13.1012.6331755.0–10.5
BTKAD025MA13_0395.5–5.6Chamosite–658.33.1012.6331249.2–23.1
BTKAD025MA13_0406.9–7.0Chamosite–522.23.1012.6332053.1–9.8
BTKAD026MA13_0672.9–3Rectorite–695.92.4584.1922560.8 
BTKAD026MA13_0695.6–5.8Chamosite–615.13.1012.6331616.0–18.7
BTKAD012MA13_0622.8–2.9Corrensite–859.64.288.7832272.1 
BTKAD012MA13_0634.3–4.4Corrensite–8011.14.288.7832073.6 
CTKAD032MA13_0756.4–6.5Nontronite–9212.32.554.0511386.7 
ETKAD003POO1211.1–11.2Chamosite–564.63.1012.6331695.4–14.1
ETKAD003POO1416.4–16.5Chamosite–574.03.1012.6331774.9–15.1
FTKAD007POO2211.2–11.3Chamosite–545.83.1012.6331526.7–11.7
FTKAD007POO2322.4–22.5Chamosite–533.33.1012.6331874.2–10.6

Fractionation factors after (1) Sheppard and Gilg (1996), (2) Savin and Lee (1988), and (3) Buatier et al. (1995); initial δ18Owater and (δDwater calculated for equilibrium temperature of 240°C; hydrogen fractionation factor for chlorite-water (α = 0.954) after Marumo et al. (1980)

Values of δDv-smow range from –92 to –86‰ for nontronite, –65 to –52‰ for chamosite, –85 to –80‰ for corrensite, and –69‰ in the single sample of rectorite. There is a general increase in δDV_SMOW with depth, consistent with the downhole increase in calculated temperature. Importantly, no accepted hydrogen isotope fractionation factors exist for 2:1 phyllosilicate minerals in order to calculate initial δDwater values. Gilg and Sheppard (1996) proposed approximating the fractionation factors of Al-rich and Fe- and Mg-free smectites using kaolinite factors; however, the smectites in this study are Fe-rich nontronite. Initial δDwater values in equilibrium with chamosite, using a chlorite-water fractionation factor of α = 0.954 from Marumo et al. (1980), are –18.7 to –9.8‰ (Table 6). However, hydrogen isotope fractionation in chlorite is strongly dependent on mineral chemistry (e.g., Marumo et al., 1980; Savin and Lee, 1988), thus adding uncertainty to the initial δDwater values.

High δ18Owater values have been reported for Kuroko-type hydrothermal fluids (up to 10.8‰; Ohmoto, 1983, 1996). The shifts in δ18Owater to lower values at higher temperatures (e.g., minimum of 0.8‰ in equilibrium with rectorite during chalcopyrite deposition) are typical of the stringer zones of ancient VMS deposits (e.g., Cathles, 1993). The 18O-rich fluids at Tinakula may be derived from isotopic exchange with 18O-rich rocks, boiling, or direct contributions from magmatic fluids. At high temperatures, water-rock interactions can shift δ18Owater but not δDwater (Huston, 1999). We interpret the combination of higher δ18Owater and lower δDwater at Tinakula as a result of subseafloor boiling, which is likely because of the shallow water depths there and the fluid inclusion trapping temperatures, which are close to the boiling curve (Fig. 22). The lack of low-pH mineral assemblages and the seawater-dominated S isotope signature (see below) suggest that there was no magmatic contribution to the hydrothermal fluids or only a limited one.

Fig. 22.

Temperature ranges for the formation of hydrothermal alteration minerals from Reyes (1990), with significant overlap shown among minerals. Temperature calculations for the Tinakula deposit using 18O isotopes (blue circles; Table 6) and primary fluid inclusions in anhydrite (orange circles; App. Table A4) are shown.

Fig. 22.

Temperature ranges for the formation of hydrothermal alteration minerals from Reyes (1990), with significant overlap shown among minerals. Temperature calculations for the Tinakula deposit using 18O isotopes (blue circles; Table 6) and primary fluid inclusions in anhydrite (orange circles; App. Table A4) are shown.

Sulfur isotopes

The δ34S values of chalcopyrite range from –0.8 to 2.0‰ and are, on average, slightly lower in the Cu-rich chimneys than in the Zn-rich chimneys (Fig. 23A). Pyrite has δ34S values of –0.4 to 3.1‰ (Fig. 23B), and sphalerite has δ34S of 0.3 to 3.4‰ (Fig. 23C), both minerals having higher values in Zn-rich chimneys. Fe-poor (light-colored) sphalerite has slightly higher δ34S values than Fe-rich (dark-colored) sphalerite from the same sample. However, there is no systematic spatial variation in δ34S of the sulfides along the main NS-trending fissure at Tinakula, suggesting that the sulfur source was uniform. Barite from chimney samples yielded δ34S of 19.3 to 22.5‰ (Fig. 23D). Anhydrite from the altered volcaniclastic rocks has δ34S of 20.6 to 22.2‰, whereas gypsum has δ34S of 21.3 to 21.8‰ (Fig. 23D), all close to that of seawater sulfate (δ34S = 21.24‰; Tostevin et al., 2014), with a slight decrease with depth (App. Table A5).

Fig. 23.

Variations in δ34S (V-CDT) values in minerals from chimneys at the Tinakula deposit: (A) chalcopyrite, (B) pyrite, (C) sphalerite (dark and light), and (D) anhydrite, gypsum, and barite. Rock types are shown on left. Composition of MORB sulfur (MS: δ34S = –0.9 to –0.6‰; Labidi et al., 2012, 2014), island-arc andesite sulfur (ARC: δ34S = 6.00‰; Alt et al., 1993), and seawater sulfate (SS: δ34S = 21.24‰; Tostevin et al., 2014) are shown by dashed lines. Complete S isotope data are in Appendix Table A5. Abbreviations: Brt = barite, Ccp = chalcopyrite, Sp = sphalerite.

Fig. 23.

Variations in δ34S (V-CDT) values in minerals from chimneys at the Tinakula deposit: (A) chalcopyrite, (B) pyrite, (C) sphalerite (dark and light), and (D) anhydrite, gypsum, and barite. Rock types are shown on left. Composition of MORB sulfur (MS: δ34S = –0.9 to –0.6‰; Labidi et al., 2012, 2014), island-arc andesite sulfur (ARC: δ34S = 6.00‰; Alt et al., 1993), and seawater sulfate (SS: δ34S = 21.24‰; Tostevin et al., 2014) are shown by dashed lines. Complete S isotope data are in Appendix Table A5. Abbreviations: Brt = barite, Ccp = chalcopyrite, Sp = sphalerite.

The isotopic composition of the sulfate minerals records predominantly seawater sulfate in pore spaces of the volcaniclastic rocks. Sulfide values are generally close to mid-ocean ridge basalt (MORB) values but lower than those of typical arc lavas (e.g., Alt et al., 1993; Labidi et al., 2012, 2014). This pattern may result from degassing of magmatic volatiles, which can produce H2S with low δ34S values from the disproportionation of SO2 (Gamo et al., 1997; Kusakabe et al., 2000). However, there is no evidence at Tinakula of significant volatile contributions that elsewhere produce advanced argillic alteration and very light sulfide sulfur with negative δ34S values and no evidence of bacterial sulfate reduction that would also produce negative δ34S values (cf. Hannington et al., 2005).

Timing of alteration and mineralization

Results of 226Ra/Ba dating of hydrothermal barite are shown in Table 7 and Figure 24. The calculated ages indicate that hydrothermal venting has been ongoing for at least 5,600 ± 240 years. Approximately 60% of the data yield ages between 792 and 2,369 years. Although there is no obvious spatial variation in the ages (App. Fig. A28), the two oldest samples (>4,000 yr, MA13-011 and MA13-030) were collected from drill core, whereas younger samples were collected primarily from chimneys. Drilling at site D intersected two discrete oxidized intervals that are interpreted as paleoseafloor horizons. The upper paleoseafloor horizon displays lateral continuity over ~20 m east-west and ~25 m north-south, but the lower horizon was not recovered in drill hole TKAD030 (App. Fig. A15); instead, an intact, buried sphalerite-rich chimney was cored. This shows that different parts of the mound were actively venting at different times. The ages of the talus samples confirm that the paleoseafloor horizon on which the talus was deposited is at least 1,000 years old and gives a minimum age for the volcaniclastic cover sequence, which is variably altered and mineralized throughout the Tinakula deposit. Because of the low solubility (Blount, 1977) and high preservation potential of barite, even the oldest barite chimneys are well preserved. One fully intact (inactive) chimney has a calculated age of 3,450 ± 229 years (Table 7). Currently, the system appears to be in a waning phase of hydrothermal activity, with widespread late-stage precipitation of low-temperature amorphous silica and Fe and Mn oxyhydroxides.

Fig. 24.

Distribution of 226Ra ages of hydrothermal barite from sulfide chimney material at the Tinakula massive sulfide deposit. Ages are relative to zero age of sample GMA13-024. Data are in Table 7.

Fig. 24.

Distribution of 226Ra ages of hydrothermal barite from sulfide chimney material at the Tinakula massive sulfide deposit. Ages are relative to zero age of sample GMA13-024. Data are in Table 7.

Table 7.

Results of 226Ra Dating of Hydrothermal Barite from the Tinakula Deposit

SiteSampleCommentLatitude (deg)Longitude (deg)Seafloor depth (m b.s.l.)Drill hole depth (m b.s.f.)Ba (wt %)226 Ra Activity (Bq/kg)226Ra/Ba (Bq/kg×wt %)Age (years)
ROV grab samples
AGMA13-040Inactive chimney–12.193892167.6148111,123.4-9.152,828 ± 1623092,760 ± 233
BGMA13-019ATalus–12.195040167.6151511,100.6-7.343,102 ± 1574232,038 ± 224
BGMA13-021Inactive chimney–12.195794167.6150091,087.5-0.0424 ± 4666988 ± 390
BGMA13-022Talus–12.195428167.6154301,090.3-14.105,162 ± 2523662,369 ± 222
BGMA13-045Talus–12.195959167.6159241,089.7-13.505,153 ± 2473822,273 ± 221
CGMA13-032Inactive chimney–12.197044167.6166831,074.3-0.46188 ± 154112,101 ± 265
CGMA13-061Talus–12.196961167.6161881,081.6-8.483,319 ± 1703912,215 ± 225
DGMA13-010Active chimney–12.197393167.6164951,074.5-2.722,593 ± 209953160 ± 267
DGMA13-024Active chimney–12.197836167.6166201,094.1-3.653,729 ± 1621,0220
DGMA13-029Inactive–12.198487167.6168181,109.7-16.909,187 ± 4075441,457 ± 217
DGMA13-031Talus–12.196567167.6161191,103.3-9.045,322 ± 2485891,273 ± 220
DGMA13-034Active chimney–12.198043167.6169411,108.3-1.851,333 ± 71720807 ± 227
DGMA13-062Talus–12.197878167.6163951,099.7-6.361,728 ± 1002723,057 ± 234
EGMA13-002Inactive chimney–12.198648167.6165131,109.4-14.503,323 ± 1812293,450 ± 229
Drill core samples
BMA13-028Drill hole TKAD013–12.196126167.6155241,090.80.6–0.70.07---
BMA13-029Drill hole TKAD013–12.196126167.6155241,090.82.2–2.34.072,950 ± 149725792 ± 224
BMA13-030Drill hole TKAD013–12.196126167.6155241,090.83.6–3.79.421,337 ± 871424,556 ± 244
BMA13-031Drill hole TKAD013–12.196126167.6155241,090.84.9–5.00.64420 ± 366591,014 ± 274
BMA13-065Drill hole TKAD026–12.196357167.6155491,089.00.5–0.78.193,691 ± 1874511,891 ± 224
DMA13-011Drill hole TKAD018–12.198291167.6169201,110.21.1–1.212.801,158 ± 73905,597 ± 240
DPOO43Drill hole TKAD031–12.198549167.6166431,110.32.9–3.01.21963 ± 56795578 ± 234
EPOO10Drill hole TKAD003–12.201212167.6170341,113.81.0–1.10.97583 ± 405991,234 ± 248
EPOO11Drill hole TKAD003–12.201212167.6170341,113.85.8–5.90.78699 ± 47896304 ± 247
SiteSampleCommentLatitude (deg)Longitude (deg)Seafloor depth (m b.s.l.)Drill hole depth (m b.s.f.)Ba (wt %)226 Ra Activity (Bq/kg)226Ra/Ba (Bq/kg×wt %)Age (years)
ROV grab samples
AGMA13-040Inactive chimney–12.193892167.6148111,123.4-9.152,828 ± 1623092,760 ± 233
BGMA13-019ATalus–12.195040167.6151511,100.6-7.343,102 ± 1574232,038 ± 224
BGMA13-021Inactive chimney–12.195794167.6150091,087.5-0.0424 ± 4666988 ± 390
BGMA13-022Talus–12.195428167.6154301,090.3-14.105,162 ± 2523662,369 ± 222
BGMA13-045Talus–12.195959167.6159241,089.7-13.505,153 ± 2473822,273 ± 221
CGMA13-032Inactive chimney–12.197044167.6166831,074.3-0.46188 ± 154112,101 ± 265
CGMA13-061Talus–12.196961167.6161881,081.6-8.483,319 ± 1703912,215 ± 225
DGMA13-010Active chimney–12.197393167.6164951,074.5-2.722,593 ± 209953160 ± 267
DGMA13-024Active chimney–12.197836167.6166201,094.1-3.653,729 ± 1621,0220
DGMA13-029Inactive–12.198487167.6168181,109.7-16.909,187 ± 4075441,457 ± 217
DGMA13-031Talus–12.196567167.6161191,103.3-9.045,322 ± 2485891,273 ± 220
DGMA13-034Active chimney–12.198043167.6169411,108.3-1.851,333 ± 71720807 ± 227
DGMA13-062Talus–12.197878167.6163951,099.7-6.361,728 ± 1002723,057 ± 234
EGMA13-002Inactive chimney–12.198648167.6165131,109.4-14.503,323 ± 1812293,450 ± 229
Drill core samples
BMA13-028Drill hole TKAD013–12.196126167.6155241,090.80.6–0.70.07---
BMA13-029Drill hole TKAD013–12.196126167.6155241,090.82.2–2.34.072,950 ± 149725792 ± 224
BMA13-030Drill hole TKAD013–12.196126167.6155241,090.83.6–3.79.421,337 ± 871424,556 ± 244
BMA13-031Drill hole TKAD013–12.196126167.6155241,090.84.9–5.00.64420 ± 366591,014 ± 274
BMA13-065Drill hole TKAD026–12.196357167.6155491,089.00.5–0.78.193,691 ± 1874511,891 ± 224
DMA13-011Drill hole TKAD018–12.198291167.6169201,110.21.1–1.212.801,158 ± 73905,597 ± 240
DPOO43Drill hole TKAD031–12.198549167.6166431,110.32.9–3.01.21963 ± 56795578 ± 234
EPOO10Drill hole TKAD003–12.201212167.6170341,113.81.0–1.10.97583 ± 405991,234 ± 248
EPOO11Drill hole TKAD003–12.201212167.6170341,113.85.8–5.90.78699 ± 47896304 ± 247

Bold values indicate the initial 226Ra/Ba ratio of the zero-age chimney; - = not applicable/not determined; sample MA13-028 contained insufficient Ba for analyses

Discussion and Conclusions

Drilling of the Tinakula deposits has revealed widespread hydrothermal alteration within ~10 m of the seafloor in the highly permeable volcaniclastic substrate. The alteration assemblages reported here are consistent with a near-neutral or slightly acidic pH (cf. Reyes, 1990; Corbett and Leach, 1998). Salinities of primary fluid inclusions (2.8–4.5 wt % NaCl equiv) and isotopic data on the sulfide, sulfate, and silicate minerals indicate seawater-dominated fluids, possibly with a component of condensed vapor-phase liquid, and abundant fluid mixing between hydrothermal fluids and ambient seawater. The calculated temperatures of clay formation range from 124° to 256°C, lower than the calculated trapping temperatures of primary fluid inclusions in anhydrite and sphalerite, suggesting that the late sulfate-sulfide veins were likely transporting hotter fluids. The highest-temperature fluids likely boiled at a depth just below the deepest drill holes. Several lines of evidence indicate increasing temperature with depth: (1) the transition from barite near the surface to anhydrite at depth, (2) the progressive transition from smectite to mixed-layer illite/smectite (with increasing proportions of illite) to corrensite to chamosite (e.g., Fig. 11A, B), and decreasing δ18O and increasing δD values of monomineralic clays. Calculated temperatures of formation increase by ~130°C from nontronite → chamosite → corrensite → rectorite. Downhole temperature gradients are as high as 80°C in just over 1 m (e.g., TKAD025 at site B). However, several drill holes show reverse temperature gradients. In particular, data for TKAD026 suggests decreasing temperatures from 256°C at 2.9 m b.s.f. to 161°C at 5.6 m b.s.f., reflecting the control of variable subseafloor permeability on lateral fluid flow and variable subsurface mixing with entrained seawater.

The alteration mineralogy at Tinakula is similar to proximal footwall alteration associated with Kuroko-type VMS deposits (e.g., Shirozu, 1974; Date et al., 1983; Shikazono et al., 1995). In these deposits, the shallow subseafloor strata surrounding the fossil upflow zones is dominated by low-temperature montmorillonite/nontronite, with higher-temperature illite/smectite + pyrite present in the immediate footwall of the massive sulfide and illite/smectite + chamosite and chamosite + corrensite at depth, similar to sites A, C, and D at Tinakula.

Sites A and D are interpreted to be at the center of the highest-temperature upflow, where there is active venting, chalcopyrite-rich chimneys, and extensive formation of illite. The alteration mineralogy at sites B, C, E, and F reflects relatively lower temperatures, consistent with the dominant sphalerite- and barite-rich chimneys.

The alteration assemblages are all associated with Si loss (Fig. 17; Table 5), which may be partially related to dilution associated with anhydrite precipitation at depth, with net mass changes of as much as 49%. Similar loss of Si has been documented in association with phyllosilicate alteration at the PACMANUS hydrothermal system in the Manus back-arc basin (Giorgetti et al., 2006; Paulick and Bach, 2006; Monecke et al., 2007). In contrast, at Tinakula Si gains occur at the tops of drill holes TKAD015, TKAD032, and TKAD033 (App. Figs. A26, A27), representing local silicification near the seafloor in the form of opaline silica filling void spaces. Large gains in Ca near the bottom of drill holes TKAD015, TKAD012, and TKAD018 likely record the deposition of anhydrite at the higher subseafloor temperatures at depth, which could reflect incomplete removal of the anhydrite-rich matrix material during clast separations. Large mass gains in K2O (and LREEs) in the illite/smectite + pyrite assemblages are similar to the observed mass changes in the feeder zones of VMS deposits hosted by felsic volcaniclastic rocks (e.g., Mattabi-type alteration; Morton and Franklin, 1987). Alteration in the immediate upflow zone at Tinakula might also be enhanced by the pyrite precipitation, which generates acid. In the surrounding permeable rocks, drawdown of seawater has added MgO. Intuitively, this process should be associated with cooling, but the higher temperatures of chlorite formation at depth in the Tinakula drill holes, based on δ18O data (Fig. 22), may also have increased the drawdown, resulting in the observed Mg addition.

In permeable volcaniclastic rocks, an important question is how the permeability becomes sealed enough to allow focusing of the hydrothermal fluids to the seafloor. Here, we propose that the pervasive alteration of the volcaniclastic units close to the seafloor coincided with a progressive downward increase in subseafloor temperatures, from initial low temperatures during precipitation of nontronite (136°–138°C), increasing to chamosite (124°–205°C), corrensite (207°–227°C), and eventually rectorite (256°C). At 150°C, the temperatures were high enough to form anhydrite at depths of only a few meters below the seafloor, sealing off some permeability and focusing the fluids into channelized pathways to the seafloor. Continued hydrothermal circulation resulted in overpressuring and fracturing of this relatively impermeable layer, precipitating late sulfate veins with minor sulfide minerals (Fig. 3E), including sphalerite, at temperatures up to 298°C. A similar model has been proposed for the anhydrite-cemented breccias in the Pacmanus deposit in the Manus basin (Binns and Scott, 1993; Binns et al., 2007) and the Suiyo seamount and Iheya North Knoll field in the Izu-Bonin-Mariana arc (Ishibashi et al., 2013, 2015). This process has been investigated using numerical modeling (Kawada and Yoshida, 2010); in ancient deposits, anhydrite is rarely preserved because of its retrograde solubility; the only evidence remaining of an anhydrite seal would be the presence of alteration minerals formed at temperatures that would have caused anhydrite deposition. Chlorite-rich alteration pipes in otherwise permeable volcaniclastic sequences may have been shielded from infiltrating seawater by anhydrite that is no longer present in the rocks. Although rarely preserved in ancient VMS deposits, anhydrite and gypsum are locally abundant in the Kuroko ores (e.g., Shikazono et al., 1983) and may have played an important role in the focusing of hydrothermal upflow in those deposits, as at Tinakula.

Barite and silica, which form close to the seafloor at much lower temperatures than anhydrite, may also have been important in forming an initial cap on the permeable units, promoting lateral flow of fluids and allowing subseafloor temperatures to reach the point of anhydrite precipitation (>150°C; Bischoff and Seyfried, 1978). Early cooling of the fluids and mixing near the seafloor was likely, given the highly permeable volcaniclastic cover. The style of hydrothermal upflow also plays an important role in the nature of the seafloor mineralization. Jamieson et al. (2016) showed a high proportion of seawater mixing in chimneys with predominant dendritic barite; similar diffuse venting is likely responsible for the formation of porous beehive-like structures and unusual bulbous chimneys observed at Tinakula (cf. Koski et al., 1994). This style of venting, which has been inferred for the origin of layered barite-sulfide mounds in the East Diamante caldera (Hein et al., 2014) and at Clark seamount (de Ronde et al., 2014), may be common where seafloor venting occurs through permeable volcaniclastic substrates. Diffuse discharge over a large area is also consistent with the bottom-hugging plume signal originally detected at Tinakula (McConachy et al., 2002).

The Tinakula SMS deposit is the first documented example of a modern SMS deposit developed mainly in bimodal-mafic volcaniclastic rocks. The particular mineralization and alteration styles recognized here are a direct result of the volcanic setting, which is dominated by a permeable substrate rather than the coherent lava flows typical of mid-ocean ridges. Morphologies of the volcanic cones at Tinakula and the volcaniclastic facies successions are consistent with extensive Strombolian-type eruptions, forming a row of cinder cones along a NS-trending fissure. The eruptions likely occurred over days to months (cf. Sigurdsson, 2000), while geothermal activity at Tinakula may have lasted thousands of years. The massive sulfide mineralization may not be closely related temporally to formation of the cinder cones, but instead the hydrothermal systems most likely exploited the same structural corridor. Lateral fluid flow and semiconformable alteration at Tinakula are similar to those described in a class of ancient VMS deposits hosted by volcaniclastic rocks (e.g., Franklin et al., 2005; Galley et al., 2007), such as the Sturgeon Lake VMS district (e.g., Morton and Franklin, 1987; Morton et al., 1991; Hudak et al., 2004; Mueller et al., 2004; Holk et al., 2008). In ancient VMS deposits of this type, semiconformable alteration is commonly observed at depths of several hundred meters below the paleoseafloor, associated with subseafloor deposition of massive sulfide (e.g., Allen et al. 2002; Doyle and Allen, 2003). What is unusual at Tinakula is that hydrothermal fluids affected such a large area so close to the seafloor: ~150,000 m2 at the seafloor and possibly three times as large below the seafloor. The shallow alteration and anhydrite likely acted as an insulating layer for hydrothermal circulation and mineralization at depth, as seen in ancient VMS deposits.

Acknowledgments

The authors gratefully acknowledge Neptune Minerals, Inc., for providing samples, data, and other in-kind contributions that made this study possible. In particular, K. Allen, J. Carpenter, and R. Sharpe are thanked for assistance with sample and data collection. We would also like to thank the following people from GEOMAR: N. Augustin, A. Jegen, and S. Lange for assistance in preparing and analyzing clays, M. Thörner for help with EPMA, and J. Heinze for assistance with XRD. Funding for this project was provided by a Natural Sciences and Engineering Research Council of Canada (NSERC) Alexander Graham Bell Canada Graduate Scholarship (Doctoral) (CGS-D) award to M. Anderson and an NSERC Discovery Grant to M. Hannington. This manuscript was greatly improved by reviews provided by J.F. Slack, R. Allen, G.J. Hudak, and S. Piercey. This is contribution MERC-ME-2019-168 to the modern-ancient crust project of the Canadian Metal Earth program.

REFERENCES

1.
Allen
,
R.L.
,
Weihed
,
P.
, and The Global VMS Research Project Team,
2002
,
Global comparisons of volcanic-associated massive sulfide deposits
, in
Blundell
,
D.J.
,
Neubauer
,
F.
, and
von Quadt
,
A.
, eds.,
The timing and location of major ore deposits in an evolving orogen
 :
Geological Society of London
, Special Publication 204, p.
13
38
.
2.
Alt
,
J.C.
,
Shanks
,
W.C.
III.
, and
Jackson
,
M.C.
,
1993
,
Cycling of sulfur in subduction zones: The geochemistry of sulfur in the Mariana island arc and back-arc trough
:
Earth and Planetary Science Letters
 , v.
119
, p.
477
494
.
3.
Anders
,
M.
,
2016
,
Fluid inclusions and mineralogy at the seafloor massive sulfide deposit, New Hebrides arc/backarc, Solomon Islands
: Unpublished M.Sc. thesis,
Kiel, Germany
,
Christian Albrechts University of Kiel
,
96
p.
4.
Anderson
,
M.O.
,
2018
,
Relationships between tectonics, volcanism, and hydrothermal venting in the New Hebrides and Mariana back-arc basins, western Pacific
: Unpublished Ph.D. thesis,
Ottawa, Canada
,
University of Ottawa
,
1030
p.
5.
Anderson
,
M.O.
,
Hannington
,
M.D.
,
Haase
,
K.M.
,
Schwarz-Schampera
,
U.
,
Augustin
,
N.
,
McConachy
,
T.F.
, and
Allen
,
K.
,
2016
,
Tectonic focusing of voluminous basaltic eruptions in magma-deficient backarc rifts
:
Earth and Planetary Science Letters
 , v.
440
, p.
43
55
.
6.
Auzende
,
J.-M.
,
Pelletier
,
B.
, and
Eissen
,
J.-P
.,
1995
,
The North Fiji basin: Geology, structure, and geodynamic evolution
, in
Taylor
,
B.
, ed.,
Backarc basins: Tectonics and magmatism
 :
New York
,
Plenum Press
, p.
139
175
.
7.
Barrett
,
T.J.
, and
MacLean
,
W.H.
,
1999
,
Volcanic sequences, litho-geochemistry, and hydrothermal alteration in some bimodal volcanic-associated massive sulfide systems
:
Reviews in Economic Geology
 , v.
8
, p.
101
131
.
8.
Barrie
,
C.T.
, and
Hannington
,
M.D.
,
1999
,
Classification of volcanic-associated massive sulfide deposits based on host-rock composition
:
Reviews in Economic Geology
 , v.
8
, p.
1
11
.
9.
Beaulieu
,
S.E.
, and
Szafranski
,
K.
,
2018
,
InterRidge Global database of active submarine hydrothermal vent fields
:
Prepared for InterRidge, Version 3.4
 , http://vents-data.interridge.org.
10.
Bergeot
,
N.
,
Bouin
,
M.N.
,
Diament
,
M.
,
Pelletier
,
B.
,
Régnier
,
M.
,
Calmant
,
S.
, and
Ballu
,
V.
,
2009
,
Horizontal and vertical interseismic velocity fields in the Vanuatu subduction zone from GPS measurements: Evidence for a central Vanuatu locked zone
:
Journal of Geophysical Research
 , v.
114
, p.
1
20
.
11.
Biagioni
,
C.
,
Dini
,
A.
,
Orlandi
,
P.
,
Moëlo
,
Y.
,
Pasero
,
M.
, and
Zaccarini
,
F.
,
2016
,
Lead-antimony sulfosalts from Tuscany (Italy). XX. Members of the jordanite-geocronite series from the Pollone mine, Valdicastello Carducci: Occurrence and crystal structures
:
Minerals
 , v.
6
, p.
1
15
.
12.
Binns
,
R.A.
, and
Scott
,
S.D.
,
1993
,
Actively forming polymetallic sulfide deposits associated with felsic volcanic rocks in the eastern Manus back-arc basin, Papua New Guinea
:
Economic Geology
 , v.
88
, p.
2226
2236
.
13.
Binns
,
R.A.
,
Barriga
,
F.J.A.S.
, and
Miller
,
D.J.
,
2007
,
Leg 193 synthesis: Anatomy of an active felsic-hosted hydrothermal system, eastern Manus basin, Papua New Guinea
, in
Barriga
,
F.J.A.S.
,
Binns
,
R.A.
,
Miller
,
D.J.
, and
Herzig
,
P.M.
, eds.,
Proceedings of the Ocean Drilling program, scientific results
, v.
193
:
College Station, Texas
,
Ocean Drilling Program
, p.
1
17
.
14.
Bischoff
,
J.L.
, and
Seyfried
,
W.E.
Jr.
,
1978
,
Hydrothermal chemistry of seawater from 25°C to 350°C
:
American Journal of Science
 , v.
278
, p.
838
860
.
15.
Blount
,
C.W.
,
1977
,
Barite solubilities and thermodynamic quantities up to 300°C and 1400 bars
:
American Mineralogist
 , v.
62
, p.
942
957
.
16.
Bodnar
,
R.
, and
Vityk
,
M.
,
1994
,
Interpretation of microthermometric data for H2O-NaCl fluid inclusions
, in
De Vivo
,
B.
, and
Frezzotti
,
M.
, eds.,
Fluid inclusions in minerals: Methods and applications
 :
Blacksburg, Virginia
,
Virginia Polytechnic and State University
, p.
117
130
.
17.
Boynton
,
W.V.
,
1984
, Cosmochemistry of the rare earth elements: Meteorite studies, in
Henderson
,
P.
, ed.,
Rare earth element geochemistry
 :
Amsterdam
,
Elsevier
, p.
63
114
.
18.
Brown
,
P.E.
,
1989
,
FLINCOR: A computer program for the reduction of and investigation of fluid inclusion data
:
American Mineralogist
 , v.
74
, p.
1390
1393
.
19.
Buatier
,
M.D.
,
Früh-Green
,
G.L.
, and
Karpoff
,
A.M.
,
1995
,
Mechanism of Mg-phyllosilicate formation in a hydrothermal system at a sedimented ridge (Middle Valley, Juan de Fuca)
:
Contributions to Mineralogy and Petrology
 , v.
122
, p.
134
151
.
20.
Canfield
,
D.E.
,
Raiswell
,
R.
,
Westrich
,
J.T.
,
Reaves
,
C.M.
, and
Berner
,
R.A.
,
1986
,
The use of chromium reduction in the analysis of reduced inorganic sulfur in sediments and shales
:
Chemical Geology
 , v.
54
, p.
149
155
.
21.
Cashman
,
K.V.
,
2004
,
Volatile controls on magma ascent and eruption
, in
Sparkes
,
R.S.J.
, and
Hawkesworth
,
C.J.
, eds.,
The state of the planet: Frontiers and challenges in geophysics
 :
Washington, D.C.
,
American Geophysical Union
, p.
109
124
.
22.
Cathles
,
L.M.
,
1993
,
Oxygen isotope alteration in the Noranda mining district, Abitibi greenstone belt, Quebec
:
Economic Geology
 , v.
88
, p.
1483
1511
.
23.
Clayton
,
R.N.
, and
Mayeda
T.K.
,
1963
,
The use of bromine pentafluoride in the extraction of oxygen from oxides and silicates for isotopic analysis
:
Geochimica et Cosmochimica Acta
 , v.
27
, p.
43
52
.
24.
Cook
,
N.
,
Ciobanu
,
C.L.
,
Pring
,
A.
,
Skinner
,
W.
,
Shimizu
,
M.
,
Danyushevsky
,
L.
,
Saini-Eidukat
,
B.
, and
Melcher
,
F.
,
2009
,
Trace and minor elements in sphalerite: A LA-ICPMS study
:
Geochimica et Cosmochimica Acta
 , v.
73
, p.
4761
4791
.
25.
Corbett
,
G.J.
, and
Leach
,
T.M.
,
1998
,
Southwest Pacific rim gold-copper systems: Structure, alteration and mineralization
:
Society of Economic Geologists
 , Special Publication 6,
238
p.
26.
Craig
,
H.
,
1961
,
Standard for reporting concentrations of deuterium and oxygen-18 in natural waters
:
Science
 , v.
133
, p.
1833
1834
.
27.
Date
,
J.
,
Watanabe
,
Y
, and
Saeki
,
Y.
,
1983
,
Zonal alteration around the Fukazawa Kuroko deposits, Akita Prefecture, northern Japan
:
Economic Geology Monograph
 
5
, p.
365
386
.
28.
de Ronde
,
C.E.J.
,
Walker
,
S.L.
,
Ditchburn
,
R.G.
,
Caratori Tontini
,
F.
,
Hannington
,
M.D.
,
Merle
,
S.G.
,
Timm
,
C.
,
Handler
,
M.R.
,
Wysoczanski
,
R.J.
,
Dekov
,
V.M.
,
Kamenov
,
G.D.
,
Baker
,
E.T.
,
Embley
,
R.W.
,
Lupton
,
J.E.
, and
Stoffers
,
P.
,
2014
,
The anatomy of a buried submarine hydrothermal system, Clark volcano, Kermadec arc, New Zealand
:
Economic Geology
 , v.
109
, p.
2261
2292
.
29.
Ditchburn
,
R.
,
Graham
,
I.
,
Barry
,
B.
, and
de Ronde
,
C.E.J.
,
2004
,
Uranium series disequilibrium dating of black smoker chimneys
:
New Zealand Science Reviews
 , v.
61
, p.
54
56
.
30.
Ditchburn
,
R.G.
, and
de Ronde
,
C.E.J.
,
2017
,
Evidence for remobilization of barite affecting radiometric dating using 228Ra, 228Th, and 226Ra/Ba values: Implications for the evolution of sea-floor volcanogenic massive sulfides
:
Economic Geology
 , v.
112
, p.
1231
1245
.
31.
Dogramaci
,
S.S.
,
Herczeg
,
A.L.
,
Schiff
,
S.L.
, and
Bone
,
Y.
,
2001
,
Controls on δ34S and δ18O of dissolved SO4 in aquifers of the Murray basin, Australia, and their use as indicators of flow processes
:
Applied Geochemistry
 , v.
16
, p.
475
488
.
32.
Doyle
,
M.G.
, and
Allen
,
R.L.
,
2003
,
Subsea-floor replacement in volcanic-hosted massive sulfide deposits
:
Ore Geology Reviews
 , v.
23
, p.
183
222
.
33.
Franklin
,
J.M.
,
Kasarda
,
J.
, and
Poulsen
,
K.H.
,
1975
,
Petrology and chemistry of the alteration zone of the Mattabi massive sulfide deposit
:
Economic Geology
 , v.
70
, p.
63
79
.
34.
Franklin
,
J.M.
,
Lydon
,
J.W.
, and
Sangster
,
D.F.
,
1981
,
Volcanic-associated massive sulfide deposits
:
Economic Geology 75th Anniversary Volume
 , p.
485
627
.
35.
Franklin
,
J.M.
,
Gibson
,
H.L.
,
Jonasson
,
I.R.
, and
Galley
,
A.G.
,
2005
,
Volcanogenic massive sulfide deposits
:
Economic Geology 100th Anniversary Volume
 , p.
523
560
.
36.
Galley
,
A.G.
,
1993
,
Characteristics of semi-conformable alteration zones associated with volcanogenic massive sulphide districts
:
Journal of Geochemical Exploration
 , v.
48
, p
175
199
.
37.
Galley
,
A.G.
,
Watkinson
,
D.H.
,
Jonasson
,
I.R.
, and
Riverin
,
G.
,
1995
,
The subsea-floor formation of volcanic-hosted massive sulfide; evidence from the Ansil deposit, Rouyn-Noranda, Canada
:
Economic Geology
 , v.
90
, p.
2006
2017
.
38.
Galley
,
A.G.
,
Hannington
,
M.D.
, and
Jonasson
,
I.
,
2007
,
Volcanogenic massive sulfide deposits
, in
Goodfellow
,
W.D.
, ed.,
Mineral deposits of Canada: A synthesis of major deposit-types, district metallogeny, the evolution of geological provinces, and exploration methods
 :
Mineral Deposits Division, Geological Association of Canada
, Special Publication 5, p.
141
161
.
39.
Gamo
,
T.
,
Okamura
,
K.
,
Charlou
,
J.-L.
,
Urabe
,
T.
,
Auzende
,
J.-M.
,
Ishibashi
,
J.
,
Shitashima
,
K.
,
Chiba
,
H.
, and
Shipboard Scientific
Party
,
1997
,
Acidic and sulfate-rich hydrothermal fluids from the Manus back-arc basin, Papua New Guinea
:
Geology
 , v.
25
, p.
139
142
.
40.
Gibson
,
H.L.
, and
Galley
,
A.G.
,
2007
,
Volcanogenic massive sulfide deposits of the Archean, Noranda district, Quebec
, in
Goodfellow
,
W.D.
, ed.,
Mineral deposits of Canada: A synthesis of major deposit-types, district metallogeny, the evolution of geological provinces, and exploration methods
 :
Mineral Deposits Division, Geological Association of Canada
, Special Publication 5, p.
533
552
.
41.
Gibson
,
H.L.
,
Watkinson
,
D.H.
, and
Comba
,
C.D.A.
,
1983
,
Silicification; hydrothermal alteration in an Archean geothermal system within the Amulet Rhyolite Formation, Noranda, Quebec
:
Economic Geology
 , v.
78
, p.
954
971
.
42.
Gibson
,
H.L.
,
Watkinson
,
D.H.
,
Watkins
,
J.J.
,
Labrie
,
M.
, and
Doiron
,
G.
,
1993
,
Volcanological reconstruction of the Corbet breccia pile and Cu-Zn massive sulphide deposit, Noranda, Quebec
:
Exploration and Mining Geology
 , v.
2
, p.
1
16
.
43.
Gifkins
,
C.C.
, and
Allen
,
R.L.
,
2001
,
Textural and chemical characteristics of diagenetic and hydrothermal alteration in glassy volcanic rocks: Examples from the Mount Read Volcanics, Tasmania
:
Economic Geology
 , v.
96
, p.
973
1002
.
44.
Gilg
,
H.A.
, and
Sheppard
,
S.M.F.
,
1996
,
Hydrogen isotope fractionation between kaolinite and water revisited
:
Geochimica et Cosmochimica Acta
 , v.
60
, p.
529
533
.
45.
Gill
,
J.
,
Torssander
,
P.
,
Lapierre
,
H.
,
Taylor
,
R.
,
Kaiho
,
K.
,
Koyama
,
M.
,
Kusakabe
,
M.
,
Aitchison
,
J.
,
Cisowski
,
S.
,
Dadey
,
K.
,
Fujioka
,
K.
,
Klaus
,
A.
,
Lovell
,
M.
,
Marsaglia
,
K.
,
Pezard
,
P.
,
Taylor
,
B.
, and
Tazaki
,
K.
,
1990
,
Explosive deep water basalt in the Sumisu backarc rift
:
Science
 , v.
248
, p.
1214
1217
.
46.
Giorgetti
,
G.
,
Monecke
,
T.
,
Kleeberg
,
R.
, and
Hannington
,
M.D.
,
2006
,
Low-temperature hydrothermal alteration of silicic glass at the Pacmanus hydrothermal vent field, Manus basin, an XRD, SEM, and AEM-TEM study
:
Clays and Clay Minerals
 , v.
54
, p.
240
251
.
47.
Hannington
,
M.D.
,
de Ronde
,
C.E.J.
, and
Petersen
,
S.
,
2005
,
Sea-floor tectonics and submarine hydrothermal systems
:
Economic Geology 100th Anniversary Volume
 , p.
111
141
.
48.
Hannington
,
M.D.
,
Jamieson
,
J.
,
Monecke
,
T.
, and
Petersen
,
S.
,
2010
,
Modern sea-floor massive sulfides and base metal resources: Toward an estimate of global sea-floor massive sulfide potential
:
Society of Economic Geologists Special Publication 15
 , v.
2
, p.
317
338
.
49.
Hannington
,
M.D.
,
Jamieson
,
J.
,
Monecke
,
T.
,
Petersen
,
S.
, and
Beaulieu
,
S.
,
2011
,
The abundance of seafloor massive sulfide deposits
:
Geology
 , v.
39
, p.
1155
1158
.
50.
Hein
,
J.R.
,
de Ronde
,
C.E.J.
,
Koski
,
R.A.
,
Ditchburn
,
R.G.
,
Mizell
,
K.
,
Tamura
,
Y.
,
Stern
,
R.J.
,
Conrad
,
T.A.
,
Ishizuka
,
O.
, and
Leybourne
,
M.I.
,
2014
,
Layered hydrothermal barite-sulfide mound field, East Diamante caldera, Mariana volcanic arc
:
Economic Geology
 , v.
109
, p.
2179
2206
.
51.
Holk
,
G.J.
,
Taylor
,
B.E.
, and
Galley
,
A.G.
,
2008
,
Oxygen isotope mapping of the Archean Sturgeon Lake caldera complex and VMS-related hydrothermal system, northwestern Ontario, Canada
:
Mineralium Deposita
 , v.
43
, p.
623
640
.
52.
Hudak
,
G.J.
,
Morton
,
R.L.
,
Franklin
,
J.M.
, and
Peterson
,
D.M.
,
2004
,
Morphology, distribution, and estimated eruption volumes for intracaldera tuffs associated with volcanic-hosted massive sulfide deposits in the Archean Sturgeon Lake caldera complex, northwestern Ontario
:
American Geophysical Union, Geophysical Monograph
 
140
, p.
345
360
.
53.
Huston
,
D.L.
,
1999
,
Stable isotopes and their significance for understanding the genesis of volcanic-hosted massive sulfide deposits: A review
:
Reviews in Economic Geology
 , v.
8
, p.
157
179
.
54.
Inoue
,
A.
,
Watanabe
,
T.
,
Kohyama
,
N.
, and
Brusewitz
,
A.M.
,
1990
,
Characterization of illitization of smectite in bentonite beds at Kinnekulle, Sweden
:
Clays and Clay Minerals
 , v.
38
, p.
241
249
.
55.
Ishibashi
,
J.-I.
,
Miyoshi
,
Y.
,
Inoue
,
H.
,
Yeats
,
C.
,
Hollis
,
S.P.
,
Corona
,
J.C.
,
Bowden
,
S.
,
Yang
,
S.
,
Southam
,
G.
,
Masaki
,
Y.
,
Hartnet
,
H.
, and
IODP Expedition 331 Scientists
,
2013
,
Subseafloor structure of a submarine hydrothermal system within volcaniclastic sediments: A modern analogue for “Kuroko-type” VMS deposits [ext. abs.]
:
Society for Geology Applied to Mineral Deposits (SGA) Biennial Meeting, 12th, Uppsala, Sweden, 2013
, Extended Abstracts, p.
542
544
.
56.
Ishibashi
,
J.-I.
,
Ikegami
,
F.
,
Tsuji
,
T.
, and
Urabe
,
T.
,
2015
,
Hydrothermal activity in the Okinawa trough back-arc basin: Geological background and hydrothermal mineralization
, in
Ishibashi
,
J.-I.
,
Okino
,
K.
, and
Sunamura
,
M.
, eds.,
Subseafloor biosphere linked to hydrothermal systems
 :
Tokyo
,
Springer
, p.
337
359
.
57.
Jamieson
,
J.W.
,
Hannington
,
M.D.
,
Clague
,
D.A.
,
Kelley
,
D.S.
,
Delaney
,
J.R.
,
Holden
,
J.F.
,
Tivey
,
M.K.
, and
Kimpe
,
L.E.
,
2013
,
Sulfide geochronology along the Endeavour segment of the Juan de Fuca Ridge
:
Geochemistry, Geophysics, Geosystems
 , v.
14
, p.
2084
2099
.
58.
Jamieson
,
J.W.
,
Hannington
,
M.D.
,
Tivey
,
M.K.
,
Hansteen
,
T.
,
Williamson
,
N. M.-B.
,
Stewart
,
M.
,
Fietzke
,
J.
,
Butterfield
,
D.
,
Frische
,
M.
,
Allen
,
L.
,
Cousens
,
B.
, and
Langer
,
J.
,
2016
,
Precipitation and growth of barite within hydrothermal vent deposits from the Endeavour segment, Juan de Fuca Ridge
:
Geochimica et Cosmochimica Acta
 , v.
173
, p.
64
85
.
59.
Jean-Baptiste
P.
,
Charlou
,
J.L.
, and
Stievenard
,
M.
,
1997
,
Oxygen isotope study of mid-ocean ridge hydrothermal fluids: Implication for the oxygen-18 budget of the oceans
:
Geochimica et Cosmochimica Acta
 , v.
61
, p.
2669
2677
.
60.
Jinfeng
,
Y.
,
Xing
,
L.
,
Liang
,
L.
,
Pan
,
J.
, and
Meng
,
T.
,
2014
,
Application of short-wave infrared (SWIR) spectroscopy in quantitative estimation of clay mineral contents
:
IOP Conference Series
,
Earth and Environmental Science
 , v.
17
,
6
p.
61.
Johnson
,
D.P.
,
Maillet
,
P.C.
, and
Rice
,
R.C.
,
1992
,
The seabed morphology of the Hazel Holme fracture zone and the New Hebrides arc, northern Vanuatu-Solomon Islands region
:
South Pacific Applied Geoscience Commission
 , Technical Report 138,
28
p.
62.
Kawada
,
Y.
, and
Yoshida
,
S.
,
2010
,
Formation of a hydrothermal reservoir due to anhydrite precipitation in an arc volcano hydrothermal system
:
Journal of Geophysical Research
 , v.
115
,
B11106
.
63.
Kerr
,
D.J.
, and
Gibson
,
H.L.
,
1993
,
A comparison of the Horne volcanogenic massive sulfide deposit and intracauldron deposits of the Mine sequence, Noranda, Quebec
:
Economic Geology
 , v.
88
, p.
1419
1442
.
64.
King
,
E.M.
,
Valley
,
J.W.
, and
Davis
,
D.W.
,
2000
,
Oxygen isotope evolution of volcanic rocks at the Sturgeon Lake volcanic complex, Ontario
:
Canadian Journal of Earth Science
 , v.
37
, p.
39
50
.
65.
Koski
,
R.A.
,
Jonasson
,
I.R.
,
Kadko
,
D.C.
,
Smith
,
V.K.
, and
Wong
,
F.L.
,
1994
,
Compositions, growth mechanisms, and temporal relations of hydrothermal sulfide-sulfate-silica chimneys at the northern Cleft segment, Juan de Fuca Ridge
:
Journal of Geophysical Research
 , v.
99
, p.
4813
4832
.
66.
Kroenke
,
L.W.
,
1995
,
A morphotectonic interpretation of SPOACMAPS 1:500,000 charts: Central Solomon Islands-Southern Tuvalu
:
South Pacific Applied Geoscience Commission
 , Technical Report 220,
42
p.
67.
Kusakabe
,
M.
,
Komoda
,
Y.
,
Takano
,
B.
, and
Abiko
,
T.
,
2000
,
Sulfur isotopic effects in the disproportionation reaction of sulfur dioxide in hydrothermal fluids: Implications for the δ34S variations of dissolved bisulfate and elemental sulfur from active crater lakes
:
Journal of Volcanology and Geothermal Research
 , v.
97
, p.
287
307
.
68.
Labidi
,
J.
,
Cartigny
,
P.
,
Birck
,
J.
,
Assayag
,
N.
, and
Bourrand
,
J.
,
2012
,
Determination of multiple sulfur isotopes in glasses: A reappraisal of the MORB δ34S
:
Chemical Geology
 , v.
334
, p.
189
198
.
69.
Labidi
,
J.
,
Cartigny
,
P.
,
Hamelin
,
C.
,
Moreira
,
M.
, and
Dosso
,
L.
,
2014
,
Sulfur isotope budget (32S, 33S, 34S and 36S) in Pacific-Antarctic ridge basalts: A record of mantle heterogeneity and hydrothermal sulfide assimilation
:
Geochimica et Cosmochimica Acta
 , v.
133
, p.
47
67
.
70.
Large
,
R.R.
,
Allen
,
R.L.
,
Blake
,
M.D.
, and
Herrmann
,
W.
,
2001
,
Hydrothermal alteration and volatile element halos for the Rosebery K Lens volcanic-hosted massive sulfide deposit, western Tasmania
:
Economic Geology
 , v.
96
, p.
1055
1072
.
71.
MacDonald
,
R.
,
Barrett
,
W.J.
, and
Sherlock
,
R.L.
,
1996
,
Geology and lithogeo-chemistry at the Hidden Creek massive sulfide deposit, Anyonx, west-central British Columbia
:
Exploration and Mining Geology
 , v.
5
, p.
369
398
.
72.
Maclean
,
W.H.
,
1990
,
Mass change calculation in altered rock series
:
Mineralium Deposita
 , v.
25
, p.
44
49
.
73.
Maclean
,
W.H.
, and
Barrett
,
T.J.
,
1993
,
Lithogeochemical techniques using immobile elements
:
Journal of Geochemical Exploration
 , v.
48
, p.
109
133
.
74.
Maillet
,
P.
,
Ruellan
,
E.
,
Gérard
,
M.
,
Person
,
A.
,
Bellon
,
H.
,
Cotton
,
J.
,
Joron
,
J.-L.
,
Nakada
,
S.
, and
Price
,
R.C.
,
1995
,
Tectonics, magmatism, and evolution of the New Hebrides back-arc troughs (southwest Pacific)
, in
Taylor
,
B.
, ed.,
Backarc basins: Tectonics and magmatism
 :
New York
,
Plenum
, p.
177
235
.
75.
Marumo
,
K.
,
Nagasawa
,
K.
, and
Kuroda
,
Y.
,
1980
,
Mineralogy and hydrogen isotope geochemistry of clay minerals in the Ohnuma geothermal area, northeastern Japan
:
Earth and Planetary Science Letters
 , v.
47
, p.
255
262
.
76.
McConachy
,
T.F.
,
Yeats
,
C.J.
, and
Shipboard Scientific
Party
,
2002
,
Final cruise FR03-2002, SOLAVENTS-2002, Solomons Australia vents expedition aboard the RV Franklin, 26 March–21 April 2002
:
Commonwealth Scientific and Industrial Research Organisation (CSIRO)
 , Exploration and Mining Report 1026F,
440
p.
77.
Monecke
,
T.
,
Giorgetti
,
G.
,
Scholtysek
,
O.
,
Kleeberg
,
R.
,
Götze
,
J.
,
Hannington
,
M.D.
, and
Petersen
,
S.
,
2007
,
Textural and mineralogical changes associated with the incipient hydrothermal alteration of glassy dacite at the submarine PACMANUS hydrothermal system, eastern Manus basin
:
Journal of Volcanology and Geothermal Research
 , v.
160
, p.
23
41
.
78.
Monecke
,
T.
,
Petersen
,
S.
,
Hannington
,
M.D.
,
Grant
,
H.
, and
Samson
,
I.
,
2016
,
The minor element endowment of modern sea-floor massive sulfides and comparison with deposits hosted in ancient volcanic successions
:
Reviews in Economic Geology
 , v.
18
, p.
245
306
.
79.
Monjaret
,
M.C.
,
Bellon
,
H.
, and
Maillet
,
P.
,
1991
,
Magmatism of the troughs behind the New Hebrides island arc (RV Jean Charcot SEAPSO 2 cruise): K-Ar geochronology and petrology
:
Journal of Volcanology and Geothermal Research
 , v.
46
, p.
265
280
.
80.
Moore
,
D.M.
, and
Reynolds
,
C.R.
Jr.
,
1997
,
X-ray diffraction and the identification and analysis of clay minerals, 2nd ed
:
New York
,
Oxford University Press
,
400
p.
81.
Morton
,
R.L.
, and
Franklin
,
J.M.
,
1987
,
Two-fold classification of Archean volcanic-associated massive sulfide deposits
:
Economic Geology
 , v.
82
, p.
1057
1063
.
82.
Morton
,
R.L.
,
Walker
,
J.S.
,
Hudak
,
G.J.
, and
Franklin
,
J.M.
,
1991
,
The early development of an Archean submarine caldera complex with emphasis on the Mattabi ash-flow tuff and its relationship to the Mattabi massive sulfide deposit
:
Economic Geology
 , v.
86
, p.
1002
1011
.
83.
Mottl
,
M.
,
Seewald
,
J.
,
Wheat
,
C.
,
Tivey
,
M.
,
Michael
,
P.
,
Proskurowski
,
G.
,
McCollom
,
T.M.
,
Reeves
,
E.
,
Sharkey
,
J.
, and
You
,
C.F.
,
2011
,
Chemistry of hot springs along the eastern Lau spreading center
:
Geochimica et Cosmochimica Acta
 , v.
75
, p.
1013
1038
.
84.
Mueller
,
W.
, and
White
,
J.D.L.
,
1992
,
Felsic fire-fountaining beneath Archean seas: Pyroclastic deposits of the 2730 Ma Hunter mine group, Quebec, Canada
:
Journal of Volcanology and Geothermal Research
 , v.
54
, p.
117
134
.
85.
Mueller
,
W.
,
White
,
J.D.L.
,
Stix
,
J.
, and
Hudak
,
G.J.
,
2004
,
Archean calderas
, in
Ericksson
,
P.
,
Altermann
,
W.
,
Nelson
,
D.
,
Mueller
,
W.
,
Catuneanu
,
O.
, and
Strand
,
K.
, eds.,
The Precambrian Earth: Developments in Precambrian geology
 , v.
12
:
Elsevier
, p.
345
356
.
86.
Ohmoto
,
H.
,
1983
,
Geologic setting of the Kuroko deposits, Japan: Part I. Geologic history of the Green Tuff region
:
Economic Geology Monograph
 
5
, p.
9
24
.
87.
Ohmoto
,
H.
,
1996
,
Formation of volcanogenic massive sulfide deposits: The Kuroko perspective
:
Ore Geology Reviews
 , v.
10
, p.
135
177
.
88.
Ono
,
S.
,
Wing
,
B.A.
,
Johnston
,
D.
,
Farquhar
,
J.
, and
Rumble
,
D.
,
2006
,
Mass-dependent fractionation of quadruple sulfur isotope system as a new tracer of sulfur biogeochemical cycles
:
Geochimica et Cosmochimica Acta
 , v.
70
, p.
2238
2252
.
89.
Paulick
,
H.
, and
Bach
,
W.
,
2006
,
Phyllosilicate alteration mineral assemblages in the active subsea-floor Pacmanus hydrothermal system, Papua New Guinea, ODP Leg 193
:
Economic Geology
 , v.
101
, p.
633
650
.
90.
Pelletier
,
B.
,
Lafoy
,
Y.
, and
Missegue
,
F.
,
1993
,
Morphostructure and magnetic fabric of the northwestern North Fiji basin
:
Geophysical Research Letters
 , v.
20
, p.
1151
1154
.
91.
Pelletier
,
B.
,
Calmant
,
S.
, and
Pillet
,
R.
,
1998
,
Current tectonics of the Tonga-New Hebrides region
:
Earth and Planetary Science Letters
 , v.
164
, p.
263
276
.
92.
Piercey
,
S.J.
,
Squires
,
G.C.
, and
Brace
,
T.D.
,
2014
,
Lithostratigraphic, hydrothermal, and tectonic setting of the Boundary volcanogenic massive sulfide deposit, Newfoundland Appalachians, Canada: Formation by subseafloor replacement in a Cambrian rifted arc
:
Economic Geology
 , v.
109
, p.
661
687
.
93.
Récy
,
J.
,
Charvis
,
P
,
Ruellan
,
E.
,
Monjaret
,
M.-C.
,
Gérard
,
M.
,
Auclair
,
G.
,
Baldassari
,
C.
,
Boirat
,
J.-M.
,
Brown
,
G.R.
,
Butscher
,
J.
,
Collot
,
J.Y.
,
Daniel
,
J.
,
Louat
,
R.
,
Monzier
,
M.
, and
Pontoise
,
B.
,
1986
,
Tectonique et volcanisme sous-marin à l’arrière de l’arc des Nouvelles-Hébrides (Vanuatu): Résultats préliminaires de la campagne SEAPSO (leg 2) du N.O. Jean Charcot
:
Comptes Rendus de d’Académie des Sciences
 , v.
303
, p.
685
690
.
94.
Reeves
,
E.P.
,
Seewald
,
J.S.
,
Saccocia
,
P.
,
Bach
,
W.
,
Craddock
,
P.R.
,
Shanks
,
W.C.
,
Sylva
,
S.P.
,
Walsh
,
E.
,
Pichler
,
T.
, and
Rosner
,
M.
,
2011
,
Geochemistry of hydrothermal fluids from the PACMANUS, northeast Pual and Vienna Woods hydrothermal fields, Manus basin, Papua New Guinea
:
Geochimica et Cosmochimica Acta
 , v.
75
, p.
1088
1123
.
95.
Reyes
,
A.G.
,
1990
,
Petrology of Philippine geothermal systems and the application of alteration mineralogy to their assessment
:
Journal of Volcanology and Geothermal Research
 , v.
42
, p.
279
309
.
96.
Rice
,
C.A.
,
Tuttle
,
M.L.
, and
Reynolds
,
R.L.
,
1993
,
The analysis of forms of sulfur in ancient sediments and sedimentary rocks: Comments and cautions
:
Science
 , v.
107
, p.
83
95
.
97.
Sandwell
,
D.T.
,
Müller
,
R.D.
,
Smith
,
W.H.F.
,
Garcia
,
E.
, and
Francis
,
R.
,
2014
,
New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure
:
Science
 , v.
346
, p.
65
67
.
98.
Savin
,
S.M.
, and
Lee
,
M.
,
1988
,
Isotopic studies of phyllosilicates
:
Reviews in Mineralogy
 , v.
19
, p.
189
223
.
99.
Scott
,
S.D.
, and
Barnes
,
H.L.
,
1971
,
Sphalerite geothermometry and geobarometry
:
Economic Geology
 , v.
66
, p.
653
669
.
100.
Shanks
,
W.C.
III
,
Böhlke
,
J.K.
, and
Seal
,
R.R.
II
,
1995
,
Stable isotopes in midocean ridge hydrothermal systems
:
Interactions between fluids, minerals, and organisms
 :
American Geophysical Union
Monograph 91, p.
194
221
.
101.
Shau
,
Y.H.
,
Peacor
,
D.R.
, and
Essene
,
E.J.
,
1990
,
Corrensite and mixed-layer chlorite-corrensite in metabasalt from northern Taiwan: TEM-AEM, EMPA, XRD, and optical studies
:
Contributions to Mineralogy and Petrology
 , v.
105
, p.
123
142
.
102.
Sheppard
,
S.M.F.
,
1986
,
Characterization and isotopic variations in natural waters
:
Reviews in Mineralogy
 , v.
16
, p.
165
183
.
103.
Sheppard
,
S.M.F.
, and
Gilg
,
H.A.
,
1996
,
Stable isotope geochemistry of clay minerals
:
Clay Minerals
 , v.
31
, p.
1
24
.
104.
Shikazono
,
N.
,
Holland
,
H.D.
, and
Quirk
,
R.F.
,
1983
,
Anhydrite in Kuroko deposits: Mode of occurrence and depositional mechanisms
:
Economic Geology Monograph 5
 , p.
329
344
.
105.
Shikazono
,
N.
,
Utada
,
M.
, and
Shimizu
,
M.
,
1995
,
Mineralogical and geochemical characteristics of hydrothermal alteration of basalt in the Kuroko mine area, Japan: Implications for the evolution of a back arc basin hydrothermal system
:
Applied Geochemistry
 , v.
10
, p.
621
642
.
106.
Shirozu
,
H.
,
1974
,
Clay minerals in altered wall rocks of the Kuroko-type deposits
:
Society of Mining Geology Special Issue
 , v.
6
, p.
303
311
.
107.
Sigurdsson
,
H.
,
2000
,
Volcanic episodes and rates of volcanism
, in
Sigurdsson
,
H.
,
Houghton
,
B.
,
Rymer
,
H.
,
Stix
,
J.
, and
McNutt
,
S.
, eds.,
Encyclopedia of volcanoes
 :
San Diego, California
,
Academic Press
, p.
271
279
.
108.
Skirrow
,
R.G.
, and
Franklin
,
J.M.
,
1994
,
Silicification and metal leaching in semiconformable alteration beneath the Chisel Lake massive sulfide deposit, Snow Lake, Manitoba
:
Economic Geology
 , v.
89
, p.
31
50
.
109.
Stoffers
,
P.
,
Worthington
,
T.J.
,
Schwarz-Schampera
,
U.
,
Hannington
,
M.D.
,
Massoth
,
G.J.
,
Hekinian
,
R.
,
Schmidt
,
M.
,
Lundsten
,
L.J.
,
Evans
,
L.J.
,
Vaiomo’unga
,
R.
, and
Kerby
,
T.
,
2006
,
Submarine volcanoes and high-temperature hydrothermal venting on the Tonga arc, southwest Pacific
:
Geology
 , v.
34
, p.
453
456
.
110.
Taylor
,
B.E.
,
1992
,
Degassing of H2O from rhyolite magma during eruption and shallow intrusion, and the isotopic composition of magmatic water in hydrothermal systems
, in
Hedenquist
,
J.W.
, ed.,
Magmatic contributions to hydrothermal systems: Geological Survey of Japan Reports
 , v.
279
, p.
190
194
.
111.
Thode
,
H.G.
, and
Monster
,
J.
,
1961
,
Sulphur isotope geochemistry
:
Geochimica et Cosmochimica Acta
 , v.
25
, p.
159
174
.
112.
Thompson
,
A.J.B.
,
Hauff
,
P.L.
, and
Robitaille
,
J.A.
,
1999
,
Alteration mapping in exploration: Application of short-wave infrared (SWIR) spectroscopy
:
Society of Economic Geologists Newsletter
 , no.
39
, p.
1
27
.
113.
Tostevin
,
R.
,
Turchyn
,
A.V.
,
Farquhar
,
J.
,
Johnston
,
D.T.
,
Eldridge
,
D.L.
,
Bishop
,
J.K.B.
, and
McIlvin
,
M.
,
2014
,
Multiple sulfur isotope constraints on the modern sulfur cycle
:
Earth and Planetary Science Letters
 , v.
396
, p.
14
21
.
114.
Ulrich
,
M.R.
, and
Bodnar
,
R.J.
,
1988
,
Systematics of stretching of fluid inclusions; II, Barite at 1 atm confining pressure
:
Economic Geology
 , v.
83
, p.
1037
1046
.
115.
Urabe
,
T.
,
Auzende
,
J.-M.
,
Deplus
,
C.
,
Eissen
,
J.-P.
,
Grimaud
,
D.
,
Huchon
,
P.
,
Ishibashi
,
J.
,
Joshima
,
M.
,
LaGabrielle
,
Y.
,
Mevel
,
C.
,
Naka
,
J.
,
Tanaka
,
T.
, and
Tanahashi
,
M.
,
1990
,
Hydrothermal deposits at the central spreading axis of the North Fiji basin (SW Pacific): Preliminary report of the Nautile dives of the STARMER cruise
:
Mining Geology
 , v.
40
, p.
117
124
(in Japanese with English abs.).
116.
Velde
,
B.
,
1980
,
Cell dimensions, polymorph type, and infrared spectra of synthetic white micas: The importance of ordering
:
American Mineralogist
 , v.
65
, p.
1277
1282
.
117.
Vergnoille
,
S.
, and
Mangan
,
M.T.
,
2000
,
Hawaiian and Strombolian eruptions
, in
Sigurdsson
,
H.
, ed.,
Encyclopedia of volcanoes
 :
San Diego
,
Academic Press
, p.
447
461
.
118.
Wienkenjohann
,
H.
,
2015
,
A sulfur isotope study of the Tinakula deposit, New Hebrides island arc, Solomon Islands
: Unpublished B.Sc. thesis,
Münster, Germany
,
Westfälische Wilhelms-Universität Münster
,
67
p.
119.
Wood
,
C.A.
,
1980
,
Morphometric evolution of cinder cones
:
Journal of Volcanology and Geothermal Research
 , v.
7
, p.
387
413
.
120.
Zhang
,
Y.
, and
Frantz
,
J.D.
,
1987
,
Determination of the homogenization temperatures and densities of supercritical fluids in the system NaCl-KCl-CaCl2-H2O using synthetic fluid inclusions
:
Chemical Geology
 , v.
64
, p.
335
350
.
121.
Zhu
,
W.
,
Tivey
,
M.K.
,
Gittings
,
H.
, and
Craddock
,
P.R.
,
2007
,
Permeability-porosity relationships in seafloor vent deposits: Dependence on pore evolution processes
:
Journal of Geophysical Research
 , v.
112
,
B05208
.

Melissa O. Anderson is a newly appointed assistant professor of economic geology at the University of Toronto, Canada. She obtained her B.Sc. (Hon.) degree at Brandon University, Canada, and an M.Sc. degree at the University of New Brunswick, Canada, under the supervision of Drs. Dave Lentz and Chris McFarlane. She completed her Ph.D. degree in 2018 at the University of Ottawa, Canada, under the supervision of Dr. Mark Hannington and spent time as a visiting researcher at GEOMAR Helmholtz Centre for Ocean Research in Kiel, Germany. Her research focuses on understanding the relationships between geodynamics and ore formation in modern seafloor environments.

Gold Open Access: this paper is published under the terms of the CC-BY 3.0 license.

Supplementary data