Local reflects global: Life stage-dependent changes in the phenology of coastal habitat use by North Sea herring.

Climate warming is affecting the suitability and utilization of coastal habitats by marine fishes around the world. Phenological changes are an important indicator of population responses to climate-induced changes but remain difficult to detect in marine fish populations. The design of large-scale monitoring surveys does not allow fine-grained temporal inference of population responses, while the responses of ecologically and economically important species groups such as small pelagic fish are particularly sensitive to temporal resolution. Here, we use the longest, highest resolution time series of species composition and abundance of marine fishes in northern Europe to detect possible phenological shifts in the small pelagic North Sea herring. We detect a clear forward temporal shift in the phenology of nearshore habitat use by small juvenile North Sea herring. This forward shift might be linked to changes in water temperatures in the North Sea. We next assessed the robustness of the effects we found with respect to monitoring design. We find that reducing the temporal resolution of our data to reflect the resolution typical of larger surveys makes it difficult to detect phenological shifts and drastically reduces the effect sizes of environmental covariates such as seawater temperature. Our study therefore shows how local, long-term, high-resolution time series of fish catches are essential to understand the general phenological responses of marine fishes to climate warming and to define ecological indicators of system-level changes.

Acclimation capacity to global warming of amphibians and freshwater fishes: Drivers, patterns, and data limitations.

Amphibians and fishes play a central role in shaping the structure and function of freshwater environments. These organisms have a limited capacity to disperse across different habitats and the thermal buffer offered by freshwater systems is small. Understanding determinants and patterns of their physiological sensitivity across life history is, therefore, imperative to predicting the impacts of climate change in freshwater systems. Based on a systematic literature review including 345 experiments with 998 estimates on 96 amphibian (Anura/Caudata) and 93 freshwater fish species (Teleostei), we conducted a quantitative synthesis to explore phylogenetic, ontogenetic, and biogeographic (thermal adaptation) patterns in upper thermal tolerance (CTmax) and thermal acclimation capacity (acclimation response ratio, ARR) as well as the influence of the methodology used to assess these thermal traits using a conditional inference tree analysis. We found globally consistent patterns in CTmax and ARR, with phylogeny (taxa/order), experimental methodology, climatic origin, and life stage as significant determinants of thermal traits. The analysis demonstrated that CTmax does not primarily depend on the climatic origin but on experimental acclimation temperature and duration, and life stage. Higher acclimation temperatures and longer acclimation times led to higher CTmax values, whereby Anuran larvae revealed a higher CTmax than older life stages. The ARR of freshwater fishes was more than twice that of amphibians. Differences in ARR between life stages were not significant. In addition to phylogenetic differences, we found that ARR also depended on acclimation duration, ramping rate, and adaptation to local temperature variability. However, the amount of data on early life stages is too small, methodologically inconsistent, and phylogenetically unbalanced to identify potential life cycle bottlenecks in thermal traits. We, therefore, propose methods to improve the robustness and comparability of CTmax/ARR data across species and life stages, which is crucial for the conservation of freshwater biodiversity under climate change.

Climate risk to European fisheries and coastal communities.

With the majority of the global human population living in coastal regions, correctly characterizing the climate risk that ocean-dependent communities and businesses are exposed to is key to prioritizing the finite resources available to support adaptation. We apply a climate risk analysis across the European fisheries sector to identify the most at-risk fishing fleets and coastal regions and then link the two analyses together. We employ an approach combining biological traits with physiological metrics to differentiate climate hazards between 556 populations of fish and use these to assess the relative climate risk for 380 fishing fleets and 105 coastal regions in Europe. Countries in southeast Europe as well as the United Kingdom have the highest risks to both fishing fleets and coastal regions overall, while in other countries, the risk-profile is greater at either the fleet level or at the regional level. European fisheries face a diversity of challenges posed by climate change; climate adaptation, therefore, needs to be tailored to each country, region, and fleet's specific situation. Our analysis supports this process by highlighting where and what adaptation measures might be needed and informing where policy and business responses could have the greatest impact.

The costs and trade-offs of optimal foraging in marine fish larvae.

In a warming world, both the metabolic rates of ectotherm predators and the phenology of their prey organisms is subject to change. Knowledge on how intrinsic and extrinsic factors govern predator-prey interactions is essential in order to understand how the environment regulates the vital rates of consumers. Controlled experiments, however, simultaneously testing behavioural and growth responses of the larvae of fish and other ectotherm organisms in different feeding regimes are scarce. Prey size (PS) selection was determined for young Atlantic herring Clupea harengus L. larvae offered 100- to 850-µm copepods Acartia tonsa at five different concentrations. In separate, 4- (13°C) or 7-day (7°C) trials, the effect of prey size on larval foraging behaviour, specific growth rate (SGR) and biochemical condition (RNA:DNA, RD, a proxy for individual instantaneous growth) was tested. Preferred (selected) PS was similar at all prey concentrations but increased from 3% to 5% predator length with increasing larval size. At various temperatures, dome-shaped relationships existed between PS and larval RD (and accordingly SGR). Compensatory changes in foraging behaviour (pause and feeding strike frequencies) existed but were not adequate to maintain positive SGR when available prey were substantially smaller than those preferred by larvae. A physiology-based model predicted that larvae depended more heavily on optimal prey sizes at the colder versus warmer temperature to grow well and that the profitable prey niche breadth (the range in prey sizes in which growth was positive) increased at warmer temperatures. Seemingly subtle match-mismatch dynamics between ectotherm predators and their preferred prey size can have large, temperature-dependent consequences for rates of growth and likely survival of the predator. To the best knowledge, this was the first study to directly quantify the "costs and trade-offs" of optimal foraging in marine fish larvae.

Effects of ocean acidification over successive generations decrease resilience of larval European sea bass to ocean acidification and warming but juveniles could benefit from higher temperatures in the NE Atlantic.

European sea bass (Dicentrarchus labrax) is a large, economically important fish species with a long generation time whose long-term resilience to ocean acidification (OA) and warming (OW) is not clear. We incubated sea bass from Brittany (France) for two generations (>5 years in total) under ambient and predicted OA conditions (PCO2: 650 and 1700 µatm) crossed with ambient and predicted OW conditions in F1 (temperature: 15-18°C and 20-23°C) to investigate the effects of climate change on larval and juvenile growth and metabolic rate. We found that in F1, OA as a single stressor at ambient temperature did not affect larval or juvenile growth and OW increased developmental time and growth rate, but OAW decreased larval size at metamorphosis. Larval routine and juvenile standard metabolic rate were significantly lower in cold compared with warm conditioned fish and also lower in F0 compared with F1 fish. We did not find any effect of OA as a single stressor on metabolic rate. Juvenile PO2,crit was not affected by OA or OAW in both generations. We discuss the potential underlying mechanisms resulting in the resilience of F0 and F1 larvae and juveniles to OA and in the beneficial effects of OW on F1 larval growth and metabolic rate, but contrastingly in the vulnerability of F1, but not F0 larvae to OAW. With regard to the ecological perspective, we conclude that recruitment of larvae and early juveniles to nursery areas might decrease under OAW conditions but individuals reaching juvenile phase might benefit from increased performance at higher temperatures.

Disentangling seasonal from maternal effects on egg characteristics in western Baltic spring-spawning herring Clupea harengus.

In marine fishes, the timing of spawning determines the environment offspring will face and, therefore, the chances of early life stage survival. Different waves of Atlantic herring Clupea harengus spawn throughout spring in the western Baltic Sea, and the survival of offspring from early in the season has been low in the most recent decade. The authors assessed changes in egg traits from early, middle and late phases of the spawning season to examine whether seasonal and/or maternal effects influenced embryo survival. At each phase, fertilized eggs of six to eight females were incubated at two temperatures (7 and 13°C), and egg size, fertilization success, mortality and time to hatch were recorded. A compilation of data from 2017 to 2020 spawning seasons indicated that mean total length of females decreased with progression of the season and increasing in situ water temperature. For the sub-set of females used in the laboratory study, early spawners were 7.6% larger and produced 14.2% larger eggs than late-spring spawners. Fertilization success was consistently high (>90%), and mortality to hatch was low (<3%). Neither the former nor latter were influenced by season, but both were influenced by maternity. This significant female effect was, however, not related to any maternal trait measured here (total length, Fulton's condition factor or age). There was no maternal effect on development rate at 7 or 13°C. The results suggest that intrinsic differences among females or among spawning waves are unlikely to markedly contribute to the poor survival observed for progeny from early in the season in this population and point towards other extrinsic factors or processes acting on eggs or early larval stages (e.g., seasonal match-mismatch dynamics with prey) as more likely causes of mortality.

Embryonic development and effect of temperature on larval growth of the Peruvian anchovy Engraulis ringens.

Understanding aspects of the biology of early life stages of marine fish is critical if one hopes to reveal the factors and processes that impact the survival and recruitment (year class) strength. The Peruvian anchovy (Engraulis ringens) is a key species in the Humboldt current system, and the present study provides the first description of the embryonic and larval development of this species reared in captivity. Embryonic and early exogenous feeding stages of larvae were illustrated in detail at 18.5°C. Hatching was completed within 42 and 48 h post-fertilization at 18.5 and 14.5°C, respectively. Mean ± 95% C.I. standard length (LS ) at hatch (3.40 ± 0.10 mm at 18.5°C and 2.76 ± 0.34 mm at 14.5°C) was significantly different between the two temperatures. Larval behaviour was assessed at 18.5°C; at the onset of exogenous feeding [3 days post-hatch (dph)], larvae were fed small, motile dinoflagellates, Akashiwo sanguinea. At 7 dph, larvae started to feed almost exclusively on zooplankton (rotifers and Artemia nauplii). Larval activity increased with age, and the first sign of schooling was noted at 31 dph (18.56 mm LS ) at 18.5°C. Temperature had a significant effect on size-at-age, but not on body shape (depth to LS ratio). The size-at-age data for larvae (this study) was used to parameterize a temperature-corrected von Bertalanffy growth function for Peruvian anchovy, the accuracy of which was assessed for juveniles and adults (literature values).

Altered thyroid hormone levels affect the capacity for temperature-induced developmental plasticity in larvae of Rana temporaria and Xenopus laevis.

Anuran larvae show phenotypic plasticity in age and size at metamorphosis as a response to temperature variation. The capacity for temperature-induced developmental plasticity is determined by the thermal adaptation of a population. Multiple factors such as physiological responses to changing environmental conditions, however, might influence this capacity as well. In anuran larvae, thyroid hormone (TH) levels control growth and developmental rate and changes in TH status are a well-known stress response to sub-optimal environmental conditions. We investigated how chemically altered TH levels affect the capacity to exhibit temperature-induced developmental plasticity in larvae of the African clawed frog (Xenopus laevis) and the common frog (Rana temporaria). In both species, TH level influenced growth and developmental rate and modified the capacity for temperature-induced developmental plasticity. High TH levels reduced thermal sensitivity of metamorphic traits up to 57% (R. temporaria) and 36% (X. laevis). Rates of growth and development were more plastic in response to temperature in X. laevis (+30%) than in R. temporaria (+6%). Plasticity in rates of growth and development is beneficial to larvae in heterogeneous habitats as it allows a more rapid transition into the juvenile stage where rates of mortality are lower. Therefore, environmental stressors that increase endogenous TH levels and reduce temperature-dependent plasticity may increase risks and the vulnerability of anuran larvae. As TH status also influences metabolism, future studies should investigate whether reductions in physiological plasticity also increases the vulnerability of tadpoles to global change.

Altered thyroid hormone levels affect body condition at metamorphosis in larvae of Xenopus laevis.

Chemical, physical and biological environmental stressors may affect the endocrine system, such as the thyroid hormone (TH) axis in larval amphibians with consequences for energy partitioning among development, growth and metabolism. We studied the effects of two TH level affecting compounds, exogenous l-thyroxine (T4 ) and sodium perchlorate (SP), on various measures of development and body condition in larvae of the African clawed frog (Xenopus laevis). We calculated the scaled mass index, hepatosomatic index and relative tail muscle mass as body condition indices to estimate fitness. Altered TH levels significantly altered the growth, development, survival and body condition in metamorphic larvae in different directions. While exogeno us T4 reduced growth and accelerated development, SP treatment increased growth but slowed down development. Altered TH levels improved body conditions in both treatments and particularly in larvae of the SP treatment but to the detriment of lower survival rates in both TH level altering treatments. The hepatosomatic index was negatively affected by exogenous T4 , but not by SP treatment indicating a lower lipid reserve in the liver in larvae of T4 treatment. These altered TH levels as caused by several environmental stressors may have an influence on individual fitness across life, as body condition at the onset of metamorphosis determines metamorphic and juvenile survival. Further research is needed to determine synergetic effects of environmental stressors on TH levels and its effects on physiological traits such as metabolic rate.

Patterns of temperature induced developmental plasticity in anuran larvae.

Anurans exhibit plasticity in the timing of metamorphosis and tadpoles show phenotypic plasticity in age and size at metamorphosis as a response to temperature variation. This developmental plasticity to changing thermal conditions is expected to be a primary factor that dictates the vulnerability of amphibians to increasing ambient temperatures such as are predicted in climate change scenarios. We analyzed the patterns of thermal effects on size and age at metamorphosis to investigate whether the intraspecific "temperature-size rule" is applicable over a broad range of anuran species by carrying out a combined analysis based on the data from 25 studies performed on 18 anuran species. Furthermore, we tested whether the thermal background of respective populations impacts the capacity for a plastic response in metamorphic traits. We could confirm this pattern for across-population comparisons. All included populations developed faster and 75% were smaller at the onset of metamorphosis when developmental temperatures were warmer, but the sensitivity of growth and developmental rate to a given temperature change was different. We found that the thermal background of a population influences the sensitivity of metamorphic traits and thus, the capacity for a plastic response in growth and developmental rate. Warm adapted populations were less sensitive to temperature variation indicating a reduced capacity for developmental plasticity and therefore, those species may be more vulnerable to the impacts of climate change. Future studies should include a broader range of rearing temperatures and temperature fluctuations to determine full knowledge of the capacity for developmental plasticity within a species-specific thermal window.

Short-term molecular and physiological responses to heat stress in neritic copepods Acartia tonsa and Eurytemora affinis.

Invertebrates inhabiting shallow water habitats represent particularly appropriate organisms for studying the acclimation potential to environmental stress, since they naturally experience large fluctuations in key abiotic factors such as temperature and salinity. We quantified the biochemical- (mRNA transcripts of 78-kDa glucose-regulated protein (grp78), 70-kDa heat shock protein (hsp70), 90-kDa heat shock protein (hsp90), protein synthesis of HSP70) and organismal- (oxygen consumption rates) level responses to acute heat stress on two neritic copepods (Acartia tonsa and Eurytemora affinis) with special emphasis on the role of short-term acclimation. Transcripts of hsp increased with increasing acute temperature exposure and protein quantities (HSP70) were detectable for 30h. In A. tonsa, HSP70 synthesis was also associated with handling stress. In E. affinis, heat-dependent responses were detected in hsp90, grp78 (mRNA) and HSP70 (protein) expression. Acclimation to a warmer temperature significantly decreased the heat stress response in both species. In A. tonsa, short-term acclimation to heat was not detected at the organismal level via metabolic rate. This study reveals interspecific differences in both the gene expression of stress molecules (e.g. hsp90) as well as the stress factors needed to evoke a stress response (heat vs. handling). We demonstrate that cellular stress markers can be useful measures of short-term thermal acclimation in copepods, which may remain undetected by organismal-level measures.

Projecting effects of climate change on marine systems: is the mean all that matters?

Studies dealing with the effects of changing global temperatures on living organisms typically concentrate on annual mean temperatures. This, however, might not be the best approach in temperate systems with large seasonality where the mean annual temperature is actually not experienced very frequently. The mean annual temperature across a 50-year, daily time series of measurements at Helgoland Roads (54.2° N, 7.9° E) is 10.1°C while seasonal data are characterized by a clear, bimodal distribution; temperatures are around 6°C in winter and 15°C in summer with rapid transitions in spring and autumn. Across those 50 years, the temperature at which growth is maximal for each single bloom event for 115 phytoplankton species (more than 6000 estimates of optimal temperature) mirrors the bimodal distribution of the in situ temperatures. Moreover, independent laboratory data on temperature optima for growth of North Sea organisms yielded similar results: a deviance from the normal distribution, with a gap close to the mean annual temperature, and more optima either above or below this temperature. We conclude that organisms, particularly those that are short-lived, are either adapted to the prevailing winter or summer temperatures in temperate areas and that few species exist with thermal optima within the periods characterized by rapid spring warming and autumn cooling.

Solutions for ecosystem-level protection of ocean systems under climate change.

The Paris Conference of Parties (COP21) agreement renewed momentum for action against climate change, creating the space for solutions for conservation of the ocean addressing two of its largest threats: climate change and ocean acidification (CCOA). Recent arguments that ocean policies disregard a mature conservation research field and that protected areas cannot address climate change may be oversimplistic at this time when dynamic solutions for the management of changing oceans are needed. We propose a novel approach, based on spatial meta-analysis of climate impact models, to improve the positioning of marine protected areas to limit CCOA impacts. We do this by estimating the vulnerability of ocean ecosystems to CCOA in a spatially explicit manner and then co-mapping human activities such as the placement of renewable energy developments and the distribution of marine protected areas. We test this approach in the NE Atlantic considering also how CCOA impacts the base of the food web which supports protected species, an aspect often neglected in conservation studies. We found that, in this case, current regional conservation plans protect areas with low ecosystem-level vulnerability to CCOA, but disregard how species may redistribute to new, suitable and productive habitats. Under current plans, these areas remain open to commercial extraction and other uses. Here, and worldwide, ocean conservation strategies under CCOA must recognize the long-term importance of these habitat refuges, and studies such as this one are needed to identify them. Protecting these areas creates adaptive, climate-ready and ecosystem-level policy options for conservation, suitable for changing oceans.

Thermal impacts on the growth, development and ontogeny of critical swimming speed in Atlantic herring larvae.

Increases in swimming ability have a profound influence on larval fish growth and survival by increasing foraging success, predator avoidance and the ability to favorably influence transport. Understanding how development and environmental factors combine to influence swimming performance in aquatic organisms is particularly important during the transition from viscous to inertial environments. We measured the growth, development and ontogenetic changes in critical swimming speed (Ucrit) in Atlantic herring (Clupea harengus) larvae reared at three temperatures (7, 11, 15°C). Temperature had a significant effect on growth rates (from 0.21 at 7°C to 0.34mm·d(-1) at 15°C), and larval morphology-at-length (increased dry weight (DW), body height and developmental rate at warmer temperatures). Temperature-dependent differences in morphology influenced swimming performance (e.g. the exponential increase in Ucrit with increasing body size was faster at warmer temperatures). Larvae entered the transition to an inertial environment (Reynolds numbers ≥300) at body lengths between 15 (15°C) and 17mm (7°C). Inter-individual differences in Ucrit were not related to nutritional condition (RNA·DNA(-1) or DNA·DW(-1)), but were negatively correlated to length-at-age, suggesting a trade-off between growth rate and locomotor activity. The Ucrit data from this and previously published studies suggest that Atlantic herring pass through four activity phases: 1) yolk-sac (<0.6cm·s(-1)), 2) pre-flexion (0.6-3.0cm·s(-1), temperature effect changes with body size), 3) post-flexion (up to 6-8cm·s(-1), Q10~1.8-2.0), 4) juvenile-adult period (20-170cm·s(-1)).

Statistically downscaled CMIP6 ocean variables for European waters.

Climate change impact studies need climate projections for different scenarios and at scales relevant to planning and management, preferably for a variety of models and realizations to capture the uncertainty in these models. To address current gaps, we statistically downscaled (SD) 3-7 CMIP6 models for five key indicators of marine habitat conditions: temperature, salinity, pH, oxygen, and chlorophyll across European waters for three climate scenarios SSP1-2.6, SSP2-4.5, and SSP5-8.5. Results provide ensemble averages and uncertainty estimates that can serve as input data for projecting the potential success of a range of Nature-based Solutions, including the restoration of habitat-forming species such as seagrass in the Mediterranean and kelp in coastal areas of Portugal and Norway. Evaluation of the ensemble with observations from four European regions (North Sea, Baltic Sea, Bay of Biscay, and Mediterranean Sea) indicates that the SD projections realistically capture the climatological conditions of the historical period 1993-2020. Model skill (Liu-mean efficiency, Pearson correlation) clearly improves for both surface temperature and oxygen across all regions with respect to the original ESMs demonstrating a higher skill for temperature compared to oxygen. Warming is evident across all areas and large differences among scenarios fully emerge from the background uncertainties related to internal variability and model differences in the second half of the century. Scenario-specific differences in acidification significantly emerge from model uncertainty and internal variability leading to distinct trajectories in surface pH starting before mid-century (in some cases starting from present day). Deoxygenation is also present across all domains, but the climate signal was significantly weaker compared to the other two indicators when compared to model uncertainty and internal variability, and the impact of different greenhouse gas trajectories is less distinct. The substantial regional and local heterogeneity in these three abiotic indicators underscores the need for highly spatially resolved physical and biogeochemical projections to understand how climate change may impact marine ecosystems.

Oyster larvae used for ecosystem restoration benefit from increased thermal fluctuation.

A bottleneck in restoring self-sustaining beds of the European oyster (Ostrea edulis) is the successful development and settlement of larvae to bottom habitats. These processes are largely governed by temperature but a mechanistic understanding of larval performance across ecologically relevant temperatures is lacking. We reared larvae at low (20-21 °C) and high (20-24 °C) fluctuating temperatures and applied short-term exposures of larvae to temperatures between 16 and 33 °C to assess vital rates and thermal coping ranges. Larval thermal preference was between 25 and 30 °C for both rearing treatments which corresponded with optimum temperatures for oxygen consumption rates and locomotion. Larvae had 5.5-fold higher settling success, however, when reared at the high compared to the low fluctuating temperatures. Higher mean and periods of increased temperature, as projected in a future climate, may therefore enhance recruitment success of O. edulis in northern European habitats.

Hidden impacts of ocean warming and acidification on biological responses of marine animals revealed through meta-analysis.

Conflicting results remain on the impacts of climate change on marine organisms, hindering our capacity to predict the future state of marine ecosystems. To account for species-specific responses and for the ambiguous relation of most metrics to fitness, we develop a meta-analytical approach based on the deviation of responses from reference values (absolute change) to complement meta-analyses of directional (relative) changes in responses. Using this approach, we evaluate responses of fish and invertebrates to warming and acidification. We find that climate drivers induce directional changes in calcification, survival, and metabolism, and significant deviations in twice as many biological responses, including physiology, reproduction, behavior, and development. Widespread deviations of responses are detected even under moderate intensity levels of warming and acidification, while directional changes are mostly limited to more severe intensity levels. Because such deviations may result in ecological shifts impacting ecosystem structures and processes, our results suggest that climate change will likely have stronger impacts than those previously predicted based on directional changes alone.

Cross-basin and cross-taxa patterns of marine community tropicalization and deborealization in warming European seas.

Ocean warming and acidification, decreases in dissolved oxygen concentrations, and changes in primary production are causing an unprecedented global redistribution of marine life. The identification of underlying ecological processes underpinning marine species turnover, particularly the prevalence of increases of warm-water species or declines of cold-water species, has been recently debated in the context of ocean warming. Here, we track changes in the mean thermal affinity of marine communities across European seas by calculating the Community Temperature Index for 65 biodiversity time series collected over four decades and containing 1,817 species from different communities (zooplankton, coastal benthos, pelagic and demersal invertebrates and fish). We show that most communities and sites have clearly responded to ongoing ocean warming via abundance increases of warm-water species (tropicalization, 54%) and decreases of cold-water species (deborealization, 18%). Tropicalization dominated Atlantic sites compared to semi-enclosed basins such as the Mediterranean and Baltic Seas, probably due to physical barrier constraints to connectivity and species colonization. Semi-enclosed basins appeared to be particularly vulnerable to ocean warming, experiencing the fastest rates of warming and biodiversity loss through deborealization.

Temperature and prey density drive growth and otolith formation of the world's most valuable fish stock.

Peruvian anchovy (Engraulis ringens) represents the largest single-species fishery worldwide. Knowledge on how temperature and prey availability influences growth and age estimation during marine fish early life stages is critical for predicting bottom-up processes impacting stock productivity under changing environmental conditions. We reared Peruvian anchovy larvae at two temperatures (14.5 and 18.5 °C) and prey concentrations [high (HF), and low (LF)] from 6 to 30 days post-hatch (dph) to measure growth rate and examine daily deposition of otolith increments. Peruvian anchovy larvae grew faster at 18.5 °C compared to 14.5 °C. Larvae reared at low prey concentration (18.5-LF) and low temperature (14.5-HF) grew 61 and 35% slower, respectively, than those at high prey and warm temperature (18.5-HF). Age and growth rates of larvae were well depicted in the otolith microstructure of well-fed larvae at 18.5 °C. However, larvae reared at 18.5-LF or 14.5-HF, had only 55 and 49% of the expected number of daily otolith increments. Our results suggest caution when attempting to explore how ocean processes regulate small pelagic stocks, the productivity of which are largely driven by changes in the survival and growth of young larvae.

Ocean acidification and seasonal temperature extremes combine to impair the thermal physiology of a sub-Antarctic fish.

To predict the potential impacts of climate change on marine organisms, it is critical to understand how multiple stressors constrain the physiology and distribution of species. We evaluated the effects of seasonal changes in seawater temperature and near-future ocean acidification (OA) on organismal and sub-organismal traits associated with the thermal performance of Eleginops maclovinus, a sub-Antarctic notothenioid species with economic importance to sport and artisanal fisheries in southern South America. Juveniles were exposed to mean winter and summer sea surface temperatures (4 and 10 °C) at present-day and near-future pCO2 levels (~500 and 1800 µatm). After a month, the Critical Thermal maximum and minimum (CTmax, CTmin) of fish were measured using the Critical Thermal Methodology and the aerobic scope of fish was measured based on the difference between their maximal and standard rates determined from intermittent flow respirometry. Lipid peroxidation and the antioxidant capacity were also quantified to estimate the oxidative damage potentially caused to gill and liver tissue. Although CTmax and CTmin were higher in individuals acclimated to summer versus winter temperatures, the increase in CTmax was minimal in juveniles exposed to the near-future compared to present-day pCO2 levels (there was a significant interaction between temperature and pCO2 on CTmax). The reduction in the thermal tolerance range under summer temperatures and near-future OA conditions was associated with a reduction in the aerobic scope observed at the elevated pCO2 level. Moreover, an oxidative stress condition was detected in the gill and liver tissues. Thus, chronic exposure to OA and the current summer temperatures pose limits to the thermal performance of juvenile E. maclovinus at the organismal and sub-organismal levels, making this species vulnerable to projected climate-driven warming.