Seismic detectability of carbonates in the deep Earth: A nuclear inelastic scattering study

Carbonates play an important role in the transport and storage of carbon in the Earth’s mantle. However, the abundance of carbon and carbonates in subduction zones is still an unknown quantity. To determine the most abundant accessory phases and how they influence the dynamical processes that operate within the Earth, investigations on the vibrational, elastic, and thermodynamic properties of these phases are crucial for interpreting seismological observations. Recently, the nuclear inelastic scattering (NIS) method has proved to be a useful tool to access information on the lattice dynamics, as well as to determine Debye sound velocities of Fe-bearing materials. Here we derive the acoustic velocities from two carbonate compositions in the FeCO3-MgCO3 binary system up to ~70 GPa using the NIS method. We conclude that more Mg-rich samples, in this case (Fe0.26Mg0.74)CO3, have ~19% higher sound velocities than the pure end-member Fe composition. In addition, we observed a significant velocity increase after the Fe2+ spin transition was complete. After laser heating of FeCO3 at lower mantle conditions, we observed a dramatic velocity drop, which is probably associated with thermal decomposition to another phase. Parallel to our NIS experiments, we conducted a single-crystal X-ray diffraction (SCXRD) study to derive the equation of states of FeCO3 and (Fe0.26Mg0.74)CO3. The combined information from NIS (i.e., Debye velocities) and SCXRD (i.e., densities and bulk moduli) experiments enabled us to derive the primary and shear wave velocities of our samples. Our results are consistent with results obtained by other methods in previous studies, including Brillouin spectroscopy, inelastic X-ray scattering, and DFT calculations, supporting NIS as a reliable alternative method for studying the elastic properties of Fe-bearing systems at high pressures and temperatures. Finally, we discuss the seismic detectability of carbonates. We determine that nearly 22 wt% CO2 must be present in the subduction slab to detect a 1% shear wave velocity decrease compared to non-carbonated lithologies at the transition zone to lower mantle boundary depths.