The Norwegian Climate Prediction Model version 1 (NorCPM1) is a new research tool for performing climate reanalyses and seasonal-to-decadal climate predictions. It combines the Norwegian Earth System Model version 1 (NorESM1) – which features interactive aerosol–cloud schemes and an isopycnic-coordinate ocean component with biogeochemistry – with anomaly assimilation of sea surface temperature (SST) and $T/S$ -profile observations using the ensemble Kalman filter (EnKF). We describe the Earth system component and the data assimilation (DA) scheme, highlighting implementation of new forcings, bug fixes, retuning and DA innovations. Notably, NorCPM1 uses two anomaly assimilation variants to assess the impact of sea ice initialization and climatological reference period: the first (i1) uses a 1980–2010 reference climatology for computing anomalies and the DA only updates the physical ocean state; the second (i2) uses a 1950–2010 reference climatology and additionally updates the sea ice state via strongly coupled DA of ocean observations. We assess the baseline, reanalysis and prediction performance with output contributed to the Decadal Climate Prediction Project (DCPP) as part of the sixth Coupled Model Intercomparison Project (CMIP6). The NorESM1 simulations exhibit a moderate historical global surface temperature evolution and tropical climate variability characteristics that compare favourably with observations. The climate biases of NorESM1 using CMIP6 external forcings are comparable to, or slightly larger than those of, the original NorESM1 CMIP5 model, with positive biases in Atlantic meridional overturning circulation (AMOC) strength and Arctic sea ice thickness, too-cold subtropical oceans and northern continents, and a too-warm North Atlantic and Southern Ocean. The biases in the assimilation experiments are mostly unchanged, except for a reduced sea ice thickness bias in i2 caused by the assimilation update of sea ice, generally confirming that the anomaly assimilation synchronizes variability without changing the climatology. The i1 and i2 reanalysis/hindcast products overall show comparable performance. The benefits of DA-assisted initialization are seen globally in the first year of the prediction over a range of variables, also in the atmosphere and over land. External forcings are the primary source of multiyear skills, while added benefit from initialization is demonstrated for the subpolar North Atlantic (SPNA) and its extension to the Arctic, and also for temperature over land if the forced signal is removed. Both products show limited success in constraining and predicting unforced surface ocean biogeochemistry variability. However, observational uncertainties and short temporal coverage make biogeochemistry evaluation uncertain, and potential predictability is found to be high. For physical climate prediction, i2 performs marginally better than i1 for a range of variables, especially in the SPNA and in the vicinity of sea ice, with notably improved sea level variability of the Southern Ocean. Despite similar skills, i1 and i2 feature very different drift behaviours, mainly due to their use of different climatologies in DA; i2 exhibits an anomalously strong AMOC that leads to forecast drift with unrealistic warming in the SPNA, whereas i1 exhibits a weaker AMOC that leads to unrealistic cooling. In polar regions, the reduction in climatological ice thickness in i2 causes additional forecast drift as the ice grows back. Posteriori lead-dependent drift correction removes most hindcast differences; applications should therefore benefit from combining the two products. The results confirm that the large-scale ocean circulation exerts strong control on North Atlantic temperature variability, implying predictive potential from better synchronization of circulation variability. Future development will therefore focus on improving the representation of mean state and variability of AMOC and its initialization, in addition to upgrades of the atmospheric component. Other efforts will be directed to refining the anomaly assimilation scheme – to better separate internal and forced signals, to include land and atmosphere initialization and new observational types – and improving biogeochemistry prediction capability. Combined with other systems, NorCPM1 may already contribute to skilful multiyear climate prediction that benefits society.

Increasing surface air temperature is a fundamental characteristic of a warming world. Rising temperatures have potential impacts on human health through heat stress. One heat stress metric is the wet-bulb globe temperature, which takes into consideration the effects of radiation, humidity, and wind speed. It also has broad health and environmental implications. This study presents wet-bulb globe temperatures calculated from the fifth-generation European Centre for Medium-Range Weather Forecasts atmospheric reanalysis and combines it with health guidelines to assess heat stress variability and the potential for reduction in labor hours over the past decade on both the continental and urban scale. Compared to 2010–2014, there was a general increase in heat stress during the period from 2015 to 2019 throughout the northern hemisphere, with the largest warming found in tropical regions, especially in the northern part of the Indian Peninsula. On the urban scale, our results suggest that heat stress might have led to a reduction in labor hours by up to ~20% in some Asian cities subject to work–rest regulations. Extremes in heat stress can be explained by changes in radiation and circulation. The resultant threat is highest in developing countries in tropical areas where workers often have limited legal protection and healthcare. The effect of heat stress exposure is therefore a collective challenge with environmental, economic, and social implications.

The Norwegian Climate Prediction Model version 1 (NorCPM1) is a new research tool for performing climate reanalyses and seasonal-to-decadal climate predictions. It combines the Norwegian Earth System Model version 1 (NorESM1) – which features interactive aerosol-cloud schemes and an isopycnic-coordinate ocean component with biogeochemistry – with anomaly assimilation of SST and T/S-profile observations using the Ensemble Kalman Filter (EnKF). We first describe the Earth system component and the data assimilation (DA) scheme, highlighting implementation of new forcings, bug-fixes, re-tuning and DA innovations. Notably, NorCPM1 uses two anomaly assimilation variants to assess the impact of sea ice initialisation and climatological reference period: The first (i1) uses a 1980–2010 reference climatology for computing anomalies and the DA only updates the physical ocean state; the second (i2) uses a 1950–2010 reference climatology and additionally updates the sea ice state via strongly coupled DA of ocean observations. We then assess the baseline, reanalysis and prediction performance with output contributed to the Decadal Climate Prediction Project (DCPP) as part of the sixth Coupled Model Intercomparison Project (CMIP6). The non-assimilation experiments exhibit a moderate historical global surface temperature evolution and tropical climate variability characteristics that compare favourably with observations. The climate biases of NorCPM1 using CMIP6 external forcings, are comparable to, or slightly larger than those of the original NorESM1 CMIP5 model, with positive biases in Atlantic meridional overturning circulation (AMOC) strength and Arctic sea ice thickness, too cold subtropical oceans and northern continents, and a too warm North Atlantic and Southern Ocean. The biases in the assimilation experiments are mostly unchanged except for a reduced sea ice thickness bias in i2 caused by the assimilation update of sea ice, generally confirming that the anomaly assimilation synchronises variability without changing the climatology. The i1 and i2 reanalysis/hindcast products overall show comparable performance. The benefits of initialisation are seen globally in the first year of the prediction over a range of variables, also in the atmosphere and over land. External forcings are the primary source of multi-year skills, while added benefit from initialisation is demonstrated for the subpolar North Atlantic (SPNA) and its extension to the Arctic. Both products show limited success in constraining and predicting surface ocean biogeochemistry variability. However, observational uncertainties and short temporal coverage make biogeochemistry evaluation uncertain while potential predictability is found to be high. For physical climate prediction, i2 performs marginally better than i1 for a range of variables, especially in the SPNA and in the vicinity of sea ice, with notably improved sea level variability of the Southern Ocean. Despite similar skills, i1 and i2 feature very different drift behaviours, mainly due to their use of different climatologies in DA; i2 exhibits an anomalously strong AMOC that leads to forecast drift with unrealistic warming in the SPNA, whereas i1 exhibits a weaker AMOC that leads to unrealistic cooling. In polar regions, the reduction in climatological ice thickness in i2 causes additional forecast drift as the ice grows back. Posteriori lead dependent drift correction removes most hindcast differences; applications should therefore benefit from combining the two products. The results confirm that the large-scale ocean circulation exerts strong control on North Atlantic temperature variability, implying predictive potential from better synchronisation of circulation variability. Future development will therefore focus on improving the representation of mean state and variability of AMOC and its initialisation. Other efforts will be directed to refining the anomaly assimilation scheme – to better separate between internal versus forced signals, to include land and atmosphere initialisation and new observational types – and improving biogeochemistry prediction capability. Combined with other models, NorCPM1 may already contribute to skilful multi-year climate prediction that benefits society.

Large-scale interaction between the three tropical ocean basins is an area of intense research that is often conducted through experimentation with numerical models. A common problem is that modeling groups use different experimental setups, which makes it difficult to compare results and delineate the role of model biases from differences in experimental setups. To address this issue, an experimental protocol for examining interaction between the tropical basins is introduced. The Tropical Basin Interaction Model Intercomparison Project (TBIMIP) consists of experiments in which sea surface temperatures (SSTs) are prescribed to follow observed values in selected basins. There are two types of experiments. One type, called standard pacemaker, consists of simulations in which SSTs are restored to observations in selected basins during a historical simulation. The other type, called pacemaker hindcast, consists of seasonal hindcast simulations in which SSTs are restored to observations during 12-month forecast periods. TBIMIP is coordinated by the Climate and Ocean – Variability, Predictability, and Change (CLIVAR) Research Focus on Tropical Basin Interaction. The datasets from the model simulations will be made available to the community to facilitate and stimulate research on tropical basin interaction and its role in seasonal-to-decadal variability and climate change.

Climate reanalyses are essential for studying climate variability, understanding climate processes, and initializing climate predictions. We present CoRea1860+ (Wang and Counillon, 2025, https://doi.org/10.11582/2025.00009), a 30-member coupled reanalysis spanning from 1860 to the present, produced using the Norwegian Climate Prediction Model (NorCPM) and assimilating sea surface temperature (SST) observations. NorCPM combines the Norwegian Earth System Model with the ensemble Kalman filter data assimilation method. SST, available throughout the entire period, serves as the primary source of instrumental oceanic measurements prior to the 1950s. CoRea1860+ belongs to the category of sparse-input reanalyses, designed to minimize artefacts arising from changes in the observation network over time. By exclusively assimilating oceanic data, this reanalysis offers valuable insights into the ocean's role in driving climate system variability, including its influence on the atmosphere and sea ice. This study first describes the numerical model, the SST dataset, and the assimilation implementation used to produce CoRea1860+. It then provides a comprehensive evaluation of the reanalysis across four key aspects, namely reliability, ocean variability, sea ice variability, and atmospheric variability, benchmarked against more than 10 independent reanalyses and observational datasets. Overall, CoRea1860+ demonstrates strong reliability, particularly in observation-rich periods, and provides a reasonable representation of climate variability. It successfully captures key features such as multi-decadal variability and long-term trends in ocean heat content, the Atlantic meridional overturning circulation, and sea ice variability in both hemispheres. Furthermore, to some extent, CoRea1860+ agrees with the reference atmospheric datasets for surface air temperature, precipitation, sea level pressure, and 500 hPa geopotential height, especially in the tropics where air–sea interactions are most pronounced.