It is well known that when wind turbines are deployed in large arrays, their efficiency decreases due to complex interactions among themselves and with the atmospheric boundary layer (ABL). For wind farms whose length exceeds the height of the ABL by over an order of magnitude, a "fully developed" flow regime can be established. In this asymptotic regime, changes in the streamwise direction can be neglected and the relevant exchanges occur in the vertical direction. Such a fully developed wind-turbine array boundary layer (WTABL) has not been studied systematically before. A suite of large eddy simulations (LES), in which wind turbines are modeled using the classical "drag disk" concept, is performed for various wind-turbine arrangements, turbine loading factors, and surface roughness values. The results are used to quantify the vertical transport of momentum and kinetic energy across the boundary layer. It is shown that the vertical fluxes of kinetic energy are of the same order of magnitude as the power extracted by the forces modeling the wind turbines. In the fully developed WTABL, the kinetic energy extracted by the wind turbines is transported into the wind-turbine region by vertical fluxes associated with turbulence. The results are also used to develop improved models for effective roughness length scales experienced by the ABL. The effective roughness scale is often used to model wind-turbine arrays in simulations of atmospheric dynamics at larger (regional and global) scales. The results from the LES are compared to several existing models for effective roughness lengths. Based on the observed trends, a modified model is proposed, showing improvement in the predicted effective roughness length.