©2019. American Geophysical Union. All Rights Reserved. For ice concentrations less than 85%, internal ice stresses in the sea ice pack are small and sea ice is said to be in free drift. The sea ice drift is then the result of a balance between Coriolis acceleration and stresses from the ocean and atmosphere. We investigate sea ice drift using data from individual drifting buoys as well as Arctic-wide gridded fields of wind, sea ice, and ocean velocity. We perform probabilistic inverse modeling of the momentum balance of free-drifting sea ice, implemented to retrieve the Nansen number, scaled Rossby number, and stress turning angles. Since this problem involves a nonlinear, underconstrained system, we used a Monte Carlo guided search scheme—the Neighborhood Algorithm—to seek optimal parameter values for multiple observation points. We retrieve optimal drag coefficients of CA=1.2×10−3 and CO=2.4×10−3 from 10-day averaged Arctic-wide data from July 2014 that agree with the AIDJEX standard, with clear temporal and spatial variations. Inverting daily averaged buoy data give parameters that, while more accurately resolved, suggest that the forward model oversimplifies the physical system at these spatial and temporal scales. Our results show the importance of the correct representation of geostrophic currents. Both atmospheric and oceanic drag coefficients are found to decrease with shorter temporal averaging period, informing the selection of drag coefficient for short timescale climate models.