A Mineral prospectivity mapping (MPM) approach using a GIS-based weighted linear combination implementation of a Multi-Criteria Evaluation approach utilising a fuzzy Analytical Hierarchy Process to elucidate expert knowledge has been implemented to analyse the spatial distribution of epithermal deposits on the Island of Milos, Greece and model their association with exploration evidence data with the aim of providing insights into the controls on ore deposition. An integrated field and Digital Airborne Imaging Spectrometer (DAIS) hyperspectral and thermal multispectral airborne remote sensing dataset supported by field mapping and laboratory analyses, has been utilised to resolve hydrothermal alteration and parameterise the MPM. This study has highlighted the intimate spatial relationship between topographic highs and locations with high grade silicified alteration at a number of locations. The ability of high spatial resolution multispectral Thermal InfraRed (TIR) remote sensing imagery, integrated with topographic data, to resolve these silicified topographic highs provides an additional tool in the exploration of epithermal deposits. The spatial relationships between silicified lithocaps, high-grade altered rocks, faulting and topographic highs were utilised in the development of the MPM model. A close association between the modelled results and the hydrothermal alteration mapped in the field supports the accuracy of this MPM approach.