Core-Mantle Interface - Thermodynamics and chemistry

A 5-year multidisciplinary project funded by the European Research Council under the European Union's Seventh Framework Programme (FP7/2007-2013) / ERC Grant agreement 240473.

CoMITAC and related software

MSAT — MATLAB Seismic Anisotropy Toolkit

Written as part of the CoMITAC project, MSAT provides a set of Matlab functions for handling, analysing and visualising elastic tensors. This software is freely redistributable for both commercial and non-commercial use, and source files, documentation and further information are available from the MSAT homepage.

Shear wave splitting surface of the average subduction zone xenolith from the olivine fabric database of Ben Ismaïl & Mainprice (Tectonophysics, 1998), plotted using the MS_sphere function of MSAT. Red regions show directions in which splitting is smallest, and blue largest. The black ticks show the orientation of the fast shear wave, S1.

FoX — Fortran 95 library for reading/writing XML, CML and KML

Originally written as part of the NERC project eMinerals, FoX is currently maintained by Andrew Walker. It provides an extensive library for reading and writing XML (and specific flavours like CML and KML). This software is freely redistributable for both commercial and non-commercial use, and source files, documentation and further information are available from the FoX homepage.

SPECFEM3D_GLOBE dCSE project

We have instigated a project, funded via the HECToR dCSE support program, to improve the stability and performance of the SPECFEM3D_GLOBE application on HECToR and similar large supercomputers. SPECFEM3D_GLOBE implements the spectral finite element to allow the efficient simulation of the full seismic waveform for a generally anisotropic, heterogeneous model of the Earth. We use this code to provide an important link between our seismic observations and our models of the deformation of the lowermost mantle but there are a large number of other potential uses. In the dCSE project we will improve the code’s efficieny and reliability on HECToR, and add more general models of the Earth’s elasticity. The project is funded for five months’ FET effort over a one-year period, and the majority of the work is being undertaken by Paul Roberts (NAG). The project is being led by Andy Nowacki with assistance from James Wookey and Andrew Walker. We aim to merge our changes back into the main SPECFEM3D_GLOBE code base and will make patches and the project’s final report available here.

P, S and surface waves propagating away from a thrust fault located at the north pole in the anisotropic 3D mantle model S362WMANI (Kustowski et al., JGR, 2008), as calculated using SPECFEM3D_GLOBE.
Waves travelling away from an event, as calculated using SPECFEM3D_GLOBE.