TOPAZ4b¶
Copernicus Marine Arctic Ocean Physics Reanalysis — TOPAZ4b (ARCTIC_MULTIYEAR_PHY_002_003)
Available
Current status
Active, extended twice a year. In June 2025 an AI-corrected sea ice thickness variable (sithick_corr) was added to remove the thin bias before 2010. The Arctic analysis and forecast system from the same centre (ARCTIC_ANALYSISFORECAST_PHY_002_001) moved to TOPAZ5 in November 2023, but the reanalysis is still TOPAZ4b. As of October 2026 we found no official plan to move this reanalysis to TOPAZ5 or TOPAZ6.
At a glance¶
| Item | Details |
|---|---|
| Section | Ocean › Reanalyses (also listed in: Sea ice › Reanalyses) |
| Type | Reanalysis |
| Domains | Ocean, Sea ice |
| Provider | Nansen Environmental and Remote Sensing Center (NERSC), Arctic Monitoring and Forecasting Centre / distributed by the Copernicus Marine Service (CMEMS) |
| Period | 1991-01 ~ present |
| Coverage | Arctic & N. Atlantic — Arctic Ocean and North Atlantic; the default product covers 50°N–90°N |
| Horizontal resolution | Model grid: HYCOM 880 × 800 Arakawa C-grid, about 12–16 km spacing (Arctic Ocean and North Atlantic). Original output grid (originalGrid): 12.5 km polar stereographic, 609 × 881 (x −38 to 38, y −55 to 55 in units of 100 km; +proj=stere +lon_0=-45 +lat_0=90 +k=1 +R=6378273, not WGS84). Default product: interpolated to a 0.125° latitude-longitude grid, 2880 × 321 (50°N–90°N) |
| Vertical resolution | The model uses 50 hybrid layers. Output is interpolated to 40 fixed depths from 0 to 4000 m (0, 2, 4, 6, 10, 15, 20, 25, 30, 40, 50, 60, 70, 80, 90, 100, 125, 150, 175, 200, 250, 300, 350, 400, 450, 500, 600, 700, 800, 900, 1000, 1200, 1400, 1600, 1800, 2000, 2500, 3000, 3500, 4000 m) |
| Temporal resolution | Daily, monthly and yearly means (surface momentum and heat fluxes — daily and monthly means) |
| Update frequency | Twice a year (extended in July and November) |
| Licence | Copernicus Marine Service licence — free of charge, including commercial use; attribution and citation of the product DOI required |
| Account required | Yes |
Note on period: 1990 is used for spin-up, and data start in 1991. As checked on 2026-10-12, the Copernicus Marine catalogue ends at 2026-06-30 for daily data and 2026-06 for monthly data. Per the PUM, the previous year is updated each July (near-real-time inputs) and again each November (reprocessed inputs), so the most recent period may later be regenerated with different inputs.
How to read the resolution
- Horizontal resolution lists both the native model grid and the grid you download, as given in the official documentation. The same dataset can come on different grids from different providers or products, so check the grid of your downloaded file.
- Grid spacing is not the effective resolution. A provided grid that is finer than the native one carries no extra information.
- Intervals and whether a value is instantaneous, mean or accumulated differ between variables — see 'Available variables' below.
- The 'provided resolution' in lists and the comparison table is a simplified km/degree value for quick comparison. Use this page for exact values.
Overview¶
A coupled ocean and sea-ice reanalysis of the Arctic Ocean and North Atlantic produced by NERSC for Copernicus Marine. It couples the HYCOM ocean model with the CICE sea-ice model and assimilates, in a weekly cycle with a 100-member Deterministic Ensemble Kalman Filter (DEnKF), satellite sea level anomalies, SST, sea ice concentration, sea ice drift and sea ice thickness (from 2010), plus in situ temperature and salinity profiles. From 1991 it provides daily, monthly and yearly temperature, salinity, currents, and sea ice concentration, thickness and drift on a 12.5 km grid, and is used to study Arctic sea-ice change, Atlantic Water inflow and freshwater storage.
Available variables¶
Variable names as given in the Copernicus Marine Product User Manual (PUM) and the dataset metadata. The daily, monthly and yearly datasets (cmems_mod_arc_phy_my_topaz4_P1D-m, _P1M, _P1Y) contain the same variables. On the original (polar stereographic) grid, vector components follow the grid x and y axes and must be rotated when converting to latitude-longitude.
Time-stamp convention — instantaneous, mean or accumulated? — click to show
All variables are means, not instantaneous values. Daily fields are daily means (centred at 12:00 UTC per the PUM), monthly fields are monthly means (centred mid-month per the PUM), and yearly fields are averages of the daily files; every value is the mean of the 100 ensemble members. Only the bathymetry (model_depth) is time-invariant.
- The PUM describes daily means as centred at 12:00 UTC and monthly means as centred mid-month, but the Copernicus Marine catalogue stamps daily data at 00:00 on the day and monthly data at 00:00 on the 1st. Read the time value only as a label for 'the mean of this day (month)'.
- Yearly means (…_P1Y) are computed from the daily files rather than the monthly ones, to avoid errors in nonlinear variables such as mixed layer depth.
- Data assimilation (DEnKF, 100 members) runs on a 7-day cycle; daily means average the model integration between analyses.
- Atmospheric forcing is 6-hourly ERA5.
instantaneous mean accumulated maximum minimum static — instantaneous: value at the time stamp · mean/accumulated/maximum/minimum: over the indicated interval · static: does not change with time
3D ocean (40 depths)¶
| Variable | Short name | Units | Time convention | Notes |
|---|---|---|---|---|
| Sea water potential temperature | thetao |
°C | day / month / year mean | |
| Sea water salinity | so |
1e-3 | day / month / year mean | |
| Eastward sea water velocity | vxo |
m/s | day / month / year mean | |
| Northward sea water velocity | vyo |
m/s | day / month / year mean |
2D ocean¶
| Variable | Short name | Units | Time convention | Notes |
|---|---|---|---|---|
| Sea surface height above geoid | zos |
m | day / month / year mean | |
| Ocean mixed layer thickness defined by sigma theta | mlotst |
m | day / month / year mean | |
| Sea water potential temperature at sea floor | bottomT |
°C | day / month / year mean | Interpolated 10 m above the seabed, so it can look discontinuous where layers meet |
| Ocean barotropic streamfunction | stfbaro |
m3 s-1 | day / month / year mean |
Sea ice¶
| Variable | Short name | Units | Time convention | Notes |
|---|---|---|---|---|
| Sea ice area fraction | siconc |
fraction, 0–1 | day / month / year mean | |
| Sea ice thickness | sithick |
m | day / month / year mean | Model value (uncorrected). Biased thin before 2010, when no thickness observations were assimilated |
| AI-corrected sea ice thickness | sithick_corr |
m | day / month / year mean | Pre-2010 period corrected by machine learning (Edel et al. 2025). Equal to sithick from 2010 onward |
| Surface snow thickness | sisnthick |
m | day / month / year mean | |
| Eastward sea ice velocity | vxsi |
m/s | day / month / year mean | |
| Northward sea ice velocity | vysi |
m/s | day / month / year mean |
Surface fluxes (separate datasets)¶
| Variable | Short name | Units | Time convention | Notes |
|---|---|---|---|---|
| Surface downward eastward and northward stress | tauuo, tauvo |
N m-2 | day / month mean | cmems_mod_arc_phy_my_mflux_P1D-m / _P1M-m (added 2024-11) |
| Net downward shortwave flux at sea water surface | rsntds |
W m-2 | day / month mean | cmems_mod_arc_phy_my_hflux_P1D-m / _P1M-m |
| Surface downward heat flux in sea water | hfds |
W m-2 | day / month mean | cmems_mod_arc_phy_my_hflux_P1D-m / _P1M-m |
Static fields¶
| Variable | Short name | Units | Time convention | Notes |
|---|---|---|---|---|
| Sea floor depth below sea level | model_depth |
m | static |
The full variable list is in the official documentation.
How to get the data¶
- Copernicus Marine Data Store — ARCTIC_MULTIYEAR_PHY_002_003 — Web download (Subset, File Browser) / Copernicus Marine Toolbox (copernicusmarine)
- Requires a free Copernicus Marine account. Dataset IDs — daily means cmems_mod_arc_phy_my_topaz4_P1D-m, monthly cmems_mod_arc_phy_my_topaz4_P1M, yearly cmems_mod_arc_phy_my_topaz4_P1Y (current version 202506); fluxes cmems_mod_arc_phy_my_mflux_P1D-m and _P1M-m, cmems_mod_arc_phy_my_hflux_P1D-m and _P1M-m. Each dataset is offered both interpolated to 0.125° latitude-longitude (default) and on the original 12.5 km polar stereographic grid (originalGrid).
- Copernicus Marine Toolbox documentation — Python / command line (CLI)
- Version 2.x is the supported release.
Example code¶
Install with pip install copernicusmarine, then run login() once to store your credentials.
import copernicusmarine
# Once only: copernicusmarine.login()
# Daily Arctic sea ice concentration and thickness (incl. AI-corrected) — default 0.125° lat-lon grid
copernicusmarine.subset(
dataset_id="cmems_mod_arc_phy_my_topaz4_P1D-m",
variables=["siconc", "sithick", "sithick_corr", "vxsi", "vysi"],
minimum_longitude=-180, maximum_longitude=180,
minimum_latitude=65, maximum_latitude=90,
start_datetime="2005-03-01", end_datetime="2005-03-31",
output_filename="topaz4b_seaice_200503.nc",
)
# For the original 12.5 km polar stereographic grid, add dataset_part="originalGrid"
# and give the area as minimum_x/maximum_x, minimum_y/maximum_y (units of 100 km).
Practical tips¶
- sithick jumps around 2010, when CryoSat-2 and SMOS thickness assimilation begins. For sea ice thickness or volume trends from 1991, use sithick_corr and state in your paper which variable you used.
- On the original grid, vectors (currents, ice velocity, wind stress) are aligned with the polar stereographic grid, and the PUM says they must be rotated when reprojecting to latitude-longitude. When loading original-grid files into GIS, convert the x and y coordinates from 100 km to metres and use the projection's sphere (radius 6378273 m).
- The default 0.125° interpolated product is slightly smoother than the original grid and has a slightly larger land mask. For straits and coastal transports, prefer the original grid.
- At about 12 km spacing the model cannot resolve Arctic mesoscale eddies (Rossby radius 1–2 km). Compare eddy-scale features with higher-resolution model output.
- Vertical interpolation (cubic splines) can overshoot where temperature and salinity gradients are strong, and a few cells near the seabed may be intermittently masked.
- GLORYS12 (global 1/12°) also covers the Arctic, but TOPAZ4b is an Arctic-specific reanalysis that assimilates sea ice drift and thickness directly. Comparing the two gives a sense of the uncertainty in sea ice analyses.
Notices and change history¶
2026-04-28 · Info · [MDSOP-377] Application of retention periods to ARCO datacubes
Copernicus Marine is applying the retention periods stated in the Product User Manuals to some ARCO datacubes. In the Subset GUI, Toolbox subset and MyOcean Viewer, the start date or temporal extent of some time series may change; original files (File Browser, copernicusmarine get) are not affected. The notice does not list the products concerned.
2025-07 (day not given) · Warning · Copernicus Marine Toolbox 1.x end of support
Toolbox v1 is no longer supported and should be replaced with v2. Version 2.4.0 dropped support for Python 3.9.
2025-06-25 · Info · New AI-corrected sea ice thickness variable (sithick_corr)
The reanalysis assimilates CryoSat-2 and SMOS sea ice thickness from 2010, which creates a discontinuity. A new variable, sithick_corr, corrects the pre-2010 thin bias with an AI technique (dataset version 202506); it equals the model-based sithick from 2010 onward.
2024-11-26 · Info · Surface momentum and heat flux datasets added
Four new datasets provide wind stress (cmems_mod_arc_phy_my_mflux_P1D-m, _P1M-m) and total heat fluxes (cmems_mod_arc_phy_my_hflux_P1D-m, _P1M-m). Pre-existing datasets remain unchanged.
2023-11-30 · Info · Arctic forecasting system upgraded to TOPAZ5 — reanalysis stays on TOPAZ4b
The analysis and forecast product (ARCTIC_ANALYSISFORECAST_PHY_002_001) was upgraded from TOPAZ4 to TOPAZ5, doubling the horizontal and vertical resolution. The reanalysis (002_003) is still produced with TOPAZ4b, so appending recent forecast-product data to it mixes two model systems.
References for your paper¶
Double-check before using in a paper
The citations below were compiled by the maintainer and may contain errors. Always double-check DOIs, authors, years and access-date formats against the provider's official dataset page and your journal's guidelines.
Dataset citation¶
Product DOI — DOI 10.48670/moi-00007
E.U. Copernicus Marine Service Information (CMEMS). Arctic Ocean Physics Reanalysis (ARCTIC_MULTIYEAR_PHY_002_003). Marine Data Store (MDS). DOI: 10.48670/moi-00007 (Accessed on DD-MMM-YYYY)
Show BibTeX
@misc{cmems_topaz4b,
author = {{E.U. Copernicus Marine Service Information (CMEMS)}},
title = {Arctic Ocean Physics Reanalysis ({ARCTIC\_MULTIYEAR\_PHY\_002\_003})},
year = {2020},
publisher = {Mercator Ocean International, Marine Data Store (MDS)},
doi = {10.48670/moi-00007},
note = {Accessed on DD-MMM-YYYY}
}
Reference papers¶
TOPAZ4 reanalysis quality assessment — DOI 10.5194/os-13-123-2017
Xie, J., Bertino, L., Counillon, F., Lisæter, K. A., Sakov, P. (2017): Quality assessment of the TOPAZ4 reanalysis in the Arctic over the period 1991–2013. Ocean Science, 13, 123–144.
Show BibTeX
@article{xie2017topaz4,
author = {Xie, Jiping and Bertino, Laurent and Counillon, Fran{\c{c}}ois and Lis{\ae}ter, Knut A. and Sakov, Pavel},
title = {Quality assessment of the {TOPAZ4} reanalysis in the {Arctic} over the period 1991--2013},
journal = {Ocean Science},
volume = {13},
pages = {123--144},
year = {2017},
doi = {10.5194/os-13-123-2017}
}
AI-corrected sea ice thickness (sithick_corr) — DOI 10.5194/tc-19-731-2025
Edel, L., Xie, J., Korosov, A., Brajard, J., Bertino, L. (2025): Reconstruction of Arctic sea ice thickness (1992–2010) based on a hybrid machine learning and data assimilation approach. The Cryosphere, 19, 731–752.
Show BibTeX
@article{edel2025sit,
author = {Edel, L{\'e}o and Xie, Jiping and Korosov, Anton and Brajard, Julien and Bertino, Laurent},
title = {Reconstruction of {Arctic} sea ice thickness (1992--2010) based on a hybrid machine learning and data assimilation approach},
journal = {The Cryosphere},
volume = {19},
pages = {731--752},
year = {2025},
doi = {10.5194/tc-19-731-2025}
}
SMOS sea surface salinity assimilation (cited in the PUM) — DOI 10.5194/os-19-269-2023
Xie, J., Raj, R. P., Bertino, L., Martínez, J., Gabarró, C., Catany, R. (2023): Assimilation of sea surface salinities from SMOS in an Arctic coupled ocean and sea ice reanalysis. Ocean Science, 19, 269–287.
Acknowledgements¶
Copernicus Marine recommends the acknowledgement "This study has been conducted using E.U. Copernicus Marine Service Information;" followed by the product DOI (10.48670/moi-00007). The PUM names no specific reference paper, so cite Xie et al. (2017) for the system and, if you used sithick_corr, Edel et al. (2025).
If this page helped, please acknowledge the Korean Climate Data Navigator
Please add the sentence below to the Acknowledgements of your paper.