The second harmonic of a Rayleigh wave passing through a minute surface crack has been numerically analyzed by semi-explicit FEM including special elements which account for a nonlinear stress-strain relation at crack surfaces. Minute cracks perpendicular to a free, flat surface close under compressive stress when width of the crack opening is less than the longitudinal amplitude of the Rayleigh wave. Thereafter, compressive and shear stresses are partially transmitted through the closed cracks, whereas tensile and shear stresses are not transmitted through cracks that remain open. This leads to marked nonlinear ultrasonic response. Calculation was performed for an aluminum block having a surface crack. The transverse component of the Rayleigh wave propagating through the cracks shows distorted waveforms, making the second harmonic amplitude clearly noticeable. In an experiment, the second harmonic component of the leaky Rayleigh wave was detected for a simple crack model consisting of two aluminum blocks, by use of a PVDF line-focused transducer. The results of the experiment show that the second harmonic amplitude is a second-order function of the fundamental wave amplitude, and is more pronounced for low compressive stress applied to close the crack surfaces.