Editorial Type:
Article
Article Type:
Research Article

Numerical Simulation of the Diurnal Evolution of Tropical Island Convection over the Maritime Continent

Kazuo Saito Meteorological Research Institute, Tsukuba, Japan

Search for other papers by Kazuo Saito in
Current site
Google Scholar
PubMed Close
,
Tom Keenan Bureau of Meteorology Research Centre, Melbourne, Australia

Search for other papers by Tom Keenan in
Current site
Google Scholar
PubMed Close
,
Greg Holland Bureau of Meteorology Research Centre, Melbourne, Australia

Search for other papers by Greg Holland in
Current site
Google Scholar
PubMed Close
, and
Kamal Puri Bureau of Meteorology Research Centre, Melbourne, Australia

Search for other papers by Kamal Puri in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Mar 2001
DOI:
https://doi.org/10.1175/1520-0493(2001)129<0378:NSOTDE>2.0.CO;2
Page(s):
378–400
Restricted access

Abstract

Numerical simulations of the diurnal evolution of tropical island convection observed during the Maritime Continent Thunderstorm Experiment (MCTEX) are performed using the Meteorological Research Institute nonhydrostatic model (MRI NHM). The MRI NHM is double-nested within a form of the Australian Bureau of Meteorology Research Centre’s Limited-Area Assimilation and Prediction System specially operated for the MCTEX period.

Excellent agreement is found between the simulation and observed evolution of the convective clouds over the Tiwi Islands on 27 November 1995. A transition from horizontal convection occurring during the morning to vertical convection in afternoon is evident.

In the morning, the sea breeze appears along the coastlines, with a clear contrast evident in structure between the windward and leeward sides. At the windward coast, the sea breeze intrudes inland more rapidly, where the larger surface heat flux modifies the lowest air mass and makes the sea breeze front (SBF) indistinct. On the other hand, at the leeward coast, the upward motion at the head of the SBF is larger and deeper. Shallow convective clouds therefore have a preference for alignment along the leeward SBF. Over the interior of the islands ahead of SBFs, shallow convective clouds corresponding to the Rayleigh–Benard convection occur at corners of open polygonal shaped cells and seem randomly distributed. Within the SBFs, organization of convection characteristic of horizontal convective rolls (HCRs) is evident. These HCRs are preferred at the windward coast and occur within cloud-free regions. Clouds associated with the SBFs appear to develop preferentially at the cross points of the SBFs and HCRs where the surface convergence is enhanced.

Following further inland propagation of SBFs, weak precipitation starts and the Rayleigh–Benard convection is disturbed by resulting outflows. At the merging stage, the clouds organize at the leeward central part of the islands in the form of an east–west line. In this convergence zone between the two SBFs, explosive growth of convection occurs and cloud top reaches the tropopause. In the case simulated here, the associated downdrafts are not strong compared with the upward motion due to a lack of the midlevel dry air necessary to enhance evaporative cooling.

The inclusion of ice phase physics in the simulation produces little qualitative difference in storm development and the associated surface rainfall distribution, but yields stronger updrafts and higher cloud-top heights. The vertical profile of the apparent heat source (Q1) in the ice phase experiment shows double peaks corresponding to the condensation and freezing levels.

Sensitivity experiments show that the orographic undulations as well as the horizontal scale of the island are important factors determining the timing of cloud merger and convective intensity. Without hills, the transition to the explosive growth in the merger stage is delayed. This results in weaker rainfall, even if the hills are relatively flat. A smaller island produced weaker convection, which means that the total rain produced by each island is not proportional to island area. These results suggest that the intensity of tropical island convection is determined not only by the convective stability of the environmental atmosphere but is influenced significantly by the island-scale circulations, that is, horizontal convection in the morning that ultimately forces the deep convection during the afternoon.

Corresponding author address: Dr. Kazuo Saito, Forecast Research Department, Meteorological Research Institute, 1-1 Nagamine, Tsukuba, Ibaraki 305-0052, Japan.

Abstract

Numerical simulations of the diurnal evolution of tropical island convection observed during the Maritime Continent Thunderstorm Experiment (MCTEX) are performed using the Meteorological Research Institute nonhydrostatic model (MRI NHM). The MRI NHM is double-nested within a form of the Australian Bureau of Meteorology Research Centre’s Limited-Area Assimilation and Prediction System specially operated for the MCTEX period.

Excellent agreement is found between the simulation and observed evolution of the convective clouds over the Tiwi Islands on 27 November 1995. A transition from horizontal convection occurring during the morning to vertical convection in afternoon is evident.

In the morning, the sea breeze appears along the coastlines, with a clear contrast evident in structure between the windward and leeward sides. At the windward coast, the sea breeze intrudes inland more rapidly, where the larger surface heat flux modifies the lowest air mass and makes the sea breeze front (SBF) indistinct. On the other hand, at the leeward coast, the upward motion at the head of the SBF is larger and deeper. Shallow convective clouds therefore have a preference for alignment along the leeward SBF. Over the interior of the islands ahead of SBFs, shallow convective clouds corresponding to the Rayleigh–Benard convection occur at corners of open polygonal shaped cells and seem randomly distributed. Within the SBFs, organization of convection characteristic of horizontal convective rolls (HCRs) is evident. These HCRs are preferred at the windward coast and occur within cloud-free regions. Clouds associated with the SBFs appear to develop preferentially at the cross points of the SBFs and HCRs where the surface convergence is enhanced.

Following further inland propagation of SBFs, weak precipitation starts and the Rayleigh–Benard convection is disturbed by resulting outflows. At the merging stage, the clouds organize at the leeward central part of the islands in the form of an east–west line. In this convergence zone between the two SBFs, explosive growth of convection occurs and cloud top reaches the tropopause. In the case simulated here, the associated downdrafts are not strong compared with the upward motion due to a lack of the midlevel dry air necessary to enhance evaporative cooling.

The inclusion of ice phase physics in the simulation produces little qualitative difference in storm development and the associated surface rainfall distribution, but yields stronger updrafts and higher cloud-top heights. The vertical profile of the apparent heat source (Q1) in the ice phase experiment shows double peaks corresponding to the condensation and freezing levels.

Sensitivity experiments show that the orographic undulations as well as the horizontal scale of the island are important factors determining the timing of cloud merger and convective intensity. Without hills, the transition to the explosive growth in the merger stage is delayed. This results in weaker rainfall, even if the hills are relatively flat. A smaller island produced weaker convection, which means that the total rain produced by each island is not proportional to island area. These results suggest that the intensity of tropical island convection is determined not only by the convective stability of the environmental atmosphere but is influenced significantly by the island-scale circulations, that is, horizontal convection in the morning that ultimately forces the deep convection during the afternoon.

Corresponding author address: Dr. Kazuo Saito, Forecast Research Department, Meteorological Research Institute, 1-1 Nagamine, Tsukuba, Ibaraki 305-0052, Japan.

Save
Cover Monthly Weather Review
Monthly Weather Review

  • Asai, T., 1970: Stability of plane parallel flow with variable vertical shear and unstable stratification. J. Meteor. Soc. Japan,48, 129–139.

  • ——, 1983: Science of Atmospheric Convection (in Japanese). Tokyo-do, 220 pp.

  • Carbone, R. E., T. D. Keenan, J. Hacker, and J. W. Wilson, 2000: Tropical island convection in the absence of significant topography. Part I: Life cycle of diurnally forced convection. Mon. Wea. Rev.,128, 3459–3480.

  • Deardorff, J. W., 1980: Stratocumulus-capped mixed layers derived from a three-dimensional model. Bound.-Layer Meteor.,18, 495–527.

  • Droegemeier, K. K., and R. B. Wilhelmson, 1985: Three-dimensional numerical modeling of convection produced by interacting thunderstorm outflows. Part I: Control simulation and low-level moisture. J. Atmos. Sci.,42, 2381–2403.

  • Elting, D., and R. A. Brown, 1980: Roll vortices in the planetary boundary layer. A review. Bound.-Layer Meteor.,65, 215–248.

  • Gal-Chen, T., and R. C. Somerville, 1975: On the use of a coordinate transformation for the solution of the Navier–Stokes equations. J. Comput. Phys.,17, 209–228.

  • Golding, B. W., 1993: A numerical investigation of tropical island thunderstorms. Mon. Wea. Rev.,121, 1417–1433.

  • Ikawa, M., and K. Saito, 1991: Description of a nonhydrostatic model developed at the Forecast Research Department of the MRI. Tech. Rep. MRI, 28, 238 pp. [Available from Meteorological Research Institute, Tsukuba, Ibaraki 305-0052, Japan.].

  • ——, H. Mizuno, T. Matsuo, M. Murakami, Y. Yamada, and K. Saito, 1991: Numerical modeling of the convective snow cloud over the Sea of Japan. Precipitation mechanism and sensitivity to ice crystal nucleation rates. J. Meteor. Soc. Japan,69, 641–667.

  • Kato, T., 1995: Box-Lagrangian rain-drop scheme. J. Meteor. Soc. Japan,73, 241–245.

  • ——, 1996: Hydrostatic and non-hydrostatic simulations of the 6 August 1993 Kagoshima torrential rain. J. Meteor. Soc. Japan,74, 355–363.

  • ——, 1998: Numerical simulation of a band-shaped torrential rain observed over southern Kyushu, Japan on 1 August 1993. J. Meteor. Soc. Japan,76, 97–128.

  • Keenan, T. D., G. J. Holland, M. J. Manton, and J. Simpson, 1989: The Island Thunderstorm Experiment (ITEX)—A study of tropical thunderstorms in the maritime continent. Bull. Amer. Meteor. Soc.,70, 152–159.

  • ——, B. R. Morton, X. S. Zhang, and K. Nyguen, 1990: Some characteristics of thunderstorms over Bathurst and Melville Island near Darwin, Australia. Quart. J. Roy. Meteor. Soc.,116, 1153–1172.

  • ——, B. Ferrier, and J. Simpson, 1994a: Development and structure of a maritime continent thunderstorm. Meteor. Atmos. Phys.,53, 185–222.

  • ——, and Coauthors, 2000: The Maritime Continent Thunderstorm Experiment (MCTEX): Overview and some results. Bull. Amer. Meteor. Soc.,81, 2433–2455.

  • Kimura, F., and M. Toda, 1995: A linear model of anabatic and katabatic flows using the terrain following coordinates system. Research Report on Heat and Moisture Transport on Complex Terrain, University of Tsukuba, 71–86. [Available from University of Tsukuba, Ibaraki, 305-0004, Japan.].

  • Kingsmill, D. E., 1995: Convection initiation associated with a sea-breeze front, a gust front, and their collision. Mon. Wea. Rev.,123, 2913–2933.

  • Klemp, J. B., and R. B. Wilhelmson, 1978: The simulation of three-dimensional convective storm dynamics. J. Atmos. Sci.,35, 1070–1096.

  • Kondo, J., 1975: Air-sea bulk transfer coefficients in diabatic conditions. Bound.-Layer Meteor.,9, 91–112.

  • ——, 1976: Heat balance of the East China Sea during the air mass transformation experiment. J. Meteor. Soc. Japan,54, 382–398.

  • Lin, Y. H., R. D. Farley, and H. D. Orville, 1983: Bulk parameterization of the snow field in a cloud model. J. Climate Appl. Meteor.,22, 1065–1092.

  • Liu, C., and M. W. Moncrieff, 1996: A numerical study of the effects of ambient flow and shear on density currents. Mon. Wea. Rev.,124, 2282–2303.

  • Mason, P. J., 1989: Large-eddy simulation of the convective boundary layer. J. Atmos. Sci.,46, 1492–1516.

  • McDonald, A., 1986: A semi-Lagrangian and semi-implicit two time-level integration scheme. Mon. Wea. Rev.,114, 824–833.

  • Ogura, M., 1969: A direct method of Poisson equation by dimension reduction method. J. Meteor. Soc. Japan,47, 319–323.

  • Ohno, H., and M. Isa, 1984: A. statistical relation between GMS-viewed cloud amount and relative humidity (in Japanese). Tenki,31, 493–495.

  • Orlanski, I., 1976: A simple boundary condition for unbounded hyperbolic flows. J. Comput. Phys.,21, 251–269.

  • Pielke, R. A., 1974: A three-dimensional numerical model of the sea breezes over south Florida. Mon. Wea. Rev.,102, 115–139.

  • Puri, K., G. S. Dietachmayer, G. A. Mills, N. E. Davidson, R. A. Bowen, and L. W. Logan, 1998: The new BMRC limited area prediction system, LAPS. Austr. Meteor. Mag.,47, 203–223.

  • Ramage, C. S., 1968: Role of a “Maritime Continent” in the atmospheric circulation. Mon. Wea. Rev.,96, 365–370.

  • Redelsperger, J. L., and Coauthors, 2000: A GCSS model intercomparison for a tropical squall line observed during TOGA-COARE. Part 1: Cloud-resolving models. Quart. J. Roy. Meteor. Soc.,126, 823–863.

  • Saito, K., 1993: A numerical study of the local downslope wind“Yamaji-kaze” in Japan. Part 2: Non-linear aspect of the 3-D flow over a mountain range with a col. J. Meteor. Soc. Japan,71, 247–271.

  • ——, 1994: A numerical study of the local downslope wind “Yamaji-kaze” in Japan. Part 3: Numerical simulation of the 27 September 1991 windstorm with a non-hydrostatic multi-nested model. J. Meteor. Soc. Japan,72, 301–329.

  • ——, 1997: Semi-implicit fully compressible version of the MRI mesoscale nonhydrostatic model—Forecast experiment of the 6 August 1993 Kagoshima Torrential rain. Geophys. Mag. Ser.2,2, 109–137.

  • ——, 1998: Recent modification of the MRI-NHM—Higher order advection schemes and adjustment of mass-flux through lateral boundaries according to the precipitation rates. CAS/JSC WGNE Research Activities in Atmospheric and Oceanic Modelling,27, 5.44–5.45.

  • ——, and M. Ikawa, 1991: A numerical study of the local downslope wind “Yamaji-kaze” in Japan. J. Meteor. Soc. Japan,69, 31–56.

  • ——, and T. Kato, 1999: The Meteorological Research Institute mesoscale nonhydrostatic model (in Japanese). Meteorology Research Note 196, 169–195. [Available from Meteorological Society of Japan, Otemachi, Tokyo 100-0004, Japan.].

  • ——, M. Murakami, T. Matsuo, and H. Mizuno, 1996: Sensitivity experiments on the orographic snowfall over the mountainous region of northern Japan. J. Meteor. Soc. Japan,74, 797–813.

  • ——, T. Kato, H. Eito, and C. Muroi, 2001: Documentation of the Meteorological Research Institute/Numerical Prediction Division unified nonhydrostatic model. Tech. Rep. MRI, in press. [Available from Meteorological Research Institute, Tsukuba, Ibaraki 305-0052, Japan.].

  • Schmidt, H., and U. Schumann, 1989: Coherent structure of the convective boundary layer derived from large-eddy simulations. J. Fluid Mech.,200, 511–562.

  • Seaman, R., W. Bourke, P. Steinle, T. Hart, M. Naughton, and L. Rikus, 1995: Evolution of the Bureau of Meteorology’s Global Assimilation and Prediction system. Part 1: Analysis and initialisation. Austr. Meteor. Mag.,44, 1–18.

  • Segami, A., K. Kurihara, H. Nakamura, M. Ueno, and I. Takano, 1989: Description of Japan Spectral Model. JMA/NPD Technical Rep., 25, 41 pp. [Available from Japan Meteorological Agency, Otemachi, Tokyo 100-8122, Japan.].

  • Simpson, J., T. D. Keenan, B. Ferrier, R. H. Simpson, and G. J. Holland, 1993: Cumulus mergers in the Maritime Continent Region. Meteor. Atmos. Phys.,51, 73–99.

  • Sommeria, G., 1976: Three-dimensional simulation of turbulent processes in an undisturbed trade wind boundary layer. J. Atmos. Sci.,33, 216–241.

  • Stommel, H., 1947: A summary of the theory of convection cells. Ann. N. Y. Acad. Sci.,48, 715–726.

  • Tao, W. K., and J. Simpson, 1984: Cloud interactions and merging: Numerical simulations. J. Atmos. Sci.,41, 2901–2917.

  • ——, and ——, 1989: A further study of cumulus interactions and mergers: Three-dimensional simulation with trajectory analyses. J. Atmos. Sci.,46, 2974–3004.

  • Wakimoto, R. M., 1998: Convectively-driven high wind events. Preprints, 19th Conf. on Severe Local Storms, Minneapolis, MN, Amer. Meteor. Soc., 435.

  • Willis, G. E., and J. W. Deardorff, 1979: Laboratory observations of turbulent penetrative convection platforms. J. Geophys. Res.,84 (C1), 295–301.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1113 790 43
PDF Downloads 314 75 5