Submesoscales are ephemeral oceanic features in the form of vortices, density fronts, and filaments, with horizontal scales ranging from 0.1 to 10 km and lifetimes of hours to days, concentrating mainly in the weakly stratified surface layer and the bottom boundary layer (BBL). Unlike mesoscale dynamics, where vertical stratification dominates and planetary rotation firmly constrains the flow, submesoscale structures are still influenced but only weakly constrained by the Earth's rotation . Dynamically, they can thus be defined as flows for which the Rossby and balanced Richardson numbers, 1Ro=ζf,Rib=f2N2|∇hb|2, are both O(1). Here, f denotes the planetary vorticity, ζ=vx-uy the vertical component of the relative vorticity, N2=∂b/∂z the square of the buoyancy frequency, where b=-g(ρ-ρ0)/ρ0 is the buoyancy based on the (potential) density ρ (g is the gravitational acceleration, and ρ0 a constant reference density), and ∇h denotes the horizontal nabla operator.

While surface submesoscales have been extensively studied and different mechanisms have been proposed for their emergence, their counterparts in the BBL remain little explored. However, recent modeling efforts , supported by scarce observational data , suggest that submesoscales are abundant in the oceanic interior, challenging the traditional view of a quiescent, stratified interior in quasigeostrophic equilibrium . These ageostrophic features, which primarily emerge from the interaction of geostrophic flows with steep topographic slopes (as explained below), are dynamically important, as they provide a direct pathway for the forward cascade of mesoscale energy toward dissipation scales .

Specifically, showed that when an interior current flows above a sloping bottom, the drag against the slope will create not only vertical but also significant horizontal velocity shear, i.e., vertical vorticity. Currents with large vertical vorticity may become unstable, separate from the slope, especially at locations of strong topographic curvature, and form turbulent topographic wakes. In the separated wake, submesoscale instabilities typically emerge, resulting in enhanced energy dissipation and the formation of submesoscale coherent vortices . This process has been shown, for example, in the context of the California Undercurrent , the Gulf Stream , and the Kuroshio–Oyashio Extension . The resulting submesoscale coherent vortices are usually long-lived, which enables them to transport their anomalous physical and biogeochemical properties over long distances away from the generation sites . More recently, observations from a deep western boundary current in the Southern Ocean have further highlighted the role of submesoscale instabilities in driving boundary–interior exchange and turbulent boundary mixing .

The aforementioned studies have explored the emergence of interior submesoscales by focusing particularly on the vicinity of strong current systems , coastal jets and dense water outflows . Here the focus is on wind-driven, non-tidal marine and inland water systems that do not include western boundary currents or other large-scale current systems. Examples include semi-enclosed marine systems like the Baltic Sea, the Black Sea, and large limnic systems like the Great Lakes in North America and the Caspian Sea. Our goal is to investigate how wind events influence the deep-water submesoscale dynamics in such wind-driven systems, complementing the study of , which focused on the effects of storms on surface-layer submesoscale dynamics. We therefore use a realistic high-resolution numerical simulation to study the emergence of these features during sequential storm events.

Our study area is the Baltic Sea, a shallow and strongly stratified brackish marine system connected to the North Sea through the Danish Straits (Fig. a). Persistent lateral density gradients and rich submesoscale activity dominate in the surface mixed layer , whereas mesoscale eddies are found at various depths, concentrating in the pycnocline . Temporal variability of currents is high in response to variable atmospheric forcing . Wind forcing plays an important role in driving the circulation, influencing not only the surface but also the deeper layers of the basin .

The response of a stratified, elongated basin to constant alongshore wind forcing was described for the Baltic Sea by and by . It involves cross-shore Ekman transport in the surface layer, leading to a sea-level rise along one coast and a corresponding sea-level fall along the opposite coast. This induces downwelling and upwelling, respectively, and associated baroclinic responses at both coasts. Consequently, coastal jets develop along both coasts in the direction of the wind, accompanied by a compensating return flow in the basin interior. Such bottom countercurrents, directed opposite to the wind, can arise as a response to the developed coastal jets or, in some cases, prevail as a result of open-sea wind-driven Ekman transport . Coastal upwelling/downwelling events and the formation of wind-driven boundary currents are commonly observed and can occur concurrently at the basin's eastern and western sides . The circulation therefore rapidly reacts to the wind forcing, with wind reversals followed by current reversals in the deeper layers occurring frequently in the basin . We thus hypothesize that such wind-driven changes in the direction and structure of the interior circulation can, in turn, influence the development of submesoscale features.

We focus particularly on the Eastern Gotland Basin (blue box in Fig. a), which is characterized by a permanent halocline, typically centered at around 70 m depth (Fig. b). The strong halocline acts as a natural lid, isolating the oxygen-rich surface water from the typically anoxic deeper layers of what is one of the globally largest anoxic marine systems. As tides are virtually absent and diapycnal mixing is weak, the primary process ventilating the oxygen-depleted bottom waters is short-lived intrusion events Major Baltic Inflows;, which import oxygenated waters from the North Sea. A tracer release experiment and turbulence microstructure observations have shown that basin-scale vertical deep-water mixing in the central Baltic Sea is largely driven by boundary processes, fueled by wind-driven topographic waves and, to a smaller extent, by near-inertial waves generated by transient wind forcing . At present, however, the processes governing the exchange between mixed boundary-layer waters and the weakly turbulent interior of the Baltic Sea remain poorly understood, and the potential role of subsurface submesoscale motions has not yet been investigated in detail.

Here, the Baltic Sea is used primarily as a natural laboratory to enable the investigation of interior submesoscale phenomena in a semi-enclosed, wind-driven basin, placed in the context of the more familiar open-ocean submesoscales that arise near stationary, large-scale boundary currents. Our goal is to show that in such basins, intense and variable winds can strongly influence submesoscale motions – a hypothesis that is likely to be broadly applicable to other wind-driven marine basins and large lakes. The manuscript is organized as follows. Section  introduces the model setup, the atmospheric conditions, and the analysis methods. In Sect. , we explore the co-existence of surface and subsurface submesoscales, immediately following a severe storm episode. Section  highlights the main characteristics of the wind-driven submesoscales emerging in the stratified interior, below the pycnocline. Section  investigates the formation of these features during sequential wind events by exploring the effects of wind magnitude and direction. In Sect. , we examine the different instability types emerging in the BBL during storms. The main findings of our study are summarized in Sect. .

At the peak of the last storm event (31 October 2017), Rossby number maps reveal a wealth of small-scale structures at all depth levels (Fig. ). Features with enhanced Ro appear not only in the surface layer (Fig. a) but also in the deeper parts of the water column (Fig. b and c), especially close to the sloping boundaries. Inside these features, we find widespread regions with Ro=O(1), suggesting that submesoscale processes form an important component of the deep-water dynamics also in this shallow wind-driven marine system.

Focusing on the surface layer, numerous fronts and filaments can be identified from the spatial distribution of the near-surface Rossby number (Fig. a). showed that these features are characterized by shallow mixed layers even during periods of strong wind forcing, due to submesoscale restratification. Their emergence was shown to be linked to a large-scale persistent frontal structure which, during autumn, occupies the entire eastern part of the basin. Typically, the surface submesoscale phenomena stay confined to the weakly stratified surface mixed layer, and therefore have no signature in the deeper layers .

Figure b, c show that submesoscale features with Ro=O(1) are ubiquitous also in the stratified interior region below the permanent pycnocline. These features are not spatially correlated to their counterparts in the surface layer, suggesting an independent generation mechanism. Significantly higher values of Ro close to the lateral boundaries, with some ageostrophic patterns extending further offshore, away from the sloping bottom, provide an initial indication of the relevance of submesoscale generation sites near the boundaries. A prominent example is the anticyclonic vortex with a strongly negative vorticity core that develops during the northerly storm (marked by the black arrow in Fig. b). Two eddy-like features with positive Ro=O(1) are also visible, particularly at 100 m depth (black arrows in Fig. c). These features were formed during the first wind episode (see also Fig.  below) and have a boundary origin, potentially arising from flow-topography interactions, as noted in previous studies .

Overall, small-scale structures with enhanced Ro=O(1) are ubiquitous at the peak of the storm, both at the surface and in the deeper layers. However, the surface and the interior regions exhibit distinct and phenomenologically disconnected patterns, especially near the lateral margins. This suggests that surface and subsurface submesoscale motions are able to coexist while remaining dynamically isolated from each other. This contrasts with previous studies e.g., showing that submesoscale motions observed in the surface layer were formed above sloping topography, and thus had an interior origin. Building on existing understanding, our investigations will therefore concentrate on the emergence of submesoscales in the deeper layers, below the main pycnocline, focusing particularly on the role of wind events in modulating these phenomena in a predominantly wind-driven, non-tidal system.

To explore the spatial distribution of the flow after the first wind event, horizontal slices of various dynamical quantities at 80 m depth are shown in Fig. , where the signatures of fronts, filaments, and spiral-like features are clearly visible. Close to the lateral margins, particularly large horizontal buoyancy gradients (Fig. c), often combined with Ro=O(1) (Fig. e), indicate submesoscale dynamics. The pronounced |∇hb|2 patterns also have a signature in the temperature field, where boluses of colder water can be spotted throughout the basin (Fig. a). Figure b shows the spatial distribution of the passive tracer, released at 80 m depth exclusively within the BBL 16 d prior to the study period (Sect. ). The tracer shows distinct features that correlate strongly with those observed in the buoyancy gradients. These features are ageostrophic, as evident from the large Ro values (Fig. e), and seemingly control the lateral dispersion of the tracer by initially trapping it and then transporting it away from the boundary, into the stratified interior.

A spatial discrepancy between the western and eastern margins of the basin is visible from Fig. e and f. In the west, thin stripes of negative vorticity appear that tend to follow isobaths, whereas in the east, positive vorticity signals dominate. The negative relative vorticity (or Ro<0, Fig. e) along the slope east of the island of Gotland is generated by lateral topographic drag on the southward flow (black arrows in Fig. b). Similarly, regions with large Ro>0 develop along the Latvian coastline as the current flows with the coast to its left, such that topographic drag now generates cyclonic shear. This is consistent with numerous studies demonstrating that flow-topography interactions can lead to the generation of either anticyclonic or cyclonic relative vorticity, depending on the direction of the interior flow with respect to the coast . However, different from these studies focusing on strong large-scale currents along steep topographic slopes, Fig.  shows (as will subsequent sections) that this mechanism is also at work in this semi-enclosed basin where ephemeral wind-driven currents largely dominate. Arrow (i) in Fig. e marks a location where a boundary current with large cyclonic ζ separates from the slope, and rolls up around an anticyclonic cold vortex. The same process most likely applies to the feature observed at the southeastern margin of the basin (Fig. e, black arrow ii). The strong ζ anomalies in the form of eddies, fronts and filaments are, therefore, transported from the boundary into the interior, rather than remaining confined to the margins.

In addition to the relative vorticity generation, a sloping turbulent BBL can also act as a source or sink of potential vorticity . The former case is visible along the eastern boundary (Fig. f), where stably stratified regions with positive EPV and Ro can clearly be identified. As EPV is conserved in the largely non-turbulent interior region, it serves as a tracer for BBL fluid that has detached from the basin slopes. One example is the positive EPV anomaly characterizing the high-Ro filament detaching from the boundary at the location marked by arrow (i) in Fig. e. Another prominent example is the structure marked by the blue box in Fig. f, which is characterized by high PV in its center. The near-zero PV at the feature's periphery results from the negative ζ. The structure, which can be seen 5 d later in the embedded panel of Fig. f, detaches from the boundary and carries anomalous PV further offshore. Subsequently, the positive EPV anomaly propagates northwards, and can be spotted during the last storm event from the Ro map shown in Fig. c (black arrow i). Overall, the eddies and filaments observed here and in Fig.  visibly originate from the lateral boundaries, potentially arising through BBL instabilities and flow detachment. They are characterized by enhanced horizontal buoyancy gradients and are associated with low and high vorticity anomalies.

On the other side of the basin, along the topographic slope east of the island of Gotland, the prevalence of negative EPV values indicates that the region is prone to instabilities (see Sect. ). The EPV < 0 stripes are confined to the boundary region with Ro<0, suggesting that centrifugal instability may occur there. However, as will be shown below, centrifugal instability, which arises under stable stratification (N2>0) when the EPV is vortically negative, is not abundant in our simulations. This is probably due to the lack of steep bathymetric slopes, which have been seen to favor its emergence in other situations .

To investigate the potential impact of storms on the BBL structure, vertical transects showing pre-storm flow conditions are contrasted to storm conditions in Fig. . The snapshots in Fig.  correspond to (see Fig. c): a period with relatively weak winds before the wind events in the last week of October (left panels), the onset of the last storm event (middle panel), and its final phase (right panel). During the storm event, pronounced northerly winds induced coastal upwelling in the eastern part of the basin.

Before the northerly storm (left panels in Fig. ), the interior flow below the mixed layer is directed southwards. Negative Ro values thus appear at the western boundary, and positive at the eastern boundary (Fig. g). Such opposing signs are also found in the EPV, which exhibits highly correlated patterns with Ro (see also Fig. ). Considering that the strong halocline effectively separates the surface-layer and deep-water processes, the vorticity seen here and the associated submesoscale structures in Fig.  are most likely generated by flow-topography interactions rather than surface processes. During these pre-storm conditions with winds from south-easterly direction, the dominance of EPV < 0, indicative of BBL overturning instabilities, implies that the western boundary is a pronounced mixing hotspot (Fig. j).

Driven by mostly northerly winds during the storm, however, the direction of the along-slope flow reverses, significantly affecting the vorticity sources and sinks along with the mixing hotspots. Alternating positive/negative BBL vorticity signals, therefore, appear close to the bottom in the snapshots shown in Fig. , with the signals intensifying during the severe storm episode. These alternating signals (both in time and between east and west) result from the variability of the along-slope flow, and are thus indirectly linked to the transient nature of the wind field, as shown also further below. At the onset of the storm (Fig. e), the strong northerly winds drive a southward flow in the surface layer, inducing a baroclinic return current in the deeper layers. The roughly two-layer flow in the interior and the emergence of coastal jets aligned with the wind direction at both coasts suggest that the flow exhibits features indicative of Ekman dynamics (see Introduction), as described by . Close to the eastern margin, this roughly two-layer flow (more pronounced at the peak of the wind event) generates vorticity with opposite signs above and below the mixed layer (Fig. h, around 20–21° E). Thus, opposing vertical vorticities are found not only between the eastern and western sides of the basin, but also between the near-surface and deeper layers below the pycnocline.

The opposing vertical vorticity signals between the eastern and western boundaries are more clearly illustrated in Fig. , which focuses on the near-boundary region. During the first half of the study period, prior to the northerly storm, the dominance of southward along-slope flow (Fig. a) results in positive Ro at the eastern boundary, while vorticity of opposite sign develops along the western boundary (Fig. b). During the prolonged final storm event, however, the along-slope currents reverse direction, accompanied by a corresponding reversal in the sign of the vorticity, consistent with the snapshots shown in Figs. d–i. Such current reversals following wind reversals in the deeper layers have been reported in previous studies, indicating that the circulation rapidly responds to the atmospheric forcing . Note that during the second storm event – characterized by south-southeasterly winds – a predominantly southward interior flow would be expected; however, northward currents are observed along the eastern boundary, likely reflecting the greater wind variability during this event compared to the more persistent northerly winds associated with the final storm. Wind reversals do not necessarily result in an immediate reversal of the current, as the pre-existing flow first needs to adjust to the new wind forcing. The response of Ro therefore depends on the pre-existing flow conditions but also on wind intensity and storm duration. Near-inertial waves are also evident (Fig. b), with previous studies demonstrating that in our study region they play a relatively minor role in boundary mixing (see Introduction).

Figure c further shows the fraction of the near-bottom region occupied by nearly zero and negative EPV. The percentage of low PV, indicative of conditions favorable for overturning instabilities , increases substantially at the eastern boundary during storm periods, reaching nearly 30 % at the peak of the final event, while simultaneously decreasing along the western boundary. Such negative PV values can arise under stable stratification when strong anticyclonic vertical vorticity dominates . Here, the occurrence of negative PV appears to be primarily associated with the direction of the along-slope flow, which, as discussed below, is favorably aligned to trigger downslope Ekman transport, thereby reducing stratification and PV. Overall, Figs.  and show that wind reversals frequently followed by current reversals see also can lead to alternate ζ and EPV signals within the BBL, suggesting that the submesoscale generation sites and mixing hotspots exhibit transient behavior. These signals are observed to alternate between the eastern and western sides of the basin and intensify locally during storms.

While Figs.  and mainly show that storms are important in regulating the alternating Ro patterns by controlling the direction of the along-slope currents, Fig.  more clearly indicates that storms also influence the magnitude of Ro. The time series display the absolute values of Ro (|Ro|) at 80 m depth, in the interior (depth > 100 m), as well as along the lateral boundaries of the basin (Fig. a). Consistent with Fig. , Fig.  exhibits higher |Ro| values during the severe wind events, suggesting that submesoscale structures are generated. When only the strong |Ro| values (|Ro| > 0.5, Fig. b) are considered, this storm-driven intensification becomes more evident. |Ro| is higher at the lateral boundaries, further supporting that these areas act as submesoscale generation sites, with the decay time scale for Ro appearing to be of the order of 1 week (see Fig.  in the Appendix). Low |Ro| is observed after the second storm event (Fig. b), when winds weakened and shifted from southward to northward (Fig. c and Fig. ). This is partly because the winds should first act to decelerate the pre-existing flow before the current can reverse, suggesting that the response to individual wind events also depends on the pre-existing flow conditions. Figure  further shows that the |Ro| intensity largely follows the wind variability, albeit with a time lag, especially after wind reversals (e.g., last storm event). By considering only the strong Ro (|Ro| > 0.5), we illustrate in Fig. c the fraction of the boundary occupied by such strong ageostrophic features. This fraction clearly increases during storm periods, reaching approximately 27 % during the final storm episode.

Focusing on the eastern boundary (red region in Fig. a) cross-correlation analysis between the along-slope velocity and Ro shows a maximum correlation at zero lag (r = -0.84), confirming that Ro is tightly linked to the current. Consistent with Fig. , the negative correlation indicates that the sign and magnitude of Ro follow the direction and strength of the current, becoming increasingly negative as the northward current strengthens and positive when the current is southward (see Appendix, Fig. ). Moreover, a maximum correlation of r=0.6 at a lag of approximately 13 h was found between Ro at the eastern boundary and the meridional wind component, suggesting that changes in wind forcing are followed by a response in Ro on a timescale comparable to the local inertial period. Note that we use the meridional wind component here because the storm events are predominantly oriented in the north–south direction, and the meridional component is approximately aligned with the alongshore direction at the eastern boundary. Focusing on the entire Eastern Gotland Basin, the wind field and the interior current were strongly correlated (r=0.69) at a lag of 14 h. Similar results were obtained when correlating the current magnitude specifically in the boundary region (blue shading in the inset panel of Fig. a) with the magnitude of the meridional wind component (r=0.8, lag = 16 h) and the magnitude of the meridional wind stress (r=0.82, lag = 15 h). Overall, our analyses indicate that storms affect the magnitude and, most importantly, the sign of vorticity – by reversing the interior currents. This reveals a clear dynamical link between wind forcing and near-bottom flow, providing new insights into the dynamics of submesoscale phenomena in wind-driven basins.

Regions with negative EPV in a sloping BBL indicate that the flow is prone to overturning instabilities. To explore the temporal evolution of the various instability types (see Sect. ), namely gravitational (GI), symmetric (SI), and centrifugal instability (CI), we focus specifically on the eastern boundary of the basin where, as seen in Fig. c, the EPV < 0 regions increase substantially during storms. Figure a illustrates the fraction of the near-bottom region that is unstable to each instability type. SI (blue line) clearly dominates throughout the entire study period, with the GI (green line) percentages increasing during the two severe wind episodes. The fraction of the domain that satisfies the SI criteria increases significantly during and after the last storm event, reaching approximately 25 %. As briefly discussed below, this effect can be largely traced to a wind-induced reversal of the along-slope currents northward (Fig. a, red line), which are downwelling-favorable (i.e., downslope Ekman transport) and therefore support the evolution of SI. The fraction of the domain satisfying the CI conditions remains small (yellow line in Fig. a). This may be due to the lack of steep bathymetry, which has been seen to favor CI emergence . It may also be due to the method used here to distinguish between the different instability types. As shown in Sect. , the applied instability criteria are rather crude and do not consider mixed modes. In the BBL, however, mixed centrifugal-symmetric modes have been seen to dominate . These instabilities trigger vigorous turbulent mixing in the BBL, with the negative PV values typically coinciding with high energy dissipation . Here, enhanced energy dissipation rates are observed along the basin boundary (Fig. b), and they intensify substantially during storms, pointing to a dynamical association with the low-PV regions.

We have used a realistic high-resolution simulation to explore the genesis of interior submesoscale phenomena in a semi-enclosed, predominantly wind-driven basin and to elucidate the impact of storm events on deep-water submesoscale dynamics. Focusing on the central basin of the Baltic Sea, we find that a wealth of filamentary and eddy-like features populate the basin's interior, below the main pycnocline. To the best of our knowledge, these interior structures have not been investigated previously in the Baltic Sea, with existing studies e.g., focusing primarily on submesoscale activity in the surface mixed layer. Our results show that surface and subsurface submesoscales can coexist while remaining largely decoupled from each other. As storm-intensified currents interact with the sloping bottom, anomalies in relative vorticity and EPV are generated through frictional and diabatic processes within the BBL . These strong vorticity anomalies are subsequently exported from the boundary into the stratified interior in the form of eddies, fronts and filaments, which are inherently submesoscales and act to facilitate the boundary-interior exchanges. The latter provides a missing link connecting the observation of enhanced boundary mixing in the Baltic Sea with the finding that boundary mixing determines basin-scale deep-water mixing, which implies a connection between the BBL and the interior region .

The most energetic submesoscale features appear during storms, with wind forcing indirectly modulating these motions by altering the interior circulation. Strong correlations are found between wind forcing and deep-water circulation, as well as between Ro and along-slope flow. Such along-slope flows may be either directly wind-driven, as described by , or indirectly mediated by near-inertial and topographic waves , or different types of seiching motions . By analyzing a series of sequential storm events with opposing winds, we also found that variability in wind direction is important in determining both the submesoscale generation sites and the mixing hotspots, which in turn exhibit transient behavior. Wind reversals reverse the currents, generating oppositely-signed ζ and EPV in the BBL, with EPV < 0 indicating conditions favorable to overturning instabilities . In addition to wake instabilities arising in the lee of topography due to flow–topography interactions , such overturning instabilities can also develop within the BBL due to cross-slope Ekman transport induced by along-slope currents . These instabilities typically occur when the along-slope flow is favorably aligned (i.e. with the coast on its right in the Northern Hemisphere) to trigger downslope Ekman transport, such that lighter fluid is advected beneath denser fluid, leading to boundary layer thickening, enhanced mixing, and reduction of PV . This is consistent with Figs. –, which show enhanced PV destruction, particularly when the storm-intensified along-boundary flow is downwelling-favorable.

One of the key findings of our study is therefore that intense and variable winds can strongly influence interior submesoscale dynamics. Storms not only modulate the strength of submesoscales by altering the intensity of along-slope currents, but they also generate alternating BBL mixing hotspots (both in time and between the eastern and western sides of the basin) by continuously shifting the direction of the interior currents. Such wind-induced current reversals , leading to alternating mixing hotspots, are likely to occur in other marine and limnic systems where circulation is primarily driven by transient wind events, but they are generally not expected in major current systems where the flow direction is predominantly steady. Nevertheless, the main dynamical mechanisms described here are broadly consistent with previous studies that examined similar submesoscale features, albeit in the open ocean, particularly in regions where energetic, steady currents flow over steep topography .

The Baltic Sea, largely unexplored in terms of interior submesoscale dynamics, provides a natural laboratory for studying these processes in a semi-enclosed, wind-driven basin. Most of the dynamical mechanisms identified here are not unique to the Baltic, suggesting that our findings may be broadly applicable to other semi-enclosed, wind-driven marine basins, such as the Black Sea and Caspian Sea, as well as large limnic systems like the Great Lakes and Lake Geneva. In Lake Geneva, submesoscale motions have recently been observed in the surface layer , while strong boundary mixing – potentially linked to submesoscale instabilities – has been detected in the interior by glider measurements .

For the Baltic Sea, future work could assess the influence of the slope Burger number (Bu) in modulating the occurrence of negative EPV along the basin's boundaries. As Bu in the Baltic Sea varies broadly and often abruptly, idealized modelling work e.g.,, including sensitivity experiments with different slopes and stratification conditions, could shed light on the nature of instabilities that emerge as wind-driven flows interact with sloping boundaries.

In conclusion, our work shows that variable wind forcing over sloping basin boundaries can generate a rich field of submesoscale motions via a breadth of dynamical pathways. These pathways include many canonical scenarios such as topographic wakes and instabilities within the BBL considered by past theoretical and idealised modeling studies, while extending the range of topographic and ambient flow conditions in which submesoscales may develop. This result suggests that intense and variable wind forcing over sloping topography can be particularly effective in promoting near-boundary mixing, energy dissipation, and boundary-interior exchange, and calls for an investigation of our findings' generality to other transient wind-driven basins.