Abstract. Variations in the atmospheric oxidative capacity, largely determined by variations in the hydroxyl radical (OH), form a key uncertainty in many greenhouse and other pollutant budgets, such as that of methane (CH4). Methyl chloroform (MCF) is an often-adopted tracer to indirectly put observational constraints on large-scale variations in OH. We investigated the budget of MCF in a 4DVAR inversion using the atmospheric transport model TM5, for the period 1998–2018, with the objective to derive information on large-scale, interannual variations in atmospheric OH concentrations. While our main inversion did not fully converge, we did derive interannual variations in the global oxidation of MCF that bring simulated mole fractions of MCF within 1 %–2 % of the assimilated observations from the NOAA-GMD surface network at most sites. Additionally, the posterior simulations better reproduce aircraft observations used for independent validation compared to the prior simulations. The derived OH variations showed robustness with respect to the prior MCF emissions and the prior OH distribution over the 1998 to 2008 period. Although we find a rapid 8 % increase in global mean OH concentrations between 2010 and 2012 that quickly declines afterwards, the derived interannual variations were typically small (< 3 %/yr), with no significant long-term trend in global mean OH concentrations. The inverse system found strong adjustments to the latitudinal distribution of OH, relative to widely used prior distributions, with systematic increases in tropical and decreases in extra-tropical OH concentrations (both up to 30 %). These spatial adjustments were driven by intrahemispheric biases in simulated MCF mole fractions, which have not been identified in previous studies. Given the large amplitude of these adjustments, which exceeds spread between literature estimates, and a residual bias in the MCF intrahemispheric gradients, we suggest a reversal in the extratropical ocean sink of MCF in response to declining atmospheric MCF abundance (as hypothesized in Wennberg et al., 2004). This ocean source provides a more realistic explanation for the biases, possibly complementary to adjustments in the OH distribution. We identified significant added value in the use of a 3D transport model, since it implicitly accounts for variable transport and optimizes the observed spatial gradients of MCF, which is not possible in simpler models. However, we also found a trade-off in computational expense and convergence problems. Despite these convergence problems, the derived OH variations do result in an improved match with MCF observations relative to an interannually repeating prior for OH. Therefore, we consider that variations in OH derived from MCF inversions with 3D models can add value to budget studies of long-lived gases like CH4.