Energy value of phytoplankton regulates the growth of higher trophic species, affecting the tropic balance and sustainability of marine food webs. Therefore, developing our capability to estimate and monitor, on a global scale, the concentrations of macromolecules that determine phytoplankton energy value, would be invaluable. Reported here are the first estimates of carbohydrate, protein, lipid, and overall energy value of phytoplankton in the world oceans, using ocean-colour data from satellites. The estimates are based on a novel bio-optical method that utilises satellite-derived bio-optical fingerprints of living phytoplankton combined with allometric relationships between phytoplankton cells and cellular macromolecular contents. The annually averaged phytoplankton energy value, per cubic metre of sub-surface ocean, varied from less than 0.1 kJ in subtropical gyres, to 0.5-1.0 kJ in parts of the equatorial, northern and southern latitudes, and rising to >10 kJ in certain coastal and optically complex waters. The annually averaged global stocks of carbohydrate, protein and lipid were 0.044, 0.17 and 0.108 gigatonnes, respectively, with monthly stocks highest in September and lowest in June, over 1997-2013. The fractional contributions of phytoplankton size classes e.g., picoplankton, nanoplankton and microplankton to surface concentrations and global stocks of macromolecules varied considerably across marine biomes classified as Longhurst provinces. Among these provinces, the highest annually averaged surface concentrations of carbohydrate, protein, and lipid were in North-East Atlantic Coastal Shelves, whereas, the lowest concentration of carbohydrate or lipid were in North Atlantic Tropical Gyral, and that of protein was in North Pacific Subtropical Gyre West. The regional accuracy of the estimates and their sensitivity to satellite inputs are quantified from the bio-optical model, which show promise for possible operational monitoring of phytoplankton energy value from satellite ocean colour. Adequate in situ measurements of macromolecules and improved retrievals of inherent optical properties from high-resolution satellite images, would be required to validate these estimates at local sites, and to further improve their accuracy in the world oceans.