A variety of tectonic models can be constructed by superposing the displacement field solution of a single force operative at a point in an elastic half-space. Piezomagnetic changes associated with the same models are given as well by the linear combination of fundamental piezomagnetic potentials, which arise from stress-induced magnetization produced by a unit single force at a point. The earth is assumed a homogeneous and isotropic elastic semi-infinite medium having a uniformly magnetized top layer. The method is adaptable to (a) surface load, (b) strain nuclei, (c) dislocation and (d) volume source problems. The prescribed potentials are obtained with the aid of a surface integral representation of the tectonomagnetic field (SASAI, 1983). As an application example, piezomagnetic field is calculated for a uniform circular load.