A 2-D matrix ultrasound array is used to monitor acoustic radiation force impulse (ARFI) induced shear wave propagation in 3-D in excised canine muscle. From a single acquisition, both the shear wave phase and group velocity can be calculated to estimate the shear wave speed (SWS) along and across the fibers, as well as the fiber orientation in 3-D. The true fiber orientation found using the 3-D radon transform on B-mode volumes of the muscle was used to verify the fiber direction estimated from shear wave data. For the simplified imaging case when the ARFI push can be oriented perpendicular to the fibers, the error in estimating the fiber orientation using phase and group velocity measurements was 3.5 ± 2.6° and 3.4 ± 1.4° (mean ± standard deviation), respectively, over six acquisitions in different muscle samples. For the more general case when the push is oblique to the fibers, the angle between the push and the fibers is found using the dominant orientation of the shear wave displacement magnitude. In 30 acquisitions on six different muscle samples with oblique push angles up to 40°, the error in the estimated fiber orientation using phase and group velocity measurements was 5.4 ± 2.9° and 5.3 ± 3.2°, respectively, after estimating and accounting for the additional unknown push angle. Either the phase or group velocity measurements can be used to estimate fiber orientation and SWS along and across the fibers. Although it is possible to perform these measurements when the push is not perpendicular to the fibers, highly oblique push angles induce lower shear wave amplitudes which can cause inaccurate SWS measurements.