1. Introduction
As we celebrate two hundred years of the Navier–Stokes equations, it is remarkable to note that hydrodynamic turbulence is still a much studied subject. Among the different fields of study, there is wave turbulence that has become important in geophysics and astrophysics where waves are omnipresent (Galtier Reference Galtier2023b). The strength of the (weak) wave turbulence theory is that it offers the possibility of a deep understanding of physical systems composed of a set of random waves interacting nonlinearly. The reason for this is, first of all, the possibility of analytically deriving a set of integro-differential equations for second-order spectral cumulants – the so-called kinetic equations – which are free from the closure problem classically encountered in eddy turbulence. Indeed, in wave turbulence there is a natural asymptotic closure rooted on the existence of a small parameter, the wave amplitude. Secondly, exact solutions (Kolmogorov–Zakharov spectra) can be found from the kinetic equations. In addition to the usual thermodynamic solutions, the kinetic equations have finite flux solutions that capture the flow of conserved densities from sources to sinks. Thirdly, these solutions correspond to power law spectra that can be compared with the data. The number of experiments, observations and diagnostics have increased considerably over the past two decades and today, thanks also to direct numerical simulations (DNS), wave turbulence has become a leading field in turbulence where new fundamental questions are being raised (Galtier & Nazarenko Reference Galtier and Nazarenko2017; Hassaini & Mordant Reference Hassaini and Mordant2017; Hassaini et al. Reference Hassaini, Mordant, Miquel, Krstulovic and Düring2019; Savaro et al. Reference Savaro, Campagne, Linares, Augier, Sommeria, Valran, Viboud and Mordant2020; Galtier & Nazarenko Reference Galtier and Nazarenko2021; Ricard & Falcon Reference Ricard and Falcon2021; David & Galtier Reference David and Galtier2022; Falcon & Mordant Reference Falcon and Mordant2022; Griffin et al. Reference Griffin, Krstulovic, L'vov and Nazarenko2022; Hrabski & Pan Reference Hrabski and Pan2022; Kochurin & Kuznetsov Reference Kochurin and Kuznetsov2022; Onorato et al. Reference Onorato, Dematteis, Proment, Pezzi, Ballarin and Rondoni2022; Rodda et al. Reference Rodda, Savaro, Davis, Reneuve, Augier, Sommeria, Valran, Viboud and Mordant2022; Zhang & Pan Reference Zhang and Pan2022; Dematteis & Lvov Reference Dematteis and Lvov2023; Galtier Reference Galtier2023a; Lanchon et al. Reference Lanchon, Mora, Monsalve and Cortet2023; Novkoski, Pham & Falcon Reference Novkoski, Pham and Falcon2023; Zhu et al. Reference Zhu, Semisalov, Krstulovic and Nazarenko2023).
Rotating fluids are one of the most studied examples in (strong/weak) wave turbulence as it involves inertial waves that are easy to excite experimentally (Hopfinger, Gagne & Browand Reference Hopfinger, Gagne and Browand1982; Jacquin et al. Reference Jacquin, Leuchter, Cambon and Mathieu1990; Morize, Moisy & Rabaud Reference Morize, Moisy and Rabaud2005). This regime is of interest in a number of fields such as geophysics where the Coriolis force is felt, for example, through large-scale atmospheric motions. Very early, it was recognized that rotating turbulence behaves differently from classical eddy turbulence with a reduction of the cascade along the axis of rotation 
${\boldsymbol \varOmega _0}$ with possibly a steeper energy spectrum than the well-known Kolmogorov spectrum (Hossain Reference Hossain1994; Zeman Reference Zeman1994; Zhou Reference Zhou1995; Cambon, Mansour & Godeferd Reference Cambon, Mansour and Godeferd1997; Smith & Waleffe Reference Smith and Waleffe1999; Baroud et al. Reference Baroud, Plapp, She and Swinney2002; Godeferd & Moisy Reference Godeferd and Moisy2015). Another remarkable feature, still not understood, is the self-similarity found in the scaling of velocity structure functions in the direction transverse to 
${\boldsymbol \varOmega }_0 = \varOmega _0 \boldsymbol {e}_\parallel$ (Baroud et al. Reference Baroud, Plapp, She and Swinney2002; van Bokhoven et al. Reference van Bokhoven, Clercx, van Heijst and Trieling2009), which is in strong contrast with intermittency observed in hydrodynamic turbulence. In recent years, the (weak) inertial wave turbulence regime has been specifically studied experimentally (Yarom & Sharon Reference Yarom and Sharon2014; Monsalve et al. Reference Monsalve, Brunet, Gallet and Cortet2020). For example, it was shown that the energy spectrum is concentrated along the dispersion relation, as expected in wave turbulence, with a scaling in agreement with the theoretical prediction. Numerical simulations (including DNS) have also been carried out, notably to study the spectral properties (Bellet et al. Reference Bellet, Godeferd, Scott and Cambon2006; Scott Reference Scott2014; Clark di Leoni & Mininni Reference Clark di Leoni and Mininni2016; Le Reun et al. Reference Le Reun, Favier, Barker and Le Bars2017; Sharma, Verma & Chakraborty Reference Sharma, Verma and Chakraborty2018; Galtier & David Reference Galtier and David2020; Le Reun, Favier & Le Bars Reference Le Reun, Favier and Le Bars2020; Yokoyama & Takaoka Reference Yokoyama and Takaoka2021). They confirm the previous experimental and theoretical results and reveal, for example, the existence of a non-stationary solution different from the stationary Kolmogorov–Zakharov spectrum.
The theory of inertial wave turbulence has been derived by Galtier (Reference Galtier2003). The kinetic equation obtained is valid in the most general case, i.e. without making the assumption of anisotropy. However, a simple argument based on the resonance condition shows that the cascade is anisotropic with a transfer mainly in the perpendicular (
$\perp$) direction to 
${\boldsymbol{\varOmega} }_0$. Using this feature, the kinetic equation was eventually reduced to the axisymmetric case from which the exact (Kolmogorov–Zakharov) energy spectrum was derived. This solution takes the form 
$E(k_{\perp },k_{\parallel }) \sim k_{\perp }^{-5/2} k_{\parallel }^{-1/2}$, with 
$\boldsymbol {k} = \boldsymbol {k}_{\perp } + k_{\parallel } \boldsymbol {e}_\parallel$ the wavevector. As recently proved by David & Galtier (Reference David and Galtier2023), this energy spectrum corresponds to a local turbulence (with an inertial range independent of the largest and smallest scale) for which we can also estimate the Kolmogorov constant. The derivation of the kinetic equation of inertial wave turbulence in the general case (without the axisymmetry assumption) is cumbersome and the use of the Hamiltonian formalism does not drastically simplify the calculation (Gelash, L'vov & Zakharov Reference Gelash, L'vov and Zakharov2017). Here, it is proposed to take another path for the derivation of such an equation by using the anisotropy assumption (
$k_{\perp } \gg k_{\parallel }$) at the beginning of the analysis. As shown in § 2, in this case the wave amplitude equation can be obtained without the introduction of a complex helicity basis, which is an interesting simplification: the velocity field is decomposed into poloidal and toroidal fields from which we can define the canonical variables. In § 3 the multiple time scale method introduced by Benney & Saffman (Reference Benney and Saffman1966) is adapted (a few points are also clarified) to this anisotropic problem and then used to derive such a kinetic equation that is the same as the original work when anisotropy is finally taken into account. This result shows the commutativity between asymptotic closure and anisotropy. Furthermore, the multiple time scale method informs us that the kinetic equation (of weak wave turbulence) can be derived without imposing restrictions on the probability distribution of the wave amplitude such as quasi-Gaussianity, or on the phase such as random phase approximation that naturally occurs dynamically. Finally, we conclude with a discussion in § 4.
2. Wave amplitude equation
2.1. Canonical variables
The Navier–Stokes equations with the Coriolis force read
\begin{equation} \frac{\partial {\boldsymbol w}}{\partial t} - 2 ({\boldsymbol{\varOmega}}_0 \boldsymbol{\cdot} \boldsymbol{\nabla}) {\boldsymbol u} = ({\boldsymbol w} \boldsymbol{\cdot} \boldsymbol{\nabla}) {\boldsymbol u} - ({\boldsymbol u} \boldsymbol{\cdot} \boldsymbol{\nabla}) {\boldsymbol w} + \nu \nabla^2 {\boldsymbol w}, \end{equation}
where 
${\boldsymbol u}$ is a solenoidal velocity (
$\boldsymbol {\nabla } \boldsymbol {\cdot } {\boldsymbol u} = 0$), 
${\boldsymbol w} = \boldsymbol {\nabla } \times {\boldsymbol u}$ the vorticity and 
${\boldsymbol{\varOmega} }_0$ a constant rotation rate. Hereafter, we neglect the viscosity. We introduce the toroidal (
$\psi$) and poloidal (
$\phi$) scalar fields in the following manner:
\begin{equation} {\boldsymbol u} = \boldsymbol{\nabla} \times (\psi \boldsymbol{e}_\parallel) + \boldsymbol{\nabla} \times (\boldsymbol{\nabla} \times (\phi \boldsymbol{e}_\parallel)). \end{equation}Here the Fourier transform writes
\begin{equation} \hat {\boldsymbol u}_k = {\rm i} \hat \psi_k \boldsymbol{k} \times \boldsymbol{e}_\parallel - \hat \phi_k \boldsymbol{k} \times (\boldsymbol{k} \times \boldsymbol{e}_\parallel ) = {\rm i} \hat \psi_k \boldsymbol{k} \times \boldsymbol{e}_\parallel + \hat \phi_k ( k^2 \boldsymbol{e}_\parallel - k_{{\parallel}} \boldsymbol{k} ), \end{equation}
from which we deduce the vorticity vector (
$\vert \boldsymbol {k} \vert =k$)
\begin{equation} \hat {\boldsymbol w}_k = \hat \psi_k ( k^2 \boldsymbol{e}_\parallel - k_{{\parallel}} \boldsymbol{k}) + {\rm i} k^2 \hat \phi_k \boldsymbol{k} \times \boldsymbol{e}_\parallel. \end{equation}It is straightforward to show in Fourier space that the linear contribution of (2.1) leads, after projection, to
\begin{gather} \frac{\partial \hat \phi_k}{\partial t} = 2{\rm i} \varOmega_0 \frac{k_{{\parallel}}}{k^2} \hat \psi_k , \end{gather}
\begin{gather}\frac{\partial \hat \psi_k}{\partial t} = 2{\rm i} \varOmega_0 k_{{\parallel}} \hat \phi_k . \end{gather}
The linear solutions are the well-known (helical) inertial waves with the angular frequency (
$\partial ^2_t = - \omega ^2_k$ can be used)
\begin{equation} \omega^2_k = 4 \varOmega^2_0 \frac{k_{{\parallel}}^2}{k^2} . \end{equation}From this property, we introduce the canonical variables
\begin{equation} A^{s}_k \equiv A^{s} (\boldsymbol{k}) = k^2 \hat \phi_k - s k \hat \psi_k , \end{equation}
with 
$s=\pm$ the directional polarity. With such a choice of canonical variables, we have
\begin{equation} \vert A^+_k \vert^2 + \vert A^-_k \vert^2 = 2 \vert \hat {\boldsymbol u}_k \vert^2 \end{equation}and at the linear level
\begin{equation} \frac{\partial A^s_k}{\partial t} + {\rm i} s \omega_k A^s_k = 0 . \end{equation}2.2. Resonance condition
The resonance condition for three-wave interactions can be written as (Galtier Reference Galtier2023b)
\begin{gather} s \omega_k + s_p \omega_p + s_q \omega_q = 0, \end{gather}
\begin{gather}\boldsymbol{k} + \boldsymbol{p} + \boldsymbol{q} = 0 . \end{gather}In the case of inertial waves, these relations are equivalent to the conditions
\begin{equation} \frac{s_q q - s_p p}{s \omega_k} = \frac{sk - s_q q}{s_p \omega_p} = \frac{s_p p - sk}{s_q \omega_q} . \end{equation}
Assuming that the system is initially excited at large scale in a narrow isotropic domain in Fourier space, a situation often considered in DNS, the dynamics will initially be dominated by local interactions such that 
$k \simeq p \simeq q$. As the locality of the interactions is a property of turbulence that is generally verified, we can extend its use beyond the initial instant. We obtain
\begin{equation} \frac{s_q -s_p}{s k_{{\parallel}}} \simeq \frac{s - s_q}{s_p p_\parallel} \simeq \frac{s_p - s}{s_q q_\parallel} . \end{equation}
From this expression, we can show that the associated cascade is necessarily anisotropic. Indeed, if 
$k_{\parallel }$ is non-zero, the left-hand term will only give a non-negligible contribution when 
$s_{p}=-s_{q}$. The immediate consequence is that either the middle or the right-hand term has its numerator that cancels (to leading order), which implies that the associated denominator must also cancel (to leading order) to satisfy the equality: for example, if 
$s=s_{p}$ then 
$q_\parallel \simeq 0$. This condition means that the transfer in the parallel direction is negligible: indeed, the integration of the wave amplitude equation in the parallel direction (see below) is then reduced to a few modes (since 
$p_\parallel \simeq k_{\parallel }$), which strongly limits the transfer between parallel modes. The cascade in the parallel direction is thus possible but relatively weak compared with that in the perpendicular direction. In the following, we take advantage of this property and consider the anisotropic limit 
$k_{\perp } \gg k_{\parallel }$ to simplify the derivation. Note that once turbulence is anisotropic, we can still use the locality condition with 
$k \sim k_\perp$; we then obtain 
$k_\perp \sim p_\perp \sim q_\perp$, whereas the parallel wavenumbers are limited to a narrow domain.
2.3. Wave amplitude equation
In the derivation of the wave amplitude equation, we consider a continuous medium that can lead to mathematical difficulties connected with infinite dimensional phase spaces. For this reason, it is preferable to assume a variable spatially periodic over a box of finite size 
$L$. However, in the derivation of the kinetic equation, the limit 
$L \to +\infty$ is finally taken (before the long-time limit, or equivalently, the limit 
$\epsilon \to 0$). As both approaches lead to the same kinetic equation, for simplicity, we anticipate this result and follow the original approach of Benney & Saffman (Reference Benney and Saffman1966). Note that the anisotropic limit (
$k_{\perp } \gg k_{\parallel }$) will also be taken before the (asymptotic) long-time limit.
The first nonlinear term of (2.1) writes
\begin{align} \widehat{({\boldsymbol w}
\boldsymbol{\cdot} \boldsymbol{\nabla}) {\boldsymbol u}}_k
&= {\rm i} \int (\hat {\boldsymbol w}_p \boldsymbol{\cdot}
\boldsymbol{q}) \hat {\boldsymbol u}_q \delta_{k,pq} \,{\rm
d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &=
{\rm i} \int [ {\rm i} \hat \phi_p \hat \phi_q p^2
( \boldsymbol{q} \boldsymbol{\cdot} (\boldsymbol{p}
\times \boldsymbol{e}_\parallel)) (q^2
\boldsymbol{e}_\parallel - q_\parallel \boldsymbol{q}) -
\hat \phi_p \hat \psi_q p^2 ( \boldsymbol{q}
\boldsymbol{\cdot} (\boldsymbol{p} \times
\boldsymbol{e}_\parallel)) (\boldsymbol{q} \times
\boldsymbol{e}_\parallel) \nonumber\\ &\quad -\, \hat \psi_p \hat \phi_q ( p_\parallel
\boldsymbol{p} \boldsymbol{\cdot} \boldsymbol{q} - p^2
q_\parallel ) (q^2 \boldsymbol{e}_\parallel -
q_\parallel \boldsymbol{q}) - {\rm i} \hat \psi_p \hat
\psi_q ( p_\parallel \boldsymbol{p} \boldsymbol{\cdot}
\boldsymbol{q} - p^2 q_\parallel ) (\boldsymbol{q}
\times \boldsymbol{e}_\parallel) ] \nonumber\\ &\quad
\times \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm
d}\boldsymbol{q},
\end{align}
with 
$\delta _{k,pq} \equiv \delta (\boldsymbol {k}-\boldsymbol {p}-\boldsymbol {q})$. In the anisotropic limit (
$k_{\perp } \gg k_{\parallel }$), the following first simplification arises
\begin{align} \widehat{({\boldsymbol w}
\boldsymbol{\cdot} \boldsymbol{\nabla}) {\boldsymbol u}}_k
&= {\rm i} \int [ {\rm i} \hat \phi_p \hat \phi_q
p_{{\perp}}^2 q_{{\perp}}^2 ( \boldsymbol{e}_\parallel
\boldsymbol{\cdot} (\boldsymbol{q_{{\perp}}} \times
\boldsymbol{p_{{\perp}}})) \boldsymbol{e}_\parallel - \hat \phi_p \hat \psi_q p_{{\perp}}^2 (
\boldsymbol{e}_\parallel \boldsymbol{\cdot}
(\boldsymbol{q_{{\perp}}} \times
\boldsymbol{p_{{\perp}}})) (\boldsymbol{q_{{\perp}}}
\times \boldsymbol{e}_\parallel) \nonumber\\ &\quad
- \hat \psi_p \hat \phi_q q_{{\perp}}^2 (
p_\parallel \boldsymbol{p_{{\perp}}} \boldsymbol{\cdot}
\boldsymbol{q_{{\perp}}} - p_{{\perp}}^2 q_\parallel
) \boldsymbol{e}_\parallel - {\rm i} \hat \psi_p
\hat \psi_q ( p_\parallel \boldsymbol{p_{{\perp}}}
\boldsymbol{\cdot} \boldsymbol{q_{{\perp}}} - p_{{\perp}}^2
q_\parallel ) (\boldsymbol{q_{{\perp}}} \times
\boldsymbol{e}_\parallel) ] \nonumber\\ &\quad \times
\delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm
d}\boldsymbol{q} .
\end{align}The second nonlinear term of (2.1) reads
\begin{align} \widehat{({\boldsymbol u}
\boldsymbol{\cdot} \boldsymbol{\nabla}) {\boldsymbol w}}_k
&= {\rm i}\int (\hat {\boldsymbol u}_p \boldsymbol{\cdot}
\boldsymbol{q}) \hat {\boldsymbol w}_q \delta_{k,pq} \,{\rm
d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &=
{\rm i} \int [ {\rm i} \hat \phi_p \hat \phi_q (
p^2 q_\parallel - p_\parallel \boldsymbol{p}
\boldsymbol{\cdot} \boldsymbol{q} ) q^2
(\boldsymbol{q} \times \boldsymbol{e}_\parallel)
\nonumber\\ &\quad + \hat \phi_p \hat \psi_q (
p^2 q_\parallel - p_\parallel \boldsymbol{p}
\boldsymbol{\cdot} \boldsymbol{q} ) (q^2
\boldsymbol{e}_\parallel - q_\parallel \boldsymbol{q})
\nonumber\\ &\quad - \hat \psi_p \hat \phi_q (
\boldsymbol{q} \boldsymbol{\cdot} (\boldsymbol{p} \times
\boldsymbol{e}_\parallel) ) q^2 (\boldsymbol{q}
\times \boldsymbol{e}_\parallel) \nonumber\\ &\quad
+ {\rm i} \hat \psi_p \hat \psi_q (
\boldsymbol{q} \boldsymbol{\cdot} (\boldsymbol{p} \times
\boldsymbol{e}_\parallel) ) (q^2
\boldsymbol{e}_\parallel - q_\parallel \boldsymbol{q})
] \nonumber\\ &\quad \times \delta_{k,pq} \,{\rm
d}\boldsymbol{p} \,{\rm d}\boldsymbol{q},
\end{align}which simplifies in the anisotropic limit to
\begin{align} \widehat{({\boldsymbol u}
\boldsymbol{\cdot} \boldsymbol{\nabla}) {\boldsymbol w}}_k
&= {\rm i} \int q_{{\perp}}^2 [ {\rm i} \hat \phi_p
\hat \phi_q ( p_{{\perp}}^2 q_\parallel - p_\parallel
\boldsymbol{p_{{\perp}}} \boldsymbol{\cdot}
\boldsymbol{q_{{\perp}}} ) (\boldsymbol{q_{{\perp}}}
\times \boldsymbol{e}_\parallel) \nonumber\\ &\quad
+ \hat \phi_p \hat \psi_q ( p_{{\perp}}^2 q_\parallel
- p_\parallel \boldsymbol{p_{{\perp}}} \boldsymbol{\cdot}
\boldsymbol{q_{{\perp}}} ) \boldsymbol{e}_\parallel
\nonumber\\ &\quad - \hat \psi_p \hat \phi_q
( \boldsymbol{e}_\parallel \boldsymbol{\cdot}
(\boldsymbol{q_{{\perp}}} \times \boldsymbol{p_{{\perp}}})
) (\boldsymbol{q_{{\perp}}} \times
\boldsymbol{e}_\parallel) + {\rm i} \hat \psi_p \hat \psi_q
( \boldsymbol{e}_\parallel \boldsymbol{\cdot}
(\boldsymbol{q_{{\perp}}} \times \boldsymbol{p_{{\perp}}})
) \boldsymbol{e}_\parallel ] \nonumber\\ &\quad
\times \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm
d}\boldsymbol{q}.
\end{align}The addition of these two nonlinear contributions leads to the simplified expression
\begin{align}
\widehat{NL}(\boldsymbol{k}) &= \widehat{({\boldsymbol w}
\boldsymbol{\cdot} \boldsymbol{\nabla}) {\boldsymbol u}}_k
- \widehat{({\boldsymbol u} \boldsymbol{\cdot}
\boldsymbol{\nabla}) {\boldsymbol w}}_k \nonumber\\ &=\int
\hat \phi_p \hat \phi_q p_{{\perp}}^2q_{{\perp}}^2 (
\boldsymbol{e}_\parallel \boldsymbol{\cdot}
(\boldsymbol{p_{{\perp}}} \times
\boldsymbol{q_{{\perp}}})) \boldsymbol{e}_\parallel
\delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm
d}\boldsymbol{q} \nonumber\\ &\quad +{\rm i} \int \hat
\phi_p \hat \psi_q p_{{\perp}}^2 (
\boldsymbol{e}_\parallel \boldsymbol{\cdot}
(\boldsymbol{p_{{\perp}}} \times
\boldsymbol{q_{{\perp}}})) (\boldsymbol{q_{{\perp}}}
\times \boldsymbol{e}_\parallel) \delta_{k,pq} \,{\rm
d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &\quad
-{\rm i} \int \hat \psi_p \hat \phi_q q_{{\perp}}^2 (
\boldsymbol{e}_\parallel \boldsymbol{\cdot}
(\boldsymbol{p_{{\perp}}} \times \boldsymbol{q_{{\perp}}})
) (\boldsymbol{q_{{\perp}}} \times
\boldsymbol{e}_\parallel) \delta_{k,pq} \,{\rm
d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &\quad
- \int \hat \psi_p \hat \psi_q q_{{\perp}}^2 (
\boldsymbol{e}_\parallel \boldsymbol{\cdot}
(\boldsymbol{p_{{\perp}}} \times \boldsymbol{q_{{\perp}}})
)\boldsymbol{e}_\parallel \delta_{k,pq} \,{\rm
d}\boldsymbol{p} \,{\rm d}\boldsymbol{q}.
\end{align}The introduction of the canonical variables
\begin{gather} \hat \psi_k = - \frac{1}{2 k_{{\perp}}} \sum_s s A_k^s , \end{gather}
\begin{gather}\hat \phi_k = \frac{1}{2 k_{{\perp}}^2} \sum_s A_k^s , \end{gather}gives
\begin{align}
\widehat{NL}(\boldsymbol{k}) &= \frac{1}{4} \sum_{s_p s_q}
\int A_p^{s_p} A_q^{s_q} ( \boldsymbol{e}_\parallel
\boldsymbol{\cdot} (\boldsymbol{p_{{\perp}}} \times
\boldsymbol{q_{{\perp}}})) \boldsymbol{e}_\parallel
\delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm
d}\boldsymbol{q} \nonumber\\ &\quad -\frac{{\rm i}}{4}
\sum_{s_p s_q} \int A_p^{s_p} A_q^{s_q}
\frac{s_q}{q_{{\perp}}} ( \boldsymbol{e}_\parallel
\boldsymbol{\cdot} (\boldsymbol{p_{{\perp}}} \times
\boldsymbol{q_{{\perp}}})) (\boldsymbol{q_{{\perp}}}
\times \boldsymbol{e}_\parallel) \delta_{k,pq} \,{\rm
d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &\quad
+ \frac{{\rm i}}{4} \sum_{s_p s_q} \int A_p^{s_p} A_q^{s_q}
\frac{s_p}{p_{{\perp}}} ( \boldsymbol{e}_\parallel
\boldsymbol{\cdot} (\boldsymbol{p_{{\perp}}} \times
\boldsymbol{q_{{\perp}}}) ) (\boldsymbol{q_{{\perp}}}
\times \boldsymbol{e}_\parallel) \delta_{k,pq} \,{\rm
d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &\quad
- \frac{1}{4} \sum_{s_p s_q} \int A_p^{s_p} A_q^{s_q} s_p
s_q \frac{q_{{\perp}}}{p_{{\perp}}} (
\boldsymbol{e}_\parallel \boldsymbol{\cdot}
(\boldsymbol{p_{{\perp}}} \times \boldsymbol{q_{{\perp}}})
) \boldsymbol{e}_\parallel \delta_{k,pq} \,{\rm
d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} .
\end{align}
The dummy variables 
$\boldsymbol {p}$, 
$\boldsymbol {q}$ and 
$s_p$, 
$s_q$ can be exchanged to symmetrise the equation; we find that
\begin{align} \widehat{NL}(\boldsymbol{k}) &= \frac{1}{8} \sum_{s_p s_q} \int A_p^{s_p} A_q^{s_q} \frac{\boldsymbol{e}_\parallel \boldsymbol{\cdot} (\boldsymbol{p_{{\perp}}} \times \boldsymbol{q_{{\perp}}})}{p_{{\perp}} q_{{\perp}}} (p_{{\perp}}^2-q_{{\perp}}^2)s_p s_q \boldsymbol{e}_\parallel \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &\quad +\frac{{\rm i}}{8} \sum_{s_p s_q} \int A_p^{s_p} A_q^{s_q} \frac{\boldsymbol{e}_\parallel \boldsymbol{\cdot} (\boldsymbol{p_{{\perp}}} \times \boldsymbol{q_{{\perp}}})}{p_{{\perp}} q_{{\perp}}} ( s_p q_{{\perp}} - s_q p_{{\perp}} )(\boldsymbol{k_{{\perp}}} \times \boldsymbol{e}_\parallel) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{align}Coming back to the wave amplitude equation, we can write
\begin{equation} \left(\frac{\partial \hat \psi_k}{\partial t} - 2 {\rm i} {\boldsymbol{\varOmega}}_0 k_{{\parallel}} \hat \phi_k \right) k_{{\perp}}^2 \boldsymbol{e}_\parallel + \left({\rm i} k_{{\perp}}^2 \frac{\partial \hat \phi_k}{\partial t} + 2 {\boldsymbol{\varOmega}}_0 k_{{\parallel}} \hat \psi_k \right) \boldsymbol{k} \times \boldsymbol{e}_\parallel = \widehat{NL}(\boldsymbol{k}); \end{equation}therefore, after projection and use of the dispersion relation, we obtain
\begin{gather} \frac{\partial \hat \psi_k}{\partial t} - {\rm i} \omega_k k_{{\perp}} \hat \phi_k = \sum_{s_p s_q} \int A_p^{s_p} A_q^{s_q} \frac{\boldsymbol{e}_\parallel \boldsymbol{\cdot} (\boldsymbol{p_{{\perp}}} \times \boldsymbol{q_{{\perp}}})}{8 k_{{\perp}}^2p_{{\perp}} q_{{\perp}}} (p_{{\perp}}^2-q_{{\perp}}^2)s_p s_q \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q}, \end{gather}
\begin{gather}\frac{\partial \hat \phi_k}{\partial t} - {\rm i} \omega_k \frac{\hat \psi_k}{k_{{\perp}}} = \sum_{s_p s_q} \int A_p^{s_p} A_q^{s_q} \frac{\boldsymbol{e}_\parallel \boldsymbol{\cdot} (\boldsymbol{p_{{\perp}}} \times \boldsymbol{q_{{\perp}}})}{8 k_{{\perp}}^2 p_{{\perp}} q_{{\perp}}} ( s_p q_{{\perp}} - s_q p_{{\perp}} ) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{gather}With the introduction of the canonical variables (2.7), the weighted addition of the previous expressions gives
\begin{align} \frac{\partial A^s_k}{\partial t} + {\rm i} s \omega_k A^s_k &= \sum_{s_p s_q} \int s s_p s_q \frac{\boldsymbol{e}_\parallel \boldsymbol{\cdot} (\boldsymbol{p_{{\perp}}} \times \boldsymbol{q_{{\perp}}})}{8 k_{{\perp}} p_{{\perp}} q_{{\perp}}} ( q_{{\perp}}^2-p_{{\perp}}^2 + s k_{{\perp}} ( s_q q_{{\perp}} - s_p p_{{\perp}} ) ) \nonumber\\ &\quad \times A_p^{s_p} A_q^{s_q} \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{align}Remarking that
\begin{equation} q_{{\perp}}^2-p_{{\perp}}^2 = (s_q q_{{\perp}} - s_p p_{{\perp}})(s_pp_{{\perp}} + s_q q_{{\perp}}) , \end{equation}we can rearrange the expression in the following manner:
\begin{align} \frac{\partial A^s_k}{\partial t} + {\rm i} s \omega_k A^s_k &= \sum_{s_p s_q} \int s s_p s_q \frac{\boldsymbol{e}_\parallel \boldsymbol{\cdot} (\boldsymbol{p_{{\perp}}} \times \boldsymbol{q_{{\perp}}})}{8 k_{{\perp}} p_{{\perp}} q_{{\perp}}} ( s_q q_{{\perp}} - s_p p_{{\perp}}) ( s k_{{\perp}} + s_p p_{{\perp}} + s_q q_{{\perp}}) \nonumber\\ &\quad \times A_p^{s_p} A_q^{s_q} \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{align}
We introduce the interaction representation for waves of weak amplitude (
$0 < \epsilon \ll 1$)
\begin{equation} A_k^s = \epsilon a_k^s \exp({-{\rm i}s\omega_k t}), \end{equation}and eventually get the wave amplitude equation after a few last manipulations,
\begin{equation} \frac{\partial a^s_k}{\partial t} = \epsilon \sum_{s_p s_q} \int L_{kpq}^{ss_ps_q} a_p^{s_p} a_q^{s_q} \exp({{\rm i} \varOmega_{k,pq} t}) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q}, \end{equation}
with 
$\varOmega _{k,pq} \equiv s \omega _k - s_p \omega _p - s_q \omega _q$ and
\begin{equation} L_{kpq}^{ss_ps_q} \equiv \omega_k \frac{\boldsymbol{e}_\parallel \boldsymbol{\cdot} (\boldsymbol{p_{{\perp}}} \times \boldsymbol{q_{{\perp}}})}{8 k_{{\perp}} p_{{\perp}} q_{{\perp}}} s_p s_q \left( \frac{s_q q_{{\perp}} - s_p p_{{\perp}}}{s \omega_k}\right) ( s k_{{\perp}} + s_p p_{{\perp}} + s_q q_{{\perp}}) . \end{equation}Expression (2.28) satisfies the following properties (relation (2.11) is used):
\begin{gather} L_{0pq}^{ss_ps_q} = 0, \end{gather}
\begin{gather}L_{kqp}^{ss_qs_p} = L_{kpq}^{ss_ps_q} , \end{gather}
\begin{gather}L_{pkq}^{s_pss_q} = - \frac{s_p \omega_p}{s \omega_k} L_{kpq}^{ss_ps_q} , \end{gather}
\begin{gather}L_{kpq}^{ - s-s_p-s_q} = - L_{kpq}^{ss_ps_q} , \end{gather}
\begin{gather}L_{ - kpq}^{ss_ps_q} = L_{kpq}^{ss_ps_q} , \end{gather}
\begin{gather}L_{ - k-p-q}^{ss_ps_q} = L_{kpq}^{ss_ps_q} . \end{gather}
The wave amplitude equation (2.27) governs the slow evolution of inertial waves of weak amplitude in the anisotropic limit. It is a quadratic nonlinear equation that corresponds to the interactions between waves propagating along 
$\boldsymbol {p}$ and 
$\boldsymbol {q}$, in the positive (
$s_{p},s_{q}>0$) or negative (
$s_{p},s_{q}<0$) direction. The multiple time scale method introduced in the next section is based on this expression. The symmetries listed above will also be used to simplify the derivation of the kinetic equation. Unlike the original derivation by Galtier (Reference Galtier2003), expression (2.27) has been derived without going through a complex helicity basis. The wave amplitude equation tells us that the nonlinear coupling between the states associated with the wavevectors 
$\boldsymbol {p_{\perp }}$ and 
$\boldsymbol {q_{\perp }}$ vanishes when these wavevectors are collinear. Moreover, we note that the nonlinear coupling disappears when the wavenumbers 
$p_{\perp }$ and 
$q_{\perp }$ are equal if their associated directional polarities, 
$s_p$ and 
$s_q$, are also equal. These are general properties for helical waves (Kraichnan Reference Kraichnan1973; Waleffe Reference Waleffe1992; Turner Reference Turner2000; Galtier Reference Galtier2003, Reference Galtier2006, Reference Galtier2014).
3. Asymptotic sequential closures
The method outlined here was first proposed by Benney & Saffman (Reference Benney and Saffman1966) for three-wave interactions, but to our knowledge it has never been explicitly applied to a physical system. In fact, originally the main motivation for such a development was four-wave interactions that describe gravity waves (Hasselmann Reference Hasselmann1962; Newell Reference Newell1968), and for which the main prediction of wave turbulence (Kolmogorov–Zakharov spectrum) is now well observed (Lenain & Melville Reference Lenain and Melville2017).
Unlike the classical perturbation analysis, the multiple time scale method is based on the existence of a sequence of time scales, 
$T_{0}$, 
$T_{1}$, 
$T_{2},\ldots$, with by definition (Nayfeh Reference Nayfeh2004)
\begin{equation} T_{0} \equiv t,\quad T_{1}\equiv \epsilon t, \quad T_{2}\equiv \epsilon^{2} t,\ldots. \end{equation}
Because of the weak dependence in 
$t$ of 
$T_{1}$, 
$T_{2},\ldots$, all these variables will be treated (it is an approximation) as independent. Furthermore, the smaller 
$\epsilon$ is, the better the approximation. The variation of the wave amplitude with 
$T_1$ and 
$T_2$ represents the slow variation that we wish to extract. (In practice, our analysis will be limited to 
${O} (\epsilon ^2)$.) Using the chain rule, we obtain (
$T_0$ being replaced by 
$t$)
\begin{equation} \left(\frac{\partial}{\partial t} + \epsilon \frac{\partial}{\partial T_{1}} + \epsilon^{2} \frac{\partial}{\partial T_{2}} + \cdots\right) a^{s}_{k} = \epsilon \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} a^{s_{p}}_{p} a^{s_{q}}_{q} \exp({{\rm i} \varOmega_{k,pq}t}) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{equation}
The variable 
$a^{s}_{k}$ must also be expanded to the power of 
$\epsilon$ to make the various scales appear in time,
\begin{equation} a^{s}_{k} = \sum_{n=0}^{+\infty} \epsilon^n a^{s}_{k,n}(t,T_{1},T_{2},\ldots) = a^{s}_{k,0} + \epsilon a^{s}_{k,1} + \epsilon^{2} a^{s}_{k,2} + \cdots . \end{equation}Expression (3.3) is then introduced into the fundamental equation (3.2); we obtain, for the first three terms,
\begin{gather} \frac{\partial a^{s}_{k,0} }{\partial t} = 0 , \end{gather}
\begin{gather}\frac{\partial a^{s}_{k,1} }{\partial t} = - \frac{\partial a^{s}_{k,0} }{\partial T_{1}} + \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} a_{p,0}^{s_{p}} a_{q,0}^{s_{q}} \exp({{\rm i} \varOmega_{k,pq}t}) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} , \end{gather}
\begin{gather}\frac{\partial a^{s}_{k,2} }{\partial t} = - \frac{\partial a^{s}_{k,1} }{\partial T_{1}} - \frac{\partial a^{s}_{k,0} }{\partial T_{2}} + \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} \left[a_{p,1}^{s_{p}} a_{q,0}^{s_{q}} + a_{p,0}^{s_{p}} a_{q,1}^{s_{q}} \right] \exp({{\rm i} \varOmega_{k,pq}t}) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{gather}
To lighten the writing, the time dependency of the variables has been omitted. After integration on 
$t$, one finds that
\begin{gather} a^{s}_{k,0} = a^{s}_{k,0}(T_1,T_2,\ldots) , \end{gather}
\begin{gather}a^{s}_{k,1} = - t \frac{\partial a^{s}_{k,0} }{\partial T_{1}} + b^{s}_{k,1} , \end{gather}
\begin{gather}a^{s}_{k,2} = \frac{t^2}{2} \frac{\partial^2 a^{s}_{k,0} }{\partial T_{1}^2} -t \frac{\partial a^{s}_{k,0} }{\partial T_{2}} - \int_0^t \frac{\partial b^{s}_{k,1} }{\partial T_{1}} {\rm d} t + \tilde b^{s}_{k,2}, \end{gather}with by definition
\begin{gather} b^{s}_{k,1} \equiv \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} a_{p,0}^{s_{p}} a_{q,0}^{s_{q}} {\rm \Delta}(\varOmega_{k,pq}) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q}, \end{gather}
\begin{gather}{\rm \Delta}(X) \equiv \int_{0}^{t} {\rm e}^{{\rm i} Xt} \,{\rm d} t = \frac{{\rm e}^{{\rm i} Xt} -1}{{\rm i} X}, \end{gather}and
\begin{equation} \tilde b^{s}_{k,2} \equiv \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} \int_{0}^{t} \left[ a_{p,1}^{s_{p}} a_{q,0}^{s_{q}} + a_{p,0}^{s_{p}} a_{q,1}^{s_{q}}\right] \exp({{\rm i} \varOmega_{k,pq}t}) \,{\rm d} t \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{equation}The previous equation is modified when expression (3.5b) is introduced; one finds that
\begin{align} \tilde b^{s}_{k,2} &= - \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} \frac{\partial (a^{s_p}_{p,0} a_{q,0}^{s_{q}})}{\partial T_{1}} \left( \int_0^t t \exp({{\rm i} \varOmega_{k,pq}t})\,{\rm d} t \right) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &\quad + \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} \left( \int_{0}^{t} (b^{s_p}_{p,1} a_{q,0}^{s_{q}} + b^{s_q}_{q,1} a_{p,0}^{s_{p}} ) \exp({{\rm i} \varOmega_{k,pq}t})\,{\rm d} t \right)\delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{align}Expression (3.5c) becomes
\begin{equation} a^{s}_{k,2} = \frac{t^2}{2} \frac{\partial^2 a^{s}_{k,0}}{\partial T_{1}^2} -t \frac{\partial a^{s}_{k,0}}{\partial T_{2}} + b^{s}_{k,2} , \end{equation}with
\begin{align} b^{s}_{k,2} &= - \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} \frac{\partial (a^{s_p}_{p,0} a_{q,0}^{s_{q}})}{\partial T_{1}} \left( \int_0^t \left[{\rm \Delta}( \varOmega_{k,pq}) + t \exp({{\rm i} \varOmega_{k,pq}t}) \right] {\rm d} t \right) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &\quad + \sum_{s_{p} s_{q} s_{p'} s_{q'} } \int 2 L^{s s_p s_q}_{kpq} L^{s_p s_{p'} s_{q'}}_{pp'q'} a_{p',0}^{s_{p'}} a_{q',0}^{s_{q'}} a_{q,0}^{s_{q}} \left( \int_{0}^{t} {\rm \Delta}( \varOmega_{p,p'q'}) \exp({{\rm i} \varOmega_{k,pq}t})\,{\rm d} t \right) \nonumber\\ &\quad \times \delta_{k,pq} \delta_{p,p'q'} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \,{\rm d}\boldsymbol{p}' \,{\rm d}\boldsymbol{q}' . \end{align}The time integrals give the relations
\begin{gather} \int_0^t \left[{\rm \Delta}( \varOmega_{k,pq}) + t \exp({{\rm i} \varOmega_{k,pq}t}) \right] {\rm d} t = t {\rm \Delta}( \varOmega_{k,pq}) , \end{gather}
\begin{gather}\int_{0}^{t} {\rm \Delta}( \varOmega_{p,p'q'}) \exp({{\rm i} \varOmega_{k,pq}t})\,{\rm d} t = \frac{{\rm \Delta}(\varOmega_{k,p'q'q}) - {\rm \Delta}(\varOmega_{k,pq})}{{\rm i} (\varOmega_{k,p'q'q}-\varOmega_{k,pq})}, \end{gather}that will be used below in the long-time limit.
3.1. First asymptotic closure at time 
$T_1$
With the previous definitions, the perturbative expansion of the second-order moment writes
\begin{align} \langle a^{s}_{k} a^{s'}_{k'} \rangle &= \langle (a^{s}_{k,0} + \epsilon a^{s}_{k,1} + \epsilon^2 a^{s}_{k,2} + \cdots)(a^{s'}_{k',0} + \epsilon a^{s'}_{k',1} + \epsilon^2 a^{s'}_{k',2} + \cdots) \rangle \nonumber\\ &= \langle a^{s}_{k,0} a^{s'}_{k',0} \rangle + \epsilon \langle a^{s}_{k,0} a^{s'}_{k',1} + a^{s}_{k,1} a^{s'}_{k',0} \rangle \nonumber\\ &\quad + \epsilon^2 \langle a^{s}_{k,0} a^{s'}_{k',2} + a^{s}_{k,1} a^{s'}_{k',1} + a^{s}_{k,2} a^{s'}_{k',0} \rangle + \cdots, \end{align}
where 
$\langle \rangle$ denotes the ensemble average. We shall assume that this turbulence is statistically homogeneous. In this case, the second-order moment can be written in terms of the second-order cumulant, 
$q^{ss'}(\boldsymbol {k},\boldsymbol {k}') \equiv q_{k}^{ss'}$, such that
\begin{equation} \langle a^{s}_{k} a^{s'}_{k'} \rangle = q_{k}^{ss'} \delta (\boldsymbol{k}+\boldsymbol{k}') , \end{equation}
where the presence of 
$\delta (\boldsymbol {k}+\boldsymbol {k}')$ is the consequence of the statistical homogeneity (Galtier Reference Galtier2023b). We have to adapt the original formalism developed for isotropic problems (Benney & Saffman Reference Benney and Saffman1966) to this anisotropic case where the dispersion relation depends not only on the wavenumber 
$k$ but also on the component 
$k_{\parallel }$. In this case, a contribution from 
$q^{ss'}$ is only relevant if 
$s=s'$, whereas it is for 
$s=-s'$ in the case of an isotropic problem. (With such conditions 
$q^{ss'}_{k}$ is real.) We assume – and this is the basic idea of the method – that the second-order moment (in fact, the coefficient 
$q^{ss'}_{k}$ in front of the delta function) on the left-hand side of expression (3.14) remains bounded at all time (Benney & Saffman Reference Benney and Saffman1966). As an example, we can think of the energy spectrum that, as we know, remains physically bounded. Therefore, the contributions on the right-hand side must also be bounded at each order in 
$\epsilon$. We will see that secular terms can appear at different orders in 
$\epsilon$; this leads to certain conditions to cancel their contributions to keep the development uniform in time. As we shall see, at order 
${O} (\epsilon ^2)$ this condition leads to the so-called kinetic equations. The main problem is therefore to account for all the secular contributions.
At order 
${O} (\epsilon ^0)$, we have the contribution of 
$\langle a^{s}_{k,0} a^{s'}_{k',0} \rangle$ that will therefore be assumed to be bounded at all times.
At order 
${O} (\epsilon ^1)$, we have the contribution
\begin{align} \langle a^{s}_{k,0} a^{s'}_{k',1} + a^{s}_{k,1} a^{s'}_{k',0} \rangle &= \left\langle a^{s}_{k,0} \left( - t \frac{\partial a^{s'}_{k',0}}{\partial T_{1}} + b^{s'}_{k',1}\right) + \left( - t \frac{\partial a^{s}_{k,0}}{\partial T_{1}} + b^{s}_{k,1}\right) a^{s'}_{k',0} \right\rangle \nonumber\\ &= - t \frac{\partial}{\partial T_{1}} \langle a^{s}_{k,0} a^{s'}_{k',0} \rangle + \langle a^{s}_{k,0} b^{s'}_{k',1} \rangle + \langle b^{s}_{k,1} a^{s'}_{k',0} \rangle . \end{align}
The first term on the right-hand side gives a secular contribution proportional to 
$t$. For the second term, we have
\begin{equation} \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} \langle a^{s}_{k,0} a_{p,0}^{s_{p}} a_{q,0}^{s_{q}} \rangle {\rm \Delta}(\varOmega_{k,pq}) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{equation}
The long-time limit (
$t \gg 1/\omega$) of this oscillating integral will be given by the Riemann–Lebesgue lemma (the proof requires the use of generalized functions)
\begin{equation} {\rm \Delta}(X) = \frac{{\rm e}^{{\rm i} Xt} -1}{{\rm i} X} \xrightarrow{\text{t $\to +\infty$}} {\rm \pi}\delta(X) + {\rm i} \mathcal{P} \left(\frac{1}{X}\right) , \end{equation}
where 
$\mathcal {P}$ is the Cauchy principal value of the integral. Therefore, the long-time limit of expression (3.17) gives no secular contribution. The same conclusion is obtained for the third term of expression (3.16). Also, the condition to cancel the unique secular term is
\begin{equation} \frac{\partial \langle a^{s}_{k,0} a^{s'}_{k',0} \rangle }{\partial T_{1}} = 0 , \end{equation}
which means that the second-order moment does not evolve over a time scale 
$T_1$. As will be seen later, a turbulent cascade in inertial wave turbulence is only expected on a time scale 
$T_2$. (Here we have a point of disagreement with expression (2.43) in Benney & Saffman (Reference Benney and Saffman1966): it is not correct, but this has no impact on the rest of the paper.)
3.2. Second asymptotic closure at time $T_2$
The analysis continues at order 
${O} (\epsilon ^2)$. With expression (3.10), the next contribution reads
\begin{align} & \langle a^{s}_{k,1}
a^{s'}_{k',1} + a^{s}_{k,0} a^{s'}_{k',2} + a^{s}_{k,2}
a^{s'}_{k',0} \rangle =\left\langle \left( - t
\frac{\partial a^{s}_{k,0}}{\partial T_{1}} +
b^{s}_{k,1}\right) \left( - t \frac{\partial
a^{s'}_{k',0}}{\partial T_{1}} + b^{s'}_{k',1}\right)
\right\rangle \nonumber\\ &\quad + \left\langle a^{s}_{k,0}
\left(\frac{t^2}{2} \frac{\partial^2
a^{s'}_{k',0}}{\partial T_{1}^2} -t \frac{\partial
a^{s'}_{k',0}}{\partial T_{2}}+ b^{s'}_{k',2} \right) +
a^{s'}_{k',0} \left(\frac{t^2}{2} \frac{\partial^2
a^{s}_{k,0}}{\partial T_{1}^2} -t \frac{\partial
a^{s}_{k,0}}{\partial T_{2}} + b^{s}_{k,2} \right)
\right\rangle,
\end{align}which gives after development and simplifications
\begin{align} \langle a^{s}_{k,1}
a^{s'}_{k',1} + a^{s}_{k,0} a^{s'}_{k',2} + a^{s}_{k,2}
a^{s'}_{k',0} \rangle &= \frac{t^2}{2}
\frac{\partial^2 \langle a^{s}_{k,0} a^{s'}_{k',0} \rangle
}{\partial T_{1}^2} - t \frac{\partial \langle a^{s}_{k,0}
a^{s'}_{k',0} \rangle }{\partial T_{2}} + \langle
b^{s}_{k,1} b^{s'}_{k',1} \rangle \nonumber\\ &\quad
-t \left\langle \frac{\partial a^{s}_{k,0}}{\partial T_{1}}
b^{s'}_{k',1} + b^{s}_{k,1} \frac{\partial
a^{s'}_{k',0}}{\partial T_{1}} \right\rangle \nonumber\\
&\quad + \langle a^{s}_{k,0} b^{s'}_{k',2} +
a^{s'}_{k',0} b^{s}_{k,2} \rangle .
\end{align}The first term on the right-hand side cancels over the long time as required by the first asymptotic closure. The second term gives a secular contribution. The other three terms can potentially give a secular contribution: it is obvious for the fourth term and non-trivial for the third and fifth terms that require further development.
The third term on the right-hand side writes
\begin{align} \left\langle b^{s}_{k,1} b^{s'}_{k',1} \right\rangle &= \sum_{s_{p} s_{q} s_{p'} s_{q'}} \int L^{s s_p s_q}_{kpq} L^{s' s_{p'} s_{q'}}_{k'p'q'} \langle a_{p,0}^{s_{p}} a_{q,0}^{s_{q}} a_{p',0}^{s_{p'}} a_{q',0}^{s_{q'}} \rangle {\rm \Delta}(\varOmega_{k,pq}) {\rm \Delta}(\varOmega_{k',p'q'}) \nonumber\\ &\quad \times \delta_{k,pq} \delta_{k',p'q'} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \,{\rm d}\boldsymbol{p}' \,{\rm d}\boldsymbol{q}' . \end{align}Here again, the theory of generalized functions gives us the long-time behaviour of this oscillating integral (with the Poincaré–Bertrand formula – see, e.g. Benney & Newell Reference Benney and Newell1969),
\begin{equation} {\rm \Delta}(X) {\rm \Delta}( - X) \xrightarrow{\text{t $\to +\infty$}} 2 {\rm \pi}t \delta(X) + 2 \mathcal{P} \left(\frac{1}{X}\right) \frac{\partial}{\partial X}. \end{equation}
Therefore, a secular contribution (proportional to 
$t$) involving a Dirac function is possible. The fourth-order moment in expression (3.22) can be decomposed into products of second-order cumulant plus a fourth-order cumulant such that (the statistical homogeneity is used as well as 
$\langle a_{k,0}^{s}\rangle = 0$)
\begin{align} & \langle a_{p,0}^{s_{p}} a_{q,0}^{s_{q}} a_{p',0}^{s_{p'}} a_{q',0}^{s_{q'}} \rangle = q^{s_{p}s_{q}s_{p'}s_{q'}}_{pqp',0} \delta(\boldsymbol{p}+\boldsymbol{q}+\boldsymbol{p}'+\boldsymbol{q}') + q^{s_{p} s_q}_{p,0} q^{s_{p'} s_{q'}}_{p',0} \delta(\boldsymbol{p}+\boldsymbol{q}) \delta(\boldsymbol{p}'+\boldsymbol{q}') \nonumber\\ &\quad + q^{s_{p} s_{p'}}_{p,0} q^{s_{q} s_{q'}}_{q,0} \delta(\boldsymbol{p}+\boldsymbol{p}') \delta(\boldsymbol{q}+\boldsymbol{q}') + q^{s_{p} s_{q'}}_{p,0} q^{s_{q} s_{p'}}_{q,0} \delta(\boldsymbol{p}+\boldsymbol{q}') \delta(\boldsymbol{q}+\boldsymbol{p}') . \end{align}
Note that according to expression (3.14) by homogeneity we also have the relation 
$\boldsymbol {k}=-\boldsymbol {k}'$. We obtain
\begin{align} & \langle b^{s}_{k,1} b^{s_{k'}}_{k',1} \rangle = \sum_{s_{p} s_{q} s_{p'} s_{q'}} \int L^{s s_p s_q}_{kpq} L^{s' s_{p'} s_{q'}}_{k'p'q'} \left[ q^{s_{p}s_{q}s_{p'}s_{q'}}_{pqp',0} \delta(\boldsymbol{p}+\boldsymbol{q}+\boldsymbol{p}'+\boldsymbol{q}') \right. \nonumber\\ &\quad + q^{s_{p} s_q}_{p,0} q^{s_{p'} s_{q'}}_{p',0} \delta(\boldsymbol{p}+\boldsymbol{q}) \delta(\boldsymbol{p}'+\boldsymbol{q}') + q^{s_{p} s_{p'}}_{p,0} q^{s_{q} s_{q'}}_{q,0} \delta(\boldsymbol{p}+\boldsymbol{p}') \delta(\boldsymbol{q}+\boldsymbol{q}') \nonumber\\ &\quad \left. +\, q^{s_{p} s_{q'}}_{p,0} q^{s_{q} s_{p'}}_{q,0} \delta(\boldsymbol{p}+\boldsymbol{q}') \delta(\boldsymbol{q}+\boldsymbol{p}') \right] {\rm \Delta}(\varOmega_{k,pq}) {\rm \Delta}(\varOmega_{k',p'q'}) \delta_{k,pq} \delta_{k',p'q'} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \,{\rm d}\boldsymbol{p}' \,{\rm d}\boldsymbol{q}' . \end{align}
We are looking for secular contributions. In the second line the first term does not contribute since it imposes 
$\boldsymbol {k}=\textbf{0}$ that cancels 
$L^{s s_p s_q}_{kpq}$, but the second term can contribute (in this anisotropic problem) when the conditions 
$s_p=s_{p'}$ and 
$s_q=s_{q'}$ are satisfied. Likewise, in the third line a contribution is possible when 
$s_p=s_{q'}$ and 
$s_q=s_{p'}$. There is no contribution from the first line, which means that the situation is the same as if the distribution were Gaussian (however, we do not make this assumption). In summary, in the long-time limit the secular contribution, written 
$\mathcal {C}_t \langle b^{s}_{k,1} b^{s'}_{k',1} \rangle$, is
\begin{align} \mathcal{C}_t \langle b^{s}_{k,1} b^{s'}_{k',1} \rangle &= 4 {\rm \pi}t \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} L^{s s_p s_q}_{k-p-q} q^{s_{p} s_{p}}_{p,0} q^{s_{q} s_{q}}_{q,0} \delta (\varOmega_{k,pq}) \delta_{k,pq} \delta_{kk'} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &= 4 {\rm \pi}t \sum_{s_{p} s_{q}} \int \vert L^{s s_p s_q}_{kpq} \vert^2 q^{s_{p} s_{p}}_{p,0} q^{s_{q} s_{q}}_{q,0} \delta (\varOmega_{k,pq}) \delta_{k,pq} \delta_{kk'} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{align} The fourth term on the right-hand side of (3.21) does not contribute over the long time because it depends on the 
$T_1$ derivative. The proof is given by a new relation involving the 
$n$-order moments that can be written as
\begin{align} \langle a_k^{s} a_{k'}^{s'} a_{k''}^{s''} \cdots\rangle &= \langle a_{k,0}^{s} a_{k',0}^{s'} a_{k'',0}^{s''} \cdots\rangle \nonumber\\ &\quad + \epsilon \langle a_{k,1}^{s} a_{k',0}^{s'} a_{k'',0}^{s''} \cdots+ a_{k,0}^{s} a_{k',1}^{s'} a_{k'',0}^{s''} \cdots+ a_{k,0}^{s} a_{k',0}^{s'} a_{k'',1}^{s''} \cdots. + \cdots.\rangle \nonumber\\ &\quad + \epsilon^2 \langle \cdots. \rangle + \cdots . \end{align}
As for the second-order moment, we demand that the moments (i.e. the coefficients in front of the delta functions) of order 
$n$ are bounded. At order 
${O} (\epsilon )$, we obtain the relation
\begin{align} & \langle a_{k,1}^{s} a_{k',0}^{s'} a_{k'',0}^{s''} \cdots+ a_{k,0}^{s} a_{k',1}^{s'} a_{k'',0}^{s''} \cdots+ a_{k,0}^{s} a_{k',0}^{s'} a_{k'',1}^{s''} \cdots+ \cdots\rangle \nonumber\\ &\quad = \left\langle \left( - t \frac{\partial a^{s}_{k,0}}{\partial T_{1}} + b^{s}_{k,1}\right) a_{k',0}^{s'} a_{k'',0}^{s''} \cdots+ a_{k,0}^{s} \left( - t \frac{\partial a^{s'}_{k',0}}{\partial T_{1}} + b^{s'}_{k',1}\right) a_{k'',0}^{s''} \cdots+ \cdots\right\rangle \nonumber\\ &\quad = - t \frac{\partial \langle a^{s}_{k,0} a^{s'}_{k',0} a_{k'',0}^{s''} \cdots\rangle}{\partial T_{1}} + \langle b^{s}_{k,1}a_{k',0}^{s'} a_{k'',0}^{s''} \cdots\rangle + \langle a^{s}_{k,0} b^{s'}_{k',1} a_{k'',0}^{s''} \cdots\rangle + \cdots. \end{align}Only the first term of the last line gives a secular contribution over the long time, which means that we have to impose the asymptotic condition
\begin{equation} \frac{\partial \langle a^{s}_{k,0} a^{s'}_{k',0} a_{k'',0}^{s''} \cdots\rangle}{\partial T_{1}} = 0, \end{equation}
at any order 
$n$. Therefore, the probability density function does not depend on 
$T_1$ and we can assume that the variable itself does not depend on 
$T_1$. (It is a mild hypothesis because it is difficult to imagine a turbulent system where everything fluctuates, and in which it would be possible to have a 
$T_1$ dependence for the amplitude whereas the probability density function has no such dependence.) This shows that the fourth term on the right-hand side of (3.21) does not contribute in the long time.
The last term of (3.21) writes
\begin{align} & \langle a^{s}_{k,0}
b^{s'}_{k',2} + a^{s'}_{k',0} b^{s}_{k,2} \rangle =
\sum_{s_{p} s_{q} s_{p'} s_{q'} } \int 2 L^{s' s_p
s_q}_{k'pq} L^{s_p s_{p'} s_{q'}}_{pp'q'} \langle
a^{s}_{k,0} a_{p',0}^{s_{p'}} a_{q',0}^{s_{q'}}
a_{q,0}^{s_{q}} \rangle \nonumber\\ &\quad \times
\left( \frac{{\rm \Delta}(\varOmega_{k',p'q'q}) -
{\rm \Delta}(\varOmega_{k',pq})}{{\rm i}
(\varOmega_{k',p'q'q}-\varOmega_{k',pq})} \right)
\delta_{k',pq} \delta_{p,p'q'} \,{\rm d}\boldsymbol{p}
\,{\rm d}\boldsymbol{q} \,{\rm d}\boldsymbol{p}' \,{\rm
d}\boldsymbol{q}' \nonumber\\ &\quad +\sum_{s_{p} s_{q}
s_{p'} s_{q'} } \int 2 L^{s s_p s_q}_{kpq} L^{s_p s_{p'}
s_{q'}}_{pp'q'} \langle a^{s'}_{k',0}
a_{p',0}^{s_{p'}} a_{q',0}^{s_{q'}} a_{q,0}^{s_{q}}
\rangle \left(\frac{{\rm \Delta}(\varOmega_{k,p'q'q}) -
{\rm \Delta}(\varOmega_{k,pq})}{{\rm i}
(\varOmega_{k,p'q'q}-\varOmega_{k,pq})} \right) \nonumber\\
&\quad \times \delta_{k,pq} \delta_{p,p'q'} \,{\rm
d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \,{\rm
d}\boldsymbol{p}' \,{\rm d}\boldsymbol{q}' .
\end{align}The secular contributions will be given by the theory of generalized functions with the relation
\begin{equation} \frac{{\rm \Delta}(X) - {\rm \Delta}(0)}{{\rm i}X} \xrightarrow{\text{t $\to +\infty$}} {\rm \pi}t \delta(X) + {\rm i}t \mathcal{P} \left(\frac{1}{X}\right). \end{equation}
We also need to use the following development (and its symmetric in 
$\boldsymbol {k}'$):
\begin{align} & \langle a_{k,0}^{s} a_{p',0}^{s_{p'}} a_{q',0}^{s_{q'}} a_{q,0}^{s_{q}} \rangle = q^{ss_{p'}s_{q'}s_{q}}_{kp'q',0} \delta(\boldsymbol{k}+\boldsymbol{p}'+\boldsymbol{q}'+\boldsymbol{q}) + q^{s s_{q}}_{k,0} q^{s_{p'} s_{q'}}_{p',0} \delta(\boldsymbol{k}+\boldsymbol{q}) \delta(\boldsymbol{p}'+\boldsymbol{q}') \nonumber\\ &\quad + q^{s s_{q'}}_{k,0} q^{s_{p'} s_{q}}_{p',0} \delta(\boldsymbol{k}+\boldsymbol{q}') \delta(\boldsymbol{p}'+\boldsymbol{q}) + q^{s s_{p'}}_{k,0} q^{s_{q'} s_{q}}_{q',0} \delta(\boldsymbol{k}+\boldsymbol{p}') \delta(\boldsymbol{q}'+\boldsymbol{q}) . \end{align}
On the right-hand side of expression (3.32) the first two terms do not give a secular contribution; however, the last two terms give a contribution when the following conditions are satisfied, namely 
$s=s_{q'}$, 
$s_{p'}=s_q$ and 
$s=s_{p'}$, 
$s_{q'}=s_q$, respectively. After substitution and simplification, we obtain the secular contributions in the long-time limit
\begin{align} & \mathcal{C}_t \langle a^{s}_{k,0} b^{s'}_{k',2} + a^{s'}_{k',0} b^{s}_{k,2} \rangle \nonumber\\ &\quad = + 4t \sum_{s_{p} s_{q}} \int L^{s' s_p s_q}_{k'pq} L^{s_p s_{q} s'}_{p-qk'} q^{s' s'}_{k',0} q^{s_{q} s_{q}}_{q,0} \left( {\rm \pi}\delta(\varOmega_{k',pq}) + {\rm i} \mathcal{P} \left(\frac{1}{\varOmega_{k',pq}}\right) \right) \delta_{k',pq} \delta_{kk'} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &\qquad +4t \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} L^{s_p s_{q} s}_{p-qk} q^{s s}_{k,0} q^{s_{q} s_{q}}_{q,0} \left( {\rm \pi}\delta(\varOmega_{k,pq}) + {\rm i} \mathcal{P} \left(\frac{1}{\varOmega_{k,pq}}\right) \right) \delta_{k,pq} \delta_{kk'} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &\quad = 8 {\rm \pi}t \sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} L^{s_p s_{q} s}_{p-qk} q^{s s}_{k,0} q^{s_{q} s_{q}}_{q,0} \delta(\varOmega_{k,pq}) \delta_{k,pq} \delta_{kk'} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q}. \end{align}
The last writing is obtained by using the general property 
$L^{s s_p s_q}_{-k-p-q} = L^{s s_p s_q}_{kpq}$ and the symmetry in 
$\boldsymbol {p} \to -\boldsymbol {p}$ and 
$\boldsymbol {q} \to -\boldsymbol {q}$.
If in expression (3.21) we impose the nullity of the sum of the different secular contributions, we find the asymptotic condition (after integration over 
$k'$)
\begin{align} \frac{\partial q^{ss'}_{k,0}}{\partial T_{2}} &= 4 {\rm \pi}\sum_{s_{p} s_{q}} \int \vert L^{s s_p s_q}_{kpq} \vert^2 q^{s_{p} s_{p}}_{p,0} q^{s_{q} s_{q}}_{q,0} \delta (\varOmega_{k,pq}) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} \nonumber\\ &\quad + 8 {\rm \pi}\sum_{s_{p} s_{q}} \int L^{s s_p s_q}_{kpq} L^{s_p s_q s}_{p-qk} q^{s s}_{k,0} q^{s_{q} s_{q}}_{q,0} \delta(\varOmega_{k,pq}) \delta_{k,pq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} . \end{align}
Introducing 
$e^{s}_k \equiv q^{ss}_{k,0}$, we end up with (some simple manipulations are also used to symmetrise the equation)
\begin{align} \frac{\partial e_k^{s}}{\partial t} &= \frac{{\rm \pi} \epsilon^2}{16 s\omega_k} \sum_{s_{p} s_{q}} \int \left(\frac{\sin \theta_k}{k_{{\perp}}} \right)^2 (s_p p_{{\perp}} - s_q q_{{\perp}})^2 ( s k_{{\perp}} + s_p p_{{\perp}} + s_q q_{{\perp}})^2 \nonumber\\ &\quad \times \left[ s\omega_k e_p^{s_p} e_q^{s_q} + s_p \omega_p e_k^{s} e_q^{s_q} + s_q \omega_q e_k^{s} e_p^{s_p} \right] \delta (\varOmega_{kpq}) \delta_{kpq} \,{\rm d}\boldsymbol{p} \,{\rm d}\boldsymbol{q} , \end{align}
with 
$\theta _k$ the opposite angle to 
$\boldsymbol {k_{\perp }}$ in the triangle 
$\boldsymbol {k_{\perp }}=\boldsymbol {p_{\perp }}+\boldsymbol {q_{\perp }}$. Expression (3.35) is the kinetic equation for inertial wave turbulence in the anisotropic limit (see (7) in Galtier (Reference Galtier2003); the small difference (sign and numerical factor) depends only on the normalisation of the canonical variables (2.7)). This result proves that the anisotropic and asymptotic limits commute. Note that the kinetic equation for inertial wave turbulence does not describe the slow mode (
$k_{\parallel }=0$) that involves strong turbulence.
4. Discussion and conclusion
Our study completes the original derivation of Galtier (Reference Galtier2003) where the delicate issue of the asymptotic limit was mentioned but not explicitly used. We obtained the kinetic equation using a multiple time scale method that leads to sequential asymptotic closures at times 
$T_1$ and 
$T_2$. This method consists of imposing the nullity of secular terms (proportional to 
$t$) that emerge over asymptotically long times in order to guarantee a bounded value for the associated moments and, thus, keep the development uniform in time. These secular terms do not involve fourth-order cumulants but only products of second-order cumulants. Consequently, the derivation of the kinetic equation (3.35) is performed in a systematic and consistent manner, and does not require any closure assumptions such as quasi-Gaussianity. In addition, the random phase approximation is natural as it arises dynamically from the separation of time scales. The latter property was first mentioned by Benney & Saffman (Reference Benney and Saffman1966), but has not always been recognized as an inherent property of wave turbulence (Zakharov, L'Vov & Falkovich Reference Zakharov, L'Vov and Falkovich1992; Nazarenko Reference Nazarenko2011). Generally, the random phase approximation is introduced (initially or subsequently) to fully justify the closure. This fundamental difference with the multiple time scale method suggests that the latter is the most natural method for deriving the kinetic equation of wave turbulence. Note, however, that the derivation made, although systematic, says nothing about the remaining terms (not used to obtain the kinetic equation) in the small 
$\epsilon$ perturbation expansion, the implicit conjecture being that they are subdominant. The proof of this conjecture remains a mathematical challenge (Deng & Zaher Reference Deng and Zaher2021).
The domain of validity of wave turbulence has already been discussed in Galtier (Reference Galtier2003). The conclusion is that we can always find a finite domain where inertial wave turbulence exists. To show this, we introduce the linear time 
$\tau _W \sim 1/ \omega _k \sim k/(k_\parallel \varOmega _0)$ and the nonlinear time 
$\tau _{NL} \sim 1/(k u)$. In the anisotropic limit we have 
$k \sim k_\perp$, and we obtain the time ratio
\begin{equation} \chi = \frac{\tau_W}{\tau_{NL}} \sim \frac{k_\perp^2 u}{k_\parallel \varOmega_0}, \end{equation}
which depends on the scale. If the system under study is initially excited locally at large scale, isotropically and with a Rossby number much smaller than 
$1$, then we can obtain 
$\chi \ll 1$ at large scale, which is the condition for having weak turbulence. As we have explained, such turbulence becomes anisotropic, with the energy cascading towards a region of Fourier space where 
$k_\perp \gg k_\parallel$. This leads to an increase of 
$\chi$ (the dependence of 
$u$ on wavenumbers does not alter this trend). Clearly, at sufficiently large values of 
$k_\perp$, we can find 
$\chi \sim 1$, which is synonymous with strong turbulence (and the critical balance regime). There is therefore a domain in Fourier space where the condition for the validity of wave turbulence can be satisfied even in the anisotropic limit.
In summary, it can be said that inertial wave turbulence for three-wave interactions is characterised by a dynamics on two time scales. On short time scales, of the order of the wave period, there is phase mixing that leads, due to the dispersive nature of the waves, to the decoupling of the correlations if they are initially present and to a statistics close to Gaussianity, as expected from the central limit theorem. This happens with a decay in 
$1/t$. On a longer time scale, the nonlinear coupling – weak at short times – becomes non-negligible due to the resonance mechanism. This coupling leads to a regeneration of the cumulants via the product of lower order cumulants. It is these terms that are at the origin of the energy transfer mechanism.
The second novelty of our study is the demonstration that anisotropy and asymptotic closure commute. (Note that this property was also found in Alfvén wave turbulence Galtier et al. Reference Galtier, Nazarenko, Newell and Pouquet2002.) Therefore, to obtain the main properties of inertial wave turbulence (exact power law solution, direction of the cascade, Kolmogorov constant, existence of an inertial range), the limit 
$k_{\perp } \gg k_{\parallel }$ can be taken before any statistical development. The study also reveals that it is not necessary to use an helicity basis, which simplifies the treatment. The third novelty is that the multiple time scale method has been generalized to an anisotropic problem involving different types of correlation in terms of directional polarity.
Note that in the present derivation, the system under study is assumed to be of infinite size and can therefore be treated as continuous. The numerical simulation with its grid of points escapes this description. Effects (freezing of the cascade) linked to the discretisation of the Fourier space can appear because the resonance conditions are a priori more difficult to satisfy (see, for example, Connaughton, Nazarenko & Pushkarev (Reference Connaughton, Nazarenko and Pushkarev2001) for capillary waves). In theory, the weaker the nonlinearities, the more important these effects are. In inertial wave turbulence, Bourouiba (Reference Bourouiba2008) has shown that discretisation effects become non-negligible when the Rossby number, 
$R_{o}$, is smaller than 
$10^{-3}$. Beyond this value, but still for a small 
$R_{o} \ll 1$, these effects are negligible because of the quasi-resonances that, with the resonances, contribute to the energy transfer.
Acknowledgements
I would like to thank and pay tribute to Professor V. Zakharov, whose contribution to wave turbulence has been profound and inspiring for me.
Funding
This research has received a grant from Institut universitaire de France.
Declaration of interests
The authors report no conflict of interest.
