A global 1/4° ocean/sea-ice 50-member ensemble simulation is used to disentangle the low-frequency imprints of the atmospherically-forced oceanic variability and of the chaotic intrinsic oceanic variability (IOV) on the large-scale (10° x10°) ocean heat content (OHC) between 1980 and 2010. The IOV explains most of the interannual-to-decadal large-scale OHC variance over substantial fractions of the global ocean area that increase with depth: 9%, 22%, and 31% in the 0-700m, 700-2000m and 2000m-bottom layers, respectively. Such areas concern principally eddy-active regions, mostly found in the Southern Ocean and in western boundary current extensions, but also concern the subtropical gyres at intermediate and deep levels. The oceanic chaos may also induce random multidecadal fluctuations so that large-scale regional OHC trends computed on the 1980-2010 period cannot be unambiguously attributed to the atmospheric forcing in several oceanic basins at various depths. These results are likely to raise detection and attribution issues from real observations.