During geomagnetic storms, enhanced solar wind and changes in the interplanetary magnetic field lead to ionisation anomalies across the polar regions. The superstorm of 20 November 2003 was one of the largest events in recent history. Numerical simulations of ionospheric dynamics during the storm are compared with plasma observations to understand the mechanisms forming the polar plasma anomalies. The results are important for understanding and forecasting space weather in polar regions.