Ferroelectric thin films often exhibit a displacement of the polarization versus the electric field hysteresis loops, particularly along the electric field axis. This shift is typically attributed to structural and electronic asymmetry of the film-electrode interfaces, asymmetric surface fields, as well as space charge regions. In this study, we analyze the effect of a spatial, continuous distribution of space charge on the hysteresis response and phase transition characteristics of epitaxial ͑001͒ PbZr 0.3 Ti 0.7 O 3 thin films sandwiched between metallic electrodes on ͑001͒ SrTiO 3 substrate. Using a nonlinear thermodynamic model, we compute numerically the internal electrical fields and polarizations for several different space charge distributions both in the presence of a triangular external electric field and as a function of temperature at zero applied field. We show that space charge accumulated near the metal-ferroelectric interfaces can dramatically displace the hysteresis along the electric field axis such that the otherwise symmetric coercive fields E C − and E C + in a bulk ferroelectric related to each other through ͉E C −͉ = E C + may shift depending on the space charge concentration in such a manner that both E C − , E C + Ͼ 0 or E C − , E C + Ͻ 0. This gives rise to a very strong imprint. Our findings reveal that the presence of space charges in ferroelectric thin films results in significant changes in the phase transition characteristics, including a reduction in the phase transition temperature, smearing of the transition over a temperature range instead of a sharp dielectric anomaly at the bulk Curie temperature, and a reduction in the dielectric response compared to defect-free ferroelectrics of the same composition.