
==== Front
Sci Adv
Sci Adv
sciadv
advances
Science Advances
2375-2548
American Association for the Advancement of Science

ado8550
10.1126/sciadv.ado8550
Research Article
Earth, Environmental, Ecological, and Space Sciences
SciAdv r-articles
Geophysics
Geophysics
The ultralow viscosity of volatile-rich kimberlite magma: Implications for the water content of primitive kimberlite melts
The ultralow viscosity of volatile-rich kimberlite magma
https://orcid.org/0000-0003-1380-3569
Hao Ming Conceptualization Formal analysis Investigation Methodology Project administration Visualization Writing - original draft Writing - review & editing 1 2 *
Zhou Wen-Yi Investigation Methodology 1 3
https://orcid.org/0000-0002-1489-2156
Hrubiak Rostislav Data curation Investigation Methodology Software 4
https://orcid.org/0000-0001-7258-8286
Kenney-Benson Curtis Resources 4
Kavanagh Janine L. Conceptualization Validation Writing - review & editing 5
Davis William Formal analysis 6
https://orcid.org/0000-0002-1936-7780
Zhang Jin S. Conceptualization Funding acquisition Investigation Methodology Project administration Resources Supervision Validation Visualization Writing - review & editing 1 3 *
1 Department of Earth and Planetary Sciences, University of New Mexico, Albuquerque, NM, USA.
2 Earth and Planets Laboratory, Carnegie Institution for Science, Washington, DC, USA.
3 Department of Geology and Geophysics, Texas A&M University, College Station, TX, USA.
4 HPCAT, Advanced Photon Source, Argonne National Laboratory, Argonne, IL, USA.
5 School of Environmental Sciences, University of Liverpool, Liverpool L69 3GP, UK.
6 Cecil H. and Ida M. Green Institute of Geophysics and Planetary Physics, Scripps Institution of Oceanography, University of California, San Diego, La Jolla, CA 92037, USA.
* Corresponding author. Email: unmminghao@gmail.com (M.H.); jinzhang@tamu.edu (J.S.Z.)
06 9 2024
06 9 2024
10 36 eado855023 2 2024
01 8 2024
Copyright © 2024 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The eruption of deeply sourced kimberlite magma offers the fastest route to bring deep-seated volatiles back to the Earth’s surface. However, the viscosity of kimberlite magma, a factor governing its migration and eruption dynamics within Earth, remains poorly constrained. We conducted synchrotron in situ falling sphere viscometry experiments to examine kimberlite magma with different volatile contents (0 to 5 wt % H2O and 2 to 8 wt % CO2) under high pressure-temperature conditions. The results reveal that the viscosity of volatile-rich kimberlite magma is ~1 to 2 orders lower than that of mid-ocean ridge basalt (MORB) and comparable to the ultramobile pure carbonate melt. Using the measured viscosity values, we simulated the ascent and eruption process of kimberlite magma. We found that a minimum content of ~0.5 wt % water in the primitive magma is necessary to allow the ultrafast eruption process of kimberlite, thereby enabling the preservation of diamonds and high-pressure mineral inclusions transported by the magma.

Kimberlite magmas with over 0.5 wt % water ascend rapidly due to their extremely low viscosity.

http://dx.doi.org/10.13039/100007179 University of New Mexico WeR1 Faculty Success Fund http://dx.doi.org/10.13039/100014013 UK Research and Innovation MR/S035141/1 Texas A&M University Start-up Fund
==== Body
pmcINTRODUCTION

Kimberlite magma, known for hosting diamonds, holds substantial importance among terrestrial magmas owing to its volatile-rich composition and deep-seated origin (1–4) (e.g., >200 km). They are excellent geochemical probes of the deep Earth (5, 6). For example, the Mg/Si ratios of kimberlite magma and the mantle xenoliths enclosed within it can help to constrain the local mantle chemical composition (5). Typically, kimberlite magma originates from partial melting of CO2-bearing peridotite (7–9) or subducted carbonated slab crust (10, 11) at depths greater than ~200 km, possibly within the mantle transition zone (2, 11, 12). Therefore, kimberlite serves as an important portal into understanding the deep volatile cycles inside Earth (6). Nevertheless, the water content of primitive kimberlite magma remains controversial, highlighted by recent conflicting studies (5, 13). The water in primitive kimberlite magma not only influences its transport properties but also affects the elemental partitioning (14). Moreover, if the source mantle of kimberlite is wet, then it could serve as an indicator for ancient subduction zones (15). Therefore, it is critical to constrain the water content in the primitive kimberlite melts.

Determining the water content of primitive kimberlite melts using the collected field samples is difficult, partially due to post-emplacement alteration processes such as serpentinization, potential surface water additions at shallower depths, and interactions between the kimberlite magma and host rocks during dike ascent (16–18). For example, Russell et al. (16, 17) suggested that during ascent, primitive kimberlite magma with a carbonate melt phase entraps mantle xenoliths and assimilates pyroxenes within these xenoliths, which leans to the increase in SiO2 content of the melt and the release of CO2, forming the observed typical kimberlite at the surface. The inferred water content of primitive kimberlite from field samples varies from <1 wt % to 6 to 9 wt % (19–21). Thus, exploring aspects beyond traditional rock sampling could offer valuable insights into the elusive water content of primitive kimberlite melts.

Viscosity is an important physical property influencing the migration and eruption dynamics of magmas. Kimberlite magma is renowned for its remarkable ascent and eruption speed (1, 16, 22). Sparks et al. (22) modeled kimberlite magma eruptions and suggested that a sufficiently low viscosity would enable the extremely rapid eruption process (e.g., ~4 to 20 m/s) under the turbulent flow regime. As an important network modifier, H2O is expected to have strong influence on the viscosity of magmas (23–25). Hence, an exploration into the volatile-dependent viscosity of kimberlite magmas holds the potential to provide additional constraints on the water content of primitive kimberlite melts, offering insights into their rapid eruption processes.

However, the viscosity of volatile-rich kimberlite magma, especially regarding the influence of water content, lacks robust experimental constraints. A recent study using quenched falling sphere experiments found that the viscosity of kimberlite magma is comparable to or even higher than that of MORB under high pressure-temperature conditions and does not depend on the H2O content (26). Such high viscosity values do not support the fast ascent of kimberlite magma proposed in previous studies (1, 22).

To address this knowledge gap, we conducted in situ falling sphere viscometry measurements on kimberlite magma with 2 to 8 wt % CO2 and 0 to 5 wt % H2O up to 5.3 GPa and 2173 K using the Paris-Edinburgh press at Sector 16-BM-B, High-Pressure Collaborative Access Team (HPCAT), Advanced Photon Source (27) and examined the chemical composition of a representative run product that survived after decompression (Materials and Methods, figs. S1 and S3, and table S1). Then, we calculated the viscosity of kimberlite magma under different pressure-temperature conditions using the terminal velocities, which are the calculated maximum sphere falling speeds via analyzing x-ray images from a high-speed camera. Using the new viscosity data in this study, we further model the ascent and eruption processes of kimberlite magma and discuss the possible range of water content in primitive kimberlite melts.

RESULTS

Figure 1 shows all the experimental results. With only 0.5 wt % H2O, the viscosity of kimberlite magma decreases by a factor of ~6, from 0.1 to 0.2 Pa·s to 0.02 to 0.03 Pa·s. With 2 wt % H2O, the viscosity of kimberlite magma decreases by more than an order of magnitude. When the H2O content reaches 5 wt %, the viscosity of kimberlite magma drops to less than 0.01 Pa·s, approaching values akin to that of pure carbonate melts (28) at similar pressure-temperature conditions and on the same order of the viscosity of liquid water at ambient condition. Addition of CO2 could also decrease the viscosity of kimberlite magma, consistent with prior studies reporting ultralow viscosity values for pure carbonate and carbonate-silicate transitional melts (28, 29). However, our results underscore the stronger impact of H2O compared to CO2. For example, with addition of only 0.5 wt % H2O and 2 wt % CO2, the viscosity of kimberlite magma closely resembles that of anhydrous carbonate-silicate transitional melts with ~22.5 wt % CO2 (29).

Fig. 1. The viscosities of kimberlite, MORB, CaCO3, and carbonate-silicate transitional melts under high pressure-temperature conditions.

The numbers before the “W” and “C” represent the wt % of the H2O and CO2 in the melts, respectively. The diamond symbols represent the viscosity values of kimberlite magmas with different CO2 and H2O contents from this study. The gray circles and squares represent the viscosity measurements of MORB (30) and carbonate-silicate transitional melts with ~22.5 wt % CO2 (29), respectively. The blue dash-dotted line represents the viscosity values of kimberlite magma from the quenched experiments (26). The green dashed line represents the viscosity of pure carbonate melts (28).

Persikov et al. (26) is the only existing experimental study regarding the viscosity of hydrous kimberlite magma at high pressure-temperature conditions. The reported viscosity values in Persikov et al. (26) for compositions similar to those in our study are ~0.1 to 0.7 Pa·s which are comparable to the measured viscosity of MORB at similar pressure-temperature conditions using in situ falling sphere viscometry (30). The reported viscosity of hydrous kimberlite magma in Persikov et al. (26) is even higher than MORB at pressures higher than ~3.5 GPa. However, our experimental data, which also used in situ falling sphere viscometry, reveals that the viscosity values of kimberlite magma with different volatile contents are always lower than MORB across the investigated pressure-temperature range (~2 to 5 GPa and 2073 to 2173 K). The notable higher viscosity values in Persikov et al. (26) may result from the lower velocities of falling spheres using the quenched method. In the quenched experiments, the sphere falling process was not directly recorded, leading to the determination of average sphere falling speeds instead of terminal velocity (Fig. 2). Consequently, the calculated viscosity values using the average sphere falling speeds in their quenched experiments likely represent upper limits rather than the true viscosity values of the liquids. In contrast, as shown in Fig. 2, we can calculate the velocity change over time throughout the entire sphere falling process using x-ray radiography images, which enables us to accurately determine the terminal velocities (31).

Fig. 2. Movements of the rhenium spheres and the velocity evolution during the falling sphere viscosity measurements.

The numbers before the W and C represent the wt % of the H2O and CO2 in the melts, respectively. All spheres reached terminal velocities during the experiments. The orange bars in the images represent 0.2 mm.

DISCUSSION

To investigate the influence of water on the ascent and eruption process of kimberlite magma, we calculated the chemical composition, xenoliths content, and bubble content of kimberlite magma change as a function of depth, drawing from established models by Sparks et al. (22) and Rusell et al. (16, 17) (Fig. 3). As shown in Fig. 3, the original kimberlite magma is carbonate melt at depth greater than 150 km. At the initial stage of ascent, the original carbonate melt gradually assimilates mantle xenoliths with limited chemical reactions, and thus only the proportion of suspended crystals increases (17). At ~100- to 150-km depth, the pyroxenes, especially orthopyroxene, within the captured mantle xenoliths react with carbonate melts (16), releasing CO2. This released CO2 undergoes a bubbling process, potentially separating from the remaining melts due to the low density of CO2 fluids and the reduced viscosity of the melts. Consequently, a volatile-rich buoyant tip emerges at this stage. Despite this, the overall crystal content remains relatively stable due to the countereffects of xenoliths assimilation and simultaneous chemical reaction. As xenoliths dissolve in the melts, the CO2 content in the melts decreases, while the SiO2 content increases. At depth less than 100 km, the ongoing release of CO2 results primarily from pressure exsolution instead of chemical reactions. In addition, olivine crystals start to crystallize within this depth range. The detailed calculation of the chemical composition evolution for the primitive kimberlite at depths is shown in Materials and Method.

Fig. 3. Kimberlite magma ascent and evolution process.

(A) Schematic illustration of the kimberlite magma ascent and eruption process. (B) The calculated compositional changes of kimberlite magma with depth. The CO2 and SiO2 contents represent the chemical variations of the melts. The crystal content shows the total volume fractions in the magma. The bubble content includes both the bubble within the magma and accumulated at the tip of the dyke.

The formation of crystals and bubbles resulting from assimilation and exsolution can substantially influence the effective viscosity of magmas (text S1) (32). As shown in Fig. 3, the released bubbles could either suspend within the melts or accumulate forming a volatile-rich tip. Since only the bubbles suspended in the melts affect the effective viscosity, we examined three different scenarios to calculate the viscosity of primitive kimberlite melts during ascent and eruption: (i) CO2 bubbles completely separate from the melt, exerting no influence on the effective viscosity; (ii) all CO2 bubbles remain suspended in the melt without any segregation, thereby excluding the formation of a volatile-rich tip; (iii) only a small amount of newly released CO2 bubbles (~5 vol %) remains suspended in the melts, while the majority segregates from the melt, forming a volatile-rich tip. Then, using the rheology of the three-phase suspension model with a high capillary number (text S1), we calculated the effects of crystals and bubbles on the viscosity of primitive kimberlite melts (32). To cover the large composition range of the kimberlite magma, we have included both the experimental results from this study and data from Stagno et al. (29) on carbonate-silicate transitional melts to model the effects of CO2 and H2O on the viscosity of kimberlite melts. In addition, we corrected for temperature effects using the predicted temperature profile outlined by Kavanagh and Sparks (33). The equations and detailed calculations are shown in text S1.

Turbulent flow regime has been suggested for kimberlite magma eruption based on petrologic textural evidence of field kimberlite samples (e.g., the milling of xenoliths and xenocrysts in kimberlite deposits) (34, 35). The ultrafast eruption of kimberlite magma, as indicated by various eruption models, also requires adoption of the turbulent flow regime (17, 22). The critical Reynolds number (Re) to form turbulent kimberlite magma flows during eruption is ~1000, corresponding to an effective viscosity of 1.24 Pa·s (22). As shown in Fig. 4, we assumed three different water contents (anhydrous, 0.5 and 2 wt %) in primitive kimberlite magma and fixed the water content at different depths to estimate the effective viscosity of kimberlite melts during the ascent and eruption. The calculated viscosity profiles at three hydration levels do not necessarily reflect the actual water content variations during its ascent but provided useful insights to the range of viscosity values we would expect during the eruption (33). At depths greater than ~150 km, because of the low viscosity of carbonate melts (28), the primitive kimberlite magmas with mantle xenoliths are all in turbulent flow regime regardless of the water content. At depths between 50 and 150 km, the exsolution of CO2 in primitive kimberlite melts gradually alters the chemical composition from carbonate to typical ultramafic silicate kimberlite melts (17). The effective viscosity gradually increases, but the hydrous kimberlite melts are still well within the turbulent flow regime. However, the effective viscosity of anhydrous kimberlite magmas is high enough causing the corresponding Re number to decrease below 1000, placing the flow regime between turbulent and steady flow. At depth less than 50 km, the continuous exsolution of CO2 and olivine crystallization further elevates the effective viscosity of kimberlite melts. With 0.5 wt % water in the melts, the Re number is still close to 1000 even with complete CO2 bubble segregation, while the Re number for anhydrous kimberlite is already approaching or even falling below 10 with partially segregated bubbles.

Fig. 4. Modeled effective viscosity of primitive kimberlite melts changes as a function of depth.

The bubbles start to show up at depth of ~150 km where orthopyroxene (opx) starts to react with carbonate melts. A Re number larger than 1000 represents turbulent flow, and a Re number smaller than 10 represents steady flow. The horizontal axe is in log scale.

Considering the natural eruption process of kimberlite magma, a volatile-rich tip is expected (1, 17). Thus, partial or complete bubble segregation is a more realistic scenario. In this case, anhydrous primitive kimberlite melts exhibit excessive viscosity at the top of upper mantle, hindering turbulent flow regimes. However, the addition of only 0.5 wt % water is enough to maintain turbulent flow in the magma. Therefore, a minimal quantity of water is essential for the fast ascent and eruption of primitive kimberlite magma in the upper mantle.

To summarize, we used falling sphere viscometry to measure the viscosity of volatile-rich kimberlite magma, revealing viscosity values 1 to 2 orders of magnitude lower than previous studies (26). On the basis of these viscosity measurements, we quantified the effect of water on the ascent and eruption dynamics of kimberlite magma. Our findings indicate that a relatively small quantity of water (~0.5 wt %) plays a crucial role in enabling the rapid ascent and eruption of primitive kimberlite magma, which could help to preserve the high-pressure mineral phases (e.g., ringwoodite) in the super-deep diamonds (36, 37) as well as the major elements, trace elements, and isotope signatures of the local mantle sampled by the xenoliths (38).

MATERIALS AND METHODS

Starting materials

It is difficult to determine a representative major element composition of kimberlite melt, since its chemical composition likely changes during its ascent process (16–18, 22), from carbonate melts at depths greater than 150 km to carbonate-silicate transitional melts and lastly the observed ultramafic carbon-bearing kimberlite melts (16, 17). Since previous studies has measured the viscosity of carbonate melts and carbonate-silicate transitional melts using in situ falling sphere viscometry (28, 29), this study focuses on the ultramafic kimberlite melts with relatively low CO2 contents (2 to 8 wt % CO2) at different hydration levels (0 to 5 wt % H2O). The major element composition of the starting material is the average of the kimberlite field samples from nine different places (table S1). Oxides, carbonates, and hydroxides with purity > 99% were ground, mixed, and pressed as starting materials (19).

Falling sphere viscometry

We conducted high pressure-temperature in situ falling sphere experiments at Sector 16-BM-B, HPCAT, Advanced Photon Source using Paris-Edinburgh cell in graphite capsules (27) with 1.5 mm in diameter and 2.0 mm in height (fig. S1). The temperatures were calculated from the power curve, which was calibrated according to the temperature at the center of the sample (27). The temperature uncertainty of each experimental run is estimated to be 100 K (27). The temperature gradient across the entire sample along the horizontal and vertical directions is <40 K at 2000 K (39), and thus the effect of temperature gradient on the obtained terminal velocity is minimal. The pressures were determined using the equation of state of MgO, with the pressure difference between MgO sleeve and the actual sample chamber taken into account (27).

For all experimental runs, once the target pressure was reached, we firstly increased temperature gradually from 300 to 1073 K, which is well below the solidus of the sample, in roughly 10 min. Then, we increased temperature abruptly from 1073 to 2073 K or 2173 K by adjusting the input power to the targeted value. The fast and complete melting was confirmed by immediate and straight descent of rhenium spheres after the power adjustment. Because of the ultrafast nature of these falling sphere experiments, volatile loss was negligible, and chemical contamination from graphite capsule was limited to the edge of the sample, as evidenced by the chemical analysis on the run product (40).

Parallel beam polychromatic x-rays were used for imaging the sample chamber (Fig. 2). A high-speed camera (Photron FASTCAM SA3) with 1000 frames per second was used to capture the sample chamber images during the sphere falling process. A tungsten carbide sphere with a 497-μm diameter was used to calibrate the size of the pixels of the camera. Rhenium spheres with diameters of ~80 to 190 μm were used to do falling sphere viscometry measurements. The solubility of rhenium is <10 parts per million in silicate melts at experimental conditions (41), and thus the sample is unlikely to be contaminated by rhenium spheres. The position of the rhenium spheres recorded by the high-speed camera was analyzed using ImageJ software to calculate sphere falling speeds. The terminal velocities (ν), which are the maximum sphere falling speeds, are calculated from the velocity-time curves (Fig. 2). The viscosity (η) was calculated with the Stokes equation with correction factors outlined in Kono et al. (27)η=gds2ρs−ρlF18νE(1)

F=1−2.104dsdl+2.09dsdl3−0.95dsdl5(2)

E=1+9ds16Z+9ds16Z2(3)

where g is gravitational acceleration, and d and ρ are diameters and densities of the sphere (s) and liquid (l), respectively. Z is the height of the sample capsule. F is the Faxén correction factor (42), and E is the end correction factor (43). Equation 1 is applicable for ds/dl ≪ 1 and Re ≪ 1. The first condition was satisfied by the experimental design, whereas post hoc validations of Re were attained by finding self-consistent solutions of the relationRe=ρldsνη(4)

giving Re < 0.2. Uncertainties in viscosity estimates are largely influenced by the uncertainty in the terminal velocity measurements (44); alternative edge effect corrections gave viscosity estimates within uncertainties (45).

Chemical analysis on quenched sample

Although most samples did not survive during the fast quench and decompression process, we were able to identify a well-preserved sample (#3-10 in table S2) which went through complete melting at 3.8 GPa 2173 K and contained ~2 wt % water and ~2 wt % of CO2. We then performed electron microprobe analysis (EPMA), energy-dispersive spectroscopy (EDS) mapping, and Fourier transform infrared spectra (FTIR) mapping to examine the chemical composition of the quenched sample as well as its homogeneity. The FTIR measurements were conducted using a Bruker HYPERION II FTIR microscope with a liquid N2 cooled mid-band mercury cadmium telluride detector at the Texas A&M University (46). The EPMA experiments were completed using a JEOL 8530F electron microprobe with an accelerating voltage of 15 kV and beam current of 30 nA at the Earth and Planets Laboratory, Carnegie Institute of Science. Wavelength dispersive spectroscopy was performed to quantify the chemical compositions with focused (~1 μm) and defocused beam (~5 μm) on random spots (table S1), and EDS mapping was performed to examine the chemical homogeneity of the entire sample (fig. S3). Experimental results are included in table S1 and fig. S3.

The composition evolution of kimberlite melts at depth

According to previous experiments and models (16, 17), the primary kimberlite melts at depth greater than 200 km are carbonate melts with ~10 wt % SiO2. During its ascent, at depth greater than 100 km, the melts gradually assimilate ~51 vol % orthopyroxene-rich mantle xenoliths (47, 48) to reach a typical SiO2 content of ~30 wt %. At 150- to 100-km depth, the orthopyroxenes in xenoliths react with the melt, gradually releasing CO2 and the CO2 content in the melts decreases to ~10 wt %. At depth less than 100 km, the CO2 content continues to decrease, reaching a level of around 5 wt % due to the pressure-induced exsolution. At the same time, ~5 wt % olivine crystallizes from the melts. The composition evolution is displayed in Fig. 3.

Calculation of Re number of kimberlite melts

The Re number for kimberlite melts in dyke can be calculated using (22)Re=2ΔPρw33hμ2(5)

where ΔP is the overpressure, w is the half dyke width, u is the velocity, ρ is the magma density, h is the dyke vertical length, and μ is the viscosity. Using typical values of ΔP = 10 MPa, w = 0.25 m, ρ = 3000 kg/m3, and h = 200 km, the Re number is ~1000 for a viscosity value of 1.24 Pa·s and the Re number is ~10 for a viscosity value of 12.4 Pa·s (22).

Acknowledgments

The experiments were performed at HPCAT (Sector 16), Advanced Photon Source (APS), Argonne National Laboratory. HPCAT operation is supported by DOE-NNSA under award no. DE-NA0001974, with partial instrumentation funding from NSF. The Advanced Photon Source is a US Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Argonne National Laboratory under contract no. DE-AC02-06CH11357.

Funding: This work was supported by WeR1 Faculty Success Fund from the University of New Mexico (J.S.Z.), Start-up fund from Texas A&M University (J.S.Z.), and UKRI Future Leaders Fellowship (MR/S035141/1; J.L.K.).

Author contributions: Conceptualization: M.H. and J.S.Z. Experiments: M.H., W.-Y.Z., R.H., C.K.-B., and J.S.Z. Models: M.H., J.L.K., W.D., and J.S.Z. Funding acquisition: J.S.Z. Writing—original draft: M.H. and J.S.Z. Writing—review and editing: M.H., W.-Y.Z., R.H., C.K.-B., J.K., W.D., and J.S.Z.

Competing interests: The authors declare that they have no competing interests.

Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials.

Supplementary Materials

This PDF file includes:

Supplementary Text S1

Figs. S1 to S3

Tables S1 and S2

References
==== Refs
REFERENCES AND NOTES

1 L. Wilson, J. W. Head III, An integrated model of kimberlite ascent and eruption. Nature 447 , 53–57 (2007).17476260
2 S. F. Foley, G. M. Yaxley, B. A. Kjarsgaard, Kimberlites from source to surface: Insights from experiments. Elements 15 , 393–398 (2019).
3 K. G. Nickel, D. H. Green, Empirical geothermobarometry for garnet peridotites and implications for the nature of the lithosphere, kimberlites and diamonds. Earth Planet. Sci. Lett. 73 , 158–170 (1985).
4 R. H. Mitchell, Kimberlites: Mineralogy, Geochemistry, and Petrology (Springer Science & Business Media, 2013).
5 D. G. Pearson, J. Woodhead, P. E. Janney, Kimberlites as geochemical probes of Earth’s mantle. Elements 15 , 387–392 (2019).
6 S. Tappe, R. L. Romer, A. Stracke, A. Steenfelt, K. A. Smart, K. Muehlenbachs, T. H. Torsvik, Sources and mobility of carbonate melts beneath cratons, with implications for deep carbon cycling, metasomatism and rift initiation. Earth Planet. Sci. Lett. 466 , 152–167 (2017).
7 R. Dasgupta, A. Mallik, K. Tsuno, A. C. Withers, G. Hirth, M. M. Hirschmann, Carbon-dioxide-rich silicate melt in the Earth’s upper mantle. Nature 493 , 211–215 (2013).23302861
8 M. Becker, A. P. L. Roex, Geochemistry of South African on- and off-craton, group I and group II kimberlites: Petrogenesis and source region evolution. J. Petrol. 47 , 673–703 (2006).
9 S. E. Price, J. K. Russell, M. G. Kopylova, Primitive magma from the Jericho Pipe, N.W.T., Canada: Constraints on primary kimberlite melt chemistry. J. Petrol. 41 , 789–808 (2000).
10 S. Tappe, D. Graham Pearson, B. A. Kjarsgaard, G. Nowell, D. Dowall, Mantle transition zone input to kimberlite magmatism near a subduction zone: Origin of anomalous Nd–Hf isotope systematics at Lac de Gras, Canada. Earth Planet. Sci. Lett. 371–372 , 235–251 (2013).
11 C. Sun, R. Dasgupta, Slab–mantle interaction, carbon transport, and kimberlite generation in the deep upper mantle. Earth Planet. Sci. Lett. 506 , 38–52 (2019).
12 E. S. Kiseeva, K. D. Litasov, G. M. Yaxley, E. Ohtani, V. S. Kamenetsky, Melting and phase relations of carbonated eclogite at 9–21 GPa and the Petrogenesis of Alkali-Rich melts in the deep mantle. J. Petrol. 54 , 1555–1583 (2013).
13 V. S. Kamenetsky, A. V. Golovin, R. Maas, A. Giuliani, M. B. Kamenetsky, Y. Weiss, Towards a new model for kimberlite petrogenesis: Evidence from unaltered kimberlites and mantle minerals. Earth Sci. Rev. 139 , 145–167 (2014).
14 A. G. Sokol, A. N. Kruk, Y. N. Palyanov, The role of water in generation of group II Kimberlite magmas: Constraints from multiple saturation experiments. Am. Mineral. 99 , 2292–2302 (2014).
15 S. Saha, R. Dasgupta, Phase relations of a depleted peridotite fluxed by a CO2-H2O fluid—Implications for the stability of partial melts versus volatile-bearing mineral phases in the cratonic mantle. J. Geophys. Res. Solid Earth 124 , 10089–10106 (2019).
16 J. K. Russell, L. A. Porritt, Y. Lavallée, D. B. Dingwell, Kimberlite ascent by assimilation-fuelled buoyancy. Nature 481 , 352–356 (2012).22258614
17 J. K. Russell, R. S. J. Sparks, J. L. Kavanagh, Kimberlite volcanology: Transport, ascent, and eruption. Elements 15 , 405–410 (2019).
18 R. S. J. Sparks, Kimberlite volcanism. Annu. Rev. Earth Planet. Sci. 41 , 497–528 (2013).
19 V. S. Kamenetsky, M. B. Kamenetsky, Y. Weiss, O. Navon, T. F. Nielsen, T. P. Mernagh, How unique is the Udachnaya-East kimberlite? Comparison with kimberlites from the Slave Craton (Canada) and SW Greenland. Lithos 112 , 334–346 (2009).
20 A. Soltys, A. Giuliani, D. Phillips, A new approach to reconstructing the composition and evolution of kimberlite melts: A case study of the archetypal Bultfontein kimberlite (Kimberley, South Africa). Lithos 304-307 , 1–15 (2018).
21 B. A. Kjarsgaard, D. G. Pearson, S. Tappe, G. M. Nowell, D. P. Dowall, Geochemistry of hypabyssal kimberlites from Lac de Gras, Canada: Comparisons to a global database and applications to the parent magma problem. Lithos 112 , 236–248 (2009).
22 R. S. J. Sparks, L. Baker, R. J. Brown, M. Field, J. Schumacher, G. Stripp, A. Walters, Dynamical constraints on kimberlite volcanism. J. Volcanol. Geotherm. Res. 155 , 18–48 (2006).
23 K.-U. Hess, D. B. Dingwell, Viscosities of hydrous leucogranitic melts: A non-Arrhenian model. Am. Mineral. 81 , 1297–1300 (1996).
24 H. R. Shaw, Viscosities of magmatic silicate liquids: An empirical method of prediction. Am. J. Sci. 272 , 870–893 (1972).
25 D. Giordano, J. K. Russell, D. B. Dingwell, Viscosity of magmatic liquids: A model. Earth Planet. Sci. Lett. 271 , 123–134 (2008).
26 E. S. Persikov, P. G. Bukhtiyarov, A. G. Sokol, Viscosity of hydrous kimberlite and basaltic melts at high pressures. Russ. Geol. Geophys. 58 , 1093–1100 (2017).
27 Y. Kono, C. Park, C. Kenney-Benson, G. Shen, Y. Wang, Toward comprehensive studies of liquids at high pressures and high temperatures: Combined structure, elastic wave velocity, and viscosity measurements in the Paris–Edinburgh cell. Phys. Earth Planet. In. 228 , 269–280 (2014).
28 Y. Kono, C. Kenney-Benson, D. Hummer, H. Ohfuji, C. Park, G. Shen, Y. Wang, A. Kavner, C. E. Manning, Ultralow viscosity of carbonate melts at high pressures. Nat. Commun. 5 , 5091 (2014).25311627
29 V. Stagno, Y. Kono, V. Stopponi, M. Masotta, P. Scarlato, C. E. Manning, “The viscosity of carbonate-silicate transitional melts at Earth’s upper mantle pressures and temperatures, determined by the in situ falling-sphere technique” in Carbon in Earth’s Interior (American Geophysical Union (AGU), 2020), pp. 223–236.
30 T. Sakamaki, A. Suzuki, E. Ohtani, H. Terasaki, S. Urakawa, Y. Katayama, K. Funakoshi, Y. Wang, J. W. Hernlund, M. D. Ballmer, Ponded melt at the boundary between the lithosphere and asthenosphere. Nat. Geosci. 6 , 1041–1044 (2013).
31 C. Liebske, B. Schmickler, H. Terasaki, B. Poe, A. Suzuki, K. Funakoshi, R. Ando, D. Rubie, Viscosity of peridotite liquid up to 13 GPa: Implications for magma ocean viscosities. Earth Planet. Sci. Lett. 240 , 589–604 (2005).
32 J. M. Truby, S. P. Mueller, E. W. Llewellin, H. M. Mader, The rheology of three-phase suspensions at low bubble capillary number. Proc. R. Soc. A Math. Phys. Eng. Sci. 471 , 20140557 (2015).
33 J. L. Kavanagh, R. S. J. Sparks, Temperature changes in ascending kimberlite magma. Earth Planet. Sci. Lett. 286 , 404–413 (2009).
34 R. C. Brett, J. K. Russell, G. D. M. Andrews, T. J. Jones, The ascent of kimberlite: Insights from olivine. Earth Planet. Sci. Lett. 424 , 119–131 (2015).
35 T. J. Jones, J. K. Russell, D. Sasse, Modification of mantle cargo by turbulent ascent of kimberlite. Front. Earth Sci. 7 , (2019).
36 D. G. Pearson, F. E. Brenker, F. Nestola, J. McNeill, L. Nasdala, M. T. Hutchison, S. Matveev, K. Mather, G. Silversmit, S. Schmitz, B. Vekemans, L. Vincze, Hydrous mantle transition zone indicated by ringwoodite included within diamond. Nature 507 , 221–224 (2014).24622201
37 M. J. Walter, S. C. Kohn, D. Araujo, G. P. Bulanova, C. B. Smith, E. Gaillou, J. Wang, A. Steele, S. B. Shirey, Deep mantle cycling of oceanic crust: Evidence from diamonds and their mineral inclusions. Science 334 , 54–57 (2011).21921159
38 A. H. Peslier, A. B. Woodland, J. A. Wolff, Fast kimberlite ascent rates estimated from hydrogen diffusion profiles in xenolithic mantle olivines from Southern Africa. Geochim. Cosmochim. Acta 72 , 2711–2722 (2008).
39 A. F. Riva, A. D. Rosa, C. Clavel, D. Sifre, O. Mathon, G. Garbarino, S. Pascarelli, Heat distribution in paris–edinburgh press assemblies through finite element simulations. High Pressure Res. 38 , 303–324 (2018).
40 W.-Y. Zhou, M. Hao, R. Hrubiak, C. Kenney-Benson, J. S. Zhang, Migration and accumulation of hydrous mantle incipient melt in the Earth’s asthenosphere: Constraints from in-situ falling sphere viscometry measurements. Earth Planet. Sci. Lett. 641 , 118833 (2024).
41 N. R. Bennett, J. M. Brenan, Controls on the solubility of rhenium in silicate melt: Implications for the osmium isotopic composition of earth’s mantle. Earth Planet. Sci. Lett. 361 , 320–332 (2013).
42 H. Faxén, Der Widerstand gegen die Bewegung einer starren Kugel in einer zähen Flüssigkeit, die zwischen zwei parallelen ebenen Wänden eingeschlossenist. Ann. Phys. 373 , 89–119 (1922).
43 A. D. Maude, End effects in a falling-sphere viscometer. Br. J. Appl. Phys. 12 , 293–295 (1961).
44 Y. Kono, C. Kenney-Benson, Y. Shibazaki, C. Park, G. Shen, Y. Wang, High-pressure viscosity of liquid Fe and FeS revisited by falling sphere viscometry using ultrafast X-ray imaging. Phys. Earth Planet. In. 241 , 57–64 (2015).
45 M. Brizard, M. Megharfi, E. Mahé, C. Verdier, Design of a high precision falling-ball viscometer. Rev. Sci. Instrum. 76 , 025109 (2005).
46 M. Mercier, A. D. Muro, N. Métrich, D. Giordano, O. Belhadj, C. W. Mandeville, Spectroscopic analysis (ftir, raman) of water in mafic and intermediate glasses and glass inclusions. Geochim. Cosmochim. Acta 74 , 5641–5656 (2010).
47 M. Gao, C. Wang, H. Xu, J. Zhang, Experiment constraints on orthopyroxene enrichment in the Kaapvaal Craton lithospheric mantle. Chem. Geol. 562 , 120046 (2021).
48 F. R. Boyd, N. P. Pokhilenko, D. G. Pearson, S. A. Mertzman, N. V. Sobolev, L. W. Finger, Composition of the Siberian cratonic mantle: Evidence from Udachnaya peridotite xenoliths. Contrib. Mineral. Petrol. 128 , 228–246 (1997).
49 E. S. Persikov, P. G. Bukhtiyarov, A. G. Sokol, Viscosity of Haplokimberlitic and basaltic melts at high pressures: Experimental and theoretical studies. Chem. Geol. 497 , 54–63 (2018).
50 E. J. Holden, S. Moss, J. K. Russell, M. C. Dentith, An image analysis method to determine crystal size distributions of Olivine in Kimberlite. Comput. Geosci. 13 , 255–268 (2009).
51 T. F. D. Nielsen, K. K. Sand, The Majuagaa Kimberlite Dike, Maniitsoq Region, West Greenland: Constraints on an Mg-rich silicocarbonatitic melt composition from groundmass mineralogy and bulk compositions. Can. Mineral. 46 , 1043–1061 (2008).
