One of the nagging problems which arises in application of discrete solution methods for wave‐propagation calculations is the presence of reflections or wraparound from the boundaries of the numerical mesh. These undesired events eventually override the actual seismic signals which propagate in the modeled region. The solution to avoiding boundary effects used to be to enlarge the numerical mesh, thus delaying the side reflections and wraparound longer than the range of times involved in the modeling. Obviously this solution considerably increases the expense of computation. More recently, nonreflecting boundary conditions were introduced for the finite‐difference method (Clayton and Enquist, 1977; Reynolds, 1978). These boundary conditions are based on replacing the wave equation in the boundary region by one‐way wave equations which do not permit energy to propagate from the boundaries into the numerical mesh. This approach has been relatively successful, except that its effectiveness degrades for event...