Persistent Urban Park Cooling Effects in Krakow: A Satellite-Based Analysis of Land Surface Temperature Patterns (1990–2018)

Urban green spaces provide measurable cooling that can mitigate urban heat islands, yet few studies have quantified these effects over multiple decades. This study analyzed Landsat imagery from four epochs (1990, 2000, 2013, 2018) to derive land surface temperature (LST) and vegetation indices—NDVI for greenness and NDMI for moisture content—for four large urban parks in Krakow. Late spring/summer LST in parks was compared with that of urban areas within 0–150 m and 150–300 m of park boundaries. Statistical significance was evaluated using bootstrapped confidence intervals, long-term trends were assessed via the Mann–Kendall test, and correlation analysis was used to examine relationships between LST and each vegetation index. Results show a persistent park cooling effect, with park interiors ~2–3 °C cooler than adjacent urban areas in all years. Despite an overall city-wide LST rise of ~5–6 °C from 1990 to 2018, the park cool island intensity (temperature difference between park and city) remained stable (no significant long-term trend, p > 0.7). Bootstrapped 95% confidence intervals confirmed that each park’s cooling effect was statistically significant in each year analyzed. NDMI (vegetation moisture content) correlated more strongly with LST (r ~ −0.90) than NDVI (r ~ −0.7 to −0.9), highlighting the importance of vegetation moisture in park cooling. These findings demonstrate that well-watered urban parks can sustain substantial cooling benefits over decades of urban development. The persistent ~2–3 °C daytime cooling observed underscores the value of water-sensitive green space planning as a long-term urban heat mitigation strategy.