Experimental data of femtosecond thick-crystal second-harmonic generation shows that when tuning away from phase matching, a dominating narrow spectral peak appears in the second harmonic that can be tuned over 100's of nm by changing the phase-mismatch parameter. Traditional theory explains this as phase matching between a sideband in the broadband pump to its second-harmonic. However, our experiment is conducted under high input intensities and instead shows excellent quantitative agreement with a nonlocal theory describing cascaded quadratic nonlinearities. This theory explains the detuned peak as a nonlocal resonance that arises due to phase-matching between the pump and a detuned second-harmonic frequency, but where in contrast to the traditional theory the pump is assumed dispersion-free. As a soliton is inherently dispersion-free, the agreement between our experiment and the nonlocal theory indirectly proves that we have observed a soliton-induced nonlocal resonance. The soliton exists in the self-defocusing regime of the cascaded nonlinear interaction and in the normal dispersion regime of the crystal, and needs high input intensities to become excited. A common observation in second-harmonic generation (SHG) of broadband laser pulses in thick crystals, is that when a phase mismatch ∆k is imposed, the second harmonic (SH) spectrum is dominated by a spectrally compressed peak that is wavelength-tunable through ∆k [1][2][3][4][5][6][7][8]. Figure 1(a) shows data from an experiment we performed using a thick β-barium borate (BBO) crystal. The input fundamental wave (FW, center frequency ω 1 ) was an intense femtosecond pulse loosely focused and collimated at the crystal entrance to avoid diffraction. The tuning around ∆k = 0 was achieved by rotating the external angle of the crystal. A striking wavelength tunability over 100's of nanometers is possible, and the peak is also strongly compressed compared to the ideal thin-crystal bandwidth (in this case about 40 nm FWHM). The compressed SH peak pertains for large negative tuning angles, while it disappears for large positive tuning angles (here +5• ). The spectral compression is traditionally explained by a phase-matched sidebands theory, which uses the classical result that the SH efficiency ∝ sinc 2 [∆kL/2]: this explains the decreasing bandwidth in a thick crystal, and the frequency dependence of ∆k(ω) = k 2 (ω) − 2k 1 (ω/2) explains how a SH sideband frequency strongly detuned from the degenerate SH frequency ω 2 = 2ω 1 can become phase matched when ∆k(ω 2 ) = 0. Figure 1(c) shows the predicted phasematching wavelengths by the phase-matched sidebands theory, and remarkably it cannot explain the experimental data for large positive tuning angles.Instead an alternative nonlocal theory, shown in Fig. 1(b), predicts "resonance" wavelengths in the nonlocal response function that up to +4• are in excellent agreement with the experimental peak positions. At +5• the phase-matching condition behind this resonance is no longer fulfilled, yielding a broadband, non-...