1 JULY 1999 GREISCHAR AND STULL 2201 Convective Transport Theory for Surface Fluxes Tested over the Western Pacific Warm Pool LAWRENCE GREISCHAR Department of Atmospheric and Oceanic Sciences, University of Wisconsin—Madison, Madison, Wisconsin ROLAND STULL Atmospheric Science Programme, Department of Geography, University of British Columbia, Vancouver, British Columbia, Canada (Manuscript received 21 July 1997, in final form 9 July 1998) ABSTRACT Turbulent flux measurements from five flights of the National Center for Atmospheric Research Electra aircraft during the Tropical Oceans and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE) are used to test convective transport theory (CTT) for a marine boundary layer. Flights during light to moderate winds and under the clearest sky conditions available were chosen. Fluxes of heat, moisture, and momentum were observed by the eddy-correlation method. Mean kinematic values for the observed sensible and latent heat fluxes and momentum flux were 0.0061 K m sϪ1 , 0.0313 g kgϪ1 m sϪ1, and 0.0195 m 2 sϪ2 , respectively. For the range of mixed-layer wind speeds (0.8–8.4 m s Ϫ1 ) studied here, the version of CTT that includes the mixed effects of buoyant and shear-driven transport give a better fit to the observations than either the COARE bulk algorithm or the pure free-convection version of CTT. This is to be expected because both of those latter parameterizations were designed for light winds (Ͻ5 m sϪ1 approximately). The CTT empirical coefficients listed in Table 3 exhibited slight sensitivity to the COARE light flux conditions, compared to their previous estimates during larger fluxes over land. For example, COARE heat fluxes were roughly 10 times smaller than previous land-based flux measurements used to calculate CTT coefficients, but the corresponding empirical mixed-layer transport coefficients were only 3% smaller. COARE momentum fluxes were also roughly 10 times smaller, but the CTT coefficients were about four times smaller. The greater variation in momentum coefficient may be due, in part, to insufficient flight-leg length used to compute momentum fluxes, to uncertainties in the effects of the ocean surface current and waves, or perhaps to roughness differences. 1. Introduction The western Pacific warm pool was identified during the Tropical Oceans and Global Atmosphere (TOGA) program as an important contributor to global climate variability. A Coupled Ocean–Atmosphere Response Experiment (COARE) was therefore designed, which included measurement of interfacial fluxes as one of the key elements. Webster and Lukas (1992) emphasized the importance of flux parameterizations especially for low-wind regimes. The convective transport theory (CTT) proposed by Stull (1994) is just such a parameterization. This theory parameterizes surface fluxes of heat, moisture, and momentum as the product of an empirical transport coef- Corresponding author address: Prof. Roland Stull, Atmospheric Science Programme, Department of Geography, 1984 West Mall, University of British Columbia, Vancouver, BC V6T 1Z2, Canada. E-mail: rstull@geog.ubc.ca ᭧ 1999 American Meteorological Society ficient, a ‘‘buoyancy velocity scale,’’ and the difference between mean surface and interior of mixed-layer values of the respective property. The COARE observations made within the intensive flux array (IFA) in conjunction with extensive surface measurements by ships and buoys contain the information necessary to test CTT under marine conditions. Thirty-two flights of the National Center for Atmospheric Research (NCAR) Electra aircraft were made during the intensive operations period (IOP) of November 1992 and February 1993, many of which had lowlevel boundary layer legs in the area of concurrent surface measurements made by ships and at the improved meteorological instrumentation (IMET) buoy at 1Њ45ЈS, 156ЊE. The published values of the empirical transport coefficients for CTT were originally estimated using airborne flux measurements from the Boundary Layer Field Experiment 1983 (BLX83) in Oklahoma (roughness length: 0.05 m; latitude: 35ЊN; vegetation: mixed 2202 JOURNAL OF THE ATMOSPHERIC SCIENCES pasture and crops; season: boreal spring). Data from another land-based experiment (Koorin) conducted under somewhat different conditions in Australia (roughness length: 0.4 m; latitude: 16ЊS; vegetation: uniform sparse trees; season: austral winter) were used to confirm the theory and the values of the coefficients. The objective of this paper is to test this theory and its limits under the marine conditions of the western Pacific warm pool (roughness length: 3 ϫ 10Ϫ5 m; latitude: 2ЊS; marine conditions; season: austral summer) using the Electra aircraft data obtained during the COARE IOP. An overview of CTT is given in section 2. The Electra aircraft data and ancillary observations are described in section 3. Section 4 describes the data processing steps, and CTT is tested against the Electra data in section 5. Concluding remarks are made in section 6. 2. Model description During conditions of free convection and calm or light winds, turbulence is generated by buoyant thermals, rather than mechanically by surface-layer wind shear (Stull 1997). Convective thermals efficiently transport heat, moisture, and momentum between the surface and the mid–mixed layer, with a rate of vertical transport proportional to the Deardorff velocity w* ϭ [(g/T )z iwЈЈ s ]1/3 , where g is gravitational acceleration, T is the mean absolute virtual temperature, z i is the mixed-layer (ML) depth, and wЈЈ s is the surface kinematic vertical flux of virtual potential temperature, which is proportional to a buoyancy flux. While traditional drag laws fail in the limit of zero mean wind, surface fluxes during free convection nonetheless can be described by CTT (Stull 1994) over a range of calm to light winds: u*2 ϭ C*D w* MML (1) wЈ Јs ϭ C*H w*⌬ (2) wЈrЈs ϭ C*E w*⌬r, (3) where u*2 , wЈЈs , and wЈrЈs , are the surface kinematic fluxes of momentum, heat, and water vapor, respectively, and u* is also known as the friction velocity. The last factor in each equation above is a difference between the surface skin condition and the mid-mixedlayer value: ⌬ ϭ skin Ϫ ML , ⌬r ϭ rskin Ϫ rML , and ⌬M ϭ MML Ϫ Mskin ϭ MML , where Mskin ϭ 0 by definition. By using skin and mid-mixed-layer values in the above formulations, it was proposed by Stull (1994) that the CTT parameterization should be independent of surface aerodynamic roughness length z o . This assumption will be tested here. The empirical mixed-layer transport coefficients for momentum, moisture, and heat were found (Stull 1994) to be C* D ϭ 0.023 Ϯ 0.007 and C* E ഠ C* H ϭ 0.0063 Ϯ 0.0016, respectively, based on measurements over VOLUME 56 land (United States and Australia), where the error tolerance is given as a standard deviation. Stull had insufficient data to calculate C* E but hypothesized that it might be equal to C* H . We will test that hypothesis here. Kustas et al. (1996) has demonstrated that surfaces with complex geometry, such as forests, require corrections to the above coefficients to compensate for shading related to solar illumination versus viewing angles. Equations (1)–(3) are in implicit form, because the fluxes on the left side of Eqs. (2) and (3) depend on w*, which itself is a function of heat and moisture fluxes via wЈЈ s ഠ wЈЈs (1 ϩ 0.61r) ϩ 0.61wЈrЈs . However, (1)–(3) can be rearranged into explicit form (Stull 1994): u*2 ϭ b D w B MML (4) wЈ Јs ϭ b H w B⌬ (5) wЈrЈs ϭ b E w B⌬r, (6) where convective transport coefficients (b D , b H , b E ) are related to the mixed-layer coefficients by b D ϭ 3/2 1/2 C* DC 1/2 *H, b H ϭ C *H, and b E ϭ C* EC *H. The buoyancy velocity scale is wB ϵ g z ⌬ T i B 1/2 . (7) The buoyancy temperature difference is ⌬ B ϵ ⌬, for  ϵ wЈЈ s /wЈЈs , which can be approximated in explicit form as ⌬ B ഠ ⌬(1 ϩ 0.61r s ) ϩ 0.61 s⌬r. (8) For convection over land (Oklahoma and Australia) it was shown (Stull 1994) both theoretically and empirically that ⌬ B ഠ ⌬ , (9) where ⌬ equals the difference of virtual potential temperatures between the surface skin and middle mixed layer. Equation (9) will be tested here for a maritime tropical atmosphere. This buoyancy velocity can also be used to parameterize the surface flux of virtual potential temperature: wЈЈ s ϭ b H w B⌬ B . (10) A mixed-layer Richardson number was defined as R* ϵ [w B /MML ] 2 , (11) where Stull (1994) suggested that even in noncalm conditions, turbulence is in quasi–free convection whenever R* Ͼ 3. This free convection condition will also be examined here for tropical maritime boundary layers. Another prediction of CTT is that a diagnostic expression for mixed-layer depth is possible: zi ϭ 1 wЈ Јs . (g/T )⌬ B b H⌬ This, too, will be tested. 2 (12) 1 JULY 1999 GREISCHAR AND STULL 2203 TABLE 1. Boundary layer and surface-flux data from five flights of the NCAR Electra aircraft during the TOGA COARE IOP. Variables include start time (Local ϭ UTC ϩ 11 h) of beginning of the approximately 15-min low-level horizontal flight leg; flight ID, flight identification of Electra aircraft; leg ID, identification for each near-surface horizontal flight leg; z leg, mean altitude of the flight leg above sea level; zi, depth of mixed layer; MML , ML , and rML , wind speed, mean potential temperature (relative to 100 kPa), and mixing ratio in the interior of the mixed layer; Psurf , Tskin , and skin , surface pressure, temperature, and potential temperature (relative to 100 kPa) of the surface radiometric skin temperature; rsat (Tskin), saturated mixing ratio at the surface skin temperature; ⌬ ϭ skin Ϫ ML ; ⌬r ϭ rsat (Tskin ) Ϫ rML ; Here, ͗wЈЈs ͘, ͗wЈЈvs͘, ͗wЈrЈs ͘, and u* 2 are the leg-averaged surface kinematic vertical fluxes of heat, virtual temperature, moisture, and momentum, respectively. Means and standard deviations of all numbers in each column are shown at the bottom. Symbol t represents . Finally, for those conditions where mechanical generation of turbulence by wind shear is great enough to be comparable to convective turbulence, a ‘‘mixed’’ mechanical-buoyant version of CTT was also proposed: u*2 ϭ (C DML MML ϩ b DML w B )MML (13) wЈ Јs ϭ (C HML MML ϩ b HML w B )⌬ (14) wЈrЈs ϭ (C BML MML ϩ b BML w B )⌬r, (15) where Stull (1994) empirically found the drag coefficient values for momentum and heat to be C DML ϭ 0.0035 and C HML ϭ 0.001, respectively, with the convective transport coefficients for this mixed-convection case of b DML ϭ 0.0007 and b HML ϭ 0.000 25, respectively. Based on (7), (11), and (14), the buoyancy velocity and Deardorff velocity scales are related by 1/3 w* ϭ (b HML ϩ C HML RϪ1/2 * ) wB , (16) which will be tested here. 3. Observations Thirty low-level flight segments (Table 1) from five flights of the Electra aircraft flown during the COARE IOP, between November 1992 and February 1993, were used in this study to estimate surface fluxes of heat, 2204 JOURNAL OF THE ATMOSPHERIC SCIENCES moisture, and momentum. These data also provided estimates of the surface and middle boundary layer parameters needed to test the theory. The flight segments were chosen to avoid precipitating conditions, gust-front boundaries, and to be as cloud free as possible. An attempt was made to use flight data over a variety of wind speeds in order to test the CTT over a range from near-calm free-convective conditions to conditions in which mechanically generated turbulence transport associated with mean wind shear contributes significantly to the still-dominant buoyant transport. Surface data from the IMET buoy (Release 1.0b of WHOI buoy data from Dr. Robert Weller, Woods Hole Oceanographic Institute) were useful for choosing appropriate flight days as were satellite and radar images and surface and upperair analyses (TCIPO 1993; Asencio et al. 1993; Bond and Alexander 1994). A corrected version of the Electra’s radiometric surface temperatures was obtained (Y. Serra 1997, personal communication). COARE shipbased current data obtained from the archive at the University of Hawaii were used to estimate the effects of surface currents on the model. Flight conditions varied from clear with light winds to partly cloudy with winds up to 8 m sϪ1 . Flights RF04, RF05, RF06, and RF25 were class 0 (described in the TOGA COARE Operations Plan as having clouds and very small, short-lived showers with areas less than 10 km 2 covering not more than 5% of the flight area; TCIPO 1992) and flight RF16 was class 1 (indication of convective cores less than 25 km 2 and lasting less than 30 min and no mesoscale organization greater than 50 km or long-lived stratiform rain areas greater than 2500 km 2 within the target domain). Rainfall data (15 min averages) from the R/V Franklin (Release 1.0 from Bradley and Coppin, CSIRO Centre for Environmental Mechanics, Canberra, Australia), which was in the vicinity of flights RF04, RF05, RF06, and RF25, indicated no precipitation at its location during these flights. Surface radar data were available for flights RF04 and RF16, which roughly delineated cloudy areas. Sensors on the Electra detected some liquid water content only during short sections of leg 2 of flight RF06. These sections were excluded from the analysis. Fast-rate measurements of the global positioning system (GPS)–corrected wind components, potential temperature, and mixing ratio were used to compute the eddy-correlation surface fluxes of heat, moisture, and momentum. These variables are described in detail in Miller and Friesen (1989). Slow-rate measurements of aircraft position and time, surface skin temperature, and an estimate of surface pressure were used to map lowlevel flight legs and determine the pressure and temperature conditions at the surface. Estimates of mixedlayer properties were made where the aircraft descended to or ascended from low-altitude flight legs, resulting in several vertical profiles through the mixed layer during each day’s flight pattern. Some flight patterns also had legs at higher levels within the mixed layer. VOLUME 56 FIG. 1. Low-level flight paths show locations of flight legs and vertical profiles (slant ascent and descent flights, shown by letter V with subscript indicating number of profiles at each approximate location) used in this study. The position of the IMET buoy and the approximate locations of the ship R/V Franklin (ship symbol) during the times of flights RF04, RF05, RF06, and RF25 are shown. Approximate wind direction during each flight is indicated by an arrow near the flight label. The boundary of the IFA is indicated by a dashed line. 4. Data processing Processing of the flight data began by producing twoand three-dimensional plots of the aircraft’s flight path in order to identify suitable low-level flight legs. Figure 1 shows the locations of the low-level legs studied from the five separate flights. The IMET buoy is shown, as are the approximate locations of the R/V Franklin during flights RF04, RF05, RF06, and RF25. An elevation versus time plot of the on-station part of a flight was also used to define the start and finish times of each flight leg, and of the slant ascent or descent sounding that often occurred between the low-level flight legs. Fast-rate values for the three wind components, mixing ratio, and potential temperature were extracted for each leg and plotted (as in Fig. 2). The following restrictions (Sun et al. 1996) were placed on the data: 1) differences of flight levels between adjacent sampling points were required to be less than 5 m; and 2) roll and pitch angles were restricted to Ϫ4Њ to ϩ4Њ and 0Њ to 3Њ, respectively. Also as in Sun et al. (1996), points in the fast-rate data that differed by more than four standard deviations from the mean along 1 km of flight were replaced by using a linear interpolation between good data values. Fluxes calculated using this ‘‘despiked’’ data were mostly changed by less than 1% with a maximum difference of 5%. The fast rate data are sampled at 20 Hz, which given an aircraft speed of approximately 100 m sϪ1 results in a short wavelength sampling limit of approximately 10 m. In order to remove mesoscale effects of diameters larger than 2500 m and corresponding to wavelengths 1 JULY 1999 GREISCHAR AND STULL FIG. 2. Times series plots of (a) potential temperature, (b) mixing ratio, (c) wind speed, and (d) wind direction for leg 1 of flight RF16. Time (in s) from 0000 UTC on 16 December 1992. Every fifth point of the 20-Hz data is plotted. greater than 5 km, a 50-s moving linear average and trend was subtracted from the fast-rate data in a manner similar to Otles and Young (1996). Eddy fluxes of heat, moisture, and momentum were then computed from 18 km or longer sections of these filtered flight legs using a program that first removed any remaining trend. Thicknesses of the mixed layer were estimated from vertical profiles of virtual potential temperature, mixing ratio, and winds measured during aircraft slant ascent and descent soundings at the ends of the low-level flight legs. One-second averages of the fast data were plotted in these vertical profiles. Radar altitudes were used up to their limit of about 750 m, with pressure altitudes approximately adjusted to match the radar data, used above this level. Only one mixed-layer thickness estimate was greater than 750 m for the cases studied here. Generally, the vertical profile of virtual potential temperature gave the best indication of mixed-layer thickness, with a nearly constant value of up to the mixedlayer depth z i . The locations of these slant soundings are shown in Fig. 1. Surface values of temperature, pressure, and moisture were estimated as follows. A corrected version of 2205 the Electra radiometric surface temperatures (Y. Serra 1997, personal communication) includes corrections for emissivity, reflectivity, transient trends, and offset were averaged to obtain an initial estimate of Tskin . An approximate correction for absorption due to moisture content between the aircraft and surface was estimated from the rate of change of the radiometric surface temperature with height during slant ascent and descent soundings (Fig. 1). This rate appeared to be fairly linear over the lowest several hundred meters of elevation and varied from ϩ0.06 K to ϩ0.14 K per 100-m decrease in elevation. A mean value of ϩ0.1 K per 100 m was used, resulting in corrections for this effect of less than 0.1 K for the 32–66-m flight elevations. A rough comparison of these corrected sea surface temperatures with concurrent ship values from the R/V Franklin showed occasional agreement but with a tendency for the ship observations to be slightly cooler. If real, such a bias would tend to make the free-convective coefficients for heat and moisture somewhat larger. Surface pressure was obtained from a data channel that estimates the surface pressure from aircraft values of pressure, virtual temperature, and elevation. Since the aircraft potential temperatures were defined with respect to a reference pressure of 100 kPa, Tskin was also converted to a potential temperature at this reference pressure using the surface pressure estimate. Saturation mixing ratio at the surface was obtained from Tskin and surface pressure using Teten’s formula (Stull 1988). As in Fairall et al. (1996) this saturation mixing ratio is multiplied by 0.98 to account for the reduction of vapor pressure caused by a typical salinity of 34 parts per thousand. Mean mixed-layer values of wind speed, potential temperature, and mixing ratio were estimated from mean values at the flight elevation. The aircraft offset table from the TOGA COARE Flux Group (1996) was used to correct observational offsets of the Electra aircraft instrumentation. GPS–corrected winds (D. Raymond 1997, personal communication) were used to estimate wind speeds. Several pairs of collinear flight lines flown within approximately a 0.6-h time period at different elevations from three of the flights showed that potential temperature and mixing ratio remained nearly constant between approximately 35- and 95-m elevation, although wind speed generally increased slightly with altitude. The vertical profiles used to estimate mixed-layer thicknesses also indicated fairly constant values of these parameters through the mixed layer. Table 1 lists the observed surface values of temperature Tskin , potential temperature skin , pressure Psurf , and mixing ratio rsat . Also listed are mixed-layer values of thickness z i , wind speed MML , potential temperature ML , and mixing ratio rML , determined for each flight-leg segment. Also in Table 1 are the surface fluxes of heat, moisture, and momentum. Surface fluxes of heat, moisture, and momentum were assumed to be reasonably represented by their respective values measured from 2206 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 56 TABLE 2. Mixed-layer and buoyancy scales. Variables are , ratio of virtual temperature flux to heat flux, both at the surface; ⌬ B , buoyancy temperature difference; ⌬ , virtual potential temperature difference between the surface skin and the mixed layer ; w*, Deardorff velocity; wB, buoyancy velocity; R*, mixed-layer Richardson number; C*D , mixed-layer transport coefficient for momentum; C*H , mixed-layer transport coefficient for heat; C*E, mixed-layer transport coefficient for moisture. Asterisks in the column labeled ‘‘Free conv.’’ indicate those data for which R* Ͼ 3.0, which were used to calculate the boldface mean values at the bottom of the table. Convective transport coefficients for momentum, heat, and moisture are indicated by bD, bH, and bE, respectively. the lowest aircraft flight level, because flux profile slopes were expected to be relatively small in the COARE region. The momentum flux components were combined to yield the square of the friction velocity [u*2 ϭ (uЈwЈ 2 ϩ ЈwЈ 2 )1/2 ]. The surface virtual temperature flux wЈЈ s was computed from the heat and moisture fluxes (Stull 1994), and this in turn was used to compute the Deardorff convective velocity w* ϭ [(g/T )z iwЈЈ s ]1/3 . Table 2 shows the following mixedlayer and buoyancy scales: ratio of kinematic virtual temperature flux to heat flux ; buoyancy temperature difference ⌬ B [Eq. (8)]; virtual temperature difference between the surface skin and the mixed layer ⌬ ; buoyancy velocity scale w B [Eq. (7)]; mixed-layer Richardson number R* [Eq. (11)]; and the mixed-layer transport coefficients for momentum C* D , heat C* H , and moisture C* E , and the corresponding convective transport coefficients (b D , b H , b E ). 5. Results Figure 2 shows the recorded potential temperature, mixing ratio, wind speed, and direction (vs time in seconds from 0000 UTC on 16 December 1992) for leg 1 1 JULY 1999 GREISCHAR AND STULL FIG. 3. Variation of mixed-layer transport coefficients for (a) heat C H , (b) moisture C E , and (c) momentum C D , with mixed-layer * * Richardson number R* . Values for the 30 COARE flight segments studied are plotted as *black squares, and for comparison the BLX83 land-based data are shown by diamonds. Approximate free convection occurs for R Ͼ 3.0, as indicated by the relatively constant values * transport coefficients for moisture and momentum. of mixed-layer Thin horizontal lines (dashed for BLX83) show mean values for the transport coefficients: C H ϭ 0.0061 (0.0063) and C E ϭ 0.0044 * (0.0063) for heat and moisture, respectively; and C *D ϭ 0.0053 * (0.0230) for momentum. (a) Shows linear fits to the semilog graphs for the COARE data (solid) and BLX83 data (dashed). (b) and (c) Thicker curves show second-order polynomial fits [in terms of log(R )] to the semilog graphs of the COARE data. * of flight RF16. A cool gust front similar to one seen in ship data (Tsukamoto and Ishida 1995) occurs near 95 200 sec. It is very evident in the temperature and wind data, and appears as an increase in variance of the mixing ratio observations. The liquid water channel on the Electra aircraft did not indicate any precipitation during this flight leg. The flight segments evaluated from this leg were chosen to avoid the discontinuities at the 2207 boundary of the gust front. Table 1 shows the changes in fluxes associated with this event (between RF16 1A and 1B on the warm side and 1C on the cool side). Consistent with the theory, heat, moisture, and momentum fluxes increased by approximately 100%, 23% and 50%, respectively, with the cooling, drying, and increased winds of the mixed layer. Figure 3 shows the variation of mixed-layer transport coefficients for heat, moisture, and momentum with the mixed-layer Richardson number. The land-based values from BLX83 are also plotted for comparison. As in Stull (1994) R* Ͼ 3.0 was the criterion adopted for free convection, indicating that buoyant production of turbulence kinetic energy is at least three times greater than shear production. As can be seen from both the BLX83 and the COARE data in Fig. 3, the choice of R* Ͼ 3 to define free convection is somewhat arbitrary, with the empirical coefficients becoming even more uniform as R* increases to 10 and beyond. Also plotted (Figs. 3b and 3c) are second-order polynomial fits [in terms of log(R*)], showing that the coefficients for moisture and momentum become approximately steady for free convection. The variation of C* H with R* shows a nearly linear variation (Fig. 3a) over the range of the semilog plot. This was also found over land (Stull 1994), which suggests that vertical heat transport continues to become more efficient as winds become lighter, perhaps caused by better organization of the coherent thermal structures. Figure 3a shows the remarkably similar linear trends for the COARE and BLX83 data with a fit to the combined datasets yielding the relation C* H ϭ 0.0091 Ϫ 0.0020 log(R*). Table 2 shows the results of averaging the 12 freeconvective cases for the proposed empirical values C* D ϭ 0.0053 (Ϯ0.0023), C* H (heat) ϭ 0.0061 (Ϯ0.0024), and C* E (moisture) ϭ 0.0044 (Ϯ0.0010). These values are plotted as horizontal solid lines in Fig. 3 and can be compared to the land-based BLX83 values (dashed lines). The very good agreement between the COARE and BLX83 values for C* H is evident in Fig. 3a. Figure 3b also has plotted the BLX83 values of C* H in order to test Stull’s (1994) assumption that C* E ഠ C* H . Although C* E is smaller than the BLX83 (and COARE) values of C* H , it is still within error limits. However, C* D (Fig. 3c) is only about 25% of the published value for nonmarine conditions. There is considerable scatter in the momentum data, which may have several causes. Some might be related to mesoscale heterogeneity caused by precipitation-induced intrusions of cool air masses (Williams et al. 1996). In addition, some may be due to insufficient leg length used to compute the momentum fluxes, to uncertainties in the mean mixedlayer wind speed and surface currents; but the difference is large enough to suggest that roughness-length effects or other factors may need to be incorporated into the CTT. The free and mixed formulations of CTT were tested. Figure 4 shows plots of predicted versus observed fluxes 2208 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 56 of heat, moisture, and momentum. Two predictions are presented: one based on free convection [(1)–(3)] for observations with R* Ͼ 3.0, and the other based on mixed convection [(13)–(15)] for all observations. Correlation coefficients and root-mean-square differences (mixed case) between the predicted and observed values shown in Fig. 4 give an estimate of performance. The root-mean-square error differences for the sensible and latent fluxes correspond to 2 W mϪ2 and 16 W mϪ2 , compared to typical magnitudes of those fluxes of 7.5 W mϪ2 and 94.3 W mϪ2 , respectively. Transport coefficients for the free-convective case have already been presented. The least squares best-fit coefficients for the mixed-convective case based on all the data in Tables 1 and 2 are for momentum fluxes C DML ϭ 0.000 38 (Ϯ0.000 12), b DML ϭ 0.000 41 (Ϯ0.000 09); for sensible heat fluxes C HML ϭ 0.000 44 (Ϯ0.000 05), b HML ϭ 0.000 30 (Ϯ0.000 04); and for latent heat fluxes C EML ϭ 0.000 62 (Ϯ0.000 19), b EML ϭ 0.000 18 (Ϯ0.000 03). The error estimates were obtained by a Monte Carlo simulation (Press et al. 1986) and allowed for random errors in observed values of up to 20% of their respective means. When compared to the previously published values of C DML ϭ 0.003 50, b DML 0.000 70, and C HML ϭ 0.001 00, b HML ϭ 0.000 25 there are considerable differences between most of the corresponding coefficients. The marine values for the mixed momentum coefficients are similar to those for mixed heat and moisture, which is consistent with the result for the convective case with C* D being more nearly equal to C* H and C* E for these marine data. Acoustic Doppler current profiler (ADCP) data from several of the ships in the IFA during the flights studied were used in an attempt to correct for the effects of surface currents on wind speed difference between the middle mixed layer and surface. These data give 0.5-h bin-averaged values centered at approximately 20-m depth and generally indicated eastward-moving currents with speeds up to 0.6 m sϪ1 but with considerable variability in space and time. These speeds and the general direction are born out by the many surface buoy trajectories plotted for December 1992 (TCIPO 1993) with the January 1993 trajectories showing more of a southward component. Since it was not possible to get an accurate estimate of the surface currents beneath each flight segment the transfer coefficients were recomputed for a Ϯ0.6 m sϪ1 adjustment to the wind speed. Table 3 shows the results of these calculations. The momentum coefficient C* D for the free-convective case is affected most since its computation involved the low- FIG. 4. Predicted and observed fluxes are compared for (a) heat, (b) moisture, and (c) momentum. Two predictions are presented: one based on free convection [(1)–(3)] with only free-convective cases plotted, and the other based on mixed convection [(13)–(15)] with all cases plotted. Coefficients determined in this study were used. ← Correlation coefficients (R f and R m for the free and mixed cases) between predicted and observed values are shown. The root-meansquare error for the mixed case (rms m ) is shown. The diagonal line representing a 1:1 perfect fit is shown for reference. 1 JULY 1999 2209 GREISCHAR AND STULL TABLE 3. Variation of CTT coefficients for a Ϯ0.6 m sϪ1 range of ocean surface currents, which could have affected the wind-speed difference between the mid-ML and the surface. Free convective Current speed ϩ0.6 m s 0 Ϫ0.6 m sϪ1 Ϫ1 C* D C* H C* E bD bH bE 0.0039 0.0053 0.0108 0.0059 0.0061 0.0065 0.0043 0.0044 0.0047 0.00029 0.00041 0.00081 0.00047 0.00050 0.00054 0.00033 0.00035 0.00038 C DML C HML C EML b DML b HML b EML 0.00043 0.00038 0.00028 0.00044 0.00044 0.00043 0.00061 0.00062 0.00062 0.00026 0.00041 0.00060 0.00026 0.00030 0.00033 0.00013 0.00018 0.00023 Mixed convective Current speed ϩ0.6 m s 0 Ϫ0.6 m sϪ1 Ϫ1 est wind speeds from 0.8–3.1 m sϪ1 . The small changes in the heat and moisture coefficients are due to fewer or more points falling into the free-convective range. All the mixed convective case coefficients are affected since they involve wind speed, though again the momentum coefficients are most affected. A comparison with the latest (version 2.0f ) COARE bulk parameterization (Fairall et al. 1996), which was designed for light wind conditions, was also made. Figure 5 shows the comparison with correlations and rms differences between the bulk and covariance (eddy correlation) fluxes for the 25 flight segments flown at elevations of less than one-tenth of the mixed-layer thickness. At larger flux values the bulk estimates are somewhat greater than the measured covariance values. The 10-m wind speeds for these data range from 0.8 to 7.8 m sϪ1 , so they represent a fair range of values for this region. The correlation coefficients are slightly better than those for CTT shown in Fig. 4 though the rms errors are larger. Confining the sample to cases with wind speeds of 5 m sϪ1 or less lowered the rms errors to 0.0017 K m sϪ1 , 0.0055 g kgϪ1 m sϪ1 , and 0.0033 m sϪ1 , respectively, for the heat, moisture, and momentum fluxes, values comparable to CTT. This is a fairly crude comparison since the aircraft was measuring at heights from 30–66 m and there is the possibility of some deviation from the surface values. Using completely independent data (eddy correlation flux observations from the R/V Moana Wave), Chang and Grossman (1999) compared CTT with bulk surface flux algorithms. For their calculations z i was assumed to be comparable to the lifting condensation level (LCL), and mid-mixed-layer values of other variables were estimated by adjustment from observation level to the top of the surface layer (ϳ0.1z i ) using Monin–Obukhov similarity theory and a linear adjustment from there to the middle of the mixed layer (0.5z i ). Based on the Moana Wave flux observations they determined an ‘‘ocean’’ set of CTT coefficients and concluded that CTT was the best algorithm for calculating latent heat, sensible heat, and momentum fluxes. The relationship between Deardorff velocity and buoyancy velocity scales given in (16) was tested. Figure 6 is a scatterplot of the two velocity scales showing the strong relationship of their magnitudes as evidenced by the least squares line passing through the origin with a slope of 1.0033. The corresponding correlation coefficient is 0.41 with an rms error of approximately 12% of the mean. The low correlation coefficient reflects the limited range of w* values for these TOGA COARE flights. Figure 7 shows a plot of buoyant versus virtual-temperature differences, which were shown to be approximately equal for convection over land. The least squares line through the origin has a slope of 0.954, a corresponding correlation coefficient of 0.82, and an rms error of approximately 14% of the mean. Thus, the approximation in (9) is also good for these more humid marine conditions. A test of the prediction of z i given in (12) was made with the limited number of flight legs adjacent to slant vertical profiles. The predicted versus observed values of z i showed considerable scatter. Since the predicted value of z i in (12) depends on the square of the observed heat flux and on ⌬Ϫ2 , and ⌬ under marine conditions is generally quite small, there is a strong dependence on the observed heat flux and on accurate temperature measurements both at the surface and mid–mixed layer. We do not recommend use of this relationship in the COARE region. The large impact of measurement errors on reducing skill to diagnose mixed-layer depths was anticipated by Stull (1994). 6. Conclusions Turbulent flux measurements made from five flights of the Electra aircraft during the TOGA COARE IOP confirmed the CTT free-convective coefficient for sensible heat flux. Averaging 12 estimates of the dimensionless mixed-layer transport coefficients gave values for C* H and C* E for heat and moisture fluxes that were not significantly different from each other and from the previously published value. However, the estimate of C* D from momentum fluxes is less than one-quarter of the published land-based value. There is considerable 2210 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 56 FIG. 6. Comparison of Deardorff velocity and buoyancy velocity scales described in (16). Slope of the least squares line passing through the origin, the corresponding correlation coefficient, and the root-mean-square error are shown. The diagonal line representing a 1:1 perfect fit is shown for reference. scatter in the momentum data that may be due to insufficient flight-leg length used to compute the momentum fluxes, and to uncertainties in the mean mixedlayer wind speed and surface currents; but the difference is large enough to suggest that roughness-length effects or other factors need to be incorporated into the CTT. FIG. 5. Comparison of fluxes predicted by the COARE bulk parameterization vs those observed by the Electra aircraft for (a) sensible heat, (b) moisture, and (c) momentum. Correlation coefficients between predicted and observed values are shown. The root-meansquare errors are listed with the value for the cases with wind speed less than 5 m sϪ1 shown in parentheses. The diagonal line representing a 1:1 perfect fit is shown for reference. FIG. 7. Comparison of ⌬ B and ⌬ V described in (10). Slope of the least squares line passing through the origin, the corresponding correlation coefficient, and the root-mean-square error are shown. The diagonal line representing a 1:1 perfect fit is shown for reference. 1 JULY 1999 2211 GREISCHAR AND STULL The mixed formulations of CTT, which include a wind speed factor, gave good results for heat, moisture, and momentum fluxes even under moderate wind conditions of 5–8 m sϪ1 ; however, the best-fit coefficients were considerably different than the published land-based values. A comparison of the observed fluxes with estimates computed from the COARE bulk algorithm was made. Heat, moisture, and momentum fluxes computed using this bulk algorithm, which was designed for light wind conditions, tended to be larger than the observed fluxes for wind speeds greater than 5 m sϪ1 . At wind speeds less than 5 m sϪ1 the agreement improved considerably, with rms errors comparable to those of the mixed CTT model. Various relations resulting from the CTT were tested for marine conditions. The theoretical value of the ratio of the Deardorff velocity and buoyancy velocity scales was confirmed. The approximate equality of the buoyancy and virtual temperature differences was found to hold even in the high-humidity conditions of the equatorial Pacific. Estimation of the mixed-layer thickness from surface heat flux and temperature difference between surface and mixed-layer showed considerable scatter, thus limiting its usefulness in the tropical marine context of the COARE region. The feasibility of extending CTT to marine conditions using the mixed formulation with adjusted coefficients has been demonstrated. Root-mean-square errors of roughly 22% (2 W mϪ2 and 16 W mϪ2 for sensible and latent heat fluxes, respectively) were obtained for this mixed-convective CTT model tuned for marine conditions. More marine studies are needed to reveal the roughness and possibly other effects on CTT, and whether this model can be used in other areas of the ocean. Acknowledgments. This work was supported mostly by National Science Foundation Grant ATM-9411467, with partial support from Canadian grants from the Atmospheric Environment Service and the Natural Sciences and Engineering Research Council. The TOGA COARE International Project Office was very helpful in providing the Electra data (processed at the NCAR ATD) and necessary documentation. Discussions with John Young, Bob Grossman, and Yolande Serra (who provided corrected radiometric SST data) are greatly appreciated. Comments from three anonymous reviewers greatly improved the manuscript. REFERENCES Asencio, N., J. P. Lafore, P. Pires, and J. L. Redelsperger, 1993: Analyses quotidiennes du Cepmmt durant la periode d’observations intensives de l’experience TOGA/COARE. METEO FRANCE Groupe de Meteorologie a Moyenne Echelle Note de Travail 12, 243 pp. [Available from La Librairie de Meteo-France, 2 Avenue Rapp, 75340 Paris Cedex 07 France.] Bond, G., and D. Alexander, 1994: TOGA COARE Meteorology Atlas. TOGA COARE International Project Office, 366 pp. Chang, H.-R., and R. L. Grossman, 1999: Evaluation of bulk surface flux algorithms for light wind conditions using data from the Coupled Atmosphere-Ocean Response Experiment (COARE). Quart. J. Roy. Meteor. Soc., in press. Fairall, C. W., E. F. Bradley, D. P. Rogers, J. B. Edson, and G. S. Young, 1996: Bulk parameterization of air-sea fluxes for Tropical Ocean-Global Atmosphere Coupled Ocean–Atmosphere Response Experiment. J. Geophys. Res., 101, 3747–3764. Kustas, W. P., T. J. Schmugge, and L. E. Hipps, 1996: On using mixedlayer transport parameterizations with radiometric surface skin temperature for computing regional scale sensible heat flux. Bound.-Layer Meteor., 80, 205–221. Miller, E. R., and R. B. Friesen, 1989: Standard output products from the NCAR Research Aviation Facility. NCAR Research Aviation Facility Bulletin 9, 70 pp. Otles, Z., and J. A. Young, 1996: Influence of shallow cumuli on subcloud turbulence fluxes analyzed from aircraft data. J. Atmos. Sci., 53, 665–676. Press, W., B. Flannery, S. Teukolsky, and W. Vetterling, 1986: Numerical Recipes. Cambridge University Press, 818 pp. Stull, R., 1988: An Introduction to Boundary Layer Meteorology. Kluwer Academic Publishers, 666 pp. , 1994: A convective transport theory for surface fluxes. J. Atmos. Sci., 51, 3–22. , 1997: Reply. J. Atmos. Sci., 54, 579. Sun, J., J. F. Howell, S. K. Esbensen, L. Mahrt, C. M. Greb, R. Grossman, and M. A. LeMone, 1996: Scale dependence of air– sea fluxes over the western equatorial Pacific. J. Atmos. Sci., 53, 2997–3012. TCIPO, 1992: TOGA COARE Operations Plan. TOGA COARE International Project Office, 391 pp. [Available from TCIPO UCAR, P.O. Box 3000, Boulder, CO 80307.] , 1993: TOGA COARE Intensive Operations Period Operations Summary. TOGA COARE International Project Office, 506 pp. [Available from TCIPO UCAR, P.O. Box 3000, Boulder, CO 80307.] TOGA COARE Flux Group, cited 1996: TOGA COARE Aircraft Offset Table. [Available online at http://www.wave.eng.uci.edu/ Projects/TogaCepex/togacoare.html.] Tsukamoto, O., and H. Ishida, 1995: Turbulent flux measurements and energy budget analysis over the eqauatorial Pacific during TOGA-COARE IOP. J. Meteor. Soc. Japan, 73, 557–568. Webster, P. J., and R. Lukas, 1992: TOGA COARE: The Coupled Ocean–Atmosphere Response Experiment. Bull. Amer. Meteor. Soc., 73, 1377–1416. Williams, A. G., H. Kraus, and J. M. Hacker, 1996: Transport processes in the tropical warm pool boundary layer. Part I: Spectral composition of fluxes. J. Atmos. Sci., 53, 1187–1202.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Faculty Research and Publications /
- Convective Transport Theory for Surface Fluxes Tested...
Open Collections
UBC Faculty Research and Publications
Convective Transport Theory for Surface Fluxes Tested over the Western Pacific Warm Pool. Greischar, Lawrence; Stull, Roland B. 1999-07-31
pdf
Page Metadata
Item Metadata
Title | Convective Transport Theory for Surface Fluxes Tested over the Western Pacific Warm Pool. |
Creator |
Greischar, Lawrence Stull, Roland B. |
Publisher | American Meteorological Society |
Date Issued | 1999-07 |
Description | Turbulent flux measurements from five flights of the National Center for Atmospheric Research Electra aircraft during the Tropical Oceans and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE) are used to test convective transport theory (CTT) for a marine boundary layer. Flights during light to moderate winds and under the clearest sky conditions available were chosen. Fluxes of heat, moisture, and momentum were observed by the eddy-correlation method. Mean kinematic values for the observed sensible and latent heat fluxes and momentum flux were 0.0061 K m s−1, 0.0313 g kg−1 m s−1, and 0.0195 m2 s−2, respectively. For the range of mixed-layer wind speeds (0.8–8.4 m s−1) studied here, the version of CTT that includes the mixed effects of buoyant and shear-driven transport give a better fit to the observations than either the COARE bulk algorithm or the pure free-convection version of CTT. This is to be expected because both of those latter parameterizations were designed for light winds (<5 m s−1 approximately). The CTT empirical coefficients listed in Table 3 exhibited slight sensitivity to the COARE light flux conditions, compared to their previous estimates during larger fluxes over land. For example, COARE heat fluxes were roughly 10 times smaller than previous land-based flux measurements used to calculate CTT coefficients, but the corresponding empirical mixed-layer transport coefficients were only 3% smaller. COARE momentum fluxes were also roughly 10 times smaller, but the CTT coefficients were about four times smaller. The greater variation in momentum coefficient may be due, in part, to insufficient flight-leg length used to compute momentum fluxes, to uncertainties in the effects of the ocean surface current and waves, or perhaps to roughness differences. Copyright 1999 American Meteorological Society (AMS). Permission to use figures, tables, and brief excerpts from this work in scientific and educational works is hereby granted provided that the source is acknowledged. Any use of material in this work that is determined to be “fair use” under Section 107 of the U.S. Copyright Act or that satisfies the conditions specified in Section 108 of the U.S. Copyright Act (17 USC §108, as revised by P.L. 94-553) does not require the AMS’s permission. Republication, systematic reproduction, posting in electronic form, such as on a web site or in a searchable database, or other uses of this material, except as exempted by the above statement, requires written permission or a license from the AMS. Additional details are provided in the AMS Copyright Policy, available on the AMS Web site located at (http://www.ametsoc.org/) or from the AMS at 617-227-2425 or copyright@ametsoc.org. |
Genre |
Article |
Type |
Text |
Language | eng |
Date Available | 2011-04-11 |
Provider | Vancouver : University of British Columbia Library |
DOI | 10.14288/1.0041842 |
URI | http://hdl.handle.net/2429/33499 |
Affiliation |
Science, Faculty of Earth and Ocean Sciences, Department of |
Citation | Greischar, Lawrence, Stull, Roland B. 1999. Convective Transport Theory for Surface Fluxes Tested over the Western Pacific Warm Pool. Journal of the Atmospheric Sciences. 56(13) 2201-2211. http://journals.ametsoc.org/doi/pdf/10.1175/1520-0469%281999%29056%3C2201%3ACTTFSF%3E2.0.CO%3B2 |
Peer Review Status | Reviewed |
Scholarly Level | Faculty |
Copyright Holder | Stull, Roland B. |
Aggregated Source Repository | DSpace |
Download
- Media
- [if-you-see-this-DO-NOT-CLICK]
- Stull_AMS_1999_JAS2201.pdf [ 238.86kB ]
- [if-you-see-this-DO-NOT-CLICK]
- Metadata
- JSON: 1.0041842.json
- JSON-LD: 1.0041842+ld.json
- RDF/XML (Pretty): 1.0041842.xml
- RDF/JSON: 1.0041842+rdf.json
- Turtle: 1.0041842+rdf-turtle.txt
- N-Triples: 1.0041842+rdf-ntriples.txt
- Original Record: 1.0041842 +original-record.json
- Full Text
- 1.0041842.txt
- Citation
- 1.0041842.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Country | Views | Downloads |
---|---|---|
United States | 3 | 9 |
City | Views | Downloads |
---|---|---|
Ashburn | 3 | 0 |
{[{ mDataHeader[type] }]} | {[{ month[type] }]} | {[{ tData[type] }]} |
Share
Embed
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"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.52383.1-0041842/manifest