Bjerknes Climate Prediction Unit

Sospedra-Alfonso, R., Merryfield, W.J., Toohey, M., Timmreck, C., Vernier, J-P., Bethke, I., Wang, Y., Bilbao, R., Donat, M.G., Ortega, P., Cole, J., Lee, W.-S., Delworth, T.L., Paynter, D., Zeng, F., Zhang, L., Khodri, M., Mignot, J., Swingedouw, D., Torres, O., Hu, S., Man, W., Zuo, M., Hermanson, L., Smith, D., Kataoka, T., Tatebe, H. 2024: Decadal prediction centers prepare for a major volcanic eruption. Bulletin of the American Meteorological Society. https://doi.org/10.1175/BAMS-D-23-0111.1 Summary: The World Meteorological Organization’s Lead Centre for Annual-to-Decadal Climate prediction issues operational forecasts annually as guidance for regional climate centers, climate outlook forums and national meteorological and hydrological services. The occurrence of a large volcanic eruption such as that of Mount Pinatubo in 1991, however, would invalidate these forecasts and prompt producers to modify their predictions. To assist and prepare decadal prediction centers for this eventuality, the Volcanic Response activities under the World Climate Research Programme’s Stratosphere-troposphere Processes And their Role in Climate (SPARC) and the Decadal Climate Prediction Project (DCPP) organized a community exercise to respond to a hypothetical large eruption occurring in April 2022. As part of this exercise, the Easy Volcanic Aerosol forcing generator was used to provide stratospheric sulfate aerosol optical properties customized to the configurations of individual decadal prediction models. Participating centers then reran forecasts for 2022-2026 from their original initialization dates and in most cases also from just before the eruption at the beginning of April 2022, according to two candidate response protocols. This article describes various aspects of this SPARC/DCPP Volcanic Response Readiness Exercise (VolRes-RE), including the hypothesized volcanic event, the modified forecasts under the two protocols from the eight contributing centers, the lessons learned during the coordination and execution of this exercise, and the recommendations to the decadal prediction community for the response to an actual eruption. Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Bushuk, M., Ali, S., Bailey, D.A., Bao, Q., Batté, L., Bhatt, U.S., Blanchard-Wrigglesworth, E., Blockley, E., Cawley, G., Chi, J., Counillon, F., et al. 2024: Predicting September Arctic Sea Ice: A Multimodel Seasonal Skill Comparison. Bull. Amer. Meteor. Soc.. https://doi.org/10.1175/BAMS-D-23-0163.1 Summary: This study quantifies the state of the art in the rapidly growing field of seasonal Arctic sea ice prediction. A novel multimodel dataset of retrospective seasonal predictions of September Arctic sea ice is created and analyzed, consisting of community contributions from 17 statistical models and 17 dynamical models. Prediction skill is compared over the period 2001–20 for predictions of pan-Arctic sea ice extent (SIE), regional SIE, and local sea ice concentration (SIC) initialized on 1 June, 1 July, 1 August, and 1 September. This diverse set of statistical and dynamical models can individually predict linearly detrended pan-Arctic SIE anomalies with skill, and a multimodel median prediction has correlation coefficients of 0.79, 0.86, 0.92, and 0.99 at these respective initialization times. Regional SIE predictions have similar skill to pan-Arctic predictions in the Alaskan and Siberian regions, whereas regional skill is lower in the Canadian, Atlantic, and central Arctic sectors. The skill of dynamical and statistical models is generally comparable for pan-Arctic SIE, whereas dynamical models outperform their statistical counterparts for regional and local predictions. The prediction systems are found to provide the most value added relative to basic reference forecasts in the extreme SIE years of 1996, 2007, and 2012. SIE prediction errors do not show clear trends over time, suggesting that there has been minimal change in inherent sea ice predictability over the satellite era. Overall, this study demonstrates that there are bright prospects for skillful operational predictions of September sea ice at least 3 months in advance. Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Richter, I., Ratnam, J.V., Martineau, P., Oettli, P., Doi, T., Ogata, T., Kataoka, T., Counillon, F. 2024: A simple statistical post-processing scheme for enhancing the skill of seasonal SST predictions in the tropics. Monthly Weather Review. https://doi.org/10.1175/MWR-D-23-0266.1 Summary: The prediction of year-to-year climate variability patterns, such as El Niño, offers potential benefits to society by aiding mitigation and adaptation efforts. Current prediction systems, however, may still have substantial room for improvement due to systematic model errors and due to imperfect initialization of the oceanic state at the start of predictions. Here we develop a statistical correction scheme to improve prediction skill after forecasts have been completed. The scheme shows some moderate success in improving the skill for predicting El Niño and similar climate patterns in seven prediction systems. Our results not only indicate a potential for improving prediction skill after the fact but also point to the importance of improving the way prediction systems are initialized. Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Garcia-Oliva, L., Counillon, F., Bethke, I., Keenlyside, N. 2024: Intercomparison of initialization methods for Seasonal-to-Decadal Climate Predictions with the NorCPM. Clim Dyn. https://doi.org/10.1007/s00382-024-07170-w Summary: Initialization is essential for accurate seasonal-to-decadal (S2D) climate predictions. The initialization schemes used differ on the component initialized, the Data Assimilation method, or the technique. We compare five popular schemes within NorCPM following the same experimental protocol: reanalysis from 1980 to 2010 and seasonal and decadal predictions initialized from the reanalysis. We compare atmospheric initialization—Newtonian relaxation (nudging)—against ocean initialization—Ensemble Kalman Filter—(ODA). On the atmosphere, we explore the benefit of full-field (NudF-UVT) or anomaly (NudA-UVT) nudging of horizontal winds and temperature (U, V, and T) observations. The scheme NudA-UV nudges horizontal winds to disentangle the role of wind-driven variability. The ODA+NudA-UV scheme is a first attempt at joint initialization of ocean and atmospheric components in NorCPM. During the reanalysis, atmospheric nudging improves the synchronization of the atmosphere and land components with the observed data. Conversely, ODA is more effective at synchronizing the ocean component with observations. The atmospheric nudging schemes are better at reproducing specific events, such as the rapid North Atlantic subpolar gyre shift. An abrupt climatological change using the NudA-UV scheme demonstrates that energy conservation is crucial when only assimilating winds. ODA outperforms atmospheric-initialized versions for S2D global predictions, while atmospheric nudging is preferable for accurately initializing phenomena in specific regions, with the technique’s benefit depending on the prediction’s temporal scale. For instance, atmospheric full-field initialization benefits the tropical Atlantic Niño at 1-month lead time, and atmospheric anomaly initialization benefits longer lead times, reducing hindcast drift. Combining atmosphere and ocean initialization yields sub-optimal results, as sustaining the ensemble’s reliability—required for ODA’s performance—is challenging with atmospheric nudging. Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Polkova, I, Swingedouw, D., Hermanson, L., Köhl, A., Stammer, D., Smith, D., Kröger, J., Bethke, I., Yang, X., Zhang, L., Nicolì, D., Athanasiadis, P., Karami, P., Pankatz, K., Pohlmann, H., Wu, B., Bilbao, R., Ortega, P., Yang, S., Sospedra-Alfonso, R., Merryfield, W., Kataoka, T., Tatebe, H., Imada, Y., Ishii, M., Matear, R. 2023: Initialization shock in the ocean circulation reduces skill in decadal predictions of the North Atlantic subpolar gyre. Front Clim. doi: https://doi.org/10.3389/fclim.2023.1273770 Summary: Due to large northward heat transport, the Atlantic meridional overturning circulation (AMOC) strongly affects the climate of various regions. Its internal variability has been shown to be predictable decades ahead within climate models, providing the hope that synchronizing ocean circulation with observations can improve decadal predictions, notably of the North Atlantic subpolar gyre (SPG). Climate predictions require a starting point which is a reconstruction of the past climate. This is usually performed with data assimilation methods that blend available observations and climate model states together. There is no unique method to derive the initial conditions. Moreover, this can be performed using full-field observations or their anomalies superimposed on the model’s climatology to avoid strong drifts in predictions. How critical ocean circulation drifts are for prediction skill has not been assessed yet. We analyze this possible connection using the dataset of 12 decadal prediction systems from the World Meteorological Organization Lead Centre for Annual-to-Decadal Climate Prediction. We find a variety of initial AMOC errors within the predictions related to a dynamically imbalanced ocean states leading to strongly displaced or multiple maxima in the overturning structures. This likely results in a blend of what is known as model drift and initial shock. We identify that the AMOC initialization influences the quality of the SPG predictions. When predictions show a large initial error in their AMOC, they usually have low skill for predicting internal variability of the SPG for a time horizon of 6-10 years. Full-field initialized predictions with low AMOC drift show better SPG skill than those with a large AMOC drift. Nevertheless, while the anomaly-initialized predictions do not experience large drifts, they show low SPG skill when skill also present in historical runs is removed using a residual correlation metric. Thus, reducing initial shock and model biases for the ocean circulation in prediction systems might help to improve their prediction for the SPG beyond 5 years. Climate predictions could also benefit from quality-check procedure for assimilation/initialization because currently the research groups only reveal the problems in initialization once the set of predictions has been completed, which is an expensive effort. Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Goncalves Dos Passos, Leilane (2023-11-03): Arctic-Atlantic Climate Variability and Predictability in Observations and in a Dynamical Prediction System. PhD thesis, University of Bergen, Bergen, Norway. https://bora.uib.no/bora-xmlui/handle/11250/3099594 Summary: The major focus of this thesis is on understanding decadal climate predictability to improve climate models and their predictions. Climate predictions show promising results but are still facing challenges, especially in connecting the ocean and atmosphere. The ocean is the main source of predictability. The ocean’s capacity to store and release heat over long periods of time makes it a thermal memory of the climate system. In the Arctic-Atlantic region, ocean currents transport heat to polar areas, and along this path, the ocean releases the heat to the atmosphere through surface fluxes. From this interaction, both the ocean and the atmosphere change. On the one hand, as the ocean releases heat into the atmosphere, it cools down, increasing its density. The denser water eventually flows southward as part of the Atlantic Meridional Overturning Circulation (AMOC). On the other hand, the atmosphere being warmed by the ocean affects nearby land areas through the winds, influencing the climate variability of Western Europe. This dynamic ocean-atmosphere interaction is a source of predictability in the Arctic-Atlantic region and is investigated here using observations and a dynamical prediction system, the Norwegian Climate Prediction Model (NorCPM). Dynamical prediction systems are useful tools for investigating and predicting climate variability on decadal timescales. Beginning their development in the early 2000s, these systems are currently the focus of significant efforts by the scientific community to provide operational decadal forecasts with reliable and accurate information. The research of this thesis is aligned with the development of NorCPM while also focusing on investigating key mechanisms that give rise to predictability in the Arctic-Atlantic region. Climate predictions are initialized in different ways, which affects their performance. The first study of the thesis investigates the best initialization method for the Arctic-Atlantic region using NorCPM. Paper I finds that employing a more complex data assimilation method leads to the improved predictive skill of temperature and salinity in the Subpolar North Atlantic (SPNA) but not in the Norwegian Sea. The loss of skill in the Norwegian Sea is found in regions characterized by intense surface heat fluxes and eddy activity, such as the Norwegian and Lofoten Basins. The warm Atlantic water moving northwards from the SPNA to the Norwegian Sea carries thermohaline anomalies, and it is transformed from light-to-dense waters by surface forcing along the path. These two mechanisms are investigated in observation-based data in Paper II. Their relationship is analyzed, focusing on the decadal timescale in the eastern SPNA. Paper II finds that warm anomalies are associated with surface-forced water mass transformation in the light-density classes, while during cold anomalies, more transformation happens in denser classes. This relationship was disrupted during the Great Salinity Anomaly events of the 70s and 90s. Furthermore, the study highlights a faster propagation of thermohaline anomalies in the SPNA compared to the Norwegian Sea, particularly regarding temperature. The influence of the ocean on the climate of Europe is investigated in Paper III. This study advances the understanding of how constrained ocean variability impacts the atmosphere of NorCPM. The results show a more realistic thermodynamic component of surface air temperature (SAT) over the ocean and some European regions. Paper III shows that there is potential to improve multi-annual to decadal predictions over Europe, which is currently challenging in prediction systems. The research presented in this Thesis enhances the understanding of climate predictability in the Arctic-Atlantic region. It provides insights into the interactions between the atmosphere and ocean and adds to the development of the Norwegian Climate Prediction Model, contributing to making this prediction system operational in the coming years. Following similar approaches as presented in this thesis for other dynamical prediction systems would be highly recommended. Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Rivas, D., Counillon, F., Keenlyside, N. 2023: On dynamical downscaling of ENSO-induced oceanic anomalies off Baja California Peninsula, Mexico: role of the air-sea heat flux. Front Mar Sci. https://doi.org/10.3389/fmars.2023.1179649 Summary: The El Niño Southern Oscillation (ENSO) phenomenon is responsible for important physical and biogeochemical anomalies in the Northeastern Pacific Ocean. The event of 1997-98 has been one of the most intense in the last decades and it had large implications for the waters off Baja California (BC) Peninsula with a pronounced warm sea surface temperature (SST) anomaly adjacent to the coast. Downscaling of reanalysis products was carried out using a mesoscale-resolving numerical ocean model to reproduce the regional SST anomalies. The nested model has a 9 km horizontal resolution that extend from Cabo Corrientes to Point Conception. A downscaling experiment that computes surface fluxes online with bulk formulae achieves a better representation of the event than a version with prescribed surface fluxes. The nested system improves the representation of the large scale warming and the localized SST anomaly adjacent to BC Peninsula compared to the reanalysis product. A sensitivity analysis shows that air temperature and to a lesser extent wind stress anomalies are the primary drivers of the formation of BC temperature anomaly. The warm air-temperature anomalies advect from the near-equatorial regions and the central north Pacific and is associated with sea-level pressure anomalies in the synoptic-scale atmospheric circulation. This regional warm pool has a pronounced signature on sea level anomaly in agreement with observations, which may have implications for biogeochemistry. Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Fransner, F., Olsen, A., Årthun, M., Counillon, F., Tjiputra, J., Samuelsen, A., Keenlyside, N. 2023: Phytoplankton abundance in the Barents Sea is predictable up to five years in advance. Commun Earth Environ. https://doi.org/10.1038/s43247-023-00791-9 Summary: The Barents Sea is a highly biologically productive Arctic shelf sea with several commercially important fish stocks. Interannual-to-decadal predictions of its ecosystem would therefore be valuable for marine resource management. Here, we demonstrate that the abundance of phytoplankton, the base of the marine food web, can be predicted up to five years in advance in the Barents Sea with the Norwegian Climate Prediction Model. We identify two different mechanisms giving rise to this predictability; 1) in the southern ice-free Atlantic Domain, skillful prediction is a result of the advection of waters with anomalous nitrate concentrations from the Subpolar North Atlantic; 2) in the northern Polar Domain, phytoplankton predictability is a result of the skillful prediction of the summer ice concentration, which influences the light availability. The skillful prediction of the phytoplankton abundance is an important step forward in the development of numerical ecosystem predictions of the Barents Sea. Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Pou, J.M.H., Langehaug, H.R. 2022: Predictive skill in the Nordic Seas. Nansen Environmental and Remote Sensing Center, NERSC Technical Report (413). You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Passos, L., Langehaug, HR., Årthun, M., Eldevik, T., Bethke, I., Kimmritz, M. 2022: Impact of initialization methods on the predictive skill in NorCPM: an Arctic–Atlantic case study. Clim Dyn. https://doi.org/10.1007/s00382-022-06437-4 Summary: The skilful prediction of climatic conditions on a forecast horizon of months to decades into the future remains a main scientific challenge of large societal benefit. Here we assess the hindcast skill of the Norwegian Climate Prediction Model (NorCPM) for sea surface temperature (SST) and sea surface salinity (SSS) in the Arctic–Atlantic region focusing on the impact of different initialization methods. We find the skill to be distinctly larger for the Subpolar North Atlantic than for the Norwegian Sea, and generally for all lead years analyzed. For the Subpolar North Atlantic, there is furthermore consistent benefit in increasing the amount of data assimilated, and also in updating the sea ice based on SST with strongly coupled data assimilation. The predictive skill is furthermore significant for at least two model versions up to 8–10 lead years with the exception for SSS at the longer lead years. For the Norwegian Sea, significant predictive skill is more rare; there is relatively higher skill with respect to SSS than for SST. A systematic benefit from more complex data assimilation approach can not be identified for this region. Somewhat surprisingly, skill deteriorates quite consistently for both the Subpolar North Atlantic and the Norwegian Sea when going from CMIP5 to corresponding CMIP6 versions. We find this to relate to change in the regional performance of the underlying physical model that dominates the benefit from initialization. Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.