Schultz, D. M., and C. A. Doswell III, 1999: Conceptual models of upper-level frontogenesis in southwesterly and northwesterly flow. Quart. J. Roy. Meteor. Soc., 125, 2535-2562. (October 1999, Part A) [QJRMS] [HTML] [PDF]


<font size="+3"><b>Conceptual Models of<br> Upper-Level Frontogenesis<br> in Southwesterly and Northwesterly Flow</font></b>

Conceptual Models of
Upper-Level Frontogenesis
in Southwesterly and Northwesterly Flow

David M. Schultz and Charles A. Doswell III

National Severe Storms Laboratory, National Oceanic and Atmospheric Administration, USA



Quarterly Journal of the Royal Meteorological Society
Submitted February 1998; Revised August 1998 and December 1998

22 December 1998



Corresponding author:
Dr. David M. Schultz
NOAA/National Severe Storms Laboratory
1313 Halley Circle
Norman, OK 73069 USA
email: david.schultz@noaa.gov


NOTE: This document does not contain 15 of the 16 figures. Download the PDF for the figures.


  Summary

The Shapiro (1982) conceptual model as it is applied to the evolution of an upper-level frontal zone within a baroclinic wave is reviewed and its limitations are investigated through previous literature and two case studies presented herewithin. The early stages in the evolutions of these two cases are used to examine specific limitations of this conceptual model: (1) upper-level frontogenesis in southwesterly flow that evolves from a state of equivalent barotropy to a state of cold advection along the front, and (2) upper-level frontogenesis in northwesterly flow with along-front variation in the sign of the thermal advection, such that warm advection occurs upstream of cold advection in the thermal trough.

Vector-frontogenesis diagnostics for the Lagrangian rate of change of the magnitude and direction of the horizontal potential temperature gradient, including tilting due to vertical motion, are derived. These diagnostics are applied to the two cases to examine the maintenance of the potential temperature gradient and the development of cold advection along each upper-level front. The upper-level front in southwesterly (northwesterly) flow was maintained primarily by deformation (tilting) frontogenesis, in agreement with previous research. The increasing cold advection along the upper-level front in both cases was related to an upstream vorticity maximum. For the case in southwesterly flow, the preexisting vorticity maximum approached a downstream equivalent-barotropic upper-level front in a manner similar to an instant occlusion, resulting in cold advection along the length of the upper-level front. For the case in northwesterly flow, an intensifying vorticity maximum concentrated the cold advection in the base of the thermal trough, as warm advection developed upstream.

These two cases are compared to upper-level fronts in previous literature and a climatology of upper-level fronts associated with landfalling cyclones over the eastern North Pacific Ocean. The results indicate that these two cases are typical of early evolutions of upper-level fronts that can occur in southwesterly and northwesterly flow. Therefore, a revised version of the Shapiro conceptual model is presented that more accurately represents the early evolutions exhibited in the present and previous studies.

1  Introduction

As noted in the review by Keyser and Shapiro (1986), baroclinic zones in the upper troposphere and lower stratosphere associated with the polar-front jet stream (known as upper-level fronts) apparently were observed first by Bjerknes and Palmén (1937). More studies followed that examined the structure and dynamics of upper-level fronts in more detail (e.g., Palmén and Nagler 1949; Berggren 1952; Newton and Carson 1953; Reed and Sanders 1953; Newton 1954; Reed 1955; Reed and Danielsen 1959; Staley 1960; Newton and Persson 1962; Danielsen 1964; Newton 1965; Bosart 1970; Shapiro 1970). The results from this early research can be summarized as follows. (1) Upper-level fronts are typically narrow ( ~ 100 km) in the cross-front direction and extend downward from the tropopause, oftentimes associated with stratospheric air, as indicated by high values of potential vorticity and ozone mixing ratio. (2) These upper-level fronts are typically related to regions of localized maxima of wind speed (i.e., jet streaks, as first defined by Newton and Carson (1953, p. 325); hence these systems are sometimes also referred to as upper-level jet-fronts). (3) In most cases, the intensification of the horizontal potential temperature gradient within these upper-level fronts (i.e., frontogenesis) and the vertical extent of the baroclinity are attributed primarily to isentropic tilting by horizontal gradients of vertical motion (i.e., subsidence maximized on the warm side of the upper-level front leads to frontogenesis). In other cases, upper-level frontogenesis can occur due to horizontal confluent deformation, with tilting being frontolytic. (4) Although examples of strong frontogenesis in southwesterly flow are known to occur (e.g., Reed 1955; Bosart 1970), the most intense upper-level frontogenesis tends to occur with systems forming in northwesterly flow, where the differential subsidence across the upper-level front is greatest. This last point was later demonstrated by modeling studies based on real and idealized upper-level fronts (e.g., Shapiro 1981; Keyser and Pecnick 1985; Keyser et al. 1986; Reeder and Keyser 1988): synoptic-scale confluence at a jet-entrance region and cold advection along the length of the front yield a gradient in vertical motion across the upper-level front of sufficient magnitude to produce frontogenesis. (The effect of both confluence and cold advection during upper-level frontogenesis has been referred to as the Shapiro effect (Rotunno et al. 1994).)

1(a)  The Shapiro conceptual model

As noted by Keyser and Shapiro (1986, p. 466), data limitations of these early observational studies were primarily responsible for previous ambiguities and inconsistencies in understanding the dynamical interactions between the primary and secondary circulations in upper-level frontal systems. In an effort to link secondary circulations calculated from the Sawyer-Eliassen equation in two-dimensional models with three-dimensional circulations calculated from observational data, Shapiro (1982, Fig. 7) created a schematic of the movement ``of an upper-tropospheric jet-front system through a midlatitude baroclinic wave'' (Shapiro and Keyser 1990, Fig. 10.8 caption), emphasizing the typical thermal-advection pattern for one such life cycle associated with these systems. Although this schematic (Fig. 1; hereafter, the Shapiro conceptual model) was not meant to address the issue of upper-level frontogenesis directly, Lagrangian frontogenesis is implied at different stages in the progression of the jet-front through the baroclinic wave.

Despite the Shapiro conceptual model not having been designed to illustrate frontogenesis directly, the schematic has become firmly entrenched in later scientific literature, attesting to its success in describing many of the general features of the migration of an upper-level jet-front system through a larger-scale baroclinic wave. The schematic has been reproduced in review articles on upper-level fronts (e.g., Keyser and Shapiro 1986, Fig. 19; Shapiro and Keyser 1990, Fig. 10.8) and emulated in textbooks (e.g., Carlson 1991, Fig. 15.6; Bluestein 1993, Fig. 2.82). Despite the successes of the Shapiro conceptual model in interpretation of the secondary circulations at different stages in the life cycle of a baroclinic wave, applications of the Shapiro schematic to frontogenesis, however, appear to be an oversimplification of the structural and kinematic complexities inherent in nature (as discussed further below). The present study endeavors to consider some of these complexities.

The Shapiro conceptual model begins with an upper-level front situated within confluent northwesterly flow (e.g., Namias and Clapp 1949), and isentropes are nearly parallel to isohypses (Fig. 1(a)). This phase of the evolution will be referred to hereafter as the equivalent-barotropic stage. Approximately 24 h later, the front evolves to a state where the trough in the isentropes lags the trough in the isohypses, implying geostrophic cold advection along the length of the front (the cold-advection stage; Fig. 1(b)). Later, the trough in the geopotential-height field closes off and becomes more symmetrical (the closed-low stage; Fig. 1(c)), as discussed by Bell and Keyser (1993) and Bell and Bosart (1994), with leading warm advection and trailing cold advection. Finally, the closed low opens up in confluent southwesterly flow and geostrophic warm advection along the length of the front is implied (the warm-advection stage; Fig. 1(d)).

Although it appears as if the intent of the Shapiro conceptual model was to illustrate the evolution of a jet-front system (e.g., ``t=t°'', ``t=t°+24 h'', etc., in Fig. 1), observational evidence suggests that this schematic evolution typically does not occur exactly in the manner illustrated by Shapiro (1982), particularly over North America. First, cold advection is depicted as maximum at the inflection point in the height field in the cold-advection stage of the Shapiro conceptual model with negligible thermal advection outside of the frontal zone. Nearly all observed upper-level fronts in northwesterly flow in the literature, however, exhibit the strongest cold advection in the base of the thermal trough (Table 1). Weak cold advection, or more typically warm advection, is present along the northwesterly extensions of the fronts (Table 1). Second, as noted by Newton and Carson (1953, p. 325), Newton and Palmén (1963), Cammas and Ramond (1989, p. 2460), Orlanski and Sheldon (1995, Fig. 3, p. 614), and Pyle (1997, pp. 97 and 133), the progression of a jet-front structure from northwesterly flow to the base of the trough to southwesterly flow (Figs. 1(b)-(d)) rarely exists. Since the wind speed within a jet streak is usually much faster than the speed of movement of the jet streak (e.g., Newton and Carson 1953, p. 331; Palmén and Newton 1969, 206-207; Cunningham and Keyser 1996, section 3), the jet-front cannot be tracked as an advective feature of the flow. Oftentimes, as a new jet-front is developing in the downstream southwesterly flow, the upstream jet-front in northwesterly flow remains. A further point is that the reorganization from large-scale diffluence during the cold-advection stage to large-scale confluence during the warm-advection stage (cf. Figs. 1(b),(d)) seems unlikely to occur within the relatively short time period depicted by the Shapiro conceptual model. Finally, the Shapiro conceptual model illustrates the evolution of an upper-level front in northwesterly flow, although as will be shown later in this paper, upper-level fronts in southwesterly flow also may display some similarities (e.g., evolution from equivalent-barotropic to cold-advection stages) to the Shapiro conceptual model at times.

1(b)  Purpose

The purpose of this paper, therefore, is to compare and contrast the early evolution of southwesterly and northwesterly flow upper-level frontogenesis, using two cases as illustrative examples. Specifically, we intend to demonstrate the following points. First, the kinematics and dynamics differ for upper-level fronts in southwesterly and northwesterly flow and have not been elucidated previously for pure southwesterly flow upper-level frontogenesis. Second, the mechanism that leads to the initiation and configuration of strong cold advection along observed upper-level fronts has not been determined. Understanding this mechanism is of importance because cold advection along upper-level fronts is concomitant with stronger tilting and upper-level frontogenesis (the Shapiro effect), processes that can be associated with surface cyclogenesis (e.g., Uccellini et al. 1985, 980-981). Finally, southwesterly flow upper-level frontogenesis appears to be common over the eastern North Pacific Ocean, in contrast to northwesterly flow frontogenesis, which appears to be common over North America.

Section 2 is a derivation of vector frontogenesis, including the vertical-velocity (tilting) terms. Section 3introduces the two cases of upper-level frontogenesis in southwesterly and northwesterly flow, respectively. Then, the changes in the magnitude and direction of the horizontal potential temperature gradient along the two upper-level fronts are diagnosed using vector frontogenesis. Section 4 is a climatology of upper-level fronts associated with landfalling cyclones on the west coast of North America, showing that evolutions similar to that occurring in southwesterly flow are not an uncommon occurrence. Finally, section 5 presents a concluding discussion and a revised conceptual model for the early evolution of southwesterly and northwesterly flow upper-level frontogenesis.

2  Derivation of the two-dimensional vector frontogenesis, including vertical-velocity terms

As originally defined by Petterssen (1936) for surface fronts, frontogenesis is the Lagrangian rate of change of the magnitude of the horizontal potential temperature (q) gradient due to the horizontal wind (VH = ui + vj). Miller (1948) extended Petterssen's definition to fronts in the free atmosphere by defining frontogenesis as the Lagrangian rate of change of the magnitude of the three-dimensional gradient of the potential temperature due to the three-dimensional wind (V = VH + wk = ui + vj + wk), written here in (x, y, p) coordinates. In this paper, we wish to understand the evolution of an upper-level front from a state of near equivalent barotropy to a state of cold advection along the front, involving the rotation of isentropes relative to isohypses. This reconfiguration of the isentropes suggests a useful diagnostic would be the vector-frontogenesis formulation of Keyser et al. (1988), which considers the Lagrangian rate of change of the magnitude and direction of the horizontal potential temperature gradient by the horizontal wind. Since contributions from vertical motion are typically small for surface fronts, but not necessarily so for upper-level fronts, we extend Keyser et al.'s (1988) methodology to include the frontogenetical effects due to vertical motion. (A similar derivation can be found in Lalaurette et al. (1994, 2005-2007).) Therefore, for the purposes of this paper, we define vector frontogenesis, F, as the Lagrangian rate of change of the magnitude and direction of the horizontal potential temperature gradient due to the three-dimensional wind:
F = d
dt
ÑH q,
(1)
where
d
dt
=
t
+u
x


p 
+ v
y


p 
+ w
p
,
ÑH = i
x


p 
+ j
y


p 
.
The subscript p indicates differentiation on an isobaric surface and hereafter will be implicit in this section. Resolving F into natural coordinates (s,n) such that the s axis is locally tangent to an isentrope and the n axis points towards colder air on a constant pressure surface (i.e., s points in the same direction as the thermal wind):
d
dt
ÑH q = n æ
ç
è
n · d
dt
ÑHq ö
÷
ø
+ s æ
ç
è
s · d
dt
ÑHq ö
÷
ø
,
(2)
where
n = - |ÑH q|-1 ÑH q,
(3a)
s = n ×k .
(3b)
Equation (1) becomes
F = Fnn + Fss ,
(4)
where Fn is referred to as the scalar frontogenesis and Fs is referred to as the rotational frontogenesis (Keyser et al. 1988, (2.3a) and (2.3b)):
Fn = - d
dt
|ÑH q| ,
(5a)
Fs = n · æ
ç
è
k × d
dt
ÑH q ö
÷
ø
.
(5b)
Therefore, positive Fn implies frontolysis, whereas positive Fs implies cyclonic rotation of the isentropes.

Conservation of potential temperature in three-dimensional adiabatic flow means that

q
t
+ V ·Ñ q = 0 ,
(6)
where
Ñ = i
x
+ j
y
+ k
p
.
Following Keyser et al. (1988, p. 764), differentiating (6) respectively by x and y yields:
d
dt
æ
ç
è
q
x
ö
÷
ø
= - u
x
q
x
- v
x
q
y
- w
x
q
p
,
(7a)
d
dt
æ
ç
è
q
y
ö
÷
ø
= - u
y
q
x
- v
y
q
y
- w
y
q
p
.
(7b)
Equations (5a) and (5b) can be rewritten:
Fn = - 1
|ÑH q|
ì
í
î
q
x
d
dt
æ
ç
è
q
x
ö
÷
ø
+ q
y
d
dt
æ
ç
è
q
y
ö
÷
ø
ü
ý
þ
,
(8a)
Fs = - 1
|ÑH q|
ì
í
î
q
y
d
dt
æ
ç
è
q
x
ö
÷
ø
- q
x
d
dt
æ
ç
è
q
y
ö
÷
ø
ü
ý
þ
.
(8b)
Combining (7a) and (7b) with (8a) and (8b) leads to the following expressions
Fn = - 1
|ÑH q|
ì
í
î


[(q)/( x)]( -[(u)/( x)][(q)/( x)]-[(v)/( x)][(q)/( y)])
n1 
+

[(q)/( y)]( -[(u)/( y)][(q)/( x)]-[(v)/( y)][(q)/( y)])
n2 
+

[(q)/( p)]( -[(w)/( x)][(q)/( x)]-[(w)/( y)][(q)/( y)])
n3 
ü
ý
þ
,
(9a)
Fs = - 1
|ÑH q|
ì
í
î


[(q)/( y)]( -[(u)/( x)][(q)/( x)]-[(v)/( x)][(q)/( y)])
s1 
-

[(q)/( x)]( -[(u)/( y)][(q)/( x)]-[(v)/( y)][(q)/( y)])
s2 
+

[(q)/( p)]( [(w)/( y)][(q)/( x)]-[(w)/( x)][(q)/( y)])
s3 
ü
ý
þ
.
(9b)

Terms n1, n2, s1, and s2 represent similar terms in Keyser et al.'s (1988) derivation. Term n3 represents tilting scalar frontogenesis due to vertical-velocity gradients acting on the horizontal potential temperature gradient (equivalent to the sum of terms 4 and 8 in Bluestein (1986, p. 181; 1993, p. 253)) and term s3 represents tilting rotational frontogenesis due to vertical-velocity gradients acting on the horizontal potential temperature gradient.1 Terms n3 and s3 can also be expressed as:

n3
=
- q
p
( ÑH w·ÑHq)
=
- q
p
æ
ç
è
w
n
ö
÷
ø
,
(10a)
s3
=
- q
p
k ·( ÑH w×ÑH q)
=
- q
p
æ
ç
è
w
s
ö
÷
ø
,
(10b)
respectively. These expressions for n3 and s3 exhibit the behavior that when ÑH w is in the same direction as ÑH q, the tilting scalar frontogenesis, n3, is maximized, but the tilting rotational frontogenesis, s3, is zero (Fig. 2(a)). In contrast, when ÑH w is perpendicular to ÑH q, n3 is zero, but |s3| is maximized (Fig. 2(b)).

Following Keyser et al. (1988, 764-765), expressions similar to their (2.10a) and (2.10b) result:

Fn = 1
2
|ÑH q| (ÑH ·VH - E cos2 b) - q
p
( ÑHw·ÑH q),
(11a)
Fs = 1
2
|ÑH q| (k ·ÑH ×VH + E sin2 b) - q
p
k ·( ÑH w×ÑH q),
(11b)
where E is the resultant deformation and b is the local angle between an isentrope and the axis of dilatation. Equation (11a) is comparable to similar expressions in Saucier (1955, 363-365). This expression for Fn comprises three terms related to divergence, deformation, and tilting, whereas (11b) for Fs comprises three terms related to relative vorticity, deformation, and tilting. These partitions will be useful for examining the evolution of the fronts in the two observed cases.2

3  Case studies

In this section, we present two case studies of the early stages in the evolution of upper-level fronts, in southwesterly and northwesterly upper-level flow, respectively. Section 3(a)describes the mesoscale-model simulation used to examine these two cases, which are introduced in section 3(b). The cases are further analyzed in section 3(c) using the vector-frontogenesis diagnostics developed in section 2.

3(a)  Mesoscale-model description

The Pennsylvania State University-National Center for Atmospheric Research (PSU-NCAR) mesoscale model version 5 (MM5), a nonhydrostatic, primitive-equation, s-coordinate model (Dudhia 1993; Grell et al. 1994), was employed to simulate the two upper-level frontogenesis events. Since both the southwesterly and northwesterly flow events occurred nearly simultaneously, one simulation was able to accommodate analysis of both. The simulation was initialized at 1200 UTC 12 December 1988, was ended at 1200 UTC 14 December 1988, and featured 23 variably spaced half-s levels in the vertical. A 30-km horizontal-resolution domain was nested within a 90-km horizontal-resolution domain using a two-way interactive mesh-refinement scheme. Precipitation processes were parameterized using an explicit-moisture scheme that includes prognostic equations for water vapor, cloud water, rain water, cloud ice, and snow (Dudhia 1989; Grell et al. 1994, section 5.3.1.1). The Kain-Fritsch cumulus parameterization (Kain and Fritsch 1993) was used to represent subgrid-scale convective precipitation. Other parameterizations included a multilevel planetary boundary layer (Zhang and Anthes 1982) and a radiative upper boundary condition (Klemp and Durran 1983).

Four-dimensional data assimilation was used throughout the simulation on the 90-km domain and during the first 12 h on the 30-km domain. The assimilation technique (Stauffer and Seaman 1990) employed Newtonian nudging to relax the model simulation to gridded 12-h upper-level and 3-h surface analyses. To create the upper-level analyses, National Meteorological Center (NMC, now known as the National Centers for Environmental Prediction) analyses were interpolated to the model grid. Surface and upper-air observations were then incorporated into the analysis using a Cressman-type analysis scheme (Benjamin and Seaman 1985). After the removal of superadiabatic lapse rates below 500 hPa, the analysis was interpolated to the s-coordinate system, and the integrated mean divergence was removed to avoid the production of spurious gravity waves. Three-hour surface analyses were generated similarly, with first-guess analysis fields provided by linear interpolation of 12-h NMC analyses. Lateral boundary conditions for the 90-km domain were generated by linear interpolation of the 12-h analyses. Throughout this paper, only output from the 90-km domain is displayed.

3(b)  Overview

The first case is the upper-level front associated with the landfalling Pacific cyclone over western North America on 12-13 December 1988. Our analyses are carried out at 500 hPa, a common midtropospheric level analyzed by synoptic meteorologists. While our results would be quantitatively different at upper-tropospheric levels (say, 400 or 300 hPa), our results do not differ qualitatively were another level presented.

Unlike other cases of upper-level fronts in southwesterly flow that apparently evolved from a northwesterly flow regime (e.g., Reed 1955; Bosart 1970; Sanders et al. 1991; Pyle 1997, p. 133), this case occurred entirely within southwesterly flow, and will hereafter be referred to as SW. SW was associated with a 500-hPa closed low that moved slowly into the Gulf of Alaska on 11 December 1988 (not shown). NMC analyses (not shown) indicate remnant baroclinity on the southeast side of the closed low in southwesterly flow, with the geopotential-height (hereafter, height) and potential-temperature contours nearly parallel from 0000 UTC 11 December 1988 to 1200 UTC 12 December 1988 (hereafter 12/12) (Fig. 3(a)). Unlike the equivalent-barotropic stage in the Shapiro conceptual model (Fig. 1(a)), however, the upper-level front was strongest in the southwesterly, not the northwesterly, geostrophic flow. The enhanced baroclinity in SW was also associated with a preexisting progressive relative-vorticity (hereafter, vorticity) maximum at the base of the shortwave trough, coincident with the thermal trough at this level (Fig. 3(a)). Twelve hours later (13/00; Fig. 3(b)), the orientation between the isentropes and isohypses implied geostrophic cold advection (potential-temperature advection due to total horizontal wind contoured in thick lines) along much of the length of the 500-hPa front and vorticity. This structure represents the cold-advection stage. Therefore, this case illustrates the same evolution from a nearly equivalent-barotropic stage to a cold-advection stage as in the Shapiro conceptual model (Figs. 1(a),(b)), except in southwesterly, instead of northwesterly, flow. It is curious to note that Hines and Mechoso (1991, Figs. 3-5) simulated a similar structure in their idealized model of primitive-equation upper-level frontogenesis: cold advection along a baroclinic zone developing in southwesterly flow. Further investigation of this case was abandoned because ``This simulated frontogenesis, therefore, has substantial differences with the observed phenomena [as characterized by the Shapiro conceptual model] and will not be analyzed further in this paper'' (Hines and Mechoso 1991, p. 1232).

Satellite imagery for SW (Fig. 4) shows that a comma cloud is associated with the advection of 500-hPa cyclonic vorticity and the polar-front cloud band is associated with the baroclinic zone in southwesterly flow (Fig. 4(a)). As the comma cloud approaches and merges with the polar-front cloud band, the clouds deform owing to the rotation associated with the vorticity maximum (Figs. 4(b),(c)). This evolution is remarkably similar to the instant-occlusion concept reviewed by Schultz and Mass (1993, section 2(d)) and described by, for example, McGinnigle et al. (1988, Fig. 7(a)), Evans et al. (1994), and Bader et al. (1995, section 4.4, Figs. 4.4.3(c) and 4.4.5(c)). In particular, the approach of a shortwave trough (comma cloud) toward a preexisting baroclinic zone in southwesterly flow (polar-front cloud band) epitomizes the instant occlusion. SW is consistent with these previous studies, in that observed instant occlusions tend to occur in confluent large-scale flows that favor merger between the comma cloud and the polar-front cloud band. As has been noted by the above authors, not all cyclone events that occur in southwesterly flow, however, are instant occlusions.

The other case is the upper-level front associated with the Experiment on Rapidly Intensifying Cyclones over the Atlantic (ERICA) Intensive Observing Period (IOP) 2 cyclone on 13-14 December 1988 over central North America. Because the upper-level front intensified in northwesterly upper-level flow, this case will hereafter be referred to as NW. The ERICA IOP 2 cyclone and its upper-level features have been discussed previously by Sanders (1990), Roebber (1993), Lackmann et al. (1997), Hakim (1997, section 5.1), and Reed and Albright (1997). At 13/00 (Fig. 5(a)), as well as 12 and 24 h previous (not shown), much of the frontal zone was characterized by weak cold advection, with the strongest cold advection on the equatorward side of the front and near the intensifying vorticity maximum. Lackmann et al. (1997) determined that tilting of horizontal vorticity into the vertical by differential subsidence across the upper-level front was responsible for the formation and intensification of this vorticity maximum, which would later be associated with the ERICA IOP 2 cyclogenesis. By 13/12, however, warm advection was occurring along much of the length of the upper-level front (Fig. 5(b)). Unlike in SW, where cold advection occurred along the length of the front, the cold advection in NW was limited to the base of the thermal trough (cf. Figs. 3(b) and 5(b)) associated with the ``compacting'' vorticity maximum (Lackmann et al. 1997). A compacting vorticity maximum is one that undergoes an increase in isotropy during its evolution (i.e., the length of the major axis of the vorticity maximum decreases relative to that of the minor axis) (Lackmann et al. 1997, p. 2731). Therefore, NW, as well as other similar events from the literature discussed previously, indicate significant along-front variability in thermal advection as the upper-level front reaches maturation, in contrast to the depicted cold advection dominating the entire length of the front in the Shapiro conceptual model (Fig. 1(b)).

It is apparent that significant differences exist between the early stages in the evolutions of the two cases of upper-level frontogenesis discussed presently, SW and NW (Figs. 3 and 5, respectively). In particular, the evolution from the equivalent-barotropic stage to cold-advection stage in SW occurred in southwesterly, not northwesterly flow. Also, NW experienced significant along-front variation in the sign of the thermal advection, with cold advection at the leading edge of the vorticity maximum and warm advection along the northwestern portion, rather than cold advection along the length of the front as portrayed in the Shapiro conceptual model (Fig. 1). It might be argued that cases exhibiting these same evolutions may be relatively uncommon, or alternatively, it may be that the Shapiro conceptual model applies best to cases in a limited geographical area (e.g., over North America). As will be shown in section 4, evolutions similar to SW are commonly associated with North Pacific cyclones landfalling on the west coast of North America, which, in turn, resemble cases throughout the Northern Hemisphere; we have already observed that many cases of upper-level fronts in the literature are similar to NW. Therefore, we believe our results based on these two cases to be of some generality.

3(c)  Vector-frontogenesis diagnostics

In this section, we further examine SW and NW using the vector-frontogenesis diagnostics derived in section 2. In order to illustrate the behavior of the vertical tilting terms in F, we present the 500-hPa vertical velocity. At 12/12, the vertical velocity for SW was characterized by descent on the west side of the trough and ascent on the east side of the trough, maximized within the region of strong potential temperature gradient (Fig. 6(a)). This suggests, therefore, that the front was in a region of southwesterly flow unfavorable for further intensification by tilting. Note that the ``four-cell'' pattern of ascent and descent associated with straight jet-front systems is masked by the vertical motions associated with the curvature of the flow (Uccellini 1990, section 6.3.1 and references therewithin). By 13/00, the pattern remained much the same (Fig. 6(b)). In contrast, descent was occurring on the warm side of the frontal zone in NW at 13/00 and 13/12 (Figs. 7(a),(b)), favoring frontogenesis by tilting.

This hypothesis is confirmed by examining 500-hPa -Fn and its components (divergence, deformation, and vertical-tilting) for SW and NW.3 At 12/12 for SW, the divergence term was consistently smaller than the deformation term, which was positive along much of the front (cf. Figs. 8(a),(b)). Figure 8(c) supports the interpretation from Fig. 6(a) that the ascent along the front was unfavorable for frontogenesis through tilting, as negative tilting frontogenesis occurred along the length of the front. As a result, the total scalar frontogenesis was positive along much of the length of the front, primarily supported by deformation (Figs. 8(b),(d)). At 13/00, the divergence term was once again small compared to the deformation term (cf. Figs. 9(a),(b)). Tilting resulted in frontolysis along the front so that the total frontogenesis was primarily a result of the deformation term (Figs. 9(b),(c),(d)), a result consistent with other cases of upper-level frontogenesis in southwesterly flow (e.g., Bosart 1970).

The situation was quite different for NW. At 13/00, the divergence term was smaller than the deformation term, and both were negative along the equatorward side of the front (Figs. 10(a),(b)). Tilting, on the other hand, was positive along much of the front, and dominant along the northwest extent of the front (Figs. 10(c),(d)) in the region of active frontogenesis. This confirms our hypothesis from Fig. 7(a) and the conclusion of Lackmann et al. (1997) that descent along the front results in frontogenesis through tilting. At 13/12, divergent and deformation frontolysis dominated along the poleward side of the front, maximized in the base of the thermal trough and along the northwest extent of the front (Figs. 11(a),(b)). Even though tilting frontogenesis was positive along the length of the front, it was not large enough to offset the negative deformation frontogenesis, with the result that the total scalar frontogenesis was negative, maximized at the base of the thermal trough (Figs. 11(c),(d)).

That the tilting frontogenesis for NW was maximum upstream of the strongest cold advection in a region of warm advection, particularly at 13/12, was noted in an earlier similar case study by Sanders et al. (1991, pp. 1339, 1364-1365). For SW, despite the presence of cold advection and confluence, frontogenetical tilting along the upper-level front did not occur (i.e., the Shapiro effect did not seem to be operating in SW). Uccellini et al. (1985, p. 978) and Pyle (1997, section 7.1) emphasized the importance of along-flow contributions to the vertical circulations around upper-level jet-fronts. This research, therefore, underscores the point made by these previous authors that two-dimensional conceptualizations of jet-front circulations can be difficult to apply to the real atmosphere in regions of curved flow.

We next examine 500-hPa Fs and its components for SW and NW. At 12/12 for SW, the vorticity term in Fs was largest along the vorticity maximum, with the deformation term negative along the northeastern extent of the front (Figs. 12(a),(b)). Positive Fs implies cyclonic rotation of the isentropes, consistent with positive vorticity in that area. Tilting led to negative Fs, anticyclonic rotation of the isentropes, which favored eastward migration of the thermal trough (Fig. 12(c)). Total Fs was a maximum in the base of the thermal trough (Fig. 12(d)), due primarily to the vorticity term (Fig. 12(a)). It was this feature that resulted in the cyclonic rotation of the isentropes relative to the flow field, thereby initiating cold advection along the front in the base of the thermal trough (see footnote 3). At 13/00 after the onset of cold advection, patterns similar to that at 12/12 existed along the front. The dominant term was the vorticity term, maximum along the length of the front (Figs. 13(a),(d)). The deformation and vertical terms were smaller and negative (Figs. 13(b),(c)). Total Fs was thus positive along most of the length of the front (Fig. 13(d)), consistent with the cyclonic rotation of the isentropes leading to cold advection along the length of the front.

For NW at 13/00, the vorticity term was positive along the front (Fig. 14(a)) and the deformation term was negative along the upstream part of the front (Fig. 14(b)), as opposed to the downstream part of the front in SW (Fig. 12(b)). With the negative contribution from the vertical term (Fig. 14(c)), the total Fs was positive at the leading downstream edge of the front (Fig. 14(d)), thereby implying that the rotation of isentropes leading to cold advection occurred here, as opposed to the upstream edge of the front in SW. This pattern was repeated and intensified at 13/12 (Fig. 15). The vorticity term dominated (Fig. 15(a)), but was offset by the deformation term along the northwest extension of the front upstream (Fig. 15(b)). Tilting was weak or negative such that the total Fs implied cyclonic rotation of the isentropes occurring in the base of the thermal trough (Figs. 15(c),(d)), in agreement with the region of cold advection in Fig. 5(b).

To examine these results further, we employed the methodology of Keyser et al. (1989) and Loughe et al. (1995) to calculate the divergent and rotational wind explicitly. Calculations of vector frontogenesis using these wind components (not shown) illustrate that the 500-hPa rotational wind was primarily responsible for the 500-hPa rotational frontogenesis, resulting in the transition to cold advection along the front. Also, the contribution of the divergent wind field to frontogenesis at 500 hPa is small. This confirms our earlier results using the vector-frontogenesis diagnostics on the total wind field.

4  Climatology of upper-level fronts over the eastern North Pacific Ocean

In order to assess the frequency of occurrence of the evolution epitomized by SW, a climatology of upper-level fronts associated with landfalling North Pacific surface cyclones was compiled. The NMC Daily Weather Map series was examined for all cases of North Pacific surface cyclones crossing the western coast of North America between 35°N and 60°N in December, January, and February for the six winters 1988-89 through 1993-94. The NMC 500-hPa analyses associated with each cyclone were then characterized according to two different subjective criteria. First, the 500-hPa flow was identified as southwesterly, northwesterly, zonal, or other (e.g., closed low). Second, the thermal evolution of the 500-hPa baroclinicity was then classified as equivalent-barotropic stage to cold-advection stage (i.e., increasing cold advection), decreasing cold advection, weak advection or equivalent-barotropic stage, warm advection, or a combination of the above. The results from this climatology are presented in Table 2.

Of the 19-30 winter (DJF) cyclones per year that made landfall on the west coast of North America during the six-year period of our climatology, 44% occurred in southwesterly flow versus 14% in northwesterly flow (Table 2), consistent with the wintertime-averaged southwesterly jet stream over the eastern North Pacific Ocean (e.g., Sanders 1988, Fig. 7; Bluestein 1993, Fig. 1.72(a); Lackmann et al. 1996, Fig. 1). With regard to the thermal evolution of the baroclinic zone, the largest percentage (49%) was characterized by weak thermal advection, indicating relatively old, nearly equivalent-barotropic systems, characteristic of cyclones at the end of their evolutions over the eastern North Pacific Ocean. The second largest group (23%) comprised those events that underwent the change from equivalent barotropic to cold advection.

Those 23% (35 events in total) that evolved from a state of equivalent barotropy to a state of cold advection were then examined separately (Table 3). Of these 35 events, 63% occurred in southwesterly flow, or 15% of the total of 149 cyclones. Of the 66 events that occurred in southwesterly flow (Table 2), 33% (22/66) evolved from equivalent barotropy to cold advection (not shown). This is in comparison to the 29% (6/21) of the northwesterly flow cases and the 28% (5/18) of the zonal flow cases that underwent this evolution (not shown). Since the evolution from equivalent barotropic to cold advection occurs at roughly the same rate regardless of flow type, it seems that the predominance of southwesterly flow fronts that evolve from equivalent barotropic to cold advection is favored in this area owing to the climatological tendency for southwesterly flow, rather than some inherent tendency for cold advection to develop from an initial state of equivalent barotropy, preferentially, in southwesterly versus northwesterly or zonal flow. Results from Tables 2 and 3 indicate, therefore, that the SW evolution displayed in Fig. 3, although not the only upper-level-frontal evolution possible in southwesterly flow, is relatively common on the west coast of North America.

It should be noted that this evolution in southwesterly flow from nearly equivalent barotropic to cold advection is not necessarily limited to over the eastern North Pacific Ocean. We would expect similar evolutions to SW to be common in regions where relatively old, nearly equivalent-barotropic systems traverse regions of climatological southwesterly flow (e.g., the eastern North Atlantic Ocean). Indeed, similar evolutions elsewhere have been documented in the literature: Bjerknes (1951, Figs. 11, 12, and 14), McGinnigle et al. (1988, Fig. 7), and Cammas and Ramond (1989, their case JS2, Figs. 2(a) and 19).

Additionally, M. Sinclair (1998, personal communication) has performed composite analyses of cyclogenesis in different large-scale flows over the western South Pacific Ocean. For a composite of cyclones that develop in the equatorward entrance region of an upper-level jet streak (his classification E; Sinclair and Revell 1999), the 500-hPa thermal evolution resembles that of SW, whereas for a composite of cyclones that develop in the diffluent exit region of an upstream upper-level jet streak (his classification U; Sinclair and Revell 1999), the 500-hPa thermal evolution resembles that of NW. These results are consistent with those presented in this paper.

5  Concluding discussion

We have presented the structure and early evolution of two case studies of upper-level frontogenesis. As discussed in sections 1(a) and 4, we believe these cases are reasonably representative of many upper-level fronts. The Shapiro conceptual model suggests that the evolution from equivalent barotropic to cold advection along the length of the front occurs in northwesterly flow. We found that this schematic evolution more closely resembles that which may occur southwesterly flow instead. Northwesterly flow fronts, however, tend to concentrate the cold advection in the base of the thermal trough, in contrast to the Shapiro conceptual model.

5(a)  Revised conceptual model

Therefore, the results of this research lead us to revise the Shapiro conceptual model. Fig. 16(a) presents the schematic evolution of an upper-level front in southwesterly flow over the eastern North Pacific Ocean. Initially, the isotherms and the isohypses are relatively parallel. The advance of an upstream vorticity maximum towards the baroclinic zone leads to the rotation of the isentropes relative to the height contours, resulting in the onset of cold advection. This revised model differs from the original Shapiro conceptual model because this evolution occurs in southwesterly, rather than northwesterly, flow and the vorticity maximum is shown to be responsible for the onset of the cold advection.

In contrast, Fig. 16(b) presents the schematic evolution of an upper-level front in northwesterly flow over North America. Initial cold advection along the length of the front becomes concentrated in the base of the thermal trough in conjunction with an intensifying and compacting vorticity maximum in northwesterly flow. Upstream of the thermal trough, warm advection is typically occurring in the northwesterly flow, indicating substantial along-front variation in thermal advection, in contrast to the Shapiro conceptual model.

5(b)  Thermal advection along upper-level fronts

During the evolution from equivalent barotropy to cold advection in SW, the rotation of the isentropes relative to the height field was initiated by the approach of an upstream vorticity maximum towards the baroclinic zone. As illustrated by the Fs diagnostics, the rotation was attributed to the vorticity contribution to Fs, the vertical and deformation terms being negative or small. Similar rotation of the isentropes was also observed in NW near the base of the thermal trough. Once the rotation of the isentropes by the approaching vorticity maximum is accomplished and thermal advection becomes significant, the dynamical processes associated with that advection likely become significant also. Therefore, it appears that vector frontogenesis is useful in a variety of large-scale flow regimes to identify when this change in structure will occur and provides operational forecasters a methodology to forecast this transition.

A mechanism for the onset of cold advection along upper-level fronts was suggested previously by Rotunno et al. (1994, 3390-3391) and discussed further in Keyser (1998). Based on their analysis of normal-mode cyclogenesis in a baroclinic primitive-equation channel model, Rotunno et al. (1994) hypothesized that the transition to along-front cold advection in the confluent jet-entrance region (the Shapiro effect) was due to the subsidence along the upper-level front bringing down higher potential-temperature air from aloft, thereby implying a horizontal rotation of the isentropes. This mechanism would appear to implicate the vertical tilting term in Fs, s3 in (11b), in contrast to our results that implicate the vorticity term in Fs in both SW and NW. In fact, the tilting term in both SW and NW acted in the opposite direction (anticyclonic rotation of the isentropes) than what was observed to occur (cyclonic rotation of the isentropes) (Figs. 12(c), 13(c), 14(c), and 15(c)).

Furthermore, since SW and NW are similar to other observed cases (as discussed in sections 1(a) and 4), the hypothesis advanced by Rotunno et al. (1994) that tilting is responsible for the onset of cold advection may not be applicable for most cases of observed upper-level frontogenesis. A possible interpretation of this discrepancy may be related to the developing midtropospheric potential-vorticity maximum (Rotunno et al. 1994, Fig. 12(b)). Their configuration is similar to that in NW, where the intensifying vorticity maximum resulted in cyclonic rotational frontogenesis. It therefore seems likely that the primitive-equation simulation presented by Rotunno et al. (1994) could be reinterpreted in terms of the diagnostics presented in this paper: frontogenetical tilting leads to midtropospheric vorticity generation which results in rotation of the isentropes relative to the flow, resulting in cold advection.

We also note that the evolution shown in Rotunno et al. (1994) more closely resembles NW than SW, a point made previously by Lackmann et al. (1997, 2755-2756). Whereas tilting of vorticity greatly intensified a weak vorticity maximum in NW, SW possessed a preexisting vorticity maximum that did not change intensity substantially. Idealized normal-mode model experiments by their nature must generate vorticity maxima from a basic state with very small initial disturbances. This comparison has already been noted by Keyser and Shapiro (1986, p. 494), who state, ``Although realistic reproductions of upper-level and surface fronts are found in the context of baroclinic instability theory, they arise as a consequence of upper-level wave amplification and low-level cyclogenesis, rather than appearing in advance of these processes.'' As we have shown here, the evolution of upper-level fronts is intimately related to vorticity maxima preexisting (SW) or developing concomitantly with the upper-level front (NW). The large-scale flow, therefore, acts to position the upper-level fronts and the vorticity maxima into an appropriate configuration.

5(c)  Future research

An extension of this study would be to perform a climatology, as in section 4, except for events occurring over north-central North America, the region characterized in the Shapiro conceptual model and where previous studies indicate frequent northwesterly flow frontogenesis (Sanders 1988; Lackmann et al. 1996). Another topic of interest is the later evolution of the cases presented in this paper. The later evolution of NW has been documented by Lackmann et al. (1997) and Hakim (1997). For SW, the eventual movement of that front over the intermountain region of the western United States and subsequent reintensification in the lee of the Rockies was complicated by the existence of the mountains.

  Acknowledgments

We are deeply indebted to the following individuals for their contributions to this work: Prof. Jim Steenburgh and Dr. David Stensrud for their assistance using MM5; Prof. Steenburgh and Matt Pyle for discussions on upper-level frontogenesis; and Prof. Daniel Keyser, Dr. Rich Rotunno, Prof. Gary Lackmann, Matt Wandishin, Dr. John Cortinas, Dr. Greg Hakim, Dr. Mel Shapiro, Dr. Conrad Ziegler, and an anonymous reviewer for their comments on an earlier version of this manuscript. Dr. Mark Sinclair graciously performed composite analyses of his cyclone dataset at our request. Dr. Hakim suggested the analysis using the Loughe et al. (1995) methodology, Dr. David Knight provided the source code, and Matt Pyle assisted with its implementation. The satellite imagery presented in Fig. 4 was obtained courtesy of the National Environmental Satellite Data and Information Service and the National Climatic Data Center. Access to the Storm Prediction Center's archived weather maps on microfilm is gratefully acknowledged. We are grateful to the Data Support Section of the Scientific Computing Division of NCAR for providing data used in this study, and the Mesoscale and Microscale Meteorology Division of NCAR for their support of MM5.

This research was conducted while the first author was a National Research Council Postdoctoral Research Associate at the National Severe Storms Laboratory.

REFERENCES



Bader, M. J., G. S. Forbes, J. R. Grant, R. B. E. Lilley, and A. J. Waters, 1995: Images in Weather Forecasting. Cambridge University Press, 499 pp.

Bell, G. D., and L. F. Bosart, 1994: Midtropospheric closed cyclone formation over the southwestern United States, the eastern United States, and the Alps. Mon. Wea. Rev., 122, 791-813.

Bell, G. D., and D. Keyser, 1993: Shear and curvature vorticity and potential-vorticity interchanges: Interpretation and application to a cutoff cyclone event. Mon. Wea. Rev., 121, 76-102.

Benjamin, S. G., and N. L. Seaman, 1985: A simple scheme for objective analysis in curved flow. Mon. Wea. Rev., 113, 1184-1198.

Berggren, R., 1952: The distribution of temperature and wind connected with active tropical air in the higher troposphere and some remarks concerning clear air turbulence at high altitude. Tellus, 4, 43-53.

Bjerknes, J., 1951: Extratropical cyclones. Compendium of Meteorology, T. F. Malone, Ed., Amer. Meteor. Soc., 577-598.

Bjerknes, J., and E. Palmén, 1937: Investigations of selected European cyclones by means of serial ascents. Geophys. Publ., 12(2), 1-62.

Bluestein, H. B., 1986: Fronts and jet streaks: A theoretical perspective. Mesoscale Meteorology and Forecasting, P. S. Ray, Ed., Amer. Meteor. Soc., 173-215.

Bluestein, H. B., 1993: Synoptic-Dynamic Meteorology in Midlatitudes. Volume II: Observations and Theory of Weather Systems. Oxford University Press, 594 pp.

Bosart, L. F., 1970: Mid-tropospheric frontogenesis. Quart. J. Roy. Meteor. Soc., 96, 442-471.

Cammas, J.-P., and D. Ramond, 1989: Analysis and diagnosis of the composition of ageostrophic circulations in jet-front systems. Mon. Wea. Rev., 117, 2447-2462.

Carlson, T. N., 1991: Mid-Latitude Weather Systems. Harper Collins, 507 pp.

Cunningham, P., and D. Keyser, 1996: Numerical modelling of jet-streak dynamics. Preprints, Seventh Conf. on Mesoscale Processes, Reading, United Kingdom, Amer. Meteor. Soc., 20-22.

Danielsen, E. F., 1964: Project Springfield report. DASA 1517, Defense Atomic Support Agency, 97 pp. [NTIS AD-607980].

Djuri\'c, D., Weather Analysis. Prentice Hall, 304 pp.

Dudhia, J., 1989: Numerical study of convection observed during the Winter Monsoon Experiment using a mesoscale two-dimensional model. J. Atmos. Sci., 46, 3077-3107.

Dudhia, J., 1993: A nonhydrostatic version of the Penn State-NCAR mesoscale model: Validation tests and simulation of an Atlantic cyclone and cold front. Mon. Wea. Rev., 121, 1493-1513.

Evans, M. S., D. Keyser, L. F. Bosart, and G. M. Lackmann, 1994: A satellite-derived classification scheme for rapid maritime cyclogenesis. Mon. Wea. Rev., 122, 1381-1416.

Grell, G. A., J. Dudhia, and D. R. Stauffer, 1994: A description of the fifth-generation Penn State/NCAR Mesoscale Model (MM5). NCAR Tech. Note NCAR/TN-398+STR, 138 pp. [Available from NCAR, P.O. Box 3000, Boulder, CO 80307-3000 USA.]

Hakim, G. J., 1997: Extratropical cyclogenesis in terms of baroclinic vortex dynamics. Ph.D. dissertation, University at Albany, State University of New York, 210 pp. [Available from Department of Earth and Atmospheric Sciences, SUNY Albany, 1400 Washington Avenue, Albany, NY 12222 USA.]

Hines, K. M., and C. R. Mechoso, 1991: Frontogenesis processes in the middle and upper troposphere. Mon. Wea. Rev., 119, 1225-1241.

Kain, J. S., and J. M. Fritsch, 1993: Convective parameterization for mesoscale models: The Kain-Fritsch scheme. The Representation of Cumulus Convection in Numerical Models. Meteor. Monogr., No. 24, Amer. Meteor. Soc., 165-170.

Keyser, D., 1998: On the representation and diagnosis of frontal circulations in two and three dimensions. The Life Cycles of Extratropical Cyclones, M. A. Shapiro, Ed., Amer. Meteor. Soc., in press.

Keyser, D., and M. J. Pecnick, 1985: A two-dimensional primitive equation model of frontogenesis forced by confluence and horizontal shear. J. Atmos. Sci., 42, 1259-1282.

Keyser, D., and M. A. Shapiro, 1986: A review of the structure and dynamics of upper-level frontal zones. Mon. Wea. Rev., 114, 452-499.

Keyser, D., M. J. Pecnick, and M. A. Shapiro, 1986: Diagnosis of the role of vertical deformation in a two-dimensional primitive equation model of upper-level frontogenesis. J. Atmos. Sci., 43, 839-850.

Keyser, D., M. J. Reeder, and R. J. Reed, 1988: A generalization of Petterssen's frontogenesis function and its relation to the forcing of vertical motion. Mon. Wea. Rev., 116, 762-780.

Keyser, D., B. D. Schmidt, and D. G. Duffy, 1989: A technique for representing three-dimensional vertical circulations in baroclinic disturbances. Mon. Wea. Rev., 117, 2463-2494.

Klemp, J. B., and D. R. Durran, 1983: An upper boundary condition permitting internal gravity wave radiation in numerical mesoscale models. Mon. Wea. Rev., 111, 430-444.

Krishnamurti, T. N., 1968: A study of a developing wave cyclone. Mon. Wea. Rev., 96, 208-217.

Lackmann, G. M., L. F. Bosart, and D. Keyser, 1996: Planetary- and synoptic-scale characteristics of explosive wintertime cyclogenesis over the western North Atlantic Ocean. Mon. Wea. Rev., 124, 2672-26702.

Lackmann, G. M., D. Keyser, and L. F. Bosart, 1997: A characteristic life cycle of upper-tropospheric cyclogenetic precursors during the Experiment on Rapidly Intensifying Cyclones over the Atlantic (ERICA). Mon. Wea. Rev., 125, 2729-2758.

Lalaurette, F., C. Fischer, and J.-P. Cammas, 1994: Location and interaction of upper- and lower-troposphere adiabatic frontogenesis. Mon. Wea. Rev., 122, 2004-2021.

Loughe, A. F., C.-C. Lai, and D. Keyser, 1995: A technique for diagnosing three-dimensional ageostrophic circulations in baroclinic disturbances on limited-area domains. Mon. Wea. Rev., 123, 1476-1504.

McGinnigle, J. B., M. V. Young, and M. J. Bader, 1988: The development of instant occlusions in the North Atlantic. Meteor. Mag., 117, 325-341.

Miller, J. E., 1948: On the concept of frontogenesis. J. Meteor., 5, 169-171.

Namias, J., and P. F. Clapp, 1949: Confluence theory of the high tropospheric jet stream. J. Meteor., 6, 330-336.

Neiman, P. J., and M. A. Shapiro, 1989: Retrieving horizontal temperature gradients and advections from single-station wind profiler observations. Wea. Forecasting, 4, 222-233.

Newton, C. W., 1954: Frontogenesis and frontolysis as a three-dimensional process. J. Meteor., 11, 449-461.

Newton, C. W., 1958: Variations in frontal structure of upper level troughs. Geophysica, 6, 357-375.

Newton, C. W., 1965: Variations in structure of subtropical current system accompanying a deep polar outbreak. Mon. Wea. Rev., 93, 101-110.

Newton, C. W., and J. E. Carson, 1953: Structure of wind field and variations of vorticity in a summer situation. Tellus, 5, 321-339.

Newton, C. W., and A. V. Persson, 1962: Structural characteristics of the subtropical jet stream and certain lower-stratospheric wind systems. Tellus, 14, 221-241.

Newton, C. W., and E. Palmén, 1963: Kinematic and thermal properties of a large-amplitude wave in the westerlies. Tellus, 15, 99-119.

Orlanski, I., and J. P. Sheldon, 1995: Stages in the energetics of baroclinic systems. Tellus, 47A, 605-628.

Palmén, E., and K. M. Nagler, 1949: The formation and structure of a large-scale disturbance in the westerlies. J. Meteor., 6, 227-242.

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

Petterssen, S., 1936: Contribution to the theory of frontogenesis. Geophys. Publ., 11(6), 1-27.

Pyle, M. E., 1997: A diagnostic study of jet streaks: Kinematic signatures and relationship to coherent tropopause disturbances. M.S. thesis, Department of Earth and Atmospheric Sciences, University at Albany, State University of New York, 169 pp. [Available from Department of Earth and Atmospheric Sciences, SUNY Albany, 1400 Washington Avenue, Albany, NY 12222 USA.]

Reed, R. J., 1955: A study of a characteristic type of upper-level frontogenesis. J. Meteor., 12, 226-237.

Reed, R. J., and E. F. Danielsen, 1959: Fronts in the vicinity of the tropopause. Arch. Meteor. Geophys. Bioklimatol., A11, 1-17.

Reed, R. J., and F. Sanders, 1953: An investigation of the development of a mid-tropospheric frontal zone and its associated vorticity field. J. Meteor., 10, 338-349.

Reed, R. J., and M. D. Albright, 1997: Frontal structure in the interior of an intense mature ocean cyclone. Wea. Forecasting, 12, 866-876.

Reeder, M. J., and D. Keyser, 1988: Balanced and unbalanced upper-level frontogenesis. J. Atmos. Sci., 45, 3366-3386.

Riehl, H., J. Badner, J. E. Hovde, N. E. La Seur, L. L. Means, W. C. Palmer, M. J. Schroeder, L. W. Snellman, and others, 1952: Forecasting in Middle Latitudes, Meteor. Monogr., 1(5), Amer. Meteor. Soc., 80 pp.

Riehl, H., M. A. Alaka, C. L. Jordan, and R. J. Renard, 1954: The Jet Stream. Meteor. Monogr., 2(7), Amer. Meteor. Soc., 100 pp.

Roebber, P. J., 1993: A diagnostic case study of self-development as an antecedent conditioning process in explosive cyclogenesis. Mon. Wea. Rev., 121, 976-1006.

Rotunno, R., W. C. Skamarock, and C. Snyder, 1994: An analysis of frontogenesis in numerical simulations of baroclinic waves. J. Atmos. Sci., 51, 3373-3398.

Sanders, F., 1988: Life history of mobile troughs in the upper westerlies. Mon. Wea. Rev., 116, 2629-2648.

Sanders, F., 1990: Surface analysis over the oceans-Searching for sea truth. Wea. Forecasting, 5, 596-612.

Sanders, F., L. F. Bosart, and C.-C. Lai, 1991: Initiation and evolution of an intense upper-level front. Mon. Wea. Rev., 119, 1337-1367.

Saucier, W. J., 1955: Principles of Meteorological Analysis. University of Chicago Press, 438 pp.

Schultz, D. M., and C. F. Mass, 1993: The occlusion process in a midlatitude cyclone over land. Mon. Wea. Rev., 121, 918-940.

Shapiro, M. A., 1970: On the applicability of the geostrophic approximation to upper-level frontal-scale motions. J. Atmos. Sci., 27, 408-420.

Shapiro, M. A., 1976: The role of turbulent heat flux in the generation of potential vorticity in the vicinity of upper-level jet stream systems. Mon. Wea. Rev., 104, 892-906.

Shapiro, M. A., 1981: Frontogenesis and geostrophically forced secondary circulations in the vicinity of jet stream-frontal zone systems. J. Atmos. Sci., 38, 954-973.

Shapiro, M. A., 1982: Mesoscale weather systems of the central United States. CIRES/NOAA Tech. Rep., University of Colorado, 78 pp. [Available from Cooperative Institute for Research in Environmental Sciences, University of Colorado/NOAA, Boulder, CO 80309 USA.]

Shapiro, M. A., and D. Keyser, 1990: Fronts, jet streams and the tropopause. Extratropical Cyclones, The Erik Palmén Memorial Volume, C. W. Newton and E. O. Holopainen, Eds., Amer. Meteor. Soc., 167-191.

Sinclair, M. R., and M. J. Revell, 1999: Classification and composite diagnosis of extratropical cyclogenesis in the southwest Pacific. Mon. Wea. Rev., 127, in press.

Staley, D. O., 1960: Evaluation of potential-vorticity changes near the tropopause and the related vertical motions, vertical advection of vorticity, and transfer of radioactive debris from stratosphere to troposphere. J. Meteor., 17, 591-620.

Stauffer, D. R., and N. L. Seaman, 1990: Use of four-dimensional data assimilation in a limited-area mesoscale model. Part I: Experiments with synoptic-scale data. Mon. Wea. Rev., 118, 1250-1277.

Uccellini, L. W., 1990: Processes contributing to the rapid development of extratropical cyclones. Extratropical Cyclones, The Erik Palmén Memorial Volume, C. W. Newton and E. O. Holopainen, Eds., Amer. Meteor. Soc., 81-105.

Uccellini, L. W., D. Keyser, K. F. Brill, and C. H. Wash, 1985: The Presidents' Day Cyclone of 18-19 February 1979: Influence of upstream trough amplification and associated tropopause folding on rapid cyclogenesis. Mon. Wea. Rev., 113, 962-988.

Zhang, D., and R. A. Anthes, 1982: A high-resolution model of the planetary boundary layer-Sensitivity tests and comparisons with SESAME-79 data. J. Appl. Meteor., 21, 1594-1609.

  Tables

Table 1: Cases of upper-level fronts that deviate from the Shapiro conceptual model with strong cold advection in the base of the trough and/or warm advection in the trailing ridge.
Palmén and Nagler (1949) Figs. 1 and 2
Reed and Sanders (1953) Figs. 1 and 2
Newton (1954) Figs. 3 and 4
Reed (1955) Figs. 5 and 6
Newton (1958) Figs. 2 and 4
Staley (1960) Figs. 2 and 3
Newton and Palmén (1963) Fig. 3
Shapiro (1976) Figs. 15 and 16
Sanders (1988) Figs. 9 and 10
Neiman and Shapiro (1989) Figs. 4, 9, and 12
Sanders et al. (1991) pp. 1364-1365
Djuri\'c (1994) Fig. 8-18
Lackmann et al. (1997) Figs. 12(a),(c)
Pyle (1997) p. 64, Fig. 23(d)

Table 2: Climatology of upper-level fronts associated with cyclones making landfall on the west coast of North America between 35°N and 60°N in December, January, and February for the six winters 1988-89 through 1993-94. Percentages may not add up to 100% due to round-off of decimal values.
SIX-YEAR YEARLY PERCENT
CATEGORY TOTAL RANGE OF TOTAL
number of cyclones 149 19-30 100
southwesterly flow 66 6-16 44
northwesterly flow 21 1-6 14
zonal flow 18 0-6 12
other flow 44 3-14 30
equivalent barotropic to cold advection 35 1-10 23
decreasing cold advection 19 0-6 13
weak advection or equivalent barotropic 73 7-21 49
warm advection 5 0-2 3
combination 17 1-6 11

Table 3: Climatology of equivalent-barotropic stage to cold-advection stage events from Table 2.
SIX-YEAR YEARLY PERCENT
CATEGORY TOTAL RANGE OF TOTAL
number of cyclones:
equivalent barotropic to cold advection 35 1-10 100
southwesterly flow 22 1-8 63
northwesterly flow 6 0-3 17
zonal flow 5 0-2 14
other flow 2 0-1 6

  Table Captions

Table 1: Cases of upper-level fronts that deviate from the Shapiro conceptual model with strong cold advection in the base of the trough and/or warm advection in the trailing ridge.

Table 2: Climatology of upper-level fronts associated with cyclones making landfall on the west coast of North America between 35°N and 60°N in December, January, and February for the six winters 1988-89 through 1993-94. Percentages may not add up to 100% due to round-off of decimal values.

Table 3: Climatology of equivalent-barotropic stage to cold-advection stage events from Table 2.

  Figure Captions

Figure 1: Idealized schematic depiction on an upper-tropospheric isobaric surface of the evolution of an upper-level jet-front system through a midlatitude baroclinic wave over a 72-h period: (a) formation of jet-front in the confluence between the mid- and high-latitude currents; (b) jet-front situated in the northwesterly flow inflection of amplifying wave; (c) jet-front at the base of the trough of fully developed wave; (d) jet-front situated in the southwesterly flow inflection of damping wave. Geopotential height contours, thick solid lines; isotachs, thick dashed lines; isentropes or isotherms, thin dashed lines. From Shapiro (1982, Fig. 7).

Figure 2: Schematic illustrating vertical tilting terms in the expression for vector frontogenesis of the horizontal potential temperature gradient. Solid lines are isentropes at two times (t = 0 and t = dt). Vectors are labeled. Insets: Solid lines are isentropes at t = 0 and dashed lines are isentropes at t = dt. Open arrows represent vertical motions. (a) Tilting scalar frontogenesis: ÑH w and ÑH q in the same direction such that n3 is maximized and ÑHw×ÑH q = 0. F is therefore in the direction of ÑH q. (b) Tilting rotational frontogenesis: ÑH w and ÑH q perpendicular such that |s3| is maximized and ÑHw·ÑH q = 0. F is therefore perpendicular to ÑH q.

Figure 3: SW 500-hPa geopotential height (thin dashed lines every 6 dam), relative vorticity of total horizontal wind (10-5 s-1-1, shaded according to scale at bottom of figure), potential temperature (thin solid lines every 2 K), and potential-temperature advection by the total horizontal wind [thick solid (dashed) lines every 12  ×10-5 K s-1, from -36 to 36 ×10-5 K s-1 represent warm (cold) advection]: (a) 1200 UTC 12 December 1988; (b) 0000 UTC 13 December 1988.

Figure 4: GOES-6 infrared satellite images; arrows point to cloud features described in text: (a) 0016 UTC 12 December 1988, (b) 1216 UTC 12 December 1988, (c) 2316 UTC 12 December 1988.

Figure 5: As in Fig. 3 except for NW: (a) 0000 UTC 13 December 1988; (b) 1200 UTC 13 December 1988.

Figure 6: SW 500-hPa omega (pressure vertical velocity) (thick lines every 2 mb s-1, from -8 to 8 mb s-1; dashed (solid) lines represent ascent (descent)), potential temperature (thin solid lines every 2 K), and relative vorticity of total horizontal wind (10-5 s-1-1, shaded according to scale at bottom of figure): (a) 1200 UTC 12 December 1988; (b) 0000 UTC 13 December 1988.

Figure 7: As in Fig. 6 except for NW: (a) 0000 UTC 13 December 1988; (b) 1200 UTC 13 December 1988.

Figure 8: SW 500-hPa relative vorticity of total horizontal wind (10-5 s-1-1, shaded according to scale at bottom of figure), and potential temperature (thin solid lines every 2 K) at 1200 UTC 12 December 1988: (a) -Fn divergence term (every 1 ×10-10 K m-1 s-1, from -6 to 6 ×10-10 K m-1 s-1); (b) -Fn deformation term (every 3×10-10 K m-1 s-1, from -18 to 18 ×10-10 K m-1 s-1) and axes of dilatation of total horizontal wind (10-5 s-1, scaled according to legend; separation between displayed axes of dilatation is 180 km (every other grid point)); (c) -Fn vertical tilting term (every 3 ×10-10 K m-1 s-1, from -18 to 18 ×10-10 K m-1 s-1); (d) total -Fn (every 3 ×10-10 K m-1 s-1, from -18 to 18×10-10 K m-1 s-1).

Figure 9: As in Fig. 8 except for 0000 UTC 13 December 1988.

Figure 10: As in Fig. 8 except for NW at 0000 UTC 13 December 1988.

Figure 11: As in Fig. 3 except for NW at 1200 UTC 13 December 1988.

Figure 12: SW 500-hPa relative vorticity of total horizontal wind (10-5 s-1-1, shaded according to scale at bottom of figure), and potential temperature (thin solid lines every 2 K) at 1200 UTC 12 December 1988: (a) Fs vorticity term (every 3 ×10-10 K m-1 s-1, from -18 to 18 ×10-10 K m-1 s-1); (b) Fs deformation term (every 3×10-10 K m-1 s-1, from -18 to 18 ×10-10 K m-1 s-1) and axes of dilatation of total horizontal wind (10-5 s-1, scaled according to legend; separation between displayed axes of dilatation is 180 km (every other grid point)); (c) Fs vertical tilting term (every 3 ×10-10 K m-1 s-1, from -18 to 18 ×10-10 K m-1 s-1); (d) total Fs (every 3 ×10-10 K m-1 s-1, from -18 to 18×10-10 K m-1 s-1).

Figure 13: As in Fig. 12 except for 0000 UTC 13 December 1988.

Figure 14: As in Fig. 12 except for NW at 0000 UTC 13 December 1988.

Figure 15: As in Fig. 12 except for NW at 1200 UTC 13 December 1988.

Figure 16: Revised conceptual model: Idealized schematic depiction on an upper-tropospheric isobaric surface of the early evolution of an upper-level jet-front system through a midlatitude baroclinic wave over a 12-24-h period: (a) southwesterly flow case; (b) northwesterly flow case. Geopotential height contours (solid grey lines) are labeled Z in (a)(i), isentropes (solid black lines) are labeled q in (a)(i), and relative vorticity (shaded) is labeled z in (a)(i).


Footnotes:

1 A similar expression for ``tilting scalar frontogenesis'' in an x-z plane, involving vag, the ageostrophic cross-front wind, can be found in Keyser et al. (1986, (2.12)).

2 Strictly speaking, the rotational frontogenesis does not describe the time rate of change of the angle between isohypses and isotherms, as the height field, as well as the thermal field, will be altered by the relative vorticity, deformation, and tilting terms. Rigorous resolution of this discrepancy, however, would involve inversion techniques, the methodology for which has not matured. Nevertheless, our results suggest that the rotation of the isohypses by the vorticity is a secondary effect, at least in the cases presented here.

3 We present -Fn rather than Fn in order that regions of positive -Fn represent regions of positive frontogenesis, according to (5a).


This page created on 12/22/98 by David Schultz.


Return to Current Research Interests.
Return to David Schultz's homepage.