A local seismic network was installed and operated in the area of Florina, northern Greece, during July 2013 – January 2014 for studying the high microseismic activity following the occurrence of an Mw=4.1 event on 17 February 2013. The recordings of the local network along with the ones of the Hellenic Unified Seismological Network (HUSN) are used for obtaining accurate locations and defining the characteristics of the seismic activity. A new velocity model is calculated for a broader area using the recordings of the HUSN stations. Relocation is performed for 1330 events recorded by the local network and 423 events recorded by HUSN using the double difference technique and cross correlation measurements. Fault plane solutions are determined for the two largest events (Mw=3.6 and Mw=4.1) using waveform inversion technique. The causative fault of the largest event (Mw=4.1) is striking almost E-W and dipping to the north. However, most of the activity is concentrated in its south and forms an almost vertical, south dipping plane, striking almost E-W. This southern cluster consists of multiplets and is located in an area of positive Coulomb stress changes, due to the coseismic slip of the Mw=4.1 earthquake. The spatio-temporal evolution of microseismicity is in most cases well explained by diffusion curves, a signature of fluid induced seismicity. However, there are cases where a triggering front could not describe the evolution of the activity resulting to an earthquake - earthquake triggering. The dependent or independent (Poisson process) earthquake occurrence was examined by searching their interevent time distribution. These characteristics are interpreted as the consequence of CO2 emission through faults which are used as pathways in an area enhanced by the coseismic Coulomb stress changes.