Artificial boundaries are necessary in seismic numerical modeling to constrain the computational limits. Special treatment is necessary to reduce spurious reflections from these artificial boundaries. For near-surface seismic modeling, boundary processing is particularly important due to the presence of surface waves. The Perfectly Matched Layer (PML) method is an efficient way of attenuating boundary reflections of body waves and surface waves.