On dynamical downscaling of ENSO-induced oceanic anomalies off Baja California Peninsula, Mexico: role of the air-sea heat flux

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.publishedVersio

Forecasting harmful algae blooms: Application to Dinophysis acuminata in northern Norway

Dinophysis acuminata produces Diarrhetic Shellfish Toxins (DST) that contaminate natural and farmed shellfish, leading to public health risks and economically impacting mussel farms. For this reason, there is a high interest in understanding and predicting D. acuminata blooms. This study assesses the environmental conditions and develops a sub-seasonal (7 - 28 days) forecast model to predict D. acuminata cells abundance in the Lyngen fjord located in northern Norway. A Support Vector Machine (SVM) model is trained to predict future D. acuminata cells abundance by using the past cell concentration, sea surface temperature (SST), Photosynthetic Active Radiation (PAR), and wind speed. Cells concentration of Dinophysis spp. are measured in-situ from 2006 to 2019, and SST, PAR, and surface wind speed are obtained by satellite remote sensing. D. acuminata only explains 40% of DST variability from 2006 to 2011, but it changes to 65% after 2011 when D. acuta prevalence reduced. The D. acuminata blooms can reach concentration up to 3954 cells l−1 and are restricted to the summer during warmer waters, varying from 7.8 to 12.7 °C. The forecast model predicts with fair accuracy the seasonal development of the blooms and the blooms amplitude, showing a coefficient of determination varying from 0.46 to 0.55. SST has been found to be a useful predictor for the seasonal development of the blooms, while the past cells abundance is needed for updating the current status and adjusting the blooms timing and amplitude. The calibrated model should be tested operationally in the future to provide an early warning of D. acuminata blooms in the Lyngen fjord. The approach can be generalized to other regions by recalibrating the model with local observations of D. acuminata blooms and remote sensing data.publishedVersio

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

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.publishedVersio

Framework for an Ocean-Connected Supermodel of the Earth System

A supermodel connects different models interactively so that their systematic errors compensate and achieve a model with superior performance. It differs from the standard non-interactive multi-model ensembles (NI), which combines model outputs a-posteriori. Supermodels with Earth system models (ESMs) has not been developed because it is technically challenging to combine models with different state space. Here, we formulate the first supermodel framework for ESMs and use data assimilation to synchronise models. The ocean of three ESMs is synchronised every month by assimilating pseudo sea surface temperature (SST) observations generated by them on a common grid to handle discrepancies in grid and resolution. We compare the performance of two supermodel approaches to that of the NI. In the first (EW), the models are connected to the equal-weight multi-model mean, while in the second (SINGLE), they are connected to a single model. Both versions achieve synchronisation in the ocean and in the atmosphere, where the ocean drives the variability. The time variability of the supermodel multi-model mean SST is reduced compared to observations, most where synchronisation is not achieved and is lower-bounded by NI. The damping is larger in EW, for which variability in the individual models is also damped. Hence, under partial synchronisation, the unsynchronized variability gets damped in the multi-model average pseudo-observations, causing a deflation during the assimilation. The SST bias in individual models of EW is reduced compared to that of NI, and so is its multi-model mean in the synchronised regions. A trained supermodel remains to be tested.publishedVersio

Propagation of Thermohaline Anomalies and their predictive potential along the Atlantic water pathway

We assess to what extent seven state-of-the-art dynamical prediction systems can retrospectively predict winter sea surface temperature (SST) in the subpolar North Atlantic and the Nordic seas in the period 1970–2005. We focus on the region where warm water flows poleward (i.e., the Atlantic water pathway to the Arctic) and on interannual-to-decadal time scales. Observational studies demonstrate predictability several years in advance in this region, but we find that SST skill is low with significant skill only at a lead time of 1–2 years. To better understand why the prediction systems have predictive skill or lack thereof, we assess the skill of the systems to reproduce a spatiotemporal SST pattern based on observations. The physical mechanism underlying this pattern is a propagation of oceanic anomalies from low to high latitudes along the major currents, the North Atlantic Current and the Norwegian Atlantic Current. We find that the prediction systems have difficulties in reproducing this pattern. To identify whether the misrepresentation is due to incorrect model physics, we assess the respective uninitialized historical simulations. These simulations also tend to misrepresent the spatiotemporal SST pattern, indicating that the physical mechanism is not properly simulated. However, the representation of the pattern is slightly degraded in the predictions compared to historical runs, which could be a result of initialization shocks and forecast drift effects. Ways to enhance predictions could include improved initialization and better simulation of poleward circulation of anomalies. This might require model resolutions in which flow over complex bathymetry and the physics of mesoscale ocean eddies and their interactions with the atmosphere are resolved.publishedVersio

NorCPM1 and its contribution to CMIP6 DCPP

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).publishedVersio

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

Under embargo until: 2022-10-01As 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.publishedVersio

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

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.publishedVersio

Super-resolution data assimilation

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% and the errors reduce by 40%, making the performance very close to that of the high-resolution system (52% of error reduction) that increases the cost by 800%. The reliability of the ensemble system is not degraded by SRDA.publishedVersio

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

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