AbstractA numerical model with a vortex force formalism is used to study the role of wind waves in the momentum budget and subtidal exchange of a shallow coastal plain estuary, Delaware Bay. Wave height and age in the bay have a spatial distribution that is controlled by bathymetry and fetch, with implications for the surface drag coefficient in young, underdeveloped seas. Inclusion of waves in the model leads to increases in the surface drag coefficient by up to 30% with respect to parameterizations in which surface drag is only a function of wind speed, in agreement with recent observations of air‐sea fluxes in estuaries. The model was modified to prevent whitecapping wave dissipation from generating breaking forces since that contribution is integrally equivalent to the wind stress. The proposed adjustment is consistent with previous studies of wave‐induced nearshore currents and with additional parameterizations for breaking forces in the model. The mean momentum balance during a simulated wind event was mainly between the pressure gradient force and surface stress, with negligible contributions by vortex, wave breaking (i.e., depth‐induced), and Stokes‐Coriolis forces. Modeled scenarios with realistic Delaware bathymetry suggest that the subtidal bay‐ocean exchange at storm time scales is sensitive to wave‐induced surface drag coefficient, wind direction, and mass transport due to the Stokes drift. Results herein are applicable to shallow coastal systems where the typical wave field is young (i.e., wind seas) and modulated by bathymetry.