Abstract. Leading edge erosion (LEE) is one of the most critical degradation mechanisms that occur with wind turbine blades (WTBs), generally starting from the tip section of the blade. A detailed understanding of the LEE process and the impact on aerodynamic performance due to the damaged leading edge (LE) is required to select the most appropriate leading edge protection (LEP) system and optimize blade maintenance. Providing accurate modeling tools is therefore essential. This paper presents a two-part study investigating computational fluid dynamics (CFD) modeling approaches for different orders of magnitudes in erosion damage. The first part details the flow transition modeling for eroded surfaces with roughness on the order of 0.1–0.2 mm, while the second part focuses on a novel study modeling high-resolution scanned LE surfaces from an actual blade with LEE damage on the order of 10–20 mm (approx. 1 % chord); 2D and 3D surface-resolved Reynolds-averaged Navier–Stokes (RANS) CFD models have been applied to investigate wind turbine blade sections in the Reynolds number (Re) range of 3–6 million. From the first part, the calibrated CFD model for modeling flow transition accounting for roughness shows good agreement of the aerodynamic forces for airfoils with leading-edge roughness heights on the order of 140–200 µm while showing poor agreement for smaller roughness heights on the order of 100 µm. Results from the second part of the study indicate that up to a 3.3 % reduction in annual energy production (AEP) can be expected when the LE shape is degraded by 0.8 % of the chord, based on the NREL5MW turbine. The results also suggest that under fully turbulent conditions, the degree of eroded LE shapes studied in this work show the minimal effect on the aerodynamic performances, which results in a negligible difference to AEP.