The boreal winter response of the ocean mixed layer to the Madden-Julian oscillation (MJO) in the IndoPacific region is determined using in situ observations from the Argo profiling float dataset. Composite averages over numerous events reveal that the MJO forces systematic variations in mixed layer depth and temperature throughout the domain. Strong MJO mixed layer depth anomalies (.15 m peak to peak) are observed in the central Indian Ocean and in the far western Pacific Ocean. The strongest mixed layer temperature variations (.0.68C peak to peak) are found in the central Indian Ocean and in the region between northwest Australia and Java. A heat budget analysis is used to evaluate which processes are responsible for mixed layer temperature variations at MJO time scales. Though uncertainties in the heat budget are on the same order as the temperature trend, the analysis nonetheless demonstrates that mixed layer temperature variations associated with the canonical MJO are driven largely by anomalous net surface heat flux. Net heat flux is dominated by anomalies in shortwave and latent heat fluxes, the relative importance of which varies between active and suppressed MJO conditions. Additionally, rapid deepening of the mixed layer in the central Indian Ocean during the onset of active MJO conditions induces significant basin-wide entrainment cooling. In the central equatorial Indian Ocean, MJO-induced variations in mixed layer depth can modulate net surface heat flux, and therefore mixed layer temperature variations, by up to ;40%. This highlights the importance of correctly representing intraseasonal mixed layer depth variations in climate models in order to accurately simulate mixed layer temperature, and thus air-sea interaction, associated with the MJO.