Ocean-bottom nodes (OBNs) are widely used because of their wide azimuth, long-offset, and low-frequency advantages. However, in the vertical component of the OBN geophone, a significant amount of S-wave induced noise may be recorded. This significantly impacts the quality of the vertical component data and may affect the follow-up merging of dual-sensor data. Some commercially available methods apply normal-moveout correction, which requires velocity data. We avoid this with an adaptive matching method, which relies solely on seismic data. Therefore, we develop a novel adaptive subtraction method for S-wave leakage suppression using the horizontal component as the noise model. Our method effectively handles the nonstationarity of the input seismic data in all time and space directions, mitigating instability caused by manual selection of the smoothing radius. A method is introduced for estimating a global nonstationary smoothing radius using the noise model. Compared with commercial and stationary smoothing methods, our method can better suppress S-wave noise while balancing residual noise and signal leakage more effectively. Synthetic and field data examples demonstrate significant improvement with our method.