Crustal wave guides are usually heterogeneous on many scales, some of which can be investigated for regional phases by one-way approximation methods. The advantage of one-way propagation methods is the greater saving of computing time and memory, often up to several orders, than the full-waveform numerical methods. By slicing a half-space crustal wave guide into a number of slabs perpendicular to the propagation direction, Wu et al. (2000) formulate a wide-angle pseudoscreen method for the SH half-space problem, which leads to a generalized screen propagator (gsp) for simulating Lg propagation. In this article we introduce a broadband constant-coefficient propagator for the SH half-space problem, which accounts for wide angles in large-contrast media while allowing implementation using Fourier transforms alone. Particular attention is paid to the first-order separation-of-variables screen propagator (svsp1) that significantly improves the split-step Fourier (ssf) method for large lateral variations at the cost of one more Fourier transform in each slab. Advancing wave fields by svsp1 is actually a linear interpolation in the wavenumber domain between two split-step terms. We benchmark the ssf, gsp, and svsp1 synthetic seismograms against the full-waveform boundary-element synthetics for flat, belling, and necking crustal wave guides, which shows that the svsp1 method can model the Lg phase and the mantle wave (head wave) quite well in both travel time, energy, and waveform for most common mantle velocity perturbations. These numerical comparisons also demonstrate some limitations (especially in waveform) of the one-way propagation methods to model the Lg code attributed by forward scatterings as a result of lateral irregularities of the Moho.