PurposeDeveloping severe load spectrum of transport aircraft structures is crucial for enhancing the fatigue damage correlation between full-scale fatigue testing results and operational service. The lack of consensus on severe spectrum development methods for transport aircraft has prompted the current research, resulting in a proposed approach for a severe gust load spectrum based on the acceleration cumulative exceedance surface.Design/methodology/approachThe measured load data were analyzed using a model based on the cumulative exceedance number surface to describe the variation in exceedance numbers. An improved sampling method based on multivariate Markov Chain Monte Carlo was employed to obtain the fleet fatigue damage distribution, enabling the determination of the severity of severe spectrum and the corresponding cumulative exceedance number surface, and a severe gust load spectrum was developed based on the surface.FindingsThe method that characterizes load spectrum variation using the cumulative exceedance surface minimizes the randomness of peak-trough pairs by incorporating the correlation of load spectrum peaks and troughs. This approach reduces the variation in fleet fatigue damage, thereby lowering the requirements for the severity of severe spectrum fatigue damage.Originality/valueThe proposed methodology extends from a one-dimensional curve to a two-dimensional surface, accounting for the correlation between peak and trough values to develop a severe spectrum. This approach more accurately describes the variation in acceleration cumulative exceedance numbers, directly benefiting fatigue damage calculation. This study provides valuable references for developing severe spectrum for transport aircraft.