Paleo-earthquakes along the Cascadia subduction zone inferred from offshore sediments and Japan coastal tsunami deposits approximated to M9+ and ruptured the entire margin. However, due to the lack of modern megathrust earthquake records and general quiescence of subduction fault seismicity, the potential megathrust rupture scenario and influence of downdip limit of the seismogenic zone are still obscure. In this study, we present a numerical simulation of Cascadia subduction zone earthquake sequences in the laboratory-derived rate-and-state friction framework to investigate the potential influence of the geodetic fault locking on the megathrust sequences. We consider the rate-state friction stability parameter constrained by geodetic fault locking models derived from decadal GPS records, tidal gauge and levelling-derived uplift rate data along the Cascadia margin. We incorporate historical coseismic subsidence inferred from coastal marine sediments to validate our coseismic rupture scenarios. Earthquake rupture pattern is strongly controlled by the downdip width of the seismogenic, velocity-weakening zone and by the earthquake nucleation zone size. In our model, along-strike heterogeneous characteristic slip distance is required to generate margin-wide ruptures that result in reasonable agreement between the synthetic and observed coastal subsidence for the AD 1700 Cascadia Mw∼9.0 megathrust rupture. Our results suggest the geodetically inferred fault locking model can provide a useful constraint on earthquake rupture scenarios in subduction zones. This article is part of the theme issue 'Fracture dynamics of solid materials: from particles to the globe'.