The three-dimensional, algebraically growing instability of a Blasius boundary layer is studied in the nonlinear regime, employing a nonparallel model based on boundary layer scalings. Adjoint-based optimization is used to determine the "optimal" steady leading-edge excitation that provides the maximum energy growth for a given initial energy. Like in the linear case, the largest transient growth is found for inlet streamwise vortices, that yield streamwise streaks downstream. Two different definitions of growth are employed, providing qualitatively similar results, although the spanwise wavenumbers of optimal growth differ by up to 20% in the two cases. The wavelength of the most amplified optimal disturbance increases with the initial amplitude. For large input amplitudes, significant deformations of the mean velocity field are found; in such cases it is reasonable to expect that nonlinear streaks may break down through a secondary instability. 2005 Elsevier SAS. All rights reserved.