AbstractHimalayan earthquakes have deep societal and economic impact. In this article, we implement a surrogate method of nowcasting (Rundle et al., 2016) to determine the current state of seismic hazard from large earthquakes in a dozen populous cities from India and Pakistan that belong to the west-northwest part of Himalayan orogeny. For this, we (1) perform statistical inference of natural times, intersperse counts of small-magnitude events between pairs of succeeding large events, based on a set of eight probability distributions; (2) compute earthquake potential score (EPS) of 14 cities from the best-fit cumulative distribution of natural times; and (3) carry out a sensitivity testing of parameters—threshold magnitude and area of city region. Formulation of natural time (Varostos et al., 2005) based on frequency–magnitude power-law statistics essentially avoids the daunting need of seismicity declustering in hazard estimation. A retrospective analysis of natural time counts corresponding to M≥6 events for the Indian cities provides an EPS (%) as New Delhi (56), Chandigarh (86), Dehradun (83), Jammu (99), Ludhiana (89), Moradabad (84), and Shimla (87), whereas the cities in Pakistan observe an EPS (%) as Islamabad (99), Faisalabad (88), Gujranwala (99), Lahore (89), Multan (98), Peshawar (38), and Rawalpindi (99). The estimated nowcast values that range from 38% to as high as 99% lead to a rapid yet useful ranking of cities in terms of their present progression to the regional earthquake cycle of magnitude ≥6.0 events. The analysis inevitably encourages scientists and engineers from governments and industry to join hands for better policymaking toward land-use planning, insurance, and disaster preparation in the west-northwest part of active Himalayan belt.