Received 10 December 2015; accepted 22 May 2016; published 25 May 2016
Recently, a diagnostic method based on the stability of the upper air flow as measured using enstrophy diagnostics  was shown to successfully identify the transition between high and low amplitude hemispheric-wide flow regimes  using data from as few as two 12-h time periods. The study of  used a quantity called the Wave Amplitude Index in order to identify these transitions using time series. Then,  used the Wave Amplitude Index to show that the spring to summer transition could be more gradual or quite abrupt in the west-central plains region. It is important to extend our knowledge in identifying these transitions, in order to improve the capability of long-range forecasting, which will in turn impact different sectors of Missouri’s culture and economy.
The detection of the transition between spring and summer can be accompanied by a change in daily temperature and precipitation patterns. Spring is known to be mild and it is the season with the greatest total precipitation in the Midwest region, with the peak in total precipitation occurring in most places during May. The onset of the summer pattern, which is sometimes abrupt, is associated with a period of consistently high temperature, high humidity and dry periods interjected with fewer episodes of convective precipitation   . Also, the 500 hPa geopotential height contours can be used as a measure of the amount of potential or in the case of kinetic energy a surrogate for energy in the atmosphere (e.g.,    ). By examining a change in the 500-mb height field over an extended period of time, periods of atmospheric excitation and relative mean energy states can be used as seasonal markers  .
Two important teleconnections which have an influence on the seasonal conditions in the study region are El Niño Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO). Their phases and interactions have been shown to have unique impacts on regional temperature and precipitation regimes for the seasonal scale (e.g.,  -  ). Teleconnections, or the relationship of general circulation and climate anomalies at distances of thousands of kilometers, were first noted by Sir Gilbert Walker in the early twentieth century  . Recently,  looked at significant rainfall patterns in South America as a function of ENSO mode, noting that different ENSO phases were associated with a quasi-symmetric rainfall pattern. Also,  argued for the expansion of the idea of ENSO teleconnections into the Intra-Americas Seas, i.e. the Gulf of Mexico and the Caribbean Sea, citing that atmospheric ENSO patterns affect the two bodies of water in opposite ways. This is especially true for tropical cyclones (e.g.   ). Also, many researchers have focused on the impact ENSO has on the yearly monsoons, including Australian and Asian monsoons (e.g.   ).
Some applicable studies to this research, though, examined ENSO effects for North America. The primary mechanism for the influence of ENSO on the weather of North America is via its impact on the typical configuration of the Pacific―North American area jet stream. For example,  and  were among the first to examine the impact of ENSO over North America and the rest of the globe. Also,  examined land falling hurricanes in the US with respect to ENSO. Then,  utilized National Centers for Environmental Prediction/ National Center for Atmospheric Research (NCEP/NCAR) reanalysis data to demonstrate the connection between strong ENSO phases and upper-air data in North America, including 500-mb geopotential heights and 850-mb specific humidity. Additionally,  proposed that ENSO impacts may be bi-modal based on the primary ENSO genesis region. Specifically,  found that the ENSO impact on North America differs for an event that is an eastern Pacific versus a central Pacific ENSO event. Also, other researchers have coupled the North American ENSO effects with the impacts of PDO. Also,  examined dry and wet conditions in the U.S. Great Plains, demonstrating that ENSO and PDO coupling can intensify ENSO phases. Then  also found this impact, but they concentrated specifically on the Midwest to outline the regional effects of ENSO/PDO coupling. Another  investigated PDO and El Nino interactions and their impact on the temperature and precipitation regimes in the Amur River basin of far southeast Russia, including summer season transitions similar to  .
The purpose of this study is based on the framework of the research published by  on the spring to summer transition in the West Central Plains and their interannual variability in the temperature regime. Their study, however, examined only 20 year period from 1981-2000 and focused on ENSO variability. This study  investigated daily temperature and precipitation patterns using data from Jefferson City, MO and also used nine cooperative National Weather Service sites as a quality control. They  also utilized 500 hPa heights information in two ways; first to identify a minimum height requirement for the spring-to-summer transition date, and second as a surrogate for energy in the atmosphere using the  wave amplitude index. This research will examine interannual variability the temperatures of the region and the spring-to-summer transition over a longer time period than  , and include an analysis of decadal variability with respect to the PDO and interactions with ENSO modes. While the relationship with ENSO seems to be relatively weak, a stronger connection is found in relation to the transition of ENSO phase  .
2. Data and Methods
The study of  examined spring-to-summer transitions using observational data over the period 1981-2000. This work expands the same data set, and the time period for this research is 1948 to 2015. The 500 hPa heights are archived at the National Centers for Environmental Prediction/National Center for Atmospheric Research (http://www.esrl.noaa.gov/psd/data/reanalysis)  . These are a netCDF-based data set that can be used to display assimilated data from 1948 to the present. Also,  identified twenty summer transition dates. The mean, median and mode all fell in June; and only one date was selected in the months of May and July each (21 May, 1991 and 01 July, 2000). Then  examined temperatures and precipitation from April through September, their result allows this work to focus on the May through July period.
The statistical tests used in this study are found in standard statistical textbooks (e.g.,   ). In order to test long-term linear trends, the analysis of variance (ANOVA) technique was used, which involves the F-test. A standard t-test was also used to test correlations for significance. Mean temperatures were tested in order to analyze and compare ENSO, and PDO related variability. A two-tailed standardized test statistic (z*) was the technique used for the comparison of the sample means. Means for the total time series studied served as the expected frequencies of occurrence. All statistical tests assumed the null hypothesis, or that no a priori relationship was present among the variables tested. Confidence levels of 90%, 95%, and 99% were all tested.
The previous study  chose Jefferson City, Missouri as the principal site for data collection, and used this as the representative site of the entire East-Central Missouri Ozarks (Figure 1). Utilizing principal component analysis from  -  , and  showed that focusing on Jefferson City is sufficient to make conclusions about interannual variability in the East-Central Missouri Ozarks region as a whole. Thus following those works, this study uses Jefferson City as the representative data site, but for a slightly different region. Here we use the NWS climate divisions including the Northwest Prairie, Northeast Prairie, West Central Plains, West Ozarks, East Ozarks and the Bootheel (Figure 2), a similar region to  . Additionally, it was originally recognized that these divisions were categorized into subgroups that possessed similar climatological and geographical characteristics (Climatological Data: Missouri, 2013. http://www.ncdc.noaa.gov). This study includes the Northeast Prairie, the Eastern Ozarks, and the part of the West Central Plains. Lastly, the Jefferson City reporting site moved three times during the study period. Nonetheless, every Jefferson City site is less than one degree of longitude and latitude from another. The Jefferson City Water Plant is the station current location and the location with the most longevity, encompassing the bulk of data used here. As in  there is no systematic impact from these moves.
2.3. ENSO and PDO
El Niño Southern Oscillation (ENSO) and Pacific Decadal Oscillation (PDO) are the two main teleconnections considered for this research  , as these have the strongest correlation with regional weather  . Other teleconnections such as the North Atlantic Oscillation and Arctic Oscillation did not correlate as strongly. The Japan Meteorological Agency (JMA) ENSO index is available through the Center for Ocean and Atmospheric Prediction Studies (COAPS) from 1868 to present (JMA sea surface temperature (SST) ENSO Index
Figure 1. From Ratley et al. (2002), the location of the study region including Jefferson City, Missouri. This area is bounded by the edge of the East Ozarks to the west, the Mississippi River to the east, the Missouri River to the north and the Missouri/Arkansas border to the south.
Figure 2. United States Department of Commerce Weather Bureau’s Climatological Data: Missouri Section, Divisions, 1957 to Present, image courtesy of the Climate Prediction Center (CPC).
http://www.coaps.fsu.edu/jma). The JMA SST ENSO Index has been widely used in other published works (e.g.,  ,  , and  are some recent examples), and is also used here (Table 1). Also,  found that, however, while the JMA index is more sensitive to La Niña events than other definitions, it is less sensitive than other indices to El Niño events. The JMA classifies ENSO phases using SST within the bounded region of 4˚S to 4˚N, 150˚W to 90˚W. The JMA defines the inception of an ENSO year as 1 October, and its conclusion on 30 September of the next year.
The Pacific Decadal Oscillation (PDO) positive and negative modes are catalogued also by the Center for Ocean-Atmospheric Prediction Studies (COAPS). The most important effect of PDO is how it interacts with ENSO during certain phases to create an enhanced effect on temperatures and precipitation variability (e.g.,   ). The characteristics of these modes are less pronounced than those for ENSO due to the fifty- to seventy-year cycle of PDO (e.g.,   ). COAPS describes the PDO phases (Table 2) as follows: the High (+) PDO phase is characterized by cold SSTs in the north central and western Pacific Ocean, warm SSTs off the western coast of North America, and a deep Aleutian low, while the Low (−) phase of the PDO is characterized by warm SSTs in the north central and western Pacific Ocean, cool SSTs off the western coast of North America, and no pronounced Aleutian low. The start date of the current PDO cycle from COAPS agrees with the independent analysis of  , based on the same SST patterns.
Table 1. Center for Ocean-Atmospheric Prediction Studies Japan Meteorological Agency El Niño Southern Oscillation Index, 1948 to present. Modes are El Niño (EL), La Niña (LA) and Neutral (NEU).
Table 2. Center for Ocean-Atmospheric Prediction Studies Pacific Decadal Oscillation Index, 1910 to present. Modes are high (positive) and low (negative).
2.4. Ratley et al. 2002 Temperature Criterion
We used the criterion of  , which is set in order to filter out synoptic-scale perturbations. The thresholds for surface temperatures were set at fifteen days. The criterion used by  for identifying the spring to summer transition date using 500 hPa heights was the start date for the first period of at least 10 consecutive days where the heights are in excess of 5820 m over the entire study region. The 5820 m contour is frequently used as an indicator of the difference between the baroclinicity of the mid-latitudes and the quasi-barotropic nature of the subtropics  . This height minimum is typically the trailing equatorward edge of tightly-packed mid-latitude height gradients, as also demonstrated in a derecho study by  , and in study of summer perturbations by  . Additionally, since the 5820 height contour is at the equatorward edge of the 500 hPa height gradients, they frequently coincide with gradients of potential temperatures near 350 K on the 2.0 Potential Vorticity Unit (PVU) (e.g,  ) in identifying the general location of the jet stream. Potential vorticity is conserved on surfaces of potential temperature, but the above argument is not strictly true as potential temperature or potential vorticity will not be strictly conserved on a height surface (e.g.,  ). Nonetheless, many papers have demonstrated the general utility of potential vorticity quantities in defining features in the jetstream. Thus we are confident that this choice of contour has dynamic relevance to the study here.
The temperature criteria for a summer start date are as follows from  , and using a 15-day window for surface based data is sufficient for removing smaller-scale impacts  :
1) The first date of a fifteen-day consecutive period where the mean temperature exceeds 70.0˚F (21.1˚C), and at least ten of those fifteen days have a temperature at, or exceeding, 75.0˚F (23.9˚C).
2) The first date of a fifteen-day consecutive period where the maximum temperature exceeds 77.0˚F (25.0˚C) and at least ten of those fifteen days have a temperature at, or exceeding, 82˚F (27.8˚C).
Temperatures are reported in the National Climatic Data Center in degrees Fahrenheit, and then converted to Celsius.
2.5. Blocking Criterion and Integrated Regional Enstrophy (IRE)
3. Climatological Analysis
3.1. Spring to Summer Transitions in the Study Region
The results of this study produced a transition date using each criterion, where data was available, and a final spring-to-summer transition date based on these three criterions for the 68 year period (1948-2015). There were a few years before 1980 in which Jefferson City temperature data were absent, including 1967. The 500 hPa height criterion provided a date for every year except 1967. In order to provide a possible date for that year, Columbia Regional airport data were used to provide continuity for some of the analyses, and the date identified was June 3 and 4 for each temperature criterion. It was also noted that the 500 hPa height criterion beginning with June 4 was only met for seven consecutive days, and at nine of the ten day period (only one time did not). During the summer as a whole, the longest consecutive period for meeting the criterion was eight days. Thus, a transition date of 4 June for 1967 seems quite plausible based on data used for this analysis.
An examination of Table 3 shows that the mean spring and summer temperatures and their standard deviations for the entire period, and these are close to the current 30 year means and standard deviations. The trend for summer temperatures would display a cooling trend using linear regression, but this trend was not statistically significant, a result identical to those of  . The spring season, however, Figure 3(a) showed a significant warming (at the 90% confidence level), and this resulted in a statistically significant decrease in the spring to summer temperature difference (at the 90%). This result is consistent with observations from other parts of the earth, and is consistent with the prevailing idea that, in a warmer world, the cold season will warm at a faster rate  .
Figure 3. (a) Mean spring season temperature (˚C) versus time (years) and (b) the final date determined (days from 1 May) for each year. The blue dotted line is the linear regression line.
Table 3. Statistics for the three criterion used to identify the spring-summer transition date, and temperature data for the spring and summer seasons and their differences from 1948-2014.
The average transition date for the entire data set, 1948-2015, is 12 June, not including 1967. The median transition date is also 12 June, and the median values for all variables in Table 3 were the same when the means were rounded. June was, by far, the most common month for final transition dates, accounting for 85.3% of all chosen values, and 78.0% of all transitions occurred within one standard deviation of the mean. Eight transitions occurred in May. These characteristics imply that the data set is close to a normal distribution, but one that is more peaked (kurtosis). Since the 1967 transition date was identified as 3 June, including it in the analysis does not change the overall statistics.
Concurrent with the decreased temperature differences between spring and summer, the trend for the final summer onset day is indicated to be occurring earlier (Figure 3), a trend that was significant at the 99% level for the 500 hPa data, and the mean temperature criterion data. The trend for the decreased temperature difference was significant at the at the 90% confidence level. Only the maximum criterion trended toward a later date, and this was significant at the 95% confidence level. These contradictory trends are possible since the summer minima trended higher, but also the dew point trend was significantly higher  . A higher dew point correlates with higher minima, and this is also consistent with the prevailing character of a warming climate  .
A further examination of the individual criterion demonstrates that the daily mean temperature and the final date chosen were similar statistically (Table 3) and the correlation between the two time series was 0.73, a value significant at the 99%. This was by far the strongest association in this data when cross correlating each time series. The 500 hPa and maximum temperature criterion were similar in their variation (standard deviation was 14 days), but the maximum temperature produced the earliest onset date and the 500 hPa the latest. While the 500 hPa criterion correlated with the final date at the 95% confidence level (0.28), the maximum criterion did not correlate with any other criterion at a statistically significant level. The 500 hPa height, the mean temperature criterion, and the final date chosen correlated positively with the summer mean temperature as might be expected, and these correlations were significant at the 99%, 90%, and 99% confidence level, respectively. There was only one other statistically significant relationship between these criterion and temperatures, and that was between the mean temperature criterion and the spring mean temperatures.
At the beginning of this work, it was proposed that the 500-mb height data would be the most useful criterion, as upper-air analyses are not subjected to surface observation errors and influences. However, there are errors inherent in the 500-mb data as well, and the years before the satellite era may possess larger error (e.g.,  ). Based on the statistical analysis above, an argument could be made that the mean temperature criterion would be the most useful for identifying the spring-summer transition date in the study region.
A comparison of this work and  (Table 4) shows a similar pattern to the statistics for the larger data set in that the 500 hPa criterion produced the latest date while the maximum temperature date was the earliest. Also, the mean temperature criterion and final date were within one day of each other in both studies. The study of  found that for 12 of the years (60%) of their 1981-2000 dataset, all three summer onset dates fell within a 10-day moving window. In the longer data set here, only 20 years (33%) were characterized in the same manner. In order to determine why this difference was so large we sampled the same 20 years as in  using our criterion. Here we found 11 years (55%) fit the same criterion, or a difference of only one year.
An examination of the two studies shows that a slightly different methodology was used to determine the onset date. The rules employed by  specified that the final date chosen must match one of the other criteria, which is more arbitrary. Here, we used similar rules except that the final date was the average date provided by the available criterion for that year, and thus, did not have to match any of the three criteria. The average of the two analyses gave final dates that were an average of 4 days apart. In this study of the period of  , five transition days (25%) fell within twenty-four hours of all criterion. Thirteen transition days (65%) fell within five days of one another, still within the time frame of the synoptic period. Seventeen days (85%) were within ten
Table 4. As in Table 3 except for relevant values from Ratley et al. (2002).
days, or planetary-scale time frame. The outcomes here chose an earlier transition date 45% of the time, and a later onset date 50% of the time, thus overall the change in assigning the final date did not result in large difference between the two studies.
3.2. Pacific Region Blocking
There has been speculation that frequency of Northern Hemisphere blocking has increased during the last two decades. Comparing our results to those of  (see Table 5 and Figure 4), indicate there has been a very sharp increases in the frequency and number of blocking events and blocking days in the Pacific region (140˚E - 100˚W) spring (significant at the 99% confidence level) and summer (significant at the 95% confidence level) seasons when the study is extended to the current time. There was, however, no significant change in the intensity of blocking in this region. The study of  covered the period 1970-2000, thus the sharp increase in blocking has occurred since this time. Also,  associated more and stronger Pacific Region summer blocking to cooler, and wetter summers in the region. Here we find a strong correlation between the occurrence of spring and summer blocking and regional summer temperatures (more blocking, cooler conditions for the study region), but not significant at the 90% and 95% confidence level as in  . This is especially true if only blocking east of dateline is included.
This and many other studies, even those going back to the 1980s (e.g.,   ) show that the occurrence of Pacific Region blocking, especially in the East Pacific is associated with cooler winter temperatures in the region. This same correlation appears to be true for the summer season as well  . The occurrence of blocking, however, did not correlate with any of the onset criterion here for the final day chosen. This may be expected since the occurrence of blocking is quasi-random. However, the association of blocking with the character of seasonal temperature and precipitation characteristics is well-documented. Thus, it would seem that the influence of blocking may be more indirect and this issue will be explored versus interannual variations.
3.3. Synoptic Climatology and Dynamic Analysis
In order to examine the effectiveness of this criterion in identifying a date for the spring-summer transition, the 850 hPa and 500 hPa heights were used to examine the lower (Figure 5) and mid-tropospheric (Figure 6) air flow for the 30 day period before and after the nominal transition date in order to identify whether a difference existed. Five particular years were chosen in order to represent an early onset (2 June 2012) a late onset (28 June
Figure 4. The frequency of (a) spring (FMAM) and (b) summer (JJA) blocking in the Pacific Region. The blue dotted line is a linear regression model.
Table 5. The mean spring and summer Pacific Region blocking (140˚E to 100˚W) from 1969-2015 as compared to Wiedenmann et al. (2002). The comparison with the earlier work is shown as a departure (current―Weidenmann et al., 2002).
(a) (b)(c) (d)(e) (f)(g) (h)(i) (j)
Figure 5. The mean 1200 UTC 500 hPa height map for the East Pacific and North America for the (a), (c), (e), (g), (i) 30 day period prior to, and (b), (d), (f), (h), (j) 30 day period following the identified onset date found in this study. The years used are (a), (b) 1967, (c), (d) 1971, (e), (f) 2004, (g), (h) 2008, and (i), (j) 2012. The identified dates for transition were, 4 June, 17 June, 28 June, 11 June, and 2 June, respectively. The contour interval is every 6 dam.
(a) (b)(c) (d)(e) (f)(g) (h)(i) (j)
Figure 6. As in Figure 5, except at the 850 hPa level and a contour interval of 30 dam.
2004), one close to the mean (11 June 2008), and two of particular interest (4 June, 1967―the date inferred above in section 3a for the year that lacked a criterion, and 17 June, 1971―Midwestern drought). Table 6 shows the summer mean temperature and total precipitation and their relative anomalies for each year relative to the 1981-2010 anomalies.
In all five cases at 500 hPa (Figure 5), the 5820 m contour was located within or equator ward of the study region in the 30 day period before the transition date. There was either split flow or ridging over the Rockies as is typical of the colder season. In the 2012 case, the flow was more zonal and the jet stream was clearly north of the region (but not the 582 m contour). In the 30 day period following the onset date, the 5820 m line and the jet stream was clearly poleward of the region, and summer season ridging was established in North America. Note that during the prior period for the year 2004 there was blocking in the Alaska Region. Similar to the case for the hot dry summer of 2012, it was noted that the years 1980 and 1954 also had earlier onset dates and a similar difference between the pre- and post-transition date (not shown). Thus, this criterion is rather effective for intensifying this change in the character of the 500 hPa flow over the region.
At the 850 hPa level (Figure 6), the post-onset dates for 1967, 1971, and 2012 show that the 850 hPa flow over the study region was more meridional and the stronger flow was to the west over Texas and Mexico, while the 2004 and 2008 flow regimes were more zonal and originally out of the gulf of Mexico. Table 6 shows that the former three years were dry over the study region, two of them extremely dry. This result for dry years is consistent with the study of  who looked at summer season dew points in the region. The latter two years were wetter than normal. During the pre-onset date period, the flow was more zonal over the study region for all the years except for the year 2012. Like the 500 hPa analysis, however, these were all different from the post- onset flow region.
Here we demonstrate using the integrated regional enstrophy (IRE) that the stability of the northern hemisphere flow (IRE) increases (decreases) in association with the time period identified for the spring-to-summer transition date. Figure 7 shows a time series of the IRE quantity from 1 May to 31 July for both 1971 and 2008. The transition date in 1971 was identified as 17 June. The IRE diagnostic for that year (Figure 7(a)) was higher for the period before 17 June, and considerably less following this date. However, the transition during May and June was more gradual. In the case of 2008, the transition date chosen was 11 June, and the IRE diagnostic increases rapidly (possibly indicating instability and flow regime transformation―  ), and then decreases to a
Table 6. The mean temperatures (˚C) and total precipitation (mm) and their anomalies for the years 1967, 1971, 2004, 2008, and 2012. The * represents years that were greater than one standard deviation from the normal (1.2˚C and 13.1 mm).
Figure 7. The Integrated Regional Enstrophy Diagnostic (IRE) for the time period 1 May to 31 July (in days from 1 May) for (a) 1971, and (b) 2008.
much lower level for the rest of the series. These examples demonstrate that the IRE diagnostic does indicate that the flow regime changes occur in the time frame identified by the criterion used here. Additionally, the IRE indicates that the summer season flow is more stable, and thus weather prediction is more favorable as shown statistically for the summer season in region  .
4. Interannual and Interdecadal Variability
In order to identify the interannual and interdecadal variability a periodogram of the summer temperatures, the final date chosen, and Pacific Region blocking were constructed using the method of cycles (e.g.,  ) (Figure 8). All of these show statistically significant interannual variability, but there is little indication of interdecadal variability, as was found for regional precipitation events (e.g.,   ), except for the blocking time series which shows significant variability beyond ten years. The interannual variability in Figure 8(a) for example can be inferred from the diagram which shows a significant peak at wave number 18, which in the 68 year data set represents variability with a period of three to four years. Interdecadal variability in the blocking occurrences (Figure 8(c)) is evidenced with a peak at wave number three and ten (46 years total) or at 15 and 5 years, respectively. Table 7 shows the number of early versus late onsets for each mode of ENSO and the PDO confirming that there is no strong interdecadal signal in the spring-to-summer onset. As may be expected for neutral years, there is no discernable tendency toward late or early onsets, and the same outcome can be inferred for years after El Niño onset, which  also found. However, for years after La Niña onsets, there is a stronger tendency for early onsets, especially during −PDO. Both La Niña years in Table 6 were early onsets. Additionally,  studied a period that contained few La Niña onsets, thus there is no basis for comparison with this study.
Table 8 shows the occurrence and of blocking within the Pacific region for each year shown in Figure 5, Figure 6. In section 3 we found a correlation between the more frequent occurrence of blocking and cooler,
Table 7. The number of early or late onset dates for summer stratified by ENSO phase for +PDO/−PDO.
Table 8. As in Table 5, except for the years in Figure 5, Figure 6 where blocking was examined.
Figure 8. The power spectra versus wave number (per 67 years) for the (a) summer season temperatures, and (b) the final date time series from Figure 3(b). The blue (green) dotted line is the 95% confidence level for a white (red) noise spectrum. (see Wilks, 2006). For (c) as in (a) and (b) except per 46 years for the occurrence of Pacific region summer blocking (from Figure 4(b)).
wetter summers. Both 2004 and 2008 were cooler, wet summers, while 2012 was the least blocked season after 2000. This season was very hot and dry. Additionally 1980 was tied for the least blocked summer season for the entire period 1969-2014, and 1971 was about average for the summer temperatures, but quite dry. The regular occurrence of blocking would allow for more variable summer weather and more frequent intrusions of cooler air in the region. Separating the occurrence of blocking by ENSO phase would demonstrate that more blocking occurred in years in which the transition was toward the El Nino (2.4 events) than those when the transition was the other way (1.75 events).
The study of  found no statistically significant differences between the summer temperatures as a function of ENSO in their 20 year study. Table 9 shows the mean seasonal temperature for each case in Table 7 using the larger data set. Similar to  or  , −PDO springs and summers show slight ENSO variation for spring and summer temperatures following El Niño onset averaging slightly warmer than summers following La Niña onset. This result matches the studies of  and similar studies from this group that there is weak ENSO variability in −PDO years. In +PDO yeas, which overlapped with most of the years in the  study, they found warmer summers following El Niño onset was statistically insignificant. Here we found that these summers were warmer at the 90% confidence level, but the spring-summer temperature difference is significantly larger at the 95% confidence level. As in Ratley et al.  , we can conclude that useful information may be found by relating the change in ENSO phase to summer temperatures, and  (  ) found a similar phenomenon for a study of the summers of eastern (western) Russia, respectively.
In examining the changes in ENSO phase, defined as the difference between the new phase of 1 October and the phase of the current summer,  examined the tendencies for certain ENSO phase changes to be associated with early or late onsets. Here, we will define changes toward the La Niña (El Niño) direction as positive (negative), and no change in phase as a null case. Then,  found that positive (negative and null) transitions were associated with late (early) onsets. The results of our study indicate no such tendency in the much longer data set (Table 10). In fact it appears that all possible transitions are equally likely to produce early and late transitions overall.
Table 9. As in Table 7 except for mean spring and summer temperatures.
Table 10. ENSO phase changes for early and late summer onset dates where a positive phase change is toward La Niña and negative is toward El Niño and null is no change in phase from one year to the next.
Finally,  found that summers associated with a transition toward La Niña or steady state La Niña conditions produced warmer, drier summers partly based on conjecture. Here, we tested all 67 seasonal transitions, and positive transitions were indeed warmest (24.6˚C), followed by null transitions (24.2˚C), and negative transitions (24.0˚C). This result is not statistically significant, but does support the assertions of  and  who also found La Niña transitions to be warmer. Narrowing the subset to include only El Niño (−) and La Niña (+) transitions and their null cases (back-to-back years) for both demonstrates that La Niña transitions and steady state La Niña years produced a mean temperature of 24.8˚C, versus 24.0˚C for El Niño onsets and steady El Niño conditions. Including only La Niña or El Niño onsets and eliminating their null cases, the average summer temperature was 25.1˚C versus 24.1˚C. Thus, the conclusion of  that transitions toward La Niña were warmer is valid, or even stronger, in the much larger dataset.
The methodology of  for identifying the spring-to-summer transition was used for a longer dataset, this one covering a 68 year period from 1948-2015. They  used surface temperatures to derive a criterion based on the daily maximum and mean temperatures, and a 500 hPa height map criterion with the goal of identifying the spring-to-summer transition date for a 20 year period. They used each criterion in order to determine a date that matched one of the three dates found by the available individual criterion. This study modified the practice and identified the mean date of the three available criteria as the spring-to-summer transition date. It was found that:
・ the statistics produced by this study were similar to those produced by the  study for the spring-to-summer transition date, in that the 500 hPa height mean transition date was the latest, and the maximum temperature mean transition date the earliest,
・ the daily mean temperature criterion was the closest to the overall date (June 12) found by compiling the three criterion, it was only one day different in both this and the earlier study. This criterion and the final date were very strongly correlated, significant at the 99% confidence level.
・ the mean spring temperatures trended upward at a significant rate, while the summer season temperatures were slightly cooler providing for a smaller annual cycle, and a significant trend (at the 99% confidence level) toward an earlier onset date,
・ when comparing the results of this study to the smaller time period of the  work, it was found that our modification for identifying the final transition date had little impact on the overall statistics,
・ Pacific region blocking showed a strong increase in occurrence and days, but not intensity, and more summer season blocking was associated with cooler summers and the transition toward El Niño,
・ the 850 hPa and 500 hPa height fields show that the 30 day period before five chosen onset days was different from the 30 day period following onset. At the 500 hPa level, the 582 contour was located over or equatorward of the region during the spring and poleward during the summer. At 850 hPa, the summer season flow was more meridional and during the drier years the strongest flow was to the west of the study region,
・ the  study found that summer seasons preceding La Niña onset were warmer and the transition date was later that the mean. Here we found that the association of La Niña onset summers being hotter summers was confirmed and the result was stronger than that of the  finding. However, this study did not find a tendency for early or late summer onsets in association with ENSO,
・ the IRE diagnostic used for example in  indicated that the dates chosen for the spring-to-summer transition were close to decreases in the northern hemisphere enstrophy demonstrating that the criterion proposed here has utility in identifying the seasonal transition,
・ and finally, this study showed only slight interdecadal variability in ENSO behavior, unlike other studies of different variables and seasons for this region, in particular it was noted that ENSO summers in the ?PDO phase were similar overall with little variation, while in the +PDO the summers following La Niña were cooler and associated with a smaller spring-to-summer temperature change.
In conclusion, these results would have applicability for seasonal or long-range forecasting in the region of study. The studies of  and  both demonstrate that SST anomalies and the occurrence of summer Pacific Region blocking correlate strongly to the summer weather and climate of the study region. These variables can then be used to predict the conditions of the summer season  . The information gained here provides another variable that can be used for long-range seasonal prediction, and this information should be on interest to the regional agriculture community.
The authors would like to thank Professor Neil I. Fox for his earlier comments on this work, as well as the anonymous reviewers who helped to make this a stronger contribution.