An Equivalent Isobaric Geopotential Height and Its Application to Synoptic Analysis and a Generalized ω Equation in σ Coordinates

Qiu-shi Chen Polar Meteorology Group, Byrd Polar Research Center, The Ohio State University, Columbus, Ohio

Search for other papers by Qiu-shi Chen in
Current site
Google Scholar
PubMed
Close
and
David H. Bromwich Polar Meteorology Group, Byrd Polar Research Center, The Ohio State University, Columbus, Ohio

Search for other papers by David H. Bromwich in
Current site
Google Scholar
PubMed
Close
Restricted access

We are aware of a technical issue preventing figures and tables from showing in some newly published articles in the full-text HTML view.
While we are resolving the problem, please use the online PDF version of these articles to view figures and tables.

Abstract

In σ coordinates, a variable ϕe(x, y, σ, t) whose horizontal gradient −ϕe is equal to the irrotational part of the horizontal pressure gradient force is referred to as an equivalent isobaric geopotential height. Its inner part can be derived from the solution of a Poisson equation with zero Dirichlet boundary value. Because −ϕ(x, y, p, t) is also the irrotational part of the horizontal pressure gradient force in p coordinates, the equivalent geopotential ϕe in σ coordinates can be used in the same way as the geopotential ϕ(x, y, p, t) used in p coordinates. In the sea level pressure (SLP) analysis over Greenland, small but strong high pressure systems often occur due to extrapolation. These artificial systems can be removed if the equivalent geopotential ϕe is used in synoptic analysis on a constant σ surface, for example, at σ = 0.995 level. The geostrophic relation between the equivalent geopotential and streamfunction at σ = 0.995 is approximately satisfied.

Because weather systems over the Tibetan Plateau are very difficult to track using routine SLP, 850-hPa, and 700-hPa analyses, equivalent isobaric geopotential analysis in σ coordinates is especially useful over this area. An example of equivalent isobaric geopotential analysis at σ = 0.995 shows that a secondary high separated from a major anticyclone over the Tibetan Plateau when cold air affected the northeastern part of the plateau, but this secondary high is hardly resolved by the SLP analysis. The early stage of a low (or vortex), called a southwest vortex due to its origin in southwest China, over the eastern flank of the Tibetan Plateau is more clearly identified by equivalent isobaric geopotential analysis at σ = 0.825 and 0.735 than by routine isobaric analysis at the 850- and 700-hPa levels. Anomalous high and low systems in the SLP analysis over the Tibetan Plateau due to extrapolation are all removed by equivalent isobaric geopotential analysis at σ = 0.995.

Use of equivalent geopotential ϕe in the vorticity and divergence equations is presented, and the equivalent geopotential equation is derived. These equations can be used in numerical models, initializations, and other dynamical studies. As an example, it is shown how these equations are used to derive a velocity potential form of the generalized ω equation in σ coordinates. As a check, retrieval of precipitation over Greenland using this ω equation shows that the computed precipitation distributions for 1987 and 1988 are in good agreement with the observed annual accumulation.

Corresponding author address: Qiu-shi Chen, Byrd Polar Research Center, The Ohio State University, 1090 Carmack Road, Columbus, OH 43210.

Abstract

In σ coordinates, a variable ϕe(x, y, σ, t) whose horizontal gradient −ϕe is equal to the irrotational part of the horizontal pressure gradient force is referred to as an equivalent isobaric geopotential height. Its inner part can be derived from the solution of a Poisson equation with zero Dirichlet boundary value. Because −ϕ(x, y, p, t) is also the irrotational part of the horizontal pressure gradient force in p coordinates, the equivalent geopotential ϕe in σ coordinates can be used in the same way as the geopotential ϕ(x, y, p, t) used in p coordinates. In the sea level pressure (SLP) analysis over Greenland, small but strong high pressure systems often occur due to extrapolation. These artificial systems can be removed if the equivalent geopotential ϕe is used in synoptic analysis on a constant σ surface, for example, at σ = 0.995 level. The geostrophic relation between the equivalent geopotential and streamfunction at σ = 0.995 is approximately satisfied.

Because weather systems over the Tibetan Plateau are very difficult to track using routine SLP, 850-hPa, and 700-hPa analyses, equivalent isobaric geopotential analysis in σ coordinates is especially useful over this area. An example of equivalent isobaric geopotential analysis at σ = 0.995 shows that a secondary high separated from a major anticyclone over the Tibetan Plateau when cold air affected the northeastern part of the plateau, but this secondary high is hardly resolved by the SLP analysis. The early stage of a low (or vortex), called a southwest vortex due to its origin in southwest China, over the eastern flank of the Tibetan Plateau is more clearly identified by equivalent isobaric geopotential analysis at σ = 0.825 and 0.735 than by routine isobaric analysis at the 850- and 700-hPa levels. Anomalous high and low systems in the SLP analysis over the Tibetan Plateau due to extrapolation are all removed by equivalent isobaric geopotential analysis at σ = 0.995.

Use of equivalent geopotential ϕe in the vorticity and divergence equations is presented, and the equivalent geopotential equation is derived. These equations can be used in numerical models, initializations, and other dynamical studies. As an example, it is shown how these equations are used to derive a velocity potential form of the generalized ω equation in σ coordinates. As a check, retrieval of precipitation over Greenland using this ω equation shows that the computed precipitation distributions for 1987 and 1988 are in good agreement with the observed annual accumulation.

Corresponding author address: Qiu-shi Chen, Byrd Polar Research Center, The Ohio State University, 1090 Carmack Road, Columbus, OH 43210.

Save
  • Alley, R. B., and Coauthors, 1993: Abrupt increase in Greenland snow accumulation at the end of the Younger Dryas event. Nature,362, 527–529.

  • Arakawa, A., and V. R. Lamb, 1977: Computational design of the basic dynamical process of the UCLA general circulation model. Methods Comput. Phys.,17, 173–265.

  • Bijlsma, S. J., L. M. Hafkensheid, and P. Lynch, 1986: Computation of the streamfunction and velocity potential and reconstruction of the wind field. Mon. Wea. Rev.,114, 1547–1551.

  • Bromwich, D. H., and F. M. Robasky, 1993: Recent precipitation variations over the polar ice sheets. Meteor. Atmos. Phys.,51, 259–274.

  • ——, ——, R. A. Keen, and J. F. Bolzan, 1993: Modeled variations of precipitation over the Greenland Ice Sheet. J. Climate,6, 1253–1268.

  • Chen, Q.-S., and Y.-H. Kuo, 1992a: A harmonic-sine series expansion and its application to the partitioning and reconstruction problem in a limited area. Mon. Wea. Rev.,120, 91–112.

  • ——, and ——, 1992b: A consistency condition for the wind field reconstruction in a limited area and a harmonic-cosine series expansion. Mon. Wea. Rev.,120, 2653–2670.

  • ——, ——, and D. H. Bromwich, 1996: A balanced ageostrophic initialization with a fixed external wind lateral boundary value for limited-area models. J. Meteor. Soc. Japan.,74, 325–342.

  • ——, L. Bai, and D. H. Bromwich, 1997a: A harmonic-Fourier spectral limited-area model with an external wind lateral boundary condition. Mon. Wea. Rev.,125, 143–167.

  • ——, D. H. Bromwich, and L. Bai, 1997b: Precipitation over Greenland retrieved by a dynamic method and its relation to cyclonic activity. J. Climate,10, 839–870.

  • Corby, G. A., A. Gilchrist, and R. L. Newson, 1972: A general circulation model of the atmosphere suitable for long period integration. Quart. J. Roy. Meteor. Soc.,98, 809–832.

  • DiMego, G. J., and L. F. Bosart, 1982: The transformation of Tropical Storm Agnes into an extratropical cyclone. Part I: The observed fields and vertical motion computations. Mon. Wea. Rev.,110, 385–411.

  • Gary, J. M., 1973: Estimate of truncation error in transformed coordinate, primitive equation atmospheric models. J. Atmos. Sci.,30, 223–233.

  • Harrison, L. P., 1957: Report on the problem of “reduction of pressure.” Report of the Working Group on Barometry to the World Meteorological Organization, Commission for Instruments and Methods of Observation (CIMO-II), 118 pp.

  • ——, 1970: Reduction of surface pressure to functions useful in analysis and forecasting. Meteorological Observations and Instrumentation, Meteor. Monogr., No. 33, Amer. Meteor. Soc., 121–136.

  • Hawkins, H. F., and S. L. Rosenthal, 1965: On the computation of streamfunction from the wind field. Mon. Wea. Rev.,93, 245–252.

  • Higgins, R. W., K. C. Mo, and S. D. Schubert, 1996: The moisture budget of the central United States in spring as evaluated in the NCEP/NCAR and the NASA/DAO reanalysis. Mon. Wea. Rev.,124, 939–963.

  • Hill, C. D., 1993: Forecast problems in the western region of the National Weather Service: An overview. Wea. Forecasting,8, 158–165.

  • Hirschberg, P. A., and J. M. Fritsch, 1991: Tropopause undulations and the development of extratropical cyclones. Part II: Diagnostic analysis and conceptual model. Mon. Wea. Rev.,119, 518–550.

  • Iversen, T., and T. E. Nordeng, 1984: A hierarchy of nonlinear filtered models—Numerical solutions. Mon. Wea. Rev.,112, 2048–2059.

  • Johnsen, S. J., and Coauthors, 1992: Irregular glacial interstadials recorded in a new Greenland ice core. Nature,359, 311–313.

  • Kalnay, E., and Coauthors, 1996: The NCEP/NCAR 40-Year Reanalysis Project. Bull. Amer. Meteor. Soc.,77, 437–471.

  • Knutson, T. R., and K. M. Weickmann, 1987: 30–60 day atmospheric oscillations: Composite life cycles of convection and circulation anomalies. Mon. Wea. Rev.,115, 1407–1436.

  • Krishnamurti, T. N., 1968a: A diagnostic balance model for studies of weather systems of low and high latitudes. Rossby number less than 1. Mon. Wea. Rev.,96, 197–207.

  • ——, 1968b: A study of a developing wave cyclone. Mon. Wea. Rev.,96, 208–217.

  • ——, 1985: Summer monsoon experiment—A review. Mon. Wea. Rev.,113, 1590–1626.

  • ——, and Y. Ramanathan, 1982: Sensitivity of monsoon onset to differential heating. J. Atmos. Sci.,39, 1290–1306.

  • Kuo, Y.-H., L. Cheng, and R. A. Anthes, 1986: Mesoscale analyses of the Sichuan flood catastrophe, 11–15 July 1981. Mon. Wea. Rev.,114, 1984–2003.

  • Lynch, P., 1988: Deducing the wind from vorticity and divergence. Mon. Wea. Rev.,116, 86–93.

  • ——, 1989: Partitioning the wind in a limited domain. Mon. Wea. Rev.,117, 1492–1500.

  • Ma, K. Y., and L. F. Bosart, 1987: A synoptic overview of a heavy rain event in southern China. Wea. Forecasting,2, 89–112.

  • Mesinger, F., Z. Janjic, S. Nickovic, D. Gavrilov, and D. G. Deaven, 1988: The step mountain coordinate: Model description and performance for cases of Alpine lee cyclogenesis and for a case of an Appalachian redevelopment. Mon. Wea. Rev.,116, 1493–1518.

  • Nakamura, H., 1978: Dynamical effects of mountains on the general circulation of atmosphere. Part I: Development of finite-difference scheme suitable for incorporating mountains. J. Meteor. Soc. Japan,56, 317–340.

  • Ohmura, A., and N. Reeh, 1991: New precipitation and accumulation maps for Greenland. J. Glaciol.,37, 140–148.

  • Palmén, E., and C. W. Newton, 1969: Atmospheric Circulation Systems. Academic Press, 419 pp.

  • Pauley, P. M., and S. J. Nieman, 1992: A comparison of quasigeostrophic and nonquasigeostrophic vertical motions for a model-simulated rapidly intensifying marine extratropical cyclone. Mon. Wea. Rev.,120, 1108–1134.

  • Phillips, N. A., 1957: A coordinate system having some special advantages for numerical forecasting. J. Meteor.,14, 184–185.

  • Raisanen, J., 1995: Factors affecting synoptic-scale vertical motion:A statistical study using a generalized omega equation. Mon. Wea. Rev.,123, 2447–2460.

  • Sangster, W. E., 1960: A method of representing the horizontal pressure force without reduction of station pressure to sea level. J. Meteor.,17, 166–176.

  • ——, 1987: An improved technique for computing the horizontal pressure-gradient force at the earth’s surface. Mon. Wea. Rev.,115, 1358–1369.

  • Shukla, J., and K. R. Saha, 1974: Computation of non-divergent streamfunction and irrotational velocity potential from the observed winds. Mon. Wea. Rev.,102, 419–425.

  • Shuman, F. G., 1957: On the problem of comparing station pressures at varying elevations. Joint Numerical Weather Prediction Unit Office Note 7, 10 pp. [Available from National Centers for Environmental Prediction, National Oceanic and Atmospheric Administration, 4700 Silver Hill Rd., Stop 9910, Washington, DC 20233.].

  • Sundqvist, H., 1976: On vertical interpolation and truncation in connection with use of sigma models. Atmosphere,14, 37–52.

  • Tarbell, T. C., T. T. Warner, and R. A. Anthes, 1981: An example of the initialization of the divergent wind component in a mesoscale numerical weather prediction model. Mon. Wea. Rev.,109, 77–95.

  • Taylor, K. C., G. W. Lamorey, G. A. Doyle, R. B. Alley, P. M. Grootes, P. A. Mayewski, J. W. C. White, and L. K. Barlow, 1993: The“flickering switch” of late Pleistocene climate change. Nature,361, 432–436.

  • Wang, W., Y.-H. Kuo, and T. T. Warner, 1993: A diabatically driven mesoscale vortex in the lee of the Tibetan Plateau. Mon. Wea. Rev.,121, 2542–2561.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1025 690 256
PDF Downloads 302 58 6