The estimation of atmospheric turbulence parameters is of relevance for: a) site evaluation & characterisation; b) prediction of the point spread function; c) live assessment of error budgets and optimisation of adaptive optics performance; d) optimisation of fringe trackers for long baseline optical interferometry. The ubiquitous deployment of Shack-Hartmann wavefront sensors in large telescopes makes them central for atmospheric turbulence parameter estimation via adaptive optics telemetry. Several methods for the estimation of the Fried parameter and outer scale have been developed, most of which are based on the fitting of Zernike polynomial coefficients variances reconstructed from the telemetry. The non-orthogonality of Zernike polynomial derivatives introduces modal cross coupling, which affects the variances. Furthermore, the finite resolution of the sensor introduces aliasing. In this article the impact of these effects on atmospheric turbulence parameter estimation is addressed with simulations. It is found that cross coupling is the dominant bias. An iterative algorithm to overcome it is presented. Simulations are conducted for typical ranges of the outer scale (4 to 32m), Fried parameter (10 cm) and noise in the variances (signal-to-noise ratio of 10 and above). It is found that, using the algorithm, both parameters are recovered with sub-percent accuracy.