Oceanic mesoscale eddies are associated with large thermodynamic anomalies, yet so far they are most commonly studied in terms of surface temperature and in the sense of composite mean. Here we employ an objective eddy identification and tracking algorithm together with a novel matching and filling procedure to more thoroughly examine eddy-induced thermodynamic anomalies in the North Pacific, their relationship with eddy amplitude (SSH), and the percentage of variability they explain on various timescales from submonthly to interannual. The thermodynamic anomalies are investigated in terms of sea surface temperature (SST), isothermal layer depth (ITD), and upper ocean heat content (HCT). Most eddies are weak in amplitude and are associated with small thermodynamic anomalies. In the sense of composite mean, anticyclonic eddies are generally warm eddies with deeper isothermal layer and larger heat content, and the reverse is true for cyclonic eddies. A small fraction of eddies, most probably subsurface eddies, exhibits the opposite polarities. Linear relationships with eddy amplitude are found for each of the thermodynamic parameters but with different level of scatter and seasonality. HCT-amplitude relation scatters the least and has the smallest seasonal difference, ITD-amplitude relation has the largest scatter and seasonality, while SST-amplitude relation is in between. For the Kuroshio and Oyashio Extension region, the most eddy-rich region in the North Pacific, eddies are responsible for over 50% of the total SSH variability up to the intra-seasonal scale, and ITD and HCT variability up to interannual. Eddy-induced SST variability is the highest along the Oyashio Extension Front on the order of 40–60% on submonthly scales. These results highlight the role of mesoscale eddies in ocean thermodynamic variability and in air-sea interaction.