Can upwelling regions be potential thermal refugia for marine fishes during climate warming?

Species are expected to migrate to higher latitudes as warming intensifies due to anthropogenic climate change since physiological mechanisms have been adapted to maximize fitness under specific temperatures. However, literature suggests that upwellings could act as thermal refugia under climate warming protecting marine ecosystem diversity. This research aimed to predict the effects of climate warming on commercial and non-commercial fish species reported in official Mexican documents (>200 species) based on their thermal niche to observe if upwellings can act as potential thermal refugia. Present (2000–2014) and Representative Concentration Pathway (6.0 and 8.5) scenarios (2040–2050 and 2090–2100) have been considered for this work. Current and future suitability patterns, species distribution, richness, and turnover were calculated using the minimum volume ellipsoids as algorithm. The results in this study highlight that beyond migration to higher latitudes, upwelling regions could protect marine fishes, although the mechanism differed between the innate characteristics of upwellings. Most modeled species (primarily tropical fishes) found refuge in the tropical upwelling in Northern Yucatan. However, the highest warming scenario overwhelmed this region. In contrast, the Baja California region lies within the Eastern Boundary Upwelling Systems. While the area experiences an increase in suitability, the northern regions have a higher upwelling intensity acting as environmental barriers for many tropical species. Conversely, in the southern regions where upwelling is weaker, species tend to congregate and persist even during elevated warming, according to the turnover analysis. These findings suggest that tropicalization in higher latitudes may not be as straightforward as previously assumed. Nevertheless, climate change affects numerous ecosystem features, such as trophic relationships, phenology, and other environmental variables not considered here. In addition, uncertainty still exists about the assumption of increasing intensity of upwelling systems.