Abstract. In this paper, we explore different prognostic methods to account for skin sea surface temperature diurnal variations in a coupled ocean–atmosphere regional model of the Mediterranean Sea. Our aim is to characterise the sensitivity of the considered methods with respect to the underlying assumption of how the solar radiation shapes the warm layer of the ocean. All existing prognostic methods truncate solar transmission coefficient at a warm-layer reference depth that is constant in space and time; instead, we implement a new scheme where this latter is estimated from a chlorophyll dataset as the e-folding depth of solar transmission, which thus allows it to vary in space and time depending on seawater's transparency conditions. Comparison against satellite data shows that our new scheme, compared to the one already implemented within the ocean model, improves the spatially averaged diurnal signal, especially during winter, and the seasonally averaged one in spring and autumn, while showing a monthly basin-wide averaged bias smaller than 0.1 K year-round. In April, when most of the drifters' measurements are available, the new scheme mitigates the bias during nighttime, keeping it positive but smaller than 0.12 K during the rest of the monthly averaged day. The new scheme implemented within the ocean model improves the old one by about 0.1 K, particularly during June. All the methods considered here showed differences with respect to objectively analysed profiles confined between 0.5 K during winter and 1 K in summer for both the eastern and the western Mediterranean regions, especially over the uppermost 60 m. The new scheme reduces the RMSE on the top 15 m in the central Mediterranean for summertime months compared to the scheme already implemented within the ocean model. Overall, the surface net total heat flux shows that the use of a skin sea surface temperature (SST) parameterisation brings the budget about 1.5 W m−2 closer to zero on an annual basis, despite all simulations showing an annual net heat loss from the ocean to the atmosphere. Our “chlorophyll-interactive” method proved to be an effective enhancement of existing methods, its strength relying on an improved physical consistency with the solar extinction implemented in the ocean component.