It is recognized that in large-area, very-high-frequency capacitively coupled plasma (VHF CCP) reactors, the higher harmonics generated by nonlinear sheath motion can lead to enhanced standing wave excitation. In this work, a self-consistent electromagnetic model, which couples a one-dimensional, radial nonlinear transmission line model with a bulk plasma fluid model, is employed to investigate the nonlinear standing wave excitation in a VHF driven, geometrically asymmetric capacitive argon discharge operated at low pressure. By considering a radially nonuniform plasma density profile (case I ) calculated self-consistently by the nonlinear electromagnetic model and the corresponding radially-averaged, uniform plasma density profile (case II ), we first examine the effect of the plasma density nonuniformity on the propagation of electromagnetic surface waves in a 3 Pa argon discharge driven at 100MHz and 90W. Compared to case II, the higher plasma density at the radial center in case I determines a higher plasma series resonance frequency, yielding stronger high-order harmonic excitations and more significant central peak in the harmonic current density Jz,n and the harmonic electron power absorption pn profiles. Therefore, under the assumption of the radially uniform plasma density in a CCP discharge, the self-excitation of higher harmonics at the radial center should be underestimated. Second, using the self-consistent electromagnetic model, the effect of the rf power on the excitation of nonlinear standing waves is investigated in a 3 Pa argon discharge driven at 100MHz. At a low power of 30W, the discharge is dominated by the first two harmonics. The higher harmonic excitations and the nonlinear standing waves are observed to be enhanced with increasing the rf power, resulting in a more pronounced central peak in the radial profiles of the total electron power absorption density pe, the electron temperature Te, and the electron density ne. For all rf powers, the calculated radial profiles of ne show good agreement with the experimental data obtained by a floating double probe.