Discovering hidden geothermal signatures using non-negative matrix factorization with customized k-means clustering. (December 2022)