We identified the nature and driving mechanisms of subinertial variability (variability at a time scale of several days) in four fjords in Southeast Greenland, in three high-resolution numerical simulations. We find two dominant frequency ranges in along-fjord velocity, volume transport of Atlantic Water, and along-fjord heat transport: one around 2–4 days and one around 10 days. The higher frequency is most prominent in the two smaller fjords (Sermilik Fjord and Kangerdlugssuaq Fjord), while the lower frequency peak dominates in the larger fjords (Scoresby Sund and King Oscar Fjord). The cross-fjord structure of variability patterns is determined by the fjord's dynamic width, while the vertical structure is determined by the stratification in the fjord. The dominant frequency range is a function of stratification and fjord length, through the travel time of resonant internal Kelvin waves. We find that the subinertial variability is the imprint of Coastal Trapped Waves, which manifest as Rossby-type waves on the continental shelf and as internal Kelvin-type waves inside the fjords. Between 50% and 80% of the variability in the fjord is directly forced by Coastal Trapped Waves propagating in from the shelf, with an additional role played by alongshore wind forcing on the shelf.

Computational oceanography is the study of ocean phenomena by numerical simulation, especially dynamical and physical phenomena. Progress in information technology has driven exponential growth in the number of global ocean observations and the fidelity of numerical simulations of the ocean in the past few decades. The growth has been exponentially faster for ocean simulations, however. We argue that this faster growth is shifting the importance of field measurements and numerical simulations for oceanographic research. It is leading to the maturation of computational oceanography as a branch of marine science on par with observational oceanography. One implication is that ultraresolved ocean simulations are only loosely constrained by observations. Another implication is that barriers to analyzing the output of such simulations should be removed. Although some specific limits and challenges exist, many opportunities are identified for the future of computational oceanography. Most important is the prospect of hybrid computational and observational approaches to advance understanding of the ocean.

Ocean currents along the Southeast Greenland Coast play an important role in North Atlantic circulation and the global climate system. They carry dense water over the Denmark Strait sill, fresh water from the Arctic and the Greenland Ice Sheet into the subpolar ocean, and warm Atlantic water into Greenland’s fjords, where it can interact with outlet glaciers. Observational evidence from the OSNAP array and other mooring records shows that the circulation in this region displays substantial subinertial variability, typically with periods of several days. For the dense water flowing over the Denmark Strait sill, this variability augments the time-mean transport; on the shelf, the variability is large enough to occasionally reverse the mean transport direction of the coastal current, highlighting the importance of characterizing this variability when interpreting synoptic surveys. In this study, we used the output of a high-resolution realistic simulation to diagnose and characterize subinertial variability in sea surface height and velocity along the coast. The results show that the subinertial signals on the shelf and along the shelf break are coherent over hundreds of kilometers, and consistent with Coastal Trapped Waves in two subinertial frequency bands—at periods of 1–3 days and 5–18 days—portraying a combination of Mode I and higher modes waves. Furthermore, we find that northeasterly barrier winds may trigger the 5–18 day shelf waves, whereas the 1–3 day variability is linked to high wind speeds over Sermilik Deep.

Ocean currents along the Southeast Greenland Coast play an important role in the climate system. They carry dense water over the Denmark Strait sill, fresh water from the Arctic and the Greenland Ice Sheet into the subpolar ocean, and warm Atlantic water into Greenland’s fjords, where it can interact with outlet glaciers. Observational evidence from moorings shows that the circulation in this region displays substantial subinertial variability (typically with periods of several days). For the dense water flowing over the Denmark Strait sill, this variability augments the time-mean transport. It has been suggested that the subinertial variability found in observations is associated with Coastal Trapped Waves, whose properties depend on bathymetry, stratification, and the mean flow. Here, we use the output of a high-resolution realistic simulation to diagnose and characterize subinertial variability in sea surface height and velocity along the coast. The results show that the subinertial signals are coherent over hundreds of kilometers along the shelf. We find Coastal Trapped Waves on the shelf and along the shelf break in two subinertial frequency bands—at periods of 1–3 days and 5–18 days—that are consistent with a combination of Mode I waves and higher modes. Furthermore, we find that northeasterly barrier winds may trigger the 5–18 day shelf waves, whereas the 1–3 day variability is linked to high wind speeds over Sermilik Deep.

The Denmark Strait Overflow (DSO) is an important contributor to the lower limb of the Atlantic meridional overturning circulation (AMOC). Determining DSO formation and its pathways is not only important for local oceanography but also critical to estimating the state and variability of the AMOC. Despite prior attempts to understand the DSO sources, its upstream pathways and circulation remain uncertain due to short-term (3–5 days) variability. This makes it challenging to study the DSO from observations. Given this complexity, this study maps the upstream pathways and along-pathway changes in its water properties, using Lagrangian backtracking of the DSO sources in a realistic numerical ocean simulation. The Lagrangian pathways confirm that several branches contribute to the DSO from the north such as the East Greenland Current (EGC), the separated EGC (sEGC), and the North Icelandic Jet (NIJ). Moreover, the model results reveal additional pathways from south of Iceland, which supplied over 16% of the DSO annually and over 25% of the DSO during winter of 2008, when the NAO index was positive. The southern contribution is about 34% by the end of March. The southern pathways mark a more direct route from the near-surface subpolar North Atlantic to the North Atlantic Deep Water (NADW), and needs to be explored further, with in situ observations.

Mesoscale features present at the Denmark Strait sill regularly enhance the volume transport of the Denmark Strait overflow (DSO). They are important for the Atlantic Meridional Overturning Circulation and ultimately, for the global climate system. Using a realistic numerical model, we find new evidence of the causal relationship between overflow surges (i.e., mesoscale features associated with high-transport events) and DSO cyclones observed downstream. Most of the cyclones form at the Denmark Strait sill during overflow surges and, because of potential vorticity conservation and stretching of the water column, grow as they move equatorward. A fraction of the cyclones form downstream of the sill, when anticyclonic vortices formed during high-transport events start collapsing. Regardless of their formation mechanism, DSO cyclones weaken starting roughly 150 km downstream of the sill, and potential vorticity is only materially conserved during the growth phase.

Multimodel Arctic Ocean “climate response function” experiments are analyzed in order to explore the effects of anomalous wind forcing over the Greenland Sea (GS) on poleward ocean heat transport, Atlantic Water (AW) pathways, and the extent of Arctic sea ice. Particular emphasis is placed on the sensitivity of the AW circulation to anomalously strong or weak GS winds in relation to natural variability, the latter manifested as part of the North Atlantic Oscillation. We find that anomalously strong (weak) GS wind forcing, comparable in strength to a strong positive (negative) North Atlantic Oscillation index, results in an intensification (weakening) of the poleward AW flow, extending from south of the North Atlantic Subpolar Gyre, through the Nordic Seas, and all the way into the Canadian Basin. Reconstructions made utilizing the calculated climate response functions explain ∼50% of the simulated AW flow variance; this is the proportion of variability that can be explained by GS wind forcing. In the Barents and Kara Seas, there is a clear relationship between the wind-driven anomalous AW inflow and the sea ice extent. Most of the anomalous AW heat is lost to the atmosphere, and loss of sea ice in the Barents Sea results in even more heat loss to the atmosphere, and thus effective ocean cooling. Release of passive tracers in a subset of the suite of models reveals differences in circulation patterns and shows that the flow of AW in the Arctic Ocean is highly dependent on the wind stress in the Nordic Seas.

Initial results are presented from a yearlong, high-resolution (~2 km) numerical simulation covering the east Greenland shelf and the Iceland and Irminger Seas. The model hydrography and circulation in the vicinity of Denmark Strait show good agreement with available observational datasets. This study focuses on the variability of the Denmark Strait overflow (DSO) by detecting and characterizing boluses and pulses, which are the two dominant mesoscale features in the strait. The authors estimate that the yearly mean southward volume flux of the DSO is about 30% greater in the presence of boluses and pulses. On average, boluses (pulses) are 57.1 (27.5) h long, occur every 3.2 (5.5) days, and are more frequent during the summer (winter). Boluses (pulses) increase (decrease) the overflow cross-sectional area, and temperatures around the overflow interface are colder (warmer) by about 2.6°C (1.8°C). The lateral extent of the boluses is much greater than that of the pulses. In both cases the along-strait equatorward flow of dense water is enhanced but more so for pulses. The sea surface height (SSH) rises by 4–10 cm during boluses and by up to 5 cm during pulses. The SSH anomaly contours form a bowl (dome) during boluses (pulses), and the two features cross the strait with a slightly different orientation. The cross streamflow changes direction; boluses (pulses) are associated with veering (backing) of the horizontal current. The model indicates that boluses and pulses play a major role in controlling the variability of the DSO transport into the Irminger Sea.

Seasonal variability in pathways of warm-water masses toward the Kangerdlugssuaq Fjord (KF)–Glacier (KG) system, southeast Greenland, is investigated by backtracking Lagrangian particles seeded at the fjord mouth in a high-resolution regional ocean model simulation in the ice-free and the ice-covered seasons. The waters at KF are a mixture of Atlantic-origin water advected from the Irminger Basin [Faxaflói (FF)], the deep waters from the Denmark Strait, and the waters from the Arctic Ocean, both represented by the Kögur section (KO). Below 200-m depth, the warm water is a mixture of FF and KO water masses and is warmer in winter than in summer. The authors find that seasonal differences in pathways double the fraction of FF particles in winter, causing the seasonal warming and salinification. Seasonal temperature variations at the upstream sections (FF and KO) have a negligible impact on temperature variations near the fjord. Successful monitoring of heat flux to the fjord therefore needs to take place close to the fjord and cannot be inferred from upstream conditions.

As numerical ocean simulations become more realistic, analysis of their output is increasingly time consuming. As part of a larger effort to make ocean model output more accessible for analysis, we have developed a fast particle-tracking algorithm for exploring the kinematics of the simulated flow. The algorithm is independent of operating system and fully vectorized and parallelized. Furthermore, it enables sliding of particles along solid boundaries according to the 3D along-boundary flow field. The new algorithm can easily simulate several million particle trajectories, thus opening the way for Lagrangian analysis of large-scale and multi-time scale oceanic phenomena.