Numerical simulations of single-cell, two-dimensional, time-dependent thermal convection in a square cross-section of fluid-saturated porous material heated uniformly from below reveal a series of transitions between distinct oscillatory dynamical regimes. With increasing Rayleigh number R, the flow first evolves from steady-state behaviour into periodic motion with a single frequency f which depends on R approximately according to $f\propto R^{\frac{7}{8}}$ the transition Rayleigh number lies between about 380 and 400. At a value of R between about 480 and 500 the flow transforms into a fluctuating state characterized by two frequencies. Soon thereafter, for R between about 500 and 520, it reverts back to single-frequency periodic behaviour with f approximately proportional to $R^{\frac{3}{2}}$. The two frequencies in the narrow transition regime may be locked to a rational ratio, in which case the flow is periodic, or they may be commensurate, in which case the flow is quasi-periodic. The spectral characteristics of numerical realizations of unsteady convection and the occurrences of transitions therein are highly dependent on truncation level in Galerkin schemes or resolution in finite-difference approaches.