The net storage heat flux (ΔQ[subscript S]) is important in the urban surface energy balance (SEB) but its determination remains a significant challenge. The hysteresis pattern of the diurnal relation between the ΔQ[subscript S] and net all-wave radiation (Q[superscript ∗]) has been captured in the Objective Hysteresis Model (OHM) parameterization of ΔQ[subscript S]. Although successfully used in urban areas, the limited availability of coefficients for OHM hampers its application. To facilitate use, and enhance physical interpretations of the OHM coefficients, an analytical solution of the one-dimensional advection–diffusion equation of coupled heat and liquid water transport in conjunction with the SEB is conducted, allowing development of AnOHM (Analytical Objective Hysteresis Model). A sensitivity test of AnOHM to surface properties and hydrometeorological forcing is presented using a stochastic approach (subset simulation). The sensitivity test suggests that the albedo, Bowen ratio and bulk transfer coefficient, solar radiation and wind speed are most critical. AnOHM, driven by local meteorological conditions at five sites with different land use, is shown to simulate the ΔQ[subscript S] flux well (RMSE values of ∼ 30 W m[superscript −2]). The intra-annual dynamics of OHM coefficients are explored. AnOHM offers significant potential to enhance modelling of the surface energy balance over a wider range of conditions and land covers.