Influence of the Lake Erie Overlake Boundary Layer on Deep Convective Storm Evolution

Transcription

1 OCTOBER 2012 W O R K O F F E T A L Influence of the Lake Erie Overlake Boundary Layer on Deep Convective Storm Evolution THOMAS E. WORKOFF* Department of Atmospheric Sciences, and Center for Atmospheric Science, Illinois State Water Survey, Prairie Research Institute, University of Illinois at Urbana Champaign, Urbana, Illinois DAVID A. R. KRISTOVICH Center for Atmospheric Science, Illinois State Water Survey, Prairie Research Institute, University of Illinois at Urbana Champaign, Urbana, Illinois NEIL F. LAIRD Hobart and William Smith Colleges, Geneva, New York ROBERT LAPLANTE NOAA/National Weather Service Forecast Office, Cleveland, Ohio DANIEL LEINS NOAA/National Weather Service Forecast Office, Phoenix, Arizona (Manuscript received 1 July 2011, in final form 6 February 2012) ABSTRACT The influence that the overlake boundary layer has on storm intensity and structure is not well understood. To improve scientists understanding of the evolution of storms crossing Lake Erie, 111 events during were examined using observations from Weather Surveillance Radar-1988 Doppler (WSR-88D), surface, buoy, and rawinsonde sites. It was found that on average, all storm modes tended to weaken over the lake; however, considerable variability in changes of storm intensity existed, with some storms exhibiting steadystate or increasing intensity in specific environments. Noteworthy changes in the storm maximum reflectivity generally occurred within 60 min after storms crossed the upwind shoreline. Isolated and cluster storm modes exhibited much greater weakening than those storms organized into lines or convective complexes. The atmospheric parameters having the greatest influence on storm intensity over Lake Erie varied by mode. Isolated and cluster storms generally weakened more rapidly with increasingly cold overlake surface air temperatures. Linear and complex systems, on the other hand, tended to exhibit constant or increasing maximum reflectivity with cooler overlake surface air temperatures. It is suggested that strongly stable conditions near the lake surface limit the amount of boundary layer air ingested into storms in these cases. * Current affiliation: Systems Research Group, Inc., NOAA/ Hydrometeorological Prediction Center, Camp Springs, Maryland. Corresponding author address: David Kristovich, ISWS, PRI, University of Illinois at Urbana Champaign, 2204 Griffith Dr., Champaign, IL dkristo@illinois.edu 1. Introduction and background It is well known that on climatic time scales the Great Lakes suppress clouds and precipitation near the coastlines during the warm season, particularly in regions downwind of the lakes (e.g., Augustine et al. 1994; Scott and Huff 1996). However, operational weather forecasts for these regions must account for the rate of DOI: /WAF-D Ó 2012 American Meteorological Society

2 1280 W E A T H E R A N D F O R E C A S T I N G VOLUME 27 change of storm intensity as well as situations where the storm evolution differs from climatology. Few published studies have examined how the lakes and their associated overlake boundary layers (OLBLs) 1 affect preexisting deep convective storms, particularly during the warm season where the water surface is predominately cooler than the surrounding land. To develop better forecasting techniques for areas in western lower Michigan, Graham et al. (2004) constructed a 5-yr climatology of the behavior of severe mesoscale convective systems (MCSs) crossing Lake Michigan. They found that MCSs responded variably to passing over the water surface, depending on time of year and structure. For example, MCSs were more likely to weaken or dissipate in the summer months (June August, JJA) than in the spring or fall. In addition, the stability of the OLBL appeared to influence storm evolution over the lake. The study showed that the convection tended to maintain intensity when the OLBL was strongly stable (lake surface $2.58C colder than the overlake air temperature), despite the reduction in surface-based instability. However, their study did not examine the rate at which storm intensity changes occurred. Other studies have used numerical modeling techniques to examine the effect of surface-based stable layers on convective systems. Parker (2008) used an idealized model to simulate squall-line interactions with such layers. In his study, a surface-based isothermal cool layer, approximately 1 km deep, was introduced to a mature squall line. He found that storm evolution is associated with the relative buoyancies of the cool layer and the cold pool produced by the convective system. When the cold pool is denser, lifting along the leading edge of the cold pool can overcome a significant amount of surface-based convective inhibition (CIN) if convective available potential energy (CAPE) is present. This, in turn, allows updrafts to maintain their intensity despite interaction with convectively unfavorable near-surface conditions. The results of Parker (2008) also showed that in instances where the surface cool layer ahead of the storm is cooler than the cold pool air, storms can become elevated atop the surface layer, with continuous updrafts being driven by bore propagation. The effect of changing low-level conditions near shorelines on mature squall lines was examined in an idealized modeling study by Lericos et al. (2007). Their study focused on squall lines moving eastward from the Gulf 1 Overlake boundary layers are often identified in the scientific literature and forecast discussions as marine internal boundary layers. OLBL is used in this article to emphasize the more limited spatial extent of lakes relative to oceans. of Mexico onto the northern Florida coast. In cases of moderate to strong ambient low-level vertical wind shear, land breezes were found to reduce the magnitude of ambient shear. This, in turn, caused the updraft of the squall line to tilt in an upshear direction, resulting in weakening of the squall line. While the Graham et al. (2004), Parker (2008), and Lericos et al. (2007) studies provide important insights into the evolution of convection interacting with varying near-surface conditions in the vicinity of coastlines, they consider only organized convective systems (e.g., MCSs), leaving the potential effect of a large lake on other convective modes unknown. In addition, the rate of change in the storm system, critical in an area of widely varying lake sizes, was not examined. The current study uses a climatological approach incorporating observations commonly available to forecasters (e.g., radar, hourly surface observations, rawinsonde observations) to examine how large lakes influence the intensity of the wide range of convective modes. 2. Data and methodology Statistical analyses were performed by season, time of day, and storm mode to identify common storm responses and the potential causes. Weather Surveillance Radar-1988 Doppler (WSR-88D) data from Cleveland, Ohio (KCLE) (Fig. 1), were used for this study. The area of study chosen was south of N over Lake Erie; this was selected to include storms within a range of ;100 km of the radar when crossing the lake. This is an area where convective storms frequently pass over Lake Erie (Laing and Fritsch 1997) and are routinely monitored by the Cleveland National Weather Service Forecast Office. a. Identification of events Appropriate events (preexisting convective storms that advected from land to over Lake Erie) were identified by a three-step process: 1) surface observations were used to identify dates when significant rain occurred, 2) composite radar reflectivity imagery was used to remove events where the observed precipitation was clearly not associated with convective storms moving over Lake Erie, and 3) more detailed analyses of WSR-88D Level II data were conducted to identify which of the remaining events were appropriate for this study. Precipitation observations at sites around Lake Erie (Fig. 1, appendix A) were used to identify potentially appropriate events. The azimuth range of with respect to the location of KCLE was chosen to capture the most likely path of storms crossing the lake while staying inside the study area, despite the lack of suitable

3 OCTOBER 2012 W O R K O F F E T A L FIG. 1. Locations of observation sites used in this study. Surface observation sites over land are denoted by black diamonds, precipitation observation sites denoted by numbered circles, buoys denoted by black ovals, radar locations denoted by, and raob sounding sites by. Radar ranges for KCLE and Detroit, MI (KDTX), are denoted by a dashed circle (50 km) and a solid circle (100 km), respectively. Dashed line denotes N, the northern boundary of the study. Names and locations of precipitation observation sites are located in appendix A. precipitation observation stations on the southern and eastern shores of Lake Erie. Days when snow occurred or surface air temperatures were,08c were excluded. Over the period , a total of 504 dates were identified with $1.27 cm (0.5 in.) of precipitation and above freezing temperatures. National composite radar reflectivity imagery obtained from the National Center for Atmospheric Research (NCAR) was examined for each of the 504 dates to identify those with convective storms developing upwind and moving over Lake Erie. While the national composite lacks the desired spatial and temporal resolutions, it allows for identification of precipitation systems that are clearly not convective. Any date that had a maximum base reflectivity,35 dbz on the national composite was eliminated. These criteria reduced the list of dates with possible convective storms moving over the Lake Erie study region to 280. It should be noted that an individual date may have included multiple convective storm events. More detailed analysis of Level II radar data from KCLE was conducted for each of the remaining 280 potential events using GR2Analyst (Gibson Ridge, version 1.60). These data were obtained from the National Climatic Data Center (NCDC). In order for the event to qualify for this study, the base reflectivity observed at the 0.58 elevation level with values.45 dbz, a common threshold for convective rainfall, needed to be observed during at least three consecutive radar scans ( min) prior to moving over the lake. If a reflectivity value of $45 dbz was not observed at least 15 min before the storm crossed the shoreline, the event was not included in the study. This was done to ensure storms included in the study were not experiencing large changes in intensity as they approached the lake; therefore, the evolution of intensity observed over the water was likely the result of OLBL interaction. Note that this criterion would remove cases where the reflectivity increased to.45 dbz over the lake; however, this was an infrequent occurrence and is not thought to change the findings of this study. In some instances, areas of.45-dbz reflectivity were embedded in areas of widespread rainfall associated with synoptic disturbances. Since the presence and

4 1282 W E A T H E R A N D F O R E C A S T I N G VOLUME 27 FIG. 2. Base reflectivity images from the Cleveland WSR-88D radar showing examples of (a) isolated, (b) cluster, (c) complex, and (d) linear storm modes used in this study. organization of convective storms is not always obvious in such situations, these events were also eliminated from the database. Finally, events when the convective storm moved north of the N boundary were removed. This resulted in a total of 111 events (on 101 days) used for analysis of storm evolution over Lake Erie. b. Determination of convective mode Each storm was classified by mode to investigate its role in convective evolution over the lake. During events when the storm mode underwent a change over the lake, the storm was classified by the mode at the time it moved over the upwind shoreline. Previous studies (e.g., Fowle and Roebber 2003; Gallus et al. 2008) tended to use three dominate convective modes (linear, isolated, multicellular); a fourth mode (complex) was added for this study due to the propensity for nonlinear mesoscale convective complex (MCC) development in the study area, as well as to allow for a more direct comparison with the findings from Graham et al. (2004). Examples of the four convective storm modes used in this study are shown in Fig. 2 and were defined as follows: Isolated (I) individual convective storms, including supercells, with a continuous area of reflectivity.45 dbz that was,500 km 2 (Fowle and Roebber 2003); storms must have been separated from other reflectivity maxima by greater than two diameters of the area containing reflectivity values of.45 dbz (Fig. 2a); Cluster (CL) areas of unorganized convection, with several (three or more) reflectivity maxima located within a distance equivalent of two diameters of the 45-dBZ reflectivity area of each individual storm (Fig. 2b); areas of.45-dbz reflectivity were generally small (,40 km 2 ) for individual storms and were separated by reflectivities,35 dbz; Complex (C) nonlinear, organized storm structure having an area of.500 km 2 with continuous reflectivity values of.45 dbz (Fig. 2c); and Linear (L) an area with reflectivity values.45 dbz organized in a curvilinear manner (Fig. 2d); storms were considered linear if they were organized in a line,50 km wide, exhibited a length width ratio of at least 3:1 (Fowle and Roebber 2003), and areas of reflectivity.45 dbz were separated by less than two of their diameters. c. Analysis of storm evolution Radar reflectivity changes were used as an indication of storm intensity variations. The maximum base reflectivity, recorded to the nearest whole decibel, at the

5 OCTOBER 2012 W O R K O F F E T A L FIG. 3. Number of occurrences (111 total) of each storm mode by month elevation angle was determined from KCLE Level II data for each volume scan (approximate time interval of 5 min). Maximum reflectivity was recorded from 30 min before (230 min) the time the storm moved over the water (TMOW) to at least 60 min following the TMOW (160 min), or until the storm arrived at the downwind shoreline. In events when the storm had not developed at 230 min, the maximum reflectivity was determined from the volume scan at the time when the maximum reflectivity exceeded 35 dbz prior to TMOW. The storm evolution over Lake Erie was determined by comparing the maximum reflectivity at TMOW with both the maximum reflectivity at 30 (130 min) and at 60 min (160 min) over the lake. For simplification of discussion, five categories were used to describe the evolution of an individual storm: major strengthening, with $ 18 dbz maximum reflectivity change, moderate strengthening, with 13 to 17 dbz maximum reflectivity change, no change, with 2 to 12 dbz maximum reflectivity change, moderate weakening, with 23 to 27 dbz maximum reflectivity change, and major weakening, with # 28 dbz maximum reflectivity change. d. Environmental and storm parameters To examine the influence of environmental conditions on storm evolution over Lake Erie, several parameters were developed to represent the atmospheric conditions that storms encountered. These parameters and their associated definitions are given in appendix B. The data used to calculate these parameters for each storm are described below. Observations of temperature, wind speed, and wind direction collected at surface sites over land near the upwind (downwind) shore of Lake Erie were used to represent the upwind (downwind) conditions encountered by each storm (Fig. 1). The upwind and downwind observation sites were relative to individual storm motion as determined from radar data. The majority (about 75%) of the storms examined had mean motions in the range from southwest (2258) through northwest (3158C). Observation sites used for each storm were those deemed closest to the storm path. Only observations taken before the storm appeared to influence conditions at each site were used for this purpose. Overlake conditions were derived from air and water temperature data collected at National Oceanic Atmospheric Administration (NOAA) buoy (Fig. 1) for each storm when the buoy was operational (generally April November). All buoy data were acquired from the National Data Buoy Center. Storms for which buoy data were not available were excluded from any analyses that included overlake conditions. Velocity azimuth display (VAD) wind analyses provided in the KCLE WSR-88D Level III dataset were used to represent surface to 3-km and surface to 6-km wind velocities. 3. Results a. Temporal distribution of storms The largest frequency (76%) of storms moving over Lake Erie occurred in the months of May August, with the maximum in July (Fig. 3). Isolated, complex, and linear storms showed similar monthly trends while cluster storms exhibited a peak in May. As expected, storms were more frequent during the afternoon and evening hours (Fig. 4) with differences in timing

6 1284 W E A T H E R A N D F O R E C A S T I N G VOLUME 27 FIG. 4. Same as Fig. 3 except by time of day (using TMOW). Time of day shown in local standard time (LST 5 UTC 2 4/5 h). between modes. Cluster and isolated storms had a relative maximum during the late afternoon and evening [ local standard time (LST)], indicative of the importance of surface buoyancy to their development. Linear storms showed a broad maximum during the late morning to late evening ( LST) and minimum near sunrise ( LST). The relative maximum for complex storms occurred during the overnight ( LST) and morning ( LST) hours, similar to the pattern of initial convective development in the late evening and maximum extent in the overnight hours observed by Ashley et al. (2003). This is consistent with the findings of Graham et al. (2004) and Parker (2008) that MCS events can maintain their evolution over shallow stable layers. b. Evolution of storms over Lake Erie The evolution of storm intensity, estimated by changes in the maximum base reflectivity, as a function of time over Lake Erie and storm mode are shown in Table 1. Regardless of storm mode, increases of $3dBZ are rare (5% of events) and no storms exhibited major strengthening ($8 dbz). One of the most notable features is the difference in storm evolution at 130 min after TMOW and 160 min after TMOW (Table 1). At 130 min, the majority of storms experienced no change in intensity (22/12 dbz), with smaller percentages of all storms experiencing moderate weakening (23/27 dbz) and significant weakening (#28 dbz). This was consistent for all storm modes. On average, the maximum base reflectivity only changed by 21 to22 dbz from the time the storm crossed the upwind shore. In contrast, excluding storms that moved onshore, by 160 min only about a third of all storms experienced no change in reflectivity and 28% of the storms experienced major weakening. Average reflectivity changes ranged from 22 to29 dbz at 160 min from TMOW. Differences in storm evolution over Lake Erie were noted between storm modes at 160 min (Table 1). Linear and complex storms tended to weaken less than isolated and cluster storms. For example, average reflectivity changes were only 22 to23 dbz for linear and complex storms, while for isolated and cluster storms the average reflectivity changes were about 27 to29 dbz. Additionally, only 10% 15% of linear and complex storms experienced major weakening at 160 min, compared to 40% of isolated and cluster storms. c. Relationships with environmental parameters Graham et al. (2004), Lericos et al. (2007), and Parker (2008) found that changes in environmental conditions associated with a land water boundary can play important roles in subsequent storm evolution. While some factors affecting the role of the boundary slowly change over long time scales (e.g., surface friction changes over land due to plant growth), forecasting applications focus on shorter time-scale changes due to variations in atmospheric conditions. This study focuses on how short-term changes in atmospheric conditions influence storm evolution over Lake Erie. To investigate which parameters or combinations of parameters are most closely related to storm evolution, stepwise linear multivariable regression (SLMR) analysis was conducted. SLMR selects only the most statistically significant parameters for inclusion in the linear regression model, which makes it useful for identifying parameters that exhibit the largest influence on storm evolution. The regression model takes the form DdBZ 5 C 1 b o V o 1 b 1 V 1 1,

7 OCTOBER 2012 W O R K O F F E T A L TABLE 1. Percentage of events of each storm mode (rows) that experienced associated changes in maximum reflectivity at (top) 30 and (bottom) 60 min over water. The bold italic values denote the highest percentage for each mode. The onshore column denotes percentage of events that moved over the downwind shoreline. The numbers in parentheses below the convective mode represent the number of events in each category. Far-right column shows the average maximum reflectivity change of each storm mode over water. #28 dbz 23/27 dbz 12/22 dbz 13/17 dbz $18 dbz Onshore Avg DdBZ 130 min Total (111) Cluster (30) Isolated (25) Linear (34) Complex (22) min Total (111) Cluster (30) Isolated (25) Linear (34) Complex (22) where DdBZ is the change in reflectivity at 160 min, C is a constant, V i are the predictor parameters, and b i are the coefficients of V i. Only parameters that have a statistical significance level (a)ofa # 0.10 are retained. In events where no parameters met the statistical significance criteria, the parameter with the highest statistical significance is shown. The SLMR used 160-min dbz change as the predictand and five predictor parameters: 1) the nearsurface air temperature over the lake (LAT); 2) OLBL stability, estimated by the difference between the water temperature and the overlake near-surface air temperature (LT-LAT, negative values imply stable conditions); 3) spatial variations in near-surface conditions, estimated by the difference between the overlake nearsurface air temperature and the upwind (relative to storm motion) overland near-surface temperature (LAT-UWT); 4) larger-scale spatial variations in nearsurface conditions not associated with the lake, estimated by the difference between the downwind and upwind near-surface equivalent potential temperature (u e ); and 5) 3-km VAD wind velocities, which may also give some indication of the ambient low-level vertical wind shear. Changes in surface temperature, low-level stability, and low-level wind speed were included, as they are factors generally cited as influencing convective structure and evolution. Table 2 shows which atmospheric parameters had the most influence on storm evolution determined by the SLMR analysis: overlake air temperature (LAT) for cluster and complex storms, horizontal temperature differences (LAT-UWT) for isolated storms, and 3-km wind speed for linear storms. The results suggest that near-surface temperature changes associated with the OLBL are the most influential parameters for all storm modes, with the exception of linear. However, the differing signs of the b 1 coefficients suggest the effect of lake-induced temperature differences depends on storm mode. The positive b 1 found for cluster and isolated modes shows a direct relationship between reduced overlake air temperature and storm intensity decrease. However, the negative b 1 value for complex storms shows that cooler overlake temperatures are associated with increasing storm intensity. It should be noted that very few events exhibited environmental or storm evolution changes far removed from all of the other events. Such events can have large influences on the relationships determined by SMLR. If statistical outliers (defined as greater than three standard deviations) are removed, the same parameters were chosen by SMLR. However, the significance of the regression model was greatly improved for the isolated (a ) mode. Another method of illustrating the effect of multiple parameters on storm evolution over the lake is flowchart analysis. Flowcharts (Fig. 5) were created to demonstrate how maximum reflectivity changed, at 60 min over the water, with regard to combinations of two parameters: overland lake temperature difference (LAT- UWT1/2) and overlake air temperature (LAT greater than or less than 228C). The value of 228C for the overlake air temperature was chosen to approximately denote transitions between cool and warm seasons. For this analysis, cluster and isolated storms were combined into an unorganized mode and linear and complex

8 1286 W E A T H E R A N D F O R E C A S T I N G VOLUME 27 TABLE 2. Summary of the results from SLMR using atmospheric parameters to predict storm intensity changes over Lake Erie (at 60 min over the water), for each storm mode. Columns 1 3 (italicized) display information about the linear regression model developed with the selected parameters. Columns 4 and 5 display information about the chosen parameters inside the model. Here, N is the number of storms included in the regression, R 2 represents the coefficient of determination (goodness of fit) of the regression model, and significance displays the statistical significance level (a) of the model. Parameters N R 2 b coefficient (standardized) Significance Mode: Cluster A. LAT Mode: Isolated A. LAT-UWT Mode: Linear A. 3-km speed Mode: Complex A. LAT storms were combined as organized modes due to their similar evolutionary patterns while traversing the lake (Table 1). Each row of the flowchart represents a way of categorizing the storms: 1) convective mode, 2) overland lake temperature difference, and 3) overlake air temperature. For example, the first two rows in Fig. 5a illustrate that cluster and isolated storms were divided into 7 events with conditions over the lake warmer than those over the upwindshorearea(lat-uwt. 0) and 36 events occurred when the overlake air temperature was cooler than at the upwind shore (LAT-UWT, 0). Each box in the flowchart gives information on the parameter used to categorize the events and the number of events (above the horizontal line) and convective trends (below the horizontal line). Convective trends are summarized by two categories: no change/moderate strengthening or major weakening (refer to section 2 for evolution definitions). Results shown in Fig. 5 suggest that cluster/isolated storms tend to be sensitive to overlake changes in nearsurface air temperature, and a greater percentage weakened as the air temperature decreased from land to lake (LAT-UWT, 0). This is the opposite of the results for linear/complex storms, which showed less of adecreaseinintensitywhentheoverlakeenvironment was cooler. 4. Discussion and conclusions Results of a climatological study examining the evolution of convective storms crossing over Lake Erie suggest that storm evolution is strongly related to storm mode, amount of time the convection has spent over the water, low-level wind speed, and near-surface temperatures over and near the lake. Most observed storms weakened after crossing the upwind shore of Lake Erie; storms rarely exhibited an increase in intensity of.3 dbz over the lake. While minor weakening of the convective storms was observed at 30 min after crossing the upwind shore, most substantial decreases in intensity were not observed until convection had spent 60 min FIG. 5. Flowcharts for (a) cluster and isolated, and (b) linear and complex storms. Numbers given in parentheses are the number of observed storm occurrences for each given parameter combination. The numbers below each parameter give the percentage of events that exhibited, after 60 min over the water, no change or moderate strengthening/major weakening. Percentages are not shown for conditions that had less than five storm occurrences.

9 OCTOBER 2012 W O R K O F F E T A L over the water. This may reflect a common time scale for atmospheric circulations associated with the convection to incorporate surface-influenced air parcels and influence precipitation rates. An important finding is that storm mode plays a critical role on the amount of weakening that storms experience while over Lake Erie. Cluster and isolated storms tended to decrease in intensity much more rapidly than linear and complex storms. This decrease in intensity was found to be most strongly correlated with the overlake near-surface air temperature (LAT) and its relation to the overland air temperature (LAT- UWT). Overlake air temperature and horizontal land lake temperature differences were diagnosed as the most influential parameters by stepwise linear multivariate regression (SLMR). It is thought that these relationships reflect the dependence of cluster and isolated convection on local near-surface stability, in agreement with their peak occurrence during times of greatest surface heating: late afternoon and evening. As a result, the majority of these storms were observed to undergo substantial weakening while over the lake, particularly when the overlake air temperature was cooler than the upwind overland air temperature (LAT-UWT, 0). Linear and complex storms were shown to be less sensitive to local lake-induced atmospheric changes and more influenced by the 3-km wind speed. This finding appears to be consistent with previous studies (Weisman et al. 1988; Rotunno et al. 1988; Weisman and Rotunno 2004; Bryan et al. 2006) that suggest dependence between linear storm intensity and low-level vertical wind shear. While the relationship between surface temperatures and intensity changes of the storms tended to be weaker than 3-km wind speed for linear and complex storm modes, surface temperatures did exhibit an interesting influence. The rate of decrease of reflectivity within linear and complex storms was inversely correlated with land lake air temperature differences. SLMR found this relationship to be nonstatistically significant. However, this relationship is supported by the flowchart analyses, which showed that organized (linear and complex) storms maintained intensity more often in events where the OLBL was cooler (relative to land). These results are in agreement with the earlier observational results of Graham et al. (2004) for Lake Michigan and the idealized modeling results of Parker (2008). It is important to note that the data used in this study do not capture potentially important features in the synoptic environment associated with the storms. Features such as elevated or midlevel instability, or the location of the convection relative to the original convective forcing, could not always be determined. Indeed, preliminary examination of storm environments suggested that complex-mode storms were generally found to form north of surface warm fronts or during the local overnight/morning hours, suggesting that a majority of them are elevated. This would likely result in minimal interactions with the OLBL. This is a topic of that should be explored in more detail in future studies. It would also be anticipated that the thermodynamic structure within the OLBL, as well as its depth, which could not be determined in this study, play significant roles in storm evolution over Lake Erie. The influences of the OLBL on convective storms are anticipated to be different depending on the size of the water body. For the Great Lakes it is possible that in events where the OLBL is highly stable, unstable air from nearby land areas may move over the OLBL and help maintain the convective storm. This seems likely for organized convective storms, where bore propagation atop a highly stable OLBL would work to lift unstable air atop the OLBL. For an oceanic coastal region, no downwind shore is available to maintain convective buoyancy. Further observational and numerical modeling studies of convective changes over lakes of various sizes and convective types are needed. To fully understand the mechanisms influencing convective storm evolution over a water body, further investigation of the relative importance of the sources and sinks of buoyant energy would be highly beneficial. Finally, smaller-scale variations in surface conditions can play important roles in storm evolution. For example, lake surface temperatures can vary rapidly with time and space over a lake. Likewise, spatial variations in atmospheric conditions, such as stability and wind characteristics, could not be investigated with currently available observational datasets. Future studies should seek to quantify the influences of such surface and atmospheric conditions on storm evolution over the Great Lakes. Acknowledgments. This research was supported by the National Science Foundation Mesoscale and Dynamic Meteorology Program (Grant ATM ) and the Cooperative Meteorology, Education, and Training (COMET) program (Grant COM UCAR S ). We appreciate the comments and suggestions from Michael Timlin and James Angel from the University of Illinois at Urbana Champaign, the anonymous reviewers, and the editor. This manuscript represents the opinions of the authors and does not necessarily reflect the views of the funding agencies or the authors affiliations.

10 1288 W E A T H E R A N D F O R E C A S T I N G VOLUME 27 APPENDIX A Observation Stations TABLE A1. The station location, station ID, and type of data collected for each of the observation stations used in this study. Stations where precip data were used are considered to be upwind of Cleveland for this analysis. Station ID in Fig. 2 Data Buffalo, NY KBUF Surface obs, raob, radar, VAD Dunkirk, NY KDKK Surface obs Erie, PA KERI Surface obs Ashtabula, OH KHZY Surface obs Cleveland, OH KCLE, 10 Surface obs, raob, radar, precip, VAD Mansfield, OH KMFD Surface obs Toledo, OH KTOL Surface obs Ypsilanti, MI KYPI, 3 Surface obs, precip Findlay, OH KFDY Surface obs Defiance, OH KDFI Surface obs Lambertville, MI KDUH Surface obs Adrian, MI KADG Surface obs Ann Arbor, MI 2 Precip Detroit, MI (metropolitan airport) KDTW, 4 Surface obs, precip Detroit, MI KDTX Surface obs, raob, radar, precip, VAD Windsor Riverside, ON 7 Precip Amherstburg, MI 5 Precip Windsor, ON 6 Precip Kingsville, ON 8 Precip Chatham, ON 9 Precip Erieau, ON CWAJ Surface obs NOAA buoy Surface obs NOAA buoy Surface obs APPENDIX B Study Parameters TABLE B1. Definitions and abbreviations of each of the parameters used for analyses in this study. Parameter Abbreviation Definition Time moved over water TMOW Time of radar scan when the geographical center of convection moved over lakeshore Starting reflectivity Start dbz Maximum base reflectivity of the convection 30 min prior to TMOW; in events where convection did not exist 30 min prior, the first radar scan with maximum reflectivity.35 dbz was used Reflectivity at the time the TMOW dbz Maximum base reflectivity of the convection at TMOW convection moved over water Reflectivity at 30 min after TMOW 130 min dbz Maximum base reflectivity of the convection at the time scan closest to 30 min after TMOW Reflectivity at 60 min after TMOW 160 min dbz Maximum base reflectivity of the convection at the time scan closest to 60 min after TMOW Delta 30-min reflectivity D30 min dbz Change in max reflectivity 30 min after TMOW: 30 min dbz 2 TMOW dbz Delta 60-min reflectivity D60 min dbz Change in max reflectivity 60 min after TMOW: 60 min dbz 2 TMOW dbz; events where the convection died or moved over the downwind shore are not included

11 OCTOBER 2012 W O R K O F F E T A L TABLE B1. (Continued) Parameter Abbreviation Definition Upwind equivalent potential temperature Downwind equivalent potential temperature UWTe (t) DWTe (t) Upwind equivalent potential temperature (K) as measured from the surface observation station nearest the storm path; data were collected from at least 24 to12 h of time convection crossed the lakeshore Downwind equivalent potential temperature (K) as measured from the surface observation station nearest the storm path; data were collected from at least 22 to14 h of time convection crossed the lakeshore Lake temperature LT (t) Temperature (8C) of Lake Erie taken from buoy 45005, when available; data were collected from 24 to14 h of TMOW Lake air temperature LAT (t) Temperature (8C) of air 4 m above Lake Erie taken from buoy 4505, when available; data were collected from 24 to14 hof TMOW Temperature difference LT 2 LAT LT-LAT The difference between the water temperature and air temperature over the lake, taken at TMOW (precold pool) Temperature Difference LAT 2 UWT Temperature difference DWTe 2 UWTe LAT-UWT DWTe-UWTe The difference between the air temperature over the lake at TMOW and the ambient upwind temperature over land (precold pool) The difference between the downwind equivalent potential temperature and the upwind equivalent potential temperature 3-km wind speed 3-km speed The wind speed (m s 21 ) determined from the VAD profile nearest the 3-km height 2 h prior to TMOW; if no 3-km VAD data were available at this time, the first 3-km wind speed value available prior to 2 h, but no longer than 4 h prior to TMOW, was used; zero surface winds were assumed 6-km wind speed 6-km speed The wind speed (m s 21 ) determined from the VAD profile nearest the 6-km height; if no 6-km VAD data were available at this time, the first 6-km wind speed value available prior to 2 h, but no longer than 4 h prior to TMOW, was used; zero surface winds were assumed Stability Various Several methods of estimating stability were investigated, such as calculating CAPE or most unstable CAPE (MUCAPE) from evening soundings, soundings closest in time to storm occurrence, and modified morning soundings REFERENCES Ashley, W. S., T. L. Mote, P. G. Dixon, S. L. Trotter, E. J. Powell, J. D. Durkee, and A. J. Grundstein, 2003: Distribution of mesoscale convective complex rainfall in the United States. Mon. Wea. Rev., 131, Augustine, J. A., W. L. Woodley, R. W. Scott, and S. A. Changnon, 1994: Using geosynchronous satellite imagery to estimate summer-season rainfall over the Great Lakes. J. Great Lakes Res., 20, Bryan, G. H., J. C. Knievel, and M. D. Parker, 2006: A multimodel assessment of RKW theory s relevance to squall-line characteristics. Mon. Wea. Rev., 134, Fowle, M. A., and P. J. Roebber, 2003: Short-range (0 48 h) numerical prediction of convective occurrence, mode, and location. Wea. Forecasting, 18, Gallus, W. A., Jr., N. A. Snook, and E. V. Johnson, 2008: Spring and summer severe weather reports over the Midwest as a function of convective mode: A preliminary study. Wea. Forecasting, 23, Graham, R., M. Bentley, J. Sparks, B. Dukesherer, and J. Evans, 2004: Lower Michigan MCS climatology: Trends, pattern types, and marine layer impacts. Preprints, 22nd Conf. on Severe Local Storms, Hyannis, MA, Amer. Meteor. Soc., 7B.6. [Available online at pdf.] Laing, A. G., and J. M. Fritsch, 1997: The global population of mesoscale convective complexes. Quart. J. Roy. Meteor. Soc., 123, Lericos, T. P., H. E. Fuelberg, M. L. Weisman, and A. I. Watson, 2007: Numerical simulations of the effects of coastlines on the evolution of strong, long-lived squall lines. Mon. Wea. Rev., 135, Parker, M. D., 2008: Response of simulated squall lines to low-level cooling. J. Atmos. Sci., 65, Rotunno, R., J. B. Klemp, and M. L. Weisman, 1988: A theory for strong, long-lived squall lines. J. Atmos. Sci., 45, Scott, R. W., and F. A. Huff, 1996: Impacts of the Great Lakes on regional climate conditions. J. Great Lakes Res., 22, Weisman, M. L., and R. Rotunno, 2004: A theory for strong, long-lived squall lines revisited. J. Atmos. Sci., 61, , J. B. Klemp, and R. Rotunno, 1988: Structure and evolution of numerically simulated squall lines. J. Atmos. Sci., 45,