Subseasonal surface wind variability strongly impacts the annual mean and subseasonal turbulent atmospheric surface fluxes. However, the impacts of subseasonal wind variability on the ocean are not fully understood. Here, we quantify the sensitivity of the ocean surface stress ( ), buoyancy flux (B), and mixed layer depth (MLD) to subseasonal wind variability in both a one-dimensional (1-D) vertical column model and a three-dimensional (3-D) global mesoscale-resolving ocean/sea ice model. The winds are smoothed by time filtering the pseudo-stresses, so the mean stress is approximately unchanged, and some important surface flux feedbacks are retained. The 1-D results quantify the sensitivities to wind variability at different time scales from 120 days to 1 day at a few sites. The 3-D results quantify the sensitivities to wind variability shorter than 60 days at all locations, and comparisons between 1-D and 3-D results highlight the importance of 3-D ocean dynamics. Globally, subseasonal winds explain virtually all of subseasonal variance, about half of subseasonal B variance but only a quarter of subseasonal MLD variance. Subseasonal winds also explain about a fifth of the annual mean MLD and a similar and spatially correlated fraction of the mean friction velocity, u * = √ | |∕ sw where sw is the density of seawater. Hence, the subseasonal MLD variance is relatively insensitive to subseasonal winds despite their strong impact on local B and variability, but the mean MLD is relatively sensitive to subseasonal winds to the extent that they modify the mean u * , and both of these sensitivities are modified by 3-D ocean dynamics. Key Points: • We quantify the global impact of subseasonal winds on ocean surface fluxes and mixed-layer depth (MLD) in a mesoscale-resolving model • Globally, subseasonal winds are responsible for about a fifth of the annual average MLD and about a quarter of the subseasonal MLD variance • The increase in the mean MLD with subseasonal winds is caused by a higher friction velocity, but other factors modify the sensitivity Supporting Information: • Supporting Information S1A few previous studies have looked into the sensitivity of the climatological MLD to subseasonal atmospheric variability in 3-D global ocean circulation models by explicitly modifying the wind forcing in ocean models using grids that do not fully resolve mesoscales. For example, Lee and Liu (2005) explore the impact of diurnal winds on the MLD and sea surface temperature using an eddy-parameterizing (nominal 1 • resolution) global ocean model. Their model is forced with daily averaged fluxes, except for the wind stresses, which are averaged over 12 hr or subsampled every 24 hr. The annual and zonal mean MLD increases by up to a maximum of about 10 m with 12-hr winds compared to subsampled 24-hr winds, and the response of the MLD is largely attributable to enhanced vertical mixing by high-frequency winds at middle-to-high latitudes. More recently, Condron and Renfrew (2013) parameterize the effects of unresolved mesoscale a...