Tag: counillon

Phytoplankton abundance in the Barents Sea is predictable up to five years in advance

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.

Framework for an Ocean-Connected Supermodel of the Earth System

Counillon, F., Keenlyside, N., Wang, S., Devilliers, M., Gupta, A., Koseki, S., Shen, M.-L. 2023: Framework for an Ocean-Connected Supermodel of the Earth System. JAMES. https://doi.org/10.1029/2022MS003310

Summary: Observed and future winter Arctic sea ice loss is strongest in the Barents Sea. However, the anthropogenic signal of the sea ice decline is superimposed by pronounced internal variability that represents a large source of uncertainty in future climate projections. A notable manifestation of internal variability is rapid ice change events (RICEs) that greatly exceed the anthropogenic trend. These RICEs are associated with large displacements of the sea ice edge which could potentially have both local and remote impacts on the climate system. In this study we present the first investigation of the frequency and drivers of RICEs in the future Barents Sea, using multi-member ensemble simulations from CMIP5 and CMIP6. A majority of RICEs are triggered by trends in ocean heat transport or surface heat fluxes. Ice loss events are associated with increasing trends in ocean heat transport and decreasing trends in surface heat loss. RICEs are a common feature of the future Barents Sea until the region becomes close to ice-free. As their evolution over time is closely tied to the average sea ice conditions, rapid ice changes in the Barents Sea may serve as a precursor for future changes in adjacent seas.

Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Benefit of vertical localization for sea surface temperature assimilation in isopycnal coordinate model

Wang, Y., Counillon, F., Barthélémy, S., Barth, A. 2022: Benefit of vertical localization for sea surface temperature assimilation in isopycnal coordinate model. Front Clim. https://doi.org/10.3389/fclim.2022.918572

Summary: Sea surface temperature (SST) observations are a critical data set for long-term climate reconstruction. However, their assimilation with an ensemble-based data assimilation method can degrade performance in the ocean interior due to spurious covariances. Assimilation in isopycnal coordinates can delay the degradation, but it remains problematic for long reanalysis. We introduce vertical localization for SST assimilation in the isopycnal coordinate. The tapering functions are formulated empirically from a large pre-industrial ensemble. We propose three schemes: 1) a step function with a small localization radius that updates layers from the surface down to the first layer for which insignificant correlation with SST is found, 2) a step function with a large localization radius that updates layers down to the last layer for which significant correlation with SST is found, and 3) a flattop smooth tapering function. These tapering functions vary spatially and with the calendar month and are applied to isopycnal temperature and salinity. The impact of vertical localization on reanalysis performance is tested in identical twin experiments within the Norwegian Climate Prediction Model (NorCPM) with SST assimilation over the period 1980–2010. The SST assimilation without vertical localization greatly enhances performance in the whole water column but introduces a weak degradation at intermediate depths (e.g., 2,000–4,000 m). Vertical localization greatly reduces the degradation and improves the overall accuracy of the reanalysis, in particular in the North Pacific and the North Atlantic. A weak degradation remains in some regions below 2,000 m in the Southern Ocean. Among the three schemes, scheme 2) outperforms schemes 1) and 3) for temperature and salinity.

Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Super-resolution data assimilation

Barthélémy, S., Brajard, J., Bertino, L., Counillon, F. 2022: Super-resolution data assimilation. Ocean Dyn. https://doi.org/10.1007/s10236-022-01523-x

Summary: Increasing model resolution can improve the performance of a data assimilation system because it reduces model error, the system can more optimally use high-resolution observations, and with an ensemble data assimilation method the forecast error covariances are improved. However, increasing the resolution scales with a cubical increase of the computational costs. A method that can more effectively improve performance is introduced here. The novel approach called “Super-resolution data assimilation” (SRDA) is inspired from super-resolution image processing techniques and brought to the data assimilation context. Starting from a low-resolution forecast, a neural network (NN) emulates the fields to high-resolution, assimilates high-resolution observations, and scales it back up to the original resolution for running the next model step. The SRDA is tested with a quasi-geostrophic model in an idealized twin experiment for configurations where the model resolution is twice and four times lower than the reference solution from which pseudo-observations are extracted. The assimilation is performed with an Ensemble Kalman Filter. We show that SRDA outperforms both the low-resolution data assimilation approach and a version of SRDA with cubic spline interpolation instead of NN. The NN’s ability to anticipate the systematic differences between low- and high-resolution model dynamics explains the enhanced performance, in particular by correcting the difference of propagation speed of eddies. With a 25-member ensemble at low resolution, the SRDA computational overhead is 55 percent and the errors reduce by 40 percent, making the performance very close to that of the high-resolution system (52 percent of error reduction) that increases the cost by 800 percent. The reliability of the ensemble system is not degraded by SRDA.

Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Editorial: Past Reconstruction of the Physical and Biogeochemical Ocean State

Masina, S., Counillon, F., Grégoire, M., Storto, A., Tsujino, H. 2022: Editorial: Past Reconstruction of the Physical and Biogeochemical Ocean State. Front Earth Sci. https://doi.org/10.3389/feart.2022.890370

Summary: Knowledge of the ocean’s physical, biogeochemical and ecosystem state and variability is crucial for understanding the evolution of our climate system and better predicting its future. However, the sparseness and inhomogeneous distribution of observations hinder the creation of sound 4-dimensional reconstructions of the past (for an overview of ocean observing systems see the Research Topic Oceanobs’19: An Ocean of Opportunity). Instead, we must rely on a combination of ocean modeling and data analysis to infer past changes. Over the last decade the quality of ocean reanalyses has improved mainly thanks to advances in data assimilation methods and more quality-controlled observation data sets. Reanalyses provide the best-possible state estimate by assimilating observations into a dynamical model (Balmaseda et al., 2015; Masina and Storto, 2017; Storto et al., 2019). In addition, advanced statistical mapping methods (e.g., objective or variational analysis) provide observation-based gridded fields whose resolution depends on the amount of available data (among many Cheng et al., 2017, Ishii et al., 2017; Boyer et al., 2018). For many variables, particularly biogeochemical, the lack of observations more strongly limits the spatial and temporal resolution of these gridded products (Fennel et al., 2019).

The Research Topic gathers contributions aiming at reconstructing the past physical, sea ice and biogeochemical state of the ocean using models in combination with data. Ocean reanalyses and observation-mapping are proposed to further our knowledge, to demonstrate their use in supporting various applications, and to increase confidence in these reconstructions within the scientific community. The products and applications described in this topic provide a foundation for their use in ecosystem-based management, policy advice to support mitigation and adaptation strategies, and in the identification of pathways towards a sustainable ocean.

Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

WMO Global Annual to Decadal Climate Update: A Prediction for 2021–25

Hermanson, L., Smith, D., Seabrook, M., Bilbao, R., Doblas-Reyes, F., Tourigny, E., Lapin, V., Kharin, V.V., Merryfield, W.J., Sospedra-Alfonso, R., Athanasiadis, P., Nicoli, D., Gualdi, S., Dunstone, N., Eade, R., Scaife, A., Collier, M., O’Kane, T., Kitsios, V., Sandery, P., Pankatz, K., Früh, B., Pohlmann, H., Müller, W., Kataoka, T., Tatebe, H., Ishii M., Imada, Y., Kruschke, T., Koenigk, T., Pasha Karami, M., Yang, S., Tian, T., Zhang, L., Delworth, T., Yang, X., Zeng, F., Wang, Y., Counillon, F., Keenlyside, N.S., Bethke, I., Lean, J., Luterbacher, J., Kumar Kolli, R., Kumar, A. 2022: WMO Global Annual to Decadal Climate Update: A Prediction for 2021–25. BAMS https://doi.org/10.1175/BAMS-D-20-0311.1 .

Summary: As climate change accelerates, societies and climate-sensitive socioeconomic sectors cannot continue to rely on the past as a guide to possible future climate hazards. Operational decadal predictions offer the potential to inform current adaptation and increase resilience by filling the important gap between seasonal forecasts and climate projections. The World Meteorological Organization (WMO) has recognized this and in 2017 established the WMO Lead Centre for Annual to Decadal Climate Predictions (shortened to “Lead Centre” below), which annually provides a large multimodel ensemble of predictions covering the next 5 years. This international collaboration produces a prediction that is more skillful and useful than any single center can achieve. One of the main outputs of the Lead Centre is the Global Annual to Decadal Climate Update (GADCU), a consensus forecast based on these predictions. This update includes maps showing key variables, discussion on forecast skill, and predictions of climate indices such as the global mean near-surface temperature and Atlantic multidecadal variability. it also estimates the probability of the global mean temperature exceeding 1.5°C above preindustrial levels for at least 1 year in the next 5 years, which helps policy-makers understand how closely the world is approaching this goal of the Paris Agreement. This paper, written by the authors of the GADCU, introduces the GADCU, presents its key outputs, and briefly discusses its role in providing vital climate information for society now and in the future..

Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Propagation of Thermohaline Anomalies and Their Predictive Potential along the Atlantic Water Pathway

Langehaug, H. R., Ortega, P., Counillon, F., Matei, D., Maroon, E., Keenlyside, N., Mignot, J., Wang, Y., Swingedouw, D., Bethke, I., Yang, S., Danabasoglu, G., Bellucci, A., Ruggieri, P., Nicolì, D., Årthun, M. 2022: Propagation of Thermohaline Anomalies and Their Predictive Potential along the Atlantic Water Pathway. J Clim. https://doi.org/10.1007/s10236-022-01523-x

Summary: In this study, we find that dynamical prediction systems and their respective climate models struggle to realistically represent ocean surface temperature variability in the eastern subpolar North Atlantic and Nordic seas on interannual-to-decadal time scales. In previous studies, ocean advection is proposed as a key mechanism in propagating temperature anomalies along the Atlantic water pathway toward the Arctic Ocean. Our analysis suggests that the predicted temperature anomalies are not properly circulated to the north; this is a result of model errors that seems to be exacerbated by the effect of initialization shocks and forecast drift. Better climate predictions in the study region will thus require improving the initialization step, as well as enhancing process representation in the climate models.

Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Estimation of Ocean Biogeochemical Parameters in an Earth System Model Using the Dual One Step Ahead Smoother: A Twin Experiment

Singh, T., Counillon, F., Tjiputra, J., Wang Y., El Gharamti, M. 2022: Estimation of Ocean Biogeochemical Parameters in an Earth System Model Using the Dual One Step Ahead Smoother: A Twin Experiment. Front. Mar. Sci. https://doi.org/10.3389/fmars.2022.775394 .

For an easy-to-understand overview of this publication, produced in collaboration with the TRIATLAS project, we recommend starting with this neat article written by Henrike Wilborn, at NERSC: “Making climate models more accurate by improving their tuning“.

Summary: Ocean biogeochemical (BGC) models utilise a large number of poorly-constrained global parameters to mimic unresolved processes and reproduce the observed complex spatio-temporal patterns. Large model errors stem primarily from inaccuracies in these parameters whose optimal values can vary both in space and time. This study aims to demonstrate the ability of ensemble data assimilation (DA) methods to provide high-quality and improved BGC parameters within an Earth system model in an idealized perfect twin experiment framework. We use the Norwegian Climate Prediction Model (NorCPM), which combines the Norwegian Earth System Model with the Dual-One-Step ahead smoothing-based Ensemble Kalman Filter (DOSA-EnKF). We aim to estimate five spatially varying BGC parameters by assimilating salinity and temperature profiles and surface BGC (Phytoplankton, Nitrate, Phosphate, Silicate, and Oxygen) observations in a strongly coupled DA framework—i.e., jointly updating ocean and BGC state-parameters during the assimilation. We show how BGC observations can effectively constrain error in the ocean physics and vice versa. The method converges quickly (less than a year) and largely reduces the errors in the BGC parameters. Some parameter error remains, but the resulting state variable error using the estimated parameters for a free ensemble run and for a reanalysis performs nearly as well as with true parameter values. Optimal parameter values can also be recovered by assimilating climatological BGC observations or sparse observational networks. The findings of this study demonstrate the applicability of the DA approach for tuning the system in a real framework.

Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

NorCPM1 and its contribution to CMIP6 DCPP

Bethke, I., Wang, Y., Counillon, F., Keenlyside, N., Kimmritz, M., Fransner, F., Samuelsen, A., Langehaug, H., Svendsen, L., Chiu, P.-G., Passos, L., Bentsen, M., Guo, C., Gupta, A., Tjiputra, J., Kirkevåg, A., Olivié, D., Seland, Ø., Solsvik Vågane, J., Fan, Y., Eldevik, T. 2021: NorCPM1 and its contribution to CMIP6 DCPP. Geosci Model Dev. https://doi.org/10.5194/gmd-14-7073-2021 .

For an easy-to-understand overview, we recommend starting with this neat article written by the Climate Futures team, a project connected to BCPU: “New Study: Decadal Climate Forecasts From The Norwegian Climate Prediction Model” (les heller på norsk).

Summary: 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 -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.

Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.

Twenty-one years of phytoplankton bloom phenology in the Barents, Norwegian and North seas

Silva, E.F.F., Counillon, F., Brajard, J., Korosov, A., Pettersson, L., Samuelsen, A., Keenlyside, N. 2021: Twenty-one years of phytoplankton bloom phenology in the Barents, Norwegian and North seas. Front Mar Sci.  https://doi.org/10.3389/fmars.2021.746327 .

For en flott oppsummering på norsk, les denne artikkelen av vår samarbeidspartner, Climate Futures.

Summary: Phytoplankton blooms provide biomass to the marine trophic web, contribute to the carbon removal from the atmosphere and can be deadly when associated with harmful species. This points to the need to understand the phenology of the blooms in the Barents, Norwegian, and North seas. We use satellite chlorophyll-a from 2000 to 2020 to assess robust climatological and the interannual trends of spring and summer blooms onset, peak day, duration and intensity. Further, we also correlate the interannual variability of the blooms with mixed layer depth (MLD), sea surface temperature (SST), wind speed and suspended particulate matter (SPM) retrieved from models and remote sensing. The climatological spring blooms start on March 10th and end on June 19th. The climatological summer blooms begin on July 13th and end on September 17th. In the Barents Sea, years of shallower mixed layer (ML) driven by both calm waters and higher freshwaters input keeps the phytoplankton in the euphotic zone, causing the spring bloom to start earlier and reach higher biomass but end sooner due to the lack of nutrients upwelling from the deep. In the Norwegian Sea, a correlation between SST and the spring blooms is found. Here, warmer waters are correlated to earlier and stronger blooms in most regions but with later and weaker blooms in the eastern Norwegian Sea. In the North Sea, years of shallower ML reduces the phytoplankton sinking below the euphotic zone and limits the SPM increase from the bed shear stress, creating an ideal environment of stratified and clear waters to develop stronger spring blooms. Last, the summer blooms onset, peak day and duration have been rapidly delaying at a rate of 1.25-day year–1, but with inconclusive causes based on the parameters assessed in this study.

Link to publication. You are most welcome to contact us or the corresponding author(s) directly, if you have questions.