Deepening the understanding of thermal comfort conditions of urban and rural residents is important for improving the outdoor thermal environment.Based on ERA5-Land reanalysis data,this study quantifies the thermal comfort along the urban-rural gradient(urban core-urban expansion-rural area)in 101 large cities in China during the summers from 2000 to 2020 by taking physiological equivalent temperature(PET)as the evaluation index,and explores the spatiotemporal characteristics of thermal comfort and its drivers along different gradients using the Sen's slope estimation,Mann-Kendall significance test and geodetector model.The results show that:(1)The average PET in the core areas(29.89 ℃)of case cities is slightly higher than that in the expansion areas(29.86 ℃)and much higher than that in the rural areas(28.94 ℃),and the cities with higher PET in each gradient are mainly located in the Beijing-Tianjin-Hebei region and the middle reaches of the Yangtze River.(2)PET shows an increasing trend from 2000 to 2020,with an increasing rate of 0.24℃/l0a in the core and expansion areas,and 0.22 ℃/10a in the rural areas.(3)Natural factors dominated by vegetation coverage and socioeconomic factors dominated by population density are key determinants of PET in the core and rural areas,respectively,while the drivers of PET in the expansion areas are unstable,with vegetation and impervious surface coverage playing an increasing role.As PET rises,the population exposed to uncomfortable environments increases,and this state of affairs poses new challenges for addressing thermal environmental issues.In the future,more studies on the thermal environment in urban and rural areas need to be conducted to provide a more comprehensive reference for the development of thermal adaptation strategies.