Percutaneous treatment of scaphoid fractures has found increasing interest in recent years as it promises to minimize soft-tissue damage, and minimizes the risk of infections and the loss of the joint stability. However, as this procedure is mostly performed on 2D fluoroscopic images, the accurate localization of the scaphoid bone for fracture fixation renders extremely challenging. In this work, we thus propose the integration of a statistical wrist model with 3D intraoperative ultrasound for accurate localization of the scaphoid bone. We utilize a previously developed statistical wrist model and register it to bone surfaces in ultrasound images using a probabilistic approach that involves expectation-maximization. We utilize local phase symmetry to detect features in noisy ultrasound images; in addition, we use shadow information in ultrasound images to enhance and set apart bone from other features. Feasibility experiments are performed by registering the wrist model to 3D ultrasound volumes of two different wrists at two different wrist positions. And the result indicates a potential of the proposed technique for localization of the scaphoid bone in ultrasound images.