AbstractA probabilistic volcanic hazard assessment (PVHA) for Ceboruco volcano (Mexico) is reported using PyBetVH, an e-tool based on the Bayesian Event Tree (BET) methodology. Like many volcanoes, Ceboruco is under-monitored. Despite several eruptions in the late Holocene and efforts by several university and government groups to create and sustain a monitoring network, this active volcano is monitored intermittently rather than continuously by dedicated groups. With no consistent monitoring data available, we look at the geology and the eruptive history to inform prior models used in the PVHA. We estimate the probability of a magmatic eruption within the next time window (1 year) of ~ 0.002. We show how the BET creates higher probabilities in the absence of monitoring data, which if available would better inform the prior distribution. That is, there is a cost in terms of higher probabilities and higher uncertainties for having not yet developed a sustained volcano monitoring network. Next, three scenarios are developed for magmatic eruptions: i) small magnitude (effusive/explosive), ii) medium magnitude (Vulcanian/sub-Plinian) and iii) large magnitude (Plinian). These scenarios are inferred from the Holocene history of the volcano, with their related hazardous phenomena: ballistics, tephra fallout, pyroclastic density currents, lahars and lava flows. We present absolute probability maps (unconditional in terms of eruption size and vent location) for a magmatic eruption at Ceboruco volcano. With PyBetVH we estimate and visualize the uncertainties associated with each probability map. Our intent is that probability maps and uncertainties will be useful to local authorities who need to understand the hazard when considering the development of long-term urban and land-use planning and short-term crisis management strategies, and to the scientific community in their efforts to sustain monitoring of this active volcano.