Marine species and assemblage change foreshadowed by their thermal bias over Early Jurassic warming

A mismatch of species’ thermal preferences to their environment may indicate how they will respond to future climate change. Averaging this mismatch across species may forewarn that some assemblages will undergo greater reorganization, extirpation, and possibly extinction, than others. Here, we examine how regional warming determines species occupancy and assemblage composition of marine bivalves, brachiopods, and gastropods over one-million-year time steps during the Early Jurassic. Thermal bias, the difference between modelled regional temperatures and species’ long-term thermal optima, predicts a gradient of species occupancy response to warming. Species that become extirpated or extinct tend to have cooler temperature preferences than immigrating species, while regionally persisting species fell midway. Larger regional changes in summer seawater temperatures (up to +10 °C) strengthen the relationship between species thermal bias and the response gradient, which is also stronger for brachiopods than for bivalves, while the relationship collapses during severe seawater deoxygenation. At +3 °C regional seawater warming, around 5 % of pre-existing benthic species in a regional assemblage are extirpated, and immigrating species comprise around one-fourth of the new assemblage. Our results validate thermal bias as an indicator of immigration, persistence, extirpation, and extinction of marine benthic species and assemblages under modern-like magnitudes of climate change.