Boundary integral equation methods are combined with a scattering matrix approach to provide an accurate model for acoustic propagation over an arbitrarily shaped boundary deformation joining two uniform waveguide sections of dissimilar depth. The method is illustrated by up- and down-slope propagation over an idealized shallow water shelf model.