We propose a procedure for selecting basis function orientation to improve the efficiency of solution methodologies that employ local plane-wave approximations. The proposed adaptive approach consists of a local wave tracking strategy. Each plane-wave basis set within considered elements of the mesh partition is individually or collectively rotated to best align one function of the set with the local propagation direction of the field. Systematic determination of the direction of the field inside the computational domain is formulated as a minimization problem. As the resultant system is nonlinear with respect to the directions of propagation, the Newton method is employed with exact characterization of the Jacobian and Hessian. To illustrate the salient features and evaluate the performance of the proposed wave tracking approach, we present error estimates as well as numerical results obtained by incorporating the procedure into a prototypical plane-wave based approach, the least-squares method (LSM) developed by Monk et al. [1]. The numerical results obtained for the case of a two-dimensional rigid scattering problem indicate that (a) convergence was achievable to a prescribed level of accuracy, even upon initial application of the tracking wave strategy outside the pre-asymptotic convergence region, and (b) the proposed approach reduced the size of the resulting system by up to two orders of magnitude, depending * Corresponding author. Address: Department of Mathematics, California State University, Northridge, 18111 Nordhoff St., Northridge, CA 91330-8131. Tel: +1-818-677-5867 E-mail address: rabia.djellouli@csun.edu Preprint submitted to Comput. Methods Appl. Mech. Engrg. March 7, 2014 on the frequency range, with respect to the size of the standard LSM system.