This version is available at https://strathprints.strath.ac.uk/52305/ Strathprints is designed to allow users to access the research output of the University of Strathclyde. Unless otherwise explicitly stated on the manuscript, Copyright © and Moral Rights for the papers on this site are retained by the individual authors and/or other copyright owners. Please check the manuscript for details of any other licences that may have been applied. You may not engage in further distribution of the material for any profitmaking activities or any commercial gain. You may freely distribute both the url (https://strathprints.strath.ac.uk/) and the content of this paper for research or private study, educational, or not-for-profit purposes without prior permission or charge.Any correspondence concerning this service should be sent to the Strathprints administrator: strathprints@strath.ac.ukThe Strathprints institutional repository (https://strathprints.strath.ac.uk) is a digital archive of University of Strathclyde research outputs. It has been developed to disseminate open access research outputs, expose data about those outputs, and enable the management and persistent access to Strathclyde's intellectual output. A systematic investigation of unstable steady-state solutions of the Darcy-OberbeckBoussinesq equations at large values of the Rayleigh number Ra is performed to gain insight into two-dimensional porous medium convection in domains of varying aspectratio L. The steady convective states are shown to transport less heat than the statistically steady 'turbulent' flow realised at the same parameter values: the Nusselt number N u ∼ Ra for turbulent porous medium convection, while N u ∼ Ra 0.6 for the maximum heat-transporting steady solutions. A key finding is that the lateral scale of the heat-fluxmaximising solutions shrinks roughly as L ∼ Ra −0.5 , reminiscent of the decrease of the mean inter-plume spacing observed in turbulent porous medium convection as the thermal forcing is increased. A spatial Floquet analysis is performed to investigate the linear stability of the fully nonlinear steady convective states, extending a recent study by Hewitt et al. (J. Fluid Mech. 737, 2013) by treating a base convective state -and secondary stability modes -that satisfy appropriate boundary conditions along plane parallel walls. As in that study, a bulk instability mode is found for sufficiently small aspect-ratio base states. However, the growth rate of this bulk mode is shown to be significantly reduced by the presence of the walls. Beyond a certain critical Ra-dependent aspect-ratio, the base state is most strongly unstable to a secondary mode that is localised near the heated and cooled walls. Direct numerical simulations, strategically initialised to investigate the fully nonlinear evolution of the most dangerous secondary instability modes, suggest that the (long time) mean inter-plume spacing in statistically-steady porous medium convection results from a balance between the competing effects of these two types of inst...