Abstract. A total of 18 high fidelity simulations of large wind farms have been performed by three different institutions using various inflow conditions and simulation setups. The setups differ in how the atmospheric turbulence, wind shear and wind turbine rotors are modelled, encompassing a wide range of commonly used modelling methods within the LES framework. Various turbine spacings, atmospheric turbulence intensity levels and incoming wind velocities are considered. The work performed is part of the International Energy Agency(IEA) wind task Wakebench, and is a continuation of previously published results on the subject. This work aims at providing a methodology for studying the general flow behavior in large wind farms in a systematic way. It seeks to investigate and further understand the global trends of wind farm performance, with a focus on variability. Parametric studies first map the effect of various parameters on large aligned wind farms, including wind turbine spacing, wind shear and atmospheric turbulence intensity. The results are then aggregated and compared to engineering models as well as LES results from other investigations to provide an overall picture of how much power can be extracted from large wind farms operating below rated level. The simple engineering models, although they cannot capture the variability features, capture the general trends well. Response surfaces are constructed based on the large amount of aggregated LES data corresponding to a wide range of large wind farm layouts. The response surfaces form a basis for mapping the inherently varying power characteristics inside very large wind farm, including how much the turbines are able to exploit the turbulent fluctuations within the wind farms and estimating the associated uncertainty, which is valuable information useful for risk mitigation.