Local transient extreme wind loads caused by group tower-related interference are among the major reasons that lead to wind-induced damage of super-large cooling towers. Four-tower arrangements are the most commonly seen patterns for super-large cooling towers. We considered five typical four-tower arrangements in engineering practice, namely, single row, rectangular, rhombic, L-shaped, and oblique L-shaped. Wind tunnel tests for rigid body were performed to determine the influence of different arrangements on static and dynamic wind loads and extreme interference effect. The most unfavorable working conditions (i.e., the largest overall wind loads) were determined based on the overall aerodynamic coefficient under different four-tower arrangements. Then we calculated the one-, two- and three-dimensional aerodynamic loads under different four-tower arrangements. Statistical analyses were performed on the wind pressure signals in the amplitude and time domains under the most unfavorable working conditions. On this basis, the non-Gaussian distribution characteristics of aerodynamic loads on the surface of the cooling towers under different four-tower arrangements were analyzed. We applied the Sadek-Simiu procedure to the calculation of two- and three-dimensional aerodynamic loads in the cooling towers under the four-tower arrangements, and the extreme wind load distribution patterns under the most unfavorable working conditions in each arrangement were compared. Finally, we proposed a uniform equation for fitting the extreme wind loads under the four-tower arrangements; the accuracy and reliability of the equation were verified. Our research findings will contribute to the optimization of the four-tower arrangements and the determination of extreme wind loads of super-large cooling towers.