A number of phenomena and processes in geosciences can be summarized by second order partial differential equations. The major numerical methods for their snlution include the classical finite difference method and the finite element method newly develope