The internal variability of the thermohaline circulation of the Mediterranean Sea is examined under contrasting extreme thermal and mass atmospheric forcing conditions. Three millennium-long numerical simulation experiments were performed under: (a) the current climatology, (b) a strong buoyancy forcing (SBF) scenario due to cold and dry conditions resembling the Younger Dryas event, and (c) a weak buoyancy forcing (WBF) scenario due to S1a sapropel deposition-like conditions (warm and wet). To isolate the inherent variability of the system, independent of interannual atmospheric forcing variability, the latter was defined as a perpetual year pertinent to each experiment. Self-diagnosed heat and salt fluxes, consistent to sea-surface charac...