Skip to main content

A single column separation method for barium isotope analysis of geologic and hydrologic materials with complex matrices


The increasing significance of barium (Ba) in environmental and geologic research in recent years has led to interest in the application of the Ba isotopic composition as a tracer for natural materials with complex matrices. Most Ba isotope measurement techniques require separation of Ba from the rest of sample prior to analysis. This paper presents a method using readily available materials and disposable columns that effectively separates Ba from a range of geologic and hydrologic materials, including carbonate minerals, silicate rocks, barite, river water, and fluids with high total dissolved solids and organic content such as oil and gas brines, rapidly and without need for an additional cleanup column. The technique involves off-the-shelf columns and cation exchange resin and a two-reagent elution that uses 2.5 N HCl followed by addition of 2.0 N HNO3. We present data to show that major matrix elements from almost any natural material are separated from Ba in a single column pass, and that the method also effectively reduces or eliminates isobaric interferences from lanthanum and cerium.


Barium (Ba) is a critical element of interest in studies of oceanic biogeochemical cycling and the diagenetic alteration of sediments (e.g., [1,2,3,4,5]), and the isotopic composition of Ba has been used to understand element cycling in marine, riverine, lacustrine and pedogenic environments [6,7,8,9,10,11,12,13,14,15,16]. More recently, Ba concentrations [17,18,19] and its isotope signatures [20] have been studied as indicators of Ba sources and mobility related to subsurface water–rock interactions and the fate and transport of contaminants in drilling mud, solid waste and produced waters related to hydraulic fracturing of oil and gas wells.

In order to obtain the Ba isotope composition of a sample via thermal ionization mass spectrometry (TIMS) or multi-collector inductively coupled plasma mass spectrometry (MC-ICPMS), the Ba must be separated from other elements in the sample. High levels of matrix elements will interfere with the ionization and transmission of Ba, and elements with directly overlapping isotopic masses (isobaric interferences, including Xe, La, and Ce) change the measured Ba isotope ratios unless they are removed prior to and/or corrected for during isotope ratio measurement [21,22,23,24,25,26,27]. While Xe is only introduced during MC-ICPMS measurement from the atmosphere and Ar gas source, all other interferents (matrix major elements, La, and Ce) must be removed via chemical processing prior to analysis by TIMS or MC-ICPMS. For samples with high ratios of Ba to matrix and/or isobaric interferents (e.g., high-Ba formation waters or samples of barite, BaSO4), the chemical separation of Ba is relatively straightforward [6, 14, 20, 24, 28, 29]. However, silicate rock samples can have high rare earth element (REE) to Ba ratios, necessitating robust procedures for removing La and Ce prior to analysis.

Various methods for the separation of Ba for isotopic analysis have been reported in detail. Miyazaki et al. [30] described an elution method using 2.5 N hydrochloric acid (HCl) to separate matrix major elements followed by 1.5 N nitric acid (HNO3) to enhance separation of Ba from REE in silicate rock samples. However, two major isotopes of Ba (136Ba and 138Ba) were excluded from the MC-ICPMS measurement protocol due to interferences from 138La, 136Ce, and 138Ce [30]. Nan et al. [23] used a similar elution method (3.0 N HCl followed by 4.0 N HNO3), using two successive columns to fully remove interferents. Elution curves for La and Ce were not reported in this study. For separation of Ba from a range of silicate and carbonate rock types, Nan et al. [31] and Zeng et al. [25] reported using up to three column passes to fully separate Ba from the REE. Tian et al. [28] described a method for separation of Ba from barite using 3 N HCl and 3 N HNO3. Because barite consists of ~ 59% Ba by weight, matrix and isobaric interferences are negligible, and the Ba could be purified in a single step.

This paper describes sample preparation and ion exchange chromatography separation of Ba using a two-reagent elution in a single off-the-shelf column that is suitable for the isotopic analysis of Ba in samples with complex matrices and a range of REE/Ba ratios that are encountered in geologic and environmental samples. In addition to barite and seawater, we demonstrate that the method is applicable to low-Ba silicates and carbonate rocks, river water, and high total dissolved solid (TDS) fluids containing organic compounds, such as oil and gas produced waters and other brines. The method uses readily available polypropylene disposable columns and cation exchange resin. We document that this method successfully separates Ba from matrix major elements, as well as La and Ce, in one step, thus necessitating minor to negligible correction for these isobaric interferences during measurement by MC-ICPMS.

Experimental methods

Materials and reagents

Experimental work was conducted at the University of Pittsburgh under clean lab conditions. Ultrapure 18.2 MΩ Milli-Q water (MQW) was used for acid dilutions and for washing solids. Teflon beakers and other sample vessels were acid washed, and most sensitive procedures (e.g., sample dissolution, spiking, evaporation, and column separations) were carried out in ULPA-filtered laminar flow hoods. Cation exchange resin (Bio-Rad AG® 50W-X8 200–400 mesh) was pre-cleaned with repeated, sequential applications of MQW, 2% HNO3, and 6 N HCl, which were added to the resin in a fluorinated ethylene propylene (FEP) bottle, shaken vigorously, and decanted after settling. Ultrapure (Fisher Optima™) HCl, HNO3, hydrofluoric acid (HF) and hydrogen peroxide (H2O2) were used for sample dissolution and preparation, resin/column preparation, elution and organic matter removal. Ultrapure (J.T. Baker™ Ultrex™) sodium carbonate (Na2CO3) was used for barite dissolution and seawater Ba precipitation. The Na2CO3 solution was further purified to remove Ba by using a modification of the method described by Foster et al. [32]; in this method, 1.15 M Na2CO3 solution is reacted with a 1.1 N calcium chloride solution (CaCl2; Thermo Fisher Alfa Aesar™ 99.99% purity), which precipitates Ba along with calcium carbonate, and results in an approximately 1 M Na2CO3 solution with low Ba levels.

To prevent sample carryover, we used disposable polypropylene gravity flow ion exchange columns (Bio-Rad Poly-Prep® Chromatography Columns) and disposed of the column and the resin after each use. The columns have a 2 mL bed volume with a 10 mL reservoir and a porous 30 µm polyethylene bed support frit (Fig. 1).

Fig. 1
figure 1

Bio-Rad Poly-Prep® Chromatography Columns used in Ba separation experiments

Samples for Ba separation experiments

Verification of method applicability was conducted using natural materials and standards which reflect a range of matrices encountered in geologic and environmental samples. Silicate rocks (shale, sandstones, igneous rocks) can be problematic for Ba isotopic analysis because of potentially high levels of REE interferents relative to Ba. Carbonate minerals such as limestone and marine shells have high Mg and Ca content, and some freshwater aragonite shells can have high REE relative to Ba [33]. Oil and gas produced waters have high TDS, can contain hydrocarbons and other organic compounds, and can vary greatly in Ba content (0.25 to > 10,000 ppm [34, 35]). For the experiments reported here, fluid samples included seawater (National Research Council Canada seawater reference material NASS-6 collected off the coast of Nova Scotia, Canada), surface water from the Ohio River in Pennsylvania, USA (OR-1204-1a), and produced water from the Marcellus Formation in Greene County, Pennsylvania (M4TFA0518; also used as an internal lab Ba isotope standard). The carbonate sample (PA-LO-102714) was a freshwater mussel shell collected in western Pennsylvania. Silicate rock samples include volcanics from southern Alaska [36], an organic-rich black shale (Marcellus Shale), and USGS standards BCR-2 basalt and AGV-1 andesite [37]. Selected data for these samples are provided in Additional file 1: Table S1.

Sample dissolution/Ba concentration methods

The carbonate mussel shell was powdered and dissolved in 2.5 N HCl. The silicate rock samples were dissolved under clean lab conditions using HF and HNO3 in a microwave digestion system (Milestone Ethos®). Following digestion at 200 °C, the samples were treated with aqua regia and concentrated HNO3 to remove fluorides.

Barite was dissolved following a method based on earlier work [14, 28, 38, 39]. In brief, powdered barite was reacted with 1 M purified Na2CO3 solution (1 mL/10 mg barite) at 95 °C in a sealed Savillex® vial, leading to exchange of SO42− with CO32−. After centrifuging, the supernatant was discarded, and the reaction was repeated twice on the remaining residue for a total of 72 h, resulting in a precipitate of barium carbonate (BaCO3). The precipitate was rinsed with MQW, evaporated to dryness, and the carbonate dissolved in 3 N HNO3. The sample was dried down, redissolved, and diluted to 2% HNO3 for elemental analysis.

Because Ba is present at low concentrations (30–175 nM) in seawater [1,2,3, 7, 9,10,11], it must be concentrated prior to column separation in order to reduce the matrix load. Barium was extracted from seawater (NASS-6) using a method modified from Foster et al. [32] and others [7, 10,11,12]. An aliquot of seawater containing ~ 2 µg Ba (~ 200 mL) was spiked with a 137Ba–135Ba solution (see next section) so that isotope mass fractionation associated with all subsequent processing could be corrected for during measurement. The aliquot was reacted with sufficient 1 M Na2CO3 solution (purified to remove Ba) to result in precipitation of calcium carbonate, which scavenges Ba from solution. This occurs at a molar ratio of CO32− to Ca2+ of approximately 6:1. Once saturation was reached, the reaction occurred nearly instantaneously, and the crystal-fluid mixture was then stirred for 5 min to ensure full removal of Ba. The sample was then decanted into 50 mL centrifuge tubes, rinsing any residue into the tubes using MQW, and centrifuged at 4000 rpm for 5 min. The supernatant was removed by pipette, preserving the precipitate (with most of the Ba) in the centrifuge tube. The precipitate was transferred to a Savillex vial using MQW and a vortex mixer as needed and evaporated to near dryness at 90 °C. The carbonate precipitate was then dissolved in 4 N HCl, evaporated to dryness, and redissolved in 0.5 mL 2.5 N HCl for cation column separation. This method can also be used for separation of Ba from other high TDS fluids and brines with a low Ba concentration and sufficient Ca for precipitation.

Barium double spike

A calibrated double Ba isotope spike was added to dissolved samples prior to loading of the columns and was used to correct for mass fractionation during MC-ICPMS analysis. Previous work [6,7,8,9, 23, 24, 30, 40,41,42,43] has shown that this can correct for mass fractionation during chemical processing, thus obviating the need for full (> 99%) recovery of Ba from the column. In addition, use of a double spike reduces the effects of non-isobaric matrix interferents [24]. Based on error analysis using the double spike tool of Rudge et al. [44], the 137Ba–135Ba isotope pair was chosen because it optimizes the precise determination of 138Ba/134Ba, whether obtained by direct measurement or calculated based on the measured fractionation in the 138Ba/136Ba ratio. The double spike was made using 135Ba- and 137Ba-enriched carbonate salts from the National Isotope Development Center at Oak Ridge National Laboratory and calibrated by running mixtures of double spike and NIST 3104a. Further details of spike calibration can be found in the Supporting Information of Tieman et al. [20].

Cation column separation procedure

For some samples with very low REE/Ba (carbonate shell PA-LO-102714 and produced water M4TFA0518), a mixed REE standard was added prior to column experiments in order to test the separation of Ba more rigorously from the isobaric interferents La and Ce. In the experiments reported here, approximately 3 mL of cleaned AG50W 200–400 mesh cation resin suspended in MQW was incrementally added by pipette to the acid cleaned columns until the resin bed was just above the base of the 10 mL column reservoir. The resin bed was acidified with 3 mL of HCl at the normality used for that elution experiment, and the resin adjusted to ensure no air bubbles were left and that the top of the resin bed was aligned with the base of the column reservoir. The resin-filled columns were then cleaned with 5 mL of 6.0 N HCl and equilibrated with 4.5 mL of HCl at the appropriate elution normality, added incrementally.

Prior to loading the sample into the columns, residual organics were removed from the sample by adding 1 mL of 50% H2O2 drop wise, followed by evaporation to dryness. The sample was then redissolved in 0.5 mL of HCl at the appropriate elution normality, allowed to equilibrate for a minimum of 2 h, and placed in an ultrasonic bath for 30 min to ensure that it was dissolved to the fullest extent possible. Samples were centrifuged in 2 mL polypropylene centrifuge tubes at 4000 rpm for 10 min to prevent any insoluble precipitates from being loaded onto the cation exchange resin bed. The supernatant was added to the prepared column using a pipette with a 1 mL polypropylene tip. The sample matrix (excluding Ba and the REE) was eluted by adding HCl of the appropriate normality in successive increments, and the eluent was collected for elemental analysis. This was followed by the addition of either HCl or HNO3, depending on the experiment, and these samples (containing primarily Ba and the REE) were collected for elemental analysis.

Analysis of column calibration cuts

Elemental chemistry of column calibration cuts was determined on a Thermo Element XR® sector field (SF) ICP-MS at the at the National Energy Technology Laboratory (NETL). The SF-ICP-MS sample introduction system used a Glass Expansion Ezylok Micromist nebulizer and Twister spray chamber. The SF-ICP-MS gas flow rates and mass spectrometer voltages were optimized daily in low resolution mode to obtain the maximum signal intensity while minimizing doubly charged and oxide species. Barium, Sr, and REE were analyzed using the low resolution slit (∆m/m = 300), while second row non-metals and Fe were analyzed using the medium resolution slit (∆m/m = 4000).

All samples were spiked with 10 μg/L indium as an internal standard prior to analysis to account for matrix differences. Samples which contained high TDS were minimally diluted with 2% HNO3 prior to analysis to minimize matrix effects affecting the ionization in the plasma. While quantitation was not required to generate the elution curves, limits of detection (LOD) were calculated for all analytes as 3 times the standard deviation of the blank from three consecutive blank runs. These detection limits are listed in Table 1. Data from these experiments are reported in Additional file 1: Tables S2–S5.

Table 1 Limits of detection (LOD) for all isotopes monitored in this work by SF-ICP-MS

Ba isotope mass spectrometry method

A subset of samples was analyzed for Ba isotopic composition on a Thermo Neptune Plus® MC-ICPMS at the University of Pittsburgh, part of a joint NETL-University of Pittsburgh facility. For these samples, the separated Ba cut was evaporated to dryness at 90 °C, dissolved in 1 mL of concentrated HNO3, and sonicated for 10 min before again being evaporated to dryness. The sample was redissolved in 2% HNO3 (2 mL per µg of Ba), sonicated in the sealed container for 60 min, and transferred to an acid cleaned 15 mL centrifuge tube for MC-ICPMS analysis.

We report the isotope ratio of Ba as δ138Ba, which is the permil deviation of the 138Ba/134Ba ratio of a sample from that of NIST Standard Reference Material 3104a:

$${\delta }^{138}Ba={10}^{3}\left[\frac{{\left({}^{138}Ba/{}^{134}Ba\right)}_{sample}}{{\left({}^{138}Ba/{}^{134}Ba\right)}_{3104a}}-1\right]$$

Samples were introduced into the MC-ICPMS using an ESI Apex® desolvating nebulizer. Barium isotopes (134Ba, 135Ba, 136Ba, 137Ba, 138Ba) were measured simultaneously on five Faraday cups using 1011 Ω resistors. Isobaric interferences from Xe, La and Ce were monitored by measurement of 131Xe, 139La and 140Ce on three additional Faraday cups. Xenon, which is present in the atmosphere and can have variable concentrations in the argon gas used in MC-ICP-MS analysis, has interfering masses at 134 and 136 (10.44% and 8.87% of total Xe). Cerium is an isobaric interferent of Ba at mass 136 (136Ce = 0.185% of total Ce), and both Ce and La are isobaric interferents at mass 138 (138Ce = 0.251% of total Ce and 138La = 0.09% of total La). Typical 138Ba intensities were 30–40 V for the Apex® desolvating nebulizer (~ 200 V/ppm Ba). No apparent systematic variation in δ138Ba was observed within this range of signal intensity. Mass fractionation of the sample was corrected for by iterative normalization to the 137Ba/135Ba of the double spike using an exponential law [45]. The detector configuration and additional details of the mass spectrometric methods can be found in Supporting Information in Tieman et al. [20].

Results and discussion

Optimal separation of Ba

Combinations of different normalities of HCl (1.5–2.5 N) and HNO3 (1–3 N) were tested to optimize the purity of the Ba cut while minimizing the volume and concentration of reagent necessary to elute. Results of selected experiments are reported in Additional file 1: Tables S2–S5; Figures S1, S2. Use of 2.5 N HCl up to the beginning of Ba elution resulted in the optimal separation of major matrix elements from Ba (Fig. 2). Experiments in which only HCl was used for elution (including of Ba) resulted in either overly large volumes of acid being required to fully capture the Ba (e.g., 2.0 N HCl; Additional file 1: Figure S1a) or significant overlap of Ba with the isobaric interferents La and Ce (e.g., 2.5 N HCl; Additional file 1: Figure S1b). van Zuilen et al. [24] report using 2.0 N or 6.4 N HCl to elute Ba, but do not report elution curves for La or Ce. The absence of these peaks during the TIMS measurement [24] may reflect the lower ionization efficiency of these elements relative to Ba rather than separation from the Ba cut. Previous work demonstrated that a better separation of Ba from the REE could be achieved by eluting with HNO3 [23, 28, 30]; therefore, we tested different normalities of HNO3 for optimal elution of Ba while minimizing overlap of La and Ce. Use of higher normality acid results in a tighter Ba elution peak but greater overlap of REE with Ba (e.g., Additional file 1: Figure S2a,b). The optimal concentration of HNO3 to elute Ba from the columns used in these experiments was 2.0 N (Fig. 2).

Fig. 2
figure 2

Column matrix separation for major elements and Ba for a 60 mg CaCO3 (freshwater mussel shell) sample. Cerium and La were not fully recovered even after 60 mL, so their % elution values are likely overestimated

The separation of Ba from REE isobaric interferents was further tested by analyzing Ba elution curves in detail using a dissolved volcanic rock with La/Ba of ~ 0.015 and Ce/Ba of ~ 0.033, and a Marcellus Shale produced water sample in which La and Ce were added to increase La/Ba and Ce/Ba to ~ 0.67 (Fig. 3). Based solely on the measured signal intensities (reported in counts per second) by SF-ICP-MS and the abundance of the measured isotope, the overlap of La with the Ba cut is 0.007% for the volcanic sample and 0.1% for the La-enriched produced water sample. Similarly, the overlap of Ce with the Ba cut is 0.03% and 1.1% for the volcanic sample and Ce-enriched produced water, respectively. In the latter case, such an overlap in Ce would necessitate a correction of ~ 0.3‰ to the δ138Ba (for the measurement method described here), which can be done accurately. However, most geologic materials have Ce/Ba ratios (Table 2) considerably lower than the artificially Ce-enriched produced water reported here (Ce/Ba ≈ 0.25). Unusual natural samples with very high Ce/Ba could be run through the same column twice to remove any excess Ce.

Fig. 3
figure 3

Separation of Ba from isobaric interferents La and Ce during elution of 2 N HNO3. The y-axis is the signal intensity (in counts per second by SF-ICP-MS) per mL of each column cut, recalculated to reflect the total element signal rather than just the measured isotope. The La and Ce signals are multiplied by the factors shown on the curves; if plotted at their actual intensities, the curves would be barely distinguishable above baseline. The La/Ba and Ce/Ba of the volcanic silicate sample (a) prior to loading was ~ 0.014 and ~ 0.033, respectively, while La and Ce were added to the produced water sample (b) prior to loading to generate high La/Ba and Ce/Ba ratios of ~ 0.67

Table 2 La/Ba and Ce/Ba ranges for various geologic and hydrologic materials

In the optimal separation procedure, a total of 20 mL 2.5 N HCl is used to elute the major matrix elements. Barium is eluted using 12 mL of 2.0 N HNO3 (Fig. 2). Major elements were effectively removed prior to the Ba cut; the results were similar for all samples tested (silicate, barite, and fluid samples). The REE isobaric interferents (La and Ce) only begin to elute from the columns after collection of the Ba cut. Additional Ce and La was added to the dissolved aliquot of carbonate (PA-LO-102714) prior to loading in the column to verify separation of Ba from these mass interferents (Fig. 2). Because HNO3 so effectively delays release of La and Ce, these elements were not fully removed from the columns even after adding 40 mL of 2 N HNO3 (following 20 mL of 2.5 N HCl). For all samples, including biogenic carbonate, silicate volcanic rock, organic-rich shales, and high TDS oil and gas produced brine, use of the method resulted in the release of major elements (e.g., Ca, Na, K, Al) from the sample matrix during the HCl elution. Barium was eluted during the nitric acid step, and La and Ce were effectively held on the resin until the bulk of the Ba was eluted. Using gravity flow, the total time from sample loading to removal of the Ba cut is 3–4 h. With a 24-sample column holder, processing 48 samples/day is easily achievable by a single user. The rate-limiting step is preparing the samples for column work (dissolution, evaporation, precipitation of Ba when needed).

The maximum procedural blank for Ba, determined by the limits of detection by SF-ICP-MS, was 2 ng of Ba, leading to a maximum δ138Ba uncertainty of ± 0.002‰ for a 2 µg sample. This assumes the maximum isotopic difference between the sample and blank of 2‰, equivalent to the entire range of δ138Ba values measured to date [20]. This is well within the in-run uncertainty of ± 0.03‰ or better for the MC-ICPMS analysis.

Verification from Ba isotope analysis

Barium isotope analysis of separated samples further demonstrates effective separation. Two silicate rock standards (BCR-2 basalt and AGV-1 andesite) and an internal produced water standard (M4TFA0518) were analyzed after separations using both our previous method [20] and the new method described here. With the previous method (smaller Teflon columns eluted with 2.0 N HCl, with silicates put through a second time to remove REE interferents), we obtained good agreement of the measured δ138Ba with published standard values for silicates [20]. The results of the current analyses are shown in Table 3. Because the isobaric interferences from 136Ce, 138Ce, and 138La are monitored continuously during the run using 139La and 140Ce, we can calculate the total correction to the δ138Ba from these interferents. As shown in Table 3, the maximum correction required for Ce is 0.002‰, which is an order of magnitude smaller than the typical in-run uncertainty. The correction required for La is three orders of magnitude smaller than the typical in-run uncertainty. The ratio of total La to total Ba transmitted into the MC-ICPMS after column separation was ~ 3 × 10–4 of the pre-column ratio, while the Ce/Ba ratio after column separation was 9–14 × 10–4 of the pre-column value. Given that the interfering isotopes only make up 0.09% and 0.49% of the interferent masses measured at 139La and 140Ce, respectively, this amounts to a negligible or reliably correctable signal.

Table 3 Interference corrections and δ138Ba values for selected standards

The δ138Ba value of USGS standard BCR-2 (Table 3) that was processed using the methods described in this paper agree well with values reported previously by Nan et al. [23, 31], An et al. [46] and Tieman et al. [20]. The AGV-1 δ138Ba value obtained here agrees with values reported by Nan et al. [23, 31], van Zuilen et al. [24], An et al. [46] and Tieman et al. [20].


We describe a method of Ba separation by cation exchange using readily available, disposable columns suitable for varied sample matrices, including silicates, carbonates, seawater, sulfates, brines, and produced waters. Eluting with 2.5 N HCl ensures the removal of major elements while Ba is still in the column, and the subsequent elution with 2.0 N HNO3 separates Ba effectively from REE mass interferents. We show that matrix and isobaric interferents can be removed from most sample types in a single step, with no need for an additional cleanup column. In addition, the use of disposable columns prevents possible cross contamination when dealing with complex sample matrices. The moderate column size and elution volumes allow for rapid (3–4 h) simultaneous separation of multiple samples.

Availability of data and materials

The dataset supporting the conclusions of this article is included within the article or in Additional file 1.


  1. Church TM, Wolgemuth (1972) Marine barite saturation. Earth Planet Sci Lett 15:35–44

    Article  Google Scholar 

  2. Bacon MP, Edmond JM (1972) Barium and GEOSECS III in the southwest Pacific. Earth Planet Sci Lett 16:66–74

    Article  Google Scholar 

  3. Chan LH, Edmond JM, Stallard RF, Broecker WS, Chung YC, Weiss RF, Ku TL (1976) Radium and barium at GEOSECS stations in the Atlantic and Pacific. Earth Planet Sci Lett 32:258–267

    Article  Google Scholar 

  4. Griffith EM, Paytan A (2012) Barite in the ocean—occurrence, geochemistry and palaeoceanographic applications. Sedimentology 59:1817–1835

    Article  Google Scholar 

  5. Carter SC, Paytan A, Griffith EM (2020) Toward an improved understanding of the marine barium cycle and the application of marine barite as a paleoproductivity proxy. Minerals 10:421

    Article  Google Scholar 

  6. von Allmen K, Böttcher ME, Samankassou E, Nägler TF (2010) Barium isotope fractionation in the global barium cycle: first evidence from barium minerals and precipitation experiments. Chem Geol 277:70–77

    Article  Google Scholar 

  7. Horner TJ, Kinsley CW, Nielsen SG (2015) Barium-isotopic fractionation in seawater mediated by barite cycling and oceanic circulation. Earth Planet Sci Lett 430:511–522

    Article  Google Scholar 

  8. Bullen T, Chadwick O (2016) Ca, Sr and Ba stable isotopes reveal the fate of soil nutrients along a tropical climosequence in Hawaii. Chem Geol 422:25–45

    Article  Google Scholar 

  9. Cao Z, Siebert C, Hathorne EC, Dai M, Frank M (2016) Constraining the oceanic barium cycle with stable barium isotopes. Earth Planet Sci Lett 434:1–9

    Article  Google Scholar 

  10. Bates SL, Hendry KR, Pryer HV, Kinsley CW, Pyle KM, Woodward EMS, Horner TJ (2017) Barium isotopes reveal role of ocean circulation on barium cycling in the Atlantic. Geochim Cosmochim Acta 204:286–299

    Article  Google Scholar 

  11. Hsieh Y-T, Henderson GM (2017) Barium stable isotopes in the global ocean: tracer of Ba inputs and utilization. Earth Planet Sci Lett 473:269–278

    Article  Google Scholar 

  12. Bridgestock L, Hsieh Y-T, Porcelli D, Homoky WB, Bryan A, Henderson GM (2018) Controls on the barium isotope compositions of marine sediments. Earth Planet Sci Lett 481:101–110

    Article  Google Scholar 

  13. Hendry KR, Pyle KM, Butler GB, Cooper A, Fransson A, Chierici M, Leng MJ, Meyer A, Dodd PA (2018) Spatiotemporal variability of barium in Arctic sea-ice and seawater. J Geophys Res Oceans 123:3507–3522

    Article  Google Scholar 

  14. Crockford PW, Wing BA, Paytan A, Hodgskiss MSW, Mayfield KK, Hayles JA, Middleton JE, Ahm A-SC, Johnston DT, Caxito F, Uhlein G, Halverson GP, Eickmann B, Torres M, Horner TJ (2019) Barium-isotopic constraints on the origin of post-Marinoan barites. Earth Planet Sci Lett 519:234–244

    Article  Google Scholar 

  15. Gou L-F, Jin Z, Galy A, Gong Y-Z, Nan X-Y, Jin C, Wang X-D, Bouchez J, Cai H-M, Chen J-B, Yu H-M, Huang F (2020) Seasonal riverine barium isotopic variation in the middle Yellow River: Sources and fractionation. Earth Planet Sci Lett 531:115990

    Article  Google Scholar 

  16. Gong Y, Zeng Z, Cheng W, Lu Y, Zhang L, Yu H, Huang F (2020) Barium isotopic fractionation during strong weathering of basalt in a tropical climate. Environ Int 143:105896

    Article  Google Scholar 

  17. Renock D, Landis JD, Sharma M (2016) Reductive weathering of black shale and release of barium during hydraulic fracturing. Appl Geochem 65:73–86

    Article  Google Scholar 

  18. Paukert Vankeuren AN, Hakala JA, Jarvis K, Moore JE (2017) Mineral reactions in shale gas reservoirs: barite scale formation from reusing produced water as hydraulic fracturing fluid. Environ Sci Technol 51:9391–9402

    Article  Google Scholar 

  19. Jew AD, Li Q, Cercone D, Maher K, Brown GE, Jr., Bargar JR (2018) Barium sources in hydraulic fracturing systems and chemical controls on its release into solution. Unconventional Resources Technology Conference (URTeC), paper 2899671.

  20. Tieman ZG, Stewart BW, Capo RC, Phan TT, Lopano CL, Hakala JA (2020) Barium isotopes track the source of dissolved solids in produced water from the unconventional Marcellus Shale gas play. Environ Sci Technol 54:4275–4285

    Article  Google Scholar 

  21. Andreasen R, Sharma M (2007) Mixing and homogenization in the early solar sytem: clues from Sr, Ba, Nd, and Sm isotopes in meteorites. Astrophys J 665:874–883

    Article  Google Scholar 

  22. Barling J, Weis D (2008) Influence of non-spectral matrix effects on the accuracy of Pb isotope ratio measurement by MC-ICP-MS: implications for the external normalization method of instrumental mass bias correction. J Anal Atom Spectrom 23:1017

    Article  Google Scholar 

  23. Nan X, Wu F, Zhang Z, Hou Z, Huang F, Yu H (2015) High-precision barium isotope measurements by MC-ICP-MS. J Anal Atom Spectrom 30:2307–2315

    Article  Google Scholar 

  24. van Zuilen K, Nägler TF, Bullen TD (2016) Barium isotopic compositions of geological reference materials. Geostand Geoanalyt Res 40:543–558

    Article  Google Scholar 

  25. Zeng Z, Li X, Liu Y, Huang F, Yu HM (2019) High-precision barium isotope measurements of carbonates by MC-ICP-MS. Geostand Geoanalyt Res 43:291–300

    Article  Google Scholar 

  26. Lin Y-B, Wei H-Z, Jiang S-Y, Hohl S, Lei H-L, Liu X, Dong G (2020) Accurate determination of barium isotopic compositions in sequentially leached phases from carbonates by double spike-thermal ionization mass spectrometry (DS-TIMS). Anal Chem 92:2417–2424

    Article  Google Scholar 

  27. Yu Y, Siebert C, Fietzke J, Goepfert T, Hathorne E, Cao Z, Frank M (2020) The impact of MC-ICP-MS plasma conditions on the accuracy and precision of stable isotope measurements evaluated for barium isotopes. Chem Geol 549:119697

    Article  Google Scholar 

  28. Tian L-L, Zeng Z, Nan X-Y, Yu H-M, Huang F (2019) Determining Ba isotopes of barite using the Na2CO3 exchange reaction and double-spike method by MC-ICP-MS. J Anal Atom Spectrom 34:1459–1467

    Article  Google Scholar 

  29. Tian L-L, Gong Y-Z, Wei W, Kang J-T, Yu H-M, Huang F (2020) Rapid determination of Ba isotope compositions for barites using a H2O-extraction method and MC-ICP-MS. J Anal Atom Spectrom 35:1566–1573

    Article  Google Scholar 

  30. Miyazaki T, Kimura J-I, Chang Q (2014) Analysis of stable isotope ratios of Ba by double-spike standard-sample bracketing using multiple-collector inductively coupled plasma mass spectrometry. J Anal Atom Spectrom 29:483

    Article  Google Scholar 

  31. Nan X-Y, Yu H-M, Rudnick RL, Gaschnig RM, Xu J, Li W-Y, Zhang Q, Jin Z-D, Li X-H, Huang F (2018) Barium isotopic composition of the upper continental crust. Geochim Cosmochim Acta 233:33–49

    Article  Google Scholar 

  32. Foster DA, Staubwasser M, Henderson GM (2004) 226Ra and Ba concentrations in the Ross Sea measured with multicollector ICP mass spectrometry. Mar Chem 87:59–71

    Article  Google Scholar 

  33. Merschel G, Bau M (2015) Rare earth elements in the aragonitic shell of freshwater mussel Corbicula fluminea and the bioavailability of anthropogenic lanthanum, samarium and gadolinium in river water. Sci Total Environ 533:91–101

    Article  Google Scholar 

  34. Barbot E, Vidic NS, Gregory KB, Vidic RD (2013) Spatial and temporal correlation of water quality parameters of produced waters from Devonian-age shale following hydraulic fracturing. Environ Sci Technol 47:2562–2569

    Article  Google Scholar 

  35. Ouyang B, Akob DM, Dunlap D, Renock D (2017) Microbially mediated barite dissolution in anoxic brines. Appl Geochem 76:51–59

    Article  Google Scholar 

  36. Cole RB, Stewart BW (2009) Continental margin volcanism at sites of spreading ridge subduction: examples from southern Alaska and western California. Tectonophysics 464:118–136

    Article  Google Scholar 

  37. Raczek I, Stoll B, Hofmann AW, Jochum KP (2001) High-precision trace element data for the USGS reference materials BCR-1, BCR-2, BHVO-1, BHVO-2, AGV-1, AGV-2, DTS-1, DTS-2, GSP-1 and GSP-2 by ID-TIMS and MIC-SSMS. Geostand Newslett 25:77–86

    Article  Google Scholar 

  38. Breit GN, Simmons EC, Goldhaber MB (1985) Dissolution of barite for the analysis os strontium isotopes and other chemical and isotopic variations using aqueous sodium carbonate. Chem Geol (Isotope Geosc Sect) 52:333–336

    Article  Google Scholar 

  39. Horner TJ, Pryer HV, Nielsen SG, Crockford PW, Gauglitz JM, Wing BA, Ricketts RD (2017) Pelagic barite precipitation at micromolar ambient sulfate. Nat Commun 8:1342

    Article  Google Scholar 

  40. Böttcher ME, Geprägs P, Neubert N, von Allmen K, Pretet C, Samankassou E, Nägler TF (2012) Barium isotope fractionation during experimental formation of the double carbonate BaMn[CO3]2 at ambient temperature. Isotopes Env Health Stud 48:457–463

    Article  Google Scholar 

  41. Mavromatis V, van Zuilen K, Purgstaller B, Baldermann A, Nägler TF, Dietzel M (2016) Barium isotope fractionation during witherite (BaCO3) dissolution, precipitation and at equilibrium. Geochim Cosmochim Acta 190:72–84

    Article  Google Scholar 

  42. Pretet C, van Zuilen K, Nägler TF, Reynaud S, Böttcher ME, Samankassou E (2016) Constraints on barium isotope fractionation during aragonite precipitation by corals. Deposit Rec 1:118–129

    Article  Google Scholar 

  43. van Zuilen K, Müller T, Nägler TF, Dietzel M, Küsters T (2016) Experimental determination of barium isotope fractionation during diffusion and adsorption processes at low temperatures. Geochim Cosmochim Acta 186:226–241

    Article  Google Scholar 

  44. Rudge JF, Reynolds BC, Bourdon B (2009) The double spike toolbox. Chem Geol 265:420–431

    Article  Google Scholar 

  45. Wasserburg GJ, Jacobsen SB, DePaolo DJ, McCulloch MT, Wen T (1981) Precise determination of Sm/Nd ratios, Sm and Nd isotopic abundances in standard solutions. Geochim Cosmochim Acta 45:2311–2323

    Article  Google Scholar 

  46. An YJ, Li X, Zhang ZF (2020) Barium isotopic compositions in thirty-four geological reference materials analysed by MC-ICP-MS. Geostand Geoanalyt Res 44:183–199

    Article  Google Scholar 

  47. McLennan SM (2001) Relationship between the trace element composition of sedimentary rocks and upper continental crust. Geochem Geophy Geosyst.

    Article  Google Scholar 

  48. Gale A, Dalton CA, Langmuir CH, Su Y, Schilling J-G (2013) The mean composition of ocean ridge basalts. Geochem Geophy Geosyst 14:489–518

    Article  Google Scholar 

  49. Niu Y (2004) Bulk-rock major and trace element compositions of abyssal peridotites: implications for mantle melting, melt extraction and post-melting processes beneath mid-ocean ridges. J Petrol 45:2423–2458

    Article  Google Scholar 

  50. Gromet LP, Dymek RF, Haskin LA, Korotev RL (1984) The “North American shale composite”: its compilation, major and trace element characteristics. Geochim Cosmochim Acta 48:2469–2482

    Article  Google Scholar 

  51. Li F, Webb GE, Algeo TJ, Kershaw S, Lu C, Oehlert AM, Gong Q, Pourmand A, Tan X (2019) Modern carbonate ooids preserve ambient aqueous REE signatures. Chem Geol 509:163–177

    Article  Google Scholar 

  52. Geyman BM, Ptacek JL, LaVigne M, Horner TJ (2019) Barium in deep-sea bamboo corals: Phase associations, barium stable isotopes, & prospects for paleoceanography. Earth Planet Sci Lett 525:115751

    Article  Google Scholar 

  53. Guichard F, Church TM, Treuil M, Jaffrezic H (1979) Rare earths in barites: distribution and effects on aqueous partitioning. Geochim Cosmochim Acta 43:983–997

    Article  Google Scholar 

  54. Nozaki Y (1997) A fresh look at element distribution in the North Pacific Ocean. Eos 78:221–223

    Article  Google Scholar 

  55. Gaillardet J, Viers J, Dupré B (2005) Trace elements in river waters. In: Drever JI (ed) Surface and Ground Water, Weathering, and Soils: Treatise on Geochemistry, vol 5. Elsevier, pp 225–272

    Google Scholar 

Download references


We thank W. Heck for assistance with column calibrations. We appreciate comments from two anonymous reviewers.


This work was supported by the U.S. Department of Energy, Office of Fossil Energy Advanced Carbon Storage and Onshore Unconventional portfolios, and by an Oak Ridge Institute for Science and Education research internship to RMM.

This work was funded by the Department of Energy, National Energy Technology Laboratory, an agency of the United States Government, through a support contract with Leidos Research Support Team (LRST). Neither the United States Government nor any agency thereof, nor any of their employees, nor LRST, nor any of their employees, makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.

Author information

Authors and Affiliations



RMM carried out the dissolution and Ba separation experiments, wrote the first draft, and edited subsequent drafts; BWS and RCC helped plan the experiments, carried out Ba isotope analyses, assisted with data interpretation, and edited the manuscript; RLT analyzed the column elution concentrations and edited the manuscript; JAH helped plan the experiments and edited the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to B. W. Stewart.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Table S1. Information on samples used for column calibration experiments. Table S2. Thermo Element® ICP-MS signal intensity (in counts per second) for a column calibration using 2.0 N HCl for the entire elution of a dissolved silicate volcanic rock. Table S3. Thermo Element® ICP-MS signal intensity (in counts per second) for a column calibration using 2.5 N HCl for the entire elution of a dissolved silicate volcanic rock. Table S4. Thermo Element® ICP-MS signal intensities (in counts per second) for a column calibration using 2.5 N HCl for the first 20 mL, and 3.0 N HNO3 for the remaining elution of a produced water sample with 0.5 μg of La and Ce added. Table S5. Thermo Element® ICP-MS signal intensities (in counts per second) for a column calibration using 2.5 N HCl for the first 20 mL, and 2.0 N HNO3 for the remaining elution of dissolved freshwater mussel shell calcium carbonate. Figure S1. Examples of column matrix separation for major elements, Ba, and isobaric interferents La and Ce for a volcanic silicate sample using the BIO-RAD Poly-Prep® gravity flow ion exchange columns described in the manuscript. Figure S2. Examples of column matrix separation of Ba from isobaric interferents La and Ce for a produced water sample (a) and calcium carbonate sample (b) using the BIO-RAD Poly-Prep® gravity flow ion exchange columns described in the manuscript.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Matecha, R.M., Capo, R.C., Stewart, B.W. et al. A single column separation method for barium isotope analysis of geologic and hydrologic materials with complex matrices. Geochem Trans 22, 4 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: