The source-receiver Green’s functions of Rayleigh and Love waves are determined using data recorded of events. In order to perform a 2-D tomography procedure, the group velocity dispersion curves (data parameter) of Rayleigh and Love waves are calculated using signals of micro-earthquakes with magnitudes less than 4. The group velocity dispersion curves are obtained for each geographic cell applying the resulting 2D maps. Then, by performing a 1-D nonlinear inversion process for the data parameter (dispersion curves in each node), the vertical and horizontal shear wave velocity is calculatedin each geographic cell. Using the vertical and horizontal shear wave velocities at each depth, radial anisotropy is calculated in each cell as a function of depth. Therefore, by adjoining/putting together the 1-Danisotropy functions, 2-D maps of anisotropy can be drawn/depicted.