Existing works have explored the anchor deployment for autonomous underwater vehicles (AUVs) localization under the assumption that the sound propagates straightly underwater at a constant speed. Considering that the underwater acoustic waves propagate along bent curves at varying speeds in practice, it becomes much more challenging to determine a proper anchor deployment configuration. In this paper, taking the practical variability of underwater sound speed into account, we investigate the anchor-AUV geometry problem in a 3-D time-of-flight (ToF) based underwater scenario from the perspective of localization accuracy. To address this problem, we first rigorously derive the Jacobian matrix of measurement errors to quantify the Cramer-Rao lower bound (CRLB) with a widely-adopted isogradient sound speed profile (SSP). We then formulate an optimization problem that minimizes the trace of the CRLB subject to the angle and range constraints to figure out the anchor-AUV geometry, which is multivariate and nonlinear and thus generally hard to handle. For mathematical tractability, by adopting tools from the estimation theory, we interestingly find that this problem can be equivalently transformed into a more explicit univariate optimization problem. By this, we obtain an easy-to-implement anchor-AUV geometry that yields satisfactory localization performance, referred to as the uniform sea-surface circumference (USC) deployment. Extensive simulation results validate our theoretical analysis and show that our proposed USC scheme outperforms both the cube and the random deployment schemes in terms of localization accuracy under the same parameter settings.