In a companion paper (Clarke 1993) we described a method for the calcu
lation of body wave synthetic seismograms, referred to as the 'Complet
e Ordered Ray Expansion'. Here we present a new method for the inversi
on of body wave seismograms, based on the CORE algorithm. This multiph
ase waveform tomography allows us to incorporate many body wave arriva
ls from each seismogram in an inversion, vastly increasing the amount
of structural information that may be extracted from a set of observed
traces. We present an example of a synthetic inversion for a low-velo
city anomaly representing a mantle plume, and show that our method all
ows us to retrieve the major features of the perturbation very success
fully. Our computational scheme does not demand excessive storage, but
the amount required does grow in proportion to the number of events a
nd stations. We thus investigate the effect of partitioning the proced
ure, by initially performing trace by trace inversions, and combining
the information in a final step. The fact that small amplitude phases
are less well constrained in the preliminary inversions leads to a som
ewhat degraded final image. Nevertheless this technique is useful in s
ituations where a single, full inversion is impractical, or where data
become available over an extended period of time.