We make a statistical study of ionospheric Joule heating with the Poynting flux method using six months of Astrid-2/EMMA electric and magnetic field data during 1999 (solar maximum year). For the background magnetic field we use the IGRF model. Our results are in agreement with earlier statistical satellite studies using both the Σ<i><sub>P</sub>E</i><sup>2</sup> method and the Poynting flux method. We present a rather comprehensive set of fitted Joule heating formulas expressing the Joule heating in given magnetic local time (MLT) and invariant latitude (ILAT) range under given solar illumination conditions as a function of the <i>K<sub>p</sub></i> index, the <i>AE</i> index, the Akasofu epsilon parameter and the solar wind kinetic energy flux. The study thus provides improved and more detailed estimates of the statistical Joule heating. Such estimates are necessary building blocks for future quantitative studies of the power budget in the magnetosphere and in the nightside auroral region.<br><br> <b>Key words.</b> Ionosphere (electric fields and currents; ionosphere-magnetosphere interactions) – Magnetospheric physics (magnetospheric configuration and dynamics)