The northwestern caldera wall of Brothers volcano in the southern Kermadec arc features several clusters of hydrothermal venting in a large area that extends from near the caldera floor (~1700 mbsl) almost up to the crater rim (~1300 mbsl). Abundant black smoker-type hydrothermal chimneys and exposed stockwork mineralization in this area provide an excellent archive of hydrothermal processes that form seafloor massive sulfide deposits. Using sulfate precipitates from chimneys and stockwork recently recovered by remotely operated vehicles, we conducted fluid inclusion microthermometry and Sr isotope studies to determine the role of phase separation and mixing between vent fluid and seawater. The variability in the vast majority of fluid inclusion salinities (i.e., 0.1–5.25 wt.% NaCl eq.) and entrapment temperatures of up to 346°C are indicative of phase-separated hydrothermal fluids. Large salinity variations in samples with entrapment temperatures mostly below the boiling temperature for the sample’s depth show that the majority of fluids ascending below the NW Caldera are phase separating in the subsurface and cooling, prior to discharge. In several samples, entrapment temperatures of over 343°C suggest that phase-separating fluids have at least sporadically exited the seafloor at the NW Caldera site. Isobaric-isenthalpic mixing trends between coexisting phase-separated vapors and brines with seawater are consistent with phase-separated fluids at near-seafloor pressures of ~170 bar and suggest that the vast majority of the ascending fluids continue to phase separate to within tens to hundreds of meters below seafloor prior to mixing with seawater. A small subset of the most saline fluid inclusions (up to 18.6 wt.% NaCl eq.) is unlikely formed by near-seafloor phase separation and is considered to be produced either by supercritical phase separation or by the contribution of a magmatic brine from near the magmatic-hydrothermal interface. 87Sr/86Sr values of sulfate samples range from 0.7049 (i.e., near hydrothermal end-member) to 0.7090 (i.e., near seawater) and show that the crystals grew from vapor- and brine-derived fluids in a hydrothermally dominated mixing regime. Our work provides new insights into mineral growth conditions, mixing regimes, and in particular, the extent and character of subseafloor phase separation during the formation of hydrothermal vents and their underlying stockwork in seawater-dominated, arc-related hydrothermal systems.

1. Introduction

Microthermometric data can be used as geothermobarometers to determine the depth of crystal formation in magmatic-hydrothermal environments, where temperatures are known [1, 2], or to find the temperature of hydrothermal mineralization at known pressure conditions [3, 4]. In seafloor hydrothermal systems, fluid inclusion studies provide information about salinity and temperature of hydrothermal fluids, as well as likely hydrothermal end-member and mixing regimes during mineralization [512].

Microthermometric studies on fluid inclusions typically comprise petrographic descriptions and visual observations of phase transitions in the trapped fluids. Measurement of ice-melting temperatures () and homogenization temperatures () of two-phase liquid-vapor inclusions enables calculation of salinity () and entrapment temperature () of the inclusions [13, 14]. Under the assumption that a fluid inclusion has not changed in volume or composition since it formed in the host crystal, these measurements represent the conditions of the hydrothermal fluid during precipitation of the crystal [15]. In hydrothermal environments, fluid inclusions can record multiple fluctuations providing a record of conditions during the process of mineralization [5, 9, 16, 17].

In addition, values of 87Sr/86Sr for sulfate host minerals provide insight into the mixing regime in hydrothermal environments, as demonstrated for the TAG [18], PACMANUS [19], and the Endeavour Segment [8] hydrothermal systems. The fractionation factor for Sr isotopes during the precipitation of sulfate crystals is neglectable [20], and thus, 87Sr/86Sr values in sulfate can be used either as a tracer for Sr in the hydrothermal fluid or as a “fingerprint” of the seawater-fluid mixture the crystal has precipitated from. The combination of fluid inclusion with temperatures deduced from 87Sr/86Sr values for fluid-seawater mixing hence may provide evidence for conductive cooling (or conductive heating) and constrain thermobarometric fluid inclusion results [12].

We combined microthermometric data of fluid inclusions with 87Sr/86Sr data of sulfate hosts to provide detailed insights into the extent of phase-separation processes and mixing regimes in hydrothermal environments. The mineralization at the NW Caldera (and Upper Caldera) vent site can be regarded to as a modern analogue for volcanogenic massive sulfides. Here, extensive sulfate mineralization occurred simultaneously to sulfide deposition and access to chimneys and their underlying stockwork provide an excellent opportunity to use this approach and investigate phase separation and seawater entrainment related to the deposition of massive sulfides.

2. Geologic Setting

The submarine Brothers volcano is part of the southern Kermadec arc and is situated ~350 km NNE of the North Island of New Zealand (Figure 1). Brothers volcano is arguably one of the most well-studied submarine arc volcanoes in the world. Various studies have focused on its volcanic structure [21, 22] and magmatic evolution [2226].

The mostly dacitic volcano [2224, 26, 27] has a wide caldera that hosts numerous vent sites that have been delineated by shipborne and/or autonomous underwater vehicle plume and magnetic surveys [2832]. Active hydrothermal venting occurs at three vent sites on the caldera walls (i.e., Upper Caldera, NW Caldera, and West Caldera (not shown), Figure 2(a)) and at two vent sites atop the resurgent cones (i.e., Upper Cone and Lower Cone, Figure 2(a)) in the south-central part of the caldera [28, 29, 33].

The nature and extent of hydrothermal activity and its circulation [28, 31, 3440] and mineralization [7, 33, 36, 4143] have been addressed in various studies.

Recently, IODP expedition 376 drilled basement in five sites at Brothers volcano, and the recovered cores reveal new details on hydrothermal mineral zonation and the links between caldera collapse and the evolution of the hydrothermal system(s) [7, 44].

The Upper Cone and Lower Cone vent sites host relatively low-temperature magmatic volatile-dominated, white smoker-type vents, whereas the most vigorous venting at Brothers occurs at the NW and Upper Caldera sites in the form of numerous clusters of black smoker chimney vents expelling high-temperature (to ~320°C) fluids affected by deep-seated seawater-rock reactions [33, 37, 47, 48]. Here, the large NW Caldera site spans an area of at least , extending from 1700 meter below sea level (mbsl) near the caldera floor up to a prominent plateau at 1520 mbsl [31]. Above this plateau, the hydrothermal activity continues to the uppermost part of the caldera rim at 1300 mbsl (the Upper Caldera). The two vent fields are characterized by extensive sulfate and sulfide mineralization manifest as chimneys and stockwork zones [6, 33, 4143] and are coincident with an area depicting a negative magnetic anomaly [29, 34], projected to extend ~1000 meters below seafloor (mbsf).

The mineralization at the NW Caldera vent site is a consequence of vigorous hydrothermal upwelling and represents a striking modern example of fluid flow beneath volcanogenic massive sulfides. Access to active chimneys and their underlying stockwork is unique and makes the NW Caldera wall an exceptional place to study recent processes related to ore formation. The stockwork zone is dominated by clasts of argillic-propyllitic altered volcanics with disseminated pyrite and massive vein networks of either pyrite ± chalcopyrite or anhydrite + barite [47]. Hydrodynamic models suggest the vent field to be controlled by subsurface permeability structures [39, 40]. Investigations of vent fluid composition [33, 47, 48] show salinity variations that point to phase-separation processes at depth.

Fluid inclusions were studied at Brothers before. Measurements of fluid inclusions in subsurface core samples during IODP expedition 376 demonstrated two types of hydrothermal fluids circulating in the subsurface: subcritically phase-separated fluids with salinities of less than 5 wt.% NaCl eq. and hypersaline brines with salinities of ~32-45 wt.% NaCl eq. that formed as exsolution from a silicate melt, or condensed from a high-temperature single-phase fluid at the magmatic-hydrothermal interface [7]. Microthermometric investigations of two seafloor samples dredged from Brothers NW Caldera did not provide conclusive evidence for phase-separation processes. Measurement of ~85 fluid inclusions in two barite-rich samples showed narrow salinity () ranges of NaCl eq. and NaCl eq., with average homogenization temperatures () of and , respectively [6]. Furthermore, no physical evidence for boiling (e.g., liquid-rich inclusions in close proximity to vapor-rich inclusions) was observed in that microthermometric study. These authors used small deviations from seawater salinity and fluid inclusion gas analyses to suggest that phase separation may have occurred in the subseafloor. The importance of phase separation in the formation of mineralized chimneys and stockwork zones and its role in the formation of this volcanogenic massive sulfide deposit are further evaluated in the present study.

3. Methods

3.1. Sample Collection and Preparation

Research expeditions SO253 (R/V Sonne and ROV Quest, Dec 2016-Jan 2017) and TN350 (R/V Thompson and ROV Jason II, Mar 2018) set out to investigate hydrothermal activity at Brothers volcano. Both expeditions investigated four of the five active vent sites at Brothers volcano (i.e., NW Caldera, Upper Caldera, Upper Cone, and Lower Cone) using seafloor observations, heat flow measurements, and sampling of hydrothermal fluids and chimneys [49, 50].

The focus of this study is the NW Caldera site and the Upper Caldera site, which are characterized by the regular occurrence of barite and anhydrite in chimneys and stockwork mineralization. Thirteen samples from these sites were collected in order to characterize the nature of the hydrothermal fluids that gave rise to the massive sulfide mineralization (Table 1); twelve samples originate from the lower part of the caldera walls situated below a prominent plateau, or caldera rim, at a water depth of ~1520 mbsl (Figure 2(b); Table 1), and one sample originates from the walls above this plateau in an area referred to as the Upper Caldera. These samples include three distinct lithologies: chimneys, stockwork, and a propylitically altered volcanic wallrock adjacent to the stockwork.

Different structural and mineralogical types of chimneys can be distinguished. A Cu-rich pyrite + chalcopyrite + anhydrite/barite ± sphalerite ± Fe-oxyhydroxide type is most common. These Cu-rich chimneys show a distinct mineralogical zonation with large inner conduits covered by mm-cm thick massive layers of chalcopyrite. This zone rapidly grades into a highly permeable inner wall composed of pyrite + barite + sphalerite. A less permeable, several mm thick outer wall is composed of massive oxyhydroxide crusts. Zn-rich chimneys are mainly composed of pyrite + sphalerite + anhydrite/barite. In some instances, barite or anhydrite occur as a monomineralic layer within these samples. Chimney sample 064-10S represents a third class of Fe-rich chimneys with pyrite + pyrrhotite ± anhydrite ± sphalerite. In contrast to the Cu-rich type, Zn- and Fe-rich chimneys often lack an inner conduit and a strong mineralogical zonation in the chimney wall. Active chimneys of this type often feature fragile beehive tops of finely intergrown anhydrite and sphalerite, through which vent fluid issues diffusely.

The stockwork samples comprise propylitic- to argillic-altered volcanic rocks with disseminated sulfate and sulfide. The rocks are cut by massive veins of sulfate (anhydrite and/or barite) or sulfide (pyrite ± chalcopyrite ± sphalerite). The vein types occur in different areas along the NW Caldera wall, with sulfate veins more prevalent in the northern outcrop of stockwork mineralization (e.g., sample 064-9R). Here, the host rocks are altered to Fe-Mg-rich, mostly chlorite-dominated propylitic assemblages including chlorite + illite + smectite [47]. Sulfide veins are found in the southernmost outcrop (e.g., samples 085-6R, 085-8R, and 1039-17-R1) where the host rocks are altered to an argillic assemblage of quartz + smectite + pyrite + illite [47].

Crystals were extracted from each of these lithologies for detailed microthermometric analyses. In the chimneys, the crystals were extracted either from the central conduit that was lined with barite and/or anhydrite crystals, or from the central wall of the chimney, where mixing of the hot hydrothermal fluids from the inner conduit and ambient seawater has taken place. In the stockwork lithologies, barite and anhydrite crystals either originate from vein material or were extracted from the altered groundmass (in samples with sulfide veins). Finally, crystals from the propylitic-altered volcanic host rock originate from voids that are lined with larger (up to mm size) quartz crystals. Sample preparation for microthermometry consisted of separating transparent minerals from the rock/chimney specimen and picking of euhedral or carefully crushed crystals (with plane crystal or cleavage faces). The shortest diameter of these crystals was chosen to be smaller than a few hundreds of μm to assure good optical conditions and to prevent thermal gradients within the crystals during microthermometric work. Independent petrographic investigation in thin sections of the corresponding samples has shown that the chimneys typically contain anhydrite, as well as barite. Hence, each individual crystal used for the microthermometric study was investigated under a petrographic microscope, and anhydrite and barite were distinguished by their difference in birefringence (i.e., barite: 0.012; anhydrite: 0.044).

3.2. Fluid Inclusion Petrography and Microthermometry

Our microthermometric study comprised the following: (1) detailed petrographic description and photo documentation of the fluid inclusions, (2) cooling the inclusions to about -60°C with subsequent heating to detect the melting temperature (), and (3) heating to high temperatures to observe the homogenization temperature () of two-phase inclusions. Petrographic description and photo documentation were followed by the classification of the inclusions as primary, pseudosecondary, or secondary [52, 53]. Additionally, we classified the inclusions according to their phase relations at room temperature.

A Linkam microthermometry heating/freezing stage attached to a Zeiss Axioskop microscope was used to measure and of these fluid inclusions. A Linkam TP92 heat controller provided dynamic temperature control to heat the samples with a defined heating rate between 0.1 and 90°C min-1. An LNP2 flow regulator fed the stage with a constant stream of nitrogen gas from a Dewar of liquid nitrogen. The stage-temperature reading of the system was calibrated using the FLINC® synthetic fluid inclusions in quartz. The freezing point of pure water (0.0°C) and the critical homogenization temperature (374.1°C) were both used to calibrate the stage’s thermocouple sensor. During the study, the calibration procedure was repeated 14 times; the freezing point was determined to be , whereas the critical temperature was determined to be . This shows that the stage provided accurate temperatures between the measured values of the synthetic fluid inclusions (; ). In each sample, several inclusions were cooled to -100°C to detect the presence of mixtures others than H2O-NaCl. Since in all but two inclusions no phase transitions except for the ones expected in the H2O-NaCl system were observed, the H2O-NaCl system was used to interpret phase transitions. The microthermometric technique to measure and and to calculate salinity () and entrapment temperature () using the H2O-NaCl binary system is described by Bodnar and Vityk [14]. Salinity was calculated directly from , and was calculated from and using a pressure correction and following an isochoric path from the pressure of to the pressure of .

Generally, reequilibration can be a problem in fluid inclusion studies [15]. Considering that the minerals anhydrite and barite have low hardness and are prone to stretching and reequilibration during the fluid inclusion measurement procedure [15], the technique described by Vanko and Bach [54] was applied. That is, their heating and freezing experiments on anhydrite-hosted fluid inclusions showed that during microthermometric measurements, hydration of anhydrite hosts may form gypsum at the inclusion’s wall under lower temperature conditions (i.e., at minus several tens °C). In subsequent measurements, the salinity of the fluid inclusion may therefore be increased. This process seems to occur reversibly, as heating the sample to 130°C converts the newly formed gypsum back to anhydrite. Furthermore, Vanko and Bach [54] found that overheating the inclusions by only ~10°C above their commonly leads to stretching or reequilibration and thus significantly changes measured of previously overheated inclusions in anhydrite. Thus, using the procedure proposed by Vanko and Bach [54], inclusions were not measured one after another, but instead, all were measured in freezing runs first. This was an attempt to retain original inclusion shapes (volume) and to prevent stretching. Thereafter, the temperature of the sample was subsequently increased for the measurements with all inclusions observed at 5°C intervals. Following this procedure, no inclusion was significantly overheated before measuring .

3.3. Analysis of 87Sr/86Sr

Ratios of 87Sr/86Sr for five hand-picked mineral separates of barite and anhydrite and six bulk sample powders were analyzed by thermal ionization mass spectrometry (TIMS). A Thermo Scientific Triton Plus mass spectrometer was used at the Isotope Geochemistry Laboratory, Center for Marine Environmental Science MARUM, University of Bremen. In total, eleven analyses (eight stockwork, two chimneys, and one propyliticly altered rock) were carried out to determine 87Sr/86Sr values on six individual sulfate samples.

The samples were digested in concentrated subboiled HNO3 at 130°C. Strontium was isolated from the matrix elements, including Ba, using SrSpec resin by a method adapted from Deniel and Pin [55]. Samples were analyzed with a Ta-emitter on rhenium filaments using the static acquisition method. Two analyses of the NIST SRM 987 reference material were conducted along with the samples, which yielded 87Sr/86Sr values of and that plot well within the long-term reproducibility value of the laboratory (; 2σ of ). This in turn compares well to average published data for NIST SRM 987 (; 2σ, , are discarded; GeoReM database [56], request September 2017).

3.4. Isenthalpic-Isobaric Mixing Models

Fluid inclusion data and 87Sr/86Sr ratios were used to develop mixing models between seawater and hydrothermal fluids. We used --dependent brine-vapor pairs formed by phase separation to account for the variation in the fluid inclusion data. Potential end-member pairs matching the fluid inclusion data were determined by the phase relations in the H2O-NaCl system as described by Driesner and Heinrich [57]. Isobaric liquid-vapor (LV) curves were calculated via the program “SOWAT” from which pressure- and temperature-dependent coexisting vapor and brine pairs can be derived.

For simulation of isenthalpic mixing between ambient seawater and phase-separated vapors and brines, we calculated the heat capacity as a function of -- following the procedure of Driesner [58]. The thermodynamic properties of pure water are given by Haar et al. [59] and were used to develop a correlation method to calculate the thermodynamic properties of binary H2O-NaCl solutions. From this data, we calculated the compositional dependency of the heat capacity of salty solutions (see Driesner [58] and references therein). For individual mixing paths of our different end-member brines and vapors with seawater, we simultaneously calculated heat capacities and temperatures of those mixtures.

4. Results

4.1. Fluid Inclusion Petrography and Classification

We studied a total of 471 fluid inclusions hosted by barite (, 59%), anhydrite (, 37%), and quartz (, 4%). The inclusions were measured in chimneys (, 50%), stockwork material (, 46%), and propylitically altered wallrock (, 4%). We separated both barite and anhydrite from chimneys and stockwork zones, while quartz represented vug infills in the one altered rock specimen. All inclusions were classified following the petrographic classification scheme outlined by Roedder [53]. He described primary inclusions to occur as single inclusions, as small groups of inclusions, or in trails that follow crystal growth faces. Secondary inclusions occur as trails parallel to cleavage faces or with no visible relation to crystal growth faces and do intersect crystal boundaries. Pseudosecondary inclusions occur like secondary inclusions, but trails do not intersect crystal boundaries. In our investigation inclusions are classified into primary (; 85%) and secondary (; 15%) inclusions. No pseudosecondary inclusions are identified. Inclusions are either euhedral, with a negative crystal shape (; 59%), round (; 25%), or irregular (; 16%) in shape. Inclusions where the largest diameter is equal to, or larger by three times the smallest diameter, are classified as elongated (; 28%). All shapes occur for primary, as well as secondary, inclusions, but secondary inclusions tend to be rather round. Fluid inclusion size is calculated using the smallest and the largest diameter to calculate the area of an ellipse. Inclusions (; 28%) with an area smaller than 5 μm2 are considered small. If the area ranges between 5 and 50 μm2, the inclusions are considered medium sized (; 49%), and the remaining inclusions (; 23%) with a μm2 are considered large.

The volume proportions of the vapor and liquid phase were estimated from the area expression of the phases in microscopic view. The area occupied by the vapor bubble was related to the total inclusion area. The majority of inclusions are liquid rich at room temperature, with most inclusions (; 46%) having a vapor fraction of 20%. In total, in 402 inclusions (87%), the vapor fraction occupies <50%, while in a smaller number of inclusions (), the vapor fraction occupies ≥50%. The remaining inclusions decrepitated before the liquid-vapor ratio was detected or were found as single-phase (liquid) inclusions. According to phase relations at room temperature, the fluid inclusions are further divided into three types. Type I inclusions are single phase and were not further investigated. Abundant () Type II inclusions are two phased and contain a vapor bubble and a liquid phase. Type II inclusions are subdivided into liquid-rich Type IIa (, 89%) and vapor-rich Type IIb (), with a division set at a liquid-to-vapor ratio of 0.5. Type IIa and Type IIb inclusions both co-occur in single trails. Rare Type III inclusions (, 2%) are three-phased and contain a liquid, a vapor bubble, and a transparent solid phase. Figure 3 provides a representative overview of the petrographic relationship of the NW Caldera fluid inclusions.

4.2. Microthermometry

Microthermometric measurements are shown in Table 2 and Figure 4, including of 429 and of 346 inclusions (Figures 4(a) and 4(b)) out of 471 investigated inclusions. Notably, the of two inclusions is +1.8°C and +3.3°C, respectively, which is higher than the for pure H2O, and hence higher than possible for H2O-NaCl binary solutions. The melting temperature in these inclusions can only be explained with the presence of solutes that shift towards higher temperatures in these inclusions. Two inclusions have of -14.0°C and -15.0°C, which are distinctly lower than in all other inclusions and result in significantly higher salinities for these inclusions. For the sake of clarity and to keep statistical integrity in our assessment, these four inclusions are excluded in the following figures and tables and will be discussed separately. For 325 inclusions, both and were measured (Figure 4(c), Table 2). For these inclusions, the corresponding and are calculated. As the pressure during homogenization of the fluid inclusions is close to the pressure of entrapment (i.e., the pressure at the seafloor), the path along the isochore is short, and hence, the correction for temperature is similarly small. Overall, the temperature correction (i.e., the difference between and ) is small, ranging between 2°C and 10°C (average 9°C).

The inability of being able to detect ice melting in 42 inclusions can be mostly attributed to optical limitations. For example, some Type II inclusions were clearly visible, but ice melting was still not observed. As these inclusions are vapor rich, they contain only a small liquid component, which made it hard to observe the process of ice melting. Some of the 134 inclusions, in which could not be determined, were simply missed during heating runs and thus have been excluded to avoid measurements of potentially stretched inclusions. Some inclusions did stretch, or decrepitate, and changed from two phases to a single phase during the microthermometric measurements, even at relatively low temperatures (<<100°C). Indeed, most of the inclusions could not be measured as a consequence of (audible) decrepitation events, which occurred in barite as well as anhydrite during measurement of . Such decrepitation events occurred regularly at °C. Forty-five inclusions were not measured due to decrepitation events. For all these inclusions, is >270°C; for 39 of them, is >300°C; for 9, is >340°C; for 5 inclusions, is >400°C; and for 2 inclusions, is >500°C. Decrepitation often occurred while heating Type IIb inclusions, which are therefore underrepresented in the statistical treatment of the fluid inclusion data.

In summary, among all inclusions that could not be homogenized, two-third is of Type IIa and one-third is of Type IIb. The difference to the overall abundance of inclusion types, in which Type IIb comprised only 9%, suggests that these inclusions commonly homogenize at higher temperatures and are prone to stretching and/or are prone to decrepitation. For these inclusions, it was noted that the vapor phase grows during heating, instead of shrinking like Type IIa. Thus, we suggest that a number of these inclusions homogenize into the vapor phase and represent low-density inclusions.

In Type III inclusions, no phase transitions were detected. For two of these inclusions, temperatures of up to 400°C during the microthermometric measurements did not lead to dissolution of the daughter mineral. The fact that no visible phase transition occurred at low temperatures suggests that these inclusions do not contain saline solutions with salinities below the eutectic composition comparable to the measured Type II inclusions. If the transparent daughter mineral is a sulfate crystal, these inclusions should show ice melting comparable to Type II inclusions. In contrast, if the Type III inclusions contain H2O-NaCl solutions with a salinity above the eutectic composition, no ice melting is to be expected at low (below 0°C), as we observed. For high-salinity inclusions, the dissolution of halite is the first phase transition to be observed at high temperature, followed by complete homogenization at even higher temperatures. Assuming that the inclusions indeed represent highly saline inclusions with halite daughter minerals and their dissolution occurs above 400°C, this means that the incorporated fluids contain NaCl equivalent.

4.3. Salinities and Entrapment Temperatures

A total of 325 inclusions have both and measurements, for which the corresponding and are calculated (Table 2). In addition, values are calculated for 427 inclusions based on thermodynamic relationships in the binary H2O-NaCl system. Two inclusions with of 17.7 and 18.6 wt.% NaCl eq., respectively, are discussed separately because these two inclusions would impede data visualization and significantly bias statistical parameters.

Calculated salinities () in the remaining 425 measured fluid inclusions range from 0.1 to almost 9 wt.% NaCl eq., with a mean value of NaCl eq. (cf. seawater is 3.2 wt.% NaCl equivalent). The data show a multimodal distribution, with two pronounced modes. The first and most prominent group includes 225 inclusions and occurs between 3.0 and 4.0 wt.% NaCl eq. (Figure 5(a)), coincident with salinities in fluid inclusions measured by de Ronde et al. [6]. The second group includes 45 fluid inclusions in a narrower range, between 1.5 and 1.8 wt.% NaCl equivalent. Finally, a number of fluid inclusions show either very low values of of 0–1 wt.% NaCl eq. () or much higher values of >4 wt.% NaCl eq. (). These fluid inclusions, with noticeably higher and lower salinities when compared to seawater, have not been previously reported in seafloor samples collected from the NW Caldera site. The of Type IIb inclusions essentially show the same data distribution as Type IIa inclusions (Figure 5(b)), despite the fact that the relative abundance of inclusions with NaCl eq. is larger than that for Type IIa inclusions. In terms of seawater salinity (), the inclusions are divided into three groups: (1) those with salinities significantly below that of seawater (<3.1 wt.% NaCl eq.), comprising of 33% of the inclusions; (2) salinities of seawater composition (3.1–3.3 wt.% NaCl eq.), comprising of 16% of the inclusions; and (3) salinities significantly higher than seawater (>3.3 wt.% NaCl eq.), comprising 51% of the inclusions. We compared the range in salinities of the fluid inclusions with the salinity variations of the vent fluids from the NW Caldera site (Figure 5(c)). The majority of vent fluid salinities (45%) are seawater like (3.1–3.3 wt.% NaCl eq.) or slightly above (>3.3 wt.% NaCl eq. and up to 4.5 wt.% NaCl eq.) seawater salinity (51%). Only two out of 49 vent fluid salinities (or 4% of the data) lie significantly below the salinity of seawater. While overall the bulk of the vent fluid salinities are similar to those found associated with the fluid inclusions, the most extreme values of NaCl eq. and NaCl eq. that represent a considerable fraction (20% of all measured ) of the fluid inclusion data are not seen in vent fluid compositions.

for the entire set of inclusions fall in the range between 135°C and 347°C, with a mean value of . This value is remarkably close to fluid inclusion results given in the study of de Ronde et al. [6]. The distribution of , like , also shows a multimodal distribution with two pronounced modes (Figure 5(d)). One mode is situated at around 250°C and coincides with the average from the work of de Ronde et al. [6]. Another one consists of inclusions with of 280–300°C, which overlaps with the interval of highest given by de Ronde et al. [6]. A third mode at of 330°C shows temperatures in excess of the aforementioned study, with the exception of one inclusion with of 322°C in de Ronde et al. [6]. The few determinations made for Type IIb inclusions (Figure 5(e)) show that the relative abundance of inclusions formed at temperatures in excess of 325°C (seven out of 18, or > 33% of all Type IIb determinations) is higher than that for Type IIa inclusions, coincident with a higher proportion of relatively high inclusions. Measured vent fluid temperatures (Figure 5(f)) agree with values for the fluid inclusions. Temperatures of a small number of vent fluids are below 250°C, but most samples sit in two groups of higher fluid temperatures. The first group comprises of about 30% of all data with fluid temperatures between 250 and 275°C and can be matched to the mode at ~250°C. The second group (43%) is given by fluid temperatures of between 290 and 310°C and can be matched to the inclusion mode at 300°C. One fluid sample with a temperature of 320°C can be matched to the inclusion mode of higher at 330°C. Higher vent fluid temperatures of have not been reported although such high fluid inclusion represents almost 10% of all measured .

4.4. Temperature and Salinity Variations with Lithology and Depth

If we consider the different lithologies (chimneys, stockwork, and propylitic wallrock), we find that each lithology and each mineral type records the entire spectrum of fluid inclusion and (Figure 6), especially inclusions with salinities similar to, or greater than, seawater that mark the most pronounced mode of the entire population are found in stockwork, chimneys, and propylitic rock. Still, there are apparent relationships between / and lithology, as well as mineral type. These are the lowest (<0.7 wt.% NaCl eq.) and of 280–320°C which are measured almost exclusively in chimney samples (Figures 6(a) and 6(b)). In terms of host minerals in chimneys, a relationship between inclusion composition and the host mineral is found. Fluid inclusions in anhydrite record proportionally lower salinity fluids and higher temperatures when compared to barite-hosted fluid inclusions (Figures 6(a) and 6(b)).

In contrast to chimneys, low-salinity inclusions in stockwork samples are mostly confined to higher , in the range of 0.7–3.0 wt.% NaCl eq., and are predominantly found at lower of ~250°C (Figures 6(c) and 6(d)). In terms of mineral types, there is no apparent relationship between / and the host mineral for stockwork mineralization (Figures 6(c) and 6(d)). Anhydrite and barite inclusions cover virtually the entire range of /. However, inclusions with the lowest (<1 wt.% NaCl eq.) and highest (>4 wt.% NaCl eq.) salinities are here exclusively found in barite crystals.

Inclusions from quartz crystals in the propylitically altered wallrock (Table 2, Figures 6(e) and 6(f)) show a narrower but still broad range of (1.1–5.0 wt.% NaCl eq.) and a narrow range of higher (306–346°C). Despite the small sample size of only 17 measured inclusions, these fluid inclusions make up almost half of the inclusions in the mode of highest (~330°C) given for the overall dataset.

The variation of and as a function of water depth is shown in Figure 7. The fluid inclusion data is plotted on a projected depth profile of the inner caldera at Brothers and is compared to seawater salinity and the seawater boiling curve. This plot shows how the modes of and in the overall dataset are related to individual samples (cf. Table 2).

The mean salinities for single samples range from (sample 064-10S) to NaCl eq. (sample 1038-R8), with medians of 0.9 and 3.4 wt.% NaCl eq., respectively. Mean values for single samples range from (i.e., sample 085-8R) to (i.e., sample 081-6R). Median values for these samples are 198°C and 329°C, respectively.

Mean values for close to seawater are characteristic of 11 of 13 samples (Figure 7(a)). Two samples give mean NaCl equivalent. Although mean are similar to seawater on average, individual inclusions in 10 of 13 samples contain a wide range of (>>3 wt.% NaCl eq.) spanning almost the entire range of the dataset. Four samples contain inclusions with very low (<0.5 wt.% NaCl eq.), and in another four samples, inclusions with the highest NaCl eq. are found. By contrast, three samples are different with a small range of in individual inclusions (1039-CH2, 081-1S, and 1041-CH1 at 1611, 1619, and 1622 mbsl, respectively). In these samples, the mean of all inclusions is close to seawater and the range of (<<2 wt.% NaCl eq.) is comparably small, which is expressed by a small standard deviation (<0.6 wt.% NaCl eq.). The distribution in these three samples is most comparable to the inclusion data of de Ronde et al. [6].

Mean in individual samples cluster into three groups (Figure 7(b)). Eight of 13 samples with a mean of 200–250°C are situated deeper than 1600 mbsl. Four samples with a distinctly higher mean of 280–300°C are situated shallower than 1600 mbsl. The mean value of 330°C for the propylitically altered wallrock is significantly higher than all other mean temperatures. In all fluid inclusion samples, the mean is lower than the depth-to-boiling-point temperature of seawater for any given depth. However, a subset of samples shows maximum very close to (i.e., within 5°C) the boiling point curve for seawater. Such high occur in four of 13 samples investigated (i.e., samples 061-17S, 064-10S, 064-9R, and 081-6R; cf. Table 2).

4.5. 87Sr/86Sr Composition of Mineral Separates and Bulk Samples

Values for 87Sr/86Sr for eleven mineral separates and bulk rock samples range from near hydrothermal end-member values established for NW Caldera vent site hydrothermal fluids (, , composition extrapolated to ) to near seawater values (0.70918; [47]) (Table 3; Figure 8). The least radiogenic sample is an anhydrite separate from sample 085-6R1 with an 87Sr/86Sr value of 0.704852. By contrast, the most radiogenic sample is a barite separate from sample 064-9R with an 87Sr/86Sr value of 0.709040. The 87Sr/86Sr value for two samples (064-9R and 085-6R1) has been measured multiple times in bulk rock samples and inclusion separates. The measurements show that the isotopic composition for these samples is highly variable. For example, in sample 064-9R, the 87Sr/86Sr values vary from 0.705260 to 0.708572 in four subsamples of one rock specimen. Similarly, for sample 085-6R1, values of four subsamples of one specimen range from 0.704852 to 0.706894. This shows the inhomogeneous nature of the sulfate crystals originating in this environment. Values of 87Sr/86Sr have been measured either from single crystals from inclusion separates or from small portions of crushed bulk rock adjacent to the inclusion samples. The measurement of six crushed bulk samples shows a narrower range of isotope compositions from 0.706130 to 0.708572 than the larger range of 0.704852 to 0.709040 for individual crystals (cf. Table 3 and Figure 8).

4.6. Isenthalpic-Isobaric Mixing Models

The fluid inclusion and 87Sr/86Sr data from this study are used to develop mixing models between seawater and hydrothermal fluids. The fluid inclusion data given in this work provide direct evidence for phase separation having occurred at the NW Caldera site, as other recent studies do [7, 33, 48].

In simulations of isenthalpic mixing between ambient seawater, phase-separated vapors, and brines, we calculate the heat capacity as a function of -- following the procedure of [58]. Isenthalpic-isobaric mixing models are then employed to examine possible mixing trends for generating the observed variability in and in fluid inclusions, as well as Sr isotope variations in the host phases. In these models, we determine the positions of phase-separated end-member fluids and their corresponding mixing trends in - space, where is temperature and is either salinity or Sr isotope composition. The computations are conducted for subcritical conditions, based on evidence from the fluid inclusion data that demonstrate near the boiling temperature at seafloor pressures and close to 0 wt.% NaCl equivalent. Further, we compute simulations for supercritical conditions to account for high-salinity inclusions (>>5 wt.% NaCl eq.). Isobaric - sections of phase relations in -- space are computed following Driesner and Heinrich [57]. Our model results predict the consequences of mixing between subcritically and supercritically phase-separated fluids and seawater in terms of , , and 87Sr/86Sr variations along isenthalpic-isobaric mixing paths. Plausible brine-vapor pairs are tabulated in Table 4.

Hydrothermal end-member fluids produced by near-seafloor phase- separation at 170 bar (~1700 mbsl or ~0–300 mbsf; dependent on the location on the caldera wall) and by subseafloor phase separation at supercritical pressures of 340 bar (~3000 mbsl or ~1300–1600 mbsf assuming a fully hydrostatic pressure regime) are considered. We compute mixing paths between these phase-separated end-member fluids with ambient seawater (5°C).

Figure 9 presents our best-fit isenthalpic mixtures at a pressure of 170 bar with an end-member temperature of 354°C. The vapor-brine phases are formed when a fluid with seawater salinity fractionates into 31 wt.% vapor ( NaCl eq.) and 69 wt.% brine ( NaCl equivalent). Phase separation of an ascending fluid terminates under equilibrium conditions at 170 bar. The vapor and brine fractions of the chosen end-member broadly resemble the overall occurrence of fluid inclusion salinities with approximately 33% of all inclusions showing lower than seawater, and 51% of inclusions showing higher than seawater (the remaining inclusions are indistinguishable from seawater salinity). Phase segregation of this brine/vapor pair can produce any composition between those two end-members, and subsequent mixing with seawater results in the - space highlighted in Figure 9(a).

The model scenario may explain the vast majority of fluid inclusion data except for two groups of inclusions. The first group involves those inclusions with lowest in combination with low . These inclusions can be explained only when significant conductive cooling (of ~50–100°C) is involved in addition to the hydrothermal fluid/seawater mixing. The second group comprises inclusions with relatively high salinities ( NaCl eq.) in combination with low (<270°C). A number of these inclusions (i.e., those with ) can also be explained by conductive cooling, but the model cannot account for the even higher salinities with the chosen end-member brine/vapor pair. A lower was measured in a number of inclusions where no could be determined. These inclusions are not plotted in Figure 9(a) but provide direct evidence for a contribution of fluids with significantly higher salinities than the range explained by this model. In total, the of nine inclusions indicate higher than 6 wt.% NaCl eq. and up to 18.6 wt.% NaCl equivalent.

Figures 9(b) and 9(c) show how formation temperatures and salinities covary with the 87Sr/86Sr ratios for the mixed fluids. The mixing models can reasonably account for the variation of inclusion with respect to Sr isotope values (Figure 9(b)) and explain the covariation of salinity and 87Sr/86Sr (Figure 9(c)). Extremely low salinities in vapor end-members are expected to be accompanied by extremely low Sr concentrations (as the partitioning behavior of Sr is assumed to be equal; cf. Table 4). Vapor end-member hydrothermal fluids that mix with seawater will therefore exhibit seawater-like isotopic signatures, even with small seawater fractions in the mix.

Mixing models involving phase separation at different pressures are applied to account for the fluid inclusion data. Applying higher pressures in the phase-separation process results in both the temperature of the brine-vapor pair, as well as the salinity of the vapor end-member increasing. These models can therefore not account for the observed variability found in the fluid inclusion data, except if extreme conductive cooling scenarios are invoked.

Moreover, the occurrence of high fluid inclusions is not likely explained by shallow (near-seafloor) phase separation, because formation of a brine with much higher salinities (i.e., >10 wt.% NaCl) would require fractionation of large amounts of vapor from the original fluid and would lead to a vapor/brine fractionation where the vapor fraction dominates.

The high-salinity inclusions, however, can be explained by scenarios involving supercritical phase separation. We calculate an isobaric-isenthalpic mixing line for a vapor/brine pair of supercritically phase-separated fluids at a pressure of 340 bar (Table 4). The high-salinity inclusions plot well along a mixing trend between seawater and a phase-separated brine at such a pressure (Figure 10). The scarcity of data due to the lack of measurements for these high-salinity inclusions and the broad and indefinite freedom to choose pressure, temperature, and fractionation between a supercritical vapor and brine do not allow us to constrain the depth or nature of this possible supercritical phase-separation process any further.

5. Discussion

5.1. Evidence for Phase Separation from Microthermometry

The vast majority (98%) of fluid inclusions associated with mineralization at the Brothers NW Caldera (including the Upper Caldera) hydrothermal vent field show between 0.1 and 5.25 wt.% NaCl equivalent. Most of these inclusions are two phased, liquid-rich at room temperature, and are classified as Type IIa, while a significant number of these inclusions (48%) fall in a range between 3.3 and 5.25 wt.% NaCl equivalent. This group of inclusions coincides with previously measured vent fluid compositions [33, 47, 48] and previous fluid inclusion studies of seafloor samples by de Ronde et al. [6]. The latter reported a mean of NaCl eq. and a mean of , which are almost identical to the results given in this study (i.e., of NaCl eq.; of ; cf. Table 3). However, the ranges in both and of fluid inclusions indicated by de Ronde et al. [6] are much narrower than that in our study. The authors found that <15% of their inclusions significantly deviated from seawater and suggest that their data shows a slight increase of with increasing . Based on their microthermometric results from two samples and from gas analyses of fluid inclusions in those same samples, de Ronde et al. [6] suggest that the fluid inclusion data indicate subseafloor phase separation of the hydrothermal fluid prior to mixing with seawater and subsequent venting. Our microthermometric data, representing a much larger sample size, are consistent with the results of de Ronde et al. [6] and reinforce the significant role phase separation plays in the hydrothermal vent field of the NW and Upper Caldera, which accounts for the vast majority of fluid inclusions.

In our microthermometric dataset, about 50% of all inclusions have trapped hydrothermal fluids with greater than seawater. Similarly, about 33% of the measured inclusions have trapped hydrothermal fluids with less than seawater. Only about 16% of all inclusions can be regarded to have seawater salinity. By assuming that this distribution is representative for the hydrothermal regime in a temporal and spatial fashion, we conclude that the hydrothermal fluids that passed through the NW Caldera site are dominantly brine-derived solutions, with a significant fraction comprising vapor-derived solutions, and only a subordinate fraction with seawater-derived fluids that may not have undergone any phase separation.

A significant number of inclusions (14%) had NaCl eq., with as low as 0.1 wt.% NaCl equivalent. The occurrence of these inclusions provides direct evidence for shallow phase separation in depths of tens to few hundreds meter below seafloor and the circulation of vapor-rich hydrothermal fluids that are barely diluted by seawater. The occurrence of such low-salinity vapor-derived fluids is reported for the first time in seafloor samples from the NW Caldera vent field.

In a small subset (2%) of all inclusions,  wt.% NaCl eq. and up to 18.6 wt.% NaCl eq. were observed. These high-salinity inclusions are abundant in Type IIb and are similarly reported for the first time in seafloor samples from the NW Caldera vent field. Highly saline brines (32–45 wt.% NaCl eq.) in fluid inclusions have recently been reported in a microthermometric study of drillcore samples from the NW Caldera and the Cone Sites at Brothers [7]. Additional evidence for the occurrence of high-salinity fluids in our seafloor samples eventually was found in some rare fluid inclusions with transparent daughter crystals. The three-phased inclusions did not show any phase transitions up to 400°C. The tendency of decrepitation events during the microthermometric investigations has prevented us from heating to higher temperatures and evaluating the phase relations at such high temperatures. If the solid phase in these inclusions represents halite, these inclusions must contain fluids with NaCl eq., as dissolution of halite would occur below temperatures of 400°C for inclusions with lower salinities. We assume that the measured salinities () and the three-phased inclusions provide a direct link between the saline brines in the drillcore samples and attest the sporadic occurrence of these brines in the shallow part of the hydrothermal system and at the seafloor.

The entrapment temperatures in this study range between 135 and 347°C with three pronounced modes found regarding fluid inclusion ; the first and second modes fall between 230–270°C and 290–310°C, respectively. These two modes are in accordance with vent fluid temperatures in the NW Caldera site [33, 47, 48] and the microthermometric fluid inclusion study of de Ronde et al. [6]. In addition, we found fluid inclusions in a less pronounced (10% of all inclusions) third mode, with temperatures of 320–346°C that occur in all mineralization types and host mineral phases, and which have higher temperatures than measured for vent fluids in the NW Caldera. Maximum in several samples plot very close to the depth-to-boiling-point curve of seawater for the corresponding seafloor pressure. In addition to high and extreme variations, we find petrographic evidence for phase separation. Inclusions with variable liquid-vapor ratios occur in single trails of fluid inclusions with similar . Such inclusions are consistent with entrapment from phase-separated fluids [4].

The proposed occurrence of subcritical phase separation at the NW Caldera vent field is not unique and has been mentioned by several authors [6, 7, 31, 33, 47, 48]. Phase separation was predicted by numerical simulations of fluid and heat transport for Brothers volcano by Gruen et al. [39] and Gruen et al. [40]. Though, our study presents a wider range of fluid salinities and a larger population of non-seawater salinities compared to previous investigations of vent fluids and fluid inclusions in seafloor samples. This indicates that phase separation plays a larger role in controlling fluid compositions and modulating metal transport to the seafloor than previously appreciated.

5.2. Extent and Depth of Phase Separation in the NW Caldera Wall Vent Field

Analyzing the fluid inclusions in individual samples shows that maximum in four of 13 samples are very close to the depth-to-boiling-point temperature at the corresponding seafloor pressure. This observation suggests that in the recent past, actively phase-separating fluids have sporadically reached the seafloor without much cooling, or mixing, in the subseafloor, and that the minerals have grown from phase-separating fluids. This process has not been observed by direct vent fluid temperature measurements during the two expeditions to the NW Caldera site in 2016/2017 and 2018, where samples for this study were recovered. All measured vent fluid temperatures, both previously published and those of the two recent cruises, were notably (by ≥15–20°C) below the boiling point of seawater at the corresponding seafloor pressures [33, 47, 48].

Fluid salinities are variable in ten out of 13 samples, ranging from ≤50% to ≥150% seawater salinity, and are consistent with phase-separation processes having occurred prior to sulfate mineral growth. Salinity is a good marker of phase separation as fluid inclusions with a large range in and correspondingly below the boiling temperature represent mixing of hot, upwelling fluids and entrained ambient seawater prior to crystal growth. Such mixing processes cease phase separation in the subsurface due to efficient cooling but preserve variations in . This is consistent with the considerable variability of vent fluid salinities in combination with vent fluid temperatures explicitly lower than the boiling temperature. We therefore suggest that phase separation is a ubiquitous feature in the NW Caldera hydrothermal system, but the amount of cooling and mixing with cold seawater en route to the seafloor is variable in space and time.

An important observation made in this study is that fluid inclusions in the majority of samples and even within single crystals have trapped almost the entire spectrum of (i.e., less than and greater than seawater) observed at the NW Caldera vent site as a whole. This suggests that not only the process of fluid-seawater mixing but also the process of phase separation itself is indeed highly variable in space and time. Time-series analyses of vent fluid compositions at MOR hydrothermal systems showed a temporal evolution of phase-separation processes in vent sites. For example, vent sites at the Juan de Fuca ridge and EPR 9°50 N were visited following magmatic events and were found to be predominantly discharging vapor phase-derived fluids [6062]. Fluids at the same sites then transitioned into seawater- or brine-derived solutions after some months to years after these magmatic events had occurred.

Generally, for the NW Caldera site, radiometric dating of barite mineralization has shown that mineralization has occurred over an extended period of time, with 226Ra/Ba ages for (sulfate-rich) mineralization of ~1200 years and evidence for even older hydrothermal activity to exceed 15,000 years [31, 63, 64]. For active chimneys, 228Th/228Ra and/or 228Ra/226Ra ages have shown that recent mineralization has predominantly occurred over the last four years with some chimneys dating back to 30–40 years, at time of sampling [31, 33]. Since barite and anhydrite analyzed in this study are often intergrown and have probably precipitated coevally, we tentatively transfer the radiometric results to this study. A timespan of several years to tens of years is sufficient for several magmatic events, including intrusions of magmatic dikes in the deep subsurface, to affect the hydrothermal system by providing additional heat and lead to an alternating entrapment of brine- and vapor-derived fluids in our sulfate-hosted fluid inclusions.

Detailed dating has shown that typically the age of mineralization along a transect of active chimneys progressively increases from the outer wall to the inner wall, if not younger chimneys have grown adjacent to older ones [33]. If the phase-separation process was controlled by magmatic events and a progression from vapor- to brine-dominated fluids at the NW Caldera site occurred, we expect that certain salinity ranges are restricted to some samples depending on the timing of growth of the individual host crystal. However, the salinities are widely distributed throughout most of the discrete samples. This requires a rather ephemeral feature controlling the occurrence of vapors and brines. We suggest that the key in this variability is the process of phase-segregation, as proposed by Butterfield et al. [65]. Phase segregation may be an ephemeral feature in the permeability-controlled regime of the NW Caldera [39, 40]. During rock alteration and mineral precipitation in the subseafloor, fluid pathways may regularly open or close with this process controlling phase segregation, leading to alternating occurrences of vapor and brine over periods of days to months, instead of months to years for any chimney. The variability of our fluid inclusion data may be a product of both a temporal evolution of the site as a whole and a dynamic permeability regime in the subseafloor. Given the wide range of occurring in the fluid inclusions of most samples, we suggest the latter to be the dominant process.

5.3. Evidence from Isenthalpic-Isobaric Mixing Models

Isobaric-isenthalpic models of mixing between shallow (170 bar) phase-separated vapor-brine end-member pairs and ambient seawater can explain virtually the whole spectrum (92%) of our fluid inclusion data.

Five percent of all inclusions contain low-salinity fluids in combination with low that fall below the temperature-salinity relation of the mixing line between phase-separated vapor and ambient seawater. These inclusions cannot be explained by mixing trends at seafloor pressures because fluids at lower temperatures do not phase separate at all. These inclusions can be interpreted as vapor end-member fluids that have conductively cooled prior to, or during, mixing with seawater. Alternatively, these inclusions may represent re-equilibrated inclusions that were originally produced by mixing of the pure vapor end-member with seawater, but which have since been offset to lower values. Applying conductive cooling or re-equilibration for these inclusions does not contradict our proposed mixing model, and hence, 98% of the inclusions can be explained.

The remaining 2% of inclusions, not adequately explained by our model of subcritical phase separation, contain high-salinity brines ( NaCl eq.) with low formation temperatures (<<300°C). Like the relatively low--low- inclusions, these brine inclusions may have been affected by conductive cooling or re-equilibration. Alternatively, these inclusions may have suffered from reaction between the fluid inclusions and their host crystal. For example, during microthermometric analysis, formation of gypsum after anhydrite on the inclusion walls may increase the salinity by up to ~50% from the original value [54]. However, a number of inclusions were measured to contain fluids with salinities of up to 18.6 wt.% NaCl eq., which cannot be produced by hydration of the host crystals with original salinities of <5.25 wt.% NaCl eq., as proposed for the subcritical region. Subcritical phase separation can account for such highly saline brines if a sufficient amount of vapor has fractionated off the hydrothermal fluid. For instance, a brine with 20 wt.% NaCl eq. can be produced if a mass fraction of 84 wt.% vapor fractionates off and separates from the remaining 16 wt.% brine. This fractionation would mean that the majority of fluids at the NW Caldera should represent vapor-derived fluids. This is clearly not the case considering the brine-dominated nature of the hydrothermal system as seen in the fluid inclusion data and the composition of vent fluids.

During microthermometric studies of fluid inclusions from IODP drillcore samples, of up to 45 wt.% NaCl eq. were measured [7]. These inclusions were considered to record a contribution of sequestered brines from near a magmatic-hydrothermal interface in the deeper parts of the hydrothermal system, originating from either supercritical phase separation or by the condensation of a single-phase magmatic fluid. Our second mixing model specifically calculates an isenthalpic-isobaric mixing path for a supercritically phase-separated brine-vapor pair formed at 340 bar. This model is consistent with supercritical phase separation accounting for the very high-salinity inclusions.

In summary, the rare high-salinity inclusions of up to 18.6 wt.% NaCl eq. are unlikely to be formed by crystals growing from fluids produced under shallow phase separation or conductive cooling conditions. Re-equilibration of the inclusions is also not considered a favorable mechanism to produce the low values in these high inclusions but suggests that these inclusions are either a result of mixing between seawater and a deep-seated fluid produced by supercritical phase separation or condensation of a magmatic brine, as proposed by de Ronde et al. [7]. Our study shows that such a deep-seated brine may occasionally reach the shallowest part of the NW Caldera hydrothermal system and contributes to the composition of fluids that pass through stockwork zones to vent at the seafloor.

5.4. Seawater Entrainment and Sulfate Precipitation

We applied 87Sr/86Sr values to evaluate conductive heating/cooling and/or seawater mixing in the hydrothermal system, similar to previous studies [8, 18, 19]. Values of 87Sr/86Sr in our sulfate separates (0.704852–0.709040) range from near hydrothermal end-member values of 0.70437 to seawater-like values of 0.70918 [47]. This range implies that the individual crystals/samples have precipitated from significantly different mixing regimes from largely hydrothermal fluids and seawater, as to be expected by mass balance equations. Multiple measurements of Sr isotopes for single samples in our study resulted in a considerable range of compositions (e.g., 085-6R1, ) that require highly variable mixing on a temporal scale at one location.

In a setting where subcritical phase separation may produce vapors with very low salinities, using a mass balance approach with specific Sr contents for the hydrothermal end-member fluids may be problematic. The fluid inclusion data for individual samples (with a wide range in ) shows that most crystals grew from vapor-derived fluids, as well as from brine-derived fluids. Vapor-derived fluids, assumed to be depleted in Sr, as well as all other major cations (according to the partitioning of salt) behave differently than seawater-like or brine-rich fluids. That is, Sr isotope mass balances show that such vapor-derived fluids, because of their very low Sr contents, will inherit a seawater isotope signature even after mixing with just a few weight percent of seawater. Consequently, crystals that grow from vapor-derived fluids that then mix with seawater will mimic seawater-derived 87Sr/86Sr signatures, although the hydrothermal component in a fluid-seawater mix may be dominating. The range of seawater-hydrothermal fluid mixing ratios derived from a simple Sr isotope mass balance may thus be overpredicted in systems in which phase separation plays an important role.

Considering the broad range in NW Caldera fluid inclusion within individual samples, we suggest that 87Sr/86Sr values for a single crystal, or even bulk samples, are controlled by what proportion of the crystal grew from brine-derived versus vapor-derived fluids, before mixing with seawater. This contrasts with the approach of using these values to determine hydrothermally dominated as opposed to seawater-dominated regimes. We therefore propose that much of the variation in NW Caldera 87Sr/86Sr values reflects the ephemeral character of phase segregation and the effects of vapor- and brine-derived fluids, rather than simple mixing between hydrothermal fluids and seawater.

For future investigations of this kind in relatively shallow island arc settings, we recommend conducting spatially resolved analytics of 87Sr/86Sr by laser-based methods to circumvent this problem and to determine if 87Sr/86Sr values are homogeneously distributed in individual crystals (cf. [8, 19]). In addition, to address the effect of hydrothermal fluid-seawater mixing versus vapor-brine contribution using Sr isotopes, 87Sr/86Sr values should also be evaluated in close relation with tracers of fluid salinities, or more specific, fluid Sr concentration.

We suggest that provides a better proxy for any mixing regime for NW Caldera fluids. The temperature evolution is essentially the same for mixtures between brine- and vapor-derived fluids, and thus, temperature is not affected by phase separation, as are 87Sr/86Sr values. The problem of using temperature as a proxy for a mixing regime is still that temperature may be affected by conductive heating (or conductive cooling) without noticeably changing any chemical parameter in the fluid. Nevertheless, using temperature to evaluate the mixing regime would lead to conclude that the sulfates precipitated in dominantly hydrothermal conditions. Applying values of inclusion samples to assess the mixing regime results in a hydrothermal component of between 0.5 and 0.9 (i.e., 50–90 wt.% hydrothermal fluid in the mixture) for all samples. By comparison, using individual inclusions results in a hydrothermal component between 0.35 and 0.95 (i.e., 35–95 wt.% hydrothermal fluid in the mixture).

In summary, the Sr data demonstrate that seawater entrainment was widespread and variable throughout the NW Caldera site, consistent with sulfate deposition in chimney walls and in the stockwork immediately below the vent field. This is consistent with other studies that demonstrate that sulfate precipitation as a consequence of mixing between hydrothermal fluids and seawater is commonplace in chimneys of the NW Caldera vent field [31, 41]. While the growth of sulfate crystals in mixing regimes with variable seawater contributions has been shown before in back-arc basins and mid-ocean ridges [8, 18, 19], this study demonstrates that these processes in shallow island arc systems with subcritical phase separation occurring at lower seafloor pressures are more complex.

6. Conclusions

High fluid inclusion entrapment temperatures () and a notable range in inclusion salinities () in hydrothermal sulfate samples provide robust evidence for the ubiquitous nature of phase separation of hydrothermal fluids that formed chimneys and a mineralized stockwork zone at the Brothers NW Caldera vent field. If in fluid inclusions represent an integrated view of fluid salinities in a temporal and spatial fashion, more than half the fluids being expelled on the seafloor are brines (i.e., greater than seawater salinity), a third of the fluids are vapor derived (i.e., less than seawater salinity), with the remainder of seawater salinity. Only rare inclusions possibly represent supercritical phase separation or condensation of a magmatic brine.

Isobaric-isenthalpic mixing models of relatively shallow phase-separated fluids can explain the variability of 98% of all fluid inclusions in our study. They correctly predict the dominant occurrence of brines over vapors in the dataset. The mixing models suggest that the majority of fluids are phase separating in only 10 s to few 100 s mbsf (at a pressure of ~170 bar). Mixing with cold seawater occurs during the very final stage of ascent of these fluids.

A new discovery of this study is that high (up to the depth-to-boiling-point curve) at seafloor pressures and extremely low salinities (down to 0.1 wt.% NaCl eq., in fluid inclusions) in the NW Caldera site are also consistent with phase-separating fluids having sporadically reached the seafloor.

Rare, high-salinity fluid inclusions (>5.25 to 18.6 wt.% NaCl eq.) and possibly halite-bearing three-phased inclusions are unlikely produced by shallow phase separation. These inclusions are thought to represent periodic injection of deep-seated brines into the shallower part of the hydrothermal system. This finding corroborates results from microthermometric studies in IODP drillcore samples and shows that deep-seated brines can be transported to the seafloor and may play a role in the mineralization of chimneys and stockwork mineralization, consistent with the findings of Berkenbosch et al. [42, 66] who used Cu isotopes and trace element mapping to show how some of the mineralization has a magmatic affinity.

Since phase-separation processes are shown to produce salt-poor vapors that result in seawater-type isotope signatures, 87Sr/86Sr mass balance calculations used to evaluate mixing regimes in shallow arc environments should be used with caution. We suggest that the application of Sr isotopes should be used with spatially resolvable methods and in close conjunction with the evaluation of Sr concentration of the mineralizing hydrothermal fluids. Based on fluid inclusion entrapment temperatures, our mixing models suggest that mineralization in the NW Caldera is occurring in a hydrothermally dominated mixing regime with less than 50% of seawater.

Data Availability

The datasets associated with this article are accessible via the information system PANGAEA (https://doi.org/10.1594/PANGAEA.920208). The datasets contain the fluid inclusion raw data (https://doi.pangaea.de/10.1594/PANGAEA.920191) and the 87Sr/86Sr values (https://doi.pangaea.de/10.1594/PANGAEA.920192).

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.


The research cruise TN350 was supported by the U.S. National Science Foundation awards (OCE-1558356) to Susan Humphris and (OCE-1558795) to Anna-Louise Reysenbach. Interridge is thanked for a travel fund that enabled A.D. to participate on cruise TN350. The Bundesministerium für Bildung und Forschung (BMBF) provided funding for cruise SO253 and this work with grant no. 03G0253B. A.D. was also partly supported by the MARUM Center for Marine Environmental Science. We acknowledge the science teams, crews, and captains of the R/V Sonne and R/V Thompson for their work. Thanks are due to the teams of ROV QUEST and Jason II for their efforts to recover rock specimens that made this study possible. Special thanks goes to Andreas Klügel for instructions in the use of the microthermometry stage, to Friedrich Lucassen for his support during lab work and TIMS analysis, Janis Thal for providing gridded multibeam data acquired during cruise SO253, and Patrick Monien for proofreading und useful comments in the early stage of this manuscript. We acknowledge the two anonymous reviewers for their useful comments that improved the overall quality of the manuscript.