Range errors constrain treatment planning by limiting choice of ion beam angles and requiring large margins. Ionoacoustic range verification requires recovering the location of an acoustic source from low frequency signals. A priori information is applied to stably overcome resolution limits of inverse acoustic source imaging in this simulation study. In particular, the accuracy and robustness of ionoacoustic range verification for lateral and oblique delivery of high-energy protons to the prostate is examined. Dose maps were computed using GEANT4 Monte Carlo simulations via the TOPAS user interface. Thermoacoustic pulses were propagated using k-Wave software, with initial pressures corresponding to instantaneous dose deposition and piecewise constant maps of tissue properties derived from the planning CT. A database of dose maps with corresponding thermoacoustic emissions and Bragg peak locations, referred to as "control points," were precomputed. Corresponding thermoacoustic emissions were also precomputed. Pulses were recorded at four coplanar locations corresponding to the outer surface of a virtual transrectal array. To model experimental beam delivery, k-Wave results were convolved in time with a Gaussian envelope to account for noninstantaneous proton delivery by a synchrocyclotron. Thermoacoustic pulses were bandlimited below 150 kHz, and amplitudes were directly proportional to charge delivered. To test robustness of our method, white noise was added. Range was estimated in a two-step process. The first step obtained a preliminary range estimate by one-way beamforming. The second step was taken using data corresponding to the "control point" nearest to the preliminary range estimate. For each receiver, the time of flight difference, ∆t, between the measured and control thermoacoustic signals were accurately estimated by applying the Fourier shift theorem. Receiver-Bragg peak distance was then estimated by adding vs ∆t to the known distance of the control point, where vs is soundspeed. A linear system of equations based upon all receiver locations and distances was solved to recover the Bragg peak location. All simulations were performed relative to the planning CT. Because ultrasound (US) images were not available, results were overlaid onto the planning CT. Beamformed estimates from noise-free data tracked all beam locations within 1 cm. Final estimates for oblique and lateral beams were accurate to within 1.0 and 1.6 mm respectively. Average errors of final range estimates for oblique beams from data with SNR = 0 dB were no greater than 2.0 mm. Ionoacoustic range verification may improve current practice. Ionoacoustic range estimates can be inherently co-registered to ultrasound images of underlying anatomy. To ensure estimates are robust in clinical practice, dose maps based upon the planning CT should be overlaid onto ultrasound volumes acquired at time of treatment and acoustic simulations re-computed to provide a database of control points and corresponding thermoacoustic emissions. Computation times for beamformed estimates are already fast enough for online range verification, but are not accurate enough for a measurement aperture limited to the surface of a transrectal ultrasound probe. Accelerated acoustic simulations will be required to enable online two-stage correction, but offline calculation is already suitable for adaptive planning.