Sea surface temperature (SST) data from the Copernicus Marine Environment Monitoring Service are assimilated into a pan-Arctic ice–ocean coupled model using the ensemble-based local singular evolutive interpolated Kalman (LSEIK) filter. This study found that the SST deviation between model hindcasts and independent SST observations is reduced by the assimilation. Compared with model results without data assimilation, the deviation between the model hindcasts and independent SST observations has decreased by up to 0.28degC at the end of summer. The strongest SST improvements are located in the Greenland Sea, the Beaufort Sea, and the Canadian Arctic Archipelago. The SST assimilation also changes the sea ice concentration (SIC). Improvements of the ice concentrations are found in the Canadian Arctic Archipelago, the Beaufort Sea, and the central Arctic basin, while negative effects occur in the west area of the eastern Siberian Sea and the Laptev Sea. Also, sea ice thickness (SIT) benefits from ensemble SST assimilation.A comparison with upward-looking sonar observations reveals that hindcasts of SIT are improved in the Beaufort Sea by assimilating reliable SST observations into light ice areas. This study illustrates the advantages of assimilating SST observations into an ice–ocean coupled model system and suggests that SST assimilation can improve SIT hindcasts regionally during the melting season.