The most common clinical phenotype of progressive supranuclear palsy is Richardson syndrome, characterized by levodopa unresponsive symmetric parkinsonism, with a vertical supranuclear gaze palsy, early falls and cognitive impairment. There is currently no detailed understanding of the full sequence of disease pathophysiology in progressive supranuclear palsy. Determining the sequence of brain atrophy in progressive supranuclear palsy could provide important insights into the mechanisms of disease progression, as well as guide patient stratification and monitoring for clinical trials. We used a probabilistic event-based model applied to cross-sectional structural MRI scans in a large international cohort, to determine the sequence of brain atrophy in clinically diagnosed progressive supranuclear palsy Richardson syndrome. A total of 341 people with Richardson syndrome (of whom 255 had 12-month follow-up imaging) and 260 controls were included in the study. We used a combination of 12-month follow-up MRI scans, and a validated clinical rating score (progressive supranuclear palsy rating scale) to demonstrate the longitudinal consistency and utility of the event-based model's staging system. The event-based model estimated that the earliest atrophy occurs in the brainstem and subcortical regions followed by progression caudally into the superior cerebellar peduncle and deep cerebellar nuclei, and rostrally to the cortex. The sequence of cortical atrophy progresses in an anterior to posterior direction, beginning in the insula and then the frontal lobe before spreading to the temporal, parietal and finally the occipital lobe. This in vivo ordering accords with the post-mortem neuropathological staging of progressive supranuclear palsy and was robust under cross-validation. Using longitudinal information from 12-month follow-up scans, we demonstrate that subjects consistently move to later stages over this time interval, supporting the validity of the model. In addition, both clinical severity (progressive supranuclear palsy rating scale) and disease duration were significantly correlated with the predicted subject event-based model stage (P < 0.01). Our results provide new insights into the sequence of atrophy progression in progressive supranuclear palsy and offer potential utility to stratify people with this disease on entry into clinical trials based on disease stage, as well as track disease progression.