We report new descriptions for the (probability) distributions of hourly maximum crest and wave height of water surface gravity waves for intermediate water depths. Estimated distributions are based on analysis of laboratory-scale measurements at the DHI wave basin. For a given sea state, the distribution of both hourly maximum crest and hourly maximum wave height, normalised by sea state significant wave height, is found to follow a generalised extreme value (GEV) distribution. Variation of the three parameters of the GEV distribution across sea states, is expressed in terms of a response surface model as a function of non-dimensional sea state Ursell number and wave steepness, and wave directional spreading angle. For inference, conventional Monte Carlo wave basin measurements are supplemented with measurements selected by means of a novel “pre-selection” sampling scheme using numerical simulations. This scheme effectively guarantees that extreme events from tails of distributions are produced, and reduces uncertainties associated with the estimated distributions. Estimation is performed using Bayesian inference, allowing uncertainties to be quantified, and providing estimates of posterior predictive tail distributions for sea states with arbitrary characteristics within the domain of sea state characteristics covered by the model.