LLUNPIY (lahar modeling by local rules based on an underlying pick of yoked processes, from the Quechua word “llunp’iy“, meaning flood) is a cellular automata (CA) model that simulates primary and secondary lahars, here applied to replicate those that occurred during the huge 1877 Cotopaxi Volcano eruption. The lahars flowing down the southwestern flanks of the volcano were already satisfactorily simulated in previous investigations of ours, assuming two possible different triggering mechanisms, i.e., the sudden and homogeneous melting of the summit ice and snow cap due to pyroclastic flows and the melting of the glacier parts hit by free-falling pyroclastic bombs after being upwardly ejected during the volcanic eruption. In a similar fashion, we apply here the CA LLUNPIY model to simulate the 1877 lahars sprawling out the Cotopaxi northern slopes and eventually impacting densely populated areas. Our preliminary results indicate that several important public infrastructures (among them the regional potable water supply system) and the Valle de Los Chillos and other Quito suburban areas might be devastated by northward-bound lahars, should a catastrophic Cotopaxi eruption comparable to the 1877 one occur in the near future.