Rapid urbanization and accelerating industrial activity have profoundly altered the thermal environments of cities, disrupting the natural surface–atmosphere energy balance that governs near-surface temperature regulation. One of the most persistent and well-documented manifestations of this transformation is the Urban Heat Island (UHI) effect, whereby urban areas exhibit systematically higher temperatures than their surrounding rural counterparts [1]. This urban–rural thermal contrast arises from a combination of reinforcing physical and anthropogenic mechanisms. Impervious construction materials such as asphalt, concrete, and roofing systems possess high thermal inertia, enabling efficient heat absorption during daytime and delayed nocturnal release. Concurrently, the reduction of vegetative cover diminishes shading and evapotranspiration, while dense urban morphology modifies airflow and suppresses convective cooling. These processes are further intensified by anthropogenic heat emissions from buildings, transportation systems, and industrial activities [2-4].

The implications of the UHI effect extend well beyond thermal discomfort. Elevated urban temperatures contribute to increased electricity demand for space cooling, exacerbate air pollution through enhanced photochemical reactions, heighten heat-related morbidity and mortality, and compound the impacts of regional and global climate change [5]. As climate warming intensifies and urban populations continue to expand, the need for robust and reliable modeling tools capable of capturing urban thermal dynamics has become increasingly urgent.

Historically, mathematical modeling of urban heat processes has relied predominantly on classical differential equations, particularly diffusion-based heat equations, energy balance formulations, and empirical or statistical approaches [6,7]. These frameworks have provided foundational insight into heat propagation, surface–atmosphere exchange, and temperature evolution in urban settings. However, classical models typically rest on assumptions of locality, integer-order dynamics, and short-term dependence, implicitly treating heat transfer as a Markovian process. In real urban environments, thermal behavior often departs from these assumptions. Built materials exhibit long-term heat storage and delayed cooling, urban surfaces respond heterogeneously to forcing due to land-cover diversity, and temperature anomalies may persist long after the removal of external drivers. Such phenomena suggest the presence of memory effects, nonlocal interactions, and anomalous diffusion that challenge the representational capacity of conventional integer-order models.

Fractional calculus provides a mathematically rigorous extension of classical calculus by generalizing differentiation and integration to non-integer orders, thereby enabling the explicit incorporation of memory and hereditary properties into governing equations [8,9]. Over       the past two decades, fractional differential equations (FDEs) have been successfully applied to a wide range of physical systems characterized by long-memory behavior, including anomalous diffusion, viscoelasticity, porous media flow, and thermal processes in heterogeneous materials [10,11]. The defining feature of fractional operators—their dependence on the entire past history of the system—renders them particularly suitable for modeling processes in which present states are strongly influenced by accumulated past dynamics.

In the context of urban thermal systems, fractional calculus offers a compelling modeling framework. Urban heat retention and release are intrinsically history-dependent, shaped by cumulative solar radiation, material properties, and repeated anthropogenic forcing. A fractional-order heat equation can naturally encode these persistent temporal correlations, while also providing a mechanism for representing spatially anomalous diffusion arising from heterogeneous urban fabrics. As such, fractional models have the potential to capture both the intensity and persistence of the UHI effect more realistically than classical diffusion-based formulations.

Despite the rapid growth of UHI research, the majority of existing modeling studies continue to rely on classical heat equations and statistical regressions that implicitly assume predominantly local behavior and limited memory. Yet urban heat dynamics are governed by long-term coupling between natural drivers—such as radiation, wind, and humidity—and anthropogenic influences, including land-use change, infrastructure expansion, and emissions growth. These coupled interactions often generate nonlinear responses and delayed system adjustment, suggesting that conventional models may underrepresent the persistence of urban heat and the lagged response of temperatures to changing boundary conditions. This limitation can undermine predictive reliability, particularly in applications where accurate estimation of heat duration and recovery time is critical.

The implications of this modeling gap are significant for urban planning and climate adaptation. Heat mitigation strategies—such as green infrastructure deployment, reflective and permeable materials, urban ventilation corridors, and land-use redesign—require predictive tools capable of assessing not only instantaneous temperature reductions but also the long-term evolution and memory of urban heat accumulation. Models that fail to account for memory-driven dynamics may underestimate the effectiveness or persistence of mitigation measures, thereby limiting their utility for evidence-based policy design.

Motivated by these considerations, the present study develops a fractional differential equation–based framework for modeling the Urban Heat Island effect and systematically compares its behavior with that of a classical integer-order counterpart. The overarching objective is to construct and analyze a fractional-order heat diffusion model that captures the nonlocal, memory-dependent, and heterogeneous nature of urban thermal processes. Specifically, the study synthesizes classical and fractional formulations of heat transfer relevant to urban environments; formulates a fractional-order UHI model under appropriate initial and boundary conditions; derives analytical or semi-analytical solutions where possible; and evaluates model performance through comparison with a classical model in terms of descriptive adequacy and predictive behavior.

The contributions of this work are multi-dimensional. From a theoretical perspective, the study advances urban climate modeling by embedding UHI dynamics within a fractional calculus framework that explicitly incorporates memory and nonlocality. From a practical standpoint, it provides an enhanced modeling tool that can support the evaluation of urban heat mitigation strategies by improving representation of persistent heat behavior. Environmentally and socially, improved prediction of UHI intensity and duration has the potential to inform interventions aimed at reducing heat-related health risks, lowering cooling-energy demand, and mitigating emissions associated with urban energy consumption. Academically, the study extends the application frontier of fractional calculus to a pressing climate–urbanization challenge, reinforcing the role of applied mathematics in sustainability research.

The scope of the study is confined to fractional differential equation modeling of the UHI effect, with emphasis on fractional-order heat equations that represent memory effects in urban heat retention and release. Empirical illustration relies on representative urban–rural temperature datasets or idealized data for comparative analysis, rather than fully coupled atmospheric simulations. While the model is formulated in a general form and is not geographically specific, illustrative examples are used to demonstrate interpretation and performance.

Finally, the study acknowledges inherent limitations. High-resolution temperature datasets suitable for detailed calibration may be unavailable in some urban contexts. Fractional-order models can be computationally demanding due to their intrinsic history dependence. Moreover, while the proposed framework captures broad-scale UHI persistence, it does not explicitly resolve micro-scale processes associated with fine urban geometry or transient localized activities. These limitations also point toward future research directions, including improved data assimilation, multiscale refinement, and coupling of fractional thermal models with comprehensive urban climate simulation systems.

For Full-length Manuscript, please go through the PDF File.