Electron paramagnetic resonance (EPR) of Fe 3+ ions in Al 2 O 3 is studied in powder samples prepared by different routes and/or modified by thermal or mechanical treatments, with different doping levels and grain sizes. The measurements are performed in various frequency bands (S, X, K, Q and W) and with bimodal detection in X-band. Simulations of the spectra are achieved with a code designed for computing EPR powder spectra described by any spin Hamiltonian including second-, fourth-and sixth-order ZFS terms (S 7/2). The linewidths, intensities and lineshapes are accounted for. The lineshape is Gaussian at low Fe 3+ concentration whereas it is Lorentzian for higher concentration. The linewidths are interpreted as the superimposition of three main contributions: intrinsic linewidth, dipolar broadening and broadening due to lattice imperfections. The latter is tentatively interpreted in terms of quadrupolar spin Hamiltonian parameter distributions treated using firstorder perturbation theory. Whatever the sample, only the b 2 2 spin Hamiltonian parameter is found to be distributed around a mean zero value which corresponds to rhombic distortions. Angle and bond length distributions are tentatively extracted from the b 2 2 distributions which gives some insight into the local order around the spin probe in relation to the preparation and treatment of the samples.