With the trend of device miniaturization, ultrathin ferroelectric films are gaining more and more attention. However, understanding ferroelectricity in this nanoscale context remains a formidable challenge, primarily due to the heightened relevance of surface effects, which often leads to the loss of net polarization. Here, the influence of surface effects on the polarization as a function of thickness in ultrathin BiFeO 3 films is investigated using phase-field simulations. The findings reveal a notable increase in ferroelectric polarization with increasing thickness, with a particularly discernible change occurring below the 10 nm threshold. Upon accounting for surface effects, the polarization is marginally lower than the case without such considerations, with the disparity becoming more pronounced at smaller thicknesses. Moreover, the hysteresis loop and butterfly loop of the ultrathin film were simulated, demonstrating that the ferroelectric properties of films remain robust even down to a thickness of 5 nm. Our investigations provide valuable insights into the significance of ferroelectric thin films in device miniaturization.