Shelterbelts are increasingly used to mitigate strong wind damage, but the complex canopy structures create challenges for numerical studies of windbreak effectiveness, such as the trade-off between computational cost and accuracy of results. To address these challenges and accurately investigate the downstream wind fields, most conventional studies represent shelterbelts as rectangular porous media with a uniformly distributed aerodynamic resistance coefficient. However, due to the vertical variation in canopy diameter and the irregular distribution of leaf density, the aerodynamic resistance of natural shelterbelts becomes nonuniform accordingly. To quantify the discrepancies arising from this simplification, this study first proposes a non-destructive approach to calculate canopy porosity profiles, which are further used to derive aerodynamic resistance at different heights. Then, by comparing the results obtained from the conventional and proposed approaches in Large-Eddy Simulations, the discrepancies caused by ignoring the vertical variation in canopy structures are analyzed. Finally, these discrepancies are further investigated for double-row shelterbelts. The results show that ignoring the vertical variation in canopy diameter leads to significant differences in windbreak effectiveness, especially for the downstream velocity and pressure fields at the top and middle heights of the canopy. The proposed approach provides a computationally efficient and more accurate representation of near-surface wind fields downstream of shelterbelts, thereby contributing to the accurate prediction of local wind fields for meteorological services.