We investigate the nonlinear optical response of a noble metal surface. We derive the components of the third-order nonlinear susceptibility and determine an absolute value of c (3) ≈ 0.2 nm 2 V −2 , a value that is more than two orders of magnitude larger than the values found for typical nonlinear laser crystals. Using nonlinear four-wave mixing (4WM) with incident laser pulses of frequencies u 1 and u 2 , we generate fields oscillating at the nonlinear frequency u 4WM = 2u 1 − u 2 . We identify and discuss three distinct regimes: (i) a regime where the 4WM field is propagating, (ii) a regime where it is evanescent, and (iii) a regime where the nonlinear response couples to surface plasmon polaritons.