Featured Application: This study provides a simplified method that can be applied to help preliminarily estimate the snow load on a complex roof by combining the snow loads on several simple roofs, which have been provided in the load code or obtained by numerical studies.Abstract: A considerable number of studies have been carried out for predicting snowdrifts on roofs over the years. However, few studies have focused on snowdrifts on complex long-span roofs, as the complex shape and fine structure pose significant challenges. In this study, to simplify the calculation requirements of snow load on such roofs, work was conducted to decompose the snowdrift on a complex roof into snowdrifts on several simple roofs. First, the snow-wind tunnel test similarity criteria were investigated based on a combined air-snow-wind experimental system. Thereafter, with reference to the validated experimental similarity criteria, a series of snow-wind tunnel tests were performed for snowdrifts on a complex long-span structure under the conditions of different inflow directions. Finally, based on empirical orthogonal function (EOF) analysis, the snowdrifts on the complex roof were decomposed into basic characteristic distribution modes, including snowdrifts caused by the local and overall roof forms. The snow distribution under a specific inflow direction could be derived from the weighted combination of the basic characteristic modes, based on the wind direction coefficients. Therefore, it is possible for the snow load on a complex roof to be estimated preliminarily based on the snow distributions on several simple roofs.