Many problems in the earth sciences, particularly in the use of paleomagnetic data in plate tectonics, require estimates of the angles between directions and confidence intervals for these angles. To give a set of interrelated simple methods, we approximate all estimated directions by appropriate concentrated Fisher distributions. In preliminary numerical experiments we verify previous results for the distribution of the means of samples drawn from Fisher distributions and study the distributions of six estimators of the Fisher concentration parameter &kgr;. We selected Fisher's original estimator for &kgr;, k=(n-1)/(n-R), for use in our subsequent simulations. We then study the bias in, and confidence intervals around, the estimator &thgr;ˆ of the angle between the means of samples drawn from two Fisher distributions (''Fisher means'').

We find an approximate expression for the geometric bias in this angle, B(&thgr;)=-tan &thgr;+2&thgr;+&thgr;2crit>1/2 where &thgr;crit=<1/(R1k1)+1/(R2k2)>1/2, and show that it works well for &thgr;>2&thgr;crit. (When, for example, Rk=250 for both samples, &thgr;crit=5.1¿.) We then show that the bias-corrected estimator of &thgr;,&thgr;*=&thgr;ˆ-B(&thgr;*), is normally distributed with mean equal to the true angle and variance &sgr;&thgr;ˆ2=&thgr;crit2. Thus one can construct a 95% confidence interval for this angle with the formula (&thgr;*-1.96&sgr;&thgr;ˆ,&thgr;*+1.96&sgr;&thgr;ˆ). We show by extensive simulation that coverage probabilities for this confidence interval are slightly conservative for &thgr;>2&thgr;crit and good for &thgr;>3&thgr;crit. We derive related results for angles between various combinations of fixed directions and Fisher means. We find that the estimator αˆ of the angle between two great circles containing a known fixed direction and two different Fisher means is unbiased and normally distributed with variance &sgr;αˆ2=1/(R1k1 sin2&thgr;1)+1/(R2k2 sin 2&thgr;2) where the &thgr;i are the angular distances between the Fisher means and the known fixed direction. Thus the confidence interval on this rotation angle is (αˆ-1.96&sgr;αˆ,αˆ+1.96&sgr;αˆ). We give examples of applying these techniques to paleomagnetic data analysis, especially for determining terrane motions. Our methods provide confidence intervals for poleward displacement. For rotation, simulation suggests that they provide results at least as accurate as previous methods <McWilliams, 1984; Demarest, 1983> over a wide range of relevant parameter values, while being simpler or more flexible. We apply these methods to estimating rates of apparent polar wander and point out an additional bias due to dating errors. ¿ American Geophysical Union 1995