The dynamic response of two-dimensional Josephson junction arrays close to, but above the Berezinskii-Kosterlitz-Thouless (BKT) transition temperature is described in terms of the vortex dielectric function ε(ω) and the flux noise spectrum Sφ(ω). They are calculated by considering both the contributions of free vortices interacting through a screened Coulomb potential and the pair motion of vortices that are closer to each other than the BKT correlation length. This procedure allows us to understand various anomalous features in ε(ω) and in Sφ(ω) that have been observed both experimentally and in dynamic simulations.