Urban heat island is an emerging environmental problem worldwide, and has been widely concerned by the interdisciplinary of geography and landscape ecology. Many previous studies have discovered that the impact mechanism of landscape pattern on urban heat/cold island varied across space. However, as far whether the urban hot/cold island has mutual effects and how the intensity of their interaction may change over space remain uncertain. By focusing on the western part of Shenzhen city as the case study area and based on land surface temperature (LST) and landscape types, this study constructed a two-dimensional matrix to identify the source/sink thermal landscape. Then, a contribution index (CI) was used to quantitatively analyze the evolution process of the source/sink thermal landscape. The Gravity Model was used to quantify the intensity of interaction between thermal source landscape and thermal sink landscape during the 1988-2019. The results showed that (1) the urban thermal sink patches originally surrounded by sporadic thermal source patches in 1988 were reversed in 2019. The area of thermal source landscape increased by 210.45%, while the thermal sink landscape showed a fragmentation trend, and decreased in area by 48.07%. The difference between regional average LST and that in thermal sink landscape increased from 0.49℃ in 1988 to 1.91℃ in 2019. (2) The CI of the thermal "source-sink" landscape had an increasing trend during the past thirty years. In more detail, the CI of thermal source landscape increased when the discrete thermal source patches became spatially more continuous, which aggregated LST. Meanwhile, the CI of thermal sink landscape is inhibited with the fragmentation of thermal sink landscapes. (3) In 1988-2019, the interaction between sources and sink thermal patches was relatively stable. However, the intensity of their interaction quantified by gravity model increased eight times after 30 years. Besides, the linking pattern of source and sink landscape from one sink to many sources changed from one source to many sinks, which led the intensity of link pairs to become stronger and stronger. Overall, our results indicated that the intensity of spatial interaction between urban thermal source-sink landscapes are closely related to the contribution degree of source/sink thermal landscape. Therefore, controlling the size of thermal source landscape and maintaining the size of thermal sink landscape are very important strategic for land policy makers to regulate urban thermal environment. It is also important to improve the quality of life of urban residents.