Browse
Abstract
Aggregation efficiency in the upper troposphere is highly uncertain because of the lack of laboratory experiments and aircraft measurements, especially at atmospheric temperatures below −30°C. Aggregation is physically broken down into collision and sticking. In this study, theory-based parameterizations for the collision efficiency and sticking efficiency are newly implemented into a double-moment bulk cloud microphysics scheme. Satellite observations of the global ice cloud distribution are used to evaluate the aggregation efficiency modeling. Sensitivity experiments of 9-day global simulations using a high-resolution climate model show that the use of collision efficiency parameterization causes a slight increase in the cloud ice amount above the freezing level over the tropics to midlatitudes and that the use of our new sticking efficiency parameterization causes a significant increase in the cloud ice amount and a slight decrease in the snow amount particularly in the upper troposphere over the tropics. The increase/decrease in the cloud ice/snow amount in the upper troposphere over the tropics is consistent with the vertical profile of radar echoes. Moreover, the ice fraction of the cloud optical thickness is still underestimated worldwide. Finally, the cloud radiative forcing increases over the tropics to reduce the bias in the radiation budget. These results indicate that our new aggregation efficiency modeling reasonably functions even at atmospheric temperatures below −30°C; however, further improvements in the ice cloud modeling are needed.
Significance Statement
Long-standing biases in the radiative budget in climate models indicate the existence of a missing mechanism to realistically represent the ice cloud growth in the upper troposphere. This study focuses on aggregation efficiency, which has been assumed to be a tuning parameter to optimize the global radiative budget. Therefore, this study employs a theory-based parameterization to calculate the aggregation efficiency. According to the parameterization, aggregation efficiency in high clouds varies by the growth stage of the individual ice particles. As a result, small ice crystals are likely to grow more slowly, and the lifetime of cirrus clouds is prolonged to enhance cloud radiative forcing, particularly over the tropics. These results are promising for reducing the biases observed in climate models.
Abstract
Aggregation efficiency in the upper troposphere is highly uncertain because of the lack of laboratory experiments and aircraft measurements, especially at atmospheric temperatures below −30°C. Aggregation is physically broken down into collision and sticking. In this study, theory-based parameterizations for the collision efficiency and sticking efficiency are newly implemented into a double-moment bulk cloud microphysics scheme. Satellite observations of the global ice cloud distribution are used to evaluate the aggregation efficiency modeling. Sensitivity experiments of 9-day global simulations using a high-resolution climate model show that the use of collision efficiency parameterization causes a slight increase in the cloud ice amount above the freezing level over the tropics to midlatitudes and that the use of our new sticking efficiency parameterization causes a significant increase in the cloud ice amount and a slight decrease in the snow amount particularly in the upper troposphere over the tropics. The increase/decrease in the cloud ice/snow amount in the upper troposphere over the tropics is consistent with the vertical profile of radar echoes. Moreover, the ice fraction of the cloud optical thickness is still underestimated worldwide. Finally, the cloud radiative forcing increases over the tropics to reduce the bias in the radiation budget. These results indicate that our new aggregation efficiency modeling reasonably functions even at atmospheric temperatures below −30°C; however, further improvements in the ice cloud modeling are needed.
Significance Statement
Long-standing biases in the radiative budget in climate models indicate the existence of a missing mechanism to realistically represent the ice cloud growth in the upper troposphere. This study focuses on aggregation efficiency, which has been assumed to be a tuning parameter to optimize the global radiative budget. Therefore, this study employs a theory-based parameterization to calculate the aggregation efficiency. According to the parameterization, aggregation efficiency in high clouds varies by the growth stage of the individual ice particles. As a result, small ice crystals are likely to grow more slowly, and the lifetime of cirrus clouds is prolonged to enhance cloud radiative forcing, particularly over the tropics. These results are promising for reducing the biases observed in climate models.
Abstract
The rapidly increasing resolution of global atmospheric reanalysis and climate model datasets necessitates finding methods for computing convective available potential energy (CAPE) both efficiently and accurately. To this end, this article compares two common methods for computing CAPE which conserve either energy or entropy. Inaccuracies in these computations arise from both physical and numerical errors. For instance, computing CAPE with entropy conserved results in physical errors from nonequilibrium phase transitions but minimizes numerical errors because solutions are analytic at each height. In contrast, computing CAPE with energy conserved avoids these physical errors, but accumulates numerical errors that are grid-resolution-dependent because the numerical integration of a differential equation is required. Analysis of CAPE computed with large databases of soundings from the tropical Amazon and midlatitude storm environments shows that physical errors from the entropy method are typically 1%–3% as large as CAPE, which is comparable to the numerical errors from conserving energy with grid spacing of 25 and 250 m using explicit first-order and second-order integration schemes, respectively. Errors in entropy-based CAPE calculations are also insensitive to vertical grid spacing, in contrast to energy-based calculations whose error strongly scales with the grid spacing. It is shown that entropy-based methods are advantageous when intercomparing datasets with differing vertical resolution because they produce accurate and reasonably fast results that are insensitive to grid resolution, whereas a second-order energy-based method is advantageous when analyzing data with a consistent vertical resolution because of its superior computational efficiency.
Significance Statement
Convective available potential energy (CAPE) is a measure of instability in the atmosphere that helps forecasters and researchers understand when and where thunderstorms will form. The purpose of this article is to identify the most efficient and accurate methods for computing CAPE. Two methods are considered here, one that relates to the entropy (a measure of thermodynamic disorder) of an air parcel and one that relates to the energy of an air parcel. Results indicate that the entropy method is most accurate and insensitive to the resolution of the data used for the calculation (which can vary considerably), whereas the energy method uses the least computation time.
Abstract
The rapidly increasing resolution of global atmospheric reanalysis and climate model datasets necessitates finding methods for computing convective available potential energy (CAPE) both efficiently and accurately. To this end, this article compares two common methods for computing CAPE which conserve either energy or entropy. Inaccuracies in these computations arise from both physical and numerical errors. For instance, computing CAPE with entropy conserved results in physical errors from nonequilibrium phase transitions but minimizes numerical errors because solutions are analytic at each height. In contrast, computing CAPE with energy conserved avoids these physical errors, but accumulates numerical errors that are grid-resolution-dependent because the numerical integration of a differential equation is required. Analysis of CAPE computed with large databases of soundings from the tropical Amazon and midlatitude storm environments shows that physical errors from the entropy method are typically 1%–3% as large as CAPE, which is comparable to the numerical errors from conserving energy with grid spacing of 25 and 250 m using explicit first-order and second-order integration schemes, respectively. Errors in entropy-based CAPE calculations are also insensitive to vertical grid spacing, in contrast to energy-based calculations whose error strongly scales with the grid spacing. It is shown that entropy-based methods are advantageous when intercomparing datasets with differing vertical resolution because they produce accurate and reasonably fast results that are insensitive to grid resolution, whereas a second-order energy-based method is advantageous when analyzing data with a consistent vertical resolution because of its superior computational efficiency.
Significance Statement
Convective available potential energy (CAPE) is a measure of instability in the atmosphere that helps forecasters and researchers understand when and where thunderstorms will form. The purpose of this article is to identify the most efficient and accurate methods for computing CAPE. Two methods are considered here, one that relates to the entropy (a measure of thermodynamic disorder) of an air parcel and one that relates to the energy of an air parcel. Results indicate that the entropy method is most accurate and insensitive to the resolution of the data used for the calculation (which can vary considerably), whereas the energy method uses the least computation time.
Abstract
Previous observational studies have shown that the intensification rate (IR) of a tropical cyclone (TC) is often correlated with its real-time size. However, no any size parameter explicitly appears in the recent time-dependent theory of TC intensification, while the theory can still well capture the intensity evolution of simulated TCs. This study provides a detailed analysis to address how TC real-time size affects its intensification and why no size parameter explicitly appears in the theory based on the results from axisymmetric numerical simulations. The results show that the overall correlation between the TC IR and real-time size as reported in previous observational studies, in terms of both the radius of maximum wind (RMW) and the radius of 17 m s−1 wind (R17), is largely related to the correlation between the IR and intensity because the size and intensity are highly interrelated. As a result, the correlation between the TC IR and size for a given intensity is rather weak. Diagnostic analysis shows that the TC real-time size (RMW and R17) has two opposing effects on intensification. A larger TC size tends to result in a higher steady-state intensity but reduce the conversion efficiency of thermodynamic energy to inner-core kinetic energy or the degree of moist neutrality of the eyewall ascent for a given intensity. The former is favorable, while the latter is unfavorable for intensification. The two effects are implicitly included in the theory and largely offset, resulting in the weak dependence of the IR on TC size for a given intensity.
Abstract
Previous observational studies have shown that the intensification rate (IR) of a tropical cyclone (TC) is often correlated with its real-time size. However, no any size parameter explicitly appears in the recent time-dependent theory of TC intensification, while the theory can still well capture the intensity evolution of simulated TCs. This study provides a detailed analysis to address how TC real-time size affects its intensification and why no size parameter explicitly appears in the theory based on the results from axisymmetric numerical simulations. The results show that the overall correlation between the TC IR and real-time size as reported in previous observational studies, in terms of both the radius of maximum wind (RMW) and the radius of 17 m s−1 wind (R17), is largely related to the correlation between the IR and intensity because the size and intensity are highly interrelated. As a result, the correlation between the TC IR and size for a given intensity is rather weak. Diagnostic analysis shows that the TC real-time size (RMW and R17) has two opposing effects on intensification. A larger TC size tends to result in a higher steady-state intensity but reduce the conversion efficiency of thermodynamic energy to inner-core kinetic energy or the degree of moist neutrality of the eyewall ascent for a given intensity. The former is favorable, while the latter is unfavorable for intensification. The two effects are implicitly included in the theory and largely offset, resulting in the weak dependence of the IR on TC size for a given intensity.
Abstract
Flight-level airborne observations have often detected gravity waves with horizontal wavelengths
Abstract
Flight-level airborne observations have often detected gravity waves with horizontal wavelengths
Abstract
This study investigates the in situ generation of planetary waves (PWs) by zonally asymmetric gravity wave drag (GWD) in the mesosphere using a fully nonlinear general circulation model extending to the lower thermosphere. To isolate the effects of GWD, we establish a highly idealized but efficient framework that excludes stationary PWs propagating from the troposphere and in situ PWs generated by barotropic and baroclinic instabilities. The GWD is prescribed in a zonally sinusoidal form with a zonal wavenumber (ZWN) of either 1 or 2 in the lower mesosphere of the Northern Hemisphere midlatitude. Our idealized simulations clearly show that zonally asymmetric GWD generates PWs by serving as a nonconservative source Z′ of linearized disturbance quasigeostrophic potential vorticity q′. While Z′ initially amplifies PWs through enhancing q′ tendency, the subsequent zonal advection of q′ gradually balances with Z′, thereby attaining steady-state PWs. The GWD-induced PWs predominantly have the same ZWN as the applied GWD with minor contributions from higher ZWN components attributed to nonlinear processes. The amplitude of the induced PWs increases in proportion with the magnitude of the peak GWD, while it decreases in proportion to the square of ZWN. Moreover, the amplitude of PWs increases as the meridional range of GWD expands and as GWD shifts toward lower latitudes. These PWs deposit substantial positive Eliassen–Palm flux divergence (EPFD) of ∼30 m s−1 day−1 at their origin and negative EPFD of 5–10 m s−1 day−1 during propagation. In addition, the in situ PWs exhibit interhemispheric propagation following westerlies that extend into the Southern Hemisphere.
Abstract
This study investigates the in situ generation of planetary waves (PWs) by zonally asymmetric gravity wave drag (GWD) in the mesosphere using a fully nonlinear general circulation model extending to the lower thermosphere. To isolate the effects of GWD, we establish a highly idealized but efficient framework that excludes stationary PWs propagating from the troposphere and in situ PWs generated by barotropic and baroclinic instabilities. The GWD is prescribed in a zonally sinusoidal form with a zonal wavenumber (ZWN) of either 1 or 2 in the lower mesosphere of the Northern Hemisphere midlatitude. Our idealized simulations clearly show that zonally asymmetric GWD generates PWs by serving as a nonconservative source Z′ of linearized disturbance quasigeostrophic potential vorticity q′. While Z′ initially amplifies PWs through enhancing q′ tendency, the subsequent zonal advection of q′ gradually balances with Z′, thereby attaining steady-state PWs. The GWD-induced PWs predominantly have the same ZWN as the applied GWD with minor contributions from higher ZWN components attributed to nonlinear processes. The amplitude of the induced PWs increases in proportion with the magnitude of the peak GWD, while it decreases in proportion to the square of ZWN. Moreover, the amplitude of PWs increases as the meridional range of GWD expands and as GWD shifts toward lower latitudes. These PWs deposit substantial positive Eliassen–Palm flux divergence (EPFD) of ∼30 m s−1 day−1 at their origin and negative EPFD of 5–10 m s−1 day−1 during propagation. In addition, the in situ PWs exhibit interhemispheric propagation following westerlies that extend into the Southern Hemisphere.
Abstract
Effects of turbulence on ice supersaturation at cirrus heights (>8 km) remain unexplored. Small-scale mixing processes become important for high Reynolds number flows, which may develop below the buoyancy length scale (10–100 m). The current study couples a stochastic turbulent mixing model with reduced dimensionality to an entraining parcel model to investigate, in large-ensemble simulations, how supersaturation evolves due to homogeneous turbulence in the stably stratified, cloud-free upper troposphere. The rising parcel is forced by a mesoscale updraft. The perturbation of an initially homogeneous vertical distribution of supersaturation is studied after a 36-m ascent in a baseline case and several sensitivity scenarios. Turbulent mixing and associated temperature fluctuations alone lead to changes in ensemble-mean distributions with standard deviations in the range 0.001–0.006, while mean values are hardly affected. Large case-to-case variability in the supersaturation field is predicted with fluctuation amplitudes of up to 0.03, although such large values are rare. A vertical gradient of supersaturation (≈10−3 m−1) is generated for high turbulence intensities due to the development of a dry-adiabatic lapse rate. Entrainment of slightly warmer (less than 0.1 K) environmental air into the parcel decreases the mean supersaturation by less than 0.01. Supersaturation fluctuations are substantially larger after entrainment events with an additional small offset in absolute humidity (by ±3.5%) between the parcel and environmental air. The predicted perturbations of ice supersaturation are significant enough to motivate studies of turbulence–ice nucleation interactions during cirrus formation that abandon the assumption of instantaneous mixing inherent to traditional parcel models.
Significance Statement
The purpose of our study is to investigate the effects of microscale turbulence on ice supersaturation in the upper troposphere. The associated variability in temperature and moisture fields is not resolved in cloud models and cannot easily be represented in terms of large-scale flow variables. We specify the conditions in which turbulent mixing and entrainment cause substantial variations in distributions of supersaturation. These include high turbulence intensity, strong atmospheric stability, and large moisture gradients. Our results suggest that turbulence may affect the strongly supersaturation-dependent ice formation processes in high-altitude clouds, pointing to the need to investigate cirrus formation in the presence of turbulence.
Abstract
Effects of turbulence on ice supersaturation at cirrus heights (>8 km) remain unexplored. Small-scale mixing processes become important for high Reynolds number flows, which may develop below the buoyancy length scale (10–100 m). The current study couples a stochastic turbulent mixing model with reduced dimensionality to an entraining parcel model to investigate, in large-ensemble simulations, how supersaturation evolves due to homogeneous turbulence in the stably stratified, cloud-free upper troposphere. The rising parcel is forced by a mesoscale updraft. The perturbation of an initially homogeneous vertical distribution of supersaturation is studied after a 36-m ascent in a baseline case and several sensitivity scenarios. Turbulent mixing and associated temperature fluctuations alone lead to changes in ensemble-mean distributions with standard deviations in the range 0.001–0.006, while mean values are hardly affected. Large case-to-case variability in the supersaturation field is predicted with fluctuation amplitudes of up to 0.03, although such large values are rare. A vertical gradient of supersaturation (≈10−3 m−1) is generated for high turbulence intensities due to the development of a dry-adiabatic lapse rate. Entrainment of slightly warmer (less than 0.1 K) environmental air into the parcel decreases the mean supersaturation by less than 0.01. Supersaturation fluctuations are substantially larger after entrainment events with an additional small offset in absolute humidity (by ±3.5%) between the parcel and environmental air. The predicted perturbations of ice supersaturation are significant enough to motivate studies of turbulence–ice nucleation interactions during cirrus formation that abandon the assumption of instantaneous mixing inherent to traditional parcel models.
Significance Statement
The purpose of our study is to investigate the effects of microscale turbulence on ice supersaturation in the upper troposphere. The associated variability in temperature and moisture fields is not resolved in cloud models and cannot easily be represented in terms of large-scale flow variables. We specify the conditions in which turbulent mixing and entrainment cause substantial variations in distributions of supersaturation. These include high turbulence intensity, strong atmospheric stability, and large moisture gradients. Our results suggest that turbulence may affect the strongly supersaturation-dependent ice formation processes in high-altitude clouds, pointing to the need to investigate cirrus formation in the presence of turbulence.
Abstract
It has been widely reported that an increase in aerosol concentration in nonprecipitating clouds leads to a decrease in their liquid water path. Here, we examine the physical mechanisms that drive this response in both subtropical and Arctic stratocumulus clouds using large-eddy simulations and mechanism suppression tests. Three processes have been previously identified to contribute to the decrease, namely, the size dependency of evaporation, sedimentation, and radiation and all act to modulate the rate of entrainment of warm, dry air at the boundary layer top. We find that the liquid water path decrease is correlated with the increased entrainment, as expected, but that decrease is enhanced by a reduction in cloud radiative cooling. The reduced cloud radiative cooling can occur even though locally at cloud top, the radiative cooling rates are stronger and helping to enhance entrainment. We find that slower droplet sedimentation contributes to the increased entrainment and decreased liquid water in both cases. Faster evaporation caused directly by smaller, more numerous droplets decreases the liquid water path but does not necessarily increase the entrainment rate. On the other hand, stronger radiative cloud-top cooling caused directly by smaller droplets increases the entrainment as much as slower sedimentation does but does not decrease the liquid water path as much. In general, processes that either directly or indirectly increase radiative cooling at cloud top are more important in the Arctic case and processes that increase the evaporation rate are more important in the subtropical case.
Abstract
It has been widely reported that an increase in aerosol concentration in nonprecipitating clouds leads to a decrease in their liquid water path. Here, we examine the physical mechanisms that drive this response in both subtropical and Arctic stratocumulus clouds using large-eddy simulations and mechanism suppression tests. Three processes have been previously identified to contribute to the decrease, namely, the size dependency of evaporation, sedimentation, and radiation and all act to modulate the rate of entrainment of warm, dry air at the boundary layer top. We find that the liquid water path decrease is correlated with the increased entrainment, as expected, but that decrease is enhanced by a reduction in cloud radiative cooling. The reduced cloud radiative cooling can occur even though locally at cloud top, the radiative cooling rates are stronger and helping to enhance entrainment. We find that slower droplet sedimentation contributes to the increased entrainment and decreased liquid water in both cases. Faster evaporation caused directly by smaller, more numerous droplets decreases the liquid water path but does not necessarily increase the entrainment rate. On the other hand, stronger radiative cloud-top cooling caused directly by smaller droplets increases the entrainment as much as slower sedimentation does but does not decrease the liquid water path as much. In general, processes that either directly or indirectly increase radiative cooling at cloud top are more important in the Arctic case and processes that increase the evaporation rate are more important in the subtropical case.
Abstract
Gravity waves dispersing upward through the tropical stratosphere during opposing phases of the QBO are investigated using ERA5 data for 1979–2019. Log–log plots of two-sided zonal wavenumber–frequency spectra of vertical velocity, and cospectra representing the vertical flux of zonal momentum in the tropical lower stratosphere, exhibit distinctive gravity wave signatures across space and time scales ranging over two orders of magnitude. Spectra of the vertical flux of momentum are indicative of a strong dissipation of westward-propagating gravity waves during the easterly phase and vice versa. This selective “wind filtering” of the waves as they disperse upward imprints the vertical structure of the zonal flow on the resolved wave spectra, characteristic of (re)analysis and/or free-running models. The three-dimensional structures of the gravity waves are documented in composites of the vertical velocity field relative to grid-resolved tropospheric downwelling events at individual reference grid points along the equator. In the absence of a background zonal flow, the waves radiate outward and upward from their respective reference grid points in concentric rings. When a zonal flow is present, the rings are displaced downstream relative to the source and they are amplified upstream of the source and attenuated downstream of it, such that instead of rings, they assume the form of arcs. The log–log spectral representation of wind filtering of equatorial waves by the zonal flow in this paper can be used to diagnose the performance of high-resolution models designed to simulate the circulation of the tropical stratosphere.
Abstract
Gravity waves dispersing upward through the tropical stratosphere during opposing phases of the QBO are investigated using ERA5 data for 1979–2019. Log–log plots of two-sided zonal wavenumber–frequency spectra of vertical velocity, and cospectra representing the vertical flux of zonal momentum in the tropical lower stratosphere, exhibit distinctive gravity wave signatures across space and time scales ranging over two orders of magnitude. Spectra of the vertical flux of momentum are indicative of a strong dissipation of westward-propagating gravity waves during the easterly phase and vice versa. This selective “wind filtering” of the waves as they disperse upward imprints the vertical structure of the zonal flow on the resolved wave spectra, characteristic of (re)analysis and/or free-running models. The three-dimensional structures of the gravity waves are documented in composites of the vertical velocity field relative to grid-resolved tropospheric downwelling events at individual reference grid points along the equator. In the absence of a background zonal flow, the waves radiate outward and upward from their respective reference grid points in concentric rings. When a zonal flow is present, the rings are displaced downstream relative to the source and they are amplified upstream of the source and attenuated downstream of it, such that instead of rings, they assume the form of arcs. The log–log spectral representation of wind filtering of equatorial waves by the zonal flow in this paper can be used to diagnose the performance of high-resolution models designed to simulate the circulation of the tropical stratosphere.
Abstract
We investigate the mechanism for eye formation in hurricane-like vortices, using a formulation adapted from Oruba, Davidson, and Dormy. Numerical simulations are performed using an axisymmetric model of dry rotating Rayleigh–Bénard convection under the Boussinesq approximation. The fluxes of heat and momentum at the sea surface are described using the bulk aerodynamic formula. A simplified model for radiative cooling is also implemented. We find that the mechanism for eye formation introduced in Oruba et al., relying on vorticity stripping from the boundary layer, is robust in dry hurricane-like vortices. Furthermore, with these boundary conditions, the structure of the flow is closer to the flow of actual tropical cyclones. The applicability of this mechanism to the moist case however remains uncertain and deserves further study. Finally, energy budgets, obtained either by a heat engine approach or by a direct estimation of the work of buoyancy forces, are investigated. They provide estimations of the surface wind speed as a function of the controlling parameters.
Significance Statement
Tropical cyclones (also known as hurricanes or typhoons, depending on their location) are life-threatening and devastating atmospheric vortices. Their impact worsens with sea level rise and increasing coastal population. Here, we explore, using idealized models, the physics behind the formation of an eye (the quiet and typically clear region at the center of these storms). We then investigate the controlling parameters for the amplitude of the strongest winds in the model, using energy constraints, and compare them to what is observed in actual tropical cyclones.
Abstract
We investigate the mechanism for eye formation in hurricane-like vortices, using a formulation adapted from Oruba, Davidson, and Dormy. Numerical simulations are performed using an axisymmetric model of dry rotating Rayleigh–Bénard convection under the Boussinesq approximation. The fluxes of heat and momentum at the sea surface are described using the bulk aerodynamic formula. A simplified model for radiative cooling is also implemented. We find that the mechanism for eye formation introduced in Oruba et al., relying on vorticity stripping from the boundary layer, is robust in dry hurricane-like vortices. Furthermore, with these boundary conditions, the structure of the flow is closer to the flow of actual tropical cyclones. The applicability of this mechanism to the moist case however remains uncertain and deserves further study. Finally, energy budgets, obtained either by a heat engine approach or by a direct estimation of the work of buoyancy forces, are investigated. They provide estimations of the surface wind speed as a function of the controlling parameters.
Significance Statement
Tropical cyclones (also known as hurricanes or typhoons, depending on their location) are life-threatening and devastating atmospheric vortices. Their impact worsens with sea level rise and increasing coastal population. Here, we explore, using idealized models, the physics behind the formation of an eye (the quiet and typically clear region at the center of these storms). We then investigate the controlling parameters for the amplitude of the strongest winds in the model, using energy constraints, and compare them to what is observed in actual tropical cyclones.
Abstract
Previous observations and modeling studies showed that tropical cyclones (TCs) in a sheared environment develop an asymmetric boundary layer (BL). While the relationship between the BL asymmetries and environmental shear has been demonstrated, the exact cause of these BL asymmetries and the phase relationship between them are less well understood. In this study, we examine the dynamical processes leading to the asymmetric structure of the TC BL in a sheared environment using idealized, convection-permitting model simulations. Our results show that the emergence of the BL asymmetries is closely linked to the TC vortex tilt and rainband processes. Specifically, stratiform diabatic processes in the downtilt-left region result in midlevel descending inflow, which brings midtropospheric, low-θE air toward the BL and forms a surface cold pool in the downtilt-left quadrant. This descending inflow also advects high absolute angular momentum inward, redistributing the vertical vorticity and causing a storm-scale tangential wind acceleration within the downtilt-left quadrant. As the BL low-θE air advances inward, it becomes supergradient and decelerates radially, forming BL outflow in the uptilt-left quadrant. The outflow advects positive relative vorticity uptilt, forming an elliptic BL vorticity and circulation structure. As the tilted TC vortex and the accompanying rainband precess cyclonically over time, the above sequence of events and the resultant BL asymmetries also precess cyclonically, maintaining a quasi-stationary configuration relative to the vortex tilt. These results suggest that the primary organizing factor of the boundary layer asymmetries is the tilted vortex structure and not strictly the environmental shear direction.
Abstract
Previous observations and modeling studies showed that tropical cyclones (TCs) in a sheared environment develop an asymmetric boundary layer (BL). While the relationship between the BL asymmetries and environmental shear has been demonstrated, the exact cause of these BL asymmetries and the phase relationship between them are less well understood. In this study, we examine the dynamical processes leading to the asymmetric structure of the TC BL in a sheared environment using idealized, convection-permitting model simulations. Our results show that the emergence of the BL asymmetries is closely linked to the TC vortex tilt and rainband processes. Specifically, stratiform diabatic processes in the downtilt-left region result in midlevel descending inflow, which brings midtropospheric, low-θE air toward the BL and forms a surface cold pool in the downtilt-left quadrant. This descending inflow also advects high absolute angular momentum inward, redistributing the vertical vorticity and causing a storm-scale tangential wind acceleration within the downtilt-left quadrant. As the BL low-θE air advances inward, it becomes supergradient and decelerates radially, forming BL outflow in the uptilt-left quadrant. The outflow advects positive relative vorticity uptilt, forming an elliptic BL vorticity and circulation structure. As the tilted TC vortex and the accompanying rainband precess cyclonically over time, the above sequence of events and the resultant BL asymmetries also precess cyclonically, maintaining a quasi-stationary configuration relative to the vortex tilt. These results suggest that the primary organizing factor of the boundary layer asymmetries is the tilted vortex structure and not strictly the environmental shear direction.
