The Princeton Ocean Model for Tropical Cyclones (POM-TC), a version of the three-dimensional primitive equation numerical ocean model known as the Princeton Ocean Model, was the ocean component of NOAA’s operational Hurricane Weather Research and Forecast Model (HWRF) from 2007 to 2013. The coupled HWRF–POM-TC system facilitates accurate tropical cyclone intensity forecasts through proper simulation of the evolving SST field under simulated tropical cyclones. In this study, the 2013 operational version of HWRF is used to analyze the POM-TC ocean temperature response in retrospective HWRF–POM-TC forecasts of Atlantic Hurricanes Earl (2010), Igor (2010), Irene (2011), Isaac (2012), and Leslie (2012) against remotely sensed and in situ SST and subsurface ocean temperature observations. The model generally underestimates the hurricane-induced upper-ocean cooling, particularly far from the storm track, as well as the upwelling and downwelling oscillation in the cold wake, compared with observations. Nonetheless, the timing of the model SST cooling is generally accurate (after accounting for along-track timing errors), and the ocean model’s vertical temperature structure is generally in good agreement with observed temperature profiles from airborne expendable bathythermographs.