Impacts of Paleoecology on the TEX86 Sea Surface Temperature Proxy in the Pliocene-Pleistocene Mediterranean Sea

The TEX86 proxy, based on the distribution of isoprenoid glycerol dialkyl glycerol tetraethers (iGDGTs) from planktonic Thaumarchaeota, is widely used to reconstruct sea surface temperature (SST). Recent observations of species-specific and regionally dependent TEX86-SST relationships in cultures and the modern ocean raise the question of whether nonthermal factors may have impacted TEX86 paleorecords. Here we evaluate the effects of ecological changes on TEX86 using one Pliocene and two Pleistocene sapropels from the Mediterranean Sea. We find that TEX86-derived SSTs deviate from U-37(k) -derived SSTs before, during, and after each sapropel event. U-37(k)-derived SSTs vary by less than 6 degrees C, while TEX86-derived SSTs vary by up to 15 degrees C within a single record. Compound-specific carbon isotope compositions indicate minimal confounding influence on TEX86 from exogenous sources. Some of the variation can be accounted for by changes in nitrogen cycling intensity affecting thaumarchaeal iGDGT biosynthesis, as demonstrated by an inverse relationship between TEX86 and delta N-15(TN). TEX86-derived SSTs also consistently show warm anomalies in the Pleistocene, while the Pliocene samples exhibit both warmer and cooler relative offsets. These anomalies result from systematic differences between Plio-Pleistocene iGDGT distributions and both modern Mediterranean and modern, globally distributed core top samples. Through characteristic GDGT distributions, we suggest the existence of three distinct endemic populations of Thaumarchaeota in the Pliocene, Pleistocene, and modern Mediterranean Sea, respectively. Importantly, these communities prevailed during both sapropel and oligotrophic conditions. Our results demonstrate that ecological and community-specific effects must be considered when applying the TEX86 proxy to paleorecords.