Coal fires pose great threats to valuable energy resources, the ecological environment, and human safety. They are one of the most persistent fires on the Earth, which can burn for an extremely long-term period from decades to hundreds or even thousands of years. Remote sensing detection of coal fires is of significance for mitigating coal fire hazards. Nevertheless, short-term or temporal discrete land surface temperature (LST) data have limited capability in characterizing the persistent coal fire. This study proposed a methodology to monitor persistent coal fires using long-term Landsat thermal images and further to analyze spatiotemporal dynamics of coal fires. A total of 1118 high-quality Landsat images (each image containing <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$446\times446$ </tex-math></inline-formula> pixels) spanning 35 years from 1986 to 2020 in the Wuda coalfield area (China) were processed to retrieve the LST. LST time series of each pixel was decomposed into the seasonal, trend, and remainder components. Coal fire areas were demarcated by using the range of the trend components. To trace the trend and change point of the LST time series, the Mann–Kendall test was applied to the trend components, and the Pettitt test was employed to the original time series vectors of those pixels located in the coal fire areas, respectively. The random sample consensus algorithm was utilized to identify the background temperature (inliers) and high temperature (outliers) and, thus, judge the coal fire burning period, and the symbolic aggregate approximation algorithm was used to evaluate the robustness of the judgment. The calibration was conducted according to the filed surveys, obtaining spatiotemporal 3-D coal fire dynamics. The proposed methodology was testified by comparisons with fieldwork and regional anomaly extractor algorithm, demonstrating good performance in comprehensive monitoring of persistent coal fires.