The Kuramoto–Sivashinsky equation in one spatial dimension (1D KSE) is one of the most well-known and well-studied partial differential equations. It exhibits spatio-temporal chaos that emerges through various bifurcations as the domain length increases. There have been several notable analytical studies aimed at understanding how this property extends to the case of two spatial dimensions. In this study, we perform an extensive numerical study of the Kuramoto–Sivashinsky equation (2D KSE) to complement this analytical work. We explore in detail the statistics of chaotic solutions and classify the solutions that arise for domain sizes where the trivial solution is unstable and the long-time dynamics are completely two-dimensional. While we find that many of the features of the 1D KSE, including how the energy scales with system size, carry over to the 2D case, we also note several differences including the various paths to chaos that are not through period doubling.