Triadic resonances in internal wave modes with background shear

Abstract In this paper, we use asymptotic theory and numerical methods to study resonant triad interactions among discrete internal wave modes at a fixed frequency ($\omega$) in a two-dimensional, uniformly stratified shear flow. Motivated by linear internal wave generation mechanisms in the ocean, we assume the primary wave field as a linear superposition of various horizontally propagating vertical modes at a fixed frequency $\omega$. The weakly nonlinear solution associated with the primary wave field is shown to comprise superharmonic (frequency $2\omega$) and zero frequency wave fields, with the focus of this study being on the former. When two interacting primary modes $m$ and $n$ are in triadic resonance with a superharmonic mode $q$, it results in the divergence of the corresponding superharmonic secondary wave amplitude. For a given modal interaction $(m, n)$, the superharmonic wave amplitude is plotted on the plane of primary wave frequency $\omega$ and Richardson number $Ri$, and the locus of divergence locations shows how the resonance locations are influenced by background shear. In the limit of weak background shear ($Ri\to \infty$), using an asymptotic theory, we show that the horizontal wavenumber condition $k_m + k_n = k_q$ is sufficient for triadic resonance, in contrast to the requirement of an additional vertical mode number condition ($q = |m-n|$) in the case of no shear. As a result, the number of resonances for an arbitrarily weak shear is significantly larger than that for no shear. The new resonances that occur in the presence of shear include the possibilities of resonance due to self-interaction and resonances that occur at the seemingly trivial limit of $\omega \approx 0$, both of which are not possible in the no shear limit. Our weak shear limit conclusions are relevant for other inhomogeneities such as non-uniformity in stratification as well, thus providing an understanding of several recent studies that have highlighted superharmonic generation in non-uniform stratifications. Extending our study to finite shear (finite $Ri$) in an ocean-like exponential shear flow profile, we show that for cograde–cograde interactions, a significant number of divergence curves that start at $Ri\to \infty$ will not extend below a cutoff $Ri$ $\sim O(1)$. In contrast, for retrograde–retrograde interactions, the divergence curves extend all the way from $Ri\to \infty$ to $Ri = 0.5$. For mixed interactions, new divergence curves appear at $\omega = 0$ for $Ri\sim O(10)$ and extend to other primary wave frequencies for smaller $Ri$. Consequently, the total ($\text {cograde} + \text {retrograde} + \text {mixed}$) number of resonant triads is of the same order for small $Ri\approx 0.5$ as in the limit of weak shear ($Ri\to \infty$), although it attains a maximum at $Ri\sim O(10)$.

In this paper, we use asymptotic theory and numerical methods to study resonant triad interactions among discrete internal wave modes at a fixed frequency (ω) in a two-dimensional, uniformly stratified shear flow. Motivated by linear internal wave generation mechanisms in the ocean, we assume the primary wave field as a linear superposition of various horizontally propagating vertical modes at a fixed frequency ω. The weakly nonlinear solution associated with the primary wave field is shown to comprise superharmonic (frequency 2ω) and zero frequency wave fields, with the focus of this study being on the former. When two interacting primary modes m and n are in triadic resonance with a superharmonic mode q, it results in the divergence of the corresponding superharmonic secondary wave amplitude. For a given modal interaction (m, n), the superharmonic wave amplitude is plotted on the plane of primary wave frequency ω and Richardson number Ri, and the locus of divergence locations shows how the resonance locations are influenced by background shear. In the limit of weak background shear (Ri → ∞), using an asymptotic theory, we show that the horizontal wavenumber condition k m + k n = k q is sufficient for triadic resonance, in contrast to the requirement of an additional vertical mode number condition (q = |m − n|) in the case of no shear. As a result, the number of resonances for an arbitrarily weak shear is significantly larger than that for no shear. The new resonances that occur in the presence of shear include the possibilities of resonance due to self-interaction and resonances that occur at the seemingly trivial limit of ω ≈ 0, both of which are not possible in the no shear limit. Our weak shear limit conclusions are relevant for other inhomogeneities such as non-uniformity in stratification as well, thus providing an understanding of several recent studies that have highlighted superharmonic generation in non-uniform stratifications. Extending our study to finite shear (finite Ri) in an ocean-like exponential shear flow profile, we show that for cograde-cograde interactions, a significant number of divergence curves that start at Ri → ∞ will not extend below a cutoff Ri ∼ O(1). In contrast, for

Introduction
Internal waves generated by tides and winds can cause intense mixing in the deep ocean (Alford 2003;Garrett 2003). Dissipation of these internal waves plays a crucial role in meridional overturning circulation (Munk & Wunsch 1998) and processes such as nutrient and plankton transport (Garrett & Munk 1979).
Triadic resonance is one of the important mechanisms leading to internal wave dissipation (Staquet & Sommeria 2002). A well-studied manifestation of triadic resonance is parametric subharmonic instability (PSI), where a primary internal wave of frequency ω and wave vector k is unstable to perturbations of two secondary waves with frequencies ω 1,2 = ω/2 and wave vectors k 1 and k 2 such that k 1 + k 2 = k. (Hasselmann 1967;Koudella & Staquet 2006). Field observations (Alford et al. 2007;MacKinnon et al. 2013), however, show much less energy transfer from internal tides to subharmonic waves than what is predicted by the theory of PSI. The effects of various realistic ocean settings such as non-uniform stratification (Young, Tsang & Balmforth 2008;Gayen & Sarkar 2013;Gururaj & Guha 2020), finite width of the wave beam (Bourget et al. 2014;Dauxois et al. 2018) and background flow (Fan & Akylas 2019 are potential reasons for the discrepancy between theory and observations of PSI. The effects of a background flow is the focus of the current study, albeit on a different manifestation of triadic resonance as described below. Another manifestation of triadic resonance occurs when two monochromatic (frequency ω) primary internal waves resonantly excite a secondary wave at superharmonic frequency 2ω. Resonant generation of superharmonic internal waves has been studied in the context of interacting internal wave beams (Teoh, Ivey & Imberger 1997;Tabaei, Akylas & Lamb 2005;Jiang & Marcus 2009). Tide-topography interaction (Lamb 2004;Korobov & Lamb 2008), internal wave beam reflection from a solid boundary (Javam, Imberger & Armfield 1999;Peacock & Tabaei 2005;Gerkema, Staquet & Bouruet-Aubertot 2006;Rodenborn et al. 2011) or a pycnocline (Thorpe 1998;Gayen & Sarkar 2013;Diamessis et al. 2014;Wunsch et al. 2014;Mercier et al. 2015) are example scenarios where interacting internal wave beams generate superharmonic internal waves. Higher harmonic generation due to surface reflection of internal tides has been observed in the ocean too (Xie et al. 2013). In a fixed-depth domain like the region between the ocean floor and surface, only discretized wavenumbers are possible and linear internal wave fields are a superposition of internal wave modes. Superharmonic generation due to modal interactions, as summarized below, has received much attention recently.
In a uniform stratification with no shear, interaction between two different internal wave modes m and n at a given frequency ω can resonantly excite a superharmonic 2ω internal wave mode |m − n| at specific values of ω if m/3 < n < 3m (Thorpe 1966). As shown in figure 1, only a fraction of the points where the horizontal resonance condition is satisfied are actual triadic resonance locations. The amplitude evolution of   Figure 1. (a) Dispersion curves (ω vs k q , grey lines) for mode numbers q = 1, 2, 3, 4, 5 in a uniform stratification with no shear. The black lines show 2ω vs k m + k n for m + n = 2 to 10, with modes m and n being at frequency ω. Horizontal resonance condition k m + k n = k q is satisfied at the points of intersection (filled and unfilled circles). Filled circles represent actual resonances. (b) Number (N H , unfilled circles) of modal interactions (between modes m and n at frequency ω, 0 < ω < 0.5) for which the horizontal resonance condition k m + k n = k q (mode q being at 2ω) is satisfied, plotted as a function of m + n. Number N R (filled circles) from the N H modal interactions that are also resonances. The symbol refers to the floor operator.
such resonant triads in a uniform stratification with no shear has recently been studied numerically (Varma, Chalamalla & Mathur 2020) and experimentally (Husseini et al. 2020). In contrast to a uniform stratification, several more triadic resonances occur in a finite-depth non-uniform stratification (Varma & Mathur 2017), including the possibility of self-interacting individual modes exciting superharmonic wave fields (Thorpe 1968;Sutherland 2016;Wunsch 2017). Varma et al. (2020) and Baker & Sutherland (2020) have studied amplitude evolution in self-interacting modes at and off resonance, respectively. In an inviscid, stratified shear flow, if the gradient Richardson number is greater than 1/4 throughout the domain, Booker & Bretherton (1967) showed that significant momentum is transferred from internal waves to the mean flow at the critical layer (where internal wave phase velocity matches the local background horizontal velocity), and strong nonlinear effects ensue. As a result, numerous studies have considered nonlinear resonant interactions near the critical layer (Brown & Stewartson 1980;Grimshaw 1988Grimshaw , 1994 in infinite-depth media. To complement the studies of Brown & Stewartson (1980, 1982a, Grimshaw (1988Grimshaw ( , 1994 derived amplitude evolution equations for primary (i.e. three interacting first-order waves) and secondary (i.e. a second-order wave interacting with two first-order waves) interactions, respectively, in a slowly varying background shear and stratification in infinite depth. While Grimshaw (1988) focused on 'weak resonance' (resonance conditions satisfied only on certain surfaces in space and time) near the critical layer, he pointed out that triadic resonance in the homogeneous flow represents the leading-order 'strong resonance' conditions in weakly inhomogeneous flow. In the specific case of stratified, anti-symmetric shear layer, Kelly (1968) analysed the second-order resonant interactions of two specific interacting singular neutral modes at constant frequency and numerically showed how the amplitude of different waves is modulated. In a finite-depth stratified shear flow, the necessary condition for an explosive interaction (i.e. finite time blow-up in the amplitude evolution) of internal wave modes is the existence of a critical layer (Becker & Grimshaw 1993;Vanneste & Vial 1994). Considering a sinusoidal background velocity profile in a uniform stratification and fixed horizontal wavenumbers (k 1 , k 2 , k 3 ) that satisfy the triadic resonance condition k 1 + k 2 + k 3 = 0, Vanneste & Vial (1994) numerically showed how different interaction coefficients vary with wave amplitudes, and thereby lead to resonance.
Realistic ocean settings include factors such as non-uniform stratification, background rotation, background shear, finite depth, excitation of a wide range of wavenumbers etc. An earlier study (Varma & Mathur 2017) has shown modal interactions, including self-interaction, can lead to resonant generation of superharmonic internal waves in a finite-depth ocean-like non-uniform stratification with background rotation. Here, we consider the effects of inhomogeneity introduced by a background flow, thus building towards a generalization of the effects of inhomogeneities on finite-depth internal wave triadic resonances. Specifically, we consider triadic resonances in a finite-depth uniform stratification in the presence of an ocean-like nonlinear background shear flow (corresponding to the wind drift layer) that monotonically increases from zero at the ocean floor to a finite value at the ocean surface. With no critical layers being present for discrete modes in a continuous shear flow (Banks, Drazin & Zaturska 1976), and motivated by forcing mechanisms being at specific frequencies in the ocean (the semi-diurnal frequency, for example), we consider triadic resonances resulting from interaction between two discrete modes at the same frequency. An analytical treatment of the weak shear limit is presented, providing insights into why inhomogeneity significantly increases the number of possible resonances. A systematic study on the effects of Richardson number, spanning weak to strong shear limits, on the occurrence of triadic resonances is then performed. Owing to the loss of symmetry about the ω = 0 axis in the dispersion curves when background flow is present, both cograde (modes that travel faster than the maximum background flow velocity) and retrograde (modes that travel slower than the minimum background flow velocity) modes are considered.
The governing equations, and weakly nonlinear solutions resulting from modal interactions, are presented in § 2.1. An analytical treatment of the weak shear limit in given in § 2.2. Section 3 discusses the results and compares the solutions from asymptotics and numerics. In § 3, a systematic study of the effects of Richardson number, including weak and strong shear limits, is presented. A brief discussion and a summary of our conclusions are provided in § 4.

Governing equations
We consider an inviscid, two-dimensional flow in a uniformly stratified fluid of depth H in the Boussinesq approximation. The base flow state is described by a stably stratified, linear density profileρ(z) and a vertically varying steady horizontal shear flow U(z)e x . The corresponding constant Brunt-Väisälä frequency is N = √ −(g/ρ * )(dρ/dz), where g = −ge z is gravitational acceleration and ρ * a reference density. The non-dimensional governing equations for the perturbation flow field are (Tsutahara 1984) All quantities (including U) in (2.1) and (2.2) are non-dimensional, using H, 1/N and ρ * as the length, time and density scales, respectively. Here, t, x and z are time, horizontal and vertical coordinates, respectively and the Jacobian is defined as J(A, B) = (∂A/∂x)(∂B/∂z) − (∂A/∂z)(∂B/∂x). The non-dimensional perturbation flow field is described by (u, w) = (−∂ψ/∂z, ∂ψ/∂x), where ψ(x, z, t) is the perturbation stream function. The buoyancy perturbation is b = −gρ/(N 2 H), where ρ(x, z, t) is the non-dimensional perturbation density. The no-penetration boundary condition is given by w(x, z = 0, t) = w(x, z = 1, t) = 0, with z = 0 and z = 1 denoting the ocean bottom and top (rigid lid), respectively. Assuming a regular perturbation expansion in , a small parameter that quantifies the relative magnitude of the nonlinear terms in the governing equations, we seek solutions for the perturbation wave field in the following form: (2. 3) The O( ) wave field (ψ 1 , b 1 ), governed by linear internal wave equations, is assumed to comprise a superposition of several modes at a frequency ω where A j , k j and c j (= ω/k j ) are the complex amplitude, horizontal wavenumber and phase velocity, respectively of mode j, with c.c. denoting complex conjugate. The vertical mode shape, φ j (z) is governed by the Taylor-Goldstein-Haurwitz equation (Kundu & Cohen 2001) along with the no-penetration boundary conditions given by φ j (z = 0) = φ j (z = 1) = 0. The wave field in (2.4) could represent a linear wave field generated by forcing at a specific frequency, internal tides generated by barotropic forcing on topography (Garrett & Kunze 2007) for example.
Here, θ j = k j (x − c j t), C mn = 2ω/(k m + k n ) and A * n is the complex conjugate of amplitude A n . The particular solution of (2.6) can now be sought in the form where h mn (z) and g mn (z) are vertical structures of the superharmonic and mean-flow terms, respectively, resulting from the interaction between modes m and n. They satisfy the following equations, obtained from (2.6): where h mn (z) = h mn (z) + h nm (z) andḡ mn (z) = g mn (z) + g nm (z). The no-penetration boundary conditions are:h mn (z = 0) =h mn (z = 1) =ḡ mn (z = 0) =ḡ mn (z = 1) = 0. The magnitude of the mean-flow term in (2.11) could potentially be influenced by a class of resonances considered by Phillips (1968), who studied the interaction between an upward and a downward propagating plane internal wave in the presence of a steady (zero frequency) shear flow with twice the vertical wavenumber of the plane waves. The focus of the current study, however, is the superharmonic part of the O( 2 ) wave field, to specifically identify the role of the background shear flow U(z). Before proceeding with a fully numerical solution of (2.12), it is instructive to analyse the asymptotic limit of weak shear.

Weak shear limit
To perform an asymptotic analysis in the weak shear limit, we define a small parameter is the Richardson number and ξ = L s /H is the ratio of shear flow length scale to the ocean depth; U s and L s are dimensional velocity and length scales of the background shear flow. In the weak shear limit (δ 1), we write where υ(z) is an O(1) function describing the vertical structure of U(z) (recall that U(z) is non-dimensional, with NH being the velocity scale used for the non-dimensionalization). We will make the reasonable assumption of L s /H < 1. This implies, in conjunction with δ 1, that the shear flow time scale L s /U s is much larger than the stratification time scale 1/N (Ri 1, in other words). We proceed to analytically solve the O( ) and O( 2 ) equations presented in § 2.1, up to O(δ 1 ), with the assumption that δ 2 δ.

The O( ) equation
We begin by seeking a regular perturbation series (in δ) solutions for the vertical mode shapes (2.5) upon substituting U(z) = δυ(z). For a fixed frequency ω, we write where φ j and k j are the mode shape and horizontal wavenumber of mode j, respectively.
At O(δ 0 ), (2.5) reduces to the linear internal wave equation in quiescent fluid with j being the mode number. In this subsection, we assume ω < 1, i.e. propagating internal waves in the zero shear limit.

The O( 2 ) equation
Similar to § 2.2.1, we seek a solution for (2.12) of the form and P mn,0 is the O(δ 0 ) term in P mn (defined in (2.9)). For two different modes (m / = n), the particular solution of (2.25) ish (2.26) . (2.27) For self-interaction (m = n), the right-hand side of (2.25) vanishes and a homogeneous solution exists if and only if (k m,0 + k n,0 ) and 2ω satisfy the linear internal wave dispersion relation with no shear. Equation (2.26) indicates thath mn,0 diverges if the denominator of I mn goes to zero, i.e. if k m + k n is the horizontal wavenumber of superharmonic mode |m − n|. In other words, modes m and n at frequency ω are in triadic resonance with mode |m − n| at frequency 2ω if its horizontal wavenumber is k m + k n . This result based on the O(δ 0 ) solution is consistent with the study of Thorpe (1966) for a uniform stratification with no shear. The requirement of the superharmonic mode number being |m − n| is the reason why only a fraction of the intersections in figure 1(a) actually represent triadic resonances in a uniform stratification. It has to be noted here that, at these resonance locations, the regular perturbation expansion with constant amplitudes breaks down, and multiple-scale analysis should be used to study the amplitude variations. Our objective in this work is to identify the resonance locations and not to study the amplitude evolution.
where A(z), B(z), C(z), D(z) and E(z) are given in Appendix B. As pointed out in the beginning of this section, we have assumed that δ 2 δ, so that O(δ 2 ) terms are larger than O( 3 ) terms in the perturbation expansion. The solution of (2.28) is (2.30) where q is an integer (γ is defined soon after (2.25)). It is noteworthy that this condition for the non-existence of a solution to (2.28) also follows from the alternative theorem for the linear differential equation of the second order (Stakgold & Holst 2011). In other words, (k m,0 + k n,0 , 2ω) satisfying the linear internal wave dispersion relation with no shear is a sufficient condition for the O( 2 ) superharmonic wave field to diverge in the presence of weak shear. Unlike the requirement for triadic resonance based on the O(δ 0 ) solution, the condition for the divergence of the O( 2 ) wave field based on the O(δ 1 ) solution does not pose any requirement on the mode number of the superharmonic internal wave, which is consistent with what is reported by Vanneste & Vial (1994). As a result, for any υ(z) considered (with υ (z) / = 0 somewhere in the domain), an important implication is that all the intersections in figure 1(a), irrespective of the superharmonic wave's mode number, represent triadic resonances if a weak shear is present. Specifically, in a uniform stratification, the 2ω vs (k m + k n ) curve has (m + n − 1)/2 intersections with the dispersion curves (figure 1(a), intersections at ω = 0 not considered; x refers to the floor operator representing the greatest integer less than or equal to x), out of which only those with q = |m − n| are triadic resonances if there is no shear (see plot of N R in figure 1b). In the presence of shear, however, all the intersections (see plot of N H in figure 1b) represent triadic resonances, which include the possibility of self-interaction (m = n) too.
A seemingly trivial limit of the horizontal resonance condition k m,0 + k n,0 = k q,0 is ω = 0. In this limit, all the frequencies (primary and superharmonic) and wavenumbers are zero, although that does not guarantee γ = qπ being satisfied. Requiring γ = qπ in the limit of ω = 0 results in the additional condition of m + n = 2q, which, when satisfied, results in the divergence ofh mn,1 (z), and hence corresponds to triadic resonance at ω = 0. As a result of these additional resonances at ω = 0 for even m + n, the number of triadic resonances in the weak shear limit increases to (m + n)/2 2 for a given m + n (note that the expression N H = (m + n)/2 (m + n − 1)/2 in figure 1 In summary, the addition of a weak shear substantially increases the number of triadic resonance interactions in finite-depth uniform stratifications. In the following section, we evaluate the solutions derived in the weak shear limit for an idealized background shear flow in the ocean and subsequently validate the same with numerical solutions. Finite shear regimes (finite Ri) are then explored numerically, with a focus on the dependence of the number of triadic resonances on the Richardson number Ri.

Results
The ocean surface boundary layer is characterized by intense mixing and homogenized properties. Formed by processes such as wind stress, oceanic circulation and wave breaking, it can extend to large depths by convection, and interaction between surface waves and ocean currents (D'Asaro 2014). The generation of surface waves was analysed initially by Miles (1957), who considered a shear flow in the air layer, overlying a quiescent water layer. Subsequent studies by Valenzuela (1976), Kawai (1979) and van Gastel, Janssen & Komen (1985) have incorporated a shear flow in the water layer also, assuming it to be setup by wind stress. This wind drift velocity profile is observed to be logarithmic in field measurements (Bye 1965;Churchill & Csanady 1983) and is characterized by the wind drift layer depth, h w (typically much smaller than the depth of the ocean) and the surface velocity, U s (∼3 %-4 % of wind speed). In the current study, we investigate the effects of the wind drift velocity profile on superharmonic generation by internal wave triadic resonances. Specifically, we consider an exponential velocity profile (Zeisel, Stiassnie & Agnon 2008;Young & Wolfe 2014); it is both amenable to analytical calculations and provides results qualitatively similar to more realistic velocity profiles (Morland, Saffman & Yuen 1991;Young & Wolfe 2014). The exponential background velocity profile in the ocean is assumed to be, where δ = U s /(NH) and ξ = L s /H are as defined in § 2.2. Here, U(z) is maximum at the free surface (z = 1) and negligible at the ocean bottom (z = 0). Unlike in § 2.2, δ is not necessarily assumed small in this section. As a result, the Richardson number Ri = N 2 L 2 s /U 2 s is allowed to assume arbitrary values. As shown in Appendix A, the dispersion curves and vertical mode shapes (governed by (2.5)) in the presence of the background flow in (3.1) can be analytically obtained. The dispersion curves in the presence of background flow are not symmetric about ω = 0 (figure 2), and the cograde (modes that travel faster than the maximum background flow velocity) and retrograde (modes that travel slower than the minimum background flow velocity) modes have to be treated separately. For the cograde modes, the phase velocity rapidly approaches the asymptotic value of U(1) as the mode number is increased, while the approach to U(0) for retrograde modes is relatively slower (figure 2b). With respect to mode shapes, weak shear has a weak influence (figure 3a), while finite shear tends to accumulate zero crossings close to the boundaries (upper boundary for cograde modes and lower boundary for retrograde modes, see figure 3b,c). Owing to this accumulation, different mode shapes tend to look similar over most of the domain except near the boundaries (cograde modes 2 and 3 in figure 3(b), for example), with similar phase speeds ( figure 2b). Consequently, the difference in the vertical structures of various high modes is not significant in most part of the domain, and resonant interaction between such modes is crucial to consider (Tung, Ko & Chang 1981).
In addition to the discrete modes shown in figure 3, there exists a singular continuous spectrum of modes whose phase speed matches with the background flow velocity at some z, namely the critical layer (Banks et al. 1976;Jose et al. 2015). In this study, we will not consider such continuous spectrum solutions, for either the primary modes at frequency ω or the superharmonic modes at frequency 2ω. Hence, owing to the non-interaction theorem (Eliassen & Palm 1961), no energy or momentum exchange between internal waves and the background flow can occur, up to at least O( 2 ) (Tung et al. 1981).

Weak shear limit
We start by evaluating the asymptotic weak shear limit solutions of § 2.2 for representative modal combinations, verify if the predicted new resonances occur in the presence of shear and finally compare the asymptotic theory with numerical solutions. The conditions for divergence of the O( 2 ) superharmonic wave field give the triadic resonance criteria for the interaction between modes m, n at frequency ω and the superharmonic wave at frequency 2ω. The weak shear asymptotic theory in § 2.2 suggests that modes m and n at frequency ω and mode q at frequency 2ω are in triadic resonance if and only if k m + k n = k q , where q = |m − n| if there is no shear, and q being any integer less than or equal to (m + n − 1)/2 in the presence of a weak shear with υ (z) / = 0 somewhere in the domain. While such triadic resonances can occur only for m/3 < n < 3m with no shear (Thorpe 1966), no such restrictions exist in the presence of a shear flow. In other words, all the intersection points in figure 1(a) for a uniform stratification become triadic resonances in the presence of a weak shear.  Figure 4 shows the amplitude of the superharmonic wave (maximum value ofh mn (z), which is governed by (2.12)) for representative modal interactions of (m, n) = (2, 3) and (2, 5), plotted as a function of the primary wave frequency ω. For (m, n) = (2, 3),h max mn becomes infinitely large only at ω ≈ 0.468 for Ri = ∞ (no shear), with the superharmonic internal wave being of mode number q = |m − n| = 1 (blue curve in figure 4a). In the presence of weak shear (Ri = 10 7 ), an additional divergence ofh max mn appears at ω ≈ 0.327, which is the location where k q = k m + k n is satisfied with q = 2. Both the weak shear asymptotic theory ( § 2.2) and fully numerical solution of (2.12) confirm that triadic resonance occurs at both ω ≈ 0.468 and 0.327. For (m, n) = (2, 5), while only one resonance location exists (ω ≈ 0.285) with no shear (blue curve in figure 4b), two additional divergence locations appear with weak shear (grey dashed curve in figure 4b).
The weak shear asymptotic theory is again shown to faithfully recover the new divergences (and hence triadic resonances) in the presence of weak shear for (m, n) being equal to (2, 5) (black dotted curve in figure 4b). In summary, figure 4 confirms the main conclusion from the weak shear asymptotic theory for two representative modal interactions: in the presence of arbitrarily weak shear, additional triadic resonance locations emerge at all the locations where the horizontal wavenumber condition is satisfied. In addition, we also verified that divergence ofh max mn due to triadic resonances resulting from self-interaction is also possible in the presence of weak shear.

Finite shear
At finite shear, while it is possible to take a semi-analytical approach to solve (2.5) and (2.12) for the exponential background velocity profile (Appendix A), we present fully numerical solutions of (2.12) in this subsection owing to the simplicity in obtaining them. Using the shooting method alongside the fourth-order Runge-Kutta scheme to march from z = 0 to z = 1, (2.5) is numerically solved to obtain the horizontal wavenumbers and the vertical mode shapes of different modes at a given primary wave frequency ω. The boundary-value problem in (2.12) is then solved using a second-order finite difference scheme to obtain the superharmonic vertical structureh mn (z) for different (m, n) interactions. In the parameter space of (ω, Ri) ∈ [0.01, 0.99] × [0.50, 10 7 ], we calculate the amplitude of the superharmonic wave (h max mn ) and identify divergences via peaks inh max mn which become stronger with finer resolution in ω. The superharmonic wave mode number (q) is calculated throughout the parameter space using the number of zero crossings in the vertical structure ofh mn (z). It is worth pointing out that for ω > 0.50, while the superharmonic frequency is evanescent for Ri = ∞, propagating superharmonic internal wave modes exist for finite Ri (Bell 1974).
The distributions of log 10 [h max mn ] on the (ω, Ri) plane for three representative cograde modal interactions (m, n) = (1, 2), (2, 3) and (2, 2) are shown in figure 5(a-c). The divergence locations in the log 10 [h max mn ] vs ω plot (like in figure 4) from different Ri form the 'divergence curves' in figure 5. All the points along these divergence curves represent triadic resonance locations. In the limit of very large Ri, the divergence curves occur at locations predicted by the weak shear asymptotic theory (see hollow circles in figure 5a-c). For (m, n) = (1, 2), the triadic resonance between modes 1 and 2 at frequency ω and mode-1 at frequency 2ω occurs at all values of Ri, with ω being at 0.395 at large Ri and increasing towards 0.482 at Ri = 0.50 (figure 5a). For (m, n) = (2, 3), two different divergence locations are predicted in the weak shear limit, and they are observed to extend as divergence curves over a wide range of Ri. Similar to what was observed for (m, n) = (1, 2), the divergence curve corresponding to a mode-2 superharmonic internal wave deviates slightly from its weak shear limit location when Ri reaches small values. In contrast, the divergence curve corresponding to the mode-1 superharmonic internal wave departs significantly from its ω value in the weak shear limit as Ri becomes small. Interestingly, it does not even seem to extend all the way to Ri = 0.50. A similar behaviour is observed in the self-interaction case presented in figure 5(c), where the divergence curve associated with a mode-1 superharmonic internal wave goes from ω ≈ 0.447 at weak shear and towards large ω at Ri ≈ 1.50 (ω ≈ 0.719). It is worth recalling from § 2.2 that resonance due to self-interaction is not possible at all if there is no shear. Figure 5(a-c) shows that the primary wave frequency at the triadic resonance locations deviates very little from its weak shear limit value if Ri is larger than approximately 10 2 . For Ri < 10 2 , the divergence curves explore a larger range of primary wave frequencies, even extending to values for which the superharmonic frequency would be evanescent if there was no shear. To closely investigate the small Ri region (0.50 Ri 5), figure 5(d-f ) shows zoomed-in views of figure 5(a-c), respectively. For (m, n) = (1, 2), triadic resonance of the primary waves with the mode-1 superharmonic wave occurs at around ω ≈ 0.482 for small Ri close to 0.50, which is not far from the resonant value of ω ≈ 0.395 for Ri → ∞. Similarly, the triadic resonance associated with the mode-2 superharmonic wave for (m, n) = (2, 3) occurs at similar ω for Ri = 0.50 and Ri → ∞. In contrast, the triadic resonance associated with the mode-1 superharmonic wave ceases to exist below a cutoff Ri of 2.94 (figure 5e). Interestingly, for a given Ri larger than (but in the vicinity of) 2.94, resonant generation of a mode-1 superharmonic internal wave occurs at two different values of ω. For example, triadic resonance with the mode-1 superharmonic internal wave occurs at ω ≈ 0.743 and 0.816 for Ri ≈ 2.98, and at ω ≈ 0.587 and 0.982 for Ri ≈ 4.89. Thus, for a finite range of Ri ∈ (2.94, 5.01), there exists two different ω values at which modes 2 and 3 at ω and mode-1 at 2ω are in triadic resonance. A similar behaviour is observed with the triadic resonance between a self-interacting mode-2 at ω and mode-1 at 2ω (figure 5f ), The corresponding finite range of Ri where two resonant ω values exist is Ri ∈ (1.49, 2.38). Finally, it is worth highlighting that weakly nonlinear effects seem to be stronger at smaller Ri in general (including the regions away from divergence curves), as evidenced by the larger magnitudes of the superharmonic wave amplitude at small Ri in figure 5(e,f ).
As pointed out in § 2.2, a triadic resonance in the vicinity of ω = 0 appears for those (m, n) for which m + n is even, if a weak shear is present. This weak shear limit is indicated by the hollow circle at ω = 0 and large Ri in figure 5(c), where m + n = 4 is even. An increase in the superharmonic amplitude as one approaches the ω = 0 axis is evident at all Ri, although we could not establish that a divergence curve exists in the neighbourhood for ω > 0. Upon further investigation, we found that the actual divergence curve occurs in the ω < 0 region for all Ri, which still represents a continuous extension of what occurs at ω = 0 in the weak shear limit. An alternate view of this observation is that the (2, 2) triadic resonance interaction at ω = 0 for Ri → ∞ extends to the finite Ri region as a (−2, −2) triadic resonance interaction at negative ω values, i.e. a retrograde-retrograde interaction. This aspect is further elucidated in figure 8.
In figure 6, we provide an interpretation of our observations for the (m, n) = (2, 3) interaction (figure 5b,e) in the context of internal wave dispersion relation in the presence of shear. Like in figure 1(a), we plot the dispersion curves for the individual modes (shown in grey) and the 2ω vs (k m + k n ) curves obtained from the dispersion curves for modes m and n at frequency ω (shown in black) for various Ri in figure 6. At large Ri (Ri = 100 in figure 6a), k m + k n = k q (q is the mode number of the superharmonic wave) is satisfied at two different locations, each corresponding to q = 1 (blue dot) and 2 (red dot). The corresponding primary wave frequencies are where divergence locations occur at Ri = 100 in figure 5(b), and are very close to those values predicted by the weak shear asymptotic theory. At Ri = 5, triadic resonance of the primary waves with the mode-1 superharmonic wave occurs at two different values of 2ω (blue dots in figure 6b), again being consistent with the divergence locations at Ri = 5 in figure 5(b). As Ri is further decreased to the previously identified cutoff value of 2.94 (figure 6c), the two blue dots from figure 6(b) converge to a single location, and the mode-1 dispersion curve is now tangential to the 2ω vs (k m + k n ) curve. At an even smaller value of Ri (=0.50 in figure 6d), triadic resonance with the mode-1 superharmonic wave is absent (no blue dots), in line with the observation from figure 5(e) that it ceases to exist below Ri ≈ 2.94. Triadic resonance with the mode-2 superharmonic wave persists at all Ri, indicated by the red dots in each of figure 6(a-d).
In summary, figure 6 establishes that divergence of the superharmonic wave amplitude (as shown in figure 5) occurs at the same locations where the horizontal wave numbers satisfy k m + k n = k q , and hence represent triadic resonances. In addition, figure 6 also establishes that horizontal wavenumber condition is sufficient to ensure triadic resonance between the primary and superharmonic waves for finite shear too ( § 3.1 showed the same result in the weak shear limit). Extending the investigations in figures 5 and 6, we find that features such as multiple resonant ω values at given Ri and superharmonic mode number, and disappearance of some resonances below a cutoff Ri are common for higher mode interactions too. These features can have significant implications for the total number of possible triadic resonances at various Ri, which we proceed to discuss in § 3.3.

Number of resonance locations
In oceanic settings like internal tide generation or scattering by ocean floor topography, a finite number of modes at the tidal frequency are excited, with the low modes often containing the most energy (Garrett & Kunze 2007). In this subsection, we consider a scenario where a finite amount of energy is present in the first few modes at a fixed frequency, and investigate the total number of resonant interactions possible amongst them. Specifically, we consider all possible cograde-cograde interactions with (m, n) (5, 5), and plot all the corresponding divergence curves (obtained as locations where k m + k n = k q is satisfied) on the ω-Ri plane (figure 7a). All the divergence curves emerge from the weak shear limit (blue dots at large Ri), and expand towards other frequencies as Ri is decreased. While some of them extend all the way to Ri = 0.5, others have a cutoff Ri below which the corresponding resonances do not exist. As pointed out in § 3.2, some of the divergence curves occur at two different primary wave frequencies at a given Ri within a particular range. In this range of Ri (roughly around 1 to 100), the divergence curves extend over a large frequency range, which is in contrast with what happens at very large and very small Ri. Accounting for primary modes up to 10 (figure 7b), it becomes evident that almost every primary wave frequency corresponds to some resonant interaction for a finite band of Ri. Owing to existence of a cutoff Ri for several of the interactions, only a few interactions are possible as Ri approaches 0.5. For a given Ri, the total number of resonant interactions (N R ) over the entire range of ω ∈ [0, 1] is calculated from the plots in figure 7(a,b), and N R is subsequently plotted as a function of Ri in the bottom row of figure 7. As already pointed out, the number of resonances in the weak shear limit (Ri → ∞) is much larger than the no shear limit. In figure 7(c), which corresponds to interactions amongst the first 5 modes, the weak shear and no shear limits are shown by the red and blue dots, respectively. Considering primary wave frequencies in the range [0, 1], N R remains at 33 for 63 Ri < ∞, and starts to increase with a further decrease in Ri. It attains a maximum of 35 in the interval of 16.26 Ri 26.38, before decreasing towards 10 at Ri ≈ 1.19. N R then remains at 10 for 0.5 Ri 1.19. Considering a larger range of ω (ω ∈ [0, 1.5], shown in black), N R increases  −1, 4), plotted on the plane of primary wave frequency ω on the x-axis and Richardson number Ri on the y-axis. Hollow circles at Ri = 10 7 indicate the locations of divergence calculated from weak shear asymptotic theory. The mode number of the superharmonic internal wave at the divergence locations is indicated by the encircled numbers next to the corresponding divergence curves. The negative sign before mode number indicates that it is a retrograde mode. from its weak shear limit at a larger Ri of 268.79, and also attains a larger peak of N R = 45 at Ri ≈ 32.01. Increasing the maximum mode number considered to 10, a similar trend in N R vs Ri is observed, but with a larger value of Ri at which N R deviates from its weak shear limit (figure 7d).
In summary of figure 7, we have identified a finite range of Ri over which resonant interactions occur at a significantly wider range of primary wave frequencies than for large and small Ri. In addition, the total number of possible resonant interactions attains a maximum for intermediate Ri of O(10), before decreasing towards relatively small values at small Ri. N R at small Ri, while being much smaller than the weak shear limit, is still larger than the number of resonances in the no shear limit.
Due to symmetry breaking by the introduction of shear, retrograde-retrograde and mixed (cograde-retrograde) interactions behave differently from cograde-cograde interactions. To elucidate this, we show the contour plots of log 10 [h max mn ] for representative modal interactions of (m, n) = (−2, −4), (3, −5) and (−1, 4) in figures 8(a)-8(c), respectively. Here, a negative mode number −m indicates a retrograde mode with mode number m. Unlike the dispersion curves of cograde modes (see figure 6, for example), the dispersion curves of retrograde modes are bounded by a maximum frequency, whose magnitude increases with Ri and approaches ω = 1 in the limit of no shear (Bell 1974). This restricts the maximum primary wave frequency at which divergence curves can occur to ω = 0.5, thus restricting the primary wave frequency axis to ω ∈ [0, 0.5] in figure 8.
In the retrograde-retrograde interaction of (m, n) = (−2, −4) (figure 8a), three divergence curves emerge at Ri → ∞ from the resonance locations predicted by the weak shear asymptotic theory (hollow circles), including one from ω = 0. It is worth highlighting that the weak shear resonances occur at locations where the horizontal resonance condition is satisfied with no shear ( § 2.2), thus giving rise to the same resonance-emerging locations for cograde-cograde and retrograde-retrograde interactions. All the three divergence curves extend all the way to Ri = 0.5, with noticeable deviations from the weak shear limit. Interestingly, further investigations revealed that all the retrograde-retrograde resonance interactions that occur in the weak shear limit extend all the way to Ri = 0.5. This is in contrast to cograde-cograde interactions, some of which were shown to have a cutoff Ri below which they do not exist. As was pointed out for even m + n in cograde-cograde interactions (figure 5c), while a signature of the weak shear resonance at ω = 0 was present in the form of enhanced superharmonic amplitudes, no divergence curve was actually detected (figure 5c). It turns out that the ω = 0 resonances from the weak shear limit become retrograde-retrograde resonances for finite Ri, which explains the occurrence of a divergence curve with a retrograde mode 3 superharmonic wave in figure 8(a).
For the cograde-retrograde mixed interaction of (m, n) = (3, −5) (figure 8b), a divergence curve (with the superharmonic wave being a retrograde mode 1) appears in the vicinity of ω = 0 at large Ri, emerging from the weak shear limit of ω = 0 that is possible due to even m + n. While the resonance at Ri = 10 7 is not very close to ω = 0, we confirmed with further investigations for Ri > 10 7 that the resonance location approaches ω = 0 as Ri is made even larger. This divergence curve extends all the way to Ri = 0.5, with its location being at ω ≈ 0.416 at Ri = 0.5. In addition, a new divergence curve (with the superharmonic wave being a retrograde mode 2) emerges from (ω, Ri) = (0, 2.42); the corresponding resonance has no weak shear counterpart. Such an occurrence of a new resonance at finite Ri and ω = 0 is observed only for mixed interactions. This new divergence curve is present at all Ri 2.42, occurring at ω ≈ 0.304 for Ri = 0.5. Finally, for the mixed interaction case of (m, n) = (−1, 4) (figure 8c), a single divergence curve is observed, and it emerges from the weak shear limit of ω ≈ 0.395. Unlike for the examples in figure 8(a,b), the resonance in figure 8(c) has a cutoff Ri of 5.01 below which it does not occur.
Having investigated representative retrograde-retrograde and mixed interactions in figure 8, and recognizing that they are quite different from cograde-cograde interactions, we proceed to reproduce the calculations in figures 7(a) and 7(c) in the presence of all interactions. Specifically, considering all possible interactions between cograde and retrograde mode numbers up to 5, we plot all the divergence curves on the ω-Ri plane, with the colour indicating the type of interaction (figure 9a). The divergence curves corresponding to cograde-cograde interactions are reproduced from figure 7(a). All the divergence curves associated with retrograde-retrograde interactions extend all the way from their weak shear limit to Ri = 0.5, with quite a few of them emerging from ω = 0 at Ri → ∞. Interestingly, the ω = 0 resonances expand towards larger values of ω as Ri is decreased, thus making a significant range of ω ∈ [0, 0.5] susceptible to retrograde-retrograde resonances at small Ri. Mixed interactions are fewer in number compared with other interactions at large Ri, but new resonances emerge at finite Ri, a feature that is absent for cograde-cograde and retrograde-retrograde interactions.
As in figure 7(c), we plot the total number of resonances (N R ) as a function of Ri, including the contributions from all types of interactions (figure 9b). While the cograde-cograde contribution is reproduced (and already discussed) from figure 7(c), the retrograde-retrograde contribution is observed to be significant and invariant with Ri (blue markers in figure 9b). The latter feature can be understood from figure 9(a), where all the retrograde-retrograde resonances extended all the way from Ri → ∞ to Ri = 0.5. The number of mixed interactions is relatively small at large Ri, but increases noticeably once Ri is decreased below a value of around 23.95. At small Ri of O(1), the retrograde-retrograde contribution dominates, while the number of mixed interactions has overtaken the number of cograde-cograde interactions. In terms of the total number of resonances, N R remains more or less constant at its weak shear limit value of 85 for Ri 63, which is much larger than the no-shear value of 14. Upon a decrease in Ri below 63, N R increases, as already observed for cograde-cograde interactions in 7(c). At small Ri, owing to the new resonances from mixed interactions, N R does not drop as significantly as the number of cograde-cograde interactions, and has a value of 73 at Ri = 0.5. Finally, increasing the considered range of primary wave frequencies to ω 1.5 noticeably increases the maximum value of N R due to its influence on the cograde-cograde interactions at finite Ri, but has little effect on N R at small Ri.
In summary, the total number of resonances with shear is significantly larger than that without shear, and attains a maximum at a Ri ≈ 17.91, Its value at small Ri is comparable to its weak shear limit owing to the emergence of new mixed interaction resonances at finite Ri. Furthermore, it is noteworthy that the range of primary wave frequencies that are susceptible to resonances becomes very large at finite Ri due to cograde-cograde interactions, and also occupies a significant range of ω ∈ [0, 0.5] at small Ri mainly due to retrograde-retrograde interactions.

Discussion and conclusions
In this study, triadic resonance resulting from interaction between discrete internal wave modes at the same frequency in a two-dimensional, uniformly stratified shear flow were considered. For a linear (primary) wave field comprising a series of modes at the same frequency ω, superharmonic (frequency 2ω) and mean-flow (frequency zero) terms constitute the weakly nonlinear solutions. At certain primary wave frequencies, the amplitude of the superharmonic term diverges as a result of triadic resonance between two primary modes and a superharmonic secondary mode. In the no shear limit, primary modes m and n at frequency ω are in triadic resonance with mode q at frequency 2ω if k m + k n = k q and q = |m − n| are satisfied, where k i denotes the horizontal wavenumber of mode i at the corresponding frequency (Thorpe 1966). Here, we developed an asymptotic theory to investigate the influence of weak shear on the aforementioned triadic resonances. We find that, unlike in the no shear limit, the horizontal wavenumber condition of k m + k n = k q alone is sufficient to ensure a resonant interaction, independent of the superharmonic wave mode number q. As a result, several more resonances, which include self interactions, occur in the presence of an arbitrarily weak shear when compared with the no shear limit. The locations of these resonances can be traced back to those primary wave frequencies at which k m + k n = k q is satisfied in the no shear limit, including some at the seemingly trivial case of ω = 0.
In § 3, we extended our investigations to finite shear by numerically solving the equations governing the weakly nonlinear superharmonic wave field. Specifically, an ocean-like exponential velocity profile was considered, and a systematic study on the effects of Richardson number Ri was performed. Here, Ri → ∞ corresponds to the weak shear limit, and our study spans a wide range of Ri all the way to Ri = 0.5. On the (ω, Ri) plane, the superharmonic wave amplitude diverges at locations of resonance, and the locus of all such points corresponding to a particular superharmonic internal wave of mode number q, is termed as a divergence curve of mode number q. For cograde-cograde interactions, a cutoff Ri exists for a large number of divergence curves and only a few among all the divergence curves extend from Ri → ∞ to Ri = 0.5. In addition, for a finite range of Ri (roughly around 1-100) above the cutoff Ri for a particular modal interaction, resonance occurs at two different primary wave frequencies at a given Ri. For retrograde-retrograde interactions, all the resonances that occur for Ri → ∞ extend all the way to Ri = 0.5 without the occurrence of any cutoff Ri. Further, all the resonances that occur at zero primary wave frequency for Ri → ∞ appear in retrograde-retrograde interactions at finite ω for finite Ri. For mixed interactions, a new feature of the emergence of new resonances at finite Ri ∼ O(10) is observed, thus contributing to an increase in the number of mixed interaction resonances at small Ri. The total (from cograde, retrograde and mixed interactions) number of resonance locations has a dramatic increase from the no shear limit to the weak shear limit (Ri → ∞), attains a maximum at moderate Ri (∼O (10)) and approaches a value that is not far from the weak shear limit value, at Ri = 0.5. This trend, particularly at small Ri, is understood as a consequence of a decrease and increase of the cograde-cograde and mixed resonant interactions, respectively, with Ri.
Our conclusions based on the weak shear limit are potentially relevant for other inhomogeneous background conditions as well. Specifically, we showed that all the locations where the horizontal wavenumbers satisfy the triadic resonance condition in the absence of shear represent actual triadic resonances in the presence of an arbitrarily weak shear. A similar conclusion will hold if, instead of a weak shear, a weak non-uniformity in the stratification is introduced. This also explains why the previous study by Varma & Mathur (2017) identified several more resonances in non-uniform stratifications compared with a uniform stratification. While Varma & Mathur (2017) attributed the sufficiency of the horizontal resonance condition k m + k n = k q to a non-orthogonality condition being satisfied by modal pairs in a non-uniform stratification, an equivalent perspective in our weak shear limit theory is that the coefficient of 'sin(γ z)/ sin γ ' in (2.29) is non-zero. In summary, a weak inhomogeneity either in terms of shear or non-uniformity in stratification reduces the dimension of the spatial triadic resonance condition, and hence substantially increases the number of resonances compared with the case of uniform stratification with no shear.
A recent study by Biswas & Shukla (2021) directly used the sufficiency of frequency and horizontal wavenumber conditions to identify internal wave resonant triads in a uniformly stratified uniform shear flow. Considering primary modes (m, n) at the same frequency ω, Biswas & Shukla (2021) investigated the stability of a few resonant triads (specifically five different (m, n) combinations at a few arbitrarily chosen values of moderate Ri) that contained the superharmonic wave. In contrast to Biswas & Shukla (2021), we consider a realistic ocean-like background shear flow, and identify all possible low-mode interactions over a continuous range of Ri spanning from weak to moderate to strong shear. An analysis of the amplitude evolution associated with all the resonant triads we identified in this study would be useful in determining their relative importance in a realistic internal wave field comprised of different modes.
Our study indicates that superharmonic generation due to triadic resonance is likely for a large range of primary wave frequencies, particularly at Ri ∼ O (1-100), if a few different primary modes are simultaneously excited. In the future, it is important to consider the parameter values (Richardson number, frequency and mode numbers) for representative oceanic regimes (internal tides in the presence of wind forced shear flow or surface gravity wave-induced mean flow, wave-mean-flow interaction etc.) to study the extent to which superharmonic generation gets modified by background shear. It would further be interesting to study if the new resonances that occur due to shear can result in stronger secondary wave growth than those that are associated with the resonances in uniform stratification with no shear. This would require the derivation of amplitude evolution equations associated with the triadic resonances identified in this study. Additionally, incorporating the presence of background shear along with one or more of non-uniform stratification, background rotation, viscous and three-dimensional effects will take us a step closer to realistic oceanic settings. In a three-dimensional domain, additional classes of resonances (wave-vortex and vortex-vortex) analysed by Lelong & Riley (1991) would also be relevant. The ideas in this paper may be relevant for triadic resonance with subharmonic waves too, including PSI i.e. apart from the effects of background shear on the already known subharmonic resonances (such as PSI), it is important to study if additional resonances occur due to background shear. Finally, it would also be interesting to study the effects of background shear on triadic resonances when critical layers, and hence the possibility of a continuous spectrum of modes, are present.
Here, C 1 and C 2 are integration constants. This solution above is valid in the neighbourhood of the regular singular point μ = 0, with radius of convergence |μ| < 1.