We study three dimensional bi-Hamiltonian systems in general and use the obtained results to classify all three dimensional Lotka-Volterra equations, which admit a bi-Hamiltonian representation.