Marine ecosystems, particularly in high‐latitude regions such as the Arctic, have been significantly affected by human activities and contributions to climate change. Evaluating how fish populations responded to past changes in their environment is helpful for evaluating their future patterns, but is often hindered by the lack of long‐term biological data available. Using otolith increments of Northeast Arctic cod (Gadus morhua) as a proxy for individual growth, we developed a century‐scale biochronology (1924–2014) based on the measurements of 3,894 fish, which revealed significant variations in cod growth over the last 91 years. We combined mixed‐effect modeling and path analysis to relate these growth variations to selected climate, population and fishing‐related factors. Cod growth was negatively related to cod population size and positively related to capelin population size, one of the most important prey items. This suggests that density‐dependent effects are the main source of growth variability due to competition for resources and cannibalism. Growth was also positively correlated with warming sea temperatures but negatively correlated with the Atlantic Multidecadal Oscillation, suggesting contrasting effects of climate warming at different spatial scales. Fishing pressure had a significant but weak negative direct impact on growth. Additionally, path analysis revealed that the selected growth factors were interrelated. Capelin biomass was positively related to sea temperature and negatively influenced by herring biomass, while cod biomass was mainly driven by fishing mortality. Together, these results give a better understanding of how multiple interacting factors have shaped cod growth throughout a century, both directly and indirectly.

Otolith shape analysis provides a robust tool for the discrimination of many fish stocks in the context of fisheries management. However, there has been little research to examine within-stock temporal stability of otolith morphology in relation to changes in the environment and stock conditions. This study investigated the variability of Northeast Arctic (NEA) cod (Gadus morhua) otolith shape between 1933 and 2015, using elliptical Fourier descriptors extracted from archived material of 2968 mature fish. Series of hierarchical multivariate models were developed to relate shape to the identified optimal windows of some environmental drivers. Differences between years accounted for <3% of the observed variation and no significant differences were found between the average cohort shapes. The models not only confirmed that fish growth was the strongest driver of shape differences, but also highlighted effects of temperature and biomass-related variables at different life stages. Extrinsic factors described only a small fraction of the observed variance, which indicates that environmental changes over time likely account for less than the natural inter-individual variability in otolith shape. These results suggest that overall shape remains relatively stable through time within NEA cod, which further contributes towards a consensus on the biological interpretation of shape differences.

The isotopic composition of inorganic carbon in otoliths (δ13Coto) can be a useful tracer of metabolic rates and a method to study ecophysiology in wild fish. We evaluated environmental and physiological sources of δ13Coto variation in Icelandic and Northeast Arctic (NEA) cod (Gadus morhua) over the years 1914–2013. Individual annual growth increments of otoliths formed at age 3 and 8 were micromilled and measured by isotope-ratio mass spectrometry. Simultaneously, all annual increment widths of the otoliths were measured providing a proxy of fish somatic growth. We hypothesized that changes in the physiological state of the organism, reflected by the isotopic composition of otoliths, can affect the growth rate. Using univariate and multivariate mixed-effects models we estimated conditional correlations between carbon isotopic composition and growth of fish at different levels (within individuals, between individuals, and between years), controlling for intrinsic and extrinsic effects on both otolith measurements. δ13Coto was correlated with growth within individuals and between years, which was attributed to the intrinsic effects (fish age or total length). There was no significant correlation between δ13Coto and growth between individuals, which suggests that caution is needed when interpreting δ13Coto signals. We found a significant decrease in δ13Coto through the century which was explained by the oceanic Suess effect-admixture of isotopically light carbon from fossil fuel. We calculated the proportion of the respired carbon in otolith carbonate (Cresp) using carbon isotopic composition in diet and dissolved inorganic carbon of the seawater. This approach allowed us to correct the values for each stock in relation to these two environmental baselines. Cresp was on average 0.275 and 0.295 in Icelandic and NEA stock, respectively. Our results provide an insight into the physiological basis for differences in growth characteristics between these two cod stocks, and how that may vary over time.