1. Introduction and science background
A large portion of the North Atlantic Deep Water (NADW) can be traced to overflows from the Nordic Seas across the Greenland–Iceland–Scotland Ridge (GISR; Fig. 1). The overflow transport of the dense water mass, defined as those with a potential density higher than 27.8 kg m−3, is about 6 Sverdrups (Sv; 1 Sv ≡ 106 m3 s−1), which includes about 3 Sv through the Denmark Strait (DS; sill depth of 620 m), 2 Sv via the Faroe Bank Channel (FBC; sill depth of 840 m), and about 1 Sv over the Iceland–Faroe Ridge (IFR; maximum sill depth of 420 m) (Hansen et al. 2008). Downstream of the sills, the overflowed water mass mixes vigorously with the ambient Atlantic Ocean water and forms the core of the NADW that flows southward as the lower limb of the Atlantic meridional overturning circulation (AMOC; Price and Baringer 1994). Both observations (Macrander et al. 2005, 2007; Eldevik et al. 2009) and numerical simulations (e.g., Köhl et al. 2007; Köhl 2010; Serra et al. 2010) indicate that the Denmark Strait overflow (DSO) varies considerably on interannual time scales. Such variations result from atmospheric forcing and internal variability along the upstream pathways.
A marginal sea overflow involves processes on both sides of the sill. This study focuses on the upstream pathways to the north of the GISR, and in particular those that are related to the DSO. In previous studies, the DSO pathways were identified mainly by comparing hydrochemical characteristics of water masses in the Nordic and Arctic basins. Direct current measurements have been scarce. Swift et al. (1980) suggested that the water mass formed in the Iceland Sea supplies the DSO. Smethie and Swift (1989) proposed that the denser part of the DSO water originates from the Greenland Sea. Aagaard et al. (1991) and Buch et al. (1996) argued that water masses from the Arctic contribute to the overflow. Mauritzen (1996a,b) proposed that the Atlantic Water that has been transformed in the Nordic and Arctic basins, the so-called returned Atlantic Water, is the main component for the Denmark Strait overflow. Rudels et al. (2002) analyzed observational data and concluded that the East Greenland Current (EGC) is the main pathway along which various components of the DSO water mass are transported (Fig. 1 for the schematics). Their conclusion is supported by analyses of hydrochemical data (Tanhua et al. 2005; Jeansson et al. 2008). Even though the EGC has been considered the main pathway in some previous studies, it has also been recognized that the water masses in the DSO consist of source waters from various Arctic and Nordic basins and that its composition has evolved with time (Rudels et al. 2003; Tanhua et al. 2005).
The prevailing view of the EGC being the dominant pathway for the DSO was questioned when measurements from moorings deployed for 3 yr (1988–91) off the northern Icelandic shelf and slope showed a persistent southwestward flow at about 500-m depth, with a mean speed of about 10 cm s−1, toward the Denmark Strait (Jónsson 1999). Jónsson and Valdimarsson (2004) confirmed the existence of this current based on shipboard observations made in November 2001 and 2002. The existence of the Icelandic branch of the DSO was further confirmed by extensive surveys along the northern Icelandic shelf and slope in 2008 and 2009 (Våge et al. 2011). They observed that this branch, which they called the north Icelandic jet (NIJ), is centered on the 600-m isobath (approximately the sill depth in the Denmark Strait) and that its origin can be traced to the region northeast of Iceland (Fig. 1). Våge et al. (2013) analyzed a new set of observations made in August 2011 and 2012 and showed that about 70% of the DSO water approaches the sill along the Iceland continental slope, and it is supplied by two separate currents, namely, the NIJ and a separated branch of the EGC. (The latter has been diverted eastward from the Greenland continental slope upstream of the sill.) Våge et al. (2013) suggested that eddy fluxes associated with a baroclinically unstable EGC are responsible for the EGC separation. Such a separation, as noted by Våge et al. (2013), also occurs in rotating hydraulics models (e.g., Helfrich and Pratt 2003; Yang and Price 2000) when an outflowing current approaches a shallow sill. The bifurcation of a southward western boundary current (WBC) in a hydraulics model is mostly determined by the layer thickness upstream of the sill. It occurs when the layer thickness of the source water h2 is affected by the shoaling bathymetry and the flow along Greenland’s coast can no longer maintain its course toward the Denmark Strait along f/h2 isolines (where f is the Coriolis frequency). If the WBC is flowing along isobaths that are deeper than the sill depth, it tends to detach from the western boundary and then approach the sill on the eastern side. The topographic slope leading to the sill makes the eastern boundary behave like a dynamical western boundary. If the approaching WBC is along isobaths shallower than the sill depth, it usually continues southward along the western boundary toward the sill. In the EGC case, it is expected, in a hydraulics model, that the portion of the flow on the inner slope along isobaths that are shallower than the sill depth would continue to approach the sill as a WBC, while the flow along deeper isobaths detach from the western boundary. It is possible that the two mechanisms (eddy and topographic beta) are linked. Våge et al. (2013) demonstrated convincingly that the NIJ is an independent pathway that coexists with the separated EGC and that the NIJ accounts for roughly half of the DSO transport. In a study using an idealized model, Våge et al. (2011) suggest that the NIJ could be fed entirely by water that has been cooled in the Iceland Sea, and that has originally made its way into the Nordic Seas via the north Icelandic Irminger Current, which flows at the surface northwestward through the DS and around Iceland. The numerical model domain consisted of only a small portion of the Nordic Seas.
In this study, we use a two-layer marginal sea overflow model (Yang and Pratt 2013) to investigate some dynamical constraints on the DSO pathways with the main focus on the existence of the NIJ. We do not address where the DSO source water is formed but rather some basic dynamical processes and balances that relate to the DSO pathways regardless of the origin of the source water. The two-layer framework is idealized, but it accounts for the basic exchange properties between the Atlantic Ocean and the Nordic Seas across the GISR (e.g., Dickson et al. 2008; Pratt and Whitehead 2008), and it provides a good setting for intuition building. The numerical model is described in section 2, and this is followed in section 3 by discussion of some simulations that lack wind forcing and produce an NIJ-like current. In section 4, we address the apparently robust nature of this result using circulation integrals. This discussion will emphasize the importance of the “island” geometry of Iceland. Wind stress forcing leads to a more realistic lower-layer circulation while preserving the NIJ, and this is described in section 5. A summary is given in section 6.
2. A two-layer marginal sea overflow model
We use a primitive equation, two-layer model (Yang and Pratt 2013) in this study. Both layers are active and so the model includes the barotropic and the first baroclinic modes. The thickness of either layer is allowed to become zero (i.e., outcropping of the lower or grounding of the upper layer). A marginal sea overflow has been considered in many studies as a two-layer exchange over a ridge between two basins. Thus, two-layer ocean models have been one of the most commonly used tools in such studies [see Pratt and Whitehead (2008) for a comprehensive review]. It is assumed that the upper layer is filled by buoyant water masses of either high temperature or low salinity, while the lower layer contains dense water masses that have been formed in the marginal sea by surface air–sea fluxes. In the Nordic Seas, a potential density of 27.8 kg m−3 is often used to separate the upper and lower layers.
While the deep Nordic Seas are filled with a rather uniformly cold, dense, water mass, the upper layer in the Nordic Seas contains a wide range of water masses, including low salinity and cold Arctic Ocean water along Greenland’s coast and warm and saline Atlantic Water off the coast of Norway. These waters are modified by air–sea interaction as they circulate. In our two-layer model, horizontal changes in water properties are contained only in variability of the upper-layer thickness and the ability of the lower layer to outcrop, so features such as the gradual densification of the Atlantic Water as it circulates around the basins cannot be captured in a realistic way. Nevertheless, we anticipate that the model will provide fundamental information and insight concerning preferred deep pathways toward the Denmark Strait. The simplicity and transparency in dynamics are its main advantage over OGCMs. The potential density σθ = 27.8 kg m−3 is often used to separate the inflowing and outflowing layers, and we think of this surface as more or less coinciding with our two-layer interface.
The model uses both idealized and realistic bathymetry (Fig. 1). In all experiments with a realistic topography, the model extends from 55° to 80°N and has a resolution of
3. The DSO pathways
Our objective is to identify factors that influence the deep pathways, including the NIJ, from the upstream source regions to the DS. For now we will refer to any modeled deep westward flow along the north Icelandic coast as the NIJ, with the caveat that particular features such as volume flux and linkage with upstream currents may differ from observations and may vary between numerical experiments. The DSO is a topographically restricted and hydraulically controlled flow (Käse and Oschlies 2000; Käse et al. 2003; Wilkenskjeld and Quadfasel 2005). The Nordic Seas overflow, including the DSO, is primarily driven by a buoyancy flux (Hansen et al. 2008). The wind stress forcing modulates overflow variability (e.g., Biastoch et al. 2003; Köhl et al. 2007; Nilsen et al. 2003; Serra et al. 2010). In this section, we will discuss first the DSO pathways in circulations that are forced by water-mass sources and sinks, with no wind.
Previous studies indicate that there are likely multiple origins of the DSO source water, and we will consider the outcome of varying the origin later in this paper. We first wish to establish a benchmark [experiment 1 (EXP-1)] in which we determine the draining pathways that arise when the dense-water source is turned off [we = 0 in (1)] and the outflow occurs as a result of the sudden release of a preexisting volume of dense water. In the experiment setup, the dense water initially fills the basin to the north of the GISR to the depth level that is 50 m below the sea surface (Fig. 2). The lighter water in the upper layer fills the rest of the model domain, that is, the entire basin to the south of the ridge and the upper 50 m in the Nordic Seas. When this volume is released, the initially piled-up dense water floods across the ridge and descends into the deep ocean. This flow adjustment process is similar to the classical lock exchange experiment in a rotating system (Pratt and Whitehead 2008). Similar numerical experiments have been used to examine the Nordic Seas overflow pathways (Käse et al. 2009) and to assess the effective capacity of the dense-water reservoir (Yang and Pratt 2013).
Figure 3a shows the sea surface height (SSH) and upper-layer velocity at the end of a 1-yr model run. The lighter water in the upper layer flows northward into the Nordic Seas through passages between Iceland and Scotland and along the eastern side of the Denmark Strait. A southward East Greenland Current flows on the western side of the Denmark Strait. A prominent feature in the surface circulation in the Nordic Seas is a basinwide cyclonic rim current. The flow is very weak in the interior. The dense water masses in the deep basins are stored within closed geostrophic contours and thus somewhat buffered from the boundary currents, a feature explored in more depth by Yang and Pratt (2013). The circulation in the lower layer (Fig. 3b) is dominated by an anticyclonic boundary current in the Nordic Seas. The circulation in the deep Greenland Sea is very weak, and there is no deep expression of the EGC. The boundary current draws water mass mainly from the Lofoten basin, which seems to be consistent with the study of Isachsen et al. (2007). The dense water moves southward along the continental slope off Norway and feeds the Faroe Bank Channel overflow. A branch of this deep current flows toward the west along the 600–800-m isobaths north of the GISR and exits through Denmark Strait. In the vicinity of Iceland, this westward flow resembles the NIJ. But it should be noted that the westward jet in EXP-1 extends far to the east of Iceland, whereas the observations of Våge et al. (2013) suggest a termination point along the north Icelandic coast, a feature supported by four independent surveys. Our model NIJ becomes more realistic when wind stress is added (to be discussed in section 5).
The dam break experiment shows that a westward boundary current along the north Icelandic coast is the dominant pathway for the DSO in a drainage overflow without the renewal of deep water. In reality, the overflow source water is replenished by winter convection, so we next consider the simplest of cases with a source of overflow water. In EXP-2, a uniform diapycnal flux, that is, we = constant > 0 (positive for a downward flux from upper to lower layer), is applied over the entire Nordic Seas to the north of the GISR (within the red contour line shown in Fig. 4a). The total diapycnal flux is specified to be 6 Sv initially (the same for all experiments that uses we). (To close the overturning circulation, a region of negative we is imposed at southwestern corner of the model domain south of the GISR.) The wind stress is turned off as in the previous experiment.
The downwelling results in outcropping of the lower layer in the interior Nordic Seas. The area within the black contour in Fig. 4a indicates the area of outcropping at the end of the 10-yr simulation. The downwelling we is turned off whenever the lower-layer outcrops. So the net downwelling is less than the initial 6 Sv. At the end of the tenth year in EXP-2, the downwelling occurs only in the area between the red and black lines shown in Fig. 4a. Figure 4b shows that the overall rim circulation in the lower layer is again anticyclonic. The southward eastern boundary current supplies dense water to both FBC and the DS. The overall circulation pattern is similar to that in the dam break experiment (Fig. 3b). A noticeable difference is the existence of a northward western boundary current in the lower layer of the Greenland basin. This deep current is opposite in flow direction to the EGC in the upper layer and is associated with an anticyclonic Greenland Sea gyre that is forced by vortex squashing induced by the downwelling we. The lower-layer circulations in the other two major deep basins (Norwegian and Lofoten) are also anticyclonic. In areas between the black and red lines in Fig. 4a, the vortex squashing forces flow across f/h contours from high to low PV. In the area where the lower-layer outcrops (within the thick black contour in Fig. 4a), we is zero and the flow is mainly along f/h contours. The flows around the GISR are similar to that in the dam break experiment. The NIJ in the lower layer in this simulation is robust and flows westward along the 500–800-m isobaths toward the Denmark Strait. In fact, the DSO transport in this experiment is supplied almost completely by the NIJ.
In reality the water-mass transformation is certainly not uniformly distributed in the Nordic Seas. The DSO water is a composition of water masses from multiple origins according to analyses of hydrochemical data (Tanhua et al. 2005). We would like to examine whether the NIJ’s existence depends on where and how the DSO source water is formed. Intuitively, one would expect that an overflow pathway is a direct conduit from the source to the sill and thus is tied to the source’s location. For instance, the EGC would be an expected pathway if the source of the DSO water is in the Greenland Sea. But as we will show next, the NIJ’s existence is rather insensitive to where the DSO source water is formed.
In the third experiment (EXP-3), a localized water-mass flux from the upper to lower layer is placed on the continental slope off the east coast of Greenland (the location of downwelling is indicated by the purple oval in Fig. 5a). As in previous experiments, the wind stress is turned off. The velocity in the lower layer is shown in Fig. 5a. The source water, which is injected into the lower layer on the continental slope of the Greenland Sea, forms a southward western boundary current within the Greenland Sea. But instead of continuing southward to the Denmark Strait, it recirculates within the Greenland Sea and flows across the Mohn Ridge into the Lofoten and Norwegian basins. There are strong anticyclonic circulations in these two basins. The flow along the Jan Mayen Ridge is northward as a part of the anticyclonic gyre in the Norwegian basin. There are flows from the northern Greenland Sea southward along the continental slope off Spitsbergen and Norway into the Lofoten and Norwegian basins. This eastern boundary current continues southward and splits, with one branch feeding the overflow through the Faroe Bank Channel and the second feeding the westward NIJ along the 600–800-m isobaths to the north of Iceland–Faroe Ridge. The NIJ is the dominant source for the Denmark Strait overflow. The East Greenland Current just north of the Denmark Strait is weak even though the source of the low-layer water is placed in the Greenland Sea. Some features of EXP-3, including the sense of deep circulation in the three major basins, differ from observations, and are the same as in the previous two experiments without wind stress forcing. The important point is that, despite these differences, the NIJ is present.
We have conducted several additional source/sink-driven experiments by moving the water-mass source to different locations or by specifying the source as a southward inflow boundary condition across the Fram Strait from the Arctic basin. The NIJ always emerges as the dominant pathway for the Denmark Strait overflow. In fact, we have not been able to simulate a Nordic Seas overflow that does not have the NIJ as its main pathway for the DSO in all of our source- and sink-driven experiments. The NIJ is apparently a robust and permanent feature, at least in model runs that lack wind forcing. We now show that its existence in the model is consistent with some constraints, one of them quite powerful, involving potential vorticity dynamics and circulation.
4. Plausible mechanisms for the NIJ
In this section, we will discuss two integral constraints that identify conditions that are favorable to the existence of the NIJ-like current. The first is a circulation integral that extends around a contour C that rims the Nordic Seas (Fig. 6). The basin may be closed to the north or may contain an opening that acts as a potential source of overflow water from the Arctic. In terms of dynamics, a key question that will arise is whether the overflow source water originates from the deeper portions of the basin, and therefore has relatively low potential vorticity (PV), or whether it originates from shallower coastal areas. The tendencies predicted by the first constraint are thus sensitive to the location of the diapycnal velocity we and especially to whether it is concentrated in the interior of the basins or near the boundaries. We will consider both situations. The second constraint is stronger and derives from the circulation around the island contour CI that circles Iceland and the Iceland–Faroe Ridge (Fig. 6). This integration constraint is not sensitive to the location of the diapycnal velocity we in the Nordic Seas. The details of these two integration constraints are discussed next.
a. A basin circulation integral constraint
In a source- and sink-driven flow, the circulation around any closed contour C is constrained by the vorticity advection across C as well as the tangential stresses acting on C. Pratt (1997), Yang and Price (2000, 2007), and Helfrich and Pratt (2003) used circulation integrals to infer flow directions around closed geostrophic contours or lateral boundaries in semienclosed marginal seas. Yang (2005) and Karcher et al. (2007) applied this constraint to explain that the circulation in the Atlantic Water layer in the Arctic is cyclonic in part because of a net positive transport of PV from sub-Arctic basins. These integral constraints can provide information about the average direction of tangential flow around C, but they do not necessarily determine the direction of flow along any a particular segments of C. The constraints tend to be simplest when applied along slippery, vertical walls since fluxes across C are then zero. For the more realistic case of sloping topography, an advantageous choice of C is a closed f/h2 contour (provided that the contour does not wander out of the domain of interest), with h2 chosen to correspond to the depth at which knowledge of the direction of circulation is desired.
The presence of average anticyclonic flow around the 620-m isobaths of the Nordic Seas does not in itself imply the presence of a westward NIJ. However, the presence of the East Greenland Current, which flows in a cyclonic sense, does strengthen the requirement for anticyclonic flow elsewhere around C in order to satisfy the integral constraint. This would strengthen the argument for westward flow along the north Icelandic slope, but this is as far as this argument can be pushed. It should also be recognized that the sinking that feeds the overflows could occur dominantly outside of C, a situation we will investigate in section 4b.
Support for the tendencies predicted by (5) can be gained from a pair of idealized experiments with simplified geometry (Fig. 7). In both experiments we use a model domain of 500 km × 800 km with a depth of 2500 m in the deep basin and of 1250 m at the sill (Fig. 7a). There is just one draining strait in either case, so we set QFBC = 0 in (4) and (5). We use a reduced-gravity model version of (1) by setting the upper-layer velocity and surface pressure gradient to be zero. The model uses Cartesian coordinates and has a resolution of 5 km. A water-mass source for the lower layer is introduced either as downwelling from the upper to lower layer in the northwestern corner (Fig. 7b) or as an inflow through the northern boundary (Fig. 7c). An outflow is specified in the southern boundary so that the mass flux into the lower layer is balanced. We are primarily interested in upstream pathways connecting the source to the DSO. In both forcing types, the source water is introduced in the northern basin and needs to overflow a shallower strait to the southern basin.
In the first idealized experiment (IDL-1), the water-mass source in the northern basin is completely due to a specified downwelling (we > 0) in (1). Thus, QFS = 0 in (4) and the overall rim circulation about the basin is predicted to be anticyclonic. As shown in Fig. 7b, the downwelling forces a strong local recirculation, and the boundary current around the northern basin is anticyclonic. The water mass approaches the shallow strait as an NIJ-like current along the shelf and slope of the model’s “Iceland.” These features are qualitatively similar to those source- and sink-driven experiments that use realistic bathymetry (Figs. 4b, 5a).
The next idealized experiment (IDL-2) explores a case in which dense water is fed into the interior through a passage that is deeper than the draining strait. In particular, the water-mass is specified as an inflow through a northern strait, analogous to the Fram Strait. Equation (5) may now be used, with QFBC = 0 and HFS < Hint, so that anticyclonic circulation around C is again predicted. Figure 7c shows that the inflow forms a cyclonic current along the northern and western boundaries. It is an analog of the East Greenland Current. The southward western boundary current, instead of continuing southward toward the sill, turns eastward, makes a U-turn, and approaches the strait as a westward coastal current along the model’s Iceland toward the connecting strait. Since there is little flow along the northern and eastern boundaries, the anticyclonic NIJ current must be sufficiently strong to overcome the cyclonic contribution from the model East Greenland Current. The flow pattern is essentially the same as that shown by Helfrich and Pratt (2003, their Fig. 7a) and by Yang and Price (2000, their Fig. 8a). It is also similar to that in EXP-3 with realistic bathymetry and a localized source on the Greenland slope (Fig. 5a). The cumulative integral of bottom friction (dashed line in Fig. 8) shows little contribution from the eastern boundary, a negative contribution from the northern and western boundaries, and a stronger positive contribution from the southern boundary or “Iceland coast” (1350–1750 km). It is interesting to note that the detachment from the southward WBC to form a westward flow along the northern coast of Iceland is also qualitatively similar to the separated EGC pathway described by Våge et al. (2013). However, their proposed detachment mechanism involves changes in the wind field along the Greenland coast, whereas detachment in the present, wind-free simulation is due to the topographic effect discussed by Helfrich and Pratt (2003) and Yang and Price (2000).
The integral of the circulation, as shown in (5), can be reversed if the source water enters the northern basin with a high PV. This can be accomplished by simply making the sill depth at the northern channel, that is, HFS in (5), shallower than that at the connecting channel. In the next experiment, we repeated IDL-2 with a small modification in bathymetry. The sill depth is set at 750 m at the northern open boundary (Fig. 9a). This localized change in bathymetry makes the right-hand-side term of (5) to be negative, a reverse from that in IDL-2. The integral of the circulation around C is cyclonic. The NIJ-like current vanishes (Fig. 9b). In fact, the western boundary current, the model’s EGC, continues southward and approaches the sill as a western boundary current. No detachment occurs.
As an integral constraint, (3), or its approximation (4), governs the total circulation around C but does not dictate the direction of flow along any subsegment. In fact, the flow along the western boundary in the second idealized experiment (IDL-2; Fig. 7c) is actually southward or cyclonic. There are scenarios, including the one discussed below, in which simulated flows satisfy (4), but lack NIJ-like currents. The assumption that the source water is formed in a deep basin and thus has a low PV is rather ad hoc and overly constraining. It is possible that a portion of the DSO source water is formed on the shallow shelf along the boundary and outside the contour C. So the integral constraint (4) or (5) is suggestive but not definitive. A stronger constraint is discussed in section 4c below, but first we consider the case in which the dense outflow is fed from shallow, coastal areas.
b. Net sinking near the boundary
Modeling work and observations over the past two decades (e.g., Böning et al. 1996; Marotzke and Scott 1999; Spall 2003, 2004; Pedlosky 2003; Straneo 2006; Våge et al. 2011; Cenedese 2012) suggest that deep convection in the interior of a high-latitude basin produces little net sinking and that the latter tends to occur close to boundaries. A cartoon expressing some of the elements of this work (in particular, Spall 2003, 2004; Våge et al. 2011; Cenedese 2012) appears in Fig. 10. A warm surface boundary current enters the marginal sea along the right-hand coast. This current could be analogous the north Icelandic Irminger Current, in which case the right-hand coast would correspond to the north coast of Iceland. The current is baroclinically unstable and sheds warm eddies that disperse into the interior, where they are cooled. This interior region experiences thermal convection down to some level, but the spatially averaged vertical velocity is zero. So there is no net sinking, and thus the deep outflow is not fed from the interior. Instead there is a return flow of cooler eddies into the boundary current. Net sinking preferentially occurs in these models near the boundary itself, where the vorticity production by sinking-induced stretching is balanced by frictional effects. This mechanism also occurs in laminar models (e.g., Pedlosky 2003), but eddies are implicit in the parameterizations of friction and diffusion. The implication is that the cold outflow from the marginal sea is fed strictly by the boundary sinking, though this pathway has not been explicitly identified in any model. [Though interior convection is primarily an overturning process, a small net vertical velocity over a large area could also contribute to the outflow. Only one such study, that of Pratt and Spall (2008), considers a case where the marginal sea outflow is hydraulically controlled and partially blocked below sill depth, but the distribution of sinking and the source of the outflow was not determined.]
In light of the work just cited, we must also consider the possibility that the overflow is fed by sinking processes occurring near the boundary, perhaps entirely outside the contour C. We therefore explore a case in which we is finite only near the coast. In interpreting the results that follow, the reader may wish to keep in mind that the cross-interface velocity is associated with both sinking and water mass transformation in a layer model. The spatial separation between thermodynamic transformations and vorticity-producing net sinking that is possible in the continuously stratified models is difficult in the two-layer formulation. The dynamical and thermodynamical consequences of we occur together, and when we specify sinking near a coast, the implication is that the water mass transformation also occurs there.
Some of the elements of coastal sinking were present in EXP-2, where a uniform we is initially applied in the whole Nordic Seas. As the flow evolves, the interface rises in the interior and the lower layer eventually outcrops over the deeper portions of the basin. The downwelling we in the outcropped region (within the black contour in Fig. 4a) is then turned off and subsequently restricted to the region between the black line and the coast. Interpretation of the integral constraints (4) or (5) is tricky because we occurs both inside and outside C and thus the direction of Qint is uncertain.
Because of this ambiguity, we conducted an additional experiment in which we is finite only outside of C. In the fourth experiment (EXP-4), a uniform downwelling we is applied only in areas where the depth is shallower than the DS sill depth (620 m; areas within the purple contours in Fig. 5b). In the deeper region inside of C, the downwelling we is set to zero. As shown in Fig. 5b, the NIJ-like current remains a robust feature along the north coast of Iceland (Fig. 5b). The PV integral is clearly not helpful in this case, and we therefore seek a stronger constraint or mechanism that promotes the NIJ’s existence.
c. Mechanism two: An island circulation integral constraint
Here, we consider the circulation integral around the contour CI corresponding to a closed isobath lying slightly above the DS sill depth (620 m) and circling the Iceland–Faroe Ridge (Fig. 6). This contour exists only when the marginal sea has two outlets. It passes unimpeded through the deeper Faroe Bank Channel (sill depth 840 m) and lies well within the observed depth range of the NIJ.
We have tested these ideas by performing two additional idealized experiments. In both, we modify the Fig. 7 model bathymetry by adding an additional strait, analogous to the Faroe Bank Channel, so that there is a closed geostrophic contour between two outflow channels (Fig. 11). In the first of the new experiments (IDL-3; Figs. 11a,b), a water-mass source is introduced as an inflow through a northern inlet, and an outflow is placed on the southern boundary, just as that in IDL-2 (Fig. 7c). The water mass must exit through one or both straits from the northern to the southern basin, but the draining pathways are determined by the model dynamics.
Figure 11b shows the layer interface height and lower-layer velocity at the end of a 10-yr model run. In the upstream (northern) basin, the water mass enters the model domain and flows cyclonically along the northern and western boundary toward the western strait, the model’s Denmark Strait. As in the second idealized experiment (IDL-2; Fig. 7c), the southward western boundary current (model’s EGC) turns eastward into the interior instead of continuing southward toward the western strait. On the eastern side of the basin, there is a southward eastern boundary current that flows directly into the eastern strait or the model’s Faroe Bank Channel. As in Fig. 7c, the overflow through the western strait is fed largely by an NIJ-like current along the north coast of the model’s Iceland. Note the westward flow along the south coast of Iceland, originating from the model Faroe Bank Channel. The frictional stress associated with this flow is balanced in (9) by that due to the NIJ. In general, the flow is very similar to that in the IDL-2 (Fig. 7c). The model NIJ can also be motivated in terms of the basin integral (5), without using the island integral (9), so we next present an experiment designed to distinguish between the consequences of the two constraints.
In the fourth idealized experiment (IDL-4), a shallow ridge is inserted between the southern boundary and the midbasin island (Fig. 11c). The northern basin is unchanged so that the application of the PV integral about the northern basin is not affected. The insertion of this ridge, however, means that the island contour CI is no longer closed and so (9) is no longer valid. Figure 11d shows the model flow at the end of a 10-yr run. The circulation is similar to the previous run (Fig. 11b) except that the NIJ vanishes completely. The overflow is mainly through the eastern strait or the model’s FBC. The transport through the western strait (or the DS) is weak. In reality, there is a ridge (the Reykjanes Ridge) that runs southwest of Iceland, but the closed 620-m contour crosses this feature. It is diverted to the southwest as it crosses, but it is not completely impeded as in IDL-4.
5. Wind stress forcing
We have examined two mechanisms, that is, the basin circulation integral [(4)] and the around island integral [(9)], that promote the existence of the NIJ in source- and sink-driven overflows without wind stress forcing. Wind-driven flows contribute to the DSO, and wind stress is a main mechanism for transport variations in both the Denmark Strait and Faroe Bank Channel (e.g., Biastoch et al. 2003; Köhl et al. 2007; Nilsen et al. 2003; Serra et al. 2010). The role of wind stress forcing in the existence of the NIJ is examined in this section. We use annual-mean wind stress climatology from the objectively analyzed air–sea heat fluxes (OAFlux; Yu et al. 2008). It represents an average between 1988 and 2008 with a resolution of 0.25°. It is linearly interpolated to the model grids. In the first wind-driven experiment (WIND-1), the diapycnal water-mass flux term we is turned off, and the model is forced solely by the surface wind stress from an initial state of rest with the layer interface set at 500 m below sea surface. The circulations in both layers at the end of fifth year are shown in Fig. 12. The flows in both the upper and lower layers in the Nordic Seas are dominated by basinwide cyclonic gyres. They are driven by a positive curl of wind stress in the subpolar North Atlantic Ocean. The transport over the GISR is near zero in the lower layer since there are no water-mass sources or sinks in either basin.
There are several notable features in the lower-layer flows. First, the flow is always eastward on the north side of the Iceland–Faroe Ridge between the DS and FBC (Fig. 12b). In all source/sink-driven experiments, however, the flow is always westward (Figs. 3–5). The observed situation is apparently more complex than either of the above. Between the DS and Jan Mayen Ridge the flow is the westward NIJ. But to the east of the Jan Mayen Ridge, the flow is eastward according to Søiland et al. (2008), who used neutrally buoyant floats to trace the flow pathway. This current is the main pathway for the Faroe Bank Channel overflow (Søiland et al. 2008). None of the experiments described so far are able to simulate these two opposite currents. The flow is in the wrong direction between DS and Jan Mayen Ridge in the wind-driven run (Fig. 12) and between the Faroe Islands and Jan Mayen Ridge in source/sink-driven experiments (Figs. 3–5). This suggests that both wind stress and buoyancy forcing could be important for the Nordic Seas overflow pathways. To explore this further we add wind stress forcing to the previous experiment with a uniform we (EXP-2W). The flow in the lower layer is shown in Fig. 13. The circulation changes significantly when the wind stress is added. The two opposite currents along the northern Iceland–Faroe shelf/slope are simulated in EXP-2W with a uniform we (Fig. 13). The eastward flow along the Iceland–Faroe Ridge appears to be a continuation from a southward transport along the Jan Mayen Ridge. This pathway for the Faroe Bank Channel overflow is consistent with observations described by Søiland et al. (2008). The eastward jet along Iceland–Faroe Ridge also emerges in other source/sink-driven experiments when wind stress is added (not shown). Along the continental shelf north of Iceland, the addition of wind stress forcing in the model weakens but does not eliminate (at least for the climatological mean wind stress) the NIJ. The model NIJ extends far less eastward when compared with Fig. 4a. In fact, the westward flow along the northeastern coast of Iceland begins around 14°–15°W, which is roughly consistent with the observations by Våge et al. (2013). The NIJ transport appears to be supplied by a cyclonic gyre in the Iceland Sea and a small anticyclonic gyre between Jan Mayen Ridge and Iceland to the south of the Iceland Sea (Fig. 13).
A second notable feature of the wind-only experiment is that the EGC pathway becomes more prominent (Fig. 12) in both layers. The EGC is weak south of the Greenland Sea in the dam break and the source/sink-driven experiments (Figs. 3–4). It is strengthened as the western boundary current of the wind-driven cyclonic gyre in the Nordic Seas (Fig. 13). This would support previous modeling studies (e.g., Köhl 2010) suggesting that a strong wind stress forcing, such as during a high North Atlantic Oscillation state, enhances the EGC transport and weakens the NIJ pathway. It is interesting to note that there are three branches of the DSO, EGC, NIJ, and separated EGC (S-EGC) when wind stress is added. Without the wind stress, both the EGC and S-EGC are basically absent (Fig. 4b). So the wind stress is essential for the existence of both the EGC and S-EGC in the model, even though the formation of the S-EGC is mainly due to topographic effect.
Finally, we note that the deep circulations in the Greenland Sea, Norwegian Sea, Iceland Sea, and Lofoten basin become cyclonic when wind stress is added. They are mostly anticyclonic in previous source/sink-driven experiments (EXP-1, EXP-2, and EXP-3) without wind stress forcing. The eastern boundary current along Norway’s coast becomes northward in these wind-driven model runs. The flow along the midocean ridge, Jan Mayen and Mohn Ridges, is southward toward the GISR, which is consistent with what was found by Søiland et al. (2008). In the source/sink-driven flows, the flow is northward along the midocean ridge. In most cases with wind stress forcing discussed here, the lower-layer outcrops in the Nordic Seas and thus is exposed to wind stress forcing. This direct wind stress forcing contributes to circulation changes in the lower layer.
The wind stress forces a cyclonic boundary current that is opposite the NIJ along the north coast of Iceland. Figure 13 shows that the climatological wind stress in the model is not strong enough to reverse the NIJ. But one would wonder how resilient the NIJ is to an abnormally strong wind stress forcing, either because of a seasonal cycle or atmospheric variability. We conducted two additional experiments (not shown here), one with wind stress increased by a factor of 2 and the other by factor of 5. The westward NIJ is still present even when the climatological wind stress is doubled. But the flow along the north Icelandic coast changes to eastward when the wind stress is increased by a factor of 5.
6. Summary
A two-layer model is used to examine mechanisms that encourage the existence of the north Icelandic jet—ostensibly a major pathway for the Denmark Strait overflow. This study identifies and tests two integral constraints that bear on the existence of the NIJ. The first is a circulation integral around a closed contour C that circles those central portions of the upstream basin for which the depth is greater than the sill depth. An outward flow from these deeper regions implies a net divergence in the lateral PV transport in the deep Nordic basin, which implies an anticyclonic acceleration of the circulation around C. In a steady state this tendency must be balanced by a cyclonic frictional torque and implies anticyclonic flow around C. A westward NIJ is consistent with this picture, though not required, and the need for such a flow is strengthened by the presence of the (cyclonic) EGC. A similar result is obtained when the overflow is fed by a deep source whose thickness exceeds that of the overflow. Although both cases are suggestive, we also find a westward NIJ when sinking occurs, mainly in the shallower regions outside of C. In this case, the circulation integral around C is uninformative. The second, stronger constraint is contained in the circulation integral around a closed geostrophic contour CI that lies slightly above the sill depth of the Denmark Strait and circles around Iceland and the Iceland–Faroe Ridge, passing unimpeded through the deeper Faroe Bank Channel. If there is no wind and no net exchange between the upper and lower layers inside CI the circulation around CI is zero and so the tangential component of the velocity about it cannot be unidirectional. The FBC overflow produces a predominantly anticyclonic tangential velocity along the eastern and southern segments of CI, and an NIJ-like western flow along the segment north of Iceland would tend to balance this.
A positive wind stress curl forces cyclonic circulations in both layers. The basinwide circulation, including flows along the midocean ridges, is strongly affected by wind stress. The wind-driven surface flow along the north Icelandic shelf is eastward and against the NIJ. The wind weakens but does not eliminate the NIJ, and it produces a more realistic eastward flow in the lower layer to the east of the Jan Mayen Ridge. The presence of the wind also produces cyclonic flow in the Lofoten, Norwegian, and Greenland basins, which is in general agreement with observations. The wind stress also promotes the EGC pathway, and it is postulated that the NIJ pathway would be more dominant when wind stress forcing is weak, such as during a low state of the North Atlantic Oscillation, while the EGC becomes more important when wind is strong.
A reviewer has asked us to comment on the fact that net downwelling can become smaller over time because of the outcropping of the lower layer. In our runs, where the initial upper-layer thickness is hundreds of meters, outcropping occurs in about 2 yr, whereas the steady-state solutions that we actually analyze are examined after 10 yr or so. So the preoutcropping phase is quite short compared to the total run time. Accordingly, we have not attempted to analyze transients and other adjustment features that are triggered when outcropping occurs and the total downwelling decreases. This is an interesting problem, however, and perhaps one of the keys to variability of the Nordic Seas outflows. It is beyond the scope of the current paper, but a good problem for future study.
Acknowledgments
This study has been supported by National Science Foundation (OCE0927017 and ARC1107412). We thank Drs. Femke de Jong, Bob Pickart, Mike Spall, and two thorough reviewers for their comments and suggestions.
REFERENCES
Aagaard, K., E. Fahrbach, J. Meincke, and J. H. Swift, 1991: Saline outflow from the Arctic Ocean: Its contribution to the deep waters of the Greenland, Norwegian, and Icelandic seas. J. Geophys. Res., 96, 20 433–20 441, doi:10.1029/91JC02013.
Biastoch, A., R. H. Käse, and D. B. Stammer, 2003: The sensitivity of the Greenland–Scotland Ridge overflow to forcing changes. J. Phys. Oceanogr., 33, 2307–2319, doi:10.1175/1520-0485(2003)033<2307:TSOTGR>2.0.CO;2.
Böning, C. W., F. O. Bryan, W. R. Holland, and R. Doescher, 1996: Deep-water formation and meridional overturning in a high-resolution model of the North Atlantic. J. Phys. Oceanogr., 26, 1142–1164, doi:10.1175/1520-0485(1996)026<1142:DWFAMO>2.0.CO;2.
Buch, E., S.-A. Malmberg, and S. S. Kristmannsson, 1996: Arctic Ocean deep water masses in the western Iceland Sea. J. Geophys. Res., 101, 11 965–11 973, doi:10.1029/95JC03869.
Cenedese, C., 2012: Downwelling in basins subject to buoyancy loss. J. Phys. Oceanogr., 42, 1817–1833, doi:10.1175/JPO-D-11-0114.1.
Dickson, R., and Coauthors, 2008: The overflow flux west of Iceland: Variability, origins and forcing. Arctic–Subarctic Ocean Flux, R. R. Dickson, J. Meincke, and P. Rhines, Eds., Springer, 443–474, doi:10.1007/978-1-4020-6774-7_20.
Eldevik, T., J. E. Ø. Nilsen, D. Iovino, K. A. Olsson, A. B. Sandø, and H. Drange, 2009: Observed sources and variability of Nordic Seas overflow. Nat. Geosci., 2, 406–410, doi:10.1038/ngeo518.
Hansen, B., S. Østerhus, W. R. Turrell, S. Jónsson, H. Valdimarsson, H. Hátún, and S. M. Olsen, 2008: The inflow of Atlantic water, heat and salt to the Nordic Seas across the Greenland–Scotland Ridge. Arctic–Subarctic Ocean Flux, R. R. Dickson, J. Meincke, and P. Rhines, Eds., Springer, 15–43.
Helfrich, K. R., and L. J. Pratt, 2003: Rotating hydraulics and the upstream basin circulation. J. Phys. Oceanogr., 33, 1651–1663, doi:10.1175/2383.1.
Isachsen, P. E., C. Mauritzen, and H. Svendsen, 2007: Dense water formation in the Nordic Seas diagnosed from sea surface buoyancy fluxes. Deep-Sea Res., 54, 22–41, doi:10.1016/j.dsr.2006.09.008.
Jeansson, E., S. Jutterström, B. Rudels, L. G. Anderson, K. A. Olsson, E. P. Jones, W. M. Smethie Jr., and J. H. Swift, 2008: Sources to the East Greenland Current and its contribution to the Denmark Strait overflow. Prog. Oceanogr., 78, 12–28, doi:10.1016/j.pocean.2007.08.031.
Jónsson, S., 1999: The circulation in the northern part of the Denmark Strait and variability. ICES CM 9, 9 pp.
Jónsson, S., and H. Valdimarsson, 2004: A new path for the Denmark Strait Overflow Water from the Iceland Sea to Denmark Strait. Geophys. Res. Lett., 31, L03305, doi:10.1029/2003GL019214.
Karcher, M., F. Kauker, R. Gerdes, E. Hunke, and J. Zhang, 2007: On the dynamics of Atlantic Water circulation in the Arctic Ocean. J. Geophys. Res., 112, C04S02, doi:10.1029/2006JC003630.
Käse, R. H., and A. Oschlies, 2000: Flow through Denmark Strait. J. Geophys. Res.,105, 28 527–28 546, doi:10.1029/2000JC900111.
Käse, R. H., J. B. Girton, and T. B. Sanford, 2003: Structure and variability of the Denmark Strait overflow: Model and observations. J. Geophys. Res., 108, 3181, doi:10.1029/2002JC001548.
Käse, R. H., N. Serra, A. Köhl, and D. Stammer, 2009: Mechanisms for the variability of dense water pathways in the Nordic Seas. J. Geophys. Res., 114, C01013, doi:10.1029/2008JC004916.
Köhl, A., 2010: Variable source regions of Denmark Strait and Faroe Bank Channel. Tellus, 62A, 551–568, doi:10.1111/j.1600-0870.2010.00454.x.
Köhl, A., R. H. Käse, D. Stammer, and N. Serra, 2007: Causes of changes in the Denmark Strait overflow. J. Phys. Oceanogr., 37, 1678–1696, doi:10.1175/JPO3080.1.
Macrander, A., U. Send, H. Valdimarsson, S. Jónsson, and R. H. Käse, 2005: Interannual changes in the overflow from the Nordic Seas into the Atlantic Ocean through Denmark Strait. Geophys. Res. Lett., 32, L06606, doi:10.1029/2004GL021463.
Macrander, A., R. H. Käse, U. Send, H. Valdimarsson, and S. Jonsson, 2007: Spatial and temporal structure of the Denmark Strait overflow by acoustic observations. Ocean Dyn., 57, 75–89, doi:10.1007/s10236-007-0101-x.
Marotzke, J., and J. R. Scott, 1999: Convective mixing and the thermohaline circulation. J. Phys. Oceanogr., 29, 2962–2970, doi:10.1175/1520-0485(1999)029<2962:CMATTC>2.0.CO;2.
Mauritzen, C., 1996a: Production of dense overflow waters feeding the North Atlantic across the Greenland-Scotland Ridge. Part 1: Evidence for a revised circulation scheme. Deep-Sea Res. I, 43, 769–806, doi:10.1016/0967-0637(96)00037-4.
Mauritzen, C., 1996b: Production of dense overflow waters feeding the North Atlantic across the Greenland Sea-Scotland Ridge. Part 2: An inverse model. Deep-Sea Res., 43, 807–835, doi:10.1016/0967-0637(96)00038-6.
Nilsen, J. E. Ø., Y. Gao, H. Drange, T. Furevik, and M. Bentsen, 2003: Simulated North Atlantic-Nordic Seas water mass exchanges in an isopycnic coordinate OGCM. Geophys. Res. Lett., 30, 1536, doi:10.1029/2002GL016597.
Pedlosky, J., 2003: Thermally driven circulation in small ocean basins. J. Phys. Oceanogr., 33, 2333–2340, doi:10.1175/1520-0485(2003)033<2333:TDCISO>2.0.CO;2.
Pratt, L. J., 1997: Hydraulically drained flows in rotating basins. Part II: Steady flow. J. Phys. Oceanogr., 27, 2522–2535, doi:10.1175/1520-0485(1997)027<2522:HDFIRB>2.0.CO;2.
Pratt, L. J., and M. Spall, 2008: Circulation and exchange in choked marginal seas. J. Phys. Oceanogr., 38, 2639–2661, doi:10.1175/2008JPO3946.1.
Pratt, L. J., and J. Whitehead, 2008: Rotating Hydraulics. Springer, 582 pp.
Price, J. F., and M. O. Baringer, 1994: Outflows and deep water production by marginal seas. Prog. Oceanogr., 33, 161–200, doi:10.1016/0079-6611(94)90027-2.
Rudels, B., E. Fahrbach, E. Meincke, G. Budeus, and P. Eriksson, 2002: The East Greenland Current and its contribution to the Denmark Strait overflow. ICES J. Mar. Sci., 59, 1133–1154, doi:10.1006/jmsc.2002.1284.
Rudels, B., P. Eriksson, E. Buch, G. Budéus, E. Fahrbach, S.-A. Malmberg, J. Meincke, and P. Mälkki, 2003: Temporal switching between sources of the Denmark Strait Overflow Water. ICES Mar. Sci. Symp., 219, 319–325.
Serra, N., R. H. Käse, A. Kohl, D. Stammer, and D. Quadfadsel, 2010: On the low-frequency phase relation between the Denmark Strait and the Faroe-Bank Channel overflows. Tellus, 62A, 530–550, doi:10.1111/j.1600-0870.2010.00445.x.
Smethie, W. M., and J. H. Swift, 1989: The tritium:krypton-85 age of Denmark Strait Overflow Water and the Gibbs Fracture Zone Water just south of Denmark Strait. J. Geophys. Res., 94, 8265–8275, doi:10.1029/JC094iC06p08265.
Søiland, H., M. D. Prater, and T. Rossby, 2008: Rigid topographic control of currents in the Nordic Seas. Geophys. Res. Lett., 35, L18607, doi:10.1029/2008GL034846.
Spall, M., 2003: On the thermohaline circulation in flat-bottom marginal seas. J. Mar. Res., 61, 1–25, doi:10.1357/002224003321586390.
Spall, M., 2004: Boundary currents and watermass transformation in marginal seas. J. Phys. Oceanogr., 34, 1197–1213, doi:10.1175/1520-0485(2004)034<1197:BCAWTI>2.0.CO;2.
Spall, M., 2010: Non-local topographic influences on deep convection: An idealized model for the Nordic Seas. Ocean Modell., 32, 72–85, doi:10.1016/j.ocemod.2009.10.009.
Straneo, F., 2006: On the connection between dense water formation, overturning, and poleward heat transport in a convective basin. J. Phys. Oceanogr., 36, 606–628, doi:10.1175/JPO2875.1.
Swift, J. H., K. Aagaard, and S. Malmberg, 1980: The contribution of Denmark Strait overflow to the deep North Atlantic. Deep-Sea Res., 27, 29–42, doi:10.1016/0198-0149(80)90070-9.
Tanhua, T., K. A. Olsson, and E. Jeansson, 2005: Formation of Denmark Strait Overflow Water and its hydro-chemical composition. J. Mar. Syst., 57, 264–288, doi:10.1016/j.jmarsys.2005.05.003.
Våge, K., R. S. Pickart, M. A. Spall, H. Valdimarsson, S. Jonsson, D. J. Torres, S. Osterhus, and T. Eldevik, 2011: Significant role of the north Icelandic jet in the formation of Denmark Strait Overflow Water. Nat. Geosci., 4, 723–727, doi:10.1038/ngeo1234.
Våge, K., R. S. Pickart, M. A. Spall, G. W. K. Moore, H. Valdimarsson, D. J. Torres, S. Y. Erofeeva, and J. E. Nilsen, 2013: Revised circulation scheme north of the Denmark Strait. Deep-Sea Res., 79, 20–39, doi:10.1016/j.dsr.2013.05.007.
Wilkenskjeld, S., and D. Quadfasel, 2005: Response of the Greenland‐Scotland overflow to changing deep water supply from the Arctic Mediterranean. Geophys. Res. Lett.,32, L21607, doi:10.1029/2005GL024140.
Yang, J., 2005: The Arctic and subarctic ocean flux of potential vorticity and the Arctic Ocean circulation. J. Phys. Oceanogr., 35, 2387–2407, doi:10.1175/JPO2819.1.
Yang, J., and J. F. Price, 2000: Water mass formation and potential vorticity balance in an abyssal ocean circulation model. J. Mar. Res., 58, 789–808, doi:10.1357/002224000321358918.
Yang, J., and J. F. Price, 2007: Potential vorticity constraint on the flow between two basins. J. Phys. Oceanogr., 37, 2251–2266, doi:10.1175/JPO3116.1.
Yang, J., and L. Pratt, 2013: On the effective capacity of the dense-water reservoir for the Nordic Seas overflow: Some effects of topography and wind stress. J. Phys. Oceanogr., 43, 418–431, doi:10.1175/JPO-D-12-087.1.
Yu, L., X. Jin, and R. A. Weller, 2008: Multidecade global flux datasets from the Objectively Analyzed Air-sea Fluxes (OAFlux) project: Latent and sensible heat fluxes, ocean evaporation, and related surface meteorological variables. Woods Hole Oceanographic Institution OAFlux Project Tech. Rep. OA-2008-01, 64 pp.
The DSO pathways, which are centered along 600–700-m isobaths, are not sensitive to excluding shelf seas that are shallower than 200 m.
Although our governing equations do not resolve bottom Ekman layers, a linear bottom drag law would produce a cross-isobath flux that could be made consistent with a bottom Ekman layer through the correct choice of the coefficient in the linear bottom drag law. This is discussed in Pratt (1997).