A three-dimensional surface wave–ocean circulation coupled model and its initial testing

A theoretical framework to include the influences of nonbreaking surface waves in ocean general circulation models is established based on Reynolds stresses and fluxes terms derived from surface wave-induced fluctuation. An expression for the wave-induced viscosity and diffusivity as a function of the wave number spectrum is derived for infinite and finite water depths; this derivation allows the coupling of ocean circulation models with a wave number spectrum numerical model. In the case of monochromatic surface wave, the wave-induced viscosity and diffusivity are functions of the Stokes drift. The influence of the wave-induced mixing scheme on global ocean circulation models was tested with the Princeton Ocean Model, indicating significant improvement in upper ocean thermal structure and mixed layer depth compared with mixing obtained by the Mellor–Yamada scheme without the wave influence. For example, the model–observation correlation coefficient of the upper 100-m temperature along 35° N increases from 0.68 without wave influence to 0.93 with wave influence. The wave-induced Reynolds stress can reach up to about 5% of the wind stress in high latitudes, and drive 2–3 Sv transport in the global ocean in the form of mesoscale eddies with diameter of 500–1,000 km. The surface wave-induced mixing is more pronounced in middle and high latitudes during the summer in the Northern Hemisphere and in middle latitudes in the Southern Hemisphere.     

Numerical study of surface water circulation around Sekisei Lagoon, southwest Japan

This paper discusses the variability of surface currents around Sekisei Lagoon using a nested grid ocean circulation model. We developed a triple-nested grid system that consists of a coarse-resolution (1/60° or ∼1.85 km) model off Taiwan, an intermediate-resolution (1/300° or ∼370 m) model around the Yaeyama Islands, and a fine-resolution (1/900° or ∼123 m) model of Sekisei Lagoon. The nested grid system was forced by wind and heat flux calculated from six-hourly atmospheric reanalysis data and integrated over the period from May to July 2003. The coarse-resolution model was driven by lateral boundary conditions calculated from daily ocean reanalysis data to include realistic variation of the Kuroshio and mesoscale eddies with spatial scales of ∼500–700 km in the open ocean. The tidal forcing was included in the intermediate-resolution model by interpolating sea level data obtained from a data-assimilative tidal model. The results were then used to drive the fine-resolution model to simulate the surface water circulation around Sekisei lagoon. Model results show that (1) currents inside the lagoon are mainly driven by tide and wind; (2) there exists a strong southwestward current along the bottom slope in the southeast portion of the lagoon; the current is mainly driven by remote mesoscale eddies and at times intensified by the local wind; (3) the flow relaxation scheme is effective in reducing biases along the open boundaries. The simulated currents were used to examine the retention and dispersion of passive particles in the surface layer. Results show that the surface dispersion in the strong open ocean current region is significantly higher than that inside the lagoon.     

Multi-nest high-resolution model of submesoscale circulation features in the Gulf of Taranto

Recent oceanographic field measurements and high-resolution numerical modelling studies have revealed intense, transient, submesoscale motions characterised by a horizontal length scale of 100–10,000 m. This submesoscale activity increases in the fall and winter when the mixed layer (ML) depth is at its maximum. In this study, the submesoscale motions associated with a large-scale anticyclonic gyre in the central Gulf of Taranto were examined using realistic submesoscale-permitting simulations. We used realistic flow field initial conditions and multiple nesting techniques to perform realistic simulations, with very-high horizontal resolutions (> 200 m) in areas with submesoscale variability. Multiple downscaling was used to increase resolution in areas where instability was active enough to develop multi-scale interactions and produce 5-km-diameter eddies. To generate a submesoscale eddy, a 200-m resolution was required. The submesoscale eddy was formed through small-scale baroclinic instability in the rim of a large-scale anticyclonic gyre leading to large vertical velocities and rapid restratification of the ML in a time-scale of days. The submesoscale eddy was confirmed by observational data from the area and we can say that for the first time we have a proof that the model reproduces a realistic submesoscale vortex, similar in shape and location to the observed one.     

Observations of submesoscale eddies using high-frequency radar-derived kinematic and dynamic quantities

The spatio-temporal variability of submesoscale eddies off southern San Diego is investigated with two-year observations of subinertial surface currents [O(1) m depth] derived from shore-based high-frequency radars. The kinematic and dynamic quantities — velocity potential, stream function, divergence, vorticity, and deformation rates — are directly estimated from radial velocity maps using optimal interpolation. For eddy detection, the winding-angle approach based on flow geometry is applied to the calculated stream function. A cluster of nearly enclosed streamlines with persistent vorticity in time is identified as an eddy. About 700 eddies were detected for each rotation (clockwise and counter-clockwise). The two rotations show similar statistics with diameters in the range of 5–25 km and Rossby number of 0.2–2. They persist for 1–7 days with weak seasonality and migrate with a translation speed of 4–15 cm s⁻¹ advected by background currents. The horizontal structure of eddies exhibits nearly symmetric tangential velocity with a maximum at the defined radius of the eddy, non-zero radial velocity due to background flows, and Gaussian vorticity with the highest value at the center. In contrast divergence has no consistent spatial shape. Two episodic events are presented with other in situ data (subsurface current and temperature profiles, and local winds) as an example of frontal-scale secondary circulation associated with drifting submesoscale eddies.     

Three-dimensional structure of mesoscale eddies in the western tropical Pacific as revealed by a high-resolution ocean simulation

The three-dimensional structure of mesoscale eddies in the western tropical Pacific(6°S–20°N, 120°E–150°E)is investigated using a high-resolution ocean model simulation. Eddy detection and eddy tracking algorithms are applied to simulated horizontal velocity vectors, and the anticyclonic and cyclonic eddies identified are composited to obtain their three-dimensional structures. The mean lifetime of all long-lived eddies is about 52 days, and their mean diameter is 147 km. Two typical characteristics of mesoscale eddies are revealed and possible dynamic explanations are analyzed. One typical characteristic is that surface eddies are generally separated from subthermocline eddies along the bifurcation latitude(~13°N) of the North Equatorial Current in the western tropical Pacific, which may be associated with different eddy energy sources and vertical eddy energy fluxes in subtropical and tropical gyres. Surface eddies have maximum swirl velocities of 8–9 cm s~(-1) and can extend to about 1500 m depth. Subthermocline eddies occur below 200 m, with their cores at about 400–600 m depth, and their maximum swirl velocities can reach 10 cm s~(-1). The other typical characteristic is that the meridional velocity component of the eddy is much larger than the zonal component. This characteristic might be due to more zonal eddy pairs(two eddies at the same latitude),which is also supported by the zonal wavelength(about 200 km) in the high-frequency meridional velocity component of the horizontal velocity.     

中美洲鹦鹉湾气旋涡的亚中尺度结构特征

海洋中普遍存在的涡旋对全球海洋热、盐通量有重要贡献.一条于2000年6月在中美洲鹦鹉湾采集的地震剖面L115捕获到了一个海水次表层透镜状反射结构.透镜状反射的宽度约150 km,厚度约500 m,深度从100 m延伸到约600 m,核心深约200 m.结合和地震采集时间近同步的再分析数据中的流速和海水温度数据以及计算的Rossby数和Okubo-Weiss参数,将其解释为气旋涡,且可能是亚中尺度相干涡(SCV).由于地震剖面并未穿过涡旋核心水,涡旋顶底与核心水外层水团的双扩散作用使得剖面上的气旋涡中心表现为上下相邻的强反射双核结构.同航次采集的与L115大致正交的其他几条地震剖面也捕捉到了这个气旋涡.这些剖面整体的反射特征一致性较好,但和L115差异较大,涡旋内部普遍表现为近水平的中等强度反射.涡旋上边界为倾斜的强反射,这是气旋涡的等温线上凸引起的,这里也是亚中尺度锋面的发育区.这些地震剖面的涡旋反射特征的差异表明了该气旋涡空间结构的不对称性和时间演化特征.     

Mechanisms controlling the seasonal mixed-layer temperature and salinity of the Indonesian seas

We examine the seasonal mixed-layer temperature (MLT) and salinity (MLS) budgets in the Banda–Arafura Seas region (120–138° E, 8–3° S) using an ECCO ocean-state estimation product. MLT in these seas is relatively high during November–May (austral spring through fall) and relatively low during June–September (austral winter and the period associated with the Asian summer monsoon). Surface heat flux makes the largest contribution to the seasonal MLT tendency, with significant reinforcement by subsurface processes, especially turbulent vertical mixing. Temperature declines (the MLT tendency is negative) in May–August when seasonal insolation is smallest and local winds are strong due to the southeast monsoon, which causes surface heat loss and cooling by vertical processes. In particular, Ekman suction induced by local wind stress curl raises the thermocline in the Arafura Sea, bringing cooler subsurface water closer to the base of the mixed layer where it is subsequently incorporated into the mixed layer through turbulent vertical mixing; this has a cooling effect. The MLT budget also has a small, but non-negligible, semi-annual component since insolation increases and winds weaken during the spring and fall monsoon transitions near the equator. This causes warming via solar heating, reduced surface heat loss, and weakened turbulent mixing compared to austral winter and, to a lesser extent, compared to austral summer. Seasonal MLS is dominated by ocean processes rather than by local freshwater flux. The contributions by horizontal advection and subsurface processes have comparable magnitudes. The results suggest that ocean dynamics play a significant part in determining both seasonal MLT and MLS in the region, such that coupled model studies of the region should use a full ocean model rather than a slab ocean mixed-layer model.     

Effects of cyclonic mesoscale eddies on the marine ecosystem in the Kuroshio Extension region using an eddy-resolving coupled physical-biological model

Effects of mesoscale eddies on the marine ecosystem in the Kuroshio Extension (KE) region are investigated using an eddy-resolving coupled physical-biological model. The model captures the seasonal and intra-seasonal variability of chlorophyll distribution associated with the mesoscale eddies, front variability, Kuroshio meanders, and upwelling. The model also reproduces the observed interannual variability of sea surface height anomaly (SSHA) in the KE region along a zonal band of 32–34°N from 2002 to 2006. The distribution of high surface chlorophyll corresponds to low SSHA. Cyclonic eddies are found to detach from the KE jet near 150°E and 158°E and propagate westward. The westward propagating cyclonic eddies lift the nutrient-rich thermocline into the euphotic zone and maintain high levels of chlorophyll in summer. In the subsurface layer, the pattern in chlorophyll is influenced by both lateral and vertical advection. In winter, convection inside the eddy entrains high levels of nutrients into the mixed layer, increasing production, and resulting in high chlorophyll concentration throughout the surface mixed layer. There is significant interannual variability in both the cyclonic eddy activity and the surface phytoplankton bloom south of the KE jet, although whether or not there is a causal link is unclear.     

Numerical simulation of the circulation within the Perth Submarine Canyon,Western Australia

Surface and sub-surface currents along the ocean boundary of Western Australia were simulated using Regional Ocean Modelling System (ROMS) to examine the circulation within the Perth Canyon. Two major current systems influenced the circulation within the canyon: (1) The Leeuwin current interacted weakly with the canyon as the majority of the canyon was below the depth of the Leeuwin current and (2) Leeuwin undercurrent interacted strongly with the canyon, forming eddies within the canyon at depths of 400–800 m. The results indicated that within the canyon, the current patterns changed continuously although there were some repeated patterns. Recurrent eddies produced regions where upwelling or downwelling dominated during the model runs. Deep upwelling was stronger within the canyon than elsewhere on the shelf, but vertical transport in the upper ocean was strong everywhere when wind forcing was applied. Upwelling alone appeared to be insufficient to transport nutrients to the euphotic zone because the canyon rims were deep. Increased upwelling, combined with entrapment within eddies and strong upwelling-favourable winds, which could assist mixing, may account for the high productivity attributed to the canyon. The Leeuwin current is otherwise a strong barrier to the upwelling of nutrients.     

The mesoscale eddies and Kuroshio transport in the western North Pacific east of Taiwan from 8-year (2003–2010) model reanalysis

The relationship between the Kuroshio volume transport east of Taiwan (~24°N) and the impinging mesoscale eddies is investigated using 8-year reanalysis of a primitive equation ocean model that assimilates satellite altimetry and SST data. The mean and fluctuations of the model Kuroshio transport agree well with the available observations. Analysis of model dynamic heights and velocity fields reveals three dominant eddy modes. The first mode describes a large eddy of ~500 km in diameter, centered at ~22° N. The second mode describes a pair of the north–south counter-rotating eddies of?~?400 km in diameter each, centered at 23° and 20° N, respectively. The third mode describes a pair of the east–west counter-rotating eddies of?~?300 km in diameter each, centered at 21° N. The associated velocity fields indicate eddies extending to 600–700 m in depth with vertical shears concentrated in the upper 400 m. All three modes and the model Kuroshio transport have similar dominant timescales of 70–150 days and generally are coherent. The decreased Kuroshio volume transports typically are associated with the impinging cyclonic eddies and the increased transports with the anticyclonic eddies. Selected drifter trajectories are presented to illustrate the three eddy modes and their correspondence with the varying Kuroshio transports.     

Impacts of the Kuroshio intrusion on the two eddies in the northern South China Sea in late spring 2016

During mid-May to early June 2016, a cold eddy and a warm eddy were captured on the continental slope of the northern South China Sea by the in situ measurements. A salty lens-shaped water mass in the subsurface layer existed in these two detected eddies, which indicated they had a Kuroshio water origin. The trajectories of the observed eddies from satellite altimeter data show that the cold eddy was generated in the central part of the Luzon Strait, while the warm eddy was formed southwest of Taiwan. The genesis of the cold eddy is related to a weak Kuroshio loop current, while that of the warm eddy is associated with a strong Kuroshio loop current. The warm eddy east of the Luzon Strait may trigger the Kuroshio from a leaping path to a looping path. During the evolution of these detected eddies, they had interactions with the Kuroshio and Luzon Gyre. Energy analysis from ocean reanalysis data showed that the baroclinic conversion between the cold eddy and the Kuroshio was stronger than that between the cold eddy and Luzon Gyre. During the eddy shedding stage, the warm eddy mainly acquired energy from the Kuroshio loop current through the baroclinic conversion.     

Rapid intensification of Super Typhoon Haiyan: the important role of a warm-core ocean eddy

Super Typhoon Haiyan devastated portions of Southeast Asia, particularly the Philippines, on November 8, 2013. In this paper, observational data are used to analyze the intensification process of Super Typhoon Haiyan. Observational data showed that Typhoon Haiyan intensified and the maximum sustained winds increased to 59 ms^?1 after it encountered a double warm-core ocean eddy, while the central pressure of the typhoon dropped from 970 hPa to 920 hPa. Numerical simulations and observational data show that the presence of the warm-core eddy combined with SST increases due to climate change led to the rapid intensification of Super Typhoon Haiyan. Comparing these two factors, the warm-core ocean eddy, which brings significantly more heat into the upper ocean, plays the leading role in the intensification, with climate warming making a lesser contribution. Moreover, due to the increased thickness of the mixed layer associated with the warm-core ocean eddy, Super Typhoon Haiyan did not significantly decrease the sea surface temperature to the east of the Philippines, as is typical of typhoons, and the largest decrease was approximately 1 °C.     

Cold anticyclonic eddies formed from cold pool water in the southern Middle Atlantic Bight

AVHRR satellite imagery of the southern Mid-Atlantic Bight during May 1993 revealed a large area of cold water over the shelf break and slope that appeared to spin up into a series of southward propagating anticyclonic eddies. The eddies had diameters of 35–45 km at the surface and moved southward at about 20 cm/sec. A radial TOYO CTD (to 50m) and ADCP velocity (to 400m) transect was conducted across the southern-most of these eddies. The upper 50 meters had minimum temperatures of less than 7°C and salinities of about 33 pss, characteristics similar to cold pool waters usually found over the continental shelf. ADCP velocity data from one of the eddies revealed anticyclonic flow extending to a depth of about 250m. The transport of cold pool water by the eddies was estimated to be 0.1 to 0.2 Sv which is of the same order as the annual mean alongshore transport of shelf water in this region. The origin of the deeper water within the eddy is unlikely to be the continental shelf because the shelf break is less than 100 m. The depth and velocity profiles along the TOYO transect were consistent with the constant potential vorticity eddy model of Flierl (1979) although the source of the eddy kinetic energy is uncertain. The cause for the exodus of cold pool water from the shelf, which extended northward to at least 38°N, is unclear but must involve the establishment of an alongshore baroclinic pressure gradient against the usual southwestward shelf flow. It is possible that the intrusion of Gulf Stream waters onto the shelf near Cape Hatteras was a precursor of this off shelf transport. The southern-most eddy was marked by high biological productivity and very high oxygen supersaturation. The phytoplankton bloom detected within the exported cold pool water, located over the continental slope, suggests a mechanism whereby production fueled by nutrients derived from the shelf can be locally exported into deep water.     

The eddy pair off eastern Vietnam: Interannual variability and impact on thermohaline structure

The main features, interannual variability of the eddy pair (an anti-cyclonic southern eddy and a cyclonic northern eddy) off eastern Vietnam and its impact on local thermohaline structure are analyzed using altimeter data, the Simple Ocean Data Assimilation (SODA) data, and the Medium-Range Weather Forecasts (ECWMF) data. The results show that the eddy pair is a seasonal phenomenon usually found in the summer to fall. It can reach downward to about 400 m depth but is strongest at the upper 100 m depth. In addition, the eddy pair displays a remarkable interannual variability with periods of 5.6 and 3.6 years. The local wind stress curl plays an important role in the interannual variability of eddy pair. Moreover, the eddy pair has an obvious impact on the thermohaline structure of the local upper ocean. The eddy pair can lead to stronger temperature variability in the upper 380 m but stronger salinity variability only in the upper 150 m depth. When the eddy pair is strong, the anti-cyclonic eddy deepens the thermocline, while the cyclonic eddy shoals the thermocline. Taking the depth of 20 °C isotherm as the thermocline depth, we find that in September, the thermocline drops downward to 110–140 m depth in the anti-cyclonic eddy but rises up to 50–80 m depth in the cyclonic eddy.     

Benchmarking the mesoscale variability in global ocean eddy-permitting numerical systems

The role of data assimilation procedures on representing ocean mesoscale variability is assessed by applying eddy statistics to a state-of-the-art global ocean reanalysis (C-GLORS), a free global ocean simulation (performed with the NEMO system) and an observation-based dataset (ARMOR3D) used as an independent benchmark. Numerical results are computed on a 1/4 ^° horizontal grid (ORCA025) and share the same resolution with ARMOR3D dataset. This “eddy-permitting” resolution is sufficient to allow ocean eddies to form. Further to assessing the eddy statistics from three different datasets, a global three-dimensional eddy detection system is implemented in order to bypass the need of regional-dependent definition of thresholds, typical of commonly adopted eddy detection algorithms. It thus provides full three-dimensional eddy statistics segmenting vertical profiles from local rotational velocities. This criterion is crucial for discerning real eddies from transient surface noise that inevitably affects any two-dimensional algorithm. Data assimilation enhances and corrects mesoscale variability on a wide range of features that cannot be well reproduced otherwise. The free simulation fairly reproduces eddies emerging from western boundary currents and deep baroclinic instabilities, while underestimates shallower vortexes that populate the full basin. The ocean reanalysis recovers most of the missing turbulence, shown by satellite products , that is not generated by the model itself and consistently projects surface variability deep into the water column. The comparison with the statistically reconstructed vertical profiles from ARMOR3D show that ocean data assimilation is able to embed variability into the model dynamics, constraining eddies with in situ and altimetry observation and generating them consistently with local environment.     

Observations of the upper ocean response to storm forcing in the South Atlantic Roaring Forties

In the austral summer of 1992–1993 the passage of a storm system drove a strong upper ocean response at 45°S in the mid-South Atlantic. Good in situ observations were obtained. CTD casts revealed that the mixed layer deepened by \sim40 m over 4 days. Wind stirring dominated over buoyancy flux-driven mixing during the onset of high winds. Doppler shear currents further reveal this to be intimately related to inertial dynamics. The penetration depth of inertial currents, which are confined to the mixed layer, increases with time after a wind event, matched by a downward propagation of low values of the Richardson number. This suggests that inertial current shear is instrumental in producing turbulence at the base of the mixed layer. Evolution of inertial transport is simulated using a time series of ship-observed wind stress. Simulated transport is only 30-50% of the observed transport, suggesting that much of the observed inertial motion was forced by an earlier (possibly remote) storm. Close proximity of the subtropical front further complicates the upper ocean response to the storm. A simple heat balance for the upper 100 m reveals that surface cooling and mixing (during the storm) can account for only a small fraction of an apparent \sim1 °C mixed layer cooling.     

Eddy resolving modelling of the Gulf of Lions and Catalan Sea

The generation process of strong long-lived eddies flowing southwestwards along the Catalan slope was revealed through numerical modelling and in situ observations. Careful analyses of a particular event in autumn 2007 demonstrated a link between a “LATEX” eddy, which remained in the southwestern corner of the Gulf of Lions and a “CATALAN” eddy, which moved along the Catalan Shelf, since the death of the former gave birth to the latter. The origin of such eddies was found to be an accumulation of potential energy in the southwestern corner of the Gulf of Lions: under the influence of the negative wind stress curl associated with the Tramontane, a warm and less dense water body can be isolated and fed by a coastal current carrying warm water from the Catalan Sea. In summer, this structure can grow and intensify to generate a strong anticyclonic eddy. After a long period of Tramontane, a burst of southeasterlies and northerlies appeared to detach the “LATEX” eddy, which flowed out of the Gulf of Lions, migrating along the Catalan continental slope and continued into the Balearic Sea as the “CATALAN” eddy.     

Origin and mechanism of Subantarctic Mode Water formation and transformation in the Southern Indian Ocean

The sources and pathways of mode waters and lower thermocline waters entering the subtropical gyre of the Indian Ocean are examined. A Lagrangian analysis is performed on an eddy-admitting simulation of the Global Ocean performed by the DRAKKAR Group (NEMO/OPA), which captures the main observed features. We trace the subducted mode water’s pathways, identify their formation regions and trace whether their source waters come from the Atlantic, Pacific or Indian sectors of the Southern Ocean. Three main sites for mode waters ventilation in the Indian sector are identified with different circulation pathways and source water masses: (a) just north of Kerguelen, where 4.2 Sv of lighter Subantarctic Mode Waters (SAMW); σ ₀ ∼ 26.5) are exported—originating in the Atlantic and Agulhas Retroflection regions; (b) SW of Australia, where 6.5 Sv of medium SAMW (σ ₀ ∼ 26.6) are ventilated—originating in the southern and denser Agulhas Retroflection region; (c) SW of Tasmania and along the South Australian coast, where 3 Sv of denser SAMW (σ ₀ ∼ 26.75) are ventilated—originating from three sources: Leeuwin Current waters, Tasman Sea (Pacific) waters and Antarctic Surface Waters. In all cases, modelled mode waters were last ventilated in the Indian Ocean just north of the deepest winter-mixed layers. For the waters subducted SW of Australia, the last ventilation site extends even further north. Waters ventilated in the deepest mixed layers north of the Subantarctic Front are then re-ventilated 5 years later southwest of Australia. The model results raise new hypotheses that revisit the classical picture of the SAMW formation and transformation, where a large homogeneous mixed layer is subducted and ‘slides’ equatorward, essentially maintaining the T/S characteristics acquired at the surface. Firstly, the last ventilation of the modelled mode waters is not in the region of the deepest mixed layers, as previously thought, but further north in regions of moderate meso-scale eddy activity. Secondly, the model shows for the first time a significant source region for Indian Ocean mode waters coming from deep winter-mixed layers along the south Australian coast. Finally, this analysis shows how the mode water characteristics are modified after subduction, due to internal eddy mixing. The simulation shows that resolved eddies have a strong impact on the mixed layer properties and that isopycnal eddy mixing also contributes to the generation of more homogeneous mode water characteristics in the interior.     

Quasi-stationary striations in basin-scale oceanic circulation: vorticity balance from observations and eddy-resolving model

Recently, prominent jet-like features of the ocean circulation, called striations, with a meridional scale of O(300–500 ;km) and extending for thousands of kilometers in length, have been detected in satellite and in situ observations and in high-resolution numerical models. In this paper, we study quasi-stationary striations, which are best seen in the multi-year time-averaged velocity fields. Using 1993–2002 mean dynamic ocean topography, satellite altimeter observations, and output of the Ocean General Circulation Model for the Earth Simulator, the dynamics of the quasi-stationary striations in the eastern parts of the subtropical North and South Pacific are examined by assessing individual terms in the time-averaged equations of relative and potential vorticity. While non-linear effects are found to be essential in the dynamics of the striations, rejecting some linear hypotheses forwarded in the earlier studies, the relevance of the Rhines mechanism is not confirmed. Eddy flux does not act as a relative vorticity source for the striations. Using the potential vorticity (PV) diagnostics, we show that the time-mean PV is not conserved along the time-mean streamlines, and on the scale of the striations these changes in PV are largely induced by the eddy flux of layer thickness. The fact that eddy fluxes contribute to the striations’ time-mean PV budget suggests that the striations are not a kinematical artifact of time-averaging of westward-propagating eddies.