Search tips
Search criteria 


Logo of ncommsLink to Publisher's site
Nat Commun. 2017; 8: 16107.
Published online 2017 July 14. doi:  10.1038/ncomms16107
PMCID: PMC5594423

Effect of water activity on rates of serpentinization of olivine


The hydrothermal alteration of mantle rocks (referred to as serpentinization) occurs in submarine environments extending from mid-ocean ridges to subduction zones. Serpentinization affects the physical and chemical properties of oceanic lithosphere, represents one of the major mechanisms driving mass exchange between the mantle and the Earth’s surface, and is central to current origin of life hypotheses as well as the search for microbial life on the icy moons of Jupiter and Saturn. In spite of increasing interest in the serpentinization process by researchers in diverse fields, the rates of serpentinization and the controlling factors are poorly understood. Here we use a novel in situ experimental method involving olivine micro-reactors and show that the rate of serpentinization is strongly controlled by the salinity (water activity) of the reacting fluid and demonstrate that the rate of serpentinization of olivine slows down as salinity increases and H2O activity decreases.

Serpentinization encompasses a series of hydration reactions that occur when ultramafic rocks are exposed to circulating aqueous fluids at temperatures lower than ~400 °C, leading to the formation of serpentine phases (lizardite and chrysotile)±brucite±talc±magnetite, among other minerals1. Serpentinization affects the chemical composition, rheology, magnetic properties, seismic structure and habitability of the shallow lithosphere at slow- and ultraslow-spreading mid-ocean ridges, continental margins and forearc settings of subduction zones2,3,4,5,6,7,8. Serpentinization also influences subduction related processes9 and the geochemical cycling of volatile species (that is, H2O, CO2 and H2S) and fluid-mobile elements10,11,12,13. Recent findings suggest that serpentinization of olivine-rich lithologies also takes place on other planetary bodies, such as the icy moons of Jupiter and Saturn, which in turn, has important implications concerning their habitability14,15.

Despite the pivotal role that serpentinization plays in a number of geological and biological processes and its central role in current origin of life hypotheses16,17,18,19, few experimental studies have attempted to determine the rates of serpentinization reactions and the rates that have been reported diverge widely20,21,22,23,24,25. Furthermore, the environmental factors that affect the reaction rates are incompletely constrained. In the present work, we used synthetic fluid inclusions (SFIs)26,27 as micro-reactors in olivine crystals to monitor the serpentinization reaction with time. We trapped fluid with different initial salinities and followed reaction progress at serpentinization conditions (280 °C). The results show that the rates of olivine serpentinization are strongly influenced by the aqueous fluid salinity. The micro-reactor technique presents several advantages and permits monitoring mineral precipitation and water activity in situ and in real time.


SFI preparation

We trapped SFIs26,27 in gem quality (inclusion-free) natural olivine at 650 °C and 5.5 kbar, conditions where olivine and H2O coexist stably (Fig. 1). The starting H2O-NaCl-MgCl2 fluids in each experiment had a Na/Mg molar ratio similar to that of seawater (8/1) and total salinities of 1, 3.5, 6 and 10 wt%. The capsules were sealed with an arc welder and then placed into cold-seal pressure vessels and held at the desired trapping conditions (650 °C and 5.5 kbar) for 30 days, which allowed the olivine host to anneal and trap some of the fluid as SFIs. The PT conditions of trapping of the SFI were chosen such that the trapped fluid would have a density ~0.8 g cm−3. Following inclusion synthesis, doubly polished wafers were prepared and individual wafers were placed into a furnace at 280 °C and ambient external pressure (~1 bar)—at this temperature, olivine is no longer stable and will react with the aqueous solution to produce serpentinization. At a temperature of 280 °C and constant volume of the SFI, the calculated pressure inside the SFI is ~500 bar.

Figure 1
Photomicrographs and schematic illustrations of the evolution of the micro-reactors through time.

The SFIs can be envisioned as micro-reactors, a reaction vessel in which material can only be exchanged between the aqueous solution and the surrounding host olivine. A major advantage of this approach over most other existing experimental techniques is that changes in mineralogy and fluid composition can be monitored simultaneously in real time as the reaction progresses.

Serpentinization experiments

Following the formation of the SFI in olivine, samples were heated to 280 °C and held at this temperature for 30 to 270 days. For each fluid composition and salinity, several SFI in each olivine wafer were selected to monitor serpentinization reaction progress (Table 1). Each inclusion was photographed and mapped, and each was regularly monitored optically, by microthermometry and by Raman spectroscopy. The polished olivine wafers were briefly (2–5 h) removed from the furnace every 5 days, examined petrographically and analysed by Raman spectroscopy, and then returned to the furnace. Figure 1 (0 days) shows a typical SFI trapped in olivine immediately before starting the heating experiments to initiate the serpentinization reaction.

Table 1
Overview of experimental conditions and alteration mineralogy.

Reaction progress

After a few days of reaction time, small clusters of secondary minerals are observed in the SFI (Fig. 1, 15 days) and the size and number of individual crystals in the SFI progressively increased with time until the inclusions were nearly filled with solids, rendering them opaque (Fig. 1, 45 days).

The onset of the serpentinization reaction varied with the salinity of the starting fluid. Reaction products were first observed after 5 days in ~50% of the SFIs with a starting salinity of 1 wt%. The proportion of the SFI showing evidence of reaction after 5 days decreased to ~20, ~5 and 0% in the SFI with starting salinities of 3.5, 6 and 10 wt%, respectively. At low salinities the proportion of SFI showing reaction products increased rapidly with time until they were ubiquitous in the SFI. For the highest salinity experiments (10 wt%), the first evidence of reaction product formation was observed only after >120 days.

In addition to the experiments in which samples were regularly removed from the furnace for observation and analysis, two samples containing SFI with salinities of 1 and 3.5 wt% were maintained at constant temperature using a Chaixmeca microscope heating stage28. For these samples, the SFI could be monitored without having to cool the sample from 280 °C to room temperature. These experiments were conducted to confirm that repeatedly cycling the samples from 280 °C to room temperature and back to 280 °C did not significantly affect the progression of the experiments. No differences were observed between the thermal cycling experiments and the constant-temperature Chaixmeca experiments.

Mineral characterization

Raman spectroscopic analyses of the SFI (Fig. 2a,b) confirmed that the reaction products consist of serpentine minerals (Srp, lizardite and chrysotile) and brucite (Brc) produced according to the simplified reaction (note that Fe can substitute for Mg in all three mineral phases, but the reaction below is written in terms of the Mg endmembers):

Figure 2
Characterization of the reaction products.
An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m1.jpg

Magnetite is a minor but common byproduct of serpentinization reactions. It forms when the Fe2+ released from olivine is oxidized in the fluid and precipitates magnetite and generates H2:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m2.jpg

In our experiments, magnetite was present in detectable amounts only after >120 days of reaction time. Small magnetite crystals, formed via reaction (2) (but not detectable by optical microscopy) should begin to precipitate inside the micro-reactors soon after reaction (3) starts and the solubility of magnetite in the solution is exceeded. However, we interpret the lack of observable magnetite in the inclusions early in the experiments to reflect the fact that the crystals are not optically resolvable (~0.5 μm) by the methods used here (including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) in Fig. 2c,d). In order for magnetite to become sufficiently large that they can be recognized either optically or during microanalysis, hydrogen must be lost from the SFI. The relatively long reaction time before magnetite is recognizable with optical microscopy in the SFI is attributed to its relatively sluggish precipitation, which may be further limited by competition for iron by serpentine and brucite, or by the slow rate at which hydrogen diffuses out of the inclusions29 when (or if) the inclusion approaches equilibrium.

Mineral composition

The chemical composition of the reaction products was inferred from elemental maps obtained from time-of-flight secondary ion mass spectrometry (TOF-SIMS) analyses of the SFI (Fig. 3 and Supplementary Tables 1 and 2). The ratio of Mg to Fe, defined as the Mg number (Mg#=An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m3.jpg), of brucite is 87–88, whereas serpentine phases have an Mg# ranging from 90 to 94. These values are similar to those of secondary minerals produced during serpentinization of olivine reported from other experimental studies25.

Figure 3
TOF-SIMS elemental maps of the SFI from Fig. 2b.


The serpentinization reaction, even in its simplest form (equation (1)), is the result of several simultaneous reactions that can be summarized as follows. First, olivine dissolves:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m4.jpg

When the concentrations of dissolved species become sufficiently high, serpentine and brucite nucleate and precipitate according to:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m5.jpg
An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m6.jpg

Reactions (4) and (5) consume water and generate hydrogen ions (H+) that in turn promote forsterite dissolution (reaction 3). Eventually, the overall process reaches a steady state whereby the rates of Mg and SiO2(aq) released from olivine dissolution (reaction 3) are balanced by their rates of consumption by serpentine (reaction 4) and brucite (reaction 5) formation.

As the formation of serpentine and brucite consumes H2O, the salinity of the aqueous solution in SFI increases as the reaction proceeds. We used the increase in salinity as a proxy for H2O consumption and reaction progress, to quantify the reaction rates (mol m−2 s−1) (Fig. 4a,b, Supplementary Table 3 and Supplementary Data 1 and 2). The reaction progress is expressed as the extent of reaction (ξ), defined as the number of moles (n) of reactant consumed or product generated per unit time and normalized by the stoichiometric coefficient (An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m7.jpg) for the phase being consumed or produced. Therefore, the reaction rate of a species An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m8.jpg is defined as the time (t) derivative of the extent of reaction:

Figure 4
Reaction progress as a function of H2O consumed through time.
An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m9.jpg

As such, the reaction rate can be defined in terms of the rate of change in the amount of any of the species that participate in the reaction, that is,

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m10.jpg

where the subscripts Fo, W, Srp and Bru correspond to the phases forsterite (olivine), water, serpentine and brucite, respectively. Figure 5a shows the rate of serpentinization as a function of the activity of water. The average rates of reaction were determined to be 2.71 × 10−8, 1.17 × 10−8, 3.37 × 10−9 and 2.86 × 10−10 mol m−2 s−1 for the SFIs with starting salinities of 1, 3.5, 6 and 10 wt%, respectively (Supplementary Table 3). The results show that serpentinization rates decrease rapidly with increasing salinity (or decreasing activity of water), and agree with previous studies that have reported that dissolution rates of forsterite at laboratory conditions (25 °C and 1 bar) decrease as An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m11.jpg decreases30. A rate equation was generated by fitting the log of the reaction rates as a function of log An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m12.jpg with a second-order polynomial:

Figure 5
Experimentally determined serpentinization rates.
An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m13.jpg

where −7.24 is the log of the rate constant at 280 °C and the coefficient 103 represents the apparent reaction order of H2O. The serpentinization rates obtained here agree with the rates of reaction from Wegner and Ernst21 after extrapolating the results to unit activity (pure water) with a log J of −7.51 (Fig. 5b).

These results document the dramatic effect of salinity, as a proxy for water activity An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m14.jpg, on the rate of serpentinization of olivine and are consistent with previous studies that have suggested that salinity can have a significant effect on dissolution kinetics31,32,33 including the dissolution of forsterite at far from equilibrium conditions30,34. The effect of water activity An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m15.jpg on the reaction rate is consistent with the generalized serpentinization reaction given by equation (1), in which the driving force for production of serpentine and brucite is inversely proportional to the cube of the activity of H2O and is also consistent with the dissolution of olivine (equation (3)) being the rate controlling step of the reaction23.

These results provide a geochemical basis to support the previously recognized notion that pervasive serpentinization of oceanic lithosphere (which is dominated by the mineral olivine) requires open-system behaviour, as imprinted by the seawater-dominated isotope signature commonly observed in serpentinized peridotites35,36,37,38. In particular, this requires continual influx of a lower salinity aqueous fluid (seawater) to dilute the serpentinization fluid and allow serpentinization of olivine to continue39. Moreover, salinity may be a major rate-limiting factor where fluid influx is restricted due to the absence of open fractures and fluid migration proceeds along grain boundaries or within nanoscale porosity40,41,42. Furthermore, the effect of salinity on the serpentinization process places constraints on the environments in which serpentinization is likely to occur on Earth and possibly on other planetary bodies where high salinity fluids are believed to exist43,44,45.

The micro-reactor technique used in this study provides a novel and promising tool to monitor fluid-rock reactions in situ and in real time and can be applied to a wide variety of host minerals, reaction products, temperatures and different starting fluid compositions (for example, CaCl2, CO2, SO4 and Al2O3).


SFI preparation

SFIs were trapped in gem quality (inclusion free) natural olivine crystals (Fo85) at 650 °C and 5.5 kbar following procedures described by Sterner and Bodnar26, and similar to Bali et al.27. The crystals were thermally fractured by heating to 700 °C in an oven at atmospheric pressure and then quenching by dropping the crystals into cold distilled water. The crystals were dried overnight in a vacuum oven at 110 °C. Next, one or two 3–5 mm in diameter crystals of olivine were loaded into a platinum capsule (2.5 × 0.5 mm) together with ~100 μl of an H2O-NaCl-MgCl2 fluid (Na/Mg weight ratio of 8/1) with salinities of either 1, 3.5, 6 or 10 wt%. The capsules were sealed with an arc welder. Platinum capsules were placed into cold-seal pressure vessels and held at 650 °C and 5.5 kbar for 30 days, which allowed the olivine host to anneal and trap some of the fluid it was in contact with as synthetic secondary fluid inclusions that were then used as micro-reactors in this study. After the fluid inclusion synthesis, the capsules were cooled to room temperature and the crystals removed from the capsules and mounted in thermal cement for cutting (~100 to 200 μm wafers) and polishing. The sizes of the SFI varied between ~1 and 40 μm, with most being ~10 μm in diameter.

Raman characterization

Raman analyses were performed using a JY Horiba LabRam HR (800 mm) spectrometer, 600 grooves per mm grating, slit width of 150 μm and the confocal aperture was set at 100 μm. A 514.529 nm (green) Laser Physics 100S-514 Ar+ laser was used with a 50 mW output power at the source and ~10 mW at the sample, focused through either a × 40 or × 100 objective, and an electronically cooled open electrode 1,024 × 256 pixel CCD (charge-coupled device) detector. Raman XY mapping was conducted using a JY Horiba LabRam HR confocal Raman spectrometer with a 473 nm laser source and a 600 grooves per mm grating and a 100 μm confocal aperture and × 100 objective.

Electron microprobe analyses

The composition of the olivine host crystal (Supplementary Table 1) used in the experiments was determined using a Cameca SX-50 electron microprobe (EMPA). Accelerating potential was 15 kV, 20 nA current and 1 μm beam size, using natural and synthetic mineral standards for calibration.

Mineral and chemical characterization

To expose the reaction products in the SFI to obtain SEM images, we used a Helios 600 NanoLab FIB-SEM. Chemical analyses of the contents of the SFI after the serpentinization experiments were conducted using TOF-SIMS with a TOF.SIMS 5 IONTOF instrument46,47,48. At each step during the analysis, the ion beam consisting of Bi1+ ion pulses (30 keV ion energy) was first set in the high current bunched mode (20 ns pulse duration, ~3.7 pA measured sample current, ~1 μm lateral resolution) to determine the masses of interest with high mass resolution. The beam was typically raster-scanned over areas of 100 × 100 μm or 200 × 200 μm depending on the region of interest. Then, the beam was set in the burst alignment mode (100 ns pulse duration, ~0.03 pA measured sample current, ~200 nm lateral resolution) to produce high spatial resolution maps of the selected regions for all masses of interest. To expose the SFI in the olivine hosts at various depths we used an O2+ sputtering ion beam (1 keV ion energy, ~300 nA measured sample current) that was raster-scanned over an area of 500 × 500 μm and centred over the region of interest. The resulting sputtering rate was calculated at 0.3 nm s−1 by using a Wyko NT 9100 optical profilometer to estimate the sputtered crater depth. The sample was maintained at a base temperature of −100 °C at all times during the analysis by inductive liquid nitrogen cooling to prevent the water from desorbing from the exposed SFI. A constant current (21 eV electron energy) electron beam was shot on the sample during data acquisition to reduce sample charging. The mass resolution was >3,000 and >100 (m dm−1) for the high current bunched and burst alignment modes, respectively. The TOF-SIMS provided ratios of the number of counts for one mass (isotope) to the number of counts for some other isotope and these counts were converted into a concentration ratio. To do this, the count ratio for Mg+ and Fe+ were converted into concentration ratios using the known Mg to Fe concentration ratio of the olivine host obtained from EMPA and the count ratio obtained by TOF-SIMS analyses for this same area of the crystal (see Supplementary Table 2).

TEM analyses

The sample was polished to a thickness of 30 μm using an Allied Multiprep with 3M Imperial diamond lapping papers of 30, 15, 6, 3 and 1 μm grit sizes. After polishing, the sample was mounted on a copper aperture grid with epoxy. Then the sample was milled in a Fischione 1,010 low-angle argon ion mill at ~10 μm h−1 for several hours until the SFI were exposed. Areas immediately adjacent to the hole became thin enough to be electron transparent. The sample was placed in a JEOL single tilt holder. TEM images were collected at Virginia Tech on a JEOL 2100 TEM operated at 200 kV, using a Gatan Orius CCD camera.

Monitoring reaction progress

The salinity of the SFIs was monitored by measuring the freezing point depression (FPD) of the aqueous phase (temperature of melting of the last ice crystal, Tm) using a Linkam THMSG 600 °C heating and freezing stage. During the serpentinization reaction H2O was removed from the solution and incorporated into the reaction products, whereas Na+ and Cl remained in the solution. Trace amounts of Cl can be incorporated into serpentine, substituting for OH (0.1 to 0.4 wt%), but this is considered to have a negligible effect on salinity in the experiments. Each measurement was corrected for analytical error by measuring the Tm of a pure H2O SFI. The correction factor ranged from 0.8 to 1.0 °C and the precision of the FPD measurements is ~±0.1 °C. The salinity (wt% eq. NaCl) as a function of the FPD was calculated using FPD salinity equation49. The triple point of pure H2O occurs at 0.01 °C and 0.006 bars, where liquid H2O, vapour H2O and solid H2O (in this case, the Ice I polymorph) occur in equilibrium. With increasing salt concentration, the triple point migrates to lower temperature and this effect is referred to as the FPD of the solution. FPD is a colligative property that depends only on the ratio of the number of moles of the solute (NaCl–MgCl2 in our experiments) to that of the solvent H2O. Thus, as H2O was transferred from the aqueous solution into product phases, the ratio of moles of solute to moles of solvent (H2O) increased and the temperature of the triple point decreased. The relationship between the salinity (moles of solute in the solution) and FPD has been determined using data for the ternary system H2O-NaCl-MgCl2 (refs 49, 50).

In our experimental method, we quantify the rates of serpentinization based on the amount of water removed from the aqueous solution and incorporated into hydrous phases. We use the FPD of the aqueous solution as a proxy for the amount of water removed from solution and the precision of our FPD measurement is ±0.1 °C. However, some amount of H2O must be removed from the solution and incorporated into hydrous phases before the salinity is increased by an amount sufficient to produce a 0.1 °C lowering of the FPD. The actual amount of H2O that must be removed from solution to produce a 0.1 °C lowering of the FPD is also a function of the starting salinity. Thus, for example, if the initial salinity is 1 wt% NaCl (FPD=−0.6 °C), ~15% of the H2O in the initial fluid must be consumed by hydration reactions to drive the FPD lower by 0.1 °C to −0.7 °C (corresponding to a salinity of ~1.2 wt% NaCl). However, the same 0.1 °C lowering of the FPD for a starting salinity of 3.5 wt% NaCl requires the consumption of 4.7% of the H2O in the initial fluid, 6 wt% NaCl initial salinity requires a consumption of 2.6% of H2O and 10 wt% NaCl initial salinity requires consumption of only 1.4% of the initial H2O, as the FPD is lowered from −6.6 °C (10 wt% NaCl) to −6.7 °C (10.1 wt% NaCl). Thus, until a sufficient amount of reaction has occurred to consume enough H2O to lower the FPD by ≥0.1 °C, in the absence of observed reaction products in the inclusions we have no way of knowing a priori that the reaction has started and, depending on the starting salinity, the amount of H2O consumed (and, therefore, the amount of reaction that must occur) could require from ~1 to 15% of the initial H2O to be removed from solution. All of the microthermometric data (FPD, salinity and H2O moles consumed) are provided in Supplementary Data 1.

Salinity changes and H2O consumed during serpentinization

To calculate the mass of H2O consumed by the serpentinization reaction as a function of the FPD, we constructed a simple mass balance model that relates the mass of H2O (moles) removed from the aqueous fluid during the reaction to the FPD. The mass of H2O (g) contained in a 10 μm diameter SFI was calculated for 0.1 °C increments of FPD as follows. The total mass of solution in the SFI is given by:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m16.jpg

where VSFI is the volume of an ideal spherical fluid inclusion (5.24·10−10 cm3) and ρSFI is the density of the solution in the SFI (in g cm−3). The density of the aqueous solution in the SFI was calculated using the HOKIEFLINCS H2O_NaCl programme51 and assuming that the system H2O-NaCl provides a reasonable approximation for the relationship between salinity and density for the fluids in the SFI. Accordingly, the densities at 650 °C and 5.5 kbar for salinities of 1, 3.5, 6 and 10 wt% NaCl are ~0.84, ~0.87, ~0.88 and ~0.90 g cm−3, respectively. The masses of NaCl and H2O in the SFI are given by:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m17.jpg
An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m18.jpg

where MNacl and An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m19.jpg are the masses of NaCl and H2O, respectively, in the SFI in grams and XNacl is the mass fraction of NaCl (wt% NaCl/100) in the solution.

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m20.jpg

where An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m21.jpg is the mass of H2O (g), MNacl is the mass of NaCl (g) and XNacl is salinity of the fluid in terms of weight fraction of NaCl.

The amount of H2O consumed by the reaction is estimated from the measured FPD as follows. First, the mass of H2O in the SFI corresponding to the initial salinity and measured FPD is calculated according to equation (12). Then, the amount of H2O in the SFI after some amount of reaction has occurred to lower the FPD by some measureable amount (corresponding to an increase in salinity) is calculated. The loss of H2O from the solution to the product phases corresponding to a given change in salinity was then calculated:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m22.jpg

where, Mw,0 is the mass of H2O in grams at t=0 and Mw is the mass of H2O obtained from equation (9). The number of moles of H2O consumed nw was then calculated as:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m23.jpg

For every initial salinity (1, 3.5, 6 and 10 wt%), the number of moles of H2O consumed was fit as a polynomial expression to determine the number of moles of H2O consumed as a function of the salinity of the SFI. For every FPD measurement in an individual SFI, the number of moles of H2O leaving the system was calculated. All calculations and procedures are provided in Supplementary Data 2.

Rates of reaction and rate equation

Reaction progress was followed in several dozen SFI by monitoring the change in the FPD through time. To follow the reaction progress it was necessary to observe the changes inside the SFI. As reaction progressed, the amount of reaction products increased rendering the inclusion opaque and challenging our ability to follow the reaction progress at high extents of reaction. Moreover, as the H2O inside the SFI is consumed by reaction (1), the internal pressure and volume, as well as the concentrations and compositions inside the inclusion, will be modified from the initial conditions. To avoid the problem of visibility and the effects of changing of the fluid compositions and pressure inside the SFI micro-reactors, our observations where analysed using the initial rate method52. Batch experiments usually require data collection that includes large extents of reaction to circumvent issues related to modifications of the initial conditions (that is, composition, pH, pressure and so on)53. Alternatively, the initial rate method uses numerous short-term experiments with a small number of concentration measurements at equally spaced times. Then, the concentration data versus time data for each experiment are fitted to a function and the slope of that function (rate) is extrapolated to time=0. In our experiments, at time zero little or no brucite+serpentine had formed and no significant water was consumed. The small volume and pressure changes associated with this small extent of reaction did not affect the result. As such, in our experiments, the number of moles of H2O consumed by the serpentinization reaction with time was found from the change in nw with time:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m24.jpg

where r is the rate of H2O consumption in mol s−1, Δnw is the number of moles of H2O consumed by the transformation reaction and Δt is the time An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m25.jpgin seconds. The number of moles of H2O removed from solution and incorporated into hydrous phases as a function of time (t) for every SFI was fit by a linear regression and the rate of reaction was taken as the slope of the regression line53,54 (see Fig. 4).

As the reaction occurs at an interface that separates two phases (solid and liquid), the rate expresses how fast a component is transferred to or from that interface and the area of the interface must be taken into account as follows:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m26.jpg

where J is the flux of H2O from the liquid phase into the product phases, defined as the rate of H2O consumption per unit surface area (mol m−2 s−1), An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m27.jpg is the rate of H2O consumption in mol s−1 and A is the surface area of a 10 μm spherical fluid inclusion (1.26 × 10−9 m−2) (Supplementary Data 2). To generate rate equations, rates obtained as described above were converted into rate equations that summarize the effects of solution compositions52. These experimental rates were fit to an equation that relates the rate to the concentration (a, activity) of H2O, i, raised to a power (n).

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m28.jpg

This equation was transformed to a second order polynomial form:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m29.jpg

The data were fit to an equation of this form using a polynomial function to determine the rate constant (k) and the apparent reaction order (n). The activities (ai) used in this study were calculated with the software code EQ3/6 (ref. 55) using the thermodynamic database for 50 MPa as described in Klein et al.56. Data for the rates, fluxes, activities and extent of reaction of all the experiments are provided in Supplementary Table 3.

Serpentinization rates from the literature

To compare the serpentinization rates obtained in this study with the rate data available in the literature, we converted the published data20,21,22,23,24,25,57 into fluxes (mol m−2 s−1) by applying the Shrinking Particle Model52,58. This model assumes that all particles are spheres and the rate of dissolution equals the rate of conversion. As such, the conversion rates (olivine to serpentine and brucite, per day from the literature) were converted into fluxes (mol m−2 s−1) by determining the rate constant following expressions:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m30.jpg

where α is the conversion rate fraction, kp is the rate constant and t is time (s).

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m31.jpg

where k+ is the flux (mol m−2 s−1), R0 is the particle size radius at time zero and Vm is the molar volume of olivine. R0 was obtained by using the averaging method of Tester et al.52,59, which finds the effective diameter of the grain sizes by the following method:

An external file that holds a picture, illustration, etc.
Object name is ncomms16107-m32.jpg

where De is the effective diameter (m), Dmax is the maximum grain size and Dmin is the minimum grain size diameter.

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Additional information

How to cite this article: Lamadrid, H. M. et al. Effect of water activity on rates of serpentinization of olivine. Nat. Commun. 8, 16107 doi: 10.1038/ncomms16107 (2017).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:
Supplementary Data 1:
Supplementary Data 2:
Peer Review File:


We thank Charles Farley for help with Raman analyses; Chris Winkler, Rui Serra Maia and Jay Tuggle at the Nanoscale Characterization and Fabrication Laboratory (NCFL) for help with TEM and FIB SEM. This material is based upon work supported by the National Science Foundation under Grant OCE-1459433 to R.J.B. and E.M.S. Consejo Nacional de Ciencia y Tecnología (CONACyT), the Virginia Tech Department of Geosciences and Virginia Tech Graduate School provided partial funding to HML during this study. F.K. was supported by The Andrew W. Mellon Foundation Endowed Fund for Innovative Research. We also acknowledge the NSF grant DMR-0923096 used to purchase the TOF-SIMS instrument at Texas Materials Institute, UT Austin.


The authors declare no competing financial interests.

Author contributions H.M.L. developed the concept of SFIs as micro-reactors and performed the experiments reported in this communication. J.D.R. assisted in the interpretation of the kinetic data and geological application of the data. E.M.S. assisted in the interpretation and geological application of the data. F.K. carried out the Raman mapping as well as assisted in the interpretation and geological application of the data. S.U. performed sample preparation and TEM analyses and interpretation of the data. A.D. performed sample preparation for TOF-SIM analyses and interpretation of the data. R.J.B assisted in the experimental design, interpretation and geological application of the data. All authors discussed the results, commented and participated on the writing of the manuscript.


  • Bach W., Garrido C. J., Paulick H., Harvey J. & Rosner M. Seawater-peridotite interactions: first insights from ODP Leg 209, MAR 15 °N. Geochem. Geophys. Geosyst. 5, Q09F26 (2004).
  • Kelley D. S. et al. . An off-axis hydrothermal vent field near the Mid-Atlantic Ridge at 30 degrees N. Nature 412, 145–149 (2001). [PubMed]
  • Früh-Green G. L. et al. . 30,000 Years of hydrothermal activity at the lost city vent field. Science 301, 495–498 (2003). [PubMed]
  • Dyment J., Arkani-Hamed J. & Ghods A. Contribution of serpentinized ultramafics to marine magnetic anomalies at slow and intermediate spreading centres: insights from the shape of the anomalies. Geophys. J . Int. 129, 691–701 (1997).
  • Escartín J., Hirth G. & Evans B. Strength of slightly serpentinized peridotites: implications for the tectonics of oceanic lithosphere. Geology 2942, 1023–1026 (2001).
  • Canales J. P., Tucholke B. E., Xu M., Collins J. A. & DuBois D. L. Seismic evidence for large-scale compositional heterogeneity of oceanic core complexes. Geochem. Geophys. Geosyst. 9, (2008).
  • Kodolanyi J., Pettke T., Spandler C., Kamber B. S. & Gmeling K. Geochemistry of ocean floor and fore-arc serpentinites: constraints on the ultramafic input to subduction zones. J. Petrol. 53, 235–270 (2012).
  • Klein F. et al. . Fluid mixing and the deep biosphere of a fossil Lost City-type hydrothermal system at the Iberia Margin. Proc. Natl Acad. Sci. USA 112, 12036–12041 (2015). [PubMed]
  • Hirth G. & Guillot S. Rheology and tectonic significance of serpentinite. Elements 9, 107–113 (2013).
  • Alt J. C. et al. . The role of serpentinites in cycling of carbon and sulfur: seafloor serpentinization and subduction metamorphism. Lithos 178, 40–54 (2013).
  • Deschamps F., Guillot S., Godard M., Andreani M. & Hattori K. Serpentinites act as sponges for fluid-mobile elements in abyssal and subduction zone environments. Terra Nova 23, 171–178 (2011).
  • Schwarzenbach E. M. et al. . Sulfur geochemistry of peridotite-hosted hydrothermal systems: comparing the Ligurian ophiolites with oceanic serpentinites. Geochim. Cosmochim. Acta 91, 283–305 (2012).
  • Schwarzenbach E. M. et al. . Sources and cycling of carbon in continental, serpentinite-hosted alkaline springs in the Voltri Massif, Italy. Lithos 177, 226–244 (2013).
  • Glein C. R., Baross J. A. & Waite J. H. The pH of Enceladus’ ocean. Geochim. Cosmochim. Acta 162, 202–219 (2015).
  • Vance S. D., Hand K. P. & Pappalardo R. T. Geophysical controls of chemical disequilibria in Europa. Geophys. Res. Lett. 43, 4871–4879 (2016).
  • McCollom T. M. Methanogenesis as a potential source of chemical energy for primary biomass production by autotrophic organisms in hydrothermal systems on Europa. J. Geophys. Res. Planets 104, 30729–30742 (1999).
  • Schulte M., Blake D., Hoehler T. & McCollom T. Serpentinization and it’s implication for life on early Earth and Mars. Astrobiology 6, 364–376 (2006). [PubMed]
  • Martin W., Baross J., Kelley D. & Russell M. J. Hydrothermal vents and the origin of life. Nat. Rev. Microbiol. 6, 805 (2008). [PubMed]
  • Russell M. J., Hall A. J. & Martin W. Serpentinization as a source of energy at the origin of life. Geobiology 8, 355–371 (2010). [PubMed]
  • Martin B. & Fyfe W. S. Some experimental and theoretical observations on the kinetics of hydration reactions with particular reference to serpentinization. Chem. Geol. 6, 185–202 (1970).
  • Wegner W. W. & Ernst W. G. Experimentally determined hydration and dehydration reaction rates in the system MgO-SiO2-H2O. Am. J. Sci. 283 A, 151–180 (1983).
  • Lafay R. et al. . Mineral replacement rate of olivine by chrysotile and brucite under high alkaline conditions. J. Cryst. Growth 347, 62–72 (2012).
  • Malvoisin B., Brunet F., Carlut J., Rouméjon S. & Cannat M. Serpentinization of oceanic peridotites: 2. Kinetics and processes of San Carlos olivine hydrothermal alteration. J. Geophys. Res. Solid Earth 117, 1–14 (2012).
  • Ogasawara Y., Okamoto A., Hirano N. & Tsuchiya N. Coupled reactions and silica diffusion during serpentinization. Geochim. Cosmochim. Acta 119, 212–230 (2013).
  • McCollom T. M. et al. . Temperature trends for reaction rates, hydrogen generation, and partitioning of iron during experimental serpentinization of olivine. Geochim. Cosmochim. Acta 181, 175–200 (2016).
  • Sterner S. M. & Bodnar R. J. Synthetic fluid inclusions in natural quartz I. Compositional types synthesized and applications to experimental geochemistry. Geochim. Cosmochim. Acta 48, 2659–2668 (1984).
  • Bali E., Audétat A. & Keppler H. Water and hydrogen are immiscible in Earth’s mantle. Nature 495, 220–222 (2013). [PubMed]
  • Poty B., Leroy J. & Jachimowicz L. Un nouvel appareil pour la mesure des températures sous le microscope: l’installation de microthermométrie Chaixmeca—Google Scholar. Bull. Soc. Fr. Minéral. Cristal 99, 182–186 (1976).
  • Demouchy S. & Mackwell S. Mechanisms of hydrogen incorporation and diffusion in iron-bearing olivine. Phys. Chem. Miner. 33, 347–355 (2006).
  • Olsen A. A., Hausrath E. M. & Rimstidt J. D. Forsterite dissolution rates in Mg-sulfate-rich Mars-analog brines and implications of the aqueous history of Mars. J. Geophys. Res. Planets 120, 388–400 (2015).
  • Schott J., Pokrovsky O. S. & Oelkers E. H. The link between mineral dissolution/precipitation kinetics and solution chemistry. Rev. Mineral. Geochemistry 70, 207–258 (2009).
  • Casey W. H. & Westrich H. R. Control of dissolution rates of orthosilicate minerals by divalent metal–oxygen bonds. Nature 355, 157–159 (1992).
  • Yardley B. W. D. & Bodnar R. J. Fluids in the continental crust. Geochem. Perspect 3, 1 1–12 (2014).
  • Morrow C. P., Olsen A. A. & Kubicki J. D. Quantum mechanical modeling of hydrolysis and H2O-exchange in Mg-, Ca-, and Nisilicate clusters: implications for dissolution mechanisms of olivine minerals. Am. Mineral. 99, 2303–2312 (2014).
  • Boschi C., Dini A., Früh-Green G. L. & Kelley D. S. Isotopic and element exchange during serpentinization and metasomatism at the Atlantis Massif (MAR 30°N): insights from B and Sr isotope data. Geochim. Cosmochim. Acta 72, 1801–1823 (2008).
  • Vils F., Tonarini S., Kalt A. & Seitz H. M. Boron, lithium and strontium isotopes as tracers of seawater-serpentinite interaction at Mid-Atlantic ridge, ODP Leg 209. Earth Planet. Sci. Lett. 286, 414–425 (2009).
  • Decitre S. et al. . Behavior of Li and its isotopes during serpentinization of oceanic peridotites. Geochem. Geophys. Geosyst. 3, 1–20 (2002).
  • Malvoisin B. Mass transfer in the oceanic lithosphere: serpentinization is not isochemical. Earth Planet. Sci. Lett. 430, 75–85 (2015).
  • Rouméjon S. & Cannat M. Serpentinization of mantle-derived peridotites at mid-ocean ridges: mesh texture development in the context of tectonic exhumation. Geochem. Geophys. Geosyst. 15, 2354–2379 (2014).
  • Plumper O., Royne A., Magraso A. & Jamtveit B. The interface-scale mechanism of reaction-induced fracturing during serpentinization. Geology 40, 1103–1106 (2012).
  • Schwarzenbach E. M., Caddick M. J., Beard J. S. & Bodnar R. J. Serpentinization, element transfer, and the progressive development of zoning in veins: evidence from a partially serpentinized harzburgite. Contrib. Miner. Petrol 171, 1–22 (2016).
  • Tutolo B. M., Mildner D. F. R., Gagnon C. V. L., Saar M. O. & Seyfried W. E. Nanoscale constraints on porosity generation and fluid flow during serpentinization. Geology 44, 103–106 (2016).
  • Zolensky M. E. et al. . Asteroidal water within fluid inclusion-bearing halite in an H5 chondrite, Monahans (1998). Science 285, 1377–1379 (1999). [PubMed]
  • McEwen A. S. et al. . Seasonal flows on warm Martian slopes. Science 333, 740–743 (2011). [PubMed]
  • Martínez G. M. & Renno N. O. Water and brines on Mars: current evidence and implications for MSL. Space Sci. Rev. 175, 29–51 (2013).
  • Chou H., Ismach A., Ghosh R., Ruoff R. S. & Dolocan A. Revealing the planar chemistry of two-dimensional heterostructures at the atomic level. Nat. Commun. 6, 7482 (2015). [PMC free article] [PubMed]
  • Zu C. et al. . Breaking down the crystallinity: the path for advanced lithium batteries. Adv. Energy Mater. 6, 1–9 (2015).
  • Colleary C. et al. . Chemical, experimental, and morphological evidence for diagenetically altered melanin in exceptionally preserved fossils. Proc. Natl Acad. Sci. USA 112, 12592–12597 (2015). [PubMed]
  • Bodnar R. J. Revised equation and table for determining the freezing point depression of H2O-NaCl solutions—Google Scholar. Geochim. Cosmochim. Acta 57, 683–684 (1993).
  • Dubois M. & Marignac C. The H2O-NaCl-MgCl2 ternary phase diagram with special application to fluid inclusion studies. Econ. Geol. 92, 114–119 (1997).
  • Steele-Macinnis M., Lecumberri-Sanchez P. & Bodnar R. J. HOKIEFLINCS_H2O—NACL: A Microsoft Excel spreadsheet for interpreting microthermometric data from fluid inclusions based on the PVTX properties of H2O—NaCl. Comput. Geosci. 49, 334–337 (2012).
  • Rimstidt J. D. Geochemical Rate Models: An Introduction to Geochemical Kinetics Cambridge Univ. Press (2013).
  • Rimstidt J. D. & Newcomb W. D. Measurement and analysis of rate data: the rate of reaction of ferric iron with pyrite. Geochim. Cosmochim. Acta 57, 1919–1934 (1993).
  • Laidler K. Theories of Chemical Reaction Rates McGraw-Hill (1969).
  • Wolery T. EQ3/6: A software package for geochemical modeling of aqueous systems: package overview and installation guide (version 7.0) (Lawrence Livermore National Laboratory, Livermore, CA, USA, 1992).
  • Klein F., Bach W. & McCollom T. M. Compositional controls on hydrogen generation during serpentinization of ultramafic rocks. Lithos 178, 55–69 (2013).
  • Andreani M., Daniel I. & Pollet-Villard M. Aluminum speeds up the hydrothermal alteration of olivine. Am. Miner. 98, 1738–1744 (2012).
  • Burkin A. R. Chemical Hydrometallurgy: Theory and Principles Imperial College Press (2001).
  • Tester J. W., Worley W. G., Robinson B. A., Grigsby C. O. & Feerer J. L. Correlating quartz dissolution kinetics in pure water from 25 to 625 °C. Geochim. Cosmochim. Acta 58, 2407–2420 (1994).

Articles from Nature Communications are provided here courtesy of Nature Publishing Group