HTML
--> --> -->Emanuel proposed a simple steady-state axisymmetric TC model (Emanuel, 1986, 1988, 1995) in which the air was assumed to be slantwise neutral and inviscid above the boundary layer. The axisymmetric, adiabatic and steady-state assumption is generally acceptable for the main circulation above the boundary layer in a mature TC. In Emanuel's model, the axisymmetric structure of the saturated entropy (or angular momentum) above the boundary layer was established. With a simple boundary layer model, the vertical velocity at the top of the boundary layer and the mean radial velocity within the boundary layer were also derived. However, the secondary circulation was not completely constructed in Emanuel's model, as the radial and vertical velocities above the boundary layer were not given. This hinders studies on the secondary circulation in the theoretical TC model. To obtain the secondary circulation with the SE equation, one has to retrieve the potential heat distribution from the entropy structure, which would make the calculation a bit fussy.
This study proposes a secondary circulation diagnosis method with saturated entropy structure in a steady-state axisymmetric TC under the adiabatic and inviscid assumptions, in order to facilitate the derivation of the secondary circulation in a theoretical or mature TC. The theoretical derivation of the diagnosis method is given in section 2. In section 3, a simulated idealized TC is used to verify the diagnosis method. The main results are summarized in section 4, along with some discussion.
2.1. Theory
Assume that the mature TC is steady, adiabatic and axisymmetric. According to the second law of thermodynamics: \begin{equation} \frac{{\rm d}S^{*}}{{\rm d}t}=Q , \ \ (1)\end{equation} where S* is the saturated moist entropy and Q is the diabatic heat source. For the adiabatic process, Q=0, and S* is conserved. For an axisymmetric and steady-state system, the partial derivatives with respect to time and azimuthal coordinate are zero. The total derivative can be written \begin{equation} \frac{{\rm d}}{{\rm d}t}=u\frac{\partial}{\partial r}+w\frac{\partial}{\partial z} , \ \ (2)\end{equation} where u and w are the radial and vertical velocity, respectively. So, Eq. (1) can be written for Q=0, \begin{equation} u\frac{\partial S^{*}}{\partial r}+w \frac{\partial S^{*}}{\partial z}=0 , \ \ (3)\end{equation} where r and z are the radius and height. And then, a relationship between u and w can be constructed: \begin{equation} \frac{u}{w}=\left(\frac{\partial r}{\partial z}\right)_{S^{*}} . \ \ (4)\end{equation}Physically, when the TC circulation is axisymmetric and adiabatic, S* is conserved and the air moves along the S* surface. The ratio of u to w equals the slope of the S* surface. Substituting Eq. (4) into Eq. (2), we have \begin{equation} \frac{{\rm d}}{{\rm d}t}=u\left(\frac{\partial}{\partial r}\right)_{S^{*}}=w\left(\frac{\partial}{\partial z}\right)_{S^{*}} . \ \ (5)\end{equation}
The continuity equation can be written as \begin{equation} \frac{1}{r} \frac{\partial(r u)}{\partial r}+\frac{\partial w}{\partial z}=-\frac{{\rm d}\ln \rho}{{\rm d}t} , \ \ (6)\end{equation} where ρ is the air density. Substituting Eqs. (4) and (5) into Eq. (6), \begin{equation} \frac{1}{r}\frac{\partial}{\partial r}\left[rw\left(\frac{\partial r}{\partial z}\right)_{S^{*}}\right]+\frac{\partial w}{\partial z}=-w\left(\frac{\partial \ln \rho}{\partial z}\right)_{S^{*}} . \ \ (7)\end{equation}
Dividing Eq. (7) by w, and expanding the first term, we get \begin{equation} \frac{\partial}{\partial r}\left(\frac{\partial r}{\partial z}\right)_{S^{*}}+\frac{\partial \ln (r w)}{\partial r}\left(\frac{\partial r}{\partial z}\right)_{S^{*}}+\frac{\partial \ln w}{\partial z}=-\left(\frac{\partial \ln \rho}{\partial z}\right)_{S^{*}} . \ \ (8)\end{equation}
For the left-hand side, we have \begin{eqnarray} &\dfrac{\partial}{\partial r}\left(\dfrac{\partial r}{\partial z}\right)_{S^{*}}=\left(\dfrac{\partial}{\partial z}\right)_{S^{*}} \ln \left(\dfrac{\partial r}{\partial S^{*}}\right)_{z} ,\ \ (9)&\\ &\dfrac{\partial \ln (r w)}{\partial r}\left(\dfrac{\partial r}{\partial z}\right)_{S^{*}}+\dfrac{\partial \ln w}{\partial z}=\left(\dfrac{\partial}{\partial z}\right)_{S^{*}} \ln (r w) .& \ \ (10)\end{eqnarray}
The derivations of Eq. (9) and Eq. (10) are shown in Appendix A and B, respectively. Substituting Eq. (9) and (10) into Eq. (8), we have \begin{equation} \left(\frac{\partial}{\partial z}\right)_{S^{*}} \ln \left[\rho r w\left(\frac{\partial r}{\partial S^{*}}\right)_{z}\right]=0 . \ \ (11)\end{equation}
This implies that the term inside the square brackets is a constant along an S* surface. Assuming that \begin{equation} \lambda=\rho r w\left(\frac{\partial r}{\partial S^{*}}\right)_{z} , \ \ (12)\end{equation} then Λ is a function of S*, which is constant along the S* surface, \begin{equation} w=\frac{\lambda}{\rho r}\left(\frac{\partial S^{*}}{\partial r}\right)_{z} , \ \ (13)\end{equation} and then \begin{equation} u=w\left(\frac{\partial r}{\partial z}\right)_{S^{*}}=-\frac{\lambda}{\rho r}\left(\frac{\partial S^{*}}{\partial z}\right)_{r} . \ \ (14)\end{equation}
Assuming that \begin{equation} \varLambda\left(S^{*}\right)=\int_{S_{0}^{*}}^{S^{*}} \lambda {\rm d}S^{*} , \ \ (15)\end{equation} then S0* is the value of an initial S* surface, e.g., the inner edge of the eyewall, or the axis of the TC. $\varLambda$ is also constant along the S* surface. Then, Eqs. (13) and (14) can be written as: \begin{eqnarray} w&=&\frac{1}{\rho r} \frac{{\rm d}\Lambda}{{\rm d}S^{*}}\left(\frac{\partial S^{*}}{\partial r}\right)_{z}=\frac{1}{\rho r} \frac{\partial \Lambda}{\partial r} , \ \ (16)\\ u&=&-\frac{1}{\rho r} \frac{{\rm d}\Lambda}{{\rm d}S^{*}}\left(\frac{\partial S^{*}}{\partial z}\right)_{r}=-\frac{1}{\rho r} \frac{\partial \Lambda}{\partial z} . \ \ (17)\end{eqnarray}
Equations (16) and (17) have the form of a streamfunction. Thus, $\varLambda$ is a streamfunction in a steady-state, adiabatic and axisymmetric TC, and the S* surfaces coincide with the streamlines. All the equations above can be easily changed to isobaric coordinates by simply replacing z and w with p and ω, which are the pressure and the vertical velocity in isobaric coordinates respectively.
With the steady-state, adiabatic and axisymmetric assumption, the radial and vertical components of the secondary circulation in TC can be expressed with a function of saturated moist entropy in simple forms. This indicates that the secondary circulation and the structure of the saturated moist entropy depend on each other. With the structure of saturated moist entropy, it is possible to diagnose the secondary circulation.
2
2.2. Specific feature of the theory
Considering the annular air column between two isentropic surfaces S1* and S2*, if the air density and vertical velocity, or the product of them (ρ w), are horizontally homogeneous, the horizontal integration of Eq. (13) from S1* to S2* is \begin{equation} \rho w\left(r_{2}^{2}-r_{1}^{2}\right)=2 \int_{S_{1}^{2}}^{S_{2}^{*}} \lambda {\rm d}S^{*}=2\left[\Lambda\left(S_{2}^{*}\right)-\Lambda\left(S_{1}^{*}\right)\right]={\rm const} , \ \ (18)\end{equation} where r1 and r2 are the radii of the isentropic surfaces S1* and S2* respectively. For a small piece of the annular air column moving between the two S* surfaces, its mass is constant when the entrainment effect is ignored: \begin{equation} m=\rho \pi\left(r_{2}^{2}-r_{1}^{2}\right)H={\rm const} , \ \ (19)\end{equation} where H is the vertical thickness of the air piece. Multiplying Eq. (18) by π/m, we get \begin{equation} \frac{w}{H}=\frac{2\pi\left[\varLambda\left(S_{2}^{*}\right)-\varLambda\left(S_{1}^{*}\right)\right]}{m}={\rm const} . \ \ (20)\end{equation}Equations (18) and (20) indicate the relationship between the variation of the vertical velocity and the volume change in vertical and horizontal components. The variation of ρ represents the volume change effect of the air piece during its vertical movement, and the variation of (r22-r12) and H represent the change of the horizontal area and vertical thickness of the air piece, which are the horizontal and vertical components of the volume change respectively. For an upward moving air piece, the volume expands as a result of pressure decline and latent heating. When the vertical velocity increases, the vertical thickness of the air piece increases, and then the horizontal area decreases or its expansion is partly cancelled. The physics is that the vertical acceleration (deceleration) leads to a vertical stretching (compression), which would affect the horizontal area when accompanying a volume change. On the other hand, when the variation of the horizontal area is constrained under some kind of dynamical balance (e.g., the slope of S* surfaces is constrained by the thermal wind balance in Emanuel's theoretical TC model), then the variation of vertical velocity is controlled by the volume change, which is determined by the gas state equation. The relationship between vertical velocity and horizontal area comes from the constraint of the continuity equation.
2
2.3. Physical meaning of Λ and $varLambda$
Equation (12) can be written as \begin{equation} \lambda=\frac{1}{2\pi}\rho w\left[\frac{\partial(\pi r^{2})}{\partial S^{*}}\right]_{z} . \ \ (12)\end{equation}In cylindrical coordinates, $[\partial(\pi r^2)/\partial S^*]_z$ means the area between two S* surfaces with difference of unit S*. So, $\rho w[\partial(\pi r^2)/\partial S^*]_z$ indicates the mass that passes through the horizontal area between the two S* surfaces with difference of unit S* in unit time, and the physical meaning of Λ is the passing mass in units of radians.
According to Eq. (15), $\varLambda$ is the integral of Λ from the S0* surface to the S* surface. So, the physical meaning of $\varLambda$ is that the mass passes through the horizontal area between two surfaces, S* and S0*, in units of time and units of radians.
2
2.4. Diagnosis method
In order to diagnose the secondary circulation of the steady axisymmetric TC, the values of Λ have to be computed for each S* surface using a reference level (e.g., the top of the boundary layer) first: \begin{equation} \lambda=\rho_{0} r_{0} w_{0}\left(\frac{\partial r}{\partial S^{*}}\right)_{z_{0}} . \ \ (22)\end{equation}The value of w on the reference level is needed here. In a theoretical TC model, this can be closed by inducing a boundary layer model, which can provide the radial profile of the vertical velocity on the top of the boundary layer (Emanuel, 1986; Kepert, 2001; Smith et al., 2018). Since ${(\partial r/ \partial S^{*})}_{z0}$ is the reciprocal of the radial gradient of S*, the level with low radial gradient of S* (i.e., the S* surfaces are almost horizontal) would induce large error to the computation of Λ. Besides, the low vertical velocity would make the value of Λ very sensitive to asymmetry. Thus, the reference level should avoid the area with low radial gradient of S* and low vertical velocity on the main channel of the secondary circulation.
When the value of Λ is computed on the reference level, the value of S* is the only clue to finding the value of Λ for points on the other levels. Thus, the distribution of S* on the reference level should be monotonic. The portion with repeated values of S* on the reference level has to be skipped. The surfaces with the same S* value but not originated from the reference level also have to be excluded because the corresponding values of Λ may be different. Usually, the main channel of the secondary circulation is left after the exclusion. With the distribution of S* in the r-z plane, the radial and vertical gradient of S* and the value of Λ for each available point can be obtained, then the radial and vertical components of the secondary circulation can be calculated with Eq. (16) and (17).
The intensity of the idealized TC increases to 95 m s-1 and the TC becomes highly axisymmetric after 96 h of integration. The azimuthal mean structure of the equivalent potential temperature (θe) in the r-z plane is shown in Fig. 1 (black solid contours), along with the vertical velocity (color fill) and the 95% relative humidity (RH, red dashed contours). The less saturated portion with RH <95% in the updraft channel implies that dry air may intrude into the updraft and evaporation may happen, or there may be downward motion when computing the azimuthal mean in the updraft channel. The former situation breaks the adiabatic assumption, while the latter breaks the axisymmetric assumption. Therefore, the less saturated portion will be excluded when diagnosing the secondary circulation. This exclusion also helps to exclude the surfaces with the same θe value but not originated from the reference level. The level of 2 km is chosen as the reference level. The values of Λ within the monotonic portion are calculated. The diagnosed radial and vertical velocities along the main updraft are shown in Fig. 2 and Fig. 3.
![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2019/8/PIC/0256-1530-36-8-804/190054fig1_mini.jpg)
![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2019/8/PIC/0256-1530-36-8-804/190054fig2_mini.jpg)
![](http://www.iapjournals.ac.cn/fileDQKXJZ/journal/article/dqkxjz/2019/8/PIC/0256-1530-36-8-804/190054fig3_mini.jpg)
In the vertical component (w), there are three peaks in the simulation, and the diagnosed result captures this feature well (Fig. 2a). The relative difference between the diagnosis and the simulation is defined as dwrela=(wdiag-w)/w for quantitative validation. The relative differences over the peaks are between -0.2 to 0.2, indicating good estimations (Fig. 2b). There are some large biases at the edge of the updraft. These bias may be caused either by the θe surfaces with different origin, which is not excluded completely with RH <95%, or by the effect of asymmetry on the radial gradient of θe. Note that there is a drastic variation in the radial gradient of θe on the inner edge of the eyewall above 4 km (Fig. 1), which makes the radial gradient of θe highly sensitive to the asymmetry of the eyewall radius. This area corresponds well with the large-error area at the upper inner side of the eyewall (Fig. 2b).
There are also three peaks for the simulated positive radial velocity (u) area (Fig. 3a). However, the diagnosis covers only two of them: the mid-level one and the low-level one. Errors are small in the two peak areas, but become large below 6 km on the inner edge of eyewall (Fig. 3b). Unlike the radial gradient of θe, the drastic variation in the vertical gradient of θe happens in the vicinity where the slopes of the θe surfaces turn from vertical to slantwise or from slantwise to vertical, e.g., the inner edge of the eyewall between 1 and 6 km (Fig. 1). Asymmetry in this area could easily induce error to the vertical gradient of θe, and then to the diagnosed result. The diagnosed values of u within the boundary layer are also shown. Note that the adiabatic assumption is not fully satisfied since there is a diabatic heat source contributed by sensible heating. However, the inflow is reproduced well by the diagnosis, except that the values are underestimated in the main channel and overestimated on the top and near the surface right below the eyewall.
Compared with the SE equation, this diagnosis method is simpler. The prescribed fields needed in the method are density, entropy and vertical velocity on the reference level, which are easy to obtain from simulations and a theoretical TC model (e.g., Emanuel, 1986). Although the prescribed fields on the reference level are limited in the monotonic portion, they cover the main secondary circulation channel in a simulated TC. The limitation is not a problem in a theoretical TC model because the entropy is monotonic at any level above the boundary layer (see Emanuel, 1986). Unlike the SE equation, this diagnosis method does not need preset proper boundary conditions or to satisfy an elliptical discriminant condition (Fudeyasu and Wang, 2011). The thermal wind balance is not necessary for the diagnosis method and the application is valid in the boundary layer. Thus, the method proposed in this study facilitates the diagnosis of the secondary circulation in a simulated mature TC, especially in a theoretical TC model.
This diagnosis method comes from constraint of the continuity equation with adiabatic, axisymmetric, steady-state assumptions. These assumptions are shared in Emanuel's theory (Emanuel, 1986). The dynamics constraining the entropy structure above the boundary layer in Emanuel's theory are the gradient wind balance and the hydrostatic balance, which came from the first order (or zero order in Willoughby, 1979) approximation of the horizontal and vertical momentum equations respectively. However, the constraint of the continuity equation is not introduced in that theory. Emanuel brought in a slantwise neutral assumption to exclude the ambient convective available potential energy, which implied that the air is neutrally buoyant when ascending along surfaces of constant angular momentum and there is no upward acceleration. However, the vertical velocity should decrease when air parcels approach the top layer of the TC (Camp and Montgomery, 2001). The diagnosis method proposed in this study provides a chance to determine how the dynamic balances constrain the secondary circulation in a theoretical TC model.
3
APPENDIX A
3Derivation of Eq. (9)
$\partial/\partial r$ and $\partial/\partial z$ can be written as $[\partial/(\partial r)]_z$ and $(\partial/\partial z)_r$. The left-hand side of Eq. (9) can transform as below: \begin{equation} \left(\frac{\partial}{\partial r}\right)_{z}\left(\frac{\partial r}{\partial z}\right)_{S^{*}} =\left(\frac{\partial S^{*}}{\partial r}\right)_{z}\left(\frac{\partial}{\partial S^{*}}\right)_{z}\left(\frac{\partial r}{\partial z}\right)_{S^{*}} . \ \ (A1)\end{equation}Since \begin{equation} \left(\frac{\partial}{\partial S^{*}}\right)_{z}\left(\frac{\partial r}{\partial z}\right)_{S^{*}}=\left(\frac{\partial}{\partial z}\right)_{S^{*}}\left(\frac{\partial r}{\partial S^{*}}\right)_{z} , \ \ (A2)\end{equation} and \begin{equation} \left(\frac{\partial S^{*}}{\partial r}\right)_{z}=\frac{1}{\left(\frac{\partial r}{\partial S^{*}}\right)_{z}} , \ \ (A3)\end{equation} Eq. (A1) can be written as \begin{equation} \left(\frac{\partial}{\partial r}\right)_{z}\left(\frac{\partial r}{\partial z}\right)_{S^{*}}=\left(\frac{\partial}{\partial z}\right)_{S^{*}} \ln \left(\frac{\partial r}{\partial S^{*}}\right)_{z} . \ \ (A4)\end{equation}
3
APPENDIX B
3Derivation of Eq. (10)
Along the r surface, $\ln r$ is constant. So, there is \begin{equation} \frac{\partial \ln r}{\partial z}=0 . \ \ (B1)\end{equation} Adding $\partial \ln w/\partial z$ on both sides, we get \begin{equation} \frac{\partial \ln r}{\partial z}+\frac{\partial \ln w}{\partial z}=\frac{\partial \ln w}{\partial z} . \ \ (B2)\end{equation} Combining the two terms on the left-hand side gives \begin{equation} \frac{\partial \ln (r w)}{\partial z}=\frac{\partial \ln w}{\partial z} . \ \ (B3)\end{equation}For a continuously differentiable variable f, there is \begin{equation} \frac{\partial f}{\partial r}\left(\frac{\partial r}{\partial z}\right)_{S^{*}}+\frac{\partial f}{\partial z}=\left(\frac{\partial f}{\partial z}\right)_{S^{*}} . \ \ (B4)\end{equation} When $f=\ln(rw)$, Eq. (B4) can be written as \begin{equation} \frac{\partial \ln (r w)}{\partial r}\left(\frac{\partial r}{\partial z}\right)_{S^{*}}+\frac{\partial \ln (r w)}{\partial z}=\left(\frac{\partial}{\partial z}\right)_{S^{*}} \ln (r w) . \ \ (B5)\end{equation} Substituting Eq. (B3) into Eq. (B5), \begin{equation} \frac{\partial \ln (r w)}{\partial r}\left(\frac{\partial r}{\partial z}\right)_{S^{*}}+\frac{\partial \ln w}{\partial z}=\left(\frac{\partial}{\partial z}\right)_{S^{*}} \ln (r w) . \ \ (B6)\end{equation}