Wave description is affected by several uncertainties, with sampling variability due to limited number of observations being one of them. Ideally, temporal/spatial wave registrations should be as large as possible to eliminate this uncertainty. This is difficult to reach in nature, where stationarity of sea states is an issue, but it can in principle be obtained in laboratory tests and numerical simulations, where initial wave conditions can be kept constant and intrinsic variability can be accounted for by changing random seeds for each run. Using linear, second-order, and third-order unidirectional numerical simulations, we compare temporal and spatial statistics of selected wave parameters and show how sampling variability affects their estimators. The JONSWAP spectrum with gamma peakedness parameters γ = 1, 3.3, and 6 is used in the analysis. The third-order wave data are simulated by a numerical solver based on the higher-order spectral method which includes the leading-order nonlinear dynamical effects. Field data support the analysis. We demonstrate that the nonlinear wave field including dynamical effects is more sensitive to sampling variability than the second-order and linear ones. Furthermore, we show that the mean values of temporal and spatial wave parameters can be equal if the number of simulations is sufficiently large. Consequences for design work are discussed.