Oceanic net primary production forms the foundation of marine ecosystems. Understanding the impact of climate change on primary production is therefore critical and we rely on Earth System Models to project future changes. Stemming from their use of different physical dynamics and biogeochemical processes, these models yield a large spread in long-term projections of change on both the global and regional scale. Here we review the key physical processes and biogeochemical parameterizations that influence the estimation of primary production in Earth System Models and synthesize the available projections of productivity in the subarctic regions of the North Atlantic. The key processes and modelling issues we focus on are mixed layer depth dynamics, model resolution and the complexity and parameterization of biogeochemistry. From the model mean of five CMIP6 models, we found a large increase in PP in areas where the sea ice retreats throughout the 21st century. Stronger stratification and declining MLD in the Nordic Seas, caused by sea ice loss and regional freshening, reduce the vertical flux of nutrients into the photic zone. Following the synthesis of the primary production among the CMIP6 models, we recommend a number of measures: constraining model hindcasts through the assimilation of high-quality long-term observational records to improve physical and biogeochemical parameterizations in models, developing better parameterizations for the sub-grid scale processes, enhancing the model resolution, downscaling and multi-model comparison exercises for improved regional projections of primary production.publishedVersio
