Abstract

A conservative scheme for a compressible nonhydrostatic model including moist processes is formulated and is tested for experiments involving a squall line. The scheme is based on the flux form equations of total density, momentum, total energy, and the densities of water substance. The time-splitting scheme is used for the temporal scheme. In the small time step integration, the horizontal components of momentum are explicitly integrated, while the vertical components of momentum, density, and total energy are implicitly integrated. In particular, the flux form equation for the total energy is used to guarantee the conservation of the total energy. The internal energy is obtained by subtracting the kinetic energy and the potential energy from the total energy. This method is advantageous for the energy budget analysis. Only the warm rain cloud process is included for cloud physics. Using the squall-line experiments, the water budget and the energy budget are diagnosed and it is confirmed that the conservation of water and total energy is well satisfied.

As a quantitative improvement, more accurate formulas are used for the thermodynamics of the moist atmosphere by taking account of the effects of specific heats of the water substance and the temperature dependency of latent heat. These effects are generally neglected in most numerical models. If accurate moist thermodynamics are used, the total rain is reduced by more than 10% in comparison to the case when the simplified thermodynamics are used. The transportations of physical quantities due to rain are also appropriately introduced using a higher-order advection scheme. In the flux-form formulation, it is found that the change in energy due to the transportation of rain cannot be neglected in general, while that in momentum could be negligible.

Abstract

A new dynamical scheme with the conservative forms of the equations of density, momentum, and internal energy is proposed for the nonhydrostatic models. With this scheme, the conservations of the mass and the total energy are satisfied within round-off errors. In particular, methods for the integration of energy are discussed in detail, and three of the approaches are compared; one is in the form of the pressure equation, the second is the integration of internal energy with corrections on the transformation terms, and the third is use of the sum of internal energy and kinetic energy as a prognostic variable. This scheme is incorporated into a nonhydrostatic model with the horizontally explicit and vertically implicit time integration scheme for sound waves, and various numerical experiments for the dry atmosphere are performed. The numerical results show that the conservative properties are well satisfied.

Abstract

A newly developed global nonhydrostatic model is used for life cycle experiments (LCEs) of baroclinic waves, and the resolution dependency of frontal structures is examined. LCEs are integrated for 12 days with horizontal grid intervals ranging from 223 to 3.5 km in a global domain. In general, fronts become sharper and corresponding vertical flow strengthens as horizontal resolution increases. However, if the ratio of vertical and horizontal grid intervals is sufficiently small compared to the frontal slope s, the overall frontal structure remains unchanged. In contrast, when the ratio of horizontal and vertical grid intervals exceeds 2s − 4s, spurious gravity waves are generated at the cold front. A linear model for mountain waves quantitatively explains the mechanism of the spurious waves. The distribution of the basic wind is the major factor that determines wave amplitude and propagation. The spurious waves propagate up to a critical level at which the basic wind speed normal to the front is equal to the propagation speed of the front. Results from the linear model suggest that an effective way to eliminate spurious waves is to choose a stretched grid with a smaller vertical grid interval in lower layers where strong horizontal winds exist.

Abstract

Super Cyclone Pam (2015) formed in the central tropical Pacific under conditions that included El Niño Modoki and the passage of a convectively enhanced phase of the Madden–Julian oscillation (MJO) in the western Pacific. This study examines the influence that sea surface temperature anomalies (SSTAs) have on the MJO and low-frequency large-scale circulation, and establishes how they modulated the genesis of Pam. Two series of numerical experiments were conducted by using a nonhydrostatic global atmospheric model with observed (OBSSST) and climatological (CLMSST) SSTs. The results suggested that low-frequency westerly winds at 850 hPa (U850) were intensified in the central tropical Pacific due to the observed SSTA. The amplitude of the MJO simulated in OBSSST was larger than in CLMSST. In addition, the experiments initialized 26 February–2 March exhibited that the phase of the MJO in OBSSST was ahead of that in CLMSST, and that the genesis location in OBSSST was ~10° to the east of that in CLMSST. An analysis of large-scale fields indicated that a positive U850 maintained by SSTAs and intensification of U850 by the MJO modified distribution of large-scale cyclonic vorticity and precipitable water. These changes in large-scale fields modified the location and timing of intensification of the disturbance that become Pam and resulted in Pam’s genesis location being 10° farther east with slight impact on its genesis probability. Additional experiments showed that SSTAs in the central tropical Pacific were the dominant cause of modifications to large-scale fields, the MJO, and Pam’s genesis location.

Abstract

In this study, the single-moment 6-class bulk cloud microphysics scheme used in the operational numerical weather prediction system at the Japan Meteorological Agency was improved using the observations of the Global Precipitation Measurement (GPM) core satellite as reference values. The original cloud microphysics scheme has the following biases: underestimation of cloud ice compared to the brightness temperature of the GPM Microwave Imager (GMI) and underestimation of the lower-troposphere rain compared to the reflectivity of GPM Dual-frequency Precipitation Radar (DPR). Furthermore, validation of the dual-frequency rate of reflectivity revealed that the dominant particles in the solid phase of simulation were graupel and deviated from DPR observation. The causes of these issues were investigated using a single-column kinematic model. The underestimation of cloud ice was caused by a high ice-to-snow conversion rate, and the underestimation of precipitation in the lower layers was caused by an excessive number of small-diameter rain particles. The parameterization of microphysics scheme is improved to eliminate the biases in the single-column model. In the forecast obtained using the improved scheme, the underestimation of cloud ice and rain is reduced. Consequently, the prediction errors of hydrometeors were reduced against the GPM satellite observations, and the atmospheric profiles and precipitation forecasts were improved.

Abstract

The Nonhydrostatic Icosahedral Atmospheric Model (NICAM), a global cloud-system-resolving model, successfully simulated the life cycle of Tropical Storm Isobel that formed over the Timor Sea in the austral summer of 2006. The multiscale interactions in the life cycle of the simulated storm were analyzed in this study. The large-scale aspects that affected Isobel’s life cycle are documented in this paper and the corresponding mesoscale processes are documented in a companion paper.

The life cycle of Isobel was largely controlled by a Madden–Julian oscillation (MJO) event and the associated westerly wind burst (WWB). The MJO was found to have both positive and negative effects on the tropical cyclone intensity depending on the location of the storm relative to the WWB center associated with the MJO. The large-scale low-level convergence and high convective available potential energy (CAPE) downwind of the WWB center provided a favorable region to the cyclogenesis and intensification, whereas the strong large-scale stretching deformation field upwind of the WWB center may weaken the storm by exciting wavenumber-2 asymmetries in the eyewall and leading to the eyewall breakdown.

Five stages are identified for the life cycle of the simulated Isobel: the initial eddy, intensifying, temporary weakening, reintensifying, and decaying stages. The initial eddy stage was featured by small-scale/mesoscale convective cyclonic vortices developed in the zonally elongated rainband organized in the preconditioned environment characterized by the WWB over the Java Sea associated with the onset of an MJO event over the East Indian Ocean. As the MJO propagated eastward and the cyclonic eddies moved southward into an environment with weak vertical shear and strong low-level cyclonic vorticity, a typical tropical cyclone structure developed over the Java Sea, namely the genesis of Isobel. Isobel experienced an eyewall breakdown and a temporary weakening when it was located upwind of the WWB center as the MJO propagated southeastward and reintensified as its eyewall reformed as a result of the axisymmetrization of an inward spiraling outer rainband that originally formed downwind of the WWB center. Finally Isobel decayed as it approached the northwest coast of Australia.

Abstract

The life cycle of Tropical Storm Isobel was simulated reasonably well in the Nonhydrostatic Icosahedral Atmospheric Model (NICAM), a global cloud-system-resolving model. The evolution of the large-scale circulation and the storm-scale structure change was discussed in . Both the mesoscale and system-scale processes in the life cycle of the simulated Isobel are documented in this paper. In the preconditioned favorable environment over the Java Sea, mesoscale convective vortices (model MCVs) developed in the mesoscale convective systems (MCSs) and convective towers with cyclonic potential vorticity (PV) anomalies throughout the troposphere [model vortical hot towers (VHTs)] appeared in the model MCVs. Multiple model VHTs strengthened cyclonic PV in the interior of the model MCV and led to the formation of an upright monolithic PV core at the center of the concentric MCV (primary vortex enhancement). As the monolithic PV core with a warm core developed near the circulation center, the intensification and the increase in horizontal size of the cyclonic PV were enhanced through the system-scale intensification (SSI) process (the secondary vortex enhancement), leading to the genesis of Isobel over the Timor Sea. The SSI process can be well explained by the balanced dynamics.

After its genesis, the subsequent evolution of the simulated Isobel was controlled by both the external influence and the internal dynamics. Under the unfavorable environmental conditions, the development of asymmetric structure reduced the axisymmetric diabatic heating in the inner core and the SSI process became ineffective and the storm weakened. Later on, as the eyewall reformed as a result of the axisymmetrization of an inward-propagating outer spiral rainband, the SSI process became effective again, leading to the reintensification of Isobel. Therefore, the large-scale environmental flow provided the precondition for the genesis of Isobel and the triggering mechanism for subsequent storm-scale structure change as discussed in . The system-scale and mesoscale processes, such as the evolution of MCVs and merging VHTs, were responsible for the genesis, while the eyewall processes were critical to the storm intensity change through the SSI process.

Abstract

The degree of gradient wind balance was investigated in a number of tropical cyclones (TCs) simulated under realistic environments. The results of global-scale numerical simulations without cumulus parameterization were used, with a horizontal mesh size of 7 km. On average, azimuthally averaged maximum tangential velocities at 850 (925) hPa in the simulated TCs were 0.72% (1.95%) faster than gradient wind–balanced tangential velocity (GWV) during quasi-steady periods. Of the simulated TCs, 75% satisfied the gradient wind balance at the radius of maximum wind speed (RMW) at 850 and at 925 hPa to within about 4.0%. These results were qualitatively similar to those obtained during the intensification phase. In contrast, averages of the maximum and minimum deviations from the GWV, in all the azimuths at the RMW, achieved up to 40% of the maximum tangential velocity. Azimuthally averaged tangential velocities exceeded the GWV (i.e., supergradient) inside the RMW in the lower troposphere, whereas the velocities were close to or slightly slower than GWV (i.e., subgradient) in the other regions. The tangential velocities at 925 hPa were faster (slower) in the right-hand (left hand) side of the TC motion. When the tangential velocities at the RMW were supergradient, the primary circulation tended to decay rapidly in the vertical direction and slowly in the radial direction, and the eyewall updraft and the RMW were at larger radii. Statistical analyses revealed that the TC with supergradient wind at the RMW at 850 hPa was characterized by stronger intensity, larger RMW, more axisymmetric structure, and an intensity stronger than potential intensity.

Abstract

A nonhydrostatic stretched-grid (SG) model is used to analyze the large-scale errors generated by stretching horizontal grids and their influence on a region of interest. Simulations by a fully compressible, nonhydrostatic global atmospheric model, the Nonhydrostatic Icosahedral Atmospheric Model (NICAM), and its SG regional model, stretched-NICAM, were performed for the months of March, April, and May of 2011 using various resolutions and stretching factors. A comparison of week-long accumulative precipitation amounts between the Tropical Rainfall Measuring Mission (TRMM) satellite data and the quasi-uniform and SG simulations showed that a stretched run better represents storms and associated precipitation because the errors generated in the outer regions with coarser grid spacing do not significantly affect the inner domain centered at the focal point. For season-long simulations, in one particular set of stretched runs with the focal point located in the eastern United States, the artificial suppression of baroclinic development of midlatitude eddies in the Southern Hemisphere weakened the eddy-driven polar-front jet (PFJ), which yielded a cold bias at mid- to high latitudes. However, in the Northern Hemisphere, in contrast, the aforementioned changes are less apparent. Therefore, for the SG runs, the mean temperature was maintained at the region of interest, and an increased amount of moderate to heavy precipitation, which is also frequently found in the TRMM data, was observed; thus, the benefits of increased resolution were realized. However, careful attention must be given when applying the SG model because a regional climate response to the change in the large-scale circulations may not be fully accounted for.

Abstract

A nonhydrostatic, regional climate limited-area model (LAM) was used to analyze lateral boundary condition (LBC) errors and their influence on the uncertainties of regional models. Simulations using the fully compressible nonhydrostatic LAM (D-NICAM) were compared against the corresponding global quasi-uniform-grid Nonhydrostatic Icosahedral Atmospheric Model (NICAM) and a stretched-grid counterpart (S-NICAM). By this approach of sharing the same dynamical core and physical schemes, possible causes of model bias and LBC errors are isolated. The simulations were performed for a 395-day period from March 2011 through March 2012 with horizontal grid intervals of 14, 28, and 56 km in the region of interest. The resulting temporal mean statistics of the temperatures at 500 hPa were generally well correlated between the global and regional simulations, indicating that LBC errors had a minor impact on the large-scale flows. However, the time-varying statistics of the surface precipitation showed that the LBC errors lead to the unpredictability of convective precipitation, which affected the mean statistics of the precipitation distributions but induced only minor influences on the large-scale systems. Specifically, extratropical cyclones and orographic precipitation are not severely affected. It was concluded that the errors of the precipitation distribution are not due to the difference of the model configurations but rather to the uncertainty of the system itself. This study suggests that applications of ensemble runs, internal nudging, or simulations with longer time scales are needed to obtain more statistically significant results of the precipitation distribution in regional climate models.