We uncover novel features of three-dimensional natural convection in porous media by investigating convection in an annular porous cavity contained between two vertical coaxial cylinders. The investigations are made using a linear stability analysis, together with high-order numerical simulations using pseudospectral methods to model the nonlinear regime. The onset of convection cells and their preferred planform are studiedthe stability of the modes with respect to different types of perturbation is investigated. We also examine how variations in the Rayleigh number affect the convection modes and their stability regimes. Compared with previously published data, we show how the problem inherits an increased complexity regarding which modes will be obtained. Some stable secondary modes or mixed modes have been identified and some overlapping stability regions for different convective modes are determined.