In the present work, we study the statistics of wavefields obtained from non-linear phase-resolved simulations. The numerical model used to generate the waves models wave-wave interactions based on the fully nonlinear Zakharov equations. We vary the simulated wavefield's input spectral properties: directional spreading function, Phillips parameter and peak shape parameter. We then investigate the relationships between a wavefield's input spectral properties and its output physical properties via statistical analysis. We investigate surface elevation distribution, wave definition methods in a nonlinear wavefield with a two-dimensional wavenumber, defined waves' distributions, and the occurrence and spacing of large wave events.