Motivation: Effective management of liver tumors through thermal ablation requires precise monitoring of the ablation zone to ensure successful treatment outcomes. Computed tomography (CT) thermometry offers a promising non-invasive solution to monitor if tumor cells have been heated to the lethal temperature threshold. However, achieving reproducible, precise, and accurate temperature measurements remains a challenge, particularly due to metal artifacts introduced by the ablation equipment. Purpose: This study investigates the applicability of spectral CT thermometry in monitoring liver microwave ablation. It compares the reproducibility, precision and accuracy of CT thermometry on attenuation value images, with CT thermometry on physical density maps using spectral CT. Furthermore, it identifies the optimal metal artifact reduction (MAR) method — among O-MAR, deep learning-MAR, spectral CT, or a combination — to reduce needle artifacts and improve CT thermometry precision. Materials and Methods: Four liver-mimicking gel phantoms embedded with temperature sensors underwent a 10-minute, 60W microwave ablation imaged by dual-layer spectral CT using a Philips CT7500 scanner. Each scan was processed to reconstruct standard 120 kVp images alongside physical density maps, which were derived from virtual monochromatic imaging (70 - 150 keV) and effective atomic number maps. During each procedure, 23 CT scans were acquired to monitor attenuation and physical density values in proximity of the ablation antenna over time. Attenuation-based and physical density-based thermometry models were tested for reproducibility (coefficient of variation) over three repetitions; a fourth repetition focused on accuracy (Bland-Altman analysis). MAR techniques were applied to a single repetition to evaluate temperature precision in artifact-corrupted slices. Results: The correlation between CT value and temperature was highly linear with an R-squared value exceeding 96% for both attenuation and physical density-based thermometry. Model parameters for attenuation-based and physical density-based thermometry were -0.38 HU/ºC and 0.00039 ºC-1, with coefficients of variation of 0.023 and 0.067, respectively, indicating a high reproducibility. CT thermometry precision increased with distance from the ablation antenna, the use of attenuation maps and deep learning-MAR. Physical density maps generated at 150 keV alone and in combination with O-MAR and deep learning-MAR reduced needle artifacts by 73% on average (p=0.003) compared to attenuation images. Bland-Altman analysis reveals limits-of-agreement of -7.7°C to 5.3°C and -9.5°C to 8.1°C for attenuation and physical density-based thermometry, respectively. Conclusion: Spectral CT has the potential to make CT thermometry more universally applicable. This study demonstrates the effectiveness of spectral CT thermometry for non-invasive temperature monitoring during liver microwave ablation. It shows that using spectral physical density maps at 150 keV, alongside deep learning-MAR and O-MAR, enhances temperature accuracy and minimizes metal artifacts. However, standardizing thermometry parameters across different patient conditions remains a challenge. Future enhancements in photon counting CT and deep learning technologies could further refine this method, ultimately reducing the risk of local tumor recurrence.