Tacora Volcano (17 degrees 43'S-69 degrees 46'W) lies at the southernmost end of a 10 km-long volcanic lineament that extends between Chile and Peru. Around Tacora volcano, thermal manifestations are two active fumarolic fields located at the western flank of the stratovolcano and at the volcano summit, indicating active magma degassing in a shallow hydrothermal system. Beneath Tacora volcano is located the NW Challaviento reverse fault that belongs to the Incapuquio - Challaviento fault system of Middle Eocene age. To complement previous exploration results and conceptual modeling developed by INFINERGEO SPA, seventeen short period seismic stations were installed around Tacora Volcano, between August and December 2014. Using the P and S wave arrival times of locally recorded seismicity, a 3D velocity model was determined through a travel time tomography. We interpreted high Vp/Vs values as water-saturated areas, corresponding to the recharge zone of Tacora hydrothermal system. In addition, low values of Delta Vp/Vp (%) and Vp/Vs ratio represent the location of a gas-saturated magmatic reservoir and circulation networks of magmatic-hydrothermal fluids. Low Vp/Vs volumes (magma reservoirs/high temperature hydrothermal fluids), the presence of fumarolic fields and surface hydrothermal alteration have a spatial correlation. This suggests a structural control of the Challaviento fault in the hydrothermal flow. Finally, we present a cluster analysis using the Delta Vp/Vp (%) parameter. Through this analysis, we found a method for the identification of a key structure in depth composed by the magma reservoir (low Vp/Vs ratios, low Delta Vp/Vp (%)), clay level areas (intermediate values of Delta Vp/Vp (%)), and degassification zones (low values of Delta Vp/Vp (%)) directly related with the surface thermal manifestations.