Take the two-layer Phillips model on the
plane as reference. Setting channel width L, we consider zonal flow without horizontal shear. For the sake of simplicity, the thickness of each layer is assumed to be D when there is no movement.
is the geostrophic stream function for each layer, where
represents the upper layer and
represents the lower layer. The nondimensional governing equations are [
8]
where
Define the nondimensional parameter
F as the ratio of channel width to deformation radius, and
is the ratio of planetary vorticity gradient to the characteristic value of relative vorticity gradient, all set to
. J is the Jacobian operator. A dissipative mechanism with a rate constant
is introduced on the right side of (1) as the deboost of potential vorticity. The basic flow of disturbance is
, The total steamfunction is
The governing equations for the perturbations are
In the absence of dissipation, if
is time-independent, then the critical value required for perturbation growth will be
[
1]. Due to the potential vorticity deboost in (6), the critical value will be slightly higher than
. We are concerned with the case where the fundamental shear changes over time and is always less than this critical value. Consider the basic form of shearing as follows:
when
the shear force at every moment is below the critical value. The baroclinic model and baroclinic model of disturbed field are used to reformulate the problem. The following definitions describe the baroclinic and baroclinic modes
It can be obtained, respectively, according to (6)
Consider the case where there is no average shear force when
. For
, the flow at every moment is much less than the critical value
[
1]. According to Flierl and Pedlosky’s research [
4], we set
. This defines the scaling speed of the shear, which is the Rossby long wave velocity
. Assume that the amplitude expansion of each flow function is
The slow time scale is
, Assuming that
and
, we define
. A geostrophic flow function can be both a function of
T and
t. Therefore, the time derivatives in (10) and (11) can be converted to
The amplitude expansion solution for the lowest-order baroclinic and baroclinic Rossby modes is
where an asterisk denotes complex conjugation, and
where the asterisk indicates complex conjugate. In order to conform the situation that there is no normal current at the channel boundary, the wave number
y needs to be an integer multiple of
. In the following, we will choose
, which is the lowest and most unstable mode. Then, the governing equation of amplitude variable has the following form
where
. Frequency
is the critical frequency that defines the parametric instability of oscillating shear resonance. In the first order of
a, these waves propagate as free Rossby waves. As we will see, parametric instability occurs at a node of frequency
, where the amplitude is still an arbitrary function of the long time variable
T. In the next order of
a, the nonlinear term applies force only in (10), and the force has no connection with
x. At this point,
and
satisfies
The solution of
can be obtained as
and
The solution can be obtained by using boundary conditions
[
8]. We corrected the oscillation frequency of the average flow to
, and no time average is included in the oscillation period. In the absence of a time average of the fundamental shear force, the shear flow does not exist in the
y direction. Then, we do not expect the thickness flux of baroclinic flow to be changed in the unstable wave, and can choose the direction of the integral flux that can produce the time-averaged flux correction. Next, we will examine the average flow correction in the presence of an average but subcritical base-state shear. When expanding to the next order of
a, we can obtain the linear problem of higher-order correction of the geostrophic flow function. The time-related questions are obtained in the form of (20), The terms on both sides of the equation oscillate at the same natural frequency. These terms will produce more and more terms in the expansion of
a; so, we need to eliminate these terms to make them invalid in other ways. Under the requirement of slow scale time
T, the amplitude equations of Rossby waves are obtained as follows:
where
The fundamental shear must oscillate at the frequency
in order to suppress terms that may resonate with the linear operator in (20). If it were not for this frequency, the second term of (26) and (27) from the interaction between oscillating shear and the lowest-order Rossby waves would not have occurred. Nonlinear systems (26) and (27) control the amplitude of the Rossby waves of (14) and (15); so, (24) and (25) are also used to describe the correction for the mean zonal flow. Where
Firstly, the nonlinear terms in (26) and (27) are temporarily ignored, and the stability of Rossby waves and their capacity for change on oscillatory shear are discussed through these linear terms. If Rossby waves are unstable, this produces the expected behavior when small amplitude perturbations begin to grow. Suppose the form of the solution is
or
. We obtain the following growth rate
Thus, oscillating shear is unstable when the amplitude exceeds the critical value below
The short-wave cut-off of the shear is consistent with the standard stability problem. The shear flow frequency of
is not required. As described below, the frequency range for which this parametric instability occurs is generally enlarged for larger
H. After a period of exponential growth, the nonlinear terms in (26) and (27) can no longer be ignored. Nevertheless, nonlinear solutions of (26) and (27) can still be solved. Their solutions are unstable and oscillate periodically with slow time scale
T. The nonlinear solution has the following form:
and
Multiplied by the complex conjugate of the positive and baroclinic amplitudes, the real and imaginary parts of (26) and (27) can be expressed as follows:
and
After substituting (36) and (37) into (26) and (27), the amplitude of the oscillation is obtained
According to (10) and (14), finite amplitude solutions can exist only when linear solutions are unstable. In addition, the predicted final state from (32) and (38) contains a frequency shift proportional to the square of the equilibrium amplitude. According to the prediction of (36)–(38), the oscillation amplitude of the solution gradually decreases, and finally, the amplitude balances to a stable value.
Figure 1 and
Figure 2 show the real and imaginary parts of the positive and baroclinic amplitudes, respectively.