This work highlights the effects of nonlinearity on the crest shape of extreme directional water wave events in infinite and finite water depths. HOS-ocean, a High-Order Spectral code developed for solving the deterministic propagation of nonlinear wavefields in open ocean [Ducrozet et al., “HOS-ocean: Open-source solver for nonlinear waves in open ocean based on High-Order Spectral method,” Comput. Phys. Commun. 203, 245–254 (2016)], is used for comparisons with linear simulations. Nonlinear wave harmonics were extracted combining and extending previous methods [Fitzgerald et al., “Phase manipulation and the harmonic components of ringing forces on a surface-piercing column,” Proc. R. Soc. A 470, 20130847 (2014) and Hann et al., “A new set of focused wave linear combinations to extract non-linear wave harmonics,” in Proceedings of 29th International Workshop on Water Waves and Floating Bodies, 2014], showing the contribution of various orders of nonlinearity. The work shows that “walls of water,” [Adcock et al., “Nonlinear dynamics of wave-groups in random seas: Unexpected walls of water in the open ocean,” Proc. R. Soc. A 471, 20150660 (2015) and Adcock et al., “On the shape of large wave-groups on deep water—The influence of bandwidth and spreading,” Phys. Fluids 28, 106601 (2016)], a significant increase in the crest width, are also formulated in finite water depths. Mainly responsible for this is the first order harmonic, with significant energy transfers away from components propagating at an angle, which are more pronounced in finite water depths. As nonlinearity increases, higher order harmonics (>second) become, as expected, more significant. In deep water, rapid energy transfers lead to a broadening of the spectrum which is reflected in the increase in the first order amplitude sum and hence increased crest elevations at the point of the large event. However, in finite water depths, energy transfers are gradually leading to a narrowing of the spectrum around the spectral peak and to a general broadening of the spectrum at the higher harmonics; reflected in the significant decrease in the first order amplitude sum and hence a decrease in crest elevations at the point of the large event, also aided by the larger increase in the out-of-phase difference terms (zeroth harmonic).