Coriolis Parameter Research Articles

The importance of eddy stirring in wind-driven coastal upwelling

Abstract Coastal upwelling systems play a key role in sustaining productive coastal ecosystems in the global ocean by transporting nutrients to surface waters. However, the fundamental mechanisms and pathways responsible for nutrient upwelling are not fully understood, largely due to the historically employed two-dimensional frameworks in which coastal upwelling systems have long been studied. Using both observations and idealized numerical simulations, we identify and quantify two primary routes of nutrient upwelling: the residual circulation, resulting from a significant cancellation between Eulerian mean and eddy-induced circulations, and along-isopycnal eddy stirring. Our analysis demonstrates that their relative contributions depend on two distinct parameters: 1) the slope Burger number S, defined here as S = αN/f, where α is the topographic slope angle and N and f are the buoyancy and Coriolis frequencies and 2) the surface nutrient uptake rate by biological activities. Specifically, we propose that wind forcing induces isopycnal tilting and surface outcropping, which creates favorable conditions for along-isopycnal nutrient gradients to develop in regions of strong biological activity at the surface. The magnitude of these gradients depends on both the slope Burger number (S), which influences the strength of the residual circulation bringing nutrients from depths, and the surface biological uptake rate, which consumes nutrients. Our diagnostics provide insights into the intricate pathways for nutrient upwelling and underscore the significance of eddy stirring in coastal upwelling systems.

Theoretical Framework for the Response of Barotropic Shelf Current to Along-Shelf Bathymetric Perturbations

Abstract The mechanism underlying the response of barotropic shelf current to along-shelf bathymetric perturbations is rationalized by a one-dimensional vorticity equation. Fourier analysis of the equation in the wavenumber domain suggests that bathymetric features with relatively short length scales in the alongshore direction drive steady cross-shelf velocity undulations when the shelf circulation is against the direction of Kelvin wave propagation (upwelling regime). In contrast, much weaker and nonoscillatory cross-shelf velocity is expected when the shelf circulation is in alignment with propagating Kelvin waves (downwelling regime). The velocity undulations in the upwelling regime result from a positive feedback mechanism associated with potential vorticity (PV) advection; however, the presence of bottom friction restricts such undulations within a few wavelengths. The characteristic wavelength of undulations is determined by when the normalized drag coefficient , where u is the along-shelf mean flow speed, Hs is the shelf depth, αs is the bottom slope, f is the Coriolis parameter, and r is the linear drag coefficient. Using this formula, a scale of the cross-shelf displacement distance Yt is proposed, aligning with one-layer barotropic simulations over a large parameter space. The study provides a framework for investigating barotropic cross-shelf exchange caused by shelf circulation over irregular bottom bathymetries.

Effect of rotation on wake vortices in stratified flow

Stratified wakes past an isolated conical seamount are simulated at a Froude number of $Fr = 0.15$ and Rossby numbers of $Ro = 0.15$ , 0.75 and $\infty$ . The wakes exhibit a Kármán vortex street, unlike their unstratified, non-rotating counterpart. Vortex structures are studied in terms of large-scale global modes, as well as spatially localised vortex evolution, with a focus on rotation effects. The global modes are extracted by spectral proper orthogonal decomposition (SPOD). For all three studied $Ro$ ranging from mesoscale, submesoscale and non-rotating cases, the frequency of the SPOD modes at different heights remains coupled as a global constant. However, the shape of the SPOD modes changes from slanted ‘tongues’ at zero rotation ( $Ro=\infty$ ) to tall hill-height columns at strong rotation ( $Ro=0.15$ ). A novel method for vortex centre tracking shows that, in all three cases, the vortices at different heights advect uniformly at approximately $0.9U_{\infty }$ beyond the near wake, consistent with the lack of variability of the global modes. Under system rotation, cyclonic vortices and anticyclonic vortices (AVs) present considerable asymmetry, especially at $Ro = 0.75$ . The vorticity distribution as well as the stability of AVs are tracked downstream using statistics conditioned to the identified vortex centres. At $Ro=0.75$ , intense AVs with relative vorticity up to $\omega _z/f_{c}=-2.4$ (where $\omega_z$ is the vertical vorticity and $f_c$ is the Coriolis frequency) are seen with small regions of instability and they maintain large $\omega _z/f_{c}$ magnitude in the far wake. Recent stability analysis that accounts for stratification and viscosity is found to improve on earlier criteria and show that these intense AVs are stable.

Modeling the Tilt of Bend‐Traversing Turbidity Currents: Implications for Sinuous Submarine Channel Development

AbstractThe controls on the development of submarine channel sinuosity are contested: slope gradient and Coriolis forcing have both been recognized as key governing factors: gradient via an inverse relationship (low sinuosity at high slope and vice versa), and Coriolis forcing through its effect on sedimentation patterns (reducing lateral bend migration, and hence sinuosity development, at high latitudes and/or in large channels). Using theoretical models to calculate the bulk properties of channelized turbidity currents, this study investigates the joint role of the Coriolis force and parameters including channel size, downchannel slope and turbidity current properties in the development of submarine channel sinuosity. Model validation is undertaken through the comparison of the calculated turbidity current tilting against the measured tilting of channel levees in the Northwest Atlantic Mid‐Ocean Channel; this approach is then used to evaluate the controls on channel sinuosity in nine other modern seafloor channels. The results indicate that the Coriolis force only becomes significant when the size of the channel, the slope gradient and flow conditions are within appropriate ranges instead of solely being dependent on latitude. Thus, thick and dense (≥1% bulk sediment concentration) flows traveling within steep‐gradient, small‐scale channels were shown to be relatively less susceptible to flow modification by Coriolis forcing even at high latitudes. On the other hand, thin and dilute (≪1% bulk sediment concentration) flows in shallow‐gradient, large‐scale channels showed susceptibility to Coriolis forcing at all latitudes. These results offer new insights into submarine channel evolution and intra‐channel sedimentation patterns.

A Framework for Assessing Ocean Mixed Layer Depth Evolution

AbstractThe ocean surface mixed layer plays a crucial role as an entry or exit point for heat, salt, momentum, and nutrients from the surface to the deep ocean. In this study, we introduce a framework to assess the evolution of the mixed layer depth (MLD) for realistic forcings and preconditioning conditions. Our approach involves a physically‐based parameter space defined by three dimensionless numbers: λs representing the relative contribution of the buoyancy flux and the wind stress at the air‐sea interface, Rh the Richardson number which characterizes the stability of the water column relative to the wind shear, and f/Nh which characterizes the importance of the Earth's rotation (ratio of the Coriolis frequency f and the pycnocline stratification Nh). Four MLD evolution regimes (“restratification,” “stable,” “deepening,” and “strong deepening”) are defined based on the values of the normalized temporal evolution of the MLD. We evaluate the 3D parameter space in the context of 1D simulations and we find that considering only the two dimensions (λs, Rh) is the best choice of 2D projection of this 3D parameter space. We then demonstrate the utility of this two‐dimensional λs − Rh parameter space to compare 3D realistic ocean simulations: we discuss the impact of the horizontal resolution (1°, 1/12°, or 1/60°) and the Gent‐McWilliams parameterization on MLD evolution regimes. Finally, a proof of concept of using observational data as a truth indicates how the parameter space could be used for model calibration.

Analyzing coherent structures in the tropical cyclone boundary layer using large eddy simulations

Turbulence within the tropical cyclone boundary layer plays a crucial role in the exchange of heat, moisture, and momentum between the surface and the atmosphere. This study investigates the characteristics of coherent structures, specifically streaks and rolls, using large eddy simulations. Our results highlight significant differences across the three radius cases, with smaller radius exhibiting more intense and organized turbulence and streak/roll structures. Our analyses reveal that thermodynamic conditions significantly impact the timing of initial streak/roll development but do not affect their intensity in the steady state. Wind structures closer to the tropical cyclone center lead to stronger and more rapidly developing streaks/rolls, indicating their critical role in determining the intensity and formation of these features. Sensitivity tests on the Coriolis parameter (f) and radial decay parameter of tangential wind (n) show minimal impact on the development of streaks/rolls, suggesting these factors are less influential compared to wind and thermodynamic conditions.

Recirculation through western boundary currents varies nonlinearly with the ocean basin's aspect ratio

Recirculation gyres adjacent to western boundary currents (WBCs) in the ocean enhance the poleward transport of these currents. While it is well-established that the WBC in a barotropic ocean strengthens with increase in basin's aspect ratio (the meridional-to-zonal extent ratio), how intensity of the recirculation through the western boundary layer varies with this parameter remains unexplored. I address this using the non-dimensional form of the nonlinear, wind-driven Stommel–Munk model of westward intensification that comprises three parameters—the aspect ratio (δ), the damping coefficient (ϵ), and the β-Rossby number (Rβ). Here, ϵ is set by the ratio of Rayleigh friction coefficient (or eddy viscosity) to the meridional gradient of the Coriolis frequency and the basin's zonal dimension, while Rβ is proportional to wind stress amplitude and quantifies the strength of nonlinearity. In the weak-to-moderate nonlinearity limit (Rβ<∼ϵ), perturbation analysis reveals that recirculation varies concavely with aspect ratio, suggesting existence of an optimal aspect ratio (δopt) for which the recirculation is maximum and for typical values of ϵ (10−3−10−2), δopt follows the power-law relation δopt=4.3ϵ. Numerical simulations further validate the existence of δopt. For large ϵ (>5×10−3), the power-law predicts δopt for the numerical solutions rather accurately, but does not hold for smaller ϵ (2×10−3) due to increased importance of nonlinear terms. Nevertheless, the nonlinear variation in recirculation through the western boundary layer with aspect ratio is observed for all ϵ values and may contribute to the heterogeneous increase in the WBC's transport across different ocean basins in a warming climate.

The Coastal Sea‐Level Response to Wind Stress in the Middle Atlantic Bight

AbstractAnalysis of 40 years of tide gauge data and reanalysis wind stresses from the Middle Atlantic Bight (MAB) indicate that along‐shelf wind stresses are a dominant driver of coastal dynamic sea level (sea level plus atmospheric pressure) variability at daily to yearly time scales. The sea‐level response to along‐shelf wind stress varies substantially along the coast and is accurately reproduced by a steady, barotropic, depth‐averaged model (Csanady, 1978, https://doi.org/10.1175/1520‐0485(1978)008<0047:tatw>2.0.co;2, Arrested Topographic Wave). The model indicates that the sea‐level response in the MAB depends primarily on the along‐shelf distribution of the along‐shelf wind stress, the Coriolis frequency, the bottom drag coefficient, and the cross‐shelf bottom slope. The along‐shelf wind stress varies along the MAB shelf due primarily to changes in the shelf orientation. The sea‐level response depends on both the local and upstream (in the sense of Kelvin wave propagation) along‐shelf wind stresses. Consequently, sea‐level variability at daily, monthly and yearly time scales along much of the central MAB coast is more strongly driven by upstream winds along the southern New England shelf than by local winds along the central MAB shelf. The residual coastal sea‐level variability, after removing the wind‐driven response and the trend, is roughly uniform along the MAB coast. The along‐coast average of the residual sea level at monthly and yearly time scales is caused by variations in shelf water densities primarily associated with the large annual cycle in water temperature and interannual variations in salinity.

Cross-Shelf Exchange in Prograde Antarctic Troughs Driven by Offshore-Propagating Dense Water Eddies

Abstract This study examines the link between near-bottom outflows of dense water formed in Antarctic coastal polynyas and onshore intrusions of Circumpolar Deep Water (CDW) through prograde troughs cutting across the continental shelf. Numerical simulations show that the dense water outflow is primarily in the form of cyclonic eddies. The trough serves as a topographic guide that organizes the offshore-moving dense water eddies into a chain pattern. The offshore migration speed of the dense water eddies is similar to the velocity of the dense water offshore flow in the trough, which scaling analysis finds to be proportional to the reduced gravity of the dense water and the slope of the trough sidewalls and to be inversely proportional to the Coriolis parameter. Our model simulations indicate that, as these cyclonic dense water eddies move across the trough mouth into the deep ocean, they entrain CDW from offshore and carry CDW clockwise along their periphery into the trough. Subsequent cyclonic dense water eddies then entrain the intruding CDW further toward the coast along the trough. This process of recurring onshore entrainment of CDW by a topographically constrained chain of offshore-flowing dense water eddies is consistent with topographic hotspots of onshore intrusion of CDW around Antarctica identified by other studies. It can bring CDW from offshore to close to the coast and thus impact the heat flux into Antarctic coastal regions, affecting interactions among ocean, sea ice, and ice shelves. Significance Statement Troughs cutting across the Antarctic continental shelf are a major conduit for the transport of dense shelf water from coastal formation regions to the shelf break. This study describes a process in which clockwise-spinning eddies moving offshore in prograde troughs successively entrain filaments of relatively warm Circumpolar Deep Water from offshore across the entire shelf and into the coastal region. This eddy-induced transport provides a new understanding of the shelf edge exchange process identified in previous studies and a mechanism for further onshore intrusion of the warm Circumpolar Deep Water over parts of the Antarctic shelf. The resultant onshore heat flux could potentially bring a substantial amount of heat from offshore into the coastal region and thus affect ice–ocean interactions through melting sea ice and ice shelves.

How Atmospheric Forcing Frequency, Horizontal and Vertical Grid Resolutions Impact Mesoscale Eddy Evolution in a Numerical Model

AbstractSeasonal evolution of both surface signature and subsurface structure of a Mediterranean mesoscale anticyclones is assessed using the Coastal and Regional Ocean Community high‐resolution numerical model with realistic background stratification and fluxes. In good agreement with remote‐sensing and in‐situ observations, our numerical simulations capture the seasonal cycle of the anomalies induced by the anticyclone, both in the sea surface temperature (SST) and in the mixed layer depth (MLD). The eddy signature on the SST shifts from warm‐core in winter to cold‐core in summer, while the MLD deepens significantly in the core of the anticyclone in late winter. Our sensitivity analysis shows that the eddy SST anomaly can be accurately reproduced only if the vertical resolution is high enough (∼4 m in near surface) and if the atmospheric forcing contains high‐frequency. In summer with this configuration, the vertical mixing parameterized by the k − ϵ closure scheme is three times higher inside the eddy than outside the eddy, and leads to an anticyclonic cold core SST anomaly. This differential mixing is explained by near‐inertial waves, triggered by the high‐frequency atmospheric forcing. Near‐inertial waves propagate more energy inside the eddy because of the lower effective Coriolis parameter in the anticyclone core. On the other hand, eddy MLD anomaly appears more sensitive to horizontal resolution, and requires SST retroaction on air‐sea fluxes. These results detail the need of high frequency forcing, high vertical and horizontal resolutions to accurately reproduce the evolution of a mesoscale eddy.

A Scale-Aware Parameterization of Restratification Effect of Turbulent Thermal Wind Balance

Abstract The vertical buoyancy flux Bf under the turbulent thermal wind (TTW) balance plays an important role in restratifying the surface mixed layer in winter. So far, most of the global ocean models are too coarse to resolve this process. In this paper, a scale-aware parameterization is proposed for and implemented in a hierarchy of regional ocean simulations over the winter Kuroshio Extension with horizontal resolutions ranging from 27 to 1 km. The parameterization depends on the Coriolis parameter, model-simulated turbulent vertical viscosity, horizontal density gradient, and a scaling relationship to adjust for the effects of model horizontal resolution on the simulated horizontal density gradient. It shows good skills in reconciling the difference between in the coarse-resolution simulations (27, 9, and 3 km) and in the 1-km simulation where is well resolved. Furthermore, implementation of the parameterization improves the simulated stratification in the surface mixed layer in coarse-resolution simulations.

Parameterization of Shear‐To‐Strain Ratio Used in Finescale Parameterization

AbstractThe Gregg‐Henyey‐Polzin (GHP) finescale parameterization is widely employed to infer turbulent mixing when microstructure measurements are unavailable. However, the strain‐only GHP scaling is the only option when shear information is lacking, in which case the shear‐to‐strain ratio (Rω) is commonly treated as a constant. Since Rω has been reported to vary spatiotemporally, using a fixed value of Rω might result in a significant bias in inferring turbulent mixing. In this study, we present a global map of Rω, based on the microstructure database contributed by the Climate Process Team on internal wave‐driven ocean mixing, which covers the Indian, Pacific, and Atlantic oceans across both high and low latitudes. Then, we propose a parameterization of Rω by considering buoyancy frequency, Coriolis frequency, and topographic features. Compared to the GHP scaling with constant Rω, the turbulent dissipation rate values inferred from the GHP scaling with parameterized Rω are more accurate. In this manner, this study provides a reference for choosing optimal Rω when using the strain‐only GHP scaling to explore turbulent mixing in the open ocean.

The different dynamic influences of Typhoon Kalmaegi on two pre-existing anticyclonic ocean eddies

Abstract. Using multi-source observational data and GLORYS12V1 reanalysis data, we conduct a comparative analysis of different responses of two warm eddies, AE1 and AE2 in the northern South China Sea, to Typhoon Kalmaegi during September 2014. The findings of our research are as follows: (1) for horizontal distribution, the area and the sea surface temperature (SST) of AE1 and AE2 decrease by about 31 % (36 %) and 0.4 °C (0.6 °C). The amplitude, Rossby number (Ro = relative vorticity / Coriolis parameter) and eddy kinetic energy (EKE) of AE1 increase by 1.3 cm (5.7 %), 1.4×10-2 (20.6 %) and 107.2 cm2 s−2 (49.2 %) after the typhoon, respectively, while AE2 weakens and the amplitude, Rossby number and EKE decrease by 3.1 cm (14.6 %), 1.6×10-2 (26.2 %) and 38.5 cm2 s−2 (20.2 %), respectively. (2) In the vertical direction, AE1 demonstrates enhanced convergence, leading to an increase in temperature and a decrease in salinity above 150 m. The response below the mixed-layer depth (MLD) is particularly prominent (1.3 °C). In contrast, AE2 experiences cooling and a decrease in salinity above the MLD. Below the MLD, it exhibits a subsurface temperature drop and salinity increase due to the upwelling of cold water induced by the suction effect of the typhoon. (3) The disparity in the responses of the two warm eddies can be attributed to their different positions relative to Typhoon Kalmaegi. Under the influence of negative wind stress curl outside the maximum wind radius (Rmax) of the typhoon, triggering negative Ekman pumping velocity (EPV) and quasi-geostrophic adjustment of the eddy, the warm eddy AE1, with its center to the left of the typhoon's path, further enhances the converging sinking of the upper warm water, resulting in its intensification. On the other hand, the warm eddy AE2, situated closer to the center of the typhoon, weakens due to the cold suction caused by the strong positive wind stress curl within the typhoon's Rmax. The same polarity eddies may have different response to typhoons. The distance between eddies and typhoons, eddy intensity, and the background field need to be considered.

How heating tracers drive self-lofting long-lived stratospheric anticyclones: simple dynamical models

Abstract. Long-lived “bubbles” of wildfire smoke or volcanic aerosol have recently been observed in the stratosphere, co-located with ozone, carbon monoxide, and water vapour anomalies. These bubbles often survive for several weeks, during which time they ascend through vertical distances of 15 km or more. Meteorological analysis data suggest that this aerosol is contained within strong, persistent anticyclonic vortices. Absorption of solar radiation by the aerosol is hypothesised to drive the ascent of the bubbles, but the dynamics of how this heating gives rise to a single-sign anticyclonic vorticity anomaly have thus far been unclear. We present a description of heating-driven stratospheric vortices, based on an axisymmetric balanced model. The simplest version of this model includes a specified localised heating moving upwards at fixed velocity and produces a steadily translating solution with a single-signed anticyclonic vortex co-located with the heating, with corresponding temperature anomalies forming a vertical dipole, matching observations. A more complex version includes the two-way interaction between a heating tracer, representing the aerosol, and the dynamics. An evolving tracer provides heating which drives a secondary circulation, and this in turn transports the tracer. Through this two-way interaction an initial distribution of tracer drives a circulation and forms a self-lofting tracer-filled anticyclonic vortex. Scaling arguments show that upward velocity is proportional to heating magnitude, but the magnitude of peak quasigeostrophic vorticity is O(f) (f is the Coriolis parameter) and independent of the heating magnitude. Estimates of vorticity from observations match our theoretical predictions. We discuss 3-D effects such as vortex stripping and dispersion of tracer outside the vortex by the large-scale flow, which cannot be captured explicitly by the axisymmetric model and are likely to be important in the real atmosphere. The large O(f) vorticity of the fully developed anticyclones explains their observed persistence and their effective confinement of tracers. To further investigate the early stages of formation of tracer-filled vortices, we consider an idealised configuration of a homogeneous tracer layer. A linearised calculation reveals that the upper part of the layer is destabilised due to the decrease in tracer concentrations with height there, which sets up a self-reinforcing effect where upward lofting of tracer results in stronger heating and hence stronger lofting. Small amplitude disturbances form isolated tracer plumes that ascend out of the initial layer, indicative of a self-organisation of the flow. The relevance of these idealised models to formation and persistence of tracer-filled vortices in the real atmosphere is discussed, and it is suggested that a key factor in their formation is the time taken to reach the fully developed stage, which is shorter for strong heating rates.

Influences of coriolis force and friction on fluid dynamics in specific paraboloid basins

Natural water features, such as rivers, lakes, bays, and straits, exhibit diverse shapes and dynamics. Among these, paraboloid basins, resembling the contours of lakes or closed bays, offer a relatively simple yet intriguing case. This study delves into the dynamics of water waves within two distinct paraboloid basin variations: the parabolic canal and the circular paraboloid. In this investigation, we analyze shoreline movement, considering influential factors such as bottom friction and the Coriolis effect through a mathematical model. The foundation of our analysis lies in the two-dimensional shallow water equations, with the Thacker assumption serving as the basis for deriving an analytical solution. To enhance our understanding, we employ a numerical solution, employing the finite volume method on a staggered grid, to simulate wave behavior within these basin variations. For validation, we compare these numerical findings with our analytical solution. Furthermore, this study conducts a sensitivity analysis of the Coriolis parameter and bottom friction parameter under various conditions. Through this exploration, we gain valuable insights into the interplay of these critical factors in shaping the oscillation of water within paraboloid basins, enriching our understanding of coastal and basin dynamics.

The effects of curvature and β on zonally invariant f-plane dynamics

The modification of the dynamics on the f-plane due to the inclusion of the latitudinal variation of the Coriolis frequency, quantified by β, is examined by comparing steady and time-dependent solutions for different values of β, including 0. The comparisons are made for analytic expressions and numerical simulations of the solutions of two physical problems: geostrophic adjustment and Ekman transport. By varying the extent of the meridional domain, L, and the value of β, while keeping the same value of the deformation radius, Rd, we find that the relative change of the Coriolis frequency across L is not a uniformly valid predictor of the effect of β on the dynamics. Instead, the magnitude of this effect varies with the physical problem, the type of dynamics, i.e., steady or time dependent, and with the particular values of β and L. In contrast to the effect of β, the effect of curvature on the zonally invariant f-plane dynamics is negligible in all cases for L≤1200 km and Rd=30 km employed here.

Model for the cyclonic bias of convective vortices in a rotating system

In convective fluids, small vortices develop between neighboring convective cells. Familiar in the atmosphere in the form of dust devils and water spouts, these convective vortices have been seen in simulations of oceanic convection where the vortices exhibit an unexplained bias towards cyclonic rotation despite having a large Rossby number. Here we use large eddy simulations (LES) of an idealized oceanic convective mixed layer and vary independently the Coriolis acceleration and the surface buoyancy flux to investigate the development of the cyclonic bias of convective vortices. While large convective vortices are biased for sufficiently large values of the Coriolis parameter, small convective vortices do not exhibit a clear bias. Using Lagrangian particles, we find that the large convective vortices develop through the merger of many small convective vortices. We propose a statistical theory to predict the cyclonic bias of large convective vortices composed of many small unbiased convective vortices and test the theory using LES results. We apply the theory to typical convective conditions and find that convective vortices in upper ocean mixed layers are expected to exhibit a bias, while convective vortices in the terrestrial and Martian atmospheres are expected to be largely unbiased. Published by the American Physical Society 2024

Quantifying the turbulent mixing driven by the Faraday instability in rotating miscible fluids

The effect of the rotation on the turbulent mixing of two miscible fluids of small contrasting density, induced by Faraday instability, is investigated using direct numerical simulations. We quantify the irreversible mixing, which depicts the conversion of the available potential energy (APE) to the background potential energy (BPE) through the irreversible mixing rate M. We demonstrate that at lower forcing amplitudes, the turbulent kinetic energy (t.k.e.) increases with an increase in the Coriolis frequency f till (f/ω)2<0.25, where ω is the forcing frequency, during the sub-harmonic instability phase. This enhancement of t.k.e. is attributed to the excitement of more unstable modes. The irreversible mixing sustains for an extended period with increasing (f/ω)2 till 0.25 owing to the prolonged sub-harmonic instability phase and eventually ceases with instability saturation. When (f/ω)2>0.25, the Coriolis force significantly delays the onset of the sub-harmonic instabilities. The strong rotational effects result in lower turbulence because the bulk of the APE expends to BPE, decreasing APE that converts back to t.k.e. reservoir for (f/ω)2>0.25. Therefore, in the subsequent oscillation, the t.k.e. available to contribute to the external energy input from periodic forcing is small. Since the instability never saturates for (f/ω)2>0.25, the conversion of APE to BPE via M continues, and we find prolonged irreversible mixing. At higher forcing amplitudes, the instability delaying effect of rotation is negligible, and the turbulence is less intense and short-lived. Therefore, the irreversible mixing phenomenon also ends quickly for (f/ω)2<0.25. However, when (f/ω)2>0.25, a continuous irreversible mixing is observed. We also examine the mixing efficiency in terms of M and find that the mixing is efficient at lower forcing amplitudes and rotation rates of (f/ω)2>0.25 because the major portion of APE expends to BPE.

The Small-Amplitude Dynamics of Spontaneous Tropical Cyclogenesis. Part I: Experiments with Amplified Longwave Radiative Feedback

Abstract What determines the vortex size at the small-amplitude stage of spontaneous tropical cyclogenesis remains unclear. A doubly periodic domain is a standard setup for numerically studying this problem, but the convectively coupled standing waves inherent to the setup could directly trigger vortices, rendering an unrealistic path for tropical cyclogenesis. We increase the Coriolis parameter to suppress the wave and double the longwave radiative feedback to make the more realistic moisture–radiation instability dominant. Experiments show that the moisture–radiation instability has a short-wavelength cutoff due to the smoothing effect of convective dynamics, which includes the nonlocal convective triggering by cold pools and the nonlocal longwave radiative effect of anvil clouds. By approximating the spread of convective activity as a Gaussian filter on the column humidity, we derive a bulk convective spreading length lcp+av to parameterize the combined effect of cold pools and anvils. Using a novel diagnostic method, lcp+av is shown to be around 10 km. The contribution of cold pools and anvil clouds to convective spreading is comparable in the doubled radiative feedback experiments. An extrapolation to the normal radiative feedback state shows the anvil clouds play a smaller yet nonnegligible role.

Observed 10–20-Day Deep-Current Variability at 5°N, 90.5°E in the Eastern Indian Ocean

Abstract In the eastern off-equatorial Indian Ocean, deep current intraseasonal variability within a typical period of 10–20 days was revealed by a mooring at 5°N, 90.5°E, accounting for over 50% of the total bottom subtidal velocity variability. The 10–20-day oscillations were more energetic in the cross-isobathic direction (STD = 3.02 cm s−1) than those in the along-isobathic direction (STD = 1.50 cm s−1). The oscillations were interpreted as topographic Rossby waves (TRWs) because they satisfied the TRWs dispersion relation that considered the smaller Coriolis parameter and stronger β effect at low latitude. Further analysis indicated significant vertical coupling between the deep cross-slope oscillations and cross-isobathic 10–20-day perturbations at the depth of 300–950 m. The 10–20-day TRWs were generated by cross-isobathic motions under the potential vorticity conservation adjustment. The Mercator Ocean output reproduced the generation of kinetic energy (KE) of deep current variability. The associated diagnostic analysis of multiscale energetics showed that the KE of TRWs was mainly supplied by vertical pressure work. In the seamount region (2°–10°N, 89°–92°E), vertical and horizontal pressure works were identified to be the dominant energy source (contributing to 94% of the total KE source) and sink (contributing to 98% of the total KE sink) of the deep current variability, transporting energy downward and redistributing energy horizontally, respectively.