Abstract. The Fuego de Colima volcano (Mexico) showed a complex eruptive behaviour with periods of rapid and slow lava dome growth, punctuated by explosive activity. We reconstructed the weekly discharge rate average between 1998 and 2018 by means of satellite thermal data integrated with published discharge rate data. By using spectral and wavelet analysis, we found a multi-year long-, multi-month intermediate-, and multi-week short-term cyclic behaviour during the period of the investigated eruptive activity, as those of many others dome-forming volcanoes. We use numerical modelling in order to investigate the non-linear cyclic eruptive behaviour considering a magma feeding system composed of a dual or a single magma chamber connected to the surface through an elastic dyke evolving into a cylinder conduit in the shallowest part. We investigated the cases in which the periodicity is controlled by i) the coupled deep-shallow magma reservoirs, ii) the single shallow chamber, and iii) the elastic shallow dyke when is fed by a fixed influx rate or a constant pressure. The model outputs indicate that the observed multi-year periodicity (1.5–2.5 years) can be described by the fluctuations controlled by a shallow magma chamber with a volume of 20–50 km3 coupled with a deep reservoir of 500 km3, connected through a deep elastic dyke. The multi-month periodicity (ca. 5–10 months) appears to be controlled by the shallow magma chamber for the same range of volumes. The short-term multi-week periodicity (ca. 2.5–5 weeks) can be reproduced considering a fixed influx rate or constant pressure at the base of the shallower dyke. This work provides new insights on the non-linear cyclic behaviour of Fuego de Colima, and a general framework for the comprehension of eruptive behaviour of andesitic volcanoes.