The interpretation of electron paramagnetic resonance spectra of polynuclear transition metal complexes in terms of individual contributions from each paramagnetic center can be greatly facilitated by the availability of theoretical methods that enable the reliable prediction of local spectroscopic parameters. In this work we report an approach that enables the application of multireference ab initio methods for the calculation of local zero field splitting tensors, one of the leading terms in the spin Hamiltonian for exchange-coupled systems of high nuclearity. The method referred to as local complete active space configuration interaction (L-CASCI) represents a multireference calculation with an active space composed of local orbitals of the center of interest. By successive permutation of the active space to include the localized orbitals corresponding to a particular center of the complex, all on-site parameters can be easily obtained at a high-level of theory with a corresponding low computational cost. Benchmark calculations on synthetic complexes confirm the validity of the approach. As an example of the applicability of the L-CASCI method to large systems, we determine the local anisotropy of the Mn(III) ion of the tetranuclear manganese cluster of photosystem II in both structural forms of its S2 state.