We present a 4-D tomography model suitable for Active Monitoring application to waveform data obtained from accurately repeatable controlled point-sources such as seismic vibrators, piezo-electric ground force sources, or electro-magnetic sources. In our mathematical seismic wave propagation model we consider the static stress gradient, anelasticity, and anisotropy as perturbations in an inhomogeneous, isotropic background medium for which the Green tensor is described by body wave ray theory. The perturbations result in contrast-in-medium source integral representations for the total particle velocity field, with separable kernels for the three perturbations. The contrasts are to be resolved from inversion of the particle velocity seismograms observed at receiver stations. Continuous time-lapse travel time and contrast source inversion tomographic monitoring of the Lithosphere may reveal emergent heterogeneities potentially representing earthquake preparation zones of which the geometry, orientation, and physics parameter evolution can be determined, thus providing reliable and accurately repeatable data for earthquake precursor and pre-earthquake process studies. In additional applications, the method may enhance earthquake engineering, earthquake and explosion relocation procedures, seismograph network calibration, Comprehensive Test Ban Treaty verification, and oil, gas, and water reservoir parameter monitoring.