# Polar Ocean Observations: A Critical Gap in the Observing System and Its Effect on Environmental Predictions From Hours to a Season

Gregory Smith, Richard Allard, Marcel Babin, Laurent Bertino, Matthieu

Chevallier, Gary Corlett, Julia Crout, Fraser Davidson, Bruno Delille, Sarah
Gille, et al.

► To cite this version:

Gregory Smith, Richard Allard, Marcel Babin, Laurent Bertino, Matthieu Chevallier, et al. Polar Ocean Observations: A Critical Gap in the Observing System and Its Effect on Environmental Predictions From Hours to a Season. Frontiers in Marine Science, 2019, 6, 10.3389/fmars.2019.00429. [https://doi.org/10.3389/fmars.2019.00429](https://doi.org/10.3389/fmars.2019.00429).

HAL Id: hal-02415050

https://hal.science/hal-02415050

Submitted on 10 Apr 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Polar Ocean Observations: A Critical Gap in the Observing System and Its Effect on Environmental Predictions From Hours to a Season

Gregory C. Smith¹*, Richard Allard², Marcel Babin³, Laurent Bertino⁴,
Matthieu Chevallier⁵,⁶, Gary Corlett⁷, Julia Crout⁸, Fraser Davidson⁹, Bruno Delille¹⁰,
Sarah T. Gille¹¹, David Hebert², Patrick Hyder¹², Janet Intrieri¹³, José Lagunas³,
Gilles Larnicol¹⁴, Thomas Kaminski¹⁵, Belinda Kater¹⁶, Frank Kauker¹⁷,¹⁸,
Claudie Marec³,¹⁹, Matthew Mazloff¹¹, E. Joseph Metzger², Calvin Mordy²⁰,
Anne O'Carroll²¹, Steffen M. Olsen²², Michael Phelps⁸, Pamela Posey⁸, Pierre Prandi¹⁴,
Eric Rehm³, Phillip Reid²³, Ignatius Rigor²⁴, Stein Sandven⁴, Matthew Shupe¹³,²⁵,
Sebastiaan Swart²⁶,²⁷, Ole Martin Smedstad⁸, Amy Solomon²⁸, Andrea Storto²⁹,
Pierre Thibaut¹⁴, John Toole³⁰, Kevin Wood²⁰, Jiping Xie⁴, Qinghua Yang³¹ and
the WWRP PPP Steering Group³²

OPEN ACCESS

**Edited by:**
Tong Lee,
NASA Jet Propulsion Laboratory
(JPL), United States

**Reviewed by:**
Agnieszka Beszczynska-Möller,
Institute of Oceanology (PAN), Poland
Gunnar Spreen,
University of Bremen, Germany

*Correspondence:
Gregory C. Smith
Gregory.Smith2@canada.ca

**Specialty section:**
This article was submitted to
Ocean Observation,
a section of the journal
Frontiers in Marine Science

**Received:** 02 November 2018

**Accepted:** 05 July 2019

**Published:** 06 August 2019

**Citation:**
Smith GC, Allard R, Babin M,
Bertino L, Chevallier M, Corlett G,
Crout J, Davidson F, Delille B, Gille ST,
Hebert D, Hyder P, Intrieri J,
Lagunas J, Larnicol G, Kaminski T,
Kater B, Kauker F, Marec C,
Mazloff M, Metzger EJ, Mordy C,
O'Carroll A, Olsen SM, Phelps M,
Posey P, Prandi P, Rehm E, Reid P,
Rigor I, Sandven S, Shupe M,
Swart S, Smedstad OM, Solomon A,
Storto A, Thibaut P, Toole J, Wood K,
Xie J, Yang Q and
the WWRP PPP Steering Group
(2019) Polar Ocean Observations:
A Critical Gap in the Observing
System and Its Effect on
Environmental Predictions From
Hours to a Season.
Front. Mar. Sci. 6:429.
doi: 10.3389/fmars.2019.00429

¹ Environmental Numerical Prediction Research Section, Meteorological Research Division, Environment and Climate Change Canada, Dorval, QC, Canada, ² Stennis Space Center, U.S. Naval Research Laboratory, Bay St. Louis, MS, United States, ³ Takuvik, UMI 3376, Université Laval-CNRS, Quebec City, QC, Canada, ⁴ Nansen Environmental and Remote Sensing Center, Bergen, Norway, ⁵ Division of Marine and Oceanography, Météo France, Toulouse, France, ⁶ CNRM, Météo France, CNRS, Université de Toulouse, Toulouse, France, ⁷ European Organisation for the Exploitation of Meteorological Satellites, Darmstadt, Germany, ⁸ Perspecta, Inc., Stennis Space Center, Bay St. Louis, MS, United States, ⁹ Northwest Atlantic Fisheries Centre, Fisheries and Oceans Canada, St. John's, NL, Canada, ¹⁰ Chemical Oceanography Unit, Université de Liège, Liège, Belgium, ¹¹ Scripps Institution of Oceanography, University of California, San Diego, La Jolla, CA, United States, ¹² Met Office, Exeter, United Kingdom, ¹³ Physical Sciences Division, NOAA Earth System Research Laboratory, Boulder, CO, United States, ¹⁴ Collecte Localisation Satellites, Toulouse, France, ¹⁵ The Inversion Lab, Hamburg, Germany, ¹⁶ Arcadis Nederland B.V., Zwolle, Netherlands, ¹⁷ Ocean Atmosphere Systems, Hamburg, Germany, ¹⁸ Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany, ¹⁹ Laboratoire d'Oceanographie Physique et Spatiale, UMR 6523, CNRS - IFREMER - IRD - UBO, Plouzané, France, ²⁰ Joint Institute for the Study of the Atmosphere and Oceans, University of Washington, Seattle, WA, United States, ²¹ European Organisation for the Exploitation of Meteorological Satellites, Darmstadt, Germany, ²² Danish Meteorological Institute, Copenhagen, Denmark, ²³ Bureau of Meteorology, Hobart, TAS, Australia, ²⁴ Polar Science Center, University of Washington, Seattle, WA, United States, ²⁵ Cooperative Institute for Research in Environmental Sciences, University of Colorado Boulder, Boulder, CO, United States, ²⁶ Department of Marine Sciences, University of Gothenburg, Gothenburg, Sweden, ²⁷ Department of Oceanography, University of Cape Town, Rondebosch, South Africa, ²⁸ Earth System Research Laboratory, National Oceanic and Atmospheric Administration, Boulder, CO, United States, ²⁹ Centre for Maritime Research and Experimentation, La Spezia, Italy, ³⁰ Woods Hole Oceanographic Institution, Woods Hole, MA, United States, ³¹ Guangdong Province Key Laboratory for Climate Change and Natural Disaster Studies, School of Atmospheric Sciences, Sun Yat-sen University, Zhuhai, China, ³² World Weather Research Programme (WWRP) Polar Prediction Project (PPP) Steering Group

There is a growing need for operational oceanographic predictions in both the Arctic and Antarctic polar regions. In the former, this is driven by a declining ice cover accompanied by an increase in maritime traffic and exploitation of marine resources. Oceanographic predictions in the Antarctic are also important, both to support Antarctic operations and also to help elucidate processes governing sea ice and ice shelf stability. However, a significant gap exists in the ocean observing system in polar regions, compared to most areas of the global ocean, hindering the reliability of ocean and sea ice forecasts. This gap can also be seen from the spread in ocean and sea ice reanalyses for polar regions which provide an estimate of their uncertainty. The reduced reliability of polar predictions may affect the quality of various applications
including search and rescue, coupling with numerical weather and seasonal predictions, historical reconstructions (reanalysis), aquaculture and environmental management including environmental emergency response. Here, we outline the status of existing near-real time ocean observational efforts in polar regions, discuss gaps, and explore perspectives for the future. Specific recommendations include a renewed call for open access to data, especially real-time data, as a critical capability for improved sea ice and weather forecasting and other environmental prediction needs. Dedicated efforts are also needed to make use of additional observations made as part of the Year of Polar Prediction (YOPP; 2017–2019) to inform optimal observing system design. To provide a polar extension to the Argo network, it is recommended that a network of ice-borne sea ice and upper-ocean observing buoys be deployed and supported operationally in ice-covered areas together with autonomous profiling floats and gliders (potentially with ice detection capability) in seasonally ice covered seas. Finally, additional efforts to better measure and parameterize surface exchanges in polar regions are much needed to improve coupled environmental prediction.

**Keywords:** polar observations, operational oceanography, ocean data assimilation, ocean modeling, forecasting, sea ice, air-sea-ice fluxes, YOPP

INTRODUCTION

Over the last 10 years, there has been a significant maturing of ocean prediction systems, led by efforts such as GODAE OceanView (GOV; Bell et al., 2015; Davidson et al., 2019) and the European Union Copernicus Marine Environmental Monitoring Service (CMEMS; Le Traon et al., 2017). Numerous operational global and regional ocean analysis and forecast systems are now in place providing services for a range of applications including search and rescue, short- and long-range atmospheric and coupled prediction systems, aquaculture, energy sector activities and environmental management including environmental emergency response.

We have also seen the implementation of high-resolution operational ice-ocean prediction systems for polar regions providing forecasts on timescales of hours to days. These include the CMEMS Arctic Monitoring and Forecasting Centre (ARC MFC), the U.S. Navy Global Ocean Forecasting System, Environment and Climate Change Canada's Global and Regional Ice Ocean Prediction Systems among others (Carrières et al., 2017). While these systems are intended to provide support for marine operations and related applications, there has been a growing acceptance of the importance of including coupled interactions across the marine surface in numerical weather prediction systems (Brassington et al., 2015). Indeed, operational medium-range weather forecasting systems in Canada and at the European Centre for Medium Range Forecasting (ECWMF) now include a dynamic coupling with ice-ocean models (Mogensen et al., 2017; Smith et al., 2018). As a result, polar ocean observations may now have impacts beyond the polar regions on mid-latitude weather predictions (Jung et al., 2015). Moreover, coupled atmosphere-ice-ocean models have been shown to provide skillful forecasts at monthly-to-seasonal time scales in the polar regions, especially for the Arctic sea ice cover (e.g., Guémas et al., 2016a; Sigmond et al., 2016). The

importance of polar sea ice initial conditions in affecting skill of monthly-to-seasonal predictions at lower latitudes in the atmosphere has also been discussed (e.g., Guémas et al., 2016b). These advances are also fostering a coupled modeling approach in the context of Earth system reanalyses for climate monitoring (e.g., Buizza et al., 2018).

Environmental prediction systems rely heavily on ocean observations from a variety of platforms (both in situ and from remote sensing). Indeed, Observing System Experiments (OSEs) have shown the delicate balance and complementarity provided by the current basket of observations (Oke et al., 2015; Fujii et al., 2019). However, polar regions present a number of unique observing challenges (see Calder et al., 2010). For short lead-times (hours to days), forecast skill depends strongly on initial conditions, which in turn depend on real-time observations ingested through data assimilation. The lack of real-time Argo profiling floats in polar regions due to the sea ice cover (Roemmich and Gilson, 2009) therefore creates a significant gap. This gap is compounded by additional errors in satellite products in polar regions, for example, due to difficulties in distinguishing snow and ice from clouds (Castro et al., 2016). The relative remoteness and harsh environmental conditions in polar regions further hinder efforts to provide *in situ* measurements. Moreover, the increasing use of coupled models noted above requires collocated observations of the atmosphere, ice and ocean, including flux estimates (Bourassa et al., 2013). This is all the more important because of the role of the Arctic region in particular in shaping the heat and freshwater transports at global scales, thus having a crucial remote impact on the mid- and low- latitude climate as well (Serreze and Barry, 2011; Jung et al., 2015).

Additionally, the rapidly receding summer ice cover in the Arctic is leading to both an increase in the demand for accurate ice-ocean and weather forecast products (Jung et al., 2016), as well as challenges in how to adapt prediction systems to
accurately forecast conditions that may not have an existing
analog in the historical record (e.g., the 2012 record low Arctic sea
ice cover; Parkinson and Comiso, 2013). Similar issues apply in
the Antarctic: rising temperatures and reports of both expanding
and reduced sea ice have led to uncertainties about changes in the
system (Maksym, 2019), while at the same time growing demands
from Antarctic tourism and research logistics have driven a need
for improved forecasts (Hoke et al., 2018).

Here, we outline the status of existing near real-time
ocean observational efforts in polar regions, discuss gaps,
and explore perspectives for the future. The first section is
dedicated to *in situ* observations of temperature and salinity
and several projects underway to further extend these efforts.
Section “Satellite Observations” focuses on satellite observations
of sea surface temperature (SST) and height (SSH), along
with information on the sea ice cover. Section “Air-Sea Flux
Measurements in Polar Oceans” discusses issues associated with
surface flux measurements in polar regions. Section “Forecasting
System Experiments” presents several impact studies in different
prediction systems demonstrating the importance of ocean and
sea ice observations to forecasting skill. Section “International
Efforts to Address Gaps in Polar Ocean Observations” discusses
several international projects working to address the gap in ocean
observations in the polar regions. Finally, recommendations
and an outlook for the future are presented in Sections
“Recommendations” and “Outlook.”

IN SITU OBSERVATIONS OF
TEMPERATURE AND SALINITY

In the polar oceans there are several methods for collecting
ocean and sea ice data and sending the data in near real-
time via the Global Telecommunications System (GTS) or other
transmission systems. In ice-free areas various methods are used,
including ships, Argo floats, surface drifters, gliders, etc., while
in ice-covered areas, only ice-borne instrument systems can
operate throughout the year and transmit data in near real-time.
Owing to their significant deployment and maintenance costs, the
number of these platforms is very limited, resulting in large gaps
in the data coverage in the ice-covered Arctic Ocean. Antarctic
sea ice is largely single-year (or first-year) ice, which makes it ill-
suited for ice based measurement systems. A significant fraction
of the ocean data from the Arctic, as well as under-ice Argo
profiles from the Antarctic seasonal ice zone, are only available in
delayed mode, because underwater platforms in ice-covered areas
cannot transmit data via satellite communication. This limits the
access to ocean data available in near real-time, which is required
by monitoring and forecasting services.

In the following section, the most common observing systems
for the Arctic Ocean and Antarctic marginal sea are described
briefly. Sections “The WHOI Ice-Tethered Profiler,” “ALAMO
and Arctic Heat,” “Deployment of Argo Floats in Canadian
Marginal Ice Zone by CONCEPTS,” “Under-ice BGC Argo
Floats in the Canadian Arctic,” and “Surface Drifters and IABP”
detail several specific observing efforts as exemplars of each
class of technology.

Overview of Current Observing System

Profiling Floats

Argo floats are the backbone of the global ocean observing
system developed over the last two decades with more than
3500 units currently in operation (e.g., Roemmich et al., 2009).
In the north, Argo floats are used in ice-free areas of the
North Atlantic/Nordic Seas, Baffin Bay and Bering Sea. In
the North Atlantic/Nordic Seas about 10–20 Argo floats are
deployed each year sustaining an operational array of about
40 – 50 Argo floats in these regions. These floats are provided
mainly by national efforts (e.g., NorArgo) and coordinated at
the European level by EuroArgo. The funding of these Argo
floats is relatively secure, and provided for by the countries that
participate in EuroArgo.

There is ongoing development to adapt Argo floats to operate
under ice, using ice-sensing algorithms as well as acoustic and
optical techniques. In the Southern Ocean, the Argo program
has demonstrated success with under-ice profiling floats that
delay data transmission until the float is in open water, typically
in summer (e.g., projects SOCLIM, RemOcean, SOCCO Bio-
argo and ReMOCA). In particular, the Southern Ocean Carbon
and Climate Observations and Modeling project (SOCCOM)
has made an effort to deploy a significant number of under-
ice Argo floats in the Southern Ocean. Moreover, RAFOS floats
(Klatt et al., 2007; Reeve et al., 2016) provide the capability
to triangulate their position using sound signals from nearby
moorings for periods when they are not able to surface. In the
Arctic, the ice-sensing algorithm is starting to become more
mature and is being complemented with additional detection
methods currently being tested by Takuvik (see Under-ice BGC
Argo Floats in the Canadian Arctic).

Ice-Borne Observing Systems

Many types of ice-borne ocean measurement systems have been
developed and fielded in the last 2–3 decades. Ice-based platforms
are the only autonomous systems that can presently deliver
near-real-time subsurface ocean observations year round from
ice-covered areas. Early examples consisted of discrete sensors
mounted on a floating mooring that drifts with the sea ice.
More recent examples employ a drifting surface buoy deployed
in the ice with a Conductivity-Temperature-Depth (CTD) sensor
that is repeatedly transported along the mooring cable over the
upper 800–1000 m of the water column. Ice-Tethered Profiler
(ITP) systems from the Woods Hole Oceanographic Institution
(WHOI, see The WHOI Ice-Tethered Profiler for details) have
been fielded since 2004. Similarly, investigators from Laboratoire
d'Océanographie et du Climat (LOCEAN) have developed and
deployed IAOOS (Ice-Atmosphere-Arctic Ocean Observation
System) buoys<fnref n="1" />. In addition to CTD profile observations, the
IAOOS buoy includes ice mass balance measurements (snow
and sea ice temperature and thickness), as well as a variety
of meteorological sensors (Provost et al., 2015). Other similar
platforms include the JAMSTEC Polar Ocean Profiling system
(POPS; Kikuchi et al., 2007), the NPS Autonomous Ocean Flux

<footnote n="1">¹See http://iaoos.ipev.fr/ for more details.</footnote>
Buoy<fnref n="2" /> (AOFB), and the UpTempO ice-tethered buoy developed by APL-UW (Castro et al., 2016). The Arctic University of Norway has recently developed and field tested a new platform called the Ice-Tethered Platform Cluster for Optical, Physical, and Ecological Sensors (ICE-POPE) that will be deployed operationally in 2019 (Berge et al., 2016). The Polar Research Institute of China is also developing ice-tethered platforms, which are being tested during Arctic expeditions of the I/B Xuelong. In comparison to Argo-type profiling floats, ice-based platforms are rather expensive (>100–300 k€ per unit). To date, funding to construct and field these systems and provide data to the community has derived from individual PI-led research projects. There are no operational programs to support the ice-tethered platforms at the moment, but there is a clear requirement from the Arctic Ocean modeling and data assimilation community that 10+ such platforms need to be deployed every year and supported operationally.

It is widely known that sea ice in the Arctic is shrinking in areal coverage, thinning, and becoming more mobile. All of these changes present complications to an ice-based observing system. Although diminished, the sea ice will remain critically important to earth's climate and ecosystems as well as transportation and tourism, making ice-following observing platforms necessary into the future. However, future ice-borne instrument systems must be able to float and demonstrate resilience during fall freeze-up. Thinner, more mobile ice can be more prone to ridging that can damage ice based buoys. It has not proven feasible to maintain the array of 20 ice-based observatory systems in the Arctic that was envisioned at the turn of the century. Nor have these technologies been used extensively in the seasonal ice zone surrounding Antarctica, where sea ice typically melts completely every summer. Nevertheless, ice-based observatories have and are continuing to return valuable year-round upper ocean data from the central Arctic. Buoy clusters sampling various elements of the atmosphere, sea ice and upper ocean have proven particularly valuable.

## Ice/Snow Surface Drifters

A reasonable number of low-cost ice buoys operate in the Arctic, providing ice motion, air pressure and surface temperature data. At present, most of these buoys are drifting in the western part of the Arctic, whereas the eastern part, including the Eurasian Basin and the Russian shelves, has very few buoys. The drifter data are transmitted via the Argos or Iridium satellite systems then posted on the WMO/IOC GTS and provide baseline data for weather and ice forecasting in the Arctic (see Section Surface Drifters and IABP for more detail). Drifters are also routinely deployed throughout the Southern Ocean, providing important information (e.g., SST and sea level pressure) to forecast systems from extremely data sparse regions.

Sea ice mass balance buoys (IMBs) provide information on temperature within snow and ice as well as ice thickness and its motion. Development and use of these buoys have continued for several decades with more than 100 IMBs deployed

from 1957 to 2014. The data are transmitted in near real-time, but processing of the temperature profiles is mainly done manually. There are ongoing efforts to develop automated methods for retrieval of snow and ice thickness (Liao et al., 2018; Zuo et al., 2018). It is a challenge to develop robust algorithms that give accurate retrievals through the melting and freezing seasons.

Currently, three different systems are deployed in the Arctic including the CRREL-Dartmouth IMB (Richter-Menge et al., 2006), the SRSL Sea Ice Mass Balance Array (SIMBA; Jackson et al., 2013) and the TUT ice-tethered buoy developed by PRIC and TUC in China (Zuo et al., 2018). Some of the Ice-tethered platforms described in Section “Ice-Borne Observing Systems” are also equipped with SIMBA instruments.

## Other Systems

(1) Ferrybox systems: The Norwegian Institute for Water Research (NIVA) operates a ferrybox between Tromsø and Longyearbyen/Ny Ålesund. Data are obtained every 2 weeks year-round. This system is fairly sustainable through support from a new infrastructure project led by NIVA.

(2) Ship-based CTD data during scientific cruises and fishery management cruises: Several research vessels in the North Atlantic/Nordic seas deliver CTD observations in near real-time. The data are provided more or less regularly, depending on the schedule of the research vessels.

(3) Profiling gliders: Several institutions conduct ocean glider experiments in the Arctic (e.g., Nordic Seas, David Strait, Fram Strait, north of Svalbard, Beaufort Sea) providing CTD data in near real-time. These experiments are mainly in the summer season, however, acoustically navigated ice-capable Seagliders (developed by APL-UW) were successfully used for the year-round missions in Davis Strait and in the Beaufort Sea including under ice measurements in winter (Curry et al., 2014; Lee et al., 2016). Glider experiments are also conducted with growing frequency in the Southern Ocean. These provide valuable data and upper ocean process understanding from winter through summer (e.g., Swart et al., 2015; Thomson and Girton, 2017; du Plessis et al., 2019 – in review).

(4) Autonomous surface vehicles: Wave Gliders have been deployed in both the Arctic and Southern Ocean, providing rare surface flux (heat, momentum and CO₂) information over multiple months at a time and with high resolution (Monteiro et al., 2015; Schmidt et al., 2017; Thomson and Girton, 2017). Many of these platforms are reporting or being adapted to report real-time surface observational data.

(5) Sea-mammal borne instruments (CTD data) are becoming the most numerous *in situ* observations at high latitude in the Antarctic and also more recently in the Arctic, and may improve forecast skill and circulation patterns (Fedak, 2013; Roquet et al., 2013; Carse et al., 2015).

(6) Acoustic tomography: Acoustic sources and receivers have been deployed in the Arctic in the past (Mikhalevsky et al., 2015) and offer the unique capability to measure

<footnote n="2">²See https://www.oc.nps.edu/~stanton/fluxbuoy/ for more details.</footnote>
large-scale changes in temperature and heat content of the Arctic Ocean over long periods (Howe et al., 2019). New deployments of this technology were conducted under the CANAPE project (led by Scripps); CAATEx programs (led by NERSC) have been approved. While no observations are currently available in real-time, nor assimilated in delayed mode, acoustic tomography nonetheless presents an interesting possibility for constraining ocean state-estimate/prediction systems.

(7) Moored Arrays: Several long-term observatories have been running in both ice-covered and open-water areas with the capacity to provide real-time data. Examples include the German long-term observatory FRAM (Frontiers in Arctic Marine Monitoring; Soltwedel et al., 2013) and the Barrow Strait Observatory (Richards et al., 2017).

## The WHOI Ice-Tethered Profiler

Historically, the Arctic Ocean has been poorly sampled relative to waters at low- and mid-latitudes, especially in winter time. To address this observing shortfall, the WHOI ITP was designed to sample the upper ocean below drifting sea ice throughout the year and to return data in near real time (Krishfield et al., 2008; Toole et al., 2011). As noted earlier, there are now a variety of systems with similar capability now being fielded; the ITP is detailed here as exemplar of the group. The expendable ITP consists of a surface buoy that supports a weighted wire-rope tether extending through the ice and down to (at most) 800 m depth (Figures 1A,C,D). The heart of the ITP system is a cylindrical vehicle fitted with sensors (similar in size and shape to an Argo float) that travels up and down the tether. ITPs are equipped with CTD instruments for observing the ocean's thermohaline stratification and may also include sensors that sample, for example, dissolved oxygen (Timmermans et al., 2010), bio-optical properties (Laney et al., 2014), and ocean currents (Ice-Tethered Profiler with Velocity - ITV-V: Thwaites et al., 2011; Cole et al., 2015). In addition, instruments measuring temperature-conductivity, pCO₂, dissolved O₂ and pH have been affixed to the tether above the profiling interval (e.g., Islam et al., 2017). Deployments may be done from ice camps (supported by fixed-wing aircraft or helicopters) or ships. The majority of deployments have been through holes augured through ice floes but a handful of systems have been installed in open water (the buoy has sufficient buoyancy to support the system); most of those have survived fall freeze-up.

The basic ITP system was designed for an operational lifetime of more than 2 years assuming approximately 1500 m of profiling per day (e.g., 2 one-way profiles of 750 m span). Actual lifetimes of the full ITP system are often less than this, **Figure 1B**. There are two major failure modes of ITPs: crushing of the surface buoy and/or breaking of the tether in ice ridging events and dragging of the tether in shallow water. As is evident in **Figure 1B**, ITP surface buoys frequently transmit position data for an extended time after communication with the underwater unit is lost.

ITP data, available from the project website<fnref n="3" />, the National Environmental Data Center and the Arctic Data Center, support a

range of scientific investigations and student projects. The basin-wide and year-round coverage facilitates studies of seasonal-to-interannual physical and biogeochemical processes (e.g., Rabe et al., 2011; McPhee, 2013; Laney et al., 2014, 2017; Islam et al., 2017) and basin-scale phenomena (e.g., Timmermans et al., 2014), as well as supporting the initialization (data assimilation) and/or validation of numerical models. Smaller-scale processes may also be investigated with ITP data, including meso- and sub-mesoscale variability (e.g., Timmermans et al., 2011; Zhao et al., 2014, 2016), near-inertial internal waves (Cole et al., 2014; Dossier et al., 2014) and double diffusion (e.g., Timmermans et al., 2008; Shibley et al., 2017). Notably, the range of sensors able to be supported on ITPs provides a wide-ranging view of the evolving Arctic Ocean system.

## ALAMO and Arctic Heat

Since 2014, NOAA Pacific Marine Environmental Laboratory (PMEL) and the University of Washington's Joint Institute for the Study of the Atmosphere and Ocean (JISAO) have been developing and field testing a series of autonomous marine systems designed for Arctic applications. This work has been carried out in collaboration with scientists, forecasters, and industry, primarily under the auspices of the Innovative Technologies for Arctic Exploration (ITAE<fnref n="4" />) program. Fifteen new types of technology have been developed or advanced for Arctic use under this program, including the Air-Launched Autonomous Micro-Observer (ALAMO), the Saildrone autonomous marine vehicle, and the Oculus Coastal Glider.

An example of a multi-platform experiment along these lines is the Arctic Heat Open Science Experiment (Wood et al., 2018). Here a NOAA Twin Otter research aircraft (NOAA-56) has been equipped with a range of weather and ocean-sensing instruments, which in 2018 included flight-level weather and radiometry, an A-size sonobuoy deployment tube used for a range of air-deployed probes – including AXBT/AXCTD/AXCP, atmospheric dropsondes, ALAMO autonomous profilers, and experimental UAS (Unmanned Aerial Systems) platforms – LIDAR and thermal imaging camera. Arctic Heat deploys an array of ALAMO floats in the Chukchi Sea between May and September; floats deployed later in the season are active into winter (Figure 2), and are used to develop an experimental seasonal freeze-up projection, in combination with satellite-derived SST and historical sea ice extent/concentration. In 2018 a combination of more than two dozen profiling floats of various types were deployed by three collaborating research groups in the Chukchi and Beaufort Seas<fnref n="5" />.

The dense array of floats in 2018 provides an opportunity to thoroughly investigate the potential gain in predictive skill provided by enhanced real-time ocean observations. Understanding how assimilating more observations will impact modeled analyses and short-term forecasts is of fundamental importance as new coupled models are developed. For example, sensitivity experiments with the NOAA-Earth System Research

<footnote n="3">³http://www.whoi.edu/itp</footnote>

<footnote n="4">⁴https://www.pmel.noaa.gov/itae/</footnote>

<footnote n="5">⁵https://www.pmel.noaa.gov/arctic-heat/data</footnote>
A

<figure><img src="image_3.png" /></figure>

C

<figure><img src="image_4.png" /></figure>

D

<figure><img src="image_5.png" /></figure>

B

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

**FIGURE 1** | (A) Schematic drawing of the WHOI Ice-Tethered Profiler (ITP) system; (B) Histogram of ITP underwater vehicle lifetimes (top) and (bottom) the periods (shown as black vertical bars) over which telemetry was received from each ITP underwater unit and from each corresponding surface buoy (black plus gray bars). The history of ITP systems deployed in the Southern Ocean and in lakes are excluded from this plot. (C) Schematic drawing of the bio-optical ITP sensor suite with CTD/O₂, chlorophyll fluorescence, CDOM, optical backscatter and PAR (the latter suite housed under a retractable shutter), and (D) installation photograph of an ITP with a Modular Acoustic Velocity Sensor (ITP-V).

Laboratory Coupled Arctic Forecast System (CAFS) (see NOAA-ESRL/CIRES Coupled Arctic Forecast System) are being pursued.

## Deployment of Argo Floats in Canadian Marginal Ice Zone by CONCEPTS

A Government of Canada initiative called CONCEPTS (Canadian Operational Network of Coupled Environmental Prediction Systems; Davidson et al., 2013) has implemented global and regional environmental prediction systems over the last few years. A particular need has been identified to complement the existing observing system with additional real-time observations to better constrain water mass properties in the ocean analyses, in particular in Canadian seasonally ice-infested seas. As a result, a project is underway to test the implementation of different observing technologies for their ability to provide reliable measurements of ocean temperature and salinity at a low cost. A first step was the deployment of a standard Argo float in the Gulf of St. Lawrence, configured with a shallow profiling and parking depth (200 m). The objective was to provide observations over the spring-to-fall ice-free period and to allow the float to remain subsurface in winter when ice forms. The use of existing technology (i.e., without any particular ice detection capability) reduced costs and increased feasibility. This initial

effort was quite successful and provided an excellent complement to existing (and more costly) moorings deployed in the Gulf of St. Lawrence.

A second effort included the deployment of an Argo float on the Labrador shelf, again with a shallow profiling depth of 200 m. This depth was greater than the local bathymetry allowing the float to rest on the sea floor and reduce lateral drift. Despite the presence of ice in winter, this float was able to breach the surface periodically (roughly every 10 days) and transmit ocean measurements. Additionally, the float observed temperature and salinity measurements during winter that deviated significantly from climatological conditions (Figure 3). As these climatological values are typically used by various applications in the absence of other data, the Argo float observations filled an important gap. Moreover, their availability in real time permits their use in data assimilation systems that can allow the detected anomaly in water mass conditions to be propagated over suitably correlated water masses (i.e., along the Labrador Shelf).

As part of the YOPP (see Year of Polar Prediction), this effort has been expanded to include the deployment of 7 Argo floats during the Arctic summer Special Observing Period (July–September 2018). If successful, and found to provide adequate benefit for the cost, this effort could become part of the Canadian Argo Program, and also contribute to the newly forming
<figure><img src="image_8.png" /><figcaption>**FIGURE 2** | Temperature and salinity profiles collected by ALAMO 9119-CTD from September 17 to December 8, 2017. The float began sampling near 167W, 70N and the last profile was near 165W, 72N. See: https://www.pmel.noaa.gov/arctic-heat/ for more information, including the float track.</figcaption></figure>

Canadian Integrated Ocean Observing System (CIOOS). The use of air deployable ALAMO floats and floats with ice-detection capabilities is also being investigated.

## Under-Ice BGC Argo Floats in the Canadian Arctic

As part of project Green Edge<fnref n="6" />, Takuvik deploys Biogeochemical Argo floats (BGC Argo; manufactured by NKE<fnref n="7" />) in the Arctic Ocean to study the dynamics of ice-edge phytoplankton spring blooms as controlled by sea ice dynamics, vertical mixing, light and nutrients. BGC Argo floats are more specialized profilers, equipped with sensors capable of sampling additional essential ocean variables (Claustre et al., 2010). The study focuses on Baffin Bay, which involves navigational challenges for floats in terms of bathymetry, ice coverage and circulation. Simulations of trajectories to choose the best dropping zones, combined with observations from ice charts (climatology and real-time charts) are necessary for safe deployments. An ice-covered ocean presents a real challenge for Argo floats that must surface for geo-localization and to use satellite networks for data transmission and command reception (Riser et al., 2016), see **Figure 4**. Therefore, some technical refinements have been needed to make the floats operational in the Arctic Ocean (Fennel and Greenan, 2017).

The development of a new generation of floats (PRO-ICE) to be operated under ice, was founded by the French project

NAOS<fnref n="8" />. If sea ice is present at the surface, Argo floats need to postpone surfacing. It will perform several consecutive profiles without sending stored data. PRO-ICE have additional non-volatile memory and are able to record data from all profiles including those performed during wintertime. Because of the need to transmit a high volume of data, communications are performed via a two-way Iridium link, which is faster than the older ARGOS link. This minimizes time at the surface and allows instructions to be transmitted to the float. In addition, floats to be deployed in Arctic seas need to be able to deal with a wide range of seawater salinity/density. Some Argo floats, like PRO-ICE, do not require pre-ballasting as a function of mean seawater density. This is a huge advantage since they can surface easily in a low-density surface seawater layer (Riser et al., 2016).

The PRO-ICE floats use a combination of three technologies to detect ice: a reversed altimeter (active acoustics), an Ice Sensing Algorithm (Klatt et al., 2007) (ISA) based on seawater freezing temperature and an optical sensor. In the Arctic Ocean, the threshold of the ISA algorithm is difficult to determine because of the high levels of variation in salinity (seasonal and regional), compared to the homogeneous salinity present in Antarctica. Takuvik gathered a substantial data base of temperature and salinity data in Baffin Bay, linked to ice cover information in order to locally adapt the threshold of ISA.

Furthermore, the resolution provided by an altimeter is not sufficient to detect a thin-ice layer. For this, Takuvik has

<footnote n="6">⁶http://www.greenedgeproject.info/</footnote>

<footnote n="7">⁷http://www.nke-instrumentation.com/</footnote>

<footnote n="8">⁸Naos, “Equipex Naos.” http://www.naos-equipex.fr/.</footnote>
<figure><img src="image_9.png" /><figcaption>**FIGURE 3** | Measurements from an Argo float deployed on the Labrador Shelf. (A) Show the location of transmissions from the float over the period 01-August-2017 to 01-June-2018. (C) Show a Synthetic Aperture Radar image from RADARSAT-2 for 28-Jan-2018 with a blue star indicating the location of the Argo float. (B,D) Present analyses (orange) and 5-day forecasts (blue) from the Global Ice Ocean Prediction System for temperature and salinity respectively. Also shown are values from the World Ocean Atlas 2013 climatology (green). The presence of sea ice is indicated by the dashed red line, with values near the bottom indicating no ice and values near the top of the panel indicating the likely presence of ice as detected by GIOPS ice analyses.</figcaption></figure>

developed an ice-detection system based on laser polarimetry (Lagunas et al., 2018). Since sea-ice is a strong light depolarizer, this characteristic was used as an indicator of the presence or absence of ice. This ice-detection system has been installed on several PRO-ICE BGC-Argo floats deployed yearly in Baffin Bay since 2016.

A potential solution to data transmission issues caused by an inability to surface under the ice cover may be provided by underwater acoustic networks. Multipurpose underwater acoustic networks can provide under-ice positioning (“underwater GPS”; Lee and Gobat, 2008) and low bit-rate

communication services (Freitag et al., 2015) to AUVs fitted with low-power acoustic receivers. Acoustic positioning systems have been deployed for under-ice navigation of Argo floats in the Weddell Sea (Klatt et al., 2007). Challenges for designing basin-scale acoustic systems include modeling the time-varying nature of the Arctic sound channel, for which measurements of annual variations of under-ice sound profiles and ice-bottom roughness are scarce (Rehm et al., 2018). The role of multipurpose acoustics networks as key components of polar observing systems is explored further in a complementary article (Howe et al., 2019).
<figure><img src="image_10.png" /><figcaption>**FIGURE 4** | A daily profile cycle of BGC-Argo floats deployed during the 2016 Green Edge scientific mission in Baffin Bay. The main goal is the understanding of the dynamics of the phytoplankton spring bloom and determine its role in the Arctic. During the spring-period the risk of colliding with sea-ice when emerging is a threat to the security of the floats. Moreover, during wintertime, geo-localization and the use satellite networks for data transmission and commands reception is not yet possible (Credit: J. Sansoulet, Takuvik).</figcaption></figure>

Data acquired by Takuvik's PRO-ICE floats are available in the Argo database at the CORIOLIS Global Data Assembly Centre<fnref n="9" /> (GDAC). The data policy set for Argo mandates free access to the data for all interested scientists, research groups and operational agencies. The data management is designed to facilitate easy and immediate access to data not only in real time, but also in delayed mode.

## Surface Drifters and IABP

Surface air pressure, temperature and ocean/ice circulation in the polar regions (Figure 5) are observed by a network of drifting buoys maintained by participants of the International Arctic Buoy Programme (IABP<fnref n="10" />, Figure 6) and the International Programme for Antarctic Buoys (IPAB<fnref n="11" />). Over the Arctic Ocean and its peripheral seas (north of 60°N), there were about 120 buoys reporting as of August 2018. Most of these buoys were located in the North American sector, with only a handful reporting in the Eurasian sector of the Arctic. This gap in the Arctic Observing Network creates a significant uncertainty in the analyzed fields of sea level pressure, temperature and winds (Inoue et al., 2009, Figure 7).

In the Southern Ocean (south of 40°S) there were about 50 buoys as of August, 2018. Although over 100 buoys were deployed during the Austral summer near the coast of Antarctica, most of

these buoys were destroyed by the sea ice during freeze up or have been blown north, away from the coast by the prevailing winds.

Maintaining the polar observing networks near the coasts around most of Antarctica, and in the Eurasian Arctic is an ongoing challenge, since the prevailing winds and ocean currents quickly transport the buoys away from the coast. Ideally, the IABP and IPAB networks should be reseeded during the winter, but it is difficult to deploy buoys from ships and aircraft during the polar night and extreme cold and harsh conditions of winter.

The participants of the IABP and IPAB strive to release their data onto the WMO/IOC GTS in near real-time by both the global research and operational weather and ice forecasting communities.

## SATELLITE OBSERVATIONS

Satellite observations of polar oceans have been acquired for more than four decades by different measuring systems. The most prominent record comes from microwave radiometers that have been documenting the reduction of sea ice extent (Cavalieri and Parkinson, 2012; Osborne et al., 2018). Remote sensing systems also provide valuable information about other parameters of the sea ice (concentration and thickness) or about the underlying ocean (temperature and circulation). However, despite continuous technological advances and progress in data processing, several weaknesses can be pointed out and should be considered in the future for improving the observing system of polar regions.

<footnote n="9">⁹CORIOLIS, “Argo float data and metadata from Global Data Assembly Centre (Argo GDAC). Ifremer.” http://doi.org/10.17882/42182.</footnote>

<footnote n="10">¹⁰http://IABP.apl.uw.edu</footnote>

<footnote n="11">¹¹http://www.ipab.aq</footnote>
<figure><img src="image_11.png" /><figcaption>**FIGURE 5** | Map of drifting buoys reporting on the WMO/IOC GTS in August 2018. Source: http://OSMC.NOAA.GOV.</figcaption></figure>

## Sea Ice Concentration

Satellite-based sea ice observations are required for assimilation to provide accurate polar environmental forecasts. These observations have been proven beneficial in improving prediction skill at different temporal scales when assimilated (e.g., Lisaeter et al., 2003; Lindsay and Zhang, 2006; Stark et al., 2008; Massonnet et al., 2013; Tietsche et al., 2013). In addition, it has also been shown that summer sea ice thickness (SIT) can be constrained to some extent by assimilating sea ice concentration (e.g., Yang et al., 2015). However, there is a significant spread in sea ice concentration products obtained through different retrieval algorithms (Ivanova et al., 2014), which affects the consistency of ocean-sea ice analyses that assimilate those products (Chevallier et al., 2016; Uotila et al., 2018), and the skill of seasonal predictions initialized from those reanalyses (e.g., Bunzel et al., 2016).

## Sea Ice Freeboard and Thickness

Observing SIT from space is a challenge (Kwok and Sulsky, 2010), and gridded data are sparse. Recently SIT has been estimated using altimetry (CryoSat-2; Kwok and Cunningham, 2015), or the L-band radiometers of SMOS (Kaleschke et al., 2012; Tian-Kunze et al., 2014) and NASA's SMAP (Pałilea et al., 2019) satellites, for thinner ice. The ICESat-2 laser altimeter satellite, which was launched in September 2018, is also starting to provide polar SIT estimates.

Sea ice thickness is typically derived from radar altimetry observations of freeboard using the hydro-static equilibrium assumption (e.g., Alexandrov et al., 2010). Both the freeboard

retrieval from altimeter measurements (e.g., Ricker et al., 2014) and the freeboard to thickness conversion are active fields of research. The freeboard to thickness conversion is constrained by the limited availability of auxiliary information (sea ice type, density and snow parameters). Snow thickness on sea ice is still poorly known and is a major source of uncertainty. Snow depth models are not yet satisfactory, and alternative strategies will need to be defined to permit improving snow thickness estimation over sea ice. The joint use of Ku and Ka frequencies may facilitate estimates of the snow load above the ice pack. Studies based on AltiKa (Ka-band about 35.7 GHz) and CryoSat-2 (Ku-band about 13.5 GHz) satellite data have shown that differences in penetration of Ka- and Ku-band are correlated with snow loading on sea ice (Armitage and Ridout, 2015; Guerreiro et al., 2016). Although AltiKa only provides measurements up to 81.5°N and thus cannot provide pan-Arctic data. The combination of laser (ICESat-2) and radar altimetry is also promising for this purpose. A better estimation of freeboard and then thickness would greatly benefit from such measurement complementarity. Data editing is also important on heterogeneous surfaces such as sea ice, where melt ponds act as bright targets in the radar footprint resulting in peaky waveforms that look very similar to returns from leads. As a result, radar altimeter ice thickness products based on CryoSat-2 (e.g., Hendricks et al., 2016) are not available in summer months.

Previous studies showed that assimilation of SMOS ice thickness significantly improves the first-year ice estimates (Yang et al., 2014; Xie et al., 2016). Furthermore, assimilating CryoSat-2 and SMOS SIT leads to a reliable pan-Arctic SIT estimates (Mu et al., 2018; Xie et al., 2018), and also has the potential to
<figure><img src="image_12.png" /><figcaption>**FIGURE 6** | Map of drifting buoys reporting in the Arctic on August 3, 2018. Source: http://IABP.apl.uw.edu.</figcaption></figure>

improve seasonal forecasts of Arctic sea ice (Chen et al., 2017; Blockley and Peterson, 2018). The Cryorad mission has been proposed by Macelloni et al. (2017) as a microwave radiometer at very low frequency (down to 500 MHz, lower than SMOS) which longer wavelengths would penetrate through thicker SIT. A critical difficulty related to this mission is the contamination by Radio-Frequency Interference.

## Sea Surface Temperature

Satellite SST retrievals in polar regions are challenging. SST is estimated from instruments operating in the infrared (IR) and microwave regions of the spectrum through the so-called atmospheric window regions. IR retrievals in polar regions are often limited by the observed abnormal (often very dry) atmospheric conditions (Vincent et al., 2008). In addition, there are issues in detecting clear-sky ocean conditions that are free from cloud and sea ice, issues which are compounded during polar nights and in areas of persistent cloud. Consequently, frequent measurements of SST at high latitudes rely on microwave imaging instruments. Although not impacted by cloud (unless precipitation), microwave retrievals are impacted by issues in detecting sea ice, especially at the edges of the instrument footprint.

Many SST retrieval algorithms for both IR and microwave rely on *in situ* data to account for both deficiencies in their calibration as well as correcting for atmospheric attenuation. The lack of *in situ* data should be addressed and new innovative approaches are needed (e.g., Castro et al., 2016). The Group for High Resolution Sea Surface Temperature (GHRSSST) coordinates the production of multi-satellite merged SST products. However, the lack of accurate satellite and *in situ* data means there is little consistency between products, especially at high latitudes (Dash et al., 2010).

The current lack of continuity of microwave imagers that can be used to derive global SST is a major concern. For polar regions this requires the inclusion of a channel around 6.9 GHz (Gentemann et al., 2010). Currently, the only approved future instrument with this capability is the Chinese Microwave Radiometer Imager (MWRI) onboard the HaiYang-2B (HY-2B). A Copernicus Imaging Microwave Radiometry (CIMR) is currently being studied by the European Space agency (ESA) and JAXA is planning a follow-on to the Advanced Microwave Scanning Radiometer (AMSR2). The AMSR3 and CIMR missions are highly complementary and in combination would provide improved coverage and sampling in polar regions.
<figure><img src="image_13.png" /></figure>

**FIGURE 7** | Standard deviation (SD) of sea level pressure measurements from various atmospheric reanalyses. The SD is low in areas where there are buoy observations (A). The spread increases to cover the whole Arctic when the observations from the buoys are removed from the reanalyses (B) (Inoue et al., 2009).

## Sea Surface Height

Satellite observations of sea level are required to constrain surface geostrophic currents in ocean forecasting systems, and several teams try to tackle the issue (Prandi et al., 2012; Andersen and Piccioni, 2016; Armitage et al., 2016). However, sea level observations from altimetry over polar regions suffer from three main issues:

1. First, the altimeter constellation has mainly been created to fulfill ocean requirements for the ice-free regions, in particular with regard to orbit coverage/inclination. With the exception of the CryoSat-2 mission, which covers the Arctic Ocean up to 88°N, altimetry missions do not cover poleward of 82°, leaving a vast region without any measurement.

2. Second, although significant progress has been made to distinguish whether measurements correspond to open ocean, ice floes or leads within the sea ice, further progress is still needed to unambiguously identify the different returns, in particular in complex mixed water/floe areas. Exploitation of close match-ups between SAR imagers and altimeter measurements as the potential for improving the identification of leads and the editing of ambiguous measurements coming from melt ponds or polynya in the melt season. An example of such collocation between a Sentinel-1 image and Sentinel-3 measurement is provided in **Figure 8** (Longépé et al., 2019).

3. Lastly, the accurate retrieval of absolute SSHs, and therefore currents, in polar regions also suffers from degraded corrections. Tide models show higher

<figure><img src="image_14.png" /></figure>

**FIGURE 8** | Collocation of one Sentinel-1 SAR image (background) and Sentinel-3 altimeter waveforms (Unfocused processing; color) over a lead in the Arctic Ocean.

inter-model variability in polar regions than anywhere else, mean sea surfaces are not as accurate. An effort to refine the geophysical corrections applied to altimeter measurements is needed to improve polar SSH accuracy levels and make them useful for assimilation in operational oceanography systems.

As a result, there is currently no dedicated operational Arctic sea level product for assimilation in models. Efforts to provide Arctic sea level information in CMEMS, following the current SL-TAC products are ongoing.
Other innovative concepts such as SKIM's (Sea surface KInematics Multiscale monitoring, Arduuin et al., 2018) rotating radars at different angles offer opportunities to monitor waves, surface currents and possibly sea ice drift. The Surface Water Ocean Topography (SWOT) mission also has the potential to provide innovative new observations for sea level and sea ice cover, although its orbit is lower than 70 degrees.

## Sea-Ice Drift

Sea ice drift data are now obtained all year round both in the Arctic and Antarctic by pattern cross-correlation of scatterometers and passive microwave images. For global information, Scanning Multichannel Microwave Radiometer (SMMR), the Special Sensor Microwave/Imager (SSM/I) and SSMI/Sounder (SSMIS), and the Advanced Very-High-Resolution Radiometer (AVHRR) are used. Currently, daily products (e.g., Lavergne et al., 2010; Tschudi et al., 2019) typically use data acquired at Day D-1, with span of 24 or 48 h. In addition, buoy observations of the International Arctic Buoy Program (IABP), and ice motion derived from NCEP/NCAR surface wind vectors can be used. Low resolution ice drift products are calculated daily from aggregate charts derived from radiometers (e.g., SSMIS, AMSR2) or scatterometers (e.g., ASCAT). The typical resolution of these input images is 12.5 km. The large acquisitions, the repetition of the acquisitions, and their independence with respect to the weather conditions allow a daily polar coverage.

Sequences of SAR images can be used to derive higher resolution drift information. Algorithms have been developed to calculate ice drift from successive pairs of SAR images covering a common area. They are generally based on a spatial correlation calculation between these images, at several resolution scales (from the coarsest to the finest). In Europe, Sentinel-1 is used by DMI to produce an operational sea ice drift product as part of CMEMS (Pedersen et al., 2015). SAR data from Sentinel-1 A/B constellation allow the derivation of daily fields of sea ice deformation at 2 km resolution (Korosov and Rampal, 2017). The algorithm developed by FMI has been operational in the Baltic Sea since early 2011 (Karvonen, 2012), using the wide Radarsat-1 ScanSAR mode and the wide swath ASAR Envisat mode data.

Revisit is the key here: higher revisit of SAR images is naturally required. Algorithms often use image tracking features between consecutive images. This type of algorithms perform better if the images are taken with the same frequency, short interval, and ideally same geometry. Joint acquisition of multi-frequency SAR would enable accurate sea ice drift products, which is not possible with stand-alone current mono-frequency SAR missions. Cross-pol channel is often preferred as it is less sensitive to incidence angle variation (at least for C-band SAR). Drift cannot be calculated in areas without characteristics (i.e., so-called "level ice" and open water areas). These are of course only estimates of the integrated ice trajectories between the time instants corresponding to the SAR acquisition times. Lagrangian drift products are typically 2-day trajectories with coarse resolution (62.5 km) and are not often assimilated even though the accuracy is satisfactory (3 km for 2-days drift).

The TOPAZ4 system does assimilate these sea ice drift operationally, although the effect is relatively weak (Sakov et al., 2012). A different sea ice rheological model, more sensitive to winds, may be more adapted to assimilate this data (Rabatel et al., 2018). More precise (500 m for 1-day drift) and detailed (10 km resolution) sea ice drift products are now obtained year-round from Sentinel-1 SAR images, which cover about 70% of the Arctic as of February 2019.

## AIR-SEA FLUX MEASUREMENTS IN POLAR OCEANS

### The Need and Challenge for Air-Sea Flux Observations in Polar Regions

Air-sea fluxes quantify the exchanges of heat, momentum, freshwater, gases and aerosols between different components of the polar climate system (i.e., atmosphere, water column, sea ice). Flux observations are essential for understanding the global energy budget (Trenberth and Fasullo, 2010), for evaluating forecasting and climate models, and for evaluating processes such as ocean heat uptake and mixed-layer temperature and SST variability (e.g., Dong et al., 2007). However, in situ air-sea or air-ice flux observations are extremely sparse in polar regions (e.g., Bourassa et al., 2013), with almost no winter observations in the Southern Ocean (e.g., Gille et al., 2016; Swart et al., 2019) or in ice-covered ocean domains in both the Arctic and Antarctic (e.g., Taylor et al., 2018). Quantifying air-sea exchange in regions with sea ice requires specialized approaches that must account for spatio-temporal heterogeneity. This has led to significant gaps in our knowledge of both air-sea and air-sea-ice exchanges.

Since there are few reliable near-surface atmospheric observations to serve as constraints, reanalyses and satellite derived surface flux products have major errors and vary considerably between products (Josey et al., 2013; Bentamy et al., 2017; Schmidt et al., 2017; Taylor et al., 2018). For example, in the Southern Ocean, atmospheric models (and reanalyses) have large air-sea heat flux biases, including substantial short-wave errors related to their inability to adequately represent super-cooled liquid cloud water (Bodas-Salcedo et al., 2014, 2016), and these errors appear to bias coupled model SST (Hyder et al., 2018), sea ice, and wind (Bracegirdle et al., 2018). Moored-buoy flux observations or year-round ice camps are required to evaluate and ultimately improve these products but to date buoys have been deployed in only two locations in the Southern Ocean (Schulz et al., 2012; Ogle et al., 2018) and in only a few instances in the Arctic, (Taylor et al., 2018). The Surface Heat Budget for the Arctic Ocean (SHEBA) program offered the only year-round ice camp in the Arctic (e.g., Persson et al., 2002), although Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC; see MOSAiC) will soon extend this. The Southern Ocean Observing System (SOOS) working group on Southern Ocean air-sea fluxes (SOFLUX; Swart et al., 2019) has been working to coordinate observing system capabilities and requirements for high-latitude air-sea fluxes.
## Impact of Flux Uncertainty on
Forecasting/Prediction

There are significant advantages to forecasting with a coupled atmosphere-ice-ocean system, especially on longer timescales. The intrinsic turbulence of the atmosphere limits predictive skill to timescales of days to weeks (Mariotti et al., 2018), whereas the large-scale ocean is predictable on monthly time scales. This provides a mechanism by which the predictability of atmosphere may be extended allowing skillful seasonal predictions.

A particular limitation in this regard is with respect to the significant uncertainty in the atmosphere-ocean boundary layer, made even more egregious when considering sea ice. Boundary layer dynamics involve vertical scales unresolved by coupled forecasting systems, which must be parameterized (e.g., Pullen et al., 2017). The exchanges between components are parameterized by so-called bulk formulae, which estimate air-sea-ice exchanges based on near-surface large-scale properties. Currently, the uncertainties in bulk formulae are a primary bottleneck to seasonal prediction (Penny and Hamill, 2017). In other words, even if we had a perfect ocean model with perfect initial conditions, the information retained in the monthly ocean prediction would be degraded when propagated to the atmosphere due to uncertainty in estimating the true exchanges (Vecchi et al., 2014). Therefore, a priority in the coming decade must be to gather *in situ* estimates of air-sea-ice exchanges in the context of large-scale properties informing how to minimize errors in the bulk formulae parameterizations. Efforts such as YOPP, MOSAiC and ONR Marginal Ice Zone (MIZ) project are examples of projects aiming to address this gap.

We recommend further research to determine how best to represent these exchanges in a coupled forecasting system, with a focus on determining what aspects of boundary layer physics need to be resolved and what can be skillfully parameterized. For parameterized physics, we require that process studies be carried out to determine parameterizations and parameterization coefficients, including identifying the observations that should be sustained to validate estimated fluxes. Addressing these areas are primary goals to enable weekly-to-seasonal skillful predictions in the polar regions.

## FORECASTING SYSTEM EXPERIMENTS

Understanding how assimilating more observations will impact modeled analyses and short-term forecasts is of fundamental importance as new coupled models are being developed. For instance, preliminary results using a regional coupled model [CAFS from the NOAA Earth System Research Laboratory (ESRL)] demonstrated that current ocean reanalyses do not include realistic representation of subsurface water masses relative to observations taken in 2015 and 2016 (including Arctic Heat data). Forecast experiments using a high-resolution fully coupled regional model show that these water masses impact sea ice evolution on synoptic time scales through upper-ocean mixing and heat flux at the ice-ocean interface. It is expected these effects will increase as sea ice continues to decline and surface heat flux processes increase. In addition, assimilating

real-time ocean observations in the initial forecast conditions allows for the identification of biases in the coupled system, which are difficult to isolate when the ocean is allowed to drift away from the observed state. Some studies (e.g., Inoue et al., 2015; Lien et al., 2016; Xie et al., 2017) have begun to show the important potential impact of assimilating enhanced observations on model-based analyses and short-term sea ice forecasts. Additional focused research in this area would allow us to explore coupled model data assimilation issues, better understand physical processes, and assess model performance in comparison to non-coupled (atmospheric) model frameworks. In the next sections, we review a few results from operational centers.

## Impacts of Arctic and Antarctic
Observations in U.S. Navy Coupled
Ice-Ocean Models

The ability to forecast sea ice conditions is of crucial importance for maritime operational planning (Greenert, 2014). The current U.S. Navy operational sea ice forecast system is the Global Ocean Forecast System (GOFS) version 3.1. GOFS utilizes the Los Alamos Community Ice CodE (CICE) version 4.0 (Hunke and Lipscomb, 2008) sea ice model, which is two-way coupled with the HYbrid Coordinate Model (HYCOM) (Metzger et al., 2014) ocean model. The grid resolution is 1/12°, with horizontal resolution approximately 3.5 km at the poles. Atmospheric forcing is provided by the NAVy Global Environment Model (NAVGEM; Hogan et al., 2014). The precursor to GOFS 3.1 was the Arctic Cap Nowcast/Forecast System (ACNFS; Posey et al., 2015). ACNFS is a coupled ice/ocean model similar to GOFS with two main differences: (1) ACNFS only covered areas north of 40°N, and (2) in ACNFS HYCOM has 32 ocean layers compared to 41 in GOFS. An important part of both GOFS and ACNFS is the assimilation of observational data, which is accomplished using the Navy Coupled Ocean Data Assimilation (NCODA) system (Cummings and Smedstad, 2014).

Assimilation of observational data is performed to reduce errors in model forecasts that can result from many factors including non-linear processes that are not deterministic responses to atmospheric forcing, poorly parameterized physical processes, limitations in numerical algorithms, and limitations in model resolution. Polar observational data assimilation is an essential part of GOFS forecasts. Sea ice concentration observations are currently assimilated from the Defense Meteorological Satellite Program (DMSP) SSMIS and AMSR2. These observations are used in conjunction with the Interactive Multi-sensor Snow and Ice Mapping System produced by the U.S. National Ice Center (NIC; Helfrich et al., 2007). A full description of the assimilation process in GOFS and ACNFS can be found in Hebert et al. (2015).

The impact of data assimilation in the U.S. Navy models is significant. One measure that the U.S. Navy uses to determine the accuracy of the modeled sea ice edge (defined in Hebert et al., 2015) is the distance between it and the sea ice edge determined daily by analysts at NIC. The sea ice edge error over six regions in each hemisphere is shown in **Figure 9** as well as the entire Arctic/Antarctic. GOFS was run for 1 year (2014) without data assimilation and compared to the current GOFS
<figure><img src="image_15.png" /><figcaption>**FIGURE 9** | Ice edge error for individual regions (km). Each region contains three numbers. First number is ice edge error without assimilation. Second bold number is error with assimilation. Third number is percent improvement with assimilation. In the Arctic the overall reduction in ice edge error with observational data assimilation is 31 km (56%); in the Antarctic the overall reduction is 28 km (37%).</figcaption></figure>

system with ocean and sea ice data assimilation, including sea ice concentration. The impact of assimilating these observations was to reduce the sea ice edge error in the Arctic by 56% (31 km) and in the Antarctic by 37% (28 km), with a reduction in sea ice edge error for each region ranging between 27 and 68%.

As model resolution increases, so does the need for high-resolution observations. SSMIS and AMSR2 are relatively coarse resolution (25 and 10 km, respectively) compared that of GOFS. Higher resolution (less than 1 km) sea ice concentration observations are available from the Suomi NPP Visible/Infrared (VIIRS). Although VIIRS observations can be obstructed by clouds, including high resolution VIIRS ice concentration into our data assimilation further reduced GOFS ice edge error by 19% (5 km) in the Arctic and 11% (4 km) in the Antarctic. This result points to the need for higher resolution sea ice concentration observations to use in model applications.

In an earlier study using ACNFS to examine the impacts of sea ice concentration observations in ship routing and planning in boreal winter (January–March), assimilating satellite ice concentration observations reduced the projected track an ice breaker would take to a ship near the sea ice edge by an average of 150 km versus not assimilating sea ice concentration observations. This improved the time for planning operations by 12 h and reduced the distance a ship needs to prepare to encounter ice by 212 km.

SIT observations are also important. Currently, pan-Arctic SIT observations on a daily basis are not available. Only limited satellite tracks per day are available that are aggregated on a monthly basis. In a recent study by Allard et al. (2018), ACNFS was re-initialized on March 15th using the March 2014 monthly CryoSat-2 thickness observations and integrated for 18 months.

It showed a reduction in SIT bias by 0.75–0.97 m compared to ACNFS SIT without CryoSat-2 initialization. The impact of this one-time re-initialization was significant and work is underway to assimilate daily satellite track SIT observations on a daily basis.

## Sensitivity of Sea Ice Forecasting Skill to Ocean Mixing Around Antarctica

The rapid evolution of the sea ice cover can have important impacts on coupled environmental predictions through a variety of processes (Smith et al., 2013). These include the formation of leads and coastal polynyas, as well as changes in the ice cover along the marginal ice zone (MIZ). In these regions, the rapid formation, melt and advection of the sea ice cover can modify atmosphere-ocean fluxes on relatively short timescales. Interestingly, small-scale ocean variability has a role to play here as the timing and intensity of changes will be sensitive to the surface ocean mixing layer depth, water mass properties and mesoscale ocean circulation (e.g., Zhang et al., 1999).

As an illustration of the sensitivity of sea ice evolution to ocean mixing, an evaluation of the skill of two sets of sea ice forecasting experiments is shown in **Figure 10**. The first set uses the standard configuration of the Global Ice-Ocean Prediction System (GIOPS: Smith et al., 2016) running operationally at the Canadian Centre for Meteorological and Environmental Prediction. GIOPS combines the Système Assimilation Mercator (SAM2) ocean analysis system with a 3DVar ice analysis (Buehner et al., 2013) to produce daily 10-day forecasts using the NEMO ocean model at 1/4° resolution coupled to the CICE ice model. The second set of experiments is identical to the first with the parameterization for surface wave breaking deactivated. **Figure 10** shows the 7-day forecast skill evaluated
<figure><img src="image_16.png" /></figure>

**FIGURE 10** | Sensitivity of sea ice forecasting skill to ocean mixing around Antarctica. Weekly 7 days sea ice forecasts from the Global Ice-Ocean Prediction System (GIOPS) running operationally at the Canadian Centre for Meteorological and Environmental Prediction are evaluated against analyses over the year 2011. The evaluation of forecast skill is restricted to points where the analysis has changed by more than 10% over the forecast period (7 days). Warmer colors indicate larger root-mean squared error (maximum of 0.3 for dark red) with zero error shown as dark blue. Panels (a) and (b) show the forecast skill for experiments without and with additional ocean mixing respectively. Adapted from Smith et al. (2013).

against 3DVar ice analyses from weekly forecasts over 2011. The verification method used here (Lemieux et al., 2016) restricts the error evaluation to areas where the ice concentration analysis has changed by more than 10% over the forecast lead time (i.e., 7 days). This verification method has the advantage that it focuses the evaluation on 'hot spots' of activity, predominantly in the marginal ice zone.

From **Figure 10** it can be seen that a small modification to the ocean vertical mixing can have a first order impact on the ice forecast errors. Interestingly, while the surface wave breaking parameterization degrades ice forecast skill, it does lead to an improvement in water mass properties over ice-free waters (as evaluated against Argo profiles; not shown). This is perhaps not surprising given that the mixing regime in polar regions is quite different from that at lower latitudes. This highlights the need for an expanded under-ice ocean monitoring program to be able to adequately model vertical mixing and constrain water mass properties and mixed layer depths.

## Impact of Temperature and Salinity Profiles in the CMEMS Arctic MFC

There have been special observing periods of the Arctic and Antarctic in the past, in particular the successive International Polar Years (IPYs), the latest of which took place in 2007–2009 (with a gradual ramping up of ocean observing systems in the preceding years). Looking back at the impact of a recent IPY in a period with similar low-ice-coverage conditions in the Arctic, expanded sea ice in the Antarctic, and similar satellite coverage as today can provide another measure of the expected impact of the YOPP Special Observing Periods.

In the Arctic, the TOPAZ4 reanalysis is based on a regional configuration of the HYCOM model coupled to an early version of the CICE model at 12 km resolution. It assimilates both satellite and *in situ* observations using an Ensemble Kalman Filter (Sakov et al., 2012). As of 2018, the TOPAZ4 reanalysis system is almost

identical to the real-time physical forecasting system used in the Copernicus Arctic MFC and operated by MET Norway (Melsom et al., 2015; Xie et al., 2017). It assimilates the same types of ocean and sea ice observations: along-track sea-level anomalies from altimeters, sea surface temperatures, sea ice concentrations and drift and *in situ* temperature and salinity profiles.

Ocean models have well-known limitations in simulating the advection of Atlantic Waters into the Arctic (Ilicak et al., 2016), which for a model like TOPAZ4 results in a typical cold bias of the 300–800 m water temperature by 0.5 K (see **Figure 11**, blue line) and regional differences by 1.5 K typically (green line). During the IPY time period, the assimilation of ITPs is able to constrain this bias down to 0.1 K and reduce the regional differences below 1 K. The end of the IPY in 2009 relaxes the constraint and the bias and RMS errors return to larger values, although slightly lower than before the IPY (i.e., there is no “dynamical shock” after the IPY stops). Even though the quantitative impact on the TOPAZ4 system is dependent on the practical setup of the model and its assimilation scheme, the qualitative behavior may apply to other forecasting systems based on similar types of models and data assimilation schemes and indicates that a density of ITP profiles at least equal to that of the IPY should be sustained continuously to constrain efficiently the Atlantic Water properties in the Arctic.

## NOAA-ESRL/CIRES Coupled Arctic Forecast System

NOAA ESRL has provided experimental, daily, 10-day forecasts of Arctic weather and sea ice evolution to stakeholders during freeze-up seasons since 2015 and daily forecasts year-round starting on February 14, 2018. CAFS produces high-resolution (10 km) regional coupled-model Arctic forecasts using global 0.5° GFS forecast fields for lateral forcing. The current configuration of the model includes the POP2 dynamical ocean model (adapted from Maslowski et al., 2012), the CICE5 sea ice model (Hunke et al., 2015), the NCAR CLM4.5 land model, and the WRF3.6 ARW atmospheric model, coupled by the NCAR CPL7 coupler. The domain is the Arctic basin and surrounding coastal regions, including Bering Strait, to provide model guidance for the National Weather Service (NWS), and Fram Strait, to include the complete planned Multidisciplinary drifting Observatory for the Study of Arctic Climate campaign (MOSAiC) domain.

The CAFS forecasts are being used by NOAA ESRL to identify sources of skill on sub-seasonal time-scales due to coupled ocean-ice-atmosphere processes and by stakeholders as model guidance for sea ice forecasts. Real-time CAFS products are made available to the community<fnref n="12" />. Figures and animations from the 10-day forecasts are provided for sea ice, atmosphere, and ocean variability, as well as, an archive of model output for users to download. These forecasts are being used for model guidance by the NWS Alaska Sea Ice Program, the NOAA Arctic Testbed, the U.S. National Ice Center, and by the U.S. Navy and NOAA for operations during Arctic campaigns.

Coupled Arctic Forecast System forecasts use initial conditions that ingest SIT measurements from ESA's CryoSat-2 and SMOS satellites and the NASA Jet Propulsion Laboratory Multi-Scale Ultra-high Resolution Sea Surface Temperature

<footnote n="12">¹²https://www.esrl.noaa.gov/psd/forecasts/seaice</footnote>
<figure><img src="image_17.png" /><figcaption>**FIGURE 11** | Time series of TOPAZ4 data assimilation diagnostics across the 24-year reanalysis for all temperature profiles in the depths 300–800 m in the whole Arctic. The blue line is the bias, the green line is the related standard deviation (Root Mean Square), the red line is the ensemble spread, and the gray line the number of temperature observations, increasing during the IPY. The other vertical lines and the bottom bars indicate changes of the other observation data sources and modifications of the data assimilation system.</figcaption></figure>

(MUR) SSTs and sea ice concentration. The MUR SSTs are used to initialize the ocean mixed layer with a mixed-layer depth diagnosed from the model. In order to identify whether using these satellite products in the initial conditions increases the skill of the 10-day forecasts, a series of 10-day hindcasts were performed for the time period of the ONR SeaState DRI, Oct. 1-Nov. 5 2015. The hindcasts are setup exactly like the real-time forecasts except the lateral boundary conditions are the GFS analyses instead of the GFS forecasts, in order to identify potential model biases. In addition, similar to the forecasts, hydrometeor mass and number are initialized with fields from the first day of the previous day's hindcast to reduce spin-up time of cloud fields.

Intensive measurements were taken of the ocean, surface, and atmospheric state during the SeaState campaign. This provides for a comprehensive observational database for model validation. **Figure 12** shows the model error of the ocean-atmosphere state at the location of the R/V Sikuliaq over the 6-week campaign at lead times of 6 h, 1 day, and 5 days. At a 6-h lead time, there are equivalent errors in the ocean and atmosphere, less than 0.5°C. Due to the initialization of the ocean mixed-layer with satellite SSTs, errors in the ocean grow slowly up to 5-day lead times. However, differently in the atmosphere, errors in the lowest 2 km grow rapidly with errors greater than +1°C in the near-surface temperatures and errors greater than -1°C at the top of the atmospheric boundary layer. This is an indication that the model

is unable to maintain the observed boundary layer stratification and rapidly evolves into a less-stable state. It would not have been possible to identify this model bias without initialization of the ocean/ice state. Process studies are currently underway to identify if this bias is due to cloud processes or boundary layer parameterizations.

## Quantitative Network Design

The Quantitative Network Design (QND) approach constitutes a computationally efficient alternative to Observing System Experiments (Fujii et al., 2019). The approach can be used to inform the design of observing networks or space missions. QND evaluates a set of observations (network) in terms of its constraint on a target quantity, i.e., a quantity of interest. This evaluation is performed in a modeling system that is capable of simulating counterparts of the observations and of the target quantity from a set of unknowns in the system (control variables). For a detailed description of the formalism we refer to Kaminski and Rayner (2017). Briefly, it proceeds in two steps: In the first step, the observational information is used to infer the uncertainty in the posterior control vector, $C(x)$, that is consistent with the observational uncertainty, $C(d)$, and the prior uncertainty of the control vector, $C(x_0)$ [Equation (1)]. In the second step the uncertainty in the control vector is mapped onto an uncertainty in the target quantity, $\sigma(y)$ [Equation (2)]. Both steps use appropriate sensitivities (linearizations/response
<figure><img src="image_18.png" /><figcaption>FIGURE 12 | Box and whisker plots of CAFS forecast temperature errors in the full atmosphere-ocean column at 6 h, 1 day, and 5 days lead times compared to radiosondes (red) and CTDs (blue) from the R/V Sikuliaq during the ONR SeaState campaign October 1 – November 5 2015. Note, the vertical scale is model levels.</figcaption></figure>

functions) of the modeling system. The first step, an inversion step, is formalized by

$$C(x)^{-1} = M'TC(d)^{-1}M' + C(x_0)^{-1}, \quad (1)$$

and the second step by

$$\sigma(y)^2 = N'C(x)N'T + \sigma(y_{\text{mod}})^2, \quad (2)$$

where $M'(N')$ denotes the sensitivity of the vector of observations (the target quantity) with respect to the control vector, as simulated by the modeling system, and $\sigma(y_{\text{mod}})$ is the residual uncertainty in the simulation of the target quantity that remains even for a perfect control vector. The approach represents a network through observational locations and times and the observational uncertainty, but does not require real observations. Consequently, it can evaluate hypothetical networks/space missions based on assumed instrumental specifications and space-time coverage. Observations can range from point-scale (in situ) to gridded data sets or level-1 satellite data, if appropriate forward models/observation operators mapping the model's state variables onto the respective data stream (Kaminski and Mathieu, 2017) are available.

Historically, QND was first applied by Hardt and Scherbaum (1994) to optimize the locations of a set of seismic sensors. Rayner et al. (1996) applied the approach to the in situ sampling network for atmospheric CO₂. For the physical sea ice-ocean system in

the Arctic domain, Kaminski et al. (2015) applied the approach to evaluate idealized flight transects derived from NASA's Operation IceBridge airborne altimeter ice surveys. Target quantities were 10-day to five-month forecasts of snow and ice volumes over areas relevant for maritime traffic (along the Alaskan coast) and offshore resource exploration (Chukchi Sea). The control vector was composed of physical constants in the model's process representations as well as initial and boundary conditions.

In an activity funded through ESA's Support to Science Element as a contribution to YOPP, Kaminski et al. (2017) constructed the Arctic Mission Benefit Analysis (ArcMBA) system. The system exploits the fact that model sensitivities at observational times and locations as well as the target quantities can be pre-computed, so that the actual assessment of a data set requires only matrix multiplications and inversions [Equation (1) and Equation (2)]. This means the assessment can be performed so fast that the ArcMBA system could be used as an interactive tool to assist decision makers, for example, in a meeting. Currently the system uses pre-computed model sensitivities $M'$ for observations of SIT, sea ice freeboard, radar freeboard, laser freeboard, and snow depth along with precomputed model sensitivities $N'$ for forecasts of snow and ice volumes for three regions along the Northern Sea Route (North-East Passage) respectively denoted by West Laptev Sea, Outer New Siberian Islands, and East Siberian Sea. These sensitivities were computed with the Max-Planck-Institute Ocean Model
(MPIOM, Jungclaus et al., 2013) in a global configuration with high resolution over the Arctic (Niederdrenk, 2013). The control vector is composed of physical constants in the model's process representations as well as initial- and boundary conditions.

Kaminski et al. (2017) applied the system to evaluate real and hypothetical remote sensing products. The real products were monthly SIT, sea ice freeboard, and radar freeboard, all derived from CryoSat-2 by AWI. These real products were complemented by two hypothetical monthly laser freeboard products (2 and 20 cm accuracy, respectively), as well as two hypothetical monthly snow depth products (2 and 15 cm accuracy, respectively). Target quantities are 4-week forecasts of snow and ice volumes over three target regions along the Northern Sea Route.

As an example, **Figure 13** shows the posterior uncertainty in snow and ice volumes over the three target regions, when the CryoSat-2 SIT product (Ricker et al., 2014) is used alone, in combination with the hypothetical snow depth product with 15 cm accuracy, or in combination with the hypothetical snow depth product with 2 cm accuracy. To sharpen the contrast between the observational scenarios, $\sigma(y_{\text{mod}})$ (which acts as an offset) is set to zero (perfect model scenario). Comparison of the top and middle panels shows the added value of the hypothetical snow product, not only for the snow volume forecast but also for the sea ice volume (SIV) forecast. SIV is sensitive not only to initial SIT but also to initial ice concentration and snow depth, which are both constrained by the snow depth product. SIV is also sensitive to some of the process parameters that are constrained by snow depth, notably the ice strength, see Kaminski et al. (2017) for details. Comparison of the middle and top panels shows the added value of a higher accuracy in the snow depth product. Increasing the accuracy from 15 to 2 cm results in a reduction in uncertainty of the SIV forecast for the East Siberian Sea target region from to 63 to 24 km³.

The ArcMBA can be extended to cover further Earth Observation (EO) products and further target variables. In the setup used here, the model can simulate a range of sea ice-ocean variables in addition to those considered in the present study (e.g., ice drift, mixed layer depth, freshwater/sea surface salinity, SST, circulation). Switching to a more comprehensive model configuration would enable the investigation of yet further variables. For example the model can be operated with its biogeochemistry module HAMOCC activated (Ilyina et al., 2013) or in a mode coupled to an atmospheric general circulation model and thus enable the analysis of biogeochemical products/target quantities or Arctic mid-latitude linkages. The extension of ArcMBA by a terrestrial biosphere component is planned, which will allow joint assessment of ocean and land observations.

## INTERNATIONAL EFFORTS TO ADDRESS GAPS IN POLAR OCEAN OBSERVATIONS

Various international efforts contribute to coordinate and support the vast and complex polar observing networks. These networks are maintained by a collection of national and international efforts and scientific projects. These include

national efforts outlined above, such as IAOOS, FRAMs and NorArgo as well as diverse multi-platform observing projects such as the US SODA (Stratified Ocean Dynamics of the Arctic) program (Lee et al., 2016) that include many of the types of platforms described in Section “In situ Observations of Temperature and Salinity” (WHOI ITPs, ALAMO floats, drifters, gliders). In the south, these efforts are coordinated through the Southern Ocean Observing

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /></figure>

**FIGURE 13** | Evaluation of CryoSat-2 sea ice thickness product alone (top), in combination with 15 cm accuracy snow depth product (middle), and in combination with 2 cm accuracy snow depth product (bottom). Prior (gray, no observations) and posterior (orange, sea ice thickness and snow depth products) uncertainties in sea ice volume (SIV) and snow volume (SNV) predictions for three regions along the Northeast Passage in km³.
System (SOOS; Newman et al., 2019) supported by the Scientific Committee on Ocean Research (SCOR) and the Scientific Committee for Antarctic Research (SCAR).

A clear need for a coordinated approach for the Arctic has arisen following the International Polar Year (IPY) and is being developed through a series of bi-annual Arctic Observing Summits (Murray et al., 2018). These summits contribute to a broad initiative by Sustaining Arctic Observing Networks (SAON) that has arisen from the International Arctic Science Committee. Recently, a framework for the development of an Arctic Region Global Ocean Observing System (ARGOOS) has been proposed by Lee et al. (2019). Additionally, the CLIVAR/CliC Northern Oceans Panel serves as an international forum for coordinating and strategizing activities on the role of the Arctic Ocean in the context of the global climate system from a coupled perspective and facilitates progress in developing new tools and methods to observe the Arctic Ocean and neighboring seas.

While these efforts aim to address a broad range of societal needs, in this section we describe several particular initiatives that work to address specific gaps in the observing system important for environmental prediction.

## INTAROS (INTegrated ARctic Observation System)

INTAROS is a research and innovation project under the European Union Horizon 2020 program, running from 2016 to 2021, with objective to build an efficient integrated Arctic Observation System in the Arctic. This requires collaborative efforts among many institutions to extend, improve and unify existing systems, which in many cases are designed and developed for specific scientific disciplines. INTAROS focuses on the *in situ* part of the observing systems, which represent the largest component of the integrated observing system. Satellite Earth Observation programs provide the most developed and operational components of the system, which are run by space agencies and satellite monitoring services such as Copernicus<fnref n="13" />. The satellite systems provide data for near real-time monitoring as well as for long-term climate observations. Validation of the satellite-derived variables is an important part of the operational services. *In situ* observations play an important role for this validation, but there is very limited access to such data from the Arctic Ocean. Furthermore, most of the ocean observations are only available in delayed mode, because they are provided by underwater moorings and seafloor observatories. Some platforms that operate at the surface can transmit data in near real-time and can therefore contribute to operational monitoring of sea ice and ocean variables.

INTAROS is multidisciplinary, implying that the observing systems encompass atmospheric, marine and terrestrial systems in the different regions of the Arctic. Marine observing systems are also divided into physical and biogeochemical components of the ocean surface (including sea ice), the water column and the seafloor. INTAROS contributes to all these components in collaboration with other observing programs and projects, by

deploying new sensors and platforms to enhance the observing capacity in different Arctic regions. The collection of new *in situ* data under INTAROS started in 2017 and will continue through 2020, coinciding with and contribute to the extensive data collection in the MOSAiC program (see MOSAiC). Examples of INTAROS supported platforms include: IMBs deployed in the Central Arctic Ocean by the Finnish Meteorological Institute during Chinese Arctic expeditions (Lei et al., 2018); a network of Argo floats in the Nordic Seas in ice-free areas; bio-Argo floats in Baffin Bay together with Takuvik; and the Ferry-box transect between Norway and Svalbard<fnref n="14" />. INTAROS also supports data collection from the cargo vessel NORBJØRN that runs between Tromsø and Longyearbyen throughout the year. Data can be transmitted to land in near realtime and is available for operational monitoring. Ship of Opportunity is promising method to collect oceanographic data in Arctic waters, since there is a growing number of ships operating in the Arctic, especially tourist ships.

Arctic observing systems will benefit greatly from collaboration with local communities (Johnson et al., 2018). INTAROS is therefore working with Community-Based Monitoring (CBM) systems which are under development in many places in North America, Greenland and Russia. CBMs can be supplementary to scientific observations when indigenous and local people collect scientifically relevant data and made them available via websites<fnref n="15" />. Community-based monitoring can also provide valuable data that cannot be obtained from normal research and monitoring programs. In the circumpolar Arctic region there are a number of observing programs addressing sea ice, oceanographic data and observations of marine mammals and fish which are very important for the communities.

## EMODNET

Observations of the ocean are usually made for specific purposes. In order to save costs and improve marine knowledge, the European Union is now moving to a new paradigm where data are collected once and then used for many purposes. This paradigm is being implemented as part of the European Marine Observation and Data Network (EMODnet) consisting of a variety of organizations working together to assemble marine data and products, and to facilitate the dissemination of these resources to both public and private users. EMODnet is currently in its third development phase with the target to be fully deployed by 2020. Six Sea Basin Checkpoint programs were initiated as part of EMODnet, including the Sea Basin Checkpoint Arctic<fnref n="16" />, with the objective to examine the current status of Arctic data and to assess how fit-for-purpose the Arctic data are.

The project (2015–2018) aimed to identify problems and knowledge gaps and was organized in 10 challenges: wind farm sitting, marine protected areas, oil platform leak, climate change, coasts, fisheries management, fisheries impact, river input, bathymetry and alien species. Within the project, each dataset found (292 in total) was assessed on its spatial and temporal

<footnote n="13">¹³https://www.copernicus.eu/en</footnote>

<footnote n="14">¹⁴https://www.niva.no/en/projectweb/norsoop</footnote>

<footnote n="15">¹⁵https://mkp28.wixsite.com/cbm-best-practice</footnote>

<footnote n="16">¹⁶http://www.emodnet-arctic.eu</footnote>
coverage, its accessibility and cost to access, the responsiveness and the temporal and vertical resolution. Also, each assessment report using a dataset (840 in total) was assessed on its adequacy for the project. For the oil spill challenge an oil accident was simulated. This demonstrated the necessity for rapid acquisition and inspection of ocean current and wind data in order to provide a reliable response capacity.

## MOSAiC

An important Arctic observing project called MOSAiC (Multidisciplinary Drifting Observatory for the Study of Arctic Climate; Shupe et al., 2018) is planned for 2019–2020. The German research icebreaker Polarstern will be frozen into the pack ice and over-winter in the Transpolar Drift to obtain measurements over a complete annual cycle. The MOSAiC sea ice platform provides an opportunity for greatly enhanced deployment of autonomous instrumentation and coordinated intensive field studies from research vessels, manned and unmanned aircraft, and distributed surface stations.

MOSAiC observations have been designed specifically to characterize the important coupled processes within the atmosphere-ice-ocean system that impact sea ice melting and freezing. These processes include heat, moisture, and momentum fluxes in the atmosphere and ocean, water vapor, clouds and aerosols, biogeochemical cycles in the ocean and ice, and many others.

The MOSAiC central observatory will be a manned, icebreaker ship-based ice camp with comprehensive instrumentation to measure coupled processes within the atmosphere, ice, and ocean. This intensive observatory will be embedded within a constellation of distributed measurements made by buoys, ice-tethered profilers, remote meteorology stations, underwater drifters, unmanned aerial systems, aircraft, additional ships, and satellites. These distributed observations will provide critical information on the spatial context and variability of key parameters, and allow for limited measurements in environments with sea ice of differing age, thickness, and concentration.

The additional observations from MOSAiC together with other sources during YOPP provide an ideal opportunity to assess the impact on forecast skill in environmental prediction systems through OSEs and process studies. Results will provide a quantitative baseline for use in decisions regarding how to configure a sustained Arctic observing system appropriate for the needs of environmental prediction.

## Year of Polar Prediction

The need for improved environmental predictions (i.e., atmosphere, ice and ocean) has motivated an international effort called the Polar Prediction Project (PPP<fnref n="17" />). The PPP was created under the auspices of the World Meteorological Organisation's World Weather Research Programme (WWRP) to “promote cooperative international research enabling development of improved weather and environmental prediction services for the polar regions, on time scales from hourly to seasonal.” YOPP is a PPP flagship activity (Jung et al., 2016) with

a core phase during 2017-2019, with the overarching goal to “Enable a significant improvement in environmental prediction capabilities for the polar regions and beyond, by coordinating a period of intensive observing, modeling, verification, user-engagement and education activities.” In particular, there are a series of Special Observing Periods for both the Arctic and Antarctic to improve the polar observing system to provide better coverage of high-quality observations in a cost effective manner, primarily by carrying out Observing System Experiments. The Southern Hemisphere SOP occurred from November 2018 to February 2019, and included deployment of extra radiosondes and drifting buoys, as well as the prolongation of the Ocean Observatories Initiative surface mooring at 55°S. The Arctic SOPs were conducted for winter (February–March, 2018) and summer (July–September, 2018). A more comprehensive list of research activities is given by PPP Steering Group (2013, 2014). Additional observations gathered through field programs will also be used to improve our understanding of key polar processes relevant for improving prediction skill. A third Arctic SOP is planned for February–March 2020 as part of YOPP to capitalize on the MOSAiC observational effort.

The additional observations from MOSAiC together with other sources during YOPP provide an ideal opportunity to assess the impact on forecast skill in environmental prediction systems through OSEs and process studies. Results will provide a quantitative baseline for use in decisions regarding how to configure a sustained Arctic observing system appropriate for the needs of Environmental Prediction.

## RECOMMENDATIONS

The scarcity of observations, the unique balance of physical processes, the key importance of sea ice, and the rapidly evolving climate of the Arctic, and the uncertainties in Antarctic sea ice trends and carbon uptake lead to a number of scientific challenges for observations in the context of a polar prediction system. Addressing these challenges motivates the following recommendations:

1. The presence of a seasonal ice cover limits the availability of real-time *in situ* data in polar regions to assist operational requirements. While several technologies have been developed (e.g., ITPs, gliders communicating via acoustic modems) a comprehensive real-time ocean observing network able to supplement Argo for polar regions has yet to be put in place, hindering progress toward coupled polar prediction. It is therefore recommended that a network of ice-borne measurement systems be deployed and supported operationally in ice-covered areas. These platforms represent well-proven technologies for year-round data collection and near-real time data transmission via satellites.

2. Antarctic measurements are also needed, in particular, to evaluate changes that could be harbingers of continental ice melt. Recent studies also highlight the importance of measuring and understanding the intra-hemispheric ocean interactions on numerical weather prediction to climate time-scales (Foppert et al., 2017). Ice-borne observing

<footnote n="17">¹⁷http://www.polarprediction.net</footnote>
systems that have proven their utility in the Arctic should now be adapted and tested in the Antarctic marginal ice zone.

3. Conditions are changing rapidly with the loss of summer sea ice extent in the Arctic and changing ice-cover patterns in Antarctic marginal seas. Phenomena long considered negligible in the Arctic and Antarctic may be becoming important (e.g., ocean waves - Cavaleri et al., 2012; Massom et al., 2018). The observing system needs to be reevaluated with this in mind. The retreat of the Arctic ice cover increases the area where open-water (or seasonally ice free) observing systems (e.g., gliders, Argo) can operate providing the possibility in future to provide a substantial amount of near-real time data for polar prediction systems.

4. Ships of opportunity present a promising method to collect oceanographic data in polar waters, since there is a growing number of ships operating in these regions. Moreover, community-based monitoring can also provide valuable data that cannot be obtained from normal research and monitoring programs. These atypical observing methods should be encouraged and enhanced to provide a low-cost expansion to the *in situ* observing system.

5. *In situ* observations are routinely used for the calibration of remote sensing products over much of the globe, with fewer such calibrations made in polar regions due to lack of *in situ* observations. Studies have shown the large benefits that such calibrations can have (e.g., Castro et al., 2016; Kwok et al., 2017). The availability of near-real time *in situ* measurements could be used to improve the quality of satellite products and thus on downstream environment predictions that assimilate them. Additional efforts involving multi-platform calibration are needed to improve the quality of remote sensing products.

6. The increasing maturity of satellite sea-ice thickness winter-time products merging several sensors (e.g., CryoSat-2 and SMOS) and its positive impact in preliminary assimilation experiments call for symmetrical efforts in the Antarctic ocean, where such products do not exist at the moment.

7. Polar surface properties are often dominated by various forms of ice that vary rapidly on small spatial scales. Some remote sensing methods of ice properties (ice cover, ice thickness, snow depth on ice, albedo, crystal structure) are not mature and offer little and/or coarse-resolution information from within the ice, whereas *in situ* methods provide poor spatial coverage. Neither is currently able to address the need for high spatial and temporal resolution observations of sea ice deformation over large regions. Observations providing information regarding ice deformation and redistribution during ridging are also lacking. There is a need for high-resolution (km-scale) remotely sensed snow and ice property data for both the Arctic and Southern Ocean with sufficient temporal resolution to address these relevant features.

8. The value of polar observations for multi-range environmental prediction emerged during the last decade from a variety of impact studies. The importance of SIT

initialization for seasonal forecasting, the significance of sub-surface initialization in ice-covered areas, are non-exhaustive examples that call for coordinated efforts (including QND, OSEs and Observing System Simulation Studies) to enhance the Arctic and Southern Ocean observing networks.

9. To date, few data withholding experiments (OSEs) or observation design experiments (OSSEs) have been undertaken for polar regions using real-time prediction systems. Performing such experiments using additional observations available during YOPP (or other periods of additional observational coverage such as IPY) is suggested to provide valuable information to guide the design of a sustainable real-time observing system for polar regions suitable for environmental prediction. In particular, multi-system exercises shall be encouraged to gain robustness in the observation impact assessments.

10. Polar environmental prediction using coupled atmosphere-ice-ocean models is strongly sensitive to errors in fluxes across the surface interface and thus requires collocated information about the state of the atmosphere, sea ice and ocean, to be used for improving interface fluxes (i.e., coupled model validation) and eventually data assimilation. Direct flux covariance measurements, in particular, would be immensely valuable in constraining bulk parameterizations used to represent fluxes in models.

11. Open access to data, especially real-time data, is a critical capability for improved sea-ice and weather forecasting and other environmental prediction needs. The optimal observing system will no doubt include a suite of different instrument types, since no single platform can be optimized for the full range of observing needs (i.e., sub-surface and lower atmosphere), and hence multi-platform coordination will be necessary, including coordination with local community-based observatories. This means that real-time dissemination of *in situ* observations in polar regions to global data assembly centers must be prioritized in order to make the observational efforts suitable for environmental prediction applications.

12. International collaboration will continue to be key for facilitating deployment of polar ocean instrument systems, including the fielding of drifting and anchored buoys, floats and gliders and free, rapid dissemination of the resulting data.

## OUTLOOK

The relative remoteness and harsh environmental conditions over polar regions will always hinder efforts to provide adequate observations for polar prediction. Over recent years, we have seen improvements in observing technology and capabilities that create new possibilities for how to construct and maintain the polar ocean observing system. The technologies make an adequate polar ocean observing system feasible, but the question remains, is it worth the cost? YOPP aims
to help address this question by coordinating international observing activities and fostering efforts to assess the impact of additional observations on environmental prediction skill, including impacts on downstream users and products. Following the YOPP core period (2017–2019), there will be a consolidation phase to assess these impacts and develop recommendations toward sustained polar observation. This effort will culminate in a YOPP Final Summit planned for summer 2022, providing an ideal opportunity for funding and implementation agencies to benefit from this community effort.

## AUTHOR CONTRIBUTIONS

GS developed the structure of the manuscript and coordinated the writing team. GS and MC contributed to section “Introduction.” SSa contributed to section “Overview of Current Observing System.” JT contributed to section “The WHOI Ice-Tethered Profiler.” KW, JI, and CM contributed to section “ALAMO and Arctic Heat.” GS and FD contributed to section “Deployment of Argo Floats in Canadian Marginal Ice Zone by CONCEPTS.” CM, JL, ER, and MB contributed to section “Under-Ice BGC Argo Floats in the Canadian Arctic.” IR contributed to section “Surface Drifters and IABP.” GC and AS contributed to section “Sea Ice Concentration.” QY, PPr, PT, and GL contributed to section “Sea Ice Freeboard and Thickness.” PT, PPr, and GL contributed to section “Sea Surface Temperature.” PT, PPr, and GL contributed to section “Sea Surface Height.” PT, PPr, GL, LB, and JX contributed to section “Sea-Ice Drift.” SSw, SG, PH, MM, and BD contributed to section “Air-Sea Flux Measurements in Polar Oceans.” RA, JC, DH, EM, MP, PPo, and OS contributed to section “Impacts of Arctic and Antarctic Observations in U.S. Navy Coupled Ice-Ocean Models.” GS contributed to section “Sensitivity of Sea Ice Forecasting Skill to Ocean Mixing Around Antarctica.” LB and JX contributed to section “Impact of Temperature and Salinity Profiles in the CMEMS Arctic MFC.” ASo contributed to section “NOAA-ESRL/CIRES Coupled Arctic Forecast System.” TK and FK contributed to section “Quantitative Network Design.” GC contributed to section “International Efforts to Address Gaps in Polar Ocean Observations.” SSa contributed to section “INTAROS (INTegrated ARctic Observation System).” BK contributed to section “EMODNET.” MS and JI contributed to section “MOSAiC.” GS, MC, SO, PR, and WWRP PPP Steering Group contributed to section “YOPP.” All coauthors contributed to section Recommendations. Revised manuscript were made by GS with particular contributions to overall editing by SG, SSw, ASt, JI, MC, and LB.

## REFERENCES

Alexandrov, V., Sandven, S., Wahlin, J., and Johannessen, O. M. (2010). The relation between sea ice thickness and freeboard in the Arctic. *Cryosphere* 4, 373–380. doi: 10.5194/tc-4-373-2010

Allard, R. A., Farrell, S. L., Hebert, D. A., Johnston, W. F., Li, L., Kurtz, N. T., et al. (2018). Utilizing CryoSat-2 sea ice thickness to initialize a coupled ice-ocean modeling system. *Adv. Space Res.* 62, 1265–1280. doi: 10.1016/j.asr.2017.12.030

## FUNDING

The development of the new generation of floats (PRO-ICE) to be operated under ice was funded by the French project NAOS. Twelve PRO-ICE were funded by NAOS and nine by the Canadian Foundation for Innovation (FCI-30124). The GreenEdge project is funded by the following French and Canadian programs and agencies: ANR (Contract #111112), CNES (project #131425), IPEV (project #1164), CSA, Fondation Total, ArcticNet, LEFE and the French Arctic Initiative (GreenEdge project). The INTAROS project has received funding from the European Union's Horizon 2020 Research and Innovation Program under grant agreement No. 727890. The setup of the ArcMBA system and the experiment described in section “Quantitative Network Design” were funded by the European Space Agency through its support to science element (contract #4000117710/16/I-NB). SSw was supported by a Wallenberg Academy Fellowship (WAF 2015.0186). The work at CLS (GL, PPr, and PT) has been funded by internal investment, in relation with on-going CNES and ESA funded studies making use of radar data over Polar regions. EMODNET (BK) is funded by the European Commission. NRL Funding (for RA, JC, DH, EM, PPo, OS) provided by NRL Research Option “Determining the Impact of Sea Ice Thickness on the Arctic’s Naturally Changing Environment (DISTANCE), ONR 6.2 Data Assimilation and under program element 0602435N (JC, RA, DH). JT’s Arctic research activities are supported by the U.S. National Science Foundation and ONR. SG was funded by NSF grants/awards PLR-1425989 and OCE 1658001. IR is funded by contributors to the US IABP (including CG, DOE, NASA, NIC, NOAA, NSF, ONR). CAFS is supported by the NOAA ESRL Physical Sciences Division (AS and JI). LB and JX are funded by CMEMS. The WWRP PPP Steering Group is funded by a WMO trust fund with support from AWI for the ICO. The publication fee is provided by ECCC.

## ACKNOWLEDGMENTS

We would like to acknowledge the valuable suggestions made by both reviewers that has improved the quality of this manuscript. Contributions to the PPP Trust Fund are acknowledged as well as AWI’s financial contribution to hosting the ICO. For work done for section “Under-ice BGC Argo Floats in the Canadian Arctic,” we thank E. Leymarie and C. Penkerc’h (LOV Laboratoire d’Océanographie de Villefranche sur Mer) for their contribution in the development of the PRO-ICE. We also thank officers and crew of CCGS Amundsen for their help during our deployments of floats.

Andersen, O. B., and Piccioni, G. (2016). Recent arctic sea level variations from satellites. *Front. Mar. Sci.* 3:76. doi: 10.3389/fmars.2016.00076

Ardhuin, F., Aksenov, Y., Benetazzo, A., Bertino, L., Brandt, P., Caubet, E., et al. (2018). Measuring currents, ice drift, and waves from space: the Sea surface Kinematics Multiscale monitoring (SKIM) concept. *Ocean Sci.* 14, 337–354. doi: 10.5194/os-14-337-2018

Armitage, T. W. K., Bacon, S., Ridout, A. L., Thomas, S. F., Aksenov, A., and Wingham, D. J. (2016). Arctic sea surface height variability and change from
satellite radar altimetry and GRACE, 2003-2013. *J. Geophys. Res. Oceans* 121, 4303–4322. doi: 10.1002/2015JC011579

Armitage, T. W. K., and Ridout, A. L. (2015). Arctic sea ice freeboard from AltiKa and comparison with Cryosat-2 and Operation IceBridge. *Geophys. Res. Lett.* 42, 6724–6731. doi: 10.1002/2015GL064823

Bell, M. J., Schiller, A., Le Traon, P.-Y., Smith, N. R., Dombrowsky, E., and Wilmer-Becker, K. (2015). An introduction to GODAE OceanView. *J. Operat. Oceanogr.* 8, s2-s11. doi: 10.1080/1755876X.2015.1022041

Bentamy, A., Piolle, J.-F., Grouazel, A., Danielson, R., Gulev, S., Paul, F., et al. (2017). Review and assessment of latent and sensible heat flux accuracy over the global oceans. *Remote Sens. Environ.* 201, 196–218. doi: 10.1016/j.rse.2017.08.016

Berge, J., Geoffroy, M., Johnsen, G., Cottier, F., Bluhm, B., and Vogedes, D. (2016). Ice-tethered observational platforms in the Arctic Ocean pack ice. *IFAC Pap. Online* 49, 494–499. doi: 10.1016/j.ifacol.2016.10.484

Blockley, E. W., and Peterson, K. A. (2018). Improving Met Office seasonal predictions of Arctic sea ice using assimilation of CryoSat-2 thickness. *Cryosphere* 12, 3419–3438. doi: 10.5194/tc-12-3419-2018

Bodas-Salcedo, A., Hill, P. G., Furtado, K., Williams, K. D., Field, P. R., Manners, J. C., et al. (2016). Large contribution of supercooled liquid clouds to the solar radiation budget of the Southern Ocean. *J. Clim.* 29, 4213–4228. doi: 10.1175/jcli-d-15-0564.1

Bodas-Salcedo, A., Williams, K. D., Ringer, M. A., Beau, I., Cole, J. N., Dufresne, J. L., et al. (2014). Origins of the solar radiation biases over the Southern Ocean in CFMIP2 models. *J. Clim.* 27, 41–56. doi: 10.1175/jcli-d-13-00169.1

Bourassa, M. A., Gille, S. T., Bitz, C., Carlson, D., Cerovecki, I., Clayson, C. A., et al. (2013). High-latitude ocean and sea ice surface fluxes: challenges for climate research. *Bull. Am. Meteorol. Soc.* 94, 403–423. doi: 10.1175/bams-d-11-00 244.1

Bracegirdle, T. J., Hyder, P., and Holmes, C. R. (2018). CMIP5 uncertainty in projected 21C change in southern westerlies significantly related to simulated historical sea ice area; strong link to strengthening and weak link to shift. *J. Clim.* 31, 195–211. doi: 10.1175/jcli-d-17-0320.1

Brassington, G. B., Martin, M. J., Tolman, H. L., Akella, S., Balmeseda, M., Chambers, C. R. S., et al. (2015). Progress and challenges in short- to medium-range coupled prediction. *J. Operat. Oceanogr.* 8, s239–s258. doi: 10.1080/ 1755876X.2015.1049875

Buehner, M., Caya, A., Pogson, L., Carrieres, T., and Pestieau, P. (2013). A new Environment Canada regional ice analysis system. *Atmos. Ocean* 51, 18–34. doi: 10.1080/07055900.2012.747171

Buizza, R., Brönnimann, S., Haimberger, L., Laloyaux, P., Martin, M. J., Fuentes, M., et al. (2018). The EU-FP7 ERA-CLIM2 project contribution to advancing science and production of earth system climate reanalyses. *Bull. Am. Meteorol. Soc.* 99, 1003–1014. doi: 10.1175/bams-d-17-0199.1

Bunzel, F., Notz, D., Baehr, J., Müller, W. A., and Fröhlich, K. (2016). Seasonal climate forecasts significantly affected by observational uncertainty of Arctic sea ice concentration. *Geophys. Res. Lett.* 43, 852–859. doi: 10.1002/2015gl06 6928

Calder, J., Proshutinsky, A., Carmack, E., Ashik, I., Loeng, H., Key, J., et al. (2010). “Community white paper: an integrated international approach to arctic ocean observations for society (a legacy of the international polar year),” in *Proceedings of OceanObs’09: Sustained Ocean Observations and Information for Society*, Vol. 2, eds J. Hall, D. E. Harrison, and D. Stammer (Venice: ESA Publication), 306. doi: 10.5270/OceanObs09.cwp.14

Carrieres, T., Caya, A., Posey, P., Metzger, E. J., Bertino, L., Melsom, A., et al. (2017). *Sea Ice Analysis and Forecasting: Towards an Increased Reliance on Automated Prediction Systems*, eds T. Carrieres, M. Buehner, J.-F. Lemieux, and L. T. Pedersen (Cambridge: Cambridge University Press), 144–173. doi: 10.1017/9781108277600

Carse, F., Martin, M. J., Sellar, A., and Blockley, E. W. (2015). Impact of assimilating temperature and salinity measurements by animal-borne sensors on FOAM ocean model fields. *Q. J. R. Meteorol. Soc.* 141, 2934–2943. doi: 10.1002/qj.2613

Castro, S. L., Wick, G. A., and Steele, M. (2016). Validation of satellite sea surface temperature analyses in the Beaufort Sea using UpTempO buoys. *Remote Sens. Environ.* 187, 458–475. doi: 10.1016/j.rse.2016.10.035

Cavaleri, L., Fox-Kemper, B., and Hemer, M. (2012). Wind waves in the coupled climate system. *Bull. Am. Meteorol. Soc.* 92, 1651–1661. doi: 10.1175/bams-d- 11-00170.1

Cavaleri, D. J., and Parkinson, C. L. (2012). Arctic sea ice variability and trends, 1979-2010. *Cryosphere* 6, 881–889. doi: 10.5194/tc-6-881-2012

Chen, Z., Liu, J., Song, M., Yang, Q., and Xu, S. (2017). Impacts of assimilating satellite sea ice concentration and thickness on Arctic sea ice prediction in the NCEP Climate Forecast System. *J. Clim.* 30, 8429–8446. doi: 10.1175/jcli-d-17- 0093.1

Chevallier, M., Smith, G., Lemieux, J.-F., Dupont, F., Forget, G., Fujii, Y., et al. (2016). Intercomparison of the Arctic sea ice cover in global ocean-sea ice reanalyses from the ORA-IP project. *Clim. Dyn. Spec. Issue Ocean Reanalysis* 49, 1107–1136. doi: 10.1007/s00382-016-2985-y

Claustre, H., Bishop, J., Boss, E., Stewart, B., Berthon, J.-F., Coatanoan, C., et al. (2010). “Bio-optical profiling floats as new observational tools for biogeochemical and ecosystem studies: potential synergies with ocean color remote sensing,” in *Proceedings of OceanObs’09: Sustained Ocean Observations and Information for Society*, Vol. 2, eds J. Hall, D. E. Harrison, and D. Stammer (Venice: ESA Publication), 306.

Cole, S. T., Thwaites, F. T., Krishfield, R. A., and Toole, J. M. (2015). “Processing of velocity observations from ice-tethered profilers,” in *Proceedings of the Oceans 2015 MTS/IEEE*, Washington, DC.

Cole, S. T., Timmermans, M. L., Toole, J. M., Krishfield, R. A., and Thwaites, F. T. (2014). Ekman veering, internal waves, and turbulence observed under Arctic sea ice. *J. Phys. Oceanogr.* 44, 1306–1328. doi: 10.1175/jpo-d-12-0191.1

Cummings, J. A., and Smedstad, O. M. (2014). Ocean data impacts in global HYCOM. *J. Atmos. Ocean. Technol.* 31, 1771–1791. doi: 10.1175/JTECH-D-14- 00011.1

Curry, B., Lee, C. M., Petrie, B., Moritz, R. E., and Kwok, R. (2014). Multiyear volume, liquid freshwater, and sea ice transports through Davis Strait, 2004–10. *J. Phys. Oceanogr.* 44, 1244–1266. doi: 10.1175/jpo-d-13-0177.1

Dash, P., Ignatov, A., Kihai, Y., and Sapper, J. (2010). The SST quality monitor (SQUAM). *J. Atmos. Ocean. Technol.* 27, 1899–1917. doi: 10.1175/ 2010JTECHO756.1

Davidson, F., Chassignet, E., Vinayachandran, P. N., Lu, Y., Smith, G. C., and Zhu, X. (2019). Synergies in operational oceanography: the intrinsic need for sustained ocean observations. *Front. Mar. Sci.* 6:450. doi: 10.3389/fmars.2019. 00450 doi: 10.1175/jtech-d-14-00011.1

Davidson, F., Smith, G. C., Lu, Y., and Woodbury, S. (2013). Operational atmosphere-ocean-ice prediction systems in Canada: providing decision- enabling marine environmental information to end users. *Can. Ocean Sci. Newsl.* 70, 2–5.

Dong, S., Gille, S. T., and Sprintall, J. (2007). An assessment of the Southern Ocean mixed-layer heat budget. *J. Clim.* 20, 4425–4442. doi: 10.1175/jcli4259.1

Dosser, H. V., Rainville, L., and Toole, J. M. (2014). Near-inertial internal wave field in the Canada basin from ice-tethered profilers. *J. Phys. Oceanogr.* 44, 413–426. doi: 10.1175/JPO-D-13-0117.1

du Plessis, M., Swart, S., Ansorge, I. J., Mahadevan, A., and Thompson, A. F. (2019). Southern Ocean seasonal stratification delayed by wind-front interaction at submesoscales. *J. Phys. Oceanogr.* 49, 1035–1053. doi: 10.1175/JPO-D-18- 0136.1

Fedak, M. A. (2013). The impact of animal platforms on polar ocean observation. *Deep Sea Res. Part II Top. Stud. Oceanogr.* 88–89, 7–13. doi: 10.1016/j.dsr2. 2012.07.007

Fennel, K., and Greenan, B. (2017). *Taking the Ocean's Pulse: A Vision for the Canadian Biogeochemical Argo Program*. Available at: https://doi.org/10.13155/ 52451

Foppert, A., Donohue, K. A., Randolph Watts, D., and Tracey, K. L. (2017). Eddy heat flux across the Antarctic Circumpolar Current estimated from sea surface height standard deviation. *J. Geophys. Res. Oceans* 122, 6947–6964. doi: 10.1002/2017JC012837

Freitag, L., Ball, K., Partan, J., Koski, P., and Singh, S. (2015). “Long range acoustic communications and navigation in the Arctic,” in *Proceedings of the OCEANS’15 Conference and Exhibition*, (Washington, DC: MTS), 1–5.

Fujii, Y., Remy, E., Zuo, H., Oke, P. R., Halliwell, G. R., Gasparin, F., et al. (2019). Observing system evaluation based on ocean data assimilation and prediction systems: on-going challenges and future vision for designing/supporting ocean observational networks. *Front. Mar. Sci.* 6:417. doi: 10.3389/fmars.2019.00417

Gentemann, C. L., Meissner, T., and Wentz, F. J. (2010). Accuracy of satellite sea surface temperatures at 7 and 11 GHz. *IEEE Trans. Geosci. Remote Sens.* 48, 1009–1018. doi: 10.1109/tgrs.2009.2030322
Gille, S., Josey, S., and Swart, S. (2016). New approaches for air-sea fluxes in the Southern Ocean. *Eos* 97. doi: 10.1029/2016EO52243

Greenert, J. W. (2014). *The United States Navy Arctic Roadmap for 2014 to 2030*. Washington, DC: Office of the Chief of Naval Operations.

Guémas, V., Blanchard-Wrigglesworth, E., Chevallier, M., Day, J., Déqué, M., Doblas-Reyes, F., et al. (2016a). A review on Arctic sea ice predictability and prediction on seasonal-to-decadal timescales. *Q. J. R. Meteorol. Soc.* 142, 546–561. doi: 10.1002/qj.2401

Guémas, V., Chevallier, M., Déqué, M., Bellprat, O., and Doblas-Reyes, F. (2016b). Impact of sea ice initialization on sea ice and atmosphere prediction skill on seasonal timescales. *Geophys. Res. Lett.* 43, 3889–3896. doi: 10.1002/2015GL066626

Guerreiro, K., Fleury, S., Zakharova, E., Rémy, F., and Maisongrande, P. (2016). Potential for estimation of snow depth on Arctic sea ice from CryoSat-2 and SARAL/AltiKa missions. *Remote Sens. Environ.* 11, 2059–2073.

Hardt, M., and Scherbaum, F. (1994). The design of optimum networks for aftershock recordings. *Geophys. J. Int.* 117, 716–726. doi: 10.1111/j.1365-246x.1994.tb02464.x

Hebert, D. A., Allard, R. A., Metzger, E. J., Posey, P. G., Preller, R. H., Wallcraft, A. J., et al. (2015). Short-term sea ice forecasting: an assessment of ice concentration and ice drift forecasts using the U.S. Navy's Arctic Cap Nowcast/Forecast System. *J. Geophys. Res. Oceans* 120, 8327–8345. doi: 10.1002/2015JC011283

Helfrich, S. R., McNamara, D., Ramsay, B. H., Baldwin, T., and Kasheta, T. (2007). Enhancements to, and forthcoming developments in the Interactive Multisensor Snow and Ice Mapping System (IMS). *Hydrol. Process.* 21, 1576–1586. doi: 10.1002/hyp.6720

Hendricks, S., Ricker, R., and Helm, V. (2016). User guide-AWI CryoSat-2 sea ice thickness data product (v1.2). Alfred Wegener Institute Technical Report.

Hogan, T. F., Liu, M., Ridout, J. A., Peng, M. S., Whitcomb, T. R., Ruston, B. C., et al. (2014). The navy global environmental model. *Oceanography* 27, 116–125.

Hoke, W., Werner, K., Goessling, H., and Jung, T. (2018). "Engaging forecast users during the year of polar prediction," in *Building Bridges at the Science-Stakeholder Interface*, ed. G. Krause (Cham: Springer), 47–53. doi: 10.1007/978-3-319-75919-7_7

Howe, B. M., Miksis-Olds, J., Rehm, E., Sagen, H., Worcester, P. F., and Haralabus, G. (2019). Observing the oceans acoustically. *Front. Mar. Sci.* 6:426. doi: 10.3389/fmars.2019.00426

Hunke, E., and Lipscomb, W. (2008). "CICE: The Los Alamos Sea Ice Model documentation and software user's manual version 4.0," in *Technical Report LA-CC-06-012, AWI CryoSat-2 Sea Ice Thickness Data Product (v1.2) User Guide*, eds S. Hendricks, R. Ricker, and V. Helm (Los Alamos, NM: Los Alamos National Laboratory).

Hunke, E. C., Lipscomb, W. H., Turner, A. K., Jeffery, N., and Elliott, S. (2015). *CICE: The Los Alamos Sea ice Model Documentation and Software User's Manual Version 5 (Tech. Rep. LA-CC-06-012)*. Los Alamos, NM: Los Alamos National Laboratory.

Hyder, P., Edwards, J. M., Allan, R. P., Hewitt, H. T., Bracegirdle, T. J., Gregory, J. M., et al. (2018). Critical Southern Ocean climate model biases traced to atmospheric model cloud errors. *Nat. Commun.* 9:3625. doi: 10.1038/s41467-018-05634-2

Ilicak, M., Drange, H., Wang, Q., Gerdes, R., Aksenov, Y., Bailey, D., et al. (2016). An assessment of the Arctic Ocean in a suite of interannual CORE-II simulations. Part III: hydrography and fluxes. *Ocean Model.* 100, 141–161. doi: 10.1016/j.ocemod.2016.02.004

Ilyina, T., Six, K. D., Segschneider, J., Maier-Reimer, E., Li, H., and NezRiboni, I. (2013). Global ocean biogeochemistry model HAMOCC: model architecture and performance as component of the MPI-Earth system model in different CMIP5 experimental realizations. *J. Adv. Model. Earth Syst.* 5, 287–315. doi: 10.1029/2012MS000178

Inoue, J., Enomoto, T., Miyoshi, T., and Yamane, S. (2009). Impact of observations from Arctic drifting buoys on the reanalysis of surface fields. *Geophys. Res. Lett.* 36:L08501.

Inoue, J., Yamazaki, A., Ono, J., Dethloff, K., Maturilli, M., Neuber, R., et al. (2015). Additional Arctic observations improve weather and sea-ice forecasts for the Northern Sea Route. *Sci. Rep.* 5:16868. doi: 10.1038/srep16868

Islam, F., DeGrandpre, M. D., Beatty, C. M., Timmermans, M. L., Krishfield, R. A., Toole, J. M., et al. (2017). Sea surface pCO₂ and O₂ dynamics in the

partially ice-covered Arctic Ocean. *J. Geophys. Res.* 122, 1425–1438. doi: 10.1002/2016JC012162

Ivanova, N., Johannessen, O. M., Pedersen, L. T., and Tonboe, R. T. (2014). Retrieval of Arctic sea ice parameters by satellite passive microwave sensors: a comparison of eleven sea ice concentration algorithms. *IEEE Trans. Geosci. Remote Sens.* 52, 7233–7246. doi: 10.1109/tgrs.2014.2310136

Jackson, K., Wilkinson, J., Maksym, T., Meldrum, D., Beckers, J., Haas, C., et al. (2013). A novel and low-cost sea ice mass balance buoy. *J. Atmos. Ocean. Technol.* 30, 2676–2688. doi: 10.1175/jtech-d-13-00058.1

Johnson, N., Fidel, M., Danielsen, F., Iversen, L., Poulsen, M. K., Hauser, D., et al. (eds). (2018). "INTAROS community-based monitoring experience exchange workshop report: Canada. Workshop organized as a contribution to INTAROS, December 11-12, 2017," in *Exchange for Local Observations and Knowledge of the Arctic (ELOKA), Nordic Foundation for Development and Ecology (NORDECO), Yukon River Inter-Tribal Watershed Council (YRITWC)*, (Québec: International Arctic Research Center), 27.

Josey, S. A., Gulev, S., and Yu, L. (2013). "Exchanges through the ocean surface," in *Ocean Circulation and Climate: A 21st Century Perspective*, eds G. Siedler, J. Church, W. J. Gould, and S. Griffies (Cambridge, MA: Academic Press), 115–140. doi: 10.1016/b978-0-12-391851-2.00005-2

Jung, T., Doblas-Reyes, F., Goessling, H., Guemas, V., Bitz, C., Buontempo, C., et al. (2015). Polar lower-latitude linkages and their role in weather and climate prediction. *Bull. Am. Meteorol. Soc.* 96, ES197–ES200. doi: 10.1175/BAMS-D-15-00121.1

Jung, T., Gordon, N. D., Bauer, P., Bromwich, D. H., Chevallier, M., Day, J. J., et al. (2016). Advancing polar prediction capabilities on daily to seasonal time scales. *Bull. Am. Meteorol. Soc.* 97, 1631–1647. doi: 10.1175/BAMS-D-14-00246.1

Jungclaus, J. H., Fischer, N., Haak, H., Lohmann, K., Marotzke, J., Matei, D., et al. (2013). Characteristics of the ocean simulations in MPIOM, the ocean component of the MPI-Earth system model. *J. Adv. Model. Earth Syst.* 5, 422–446. doi: 10.1002/jame.20023

Kaleschke, L., Tian-Kunze, X., Maß, N., Mäkynen, M., and Drusch, M. (2012). Sea ice thickness retrieval from SMOS brightness temperatures during the Arctic freeze-up period. *Geophys. Res. Lett.* 39:L05501.

Kaminski, T., Kauker, F., Eicken, H., and Karcher, M. (2015). Exploring the utility of quantitative network design in evaluating Arctic sea ice thickness sampling strategies. *Cryosphere* 9, 1721–1733. doi: 10.5194/tc-9-1721-2015

Kaminski, T., Kauker, F., Toudal Pedersen, L., Vobeck, M., Haak, H., Niederdrenk, L., et al. (2017). Arctic mission benefit analysis: impact of sea ice thickness, freeboard, and snow depth products on sea ice forecast performance. *Cryosphere Discuss.* 2017, 1–39. doi: 10.5194/tc-2017-249

Kaminski, T., and Mathieu, P.-P. (2017). Reviews and syntheses: flying the satellite into your model: on the role of observation operators in constraining models of the Earth system and the carbon cycle. *Biogeosciences* 14, 2343–2357. doi: 10.5194/bg-14-2343-2017

Kaminski, T., and Rayner, P. J. (2017). Assisting the evolution of the observing system for the carbon cycle through quantitative network design. *Biogeosci. Discuss.* 14, 1–22. doi: 10.5194/bg-2017-168

Karvonen, J. (2012). Operational SAR-based sea ice drift monitoring over the Baltic Sea. *Ocean Sci.* 8, 473–483. doi: 10.5194/os-8-473-2012

Kikuchi, T., Inoue, J., and Langevin, D. (2007). Argo-type profiling float observations under the Arctic multiyear ice. *Deep Sea Res. Part I Oceanogr. Res. Pap.* 54, 1675–1686. doi: 10.1016/j.dsr.2007.05.011

Klatt, O., Boebel, O., and Fahrbach, E. (2007). A profiling float's sense of ice. *J. Atmos. Ocean. Technol.* 24, 1301–1308. doi: 10.1175/jtech2026.1

Korosov, A., and Rampal, P. (2017). A combination of feature tracking and pattern matching with optimal parametrization for sea ice drift retrieval from SAR data. *Remote Sens.* 9:258. doi: 10.3390/rs9030258

Krishfield, R., Toole, J., Proshutinsky, A., and Timmermans, M.-L. (2008). Automated Ice-Tethered Profilers for seawater observations under pack ice in all seasons. *J. Atmos. Ocean. Technol.* 25, 2091–2105. doi: 10.1175/2008jtech0587.1

Kwok, R., and Cunningham, G. F. (2015). Variability of Arctic sea ice thickness and volume from CryoSat-2. *Philos. Trans. R. Soc. A Math. Phys. Eng. Sci.* 373:20140157. doi: 10.1098/rsta.2014.0157

Kwok, R., Kurtz, N. T., Brucker, L., Ivanoff, A., Newman, T., Farrell, S. L., et al. (2017). Intercomparison of snow depth retrievals over Arctic sea ice from
radar data acquired by Operation IceBridge. *Cryosphere* 11, 2571–2593. doi: 10.5194/tc-11-2571-2017

Kwok, R., and Sulsky, D. (2010). Arctic Ocean sea ice thickness and kinematics: satellite retrievals and modeling. *Oceanography* 23, 134–143. doi: 10.5670/oceanog.2010.11

Lagunas, J., Marec, C., Leymarie, É, Penkerc'h, C., Rehm, E., Desaulniers, P., et al. (2018). Sea-ice detection for autonomous underwater vehicles and oceanographic lagrangian platforms by continuous-wave laser polarimetry. *Ocean Sens. Monitor.* X 10631:106310W. doi: 10.1117/12.230 9571

Laney, S. R., Krishfield, R. A., and Toole, J. M. (2017). The euphotic zone under Arctic Ocean sea ice: vertical extents and seasonal trends. *Limnol. Oceanogr.* 62, 1910–1934. doi: 10.1002/lno.10543

Laney, S. R., Krishfield, R. A., Toole, J. M., Hammar, T. R., Ashjian, C. J., and Timmermans, M. L. (2014). Assessing algal biomass and bio-optical distributions in perennially ice-covered polar ocean ecosystems. *Polar Sci.* 8, 73–85. doi: 10.1016/j.polar.2013.12.003

Lavergne, T., Eastwood, S., Teffah, Z., Schyberg, H., and Breivik, L.-A. (2010). Sea ice motion from low-resolution satellite sensors: an alternative method and its validation in the Arctic. *J. Geophys. Res. Atmos.* 115:C10032. doi: 10.1029/2009JC005958

Le Traon, P. Y., Ali, A., Alvarez Fanjul, E., Aouf, L., Axell, L., and Aznar, R. (2017). The copernicus marine environmental monitoring service: main scientific achievements and future prospects. *Spec. Issue Mercatorol. Ocean J.* 56, 1–101.

Lee, C. M., Cole, S., Doble, M., Guthrie, J. D., Harper, S., MacKinnon, J., et al. (2016). “Stratified ocean dynamics in the arctic: science and experiment plan. APL-UW TR 1601,” in *Technical Report, Applied Physics Laboratory*, (Seattle, WA: University of Washington), 46.

Lee, C. M., and Gobat, J. (2008). *Acoustic Navigation and Communications for High-Latitude Ocean Research*. Seattle, WA: University of Washington.

Lee, C. M., Starkweather, S., Eicken, H., Timmermans, M.-L., Wilkinson, J., Sandven, S., et al. (2019). A framework for the development, design and implementation of a sustained arctic ocean observing system. *Front. Mar. Sci.* 6:451. doi: 10.3389/fmars.2019.00451

Lei, R., Cheng, B., Heil, P., Vihma, T., Wang, J., Ji, Q., et al. (2018). Seasonal and interannual variations of sea ice mass balance from the central arctic to the Greenland sea. *J. Geophys. Res. Oceans* 123, 2422–2439. doi: 10.1002/2017jc013548

Lemieux, J. F., Beaudoin, C., Dupont, F., Roy, F., Smith, G. C., Shlyaeva, A., et al. (2016). The Regional Ice Prediction System (RIPS): verification of forecast sea ice concentration. *Q. J. R. Meteorol. Soc.* 142, 632–643. doi: 10.1002/qj.2526

Liao, Z., Cheng, B., Zhao, J., Vihma, T., Jackson, K., Yang, Q., et al. (2018). Snow depth and ice thickness derived from SIMBA ice mass balance buoy data using an automated algorithm. *Int. J. Digit. Earth* 12, 1–18. doi: 10.1080/17538947. 2018.1545877

Lien, V. S., Hjøllo, S. S., Skogen, M. D., Svendsen, E., Wehde, H., Bertino, L., et al. (2016). An assessment of the added value from data assimilation on modelled Nordic Seas hydrography and ocean transports. *Ocean Model.* 99, 43–59. doi: 10.1016/j.ocemod.2015.12.010

Lindsay, R. W., and Zhang, J. (2006). Assimilation of ice concentration in an ice-ocean model. *J. Atmos. Ocean. Technol.* 23, 742–749. doi: 10.1175/jtech1 871.1

Lisaeter, K. A., Rosanova, J., and Evensen, G. (2003). Assimilation of ice concentration in a coupled ice-ocean model, using the Ensemble Kalman filter. *Ocean Dyn.* 53, 368–388. doi: 10.1007/s10236-003-0049-4

Longépé, N., Thibaut, P., Vadaine, R., Poisson, J. C., Guillot, A., Boy, F., et al. (2019). Comparative evaluation of sea ice lead detection based on SAR imagery and altimeter data. *IEEE Trans. Geosci. Remote Sens.* 57, 4050–4061. doi: 10.1109/tgrs.2018.2889519

Macelloni, G., Brogioni, M., Montomoli, F., Leduc-Leballeur, M., De Carolis, G., Kaleschke, L., et al. (2017). “Preliminary study for a spaceborne ultrawideband microwave radiometer for the monitoring of cryosphere elements: the Cryorad project,” in *Proceedings of the 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS)*, Fort Worth, TX, 1185–1188.

Maksym, T. (2019). Arctic and Antarctic Sea ice change: contrasts, commonalities, and causes. *Annu. Rev. Mar. Sci.* 11, 187–213. doi: 10.1146/annurev-marine- 010816-060610

Mariotti, A., Ruti, P. M., and Rixen, M. (2018). Progress in subseasonal to seasonal prediction through a joint weather and climate community effort. *NPJ Clim. Atmos. Sci.* 1:4.

Maslowski, W., Clement Kinney, J., Higgins, M., and Roberts, A. (2012). The future of Arctic sea ice. *Annu. Rev. Earth Planet. Sci.* 40, 625–654.

Massom, R. A., Scambos, T. A., Bennetts, L. G., Reid, P., Squire, V. A., and Stammerjohn, S. E. (2018). Antarctic Ice shelf disintegration triggered by sea ice loss and ocean swell. *Nature* 558, 383–389. doi: 10.1038/s41586-018-0212-1

Massonnet, F., Mathiot, P., Fichefet, T., Goosse, H., Beatty, C. K., Vancoppenolle, M., et al. (2013). A model reconstruction of the Antarctic sea ice thickness and volume changes over 1980–2008 using data assimilation. *Ocean Model.* 64, 67–75. doi: 10.1016/j.ocemod.2013.01.003

McPhee, M. G. (2013). Intensification of geostrophic currents in the Canada basin, Arctic ocean. *J. Clim.* 26, 3130–3138. doi: 10.1175/JCLI-D-12-00289.1

Melsom, A., Bertino, L., and Sutherland, G. (2015). *Quality Information Document for the Arctic Ocean Physical Analysis and Forecast Product ARCTIC ANALYSIS FORECAST PHYS 002 001 A*. Available at: http://marine.copernicus.eu/ documents/QUID/CMEMS-ARC-QUID-002-001a.pdf

Metzger, E. J., Smedstad, O. M., Thoppil, P. G., Hurlburt, H. E., Cummings, J. A., Wallcraft, A. J., et al. (2014). US Navy operational global ocean and Arctic ice prediction systems. *Oceanography* 27, 32–43. doi: 10.5670/oceanog.2014.66

Mikhalevsky, P. N., Sagen, H., Worcester, P. F., Baggeroer, A. B., Orcutt, J., Moore, S. E., et al. (2015). Multipurpose acoustic networks in the integrated arctic ocean observing system. *Arctic* 68, 11–27.

Mogensen, K. S., Magnusson, L., and Bidlot, J. R. (2017). Tropical cyclone sensitivity to ocean coupling in the ECMWF coupled model. *J. Geophys. Res. Oceans* 122, 4392–4412. doi: 10.1002/2017jc012753

Monteiro, P. M. S., Gregor, L., Lévy, M., Maener, S., Sabine, C. L., and Swart, S. (2015). Seasonal-scale robotics experiment reveals the contribution made by fine scale dynamics to the uncertainties and biases in the mean seasonal CO₂ flux in the Southern Ocean. *Geophys. Res. Lett.* 42, 8507–8514. doi: 10.1002/ 2015GL066009

Mu, L., Yang, Q., Losch, M., Losa, S. N., Ricker, R., Nerger, L., et al. (2018). Improving sea ice thickness estimates by assimilating CryoSat-2 and SMOS sea ice thickness data simultaneously. *Q. J. R. Meteorol. Soc.* 144, 529–538. doi: 10.1002/qj.3225

Murray, M. S., Sankar, R. D., and Ibarguchi, G. (2018). *The Arctic Observing Summit, Background and Synthesis of Outcomes 2013-2016*. Calgary: International Study of Arctic Change.

Newman, L., Heil, P., Trebilco, R., Katsumata, K., Constable, A. J., van Wijk, E., et al. (2019). Delivering sustained, coordinated and integrated observations of the southern ocean for global impact. *Front. Mar. Sci.* 6:433. doi: 10.3389/fmars. 2019.00433

Niederdrenk, L. (2013). *The Arctic Hydrologic Cycle and its Variability in a Regional Coupled Climate Model*. Ph.D. thesis, University of Hamburg, Hamburg, 1–186.

Ogle, S. E., Tamsitt, V., Josey, S. A., Gille, S. T., Ceroveëki, I., Talley, L. D., et al. (2018). Episodic Southern Ocean heat loss and its mixed layer impacts revealed by the furthest south multi-year surface flux mooring. *Geophys. Res. Lett.* 45, 5002–5010. doi: 10.1029/2017GL076909

Oke, P. R., Larnicol, G., Fujii, Y., Smith, G. C., Lea, D. J., Guinehut, S., et al. (2015). Assessing the impact of observations on ocean forecasts and reanalyses: part 1, Global studies. *J. Operat. Oceanogr.* 8, s49–s62. doi: 10.1080/1755876X.2015. 1022067

Osborne, E., Richter-Menge, J., and Jeffries, M. (eds). (2018). *Arctic Report Card 2018*. Available at: https://www.arctic.noaa.gov/Report-Card (accessed July 26, 2019).

Parkinson, C. L., and Comiso, J. C. (2013). On the 2012 record low Arctic sea ice cover: combined impact of preconditioning and an August storm. *Geophys. Res. Lett.* 40, 1356–1361. doi: 10.1002/grl.50349

Patilea, C., Heygster, G., Huntemann, M., and Sreen, G. (2019). Combined SMAP-SMOS thin sea ice thickness retrieval. *Cryosphere* 13, 675–691. doi: 10.5194/tc-13-675-2019

Pedersen, L. T., Saldo, R., and Fenger-Nielsen, R. (2015). “Sentinel-1 results: sea ice operational monitoring,” in *Proceedings of the 2015 IEEE International Geoscience and Remote Sensing Symposium (IGARSS)*, Milan, 2828–2831. doi: 10.1109/IGARSS.2015.7326403

Penny, S. G., and Hamill, T. M. (2017). Coupled data assimilation for integrated earth system analysis and prediction. *Bull. Am. Meteorol. Soc.* 97, ES169–ES172.
Persson, P. O. G., Fairall, C. W., Andreas, E. L., Guest, P. S., and Perovich, D. K. (2002). Measurements near the Atmospheric Surface Flux Group tower at SHEBA: near-surface conditions and surface energy budget. *J. Geophys. Res.* 107:8045.

Posey, P. G., Metzger, E. J., Wallcraft, A. J., Hebert, D. A., Allard, R. A., Smedstad, O. M., et al. (2015). Assimilating high horizontal resolution sea ice concentration data into the US Navy's ice forecast systems: arctic Cap Nowcast/Forecast System (ACNFS) and the Global Ocean Forecast System (GOFS 3.1). *Cryosphere Discuss.* 9, 2339–2365. doi: 10.5194/tcd-9-2339-2015

PPP Steering Group (2013). *WWRP Polar Prediction Project Implementation Plan*. WMO Report No. WWRP/PPP 2-2013, 72.

PPP Steering Group (2014). *WWRP Polar Prediction Project Implementation Plan for the Year of Polar Prediction (YOPP)*. WMO Report No. WWRP/PPP 3-2014, 53.

Prandi, P., Ablain, M., Cazenave, A., and Picot, N. (2012). A new estimation of mean sea level in the Arctic Ocean from satellite altimetry. *Mar. Geod.* 35, 61–81. doi: 10.1080/01490419.2012.718222

Provost, C., Pelon, J., Sennéchael, N., Calzas, M., Blouzon, F., Desautez, A., et al. (2015). IAOOS (ice-atmosphere-arctic ocean observing system, 2011-2019). *Mercator Ocean Q. Newsl.* 51, 13–15.

Pullen, J., Allard, R., Seo, H., Miller, A. J., Chen, S., Pezzi, L. P., et al. (2017). Coupled ocean-atmosphere forecasting at short and medium time scales. *J. Mar. Res.* 75, 877–921. doi: 10.1357/002224017823523991

Rabatel, M., Rampal, P., Carrassi, A., Bertino, L., and Jones, C. K. (2018). Impact of rheology on probabilistic forecasts of sea ice trajectories: application for search and rescue operations in the Arctic. *Cryosphere* 12, 935–953. doi: 10.5194/tc-12-935-2018

Rabe, B., Karcher, M., Schauer, U., Toole, J. M., Krishfield, R. A., Pisarev, S., et al. (2011). An assessment of pan-Arctic Ocean freshwater content changes from the 1990s to the IPY period. *Deep Sea Res. I* 58, 173–185. doi: 10.1016/j.dsr.2010.12.002

Rayner, P. J., Enting, I. G., and Trudinger, C. M. (1996). Optimizing the CO₂ observing network for constraining sources and sinks. *Tellus* 48B, 433–444. doi: 10.3402/tellusb.v48i4.15924

Revee, K., Boebel, O., Kanzow, T., Strass, V., Rohardt, G., and Fahrbach, E. (2016). A gridded data set of upper-ocean hydrographic properties in the Weddell Gyre obtained by objective mapping of Argo float measurements. *Earth Syst. Sci. Data* 8, 15–40. doi: 10.5194/essd-8-15-2016

Rehm, E., Dushaw, B. D., Freitag, L., Heaney, K. D., Carr, S., Dakin, T., et al. (2018). *BBANC: A Baffin Bay Acoustic Navigation and Communication System Feasibility Study*. Available at: https://gitlab.com/Takuvik/BBANC/wikis/home (accessed June 1, 2018).

Richards, C., Pittman, M., Phelan, K., Nudds, S., and Hamilton, J. (2017). "The barrow strait real time observatory: under-ice monitoring in the Canadian high arctic," in *Proceedings of the International Conference on Underwater Networks & Systems*, (New York, NY: ACM), 8.

Richter-Menge, J. A., Perovich, D. K., Elder, B. C., Claffey, K., Rigor, I., and Ortmeyer, M. (2006). Ice mass-balance buoys: a tool for measuring and attributing changes in the thickness of the Arctic sea-ice cover. *Ann. Glaciol.* 44, 205–210. doi: 10.3189/172756406781811727

Ricker, R., Hendricks, S., Helm, V., Skourup, H., and Davidson, M. (2014). Sensitivity of CryoSat-2 Arctic sea-ice freeboard and thickness on radar-waveform interpretation. *Cryosphere* 8, 1607–1622. doi: 10.5194/tc-8-1607-2014

Riser, S. C., Freeland, H. J., Roemmich, D., Wijffels, S., Troisi, A., Belbéoch, M., et al. (2016). Fifteen years of ocean observations with the global Argo array. *Nat. Clim. Change* 6, 145–153.

Roemmich, D., and Gilson, J. (2009). The 2004–2008 mean and annual cycle of temperature, salinity, and steric height in the global ocean from the Argo Program. *Prog. Oceanogr.* 82, 81–100. doi: 10.1016/j.pocean.2009.03.004

Roemmich, D., Johnson, G. C., Riser, S., Davis, R., Gilson, J., Owens, W. B., et al. (2009). The Argo Program: observing the global ocean with profiling floats. *Oceanography* 22, 34–43. doi: 10.5670/oceanog.2009.36

Roquet, F., Wunsch, C., Forget, G., Heimbach, P., Guinet, C., Reverdin, G., et al. (2013). Estimates of the Southern Ocean general circulation improved by animal-borne instruments. *Geophys. Res. Lett.* 40, 6176–6180. doi: 10.1002/2013GL058304

Sakov, P., Counillon, F., Bertino, L., Lisaeter, K. A., Oke, P. R., and Korablev, A. (2012). TOPAZ4: an ocean-sea ice data assimilation system for the North Atlantic and Arctic. *Ocean Sci.* 8, 633–656. doi: 10.5194/os-8-633-2012

Schmidt, K., Swart, S., Reason, C., and Nicholson, S. (2017). Evaluation of satellite and reanalysis wind products with in situ Wave Glider wind observations in the Southern Ocean. *J. Atmos. Ocean. Technol.* 34, 2551–2568. doi: 10.1175/JTECH-D-17-0079.1

Schulz, E. W., Josey, S. A., and Verein, R. (2012). First air-sea flux mooring measurements in the Southern Ocean. *Geophys. Res. Lett.* 39:L16606.

Serreze, M. C., and Barry, R. G. (2011). Processes and impacts of Arctic amplification: a research synthesis. *Glob. Planet. Change* 77, 85–96. doi: 10.1016/j.gloplacha.2011.03.004

Shibley, N. C., Timmermans, M. L., Carpenter, J. R., and Toole, J. M. (2017). Spatial variability of the Arctic Ocean's double-diffusive staircase. *J. Geophys. Res. Oceans* 122, 980–994. doi: 10.1002/2016jc012419

Shupe, M., de Boer, G., Dethloff, K., Hunke, E., Maslowski, W., McComiskey, A., et al. (2018). *The Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAIC) Atmosphere Science Plan*. DOE Report No. DOE/SC-ARM-18-005.

Sigmond, M., Reader, M. C., Flato, G. M., Merryfield, W. J., and Tivy, A. (2016). Skillful seasonal forecasts of Arctic sea ice retreat and advance dates in a dynamical forecast system. *Geophys. Res. Lett.* 43, 12457–12465.

Smith, G. C., Bélanger, J. M., Roy, F., Pellerin, P., Ritchie, H., Onu, K., et al. (2018). Impact of coupling with an ice-ocean model on global medium-range NWP forecast skill. *Mon. Weather Rev.* 146, 1157–1180. doi: 10.1175/MWR-D-17-0157.1

Smith, G. C., Roy, F., Belanger, J. M., Dupont, F., Lemieux, J. F., Beaudoin, C., et al. (2013). "Small-scale ice-ocean-wave processes and their impact on coupled environmental polar prediction," in *Proceedings of the ECMWF-WWRP/THORPEX Polar Prediction Workshop, 24-27 June 2013*, (Reading: ECMWF).

Smith, G. C., Roy, F., Reszka, M., Surcel Colan, D., He, Z., Deacu, D., et al. (2016). Sea ice forecast verification in the Canadian global ice ocean prediction system. *Q. J. R. Meteorol. Soc.* 142, 659–671. doi: 10.1002/qj.2555

Soltwedel, T., Schauer, U., Boebel, O., Nöthig, E. M., Bracher, A., Metfies, K., et al. (2013). "FRAM - frontiers in arctic marine monitoring visions for permanent observations in a gateway to the arctic ocean," in *Proceedings of the 2013 MTS/IEEE*, Bergen. doi: 10.1109/OCEANS-Bergen.2013.6608008

Stark, J. D., Ridley, J., Martin, M., and Hines, A. (2008). Sea ice concentration and motion assimilation in a sea ice-ocean model. *J. Geophys. Res. Oceans* 113:C05S91.

Swart, S., Gille, S. T., Delille, B., Josey, S., Mazloff, M., Newman, L., et al. (2019). Constraining Southern Ocean air-sea-ice fluxes through enhanced observations. *Front. Mar. Sci.* 6:421. doi: 10.3389/fmar.2019.00421

Swart, S., Thomalla, S. J., and Monteiro, P. M. S. (2015). The seasonal cycle of mixed layer dynamics and phytoplankton biomass in the Sub-Antarctic Zone: a high-resolution glider experiment. *J. Mar. Syst.* 147, 103–115. doi: 10.1016/j.jmarsys.2014.06.002

Taylor, P., Hegyi, B., Boeke, R., and Boisvert, L. (2018). On the increasing importance of air-sea exchanges in a thawing Arctic: a review. *Atmosphere* 9:41. doi: 10.3390/atmos9020041

Thomson, J., and Girton, J. (2017). Sustained measurements of Southern Ocean air-sea coupling from a wave glider autonomous surface vehicle. *Oceanography* 30, 104–109. doi: 10.5670/oceanog.2017.228

Thwaites, F. T., Krishfield, R., Timmermans, M. L., Toole, J. M., and Williams, A. J. (2011). "Flux measurements from an ice-tethered profiler: first look," in *Proceedings of the Oceans 2011, 6-9 June 2011, IEEE/OES*, (Santander: IEEE), 6.

Tian-Kunze, X., Kaleschke, L., Maaß, N., Mäkynen, M., Serra, N., Drusch, M., et al. (2014). SMOS-derived thin sea ice thickness: algorithm baseline, product specifications and initial verification. *Cryosphere* 8, 997–1018. doi: 10.5194/tc-8-997-2014

Tietsche, S., Notz, D., Jungclaus, J. H., and Marotzke, J. (2013). Assimilation of sea-ice concentration in a global climate model-physical and statistical aspects. *Ocean Sci.* 9, 19–36. doi: 10.5194/os-9-19-2013

Timmermans, M.-L., Cole, S. T., and Toole, J. M. (2011). Horizontal density structure and restratification of the Arctic ocean surface layer. *J. Phys. Oceanogr.* 42, 659–668. doi: 10.1175/JPO-D-11-0125.1
Timmermans, M. L., Krishfield, R., Laney, S., and Toole, J. (2010). Ice-Tethered Profiler measurements of dissolved oxygen under permanent ice cover in the Arctic Ocean. *J. Atmos. Ocean. Technol.* 27, 1936–1949. doi: 10.1175/2010JTECHO772.1

Timmermans, M. L., Proshutinsky, A., Golubeva, E., Jackson, J. M., Krishfield, R., McCall, M., et al. (2014). Mechanisms of pacific summer water variability in the Arctic's central Canada basin. *J. Geophys. Res.* 119, 7523–7548. doi: 10.1002/2014JC010273

Timmermans, M. L., Toole, J., Krishfield, R., and Winsor, P. (2008). Ice-Tethered Profiler observations of the double-diffusive staircase in the Canada Basin thermocline. *J. Geophys. Res.* 113:C00A02. doi: 10.1029/2008JC004829

Toole, J. M., Krishfield, R. A., Timmermans, M.-L., and Proshutinsky, A. (2011). The ice-tethered profiler: Argo of the arctic. *Oceanography* 24, 126–135. doi: 10.5670/oceanog.2011.64

Trenberth, K. E., and Fasullo, J. T. (2010). Simulation of present day and twenty first century energy budgets of the Southern Oceans. *J. Clim.* 23, 440–454. doi: 10.1175/2009jcli3152.1

Tschudi, M., Meier, W. N., Stewart, J. S., Fowler, C., and Maslanik, J. (2019). *Polar Pathfinder Daily 25 km EASE-Grid Sea Ice Motion Vectors, Version 4*. Boulder, CO: NASA National Snow and Ice Data Center Distributed Active Archive Center.

Uotila, P., Goosse, H., Haines, K., Chevallier, M., Barthélémy, A., Bricaud, C., et al. (2018). An assessment of ten ocean reanalyses in the polar regions. *Clim. Dyn.* 52, 1613–1650. doi: 10.1007/s00382-018-4242-z

Vecchi, G. A., Delworth, T., Gudgel, R., Kapnick, S., Rosati, A., Wittenberg, A. T., et al. (2014). On the seasonal forecasting of regional tropical cyclone activity. *J. Clim.* 27, 7994–8016.

Vincent, R. F., Marsden, R. F., Minnett, P. J., Creber, K. A. M., and Buckley, J. R. (2008). Arctic waters and marginal ice zones: a composite Arctic sea surface temperature algorithm using satellite thermal data. *J. Geophys. Res.* 113:C04021. doi: 10.1029/2007JC004353

Wood, K. R., Jayne, S. R., Mordy, C. W., Bond, N., Overland, J. E., Ladd, C., et al. (2018). Results of the first arctic heat open science experiment. *Bull. Am. Meteorol. Soc.* 99, 513–520. doi: 10.1093/aob/mcw237

Xie, J., Bertino, L., Counillon, F., Lisaeter, K. A., and Sakov, P. (2017). Quality assessment of the TOPAZ4 reanalysis in the Arctic over the period 1991-2013. *Ocean Sci.* 13, 123–144. doi: 10.5194/os-13-123-2017

Xie, J., Counillon, F., and Bertino, L. (2018). Impact of assimilating a merged sea-ice thickness from CryoSat-2 and SMOS in the Arctic

reanalysis. *Cryosphere* 12, 3671–3691. doi: 10.5194/tc-12-3671-2018

Xie, J., Counillon, F., Bertino, L., Tian-Kunze, X., and Kaleschke, L. (2016). Benefits of assimilating thin sea ice thickness from SMOS into the TOPAZ system. *Cryosphere* 10, 2745–2761. doi: 10.5194/tc-10-2745-2016

Yang, Q., Losa, S. N., Losch, M., Liu, J., Zhang, Z., Nerger, L., et al. (2015). Assimilating summer sea-ice concentration into a coupled ice-ocean model using a LSEIK filter. *Ann. Glaciol.* 56, 38–44. doi: 10.3189/2015aog69a740

Yang, Q., Losa, S. N., Losch, M., Tian-Kunze, X., Nerger, L., Liu, J., et al. (2014). Assimilating SMOS sea ice thickness into a coupled ice-ocean model using a local SEIK filter. *J. Geophys. Res. Oceans* 119, 6680–6692. doi: 10.1002/2014jc009963

Zhang, Y., Maslowski, W., and Semtner, A. J. (1999). Impact of mesoscale ocean currents on sea ice in high-resolution Arctic ice and ocean simulations. *J. Geophys. Res.* 104, 18409–18429. doi: 10.1029/1999jc900158

Zhao, M., Timmermans, M. L., Cole, S., Krishfield, R., Proshutinsky, A., and Toole, J. (2014). Characterizing the eddy field in the Arctic Ocean halocline. *J. Geophys. Res.* 119, 8800–8817. doi: 10.1002/2014JC010488

Zhao, M., Timmermans, M. L., Cole, S., Krishfield, R., and Toole, J. (2016). Evolution of the eddy field in the Arctic ocean's Canada basin, 2005 - 2015. *Geophys. Res. Lett.* 43, 8106–8114. doi: 10.1002/2016FL069671

Zuo, G., Dou, Y., and Lei, R. (2018). Discrimination algorithm and procedure of snow depth and sea ice thickness determination using measurements of the vertical ice temperature profile by the ice-tethered buoys. *Sensors* 18:4162. doi: 10.3390/s18124162

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2019 Smith, Allard, Babin, Bertino, Chevallier, Corlett, Crout, Davidson, Delille, Gille, Hebert, Hyder, Intrieri, Lagunas, Larnicol, Kaminski, Kater, Kauker, Marec, Mazloff, Metzger, Mordy, O'Carroll, Olsen, Phelps, Posey, Prandi, Rehm, Reid, Rigor, Sandven, Shupe, Swart, Smedstad, Solomon, Storto, Thibaut, Toole, Wood, Xie, Yang and the WWRP PPP Steering Group. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.