UBC Faculty Research and Publications

Parameterization of Joint Frequency Distributions of Potential Temperature and Water Vapor Mixing Ratio… Berg, Larry K.; Stull, Roland B. Apr 30, 2004

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata


52383-Stull_AMS_2004_JAS813pdf.pdf [ 587.39kB ]
JSON: 52383-1.0041849.json
JSON-LD: 52383-1.0041849-ld.json
RDF/XML (Pretty): 52383-1.0041849-rdf.xml
RDF/JSON: 52383-1.0041849-rdf.json
Turtle: 52383-1.0041849-turtle.txt
N-Triples: 52383-1.0041849-rdf-ntriples.txt
Original Record: 52383-1.0041849-source.json
Full Text

Full Text

1 APRIL 2004  BERG AND STULL  813  Parameterization of Joint Frequency Distributions of Potential Temperature and Water Vapor Mixing Ratio in the Daytime Convective Boundary Layer LARRY K. BERG Pacific Northwest National Laboratory, Richland, Washington  ROLAND B. STULL Atmospheric Science Programme, Department of Earth and Ocean Sciences, The University of British Columbia, Vancouver, British Columbia, Canada (Manuscript received 16 September 2002, in final form 3 November 2003) ABSTRACT Joint frequency distributions (JFDs) of potential temperature ( ␪) versus water vapor mixing ratio (r) within the convective boundary layer were measured during a new field experiment named Boundary Layer Experiment 1996 (BLX96). These JFDs were found to be tilted, with the tilt a function of both height and boundary layer dynamics. These distributions are also skewed and more peaked than a joint Gaussian distribution. Three different methods are used to generate joint probability density functions (JPDFs) that approximate observed JFDs. Two classical methods, one based on a Gaussian fit and another based on surface-layer processes, are reviewed. A new method is devised, which treats the observed JFD as a mixing diagram. In the absence of advection, the only source regions for air in the mixing diagram are the surface and the entrainment zone. Thus, the tilt of the JFD can be explained by various mixtures from these two source regions. Methods that can be used to parameterize the mixing JPDF are presented. The primary advantage of this method is that the tilt is determined explicitly from properties of the surface, mixed layer, and entrainment zone. Similarity methods are used to parameterize all variables needed by the Gaussian model. The Bowen ratio and the total energy input are used to parameterize the tilt of the surface energy budget JPDF, while similarity methods are used to define the spread of the JPDF along the two axes. Relationships between the surface and mixed layer, and the mixed layer and free atmosphere are used to tilt the mixing diagram JPDF, while similarity methods are used to estimate the spread of the JPDF. The parameterizations are developed using a ‘‘calibration’’ subset of data acquired during BLX96. A ‘‘verification’’ subset of data, also acquired during BLX96, is used to show that the parameterized mixing diagram method is superior to the other two methods, because it has either a smaller error or is less sensitive to the value of the correlation between ␪ and r.  1. Introduction Joint frequency distributions (JFDs) provide an alternative method to investigate physical processes in the convective boundary layer. In addition, discrete JFDs provide data reduction of long time series by combining many observations into finite-sized bins. Holland (1973) used JFDs to investigate wave shapes. Mahrt and Paumier (1982) used JFDs of specific humidity and potential temperature (␪) to look at the effects of entrainment and mixing in a stratocumulus-topped boundary layer. Grossman (1984) used JFDs of vertical velocity (w) and water vapor density to perform conditional sampling of data measured in a maritime boundary layer. Using data from the Air Mass Transformation Experiment, Mahrt Corresponding author address: Dr. Larry K. Berg, Pacific Northwest National Laboratory, P.O. Box 999 MSIN K9-30, Richland, WA 99352. E-mail: Larry.Berg@pnl.gov  ᭧ 2004 American Meteorological Society  and Paumier (1984) constructed JFDs of w and ␪. Deardorff and Willis (1985) examined JFDs of w and temperature generated from their laboratory tank model. They also used JFDs to perform conditional sampling, looking for coherent structures in the boundary layer. Other authors (Lewellen and Yoh 1993; Xu and Randall 1996; Wang and Stevens 2000) have used finescale models to generate JFDs for cloudy boundary layers. Larson et al. (2001) examined distributions of conserved variables in cloudy boundary layers. Berg (2002) developed a cumulus parameterization that used subcloud distributions of ␪ versus r to predict cloud cover and cloud size distributions. A number of authors have examined the general shape of JFDs in the boundary layer. Lenschow (1970) suggested that observations of ␪ in the boundary layer are normally distributed, while Manton (1975) suggested that distributions of temperature are not normally distributed. Later work (e.g., Lenschow and Stephens 1980; Deardorff and Willis 1985; Cuijpers and Holtslag 1998;  814  JOURNAL OF THE ATMOSPHERIC SCIENCES  Wang and Stevens 2000) supported the contention by Manton that distributions of variables in the boundary layer are not normally distributed. While JFDs have been used as tools to examine the convective boundary layer, there have been few attempts to quantify the shape of JFDs. There are two strategies that can be used to account for the tilt of the JFDs. The cross correlation between the two variables can be used to account for the tilt, or a coordinate transform can be used to describe the tilt of the JFD and eliminate the cross correlation. Wyngaard and Moeng (1992, hereafter WM92) suggested a statistical approach, using the correlation between w and any arbitrary scalar to describe the tilt of the JFDs of w and that scalar. They compared their theoretical joint probability density functions (JPDFs) to JFDs generated by a large eddy simulation (LES). Schrieber et al. (1996, hereafter SSZ) examined distributions of virtual potential temperature (␪ ␷ ) and moisture [their moisture variable was the height of the lifting condensation level (zLCL )]. They took a different tack than WM92, developing a physically based method to perform a coordinate transformation to tilt the JPDFs and to eliminate the cross correlation observed in the surface layer. A new conceptual model that can describe the tilt of JFDs of ␪ versus water vapor mixing ratio (r) is proposed in this paper. Like the work of SSZ, this method is based on physical processes in the boundary layer. This approach views mixed-layer turbulence as a mixing process, and treats the JFD of ␪ versus r as a mixing diagram (Paluch 1979; Hanson 1984; Betts 1984, 1985; Betts and Albrecht 1987). When advection is small, the mixed-layer distribution is a mixture of air from two different source regions, the surface layer and the free atmosphere. Thus, the tilt of the JFD can be determined if the thermodynamic properties of the surface, mixed layer, and entrainment zone are known. In this work,  f (␪, r) ϭ  Ά  the properties of the mixed layer, the surface, and the free atmosphere are used to define a coordinate transformation to describe the tilt of Gaussian JPDFs. Data from Boundary Layer Experiment 1996 (BLX96) will be used to demonstrate that the new methods can successfully describe the size and shape of the JFDs in the convective mixed layer. First, some definitions associated with the JFDs and JPDFs will be presented. Some relevant aspects of BLX96 will be reported in section 3. Observed JFDs measured during BLX96 are discussed next. The methods proposed by WM92 and SSZ and the new method are described in section 5. The methods used to determine the parameters derived from similarity theory are presented in section 6. Finally, in section 7, each of the three methods will be compared to a subset of BLX96 case study days. 2. Definitions Joint frequency distributions count the occurrences of observations, ␪ˆ and rˆ, that lie within an interval ⌬ ␪ and ⌬r: F(␪, r, ⌬␪, ⌬r) ϭ count(␪ Ͻ ␪ˆ Ͻ ␪ ϩ ⌬␪, r Ͻ rˆ Ͻ r ϩ ⌬r).  3. Boundary Layer Experiment 1996 BLX96 was conducted between 15 July and 13 August 1996 over regions within the U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) Clouds and Radiation Test-  (1)  Since every observation, ␪ˆ and rˆ, is somewhere within the whole ␪ – r ‘‘space,’’ if there are N observations, then 1 N  ͸  F(␪, r, ⌬␪, ⌬r)⌬␪⌬r ϵ 1.  (2)  all observations  Sometimes 1/N is included with F to make it a relative JFD, and sometimes F is written as a frequency density that must be integrated over ␪ˆ and rˆ to get the total count. If each of the variables has a normal distribution, then their joint normal distribution (written here as a JPDF) can be expressed as  [΂  1 Ϫ1 ␪Ϫ␪ exp 2 Ϫ1/ 2 2 2␲␴␪ ␴ r (1 Ϫ ␳␪,r ) 2(1 Ϫ ␳␪,r ) ␴␪  where ␳ ␪,r is the correlation between ␪ and r, ␴ ␪ is the standard deviation of ␪, ␴ r is the standard deviation of r, and the overbar represents a mean quantity. As indicated by (3), the interplay of the three parameters: ␳ ␪,r , ␴ ␪ , and ␴ r define the tilt of, and spread along, the axes of the JPDF.  VOLUME 61  ΃  2  Ϫ  ΂  2␳␪,r (␪ Ϫ ␪ )(r Ϫ r) rϪr ϩ ␴␪ ␴ r ␴r  ΃ ]· , 2  (3)  bed (CART; Stokes and Schwartz 1994). This region includes parts of north-central Oklahoma and south-central Kansas. Details of BLX96 were reported by Stull et al. (1997) and Berg et al. (1997) and can be found online (see http://weather.eos.ubc.ca). A total of 12 research flights were flown during BLX96 using the University of Wyoming King Air aircraft. Each flight started between 1000 and 1030 LST, and lasted 4 to 4.5 h. Only 10 flights were used for this study: 15, 16, 22, 23, 25, 27, 28, and 31 July and 2 and 13 August. The flight on 21 July was excluded because the tube leading to the fast response humidity sensor was not connected properly during the flight, 4 August was excluded because  1 APRIL 2004  BERG AND STULL  wind speeds (about 8 m s Ϫ1 ) led to forced, rather than free convection (Obukhov length of Ϫ257.8 m). On 16 July, the second half of the pattern was excluded because the measured moisture flux was suspect during the last near-surface leg (Berg 2002). Three different locations within ARM SGP CART were chosen for BLX96 flights. The three legs were chosen to be over regions of different land use and to be close to at least one ARM SGP CART surface station, and are named after a nearby hamlet or town. The Lamont track was over flat terrain devoted to both wheat farming and hay production. Prior to BLX96 the wheat had been cut, and the wheat fields were bare. The ARM SGP CART central facility was located at the east end of this track. The Meeker track had many rolling hills and much forest coverage. The forest was not dense, and most of the trees were less than 10 m high. Some of the land was devoted to agriculture, either crops or hay production. The Winfield leg was over a region with some rolling hills. Most of the leg was over pasture, but there was some sparse forest near the east end of the track. A unique flight pattern was designed for BLX96 to simultaneously satisfy three different experiments: a radix layer experiment (Santoso and Stull 1998, Santoso and Stull 2001), a convective transport theory experiment (Stull 1994), and work related to the JFDs and boundary layer cloud cover (Berg 2002). In each case, the pattern was oriented approximately perpendicular to the mean boundary layer flow. For all of the BLX96 flights, this resulted in tracks aligned in an east–west direction. Soundings were flown from near the surface to above cloud top 3 times during each flight; near the start, middle, and end of the pattern. Level horizontal legs about 72 km long were flown at a range of heights within the boundary layer. The highest horizontal leg was flown approximately 1000 m above ground level (AGL) or just below cloud base, whichever was lower. Other horizontal legs were flown at approximately 750, 500, and 250 m AGL. During each flight, three approximately terrain following near-surface legs were flown 30 to 60 m above the surface. Fluxes measured from these legs were used to compute standard boundary layer scales. The heights of all legs correspond to heights ranging from about 0.02z i to 0.8z i , where z i is average mixed-layer depth. A discrete Fourier transform (DFT) was used to filter the data collected from horizontal legs. A DFT was used, rather than a fast Fourier transform (FFT), to eliminate the need for zero padding or truncation of the data from each flight leg. Wavelengths less than approximately 20 m were removed from the data, primarily because the instruments on the aircraft were not collocated. Wavelengths greater than 5 km were removed to eliminate mesoscale effects that might be present, but that would have been poorly sampled. Variances and turbulent fluxes were computed using data with different long wavelength cutoffs, and the contribution by scales greater  815  FIG. 1. Variation with dimensionless height of the correlation of ␪ and r for BLX96 flight legs (crosses) and observations from Wyngaard et al. (1978; triangles) over the ocean.  than 5 km was generally small. DFTs or FFTs require that all data points be equally spaced and that there are no missing data. Replacing bad data is a challenge. In our case, the points identified as ‘‘bad’’ were not noticeably different than their neighbors; so they were not removed for the DFT filtering. However, these points were excluded for all subsequent calculations, for example calculations to determine the eddy-correlation flux or JFDs. To construct JFDs, each observation was placed into a ␪ versus r bin 0.1 K by 0.1 g kg Ϫ1 wide. 4. Observations While SSZ were the first to publish JFDs of ␪ ␷ versus zLCL , their data were limited to observations made approximately 100 m above the surface (roughly 0.05z i to 0.10z i ), and their flight legs averaged roughly 12 km in length. During BLX96, JFDs of ␪ versus r were observed at heights ranging from 0.02z i to just above 0.8z i , with most flight legs below 0.6z i . Each BLX96 flight leg was approximately 70 km long. SSZ found a number of attributes were common to all JFDs; namely, they are bounded, they are sharply peaked, the spread is skewed along any axis of the distribution, they have a tilted axes, and a few JFDs have multiple peaks. The JFDs that were observed aloft during BLX96 share the key attributes described by SSZ for JFDs near the surface. However, there are some important differences. All the JFDs observed below 0.1z i have a positive correlation between ␪ and r (␳ ␪,r Ͼ 0), in other words, there is a positive slope to the JFD (Fig. 1). There is also a large range in the observed ␳ ␪,r . Between 0.2z i and 0.4z i there is a transition zone where some, but not all the  816  JOURNAL OF THE ATMOSPHERIC SCIENCES  VOLUME 61  FIG. 2. Sequence of thee JFDs of ␪ vs r at three different heights: (a) 0.03 z i , (b) 0.30 z i and (c) 0.74 z i measured on 16 Jul 1996. The first contour indicates a normalized frequency of 0.005, and the contour interval is 0.01.  JFDs have a negative ␳ ␪,r , or in other words, a negative tilt. Above 0.4z i , all of the observed JFDs have a negative ␳ ␪,r . The change of ␳ ␪,r , or the change of tilt of the JFDs with height can be explained. Solar energy is absorbed at the earth’s surface. Surface-layer processes cause air near the surface to warm and become moister. As these warm and moist parcels rise above the surface they cause the observed positive ␳ ␪,r . Near z i , entrainment brings relatively warm and dry parcels from the free atmosphere down into the mixed layer. These warm and dry parcels lead to the observed negative ␳ ␪,r . Any method used to describe JFDs through the depth of the boundary layer must account for this behavior. These measured values of ␳ ␪,r can be compared to those reported by Wyngaard et al. (1978; Fig. 1). It appears that ␳ ␪,r measured by Wyngaard et al. (1978) are larger near the surface than those measured during BLX96. Their measurements were from over the ocean where vigorous convection was driven by cold air flowing over warm water. In this scenario ␳ ␪,r would be quite large because the surface layer is very moist. Over the middle of a continent, one would expect some of the thermals to be relatively dry and a smaller value of ␳ ␪,r is expected. The observed differences between the results of Wyngaard et al. (1978) and BLX96 highlights that ␳ ␪,r is not universal, but changes from time to time and place to place. A sequence of JFDs observed on 16 July 1996 highlight this change of tilt (Fig. 2). Each JFD was separated in time by about 15 min and each leg took about 15 min to fly. No corrections for nonstationarity of the boundary layer between flight legs have been applied. This example is typical of JFDs measured at various heights during BLX96, the positive slope for the leg near the surface is obvious, as is the negative slope near the mixed-layer top. In this case, the middle leg seems to still have a positive tilt, but there is some evidence  of a tail reaching toward warmer and dryer conditions. All of the JFDs observed during BLX96 were skewed. The average skewness of ␪ was 0.75, while the average skewness of r was Ϫ0.17. The skewness was also a function of height. Each distribution was also more peaked than a Gaussian distribution. Figure 2 suggests that the lower-altitude JFDs might also have a smaller kurtosis1 than the upper-level legs. However, this trend disappears when all BLX96 case study days are considered. a. Identification of coherent drafts Coherent drafts will be used for the new method for describing the JFDs. A number of other researchers have investigated the convective boundary layer by locating coherent structures (e.g., Manton 1975; Coulman 1978; Lenschow and Stephens 1980; Nicholls and LeMone 1980; Greenhut and Khalsa 1982, 1987; Grossman 1984; Khalsa and Greenhut 1985; Young 1988; Williams and Hacker 1992; Schumann and Moeng 1991; Siebesma and Cuijpers 1995; Wang and Stevens 2000). Three different criteria are needed to identify the drafts: an indicator variable to define the thermal, some threshold value of the indicator variable, and some length scale. Unfortunately, all of these criteria are arbitrary. Either ␪ ␷ , T, r, or w have been used by various investigators as the indicator variable (Table 1). Crum and Stull (1987) tried to determine which indicator variable works the best. Based on observations of thermals made using a research aircraft and coincident lidar observations over land and near the mixed-layer top, they suggested that r is the best choice to identify thermals. In 1 Kurtosis has been defined relative to a Gaussian distribution (e.g., Press et al. 1992), so that a Guassian distribution would have a kurtosis of 0.  1 APRIL 2004  817  BERG AND STULL  TABLE 1. List of indicator variables used in different field studies, the physical reasoning behing the choice of the variable, some potential shortfalls associated with that variable, and the researchers that have used that variable. Indicator variable  Physical reasoning  w  Direct measure of thermal velocity  Noisy signal  Potential shortfalls  ␪␷  Buoyancy is the dominant force driving the thermal  r  Thermals tend to be more moist than the environment  Other forces and inertia can be important, especially high in the mixed layer where the thermal is often less buoyant than its surroundings May not always be true over land; r is only weakly related to buoyancy  T  T contributes substantially to thermal buoyancy  r also contributes to thermal buoyancy, other forces can be important (see discussion for ␪␷ )  Researchers Nicholls and LeMone (1980); Greenhut and Khalsa (1982, 1987); Grossman (1984); Khalsa and Greenhut (1985); Young (1988); Siebesma and Cuijpers (1995); Wang and Stevens (2000) Williams and Hacker (1992)  Lenshow and Stephens (1980); Nicholls and LeMone (1980); Crum and Stull (1987); Siebesma and Cuijpers* (1995); Wang and Stevens* (2000) Manton (1977); Coulman (1978)  * Siebesma and Cuijpers (1995) and Wang and Stevens (2000) used the liquid water content to identify updrafts, rather than mixing ratio.  contrast, Schumann and Moeng (1991) suggested that using w to locate the thermals results in the most consistent results. Siebesma and Cuijpers (1995) used both r and w to define thermals in the cloudy boundary layer. In this study, w will be used as an indicator variable because there was less scatter in the thermal statistics when using w rather than r. One common way to define the threshold is ␺Ϯthresh ϭ ͙mthresh ␴ Ϯ2 ␺, where mthresh is called the threshold multiplier (dimensionless) and ␺ is the indicator variable. The value of ␴ Ϯ2 ␺ is the variance of ␺ calculated using only the observations of ␺ greater than the mean value of ␺ (i.e., ␴ ϩ2 ␺), and the value of ␴ Ϫ2 ␺ is the variance of ␺ calculated using only the observations of ␺ less than the mean value of ␺ (i.e., ␴ Ϫ2 ␺). Following the lead of Greenhut and Khalsa (1982), the threshold value in this 2 study was defined as wϮthresh ϭ ͙␴ Ϯw /2 . One advantage of this choice of wϮthresh is that the upward and downward mass flux are approximately the same (Greenhut and Khalsa 1982). To be called a coherent draft, a structure must be horizontally wider than a specified length scale. In addition, if a distance less than the length scale separates two drafts, they are assumed to be one structure. Coulman (1978) used a length scale of zero, his only criterion was the T threshold. Some authors (Greenhut and Khalsa 1982, 1987; Khalsa and Greenhut 1985) used a constant length scale. Young (1988) used a length scale of 0.1z i . He chose this value because, in general, the smaller turbulence scales are within the inertial cascade. In our work, we will use a slightly smaller length scale of 0.05z i . A number of BLX96 spectra, particularly near the surface, have a peak very close to 0.1z i , and drop off at shorter scales (H. Modzelewski 2002, personal communication). Therefore, choosing a cutoff of 0.1z i might inappropriately remove some of the coherent drafts.  b. Mixing diagrams JFDs of ␪ versus r can be viewed as a mixing diagram because both ␪ and r are conserved in the dry boundary layer. Mixing diagrams are useful because mixtures of air from two source regions lie on a straight line connecting the conservative properties of the two regions (Paluch 1979; Hanson 1984; Betts 1982, 1984, 1985; Betts and Albrecht 1987). As an example, mixing diagrams of data collected from the slant ascent/descent soundings flown during BLX96 can be examined. Figure 3 shows three sequential soundings measured on 27 July. A similar plot was presented by Betts (1992). This plot is typical of mixing diagrams generated from the BLX96 soundings. Three different families of points are highlighted: mixed-layer points, entrainment-zone points and freeatmosphere points. On this particular day, ␪ and r of the free atmosphere are found in the lower-right quadrant of the diagram. For this free-atmosphere air, the different ␪ values correspond to quasi-stationary conditions at different heights above ground. The mixed layer warmed with time (as shown by the shift to the right of the mixed-layer points in the top-left quadrant), and dried throughout the day (the mixed-layer points shifted down). The entrainment zone was determined from the slant ascent/descent soundings by locating the heights at which there were large jumps in the value of r. These large changes in r are caused by dry freeatmosphere air intruding down into the moist convective boundary layer. This type of identification is relatively easy with aircraft slant soundings, which pass through a large horizontal distance (approximately 70 km during BLX96) and thus intersect many thermal circulations. It would be more difficult with balloon soundings that may drift along with a single thermal, and thus give effectively only a point sample. The points in the entrainment zone should lie on a mixing line between the  818  JOURNAL OF THE ATMOSPHERIC SCIENCES  VOLUME 61  FIG. 3. Mixing diagrams for the three sequential soundings measured on 27 Jul. Each symbol represents a block average over a depth of 20 m. The open symbols correspond to the mixed-layer values, circles with crosses, squares with crosses, and asterisks correspond to the entrainment-zone values, and the filled symbols correspond to the free atmosphere. Shading highlights observations that are part of the entrainment zone.  free-atmosphere values just above z i , and the mixedlayer mean. That was the case for nearly all of the BLX96 soundings. Surface-layer points are not highlighted in Fig. 3. For flight safety reasons, the soundings flown during BLX96 did not go low enough to capture much of the superadiabatic surface layer, so surface-layer properties were inferred from the updraft properties observed during any given leg in the following way. We assume that the mean properties of the observed updrafts lie on the mixing line that connects the mixed-layer point and the surface point on the mixing diagram. Therefore, these updraft properties can serve as a proxy for the surface-layer properties. The mean thermodynamic properties of the updrafts, downdrafts, and environment can also be examined using a mixing diagram. An example, measured on 25 July, is plotted in Fig. 4. In this example, the mean updraft, downdraft, and environment air at any one time lie on a straight line or mixing line in ␪ versus r space. This occurred for nearly all of the flight legs during BLX96. The only exceptions occurred on a few legs near the top of the mixed layer. If the mean thermodynamic properties of the updrafts, downdrafts and environment lie on a single mixing line, then (␪ u Ϫ ␪ )/ (r u Ϫ r ) ഠ (␪ u Ϫ ␪ d )/(r u Ϫ r d ), where u indicates an updraft mean, d indicates a downdraft mean, the overbar indicates a mixed-layer average. This information can  be used in a new parameterization for JPDFs, as described in the next section. 5. JPDF parameterizations a. Classical methods Two methods have been suggested in the literature to parameterize JFDs in the boundary layer. WM92 fit joint Gaussian JPDFs to JFDs of an arbitrary scalar and w generated by a LES model. They used three parameters to describe their JPDFs: the measured correlation between the scalar and w, the standard deviation of the scalar, and the standard deviation of w. Their method is statistical and did not offer a physical explanation for the tilt of the JFDs. While they constructed JPDFs of a scalar versus w, their mathematical techniques are also valid for JPDFs of ␪ versus r. Throughout the rest of this paper, this method will be called the statistical method. Three variables are needed to generate JPDFs using this method: ␴ ␪ , ␴ r , and ␳ ␪,r . SSZ introduced a physically based method to fit joint Gaussian JPDFs to JFDs of ␪ versus r measured in the surface layer. They attempted to fit JPDFs by introducing two new, rotated axes to remove the correlation between ␪ and r. They related the slope of one axis to the Bowen ratio (B), and the slope of the other axis to a variable they called the solar forcing (␪ F ). While these  1 APRIL 2004  BERG AND STULL  819  FIG. 4. Plot showing updraft (upward-pointing triangle), downdraft (downward-pointing triangle), and environment (circle) values of ␪ vs r for six horizontal flight legs on 25 Jul. The error bars represent the standard deviations. Note the different ␪ scales between used for the left-hand and right-hand sides.  new axes are not geometrically orthogonal in ␪ versus r space, changes along the B axis are independent of changes along ␪ F axes. This method will be called the B–␪ F method. Four variables are needed to generate JPDFs using this method, the slopes of the B and ␪ F axes, the standard deviation along the B axis (␴ B ), and the standard deviation along the ␪ F axis (␴␪ F). In the appendix, SSZ’s equations are converted from their ␪ ␷ versus zLCL framework to the ␪ versus r framework as used here. b. Mixing diagram method In the absence of advection, JFDs computed from observations made in the convective boundary layer are  composed of a large number of mixtures of air originating at the surface, air that has been entrained from the free atmosphere, and air that has resided in the boundary layer for some time. Thus, the JFD should stretch along a mixing line toward the thermodynamic properties of both the surface and the entrainment zone. If advection is also considered, then the mixed-layer mixture is composed of air with three different origins (Fig. 5), which tend to ‘‘pull’’ the mixture off of a straight line between surface and entrainment-zone air. The primary advantage of this method is its basis on physical processes in the mixed layer. The primary disadvantage, compared to the statistical method, is the requirement of an additional parameter to compute the fit.  820  JOURNAL OF THE ATMOSPHERIC SCIENCES  FIG. 5. Mixing diagram for ␪ vs r. The large circle represents the mixed-layer mean, the small circles represent the surface (Sfc) and entrainment-zone properties (EZ), and the square represents the advection (Adv). The large arrows indicate the ‘‘direction’’ that the different processes pull on the mixed-layer mean.  A coordinate transform is used to fit the mixing-diagram JPDFs to the JFDs and remove the correlation between ␪ and r. SSZ defined a transform to convert from ␪ and r to a coordinate system defined using the slopes of their B and ␪ F axes. The transform is used to create a JPDF with two different axes that are tilted relative to the ␪ and r axes. SSZ’s definitions can be used to convert ␪ and r data into a coordinate system, c and m, based on slope of the surface axis (Ssfc ) and the slope of the entrainment-zone axis (Sez ): c ϭ Sez  ΂S  ΃  ␪ Ϫ Ssfc r ez Ϫ Ssfc  and m ϭ  ␪ Ϫ Sez r , Ssfc Ϫ Sez  (4)  where c has units of temperature and m has units of mixing ratio. The values of c and m can then be projected from the surface and entrainment zone axes onto the ␪ and r axes. A joint Gaussian distribution, framed in c versus m coordinate system, takes the form f (m, c) ϭ  Ά [΂ ΃ ΂ ΃ ]·  1 1 m exp Ϫ 2␲␴ m ␴ c 2 ␴m  2  ϩ  c ␴c  2  ,  (5)  where the correlation between c and m has been ignored. Before integrating (5) over ␪ and r, (5) should be multiplied by the Jacobian determinate, which is defined as |J(c, m)| ϭ  ΂‫ץ ␪ץ‬r Ϫ ‫ץ ␪ץ‬r΃ . ‫ץ‬c ‫ץ‬m  ‫ץ‬m ‫ץ‬c  (6)  Using (4), it can be shown that | J(c, m) | ϭ SEZ (SEZ Ϫ Ssfc ) Ϫ1 . Multiplication by | J(c, m) | insures that the integral of the JPDF is equal to 1. To use (5) with c and m defined by (4), one must first determine the axes slopes. The surface axis slope Ssfc is defined as Ssfc ϭ (␪sfc Ϫ ␪ )/(rsfc Ϫ r),  (7)  where the overbar indicates a mean mixed-layer value  VOLUME 61  and the subscript sfc indicates a surface value. Equation (4) admits a special case (if Ssfc and Sez are approximately equal, causing divide by zero) in which the mixing diagram method will fail because the two axes coincide as a single axis. In model forecasts, Ssfc could be directly calculated using (7). An alternative expression for Ssfc , which will be used here to parameterize the JPDFs, is based on the properties of the updrafts and the mixedlayer mean such that, Ssfc ϭ (␪up Ϫ ␪ )/(rup Ϫ r ). It is tempting to define the slope of the entrainment zone axis in the same way as Ssfc . However, as shown in Fig. 3, the entrainment zone cannot be adequately described using a single point. Rather Sez is defined using slope of a linear fit to observations of ␪ and r made in the entrainment zone, which are highlighted in Fig. 3. In model forecasts, Sez could be inferred by the difference in ␪ and r across the mixed-layer top. An alternative expression for Ssfc , based on the Bowen ratio (B), can be derived (SSZ). When B is large, the sensible heat flux is large and the moisture flux is small; variations of ␪ tend to be large relative to variations of r, and Ssfc is large. When B is small, the heat flux is small relative to the latent heat flux, the variations of ␪ tend to be smaller than the variations of r, and Ssfc is small. To investigate this approach further, it is convenient to define a new B to be a function of height: Bz ϭ  C p wЈ␪Ј(z) , L␷ wЈrЈ(z)  (8)  where the subscript z indicates B is a function of height, and where C p is the specific heat of air and L ␷ is the latent heat of vaporization. Using the mass-flux approximation (e.g., Betts 1975, 1976; Manton 1975; Greenhut and Khalsa 1982, 1987; Nicholls and LeMone 1980; Randall et al. 1992; Siebesma and Cuijpers 1995) and ignoring the environmental flux, (8) can be written as Bz ഠ  C p (␦w u ␦␪ u ␣ u ϩ ␦w d ␦␪ d ␣ d ) , L␷ (␦w u ␦r u ␣ u ϩ ␦w d ␦r d ␣ d )  (9)  where ␦w is the average draft velocity, ␦␪ is the average draft perturbation value of ␪, ␦r is the average draft perturbation value of r, ␣ is the fraction of legs covered by drafts, the subscript u indicates updrafts, and the subscript d indicates downdrafts. Equation (9) can be simplified because of mass continuity in the boundary layer, ␦w u ␣ u ϭ ␦w d ␣ d , and rearranged to yield Ssfc : Ssfc ϭ  ␦␪ u Ϫ ␦␪ d LB ϭ ␷ z, ␦r u ϩ ␦r d Cp  (10)  which is similar to equations found by Betts (1984) and SSZ. Figure 6 shows the observed slopes calculated from the observed updraft and downdraft properties observed during BLX96 and the slopes predicted by (10) using the value of B z calculated from the horizontal flight legs  1 APRIL 2004  821  BERG AND STULL  TABLE 3. Calibration subset of BLX96 case-study days used to perform the statistical fits, and the associated location and cloudcover amount. Date 15 22 25 27 28  FIG. 6. Predicted surface axis slope using (10) (circles) vs observed surface axis slope for BLX96 case study days. The solid line is the 1:1 line.  Jul Jul Jul Jul Jul  Location Winfield Winfield Winfield Lamont Meeker  Cu coverage Scattered None Scattered Scattered Scattered  Cu Cu Cu Cu  the results. All of our observations of ␴ ␪ /␪* and ␴ r /r*, where ␪* is the Deardorff convective velocity scale and r* is a moisture scale, fall onto universal curves. Therefore, our results should not be sensitive to the choice of calibration and verification subsets. In addition, our research flights were conducted at three different locations each with different land-use characteristics, which makes our findings more robust. a. Statistical method  over a range of heights (Fig. 6). A further simplification can be made to (10). In section 4b, it was shown that, (␪ u Ϫ ␪ )/(r u Ϫ r ) ഠ (␪ u Ϫ ␪ d )/(r u Ϫ r d ). This relationship can be combined with (10) so that the downdraft properties in (10) are replaced with the leg mean properties. There were some legs for which (10) do not do a good job. These cases were near the mixed-layer top where the mass-flux approximation did not do a good job representing the flux.  Expressions are needed for ␴ ␪ , ␴ r , and ␳ ␪,r used by the statistical method. Many authors have shown that ␴ ␪ /␪*, and to a lesser extent ␴ r /r* are well-behaved functions of height, when normalized by the appropriate scaling variable (e.g., Willis and Deardorff 1974; Lenschow and Stephens 1980; Khalsa and Greenhut 1987). Berg (2002) introduced a new moisture scaling variable r*,new , which is based on the moisture fluxes at both the top and bottom of the mixed layer: r*,new (z) ϭ  6. Determination of parameters The variables that must be parameterized for each method are listed in Table 2. Those variables for which a physically based parameterization or a similarity parameterization will be developed are identified. The new parameterizations are derived for a subset of BLX96 case-study days chosen at random (Table 3). These ‘‘calibration subset’’ days happened to include legs from each location, although three of the five were the Winfield leg. Only one of the chosen days was clear. This random subset of days was chosen to allow the remaining case-study days to be used for verification of  wЈrЈsfc (1 Ϫ z/z i ) ϩ wЈrЈzi (z/z i ) , w*  (11)  where the subscript sfc indicates a surface value, and z i is the mixed-layer depth. It was shown in section 4, that the correlation of ␪ and r is also a function of height for BLX96. This behavior is consistent with the results reported by both Wyngaard et al (1978) and Wyngaard and LeMone (1980), although this relation may not be universal. These expressions can be used to parameterize the statistical method. Attempts to find an expression relating the correlation between ␪ and r to Bowen ratio (B) were unsuccessful.  TABLE 2. List of methods used to model JFDs, the parameters needed, the basis for the parameterization (physical or similarity theory), and an associated equation number. The symbols have the following definitions: S indicates the slope of an axis indicated by the subscript, ␴ the spread along an axis indicated by the subscript, ␳ a correlation between two variables indicated by the subscripts, a subscript ‘‘sfc’’ indicates surface axis, a subscript ‘‘ez’’ indicates entrainment-zone axis, a subscript B indicates Bowen ratio axis, and a subscript ␪ F indicates solar-forcing axis. Method  Parameters needed  Physically based  Similarity based  Statistical  ␪ , r, ␴␪ , ␴ r , and ␳␪,r  ␴␪ (12), ␴ r (13), and ␳␪,r (14)  B–␪ F  ␪ , r, ␴ B , ␴␪F , S B , and S␪F  S B and S␪F  ␴ B [(15a) or (16a)] and ␴␪F [(15b) or (16b)]  Mixing diagram  ␪ , r, ␴sfc , ␴ez , Ssfc , and Sez  Ssfc [(7) or (10)] and Sez (computed from sounding)  ␴sfc [(17a) or (18a)] and ␴ez [(17b) or (18b)]  822  JOURNAL OF THE ATMOSPHERIC SCIENCES  VOLUME 61  relationship, which is not adjusted to fit BLX96 data, is used for all subsequent calculations of ␴ ␪ /␪*. Sorbjan (1991) also derived an expression for ␴ r /r*. However, as was shown by Berg (2002), the traditional value of r* is not the best choice, so the new expression (11) was derived. Berg (2002) shows that the new scaling does a good job collapsing the BLX96 moisture data to a universal curve (Fig. 8). For the calibration subset of BLX96 data, the least squares best-fit expression for the curve ␴ r /r*,new is  ␴ r /r*,new ϭ 2.8 ϩ 0.25(z/z i ) Ϫ0.81 .  (13)  ␳ ␪,r ϭ 0.39 Ϫ 1.4(z/z i ).  (14)  Expressions are also needed for ␳ ␪,r . The data from BLX96 show a linear variation of ␳ ␪,r with height, so a straight line was fit to the observations (Fig. 1). The best-fit line for the BLX96 data was found to be FIG. 7. Sorbjan’s expression for ␴␪ /␪ (line) and observed values of * of BLX96 case-study days. ␴␪ /␪ (symbols) vs z/zi for the subset *  Sorbjan (1991) suggested an expression for ␴ ␪ /␪* [his Eq. (6)]:  [  ]  ␴␪ (1 Ϫ z/z i ) 4/ 3 (z/z i ) 4/ 3 ϭ c2b ϩ c2t R 4/ 3 , (12) 2/ 3 ␪* (z/z i ) (1 Ϫ z/z i ϩ D) 2/ 3 where R is the Ball ratio (Ball 1960) of entrainmentzone heat flux to surface heat flux, D is the ratio of the entrainment zone depth to z i , and c 2b and c 2t are empirical constants. Sorbjan suggests that c 2b ϭ 2.2 and c 2t ϭ 6. For the calibration subset of BLX96 days, the average value of R was 0.46 and the value of D was 0.4. Figure 7 shows the results from (12) and the BLX96 data using these average values. Sorbjan’s expression does a good job representing BLX96 data, although it may underestimate ␴ ␪ /␪* higher in the mixed layer. This  As pointed out in section 4, and shown in Fig. 3, it appears that this curve is not universal. Wyngaard et al. (1978) reported values of ␳␪,r measured in the maritime boundary layer that were larger than those found during BLX96. Some implications of the variation of ␳␪,r for different parameterizations will be discussed in section 7. b. B–␪F method Expressions are needed for tilt of the axes, and for the spread along the B and ␪ F axes [␴ B (g kg Ϫ1 ) and ␴ ␪ F (K ), respectively] used by the B–␪ F method. Data from the calibration subset of BLX96 flights was used to calculate ␴ B and ␴ ␪ F. Figure 9 shows that, when normalized by ␪* or r*,new , the observed values of ␴ B and ␴ ␪ F each collapse toward a single curve. Using a least squares technique, the best-fit power law was found to be  FIG. 8. Sorbjan’s expression for ␴ r /r (thin line), BLX96 data scaled with standard Deardorff * r scaling, (circles), BLX96 data scaled with (crosses), and BLX96 best fit (13) (thick line) *,new vs z/z i .  1 APRIL 2004  BERG AND STULL  823  FIG. 9. Observed values of (left) ␴␪F/␪ (circles), (right) ␴ B /r ,new (circles), and best-fit curves [(15a) and (b)] * for the BLX96 control * days vs z/z . (lines) i  ␴ B /r*,new ϭ 2.53 ϩ 0.57(z/z i )Ϫ0.65 , and ␴␪ F /␪* ϭ 1.56 ϩ 0.47(z/z i )Ϫ0.56 .  (15a) (15b)  Alternatively, equations for ␴ B and ␴␪ F were derived by SSZ that were based on ␴ ␪ , ␴ r , and ␳ ␪,r , yielding  ΂  ΃  2  ΂S  ΃  2  1 ␴ 2B ϭ (␴ 2␪ Ϫ 2S␪ F ␳␪,r ␴␪ ␴ r ϩ S ␪2 F ␴ 2r ), S B Ϫ S␪ F  and (16a)  ␴ 2␪ F ϭ  B  S␪ F (␴ 2␪ Ϫ 2S B ␳␪,r ␴␪ ␴ r ϩ S 2B ␴ 2r ), (16b) Ϫ S␪ F  where S B is the slope of the B axis and S␪ F is the slope of the ␪ F axis. Equations (16a) and (16b) were used to generate the values of ␴ B and ␴␪ F used in this study, rather than (15a) and (15b), so that results from the B– ␪ F and the statistical method are directly comparable. In other words, differences in the fits are due to differences in theories, rather than differences in the predicted value of the standard deviations. Values of ␴ ␪ , ␴ r and  ␳ ␪,r needed by (16a) and (16b) were calculated from (12), (13), and (14). c. Mixing diagram method The mixing-diagram method requires parameterization of the following: the slope of the surface axis (Ssfc ), the slope of the entrainment-zone axis (Sez ), the standard deviation along the surface axis [␴sfc , (g kg Ϫ1 )], and the standard deviation along the entrainment-zone axis [␴ez , (K)]. Only ␴sfc and ␴ez need to be parameterized using similarity methods because Ssfc and Sez are either known from theory or model output. Using data from the BLX96 calibration subset days, least-square best-fit power-law functions can be fit to the curves of ␴ez /␪* and ␴sfc /r*,new :  ␴ez /␪* ϭ 2.0 ϩ 8.2 ϫ 10Ϫ4 (z/z i )Ϫ1.8 and ␴sfc /r*,new ϭ 2.3 ϩ 1.1 ϫ 10Ϫ2 (z/z i )Ϫ1.6 .  (17a) (17b)  Figure 10 shows (17a) and (17b) plotted as a function  FIG. 10. Observed (a) ␴ez /␪ (17a) and (b) ␴sfc /r ,new (17b) on the BLX96 case study days vs z/z i . * *  824  JOURNAL OF THE ATMOSPHERIC SCIENCES  TABLE 4. Verification subset of BLX96 case-study days used to perform the statistical fits, and the associated location and cloudcover amount. Date 16 23 31 2 13  Jul Jul Jul Aug Aug  Location  Cu coverage  Meeker Winfield Winfield Meeker Lamont  Scattered Cu Scattered Cu Scattered Cu None None  of z/z i . The agreement seems to be good for the range of z/z i observed during BLX96. However, it is not clear that ␪* and r*,new are the relevant scales. While the value of ␴ez does have units of temperature, it represents a relationship that is dependant on both ␪ and r. Likewise, ␴sfc is also a function of both ␪ and r. Efforts were made to define two new scaling variables using the projected variables c and m. First, the observed heat and moisture fluxes were converted to fluxes of c and m using SSZ’s (14a) and (14b). The fluxes of c and m, along with w* were used to define c* and m* using definitions analogous to those used by standard Deardorff scaling. Unfortunately, this scaling did not improve the results. As an alternative to (17a) and (17b), ␴sfc and ␴ez can be defined in terms of ␴ ␪ , ␴ r , and ␳ ␪,r using definitions derived by SSZ as  ␴ 2ez ϭ  ΂  ΃  2  Ssfc 2 (␴ 2␪ Ϫ 2Ssfc ␳␪,r ␴␪ ␴ r ϩ S sfc ␴ 2r ), Sez Ϫ Ssfc  and (18a)  ␴  2 sfc  ΂  ΃  2  1 2 ϭ (␴ 2␪ Ϫ 2Sez ␳␪,r ␴␪ ␴ r ϩ S ez ␴ 2r ). (18b) Sez Ϫ Ssfc  Thus, the value of ␴sfc , ␴ez , and ␳ ␪,r can be found from a combination of (18) with (12), (13), and (14) to account for the influences of both ␪ and r in each standard deviation. As described in the previous section, (18a) and (18b) give results that are more consistent with the methods used to fit the statistical method. 7. Verification against independent data A ‘‘verification subset’’ was created from all the remaining BLX96 cases that were not previously chosen at random to help with calibration. The verification subset is listed in Table 4. JPDFs were created for the verification subset of BLX96 cases using theories presented in section 5 and the equations presented in section 6. Namely, values of ␴ ␪ , ␴ r , and ␳ ␪,r needed by the statistical method were determined from (12), (13), and (14), respectively. SSZ’s theories were used to predict S B and S ␪ F for the B–␪ F method. Equations (16a) and (16b) were used to predict ␴ B and ␴␪ F, respectively. For the mixing diagram method, the value of Ssfc was predicted from (10), while ␴sfc and ␴ez were predicted from Eqs. (18a) and (18b). Relevant soundings were used to  VOLUME 61  calculate Sez by fitting a line to the entrainment zone observations, and assuming Sez varied linearly in time between the soundings. The similarity relationships derived in section 6 [Eqs. (15a), (15b), (17a), and (17b)] could have been used to find ␴ B , ␴ ␪ F, ␴ez , and ␴sfc based on observations from BLX96. However, using (16a), (16b), (18a), and (18b), which are formulated in terms of ␴ ␪ , ␴ r , and ␳ ␪ ,r , to predict the values of ␴ B , ␴ ␪ F, ␴ez , and ␴sfc has an advantage: results from each method are more consistent with each other. In other words, differences in the fits are due to differences in the theories, rather than differences in the predicted value of the standard deviations. Figure 11 shows sample fits for a single leg flown on 27 July 1996. This leg was flown at a height of 0.62z i . In this case, the tilt of the JFD is negative. The B–␪ F method is unable to capture this behavior, while both the statistical and mixing diagram methods can. For this case the difference between the statistical method and the mixing diagram method is very small. The peak frequency of all methods is less than the observed JFD, although in this example the differences are small. The observed peak frequency was 0.083, while the predicted frequencies for the statistical, B–␪ F , and mixing diagram method were 0.079, 0.060, and 0.079, respectively. Similar results were found for all BLX96 JPDFs; namely, each method underestimated the peak frequency. JPDFs were computed for all the BLX96 verification subset cases, and the error of each method was computed using the method of SSZ: err ϭ  ͸ | fˆ  obs, j  Ϫ fˆparam, j |,  (19)  where ˆfobsj is the normalized frequency observed for the jth bin, and ˆfparam.j is the normalized frequency from the jth bin of the parameterization. The results from each method were complied and are shown in Fig. 12. It appears that both the statistical method and the mixing diagram method work better than the B–␪ F method. There is only one case where the error from the B–␪ F method is close to the one-to-one line. There is little difference in the error predicted by the statistical method and the mixing diagram method; most of the points are very close to the 1:1 line. As pointed out by SSZ, an error of 1.0 shows no skill. Both the statistical and mixing diagram method are better than no fit for most of the flight legs; only one flight leg had an error in the statistical or mixing diagram method that was greater than 1.0. Many of the B–␪ F JPDFs had errors larger than 1.0, Berg (2002) attributed these errors to shortcomings associated with the prediction of the slopes S B and S ␪ F. The good agreement between the statistical method and the mixing diagram method is encouraging. For BLX96 case-study days the statistical method should do a good job because each of the parameters—␴ ␪ , ␴ r , and ␳ ␪,r—are well behaved. However, the mixing diagram method has several advantages over the statistical method. First, the mixing diagram method is more robust  1 APRIL 2004  BERG AND STULL  FIG. 11. Observed (broken lines) JFDs and parameterized (solid lines) JPDFs using (a) statistical, (b) B–␪ F , and (c) mixing diagram methods for 27 Jul 1996, where z/z i ϭ 0.62 at 1211 LST.  FIG. 12. Errors associated with the parameterized statistical method JPDF vs error associated with either the mixing diagram method (circles) or B–␪ F (crosses) JPDFs. Points to the right of the one-to-one line have a smaller error than the statistical method.  825  826  JOURNAL OF THE ATMOSPHERIC SCIENCES  FIG. 13. Mixing diagram predicted ␳ ␪,r vs observed ␳ ␪,r for BLX96 flight legs. The solid line is the 1:1 line.  than the statistical method because the tilt of the JFD is a function of the slopes of the surface and entrainment-zone axes. Second, only ␴ez and ␴sfc must be parameterized in the mixing diagram method, compared to ␴ ␪ , ␴ r , and ␳ ␪,r that must be parameterized in the statistical method. But, if (18a) and (18b) are used to estimate ␴ez and ␴sfc , then a third variable (␳ ␪,r ) is needed and both the statistical method and the mixing diagram method depend on ␴ ␪ , ␴ r , and ␳ ␪,r . At first glance, inclusion of ␳ ␪,r would seem to reduce the advantage of the mixing diagram method over the statistical method. Closer examination of the terms in the equations for 2 ␴ ez2 [(18a)] and ␴ sfc [(18b)] indicate that, for BLX96, the term that includes ␳ ␪,r is nearly an order of magnitude smaller than the other terms, so the dependence of ␴ ez2 2 and ␴ sfc on ␳ ␪,r is weak and the mixing diagram method is not sensitive to the observed value of ␳ ␪,r . 2 Given the weak dependence of ␴ ez2 and ␴ sfc on the observed value of ␳ ␪,r , the JPDFs found using the mixing diagram method can be used to predict ␳ ␪,r . There is good agreement between ␳ ␪,r predicted from the mixing diagram method JPDF and ␳ ␪,r observed during BLX96 (Fig. 13). This results shows that the axes selected by the mixing diagram method account for a significant amount of the observed correlation between ␪ and r. For cases where ␳ ␪,r calculated from the BLX96 observations and ␳ ␪,r predicted from the mixing diagram JPDFs are both negative, ␳ ␪,r predicted from the mixing diagram JPDFs is slightly smaller than the observed. It is not clear why this is the case, but a number of these observations occur for legs near the top of the boundary layer. This could be related to errors in the entrainmentzone slope, which are interpolated to the time of the relevant flight leg. The skill of the mixing diagram predictions of ␳ ␪,r is a significant advantage because ␳ ␪,r may not be universal. 8. Conclusions A new field experiment was conducted in which JFDs of ␪ versus r were measured at heights ranging from  VOLUME 61  0.02 to 0.8z i in the daytime convective boundary layer. Results from BLX96 were consistent with findings of SSZ. It was found that JFDs near the surface had a positive tilt, in other words, warm parcels also tended to be relatively moist, while those JFDs measured higher in the mixed layer had a negative tilt because warm (entrained) parcels tended to be relatively dry. Each distribution was found to be skewed, and more peaked than a Gaussian distribution. A new parameterization to describe the tilt of the JFDs was introduced. This new parameterization describes the mixed-layer JFDs as a mixture of air originating from near the surface, from the entrainment zone and from the environment. The relationship of these three source regions was used to tilt the JFDs. Similarity relationships were derived to describe the spread of the JFD along each axis. Using a calibration subset of observations from BLX96 the derived parameterizations were compared to a verification subset, also from BLX96. Fortunately, BLX96 flights were flown at a three different locations, so our findings are more likely applicable to other relatively flat geographic locations. This new parameterization was compared to two classic parameterizations. The new parameterization was more successful than the method of SSZ, and about as successful as that of WM92. The primary benefit of the new parameterization are the physical relationships used to predict the tilt of the JFD and skill predicting ␳ ␪,r based only on the slopes of the two new axes. The primary disadvantage is the requirement of the two additional variables: the slope of the surface and entrainment zone axes. Acknowledgments. This work was completed as part of the lead authors doctoral research at the University of British Columbia (UBC), where funding was provided by a UBC Graduate Fellowship and the Canadian Climate Research Network through grants from the Natural Science and Engineering Research Council (NSERC), Meteorological Service of Canada (MSC), and the Canadian Foundation for Climate and Atmospheric Science. Additional support was provided by the Geophysical Disaster Computational Fluid Dynamics Centre, other grants from NSERC and Environment Canada, and a Laboratory Directed Research and Development Program of the U.S. Department of Energy under Contract DE-AC06-76RLO at the Pacific Northwest National Laboratory, which is operated for the U.S. Department of Energy by Battelle Memorial Institute. The BLX96 field program was funded by the U.S. National Science Foundation (NSF) under Grant ATM9411467. E. Santoso and J. Hacker are thanked for their work as airborne scientists during BLX96. The staff and flight crew of the University of Wyoming King Air aircraft (sponsored by NSF) helped to make BLX96 a success. The comments of three anonymous reviewers greatly improved this manuscript.  1 APRIL 2004  827  BERG AND STULL  APPENDIX  REFERENCES  Slopes Used in the B–␪ F Method  Ball, F. K., 1960: Control of inversion height by surface heating. Quart. J. Roy. Meteor. Soc., 86, 483–494. Berg, L. K., 2002: A simple parameterization coupling the convective daytime boundary layer and fair-weather cumuli. Ph.D. dissertation. Atmospheric Science Programme, Department of Earth Sciences, The University of British Columbia, 166 pp. [Available from The University of British Columbia, 6339 Stores Rd., Vancouver, BC, V6T 1Z4, Canada.] ——, R. B. Stull, E. Santoso, and J. P. Hacker, 1997: Boundary layer experiment–1996, airborne scientist flight log. Bound. Layer Tech. Rep. 97-1, 116 pp. [Available from Roland Stull, Atmospheric Science Programme, Department of Earth and Ocean Sciences, University of British Columbia, 6339 Stores Rd., Vancouver, BC V6T-1Z4, Canada.] Betts, A. K., 1975: Parametric interpretation of trade-wind cumulus budget studies. J. Atmos. Sci., 32, 1934–1945. ——, 1976: Modeling subcloud layer structure and interaction with a shallow cumulus layer. J. Atmos. Sci., 33, 2363–2382. ——, 1982: Saturation point analysis of moist convective overturning. J. Atmos. Sci., 39, 1484–1505. ——, 1984: Boundary layer thermodynamics of a high plains severe storm. Mon. Wea. Rev., 112, 2199–2211. ——, 1985: Mixing line analysis of clouds and cloudy boundary layers. J. Atmos. Sci., 42, 2751–2763. ——, 1992: FIFE atmospheric boundary layer budget methods. J. Geophys. Res., 97, 18 523–18 531. ——, and B. A. Albrecht, 1987: Conserved variable analysis of the convective boundary layer thermodynamic structure over tropical oceans. J. Atmos. Sci., 44, 83–99. Coulman, C. E., 1978: Boundary layer evolution and nocturnal inversion dispersal—Part II. Bound.-Layer Meteor., 14, 493–513. Crum, T. D., and R. B. Stull, 1987: Field measurements of the amount of surface layer air versus height in the entrainment zone. J. Atmos. Sci., 44, 2743–2753. Cuijpers, J. W. M., and A. A. M. Holtslag, 1998: Impact of skewness and nonlocal effects on scalar buoyancy fluxes in convective boundary layers. J. Atmos. Sci., 55, 151–162. Deardorff, J. W., and G. E. Willis, 1985: Further results from a laboratory model of the convective planetary boundary layer. Bound.-Layer Meteor., 32, 205–236. Greenhut, G. K., and S. J. S. Khalsa, 1982: Updraft and downdraft  Schrieber et al. (1996) derived equations for ␪Ј and rЈ, the perturbation values from the means, as [their (4) and (5)]  ␪Ј ϭ  B␪ F , 1ϩB  rЈ ϭ  ␥␪ F , 1ϩB  and  (A1) (A2)  where ␥ ϭ C p /L ␷ is the psychrometric constant and ␪ F is the solar forcing (defined in section 5). They modified these equations to give slopes in terms of ␪ ␷ and zLCL . However, (A1) and (A2) can be combined to give the slopes in terms of ␪ and r. Like the work of Schrieber et al. (1996), lines of constant B can be found by dividing (A2) by (A1): rЈ/␪Ј ϭ ␥ /B.  (A3)  Lines of constant solar forcing can also be found from (A1) and (A2). First, (A2) is solved in terms of B. This value is substituted into (A1) leaving rЈ ϭ Ϫ␥␪Ј ϩ ␥␪ F ,  (A4)  which is the equation for the ␪ F lines. In ␪ versus r space all lines of constant solar forcing have the same slope Ϫ␥ . This result is different than lines of constant solar forcing in ␪ ␷ versus zLCL space, which are a function of the reference state. Figure A1 shows an example of (A3) and (A4) plotted in ␪ versus r space. Figure A1 appears to be a mirror image of Schrieber et al.’s Fig. 2 because large values of zLCL correspond to small values of r.  FIG. A1. Axes of B (dashed lines) and ␪ F (solid lines) for a range of values for a reference state defined as ␪ ϭ 300 K and r ϭ 11 g kg Ϫ1 , after Schrieber et al. (1996).  828  JOURNAL OF THE ATMOSPHERIC SCIENCES  events in the atmospheric boundary layer over the equatorial Pacific Ocean. J. Atmos. Sci., 39, 1803–1818. ——, and ——, 1987: Convective elements in the marine atmospheric boundary layer. Part I: Conditional sampling statistics. J. Climate Appl. Meteor., 26, 813–822. Grossman, R. L., 1984: Bivariate conditional sampling of moisture flux over a tropical ocean. J. Atmos. Sci., 41, 3238–3252. Hanson, H. P., 1984: On mixed layer modeling of the stratocumulustopped marine boundary layer. J. Atmos. Sci., 41, 1226–1234. Holland, J. Z., 1973: A statistical method for analyzing wave shapes and phase relationships of fluctuationg geophysical variables. J. Phys. Oceanogr., 3, 139–155. Khalsa, S. J. S., and G. K. Greenhut, 1985: Conditional sampling of updrafts and downdrafts in the marine atmospheric boundary layer. J. Atmos. Sci., 42, 2550–2562. Larson, V. E., R. Wood, P. R. Field, J.-C. Golaz, T. H. Vonder Haar, and W. R. Cotton, 2001: Small-scale and mesoscale variability of scalars in cloudy boundary layers: One-dimensional probability density functions. J. Atmos. Sci., 58, 1978–1994. Lenschow, D. H., 1970: Airplane measurements of the planetary boundary layer structure. J. Appl. Meteor., 9, 874–884. ——, and P. L. Stephens, 1980: The role of thermals in the convective boundary layer. Bound.-Layer Meteor., 19, 509–532. Lewellen, W. S., and S. Yoh, 1993: Binormal model ensemble partial cloudiness. J. Atmos. Sci., 50, 1228–1237. Mahrt, L., and J. Paumier, 1982: Cloud-top entrainment instability observed in AMTEX. J. Atmos. Sci., 39, 622–634. ——, and ——, 1984: Heat transport in the atmospheric boundary layer. J. Atmos. Sci., 41, 3061–3075. Manton, J. R., 1975: Penetrative convection due to a field of thermals. J. Atmos. Sci., 32, 2272–2277. Nicholls, S., and M. A. LeMone, 1980: The fair weather boundary layer in GATE: The relationship of subcloud fluxes and structure to the distribution and enhancement of cumulus clouds. J. Atmos. Sci., 37, 2051–2067. Paluch, I. R., 1979: Entrainment mechanism in Colorado cumuli. J. Atmos. Sci., 36, 2467–2478. Press, W. H., S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, 1992: Numerical Recipes in C: The Art of Scientific Computing. Cambridge University Press, 994 pp. Randall, D. A., Q. Shao, and C.-H. Moeng, 1992: A second-order bulk boundary-layer model. J. Atmos. Sci., 49, 1903–1923. Santoso, E., and R. B. Stull, 1998: Wind and temperature profiles in the radix layer, the bottom fifth of the convective boundary layer. J. Appl. Meteor., 37, 545–558. ——, and ——, 2001: Similarity equations for wind and temperature  VOLUME 61  profiles in the radix layer, at the bottom of the convective boundary layer. J. Atmos. Sci., 58, 1446–1464. Schrieber, K., R. Stull, and Q. Zhang, 1996: Distributions of surfacelayer buoyancy versus lifting condensation level over a heterogeneous land surface. J. Atmos. Sci., 53, 1086–1107. Schumann, U., and C.-H. Moeng, 1991: Plume fluxes in the clear and cloudy convective boundary layers. J. Atmos. Sci., 48, 1746– 1757. Siebesma, A. P., and J. W. M. Cuijpers, 1995: Evaluation of parametric assumptions for shallow cumulus convection. J. Atmos. Sci., 52, 650–666. Sorbjan, Z., 1991: Evaluation of local similarity functions in the convective boundary layer. J. Appl. Meteor., 30, 1565–1583. Stokes, G. M., and S. E. Swartz, 1994: The atmospheric radiation measurement (ARM) program: Programmatic background and design of the cloud and radiation testbed. Bull. Amer. Meteor. Soc., 75, 1201–1221. Stull, R. B., 1994: A convective transport theory for surface fluxes. J. Atmos. Sci., 51, 3–22. ——, E. Santoso, L. Berg, and J. Hacker, 1997: Boundary Layer Experiment 1996 (BLX96). Bull. Amer. Meteor. Soc., 78, 1149– 1158. Wang, S., and B. Stevens, 2000: Top-hat representation of turbulence statistics in cloud-topped boundary layers: A large eddy simulation study. J. Atmos. Sci., 57, 423–441. Williams, A. G., and J. M. Hacker, 1992: The composite shape and structure of coherent eddies in the convective mixed layer. Bound.-Layer Meteor., 61, 213–245. Willis, G. E., and J. W. Deardorff, 1974: A laboratory model of the unstable planetary boundary layer. J. Atmos. Sci., 31, 1297– 1307. Wyngaard, J. C., and M. A. LeMone, 1980: Behavior of the refractive index structure parameter in the entraining convective boundary layer. J. Atmos. Sci., 37, 1573–1585. ——, and C.-H. Moeng, 1992: Parameterizing turbulent diffusion through the joint probability density. Bound.-Layer Meteor., 60, 1–13. ——, W. T. Pennell, D. H. Lenschow, and M. A. LeMone, 1978: The temperature–humidity covariance budget in the convective boundary layer. J. Atmos. Sci., 35, 47–58. Xu, K.-M., and D. A. Randall, 1996: Evaluation of statistically based cloudiness parameterizations used in climate models. J. Atmos. Sci., 53, 3103–3119. Young, G. S., 1988: Turbulence structure of the convective boundary layer. Part II: Phoenix 78 aircraft observations of thermals and their environment. J. Atmos. Sci., 45, 727–735.  


Citation Scheme:


Citations by CSL (citeproc-js)

Usage Statistics



Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items