Electricity consumption is contributed by many types of human activities, such as heating, air conditioning, lighting in both business and residential sectors, and major contributions come from operating equipment in industrial sectors. Whilst lighting and operating equipment might not be directly linked to climate change, heating and air conditioning have a direct impact on air temperature  . All the climate-change-related impacts on electricity demand and supply can be easily observed from the quantifiable effects of temperature on the use of heating and air conditioning, and these numbers are usually described by different measurements based on the concept of heating degree days (HDDs) and cooling degree days (CDDs).
HDDs is defined as the sum of negative deviations from the actually measured temperatures to the reference temperature (or base temperature) over a given time period; in contrast, CDDs indicates the sum of positive deviations from the average temperatures to the reference temperature over a given time period. The data frequency of the given time period is usually daily, weekly or monthly. The reference temperature is defined by the temperature level without additionally using electricity for heating or cooling. That is, if the air temperature is comfortable for humans, there will be less electricity consumption for heating or cooling.
The reference temperature can be generally considered to be 18.3˚C (65˚F)  . However, Parkpoom and Harrison  used 11.7˚C (53˚F) to be the reference temperature in Thailand; Howden and Crimp  determined 17.5˚C (63.5˚F) to be the reference temperature for Sydney; Ahmed et al.  proposed 14.3˚C (57.7˚F) as the reference temperature for the State of New South Wales in Australia after their calculation; Zachariadis and Hadjinicolaou  employed 18˚C (64.4˚F) and 22˚C (71.6˚F) respectively to be the reference temperature of HDDs and CDDs for the area of Mediterranean Europe. In sum, there could be different reference temperatures within different geographical regions.
Global warming could lead to increases in CDDs and decreases in HDDs, concluded by Benestad , whose report indicates that climate change could trigger more energy consumption due to air conditioning in the hot areas. De Cian et al.  used the panel data from 31 countries to investigate the relationship between energy consumption and variations in temperature. Their empirical results suggest that higher average temperature leads to more energy consumption during hot seasons in the warmer countries, but less energy is consumed during cold seasons in the colder countries.
Hekkenberg et al.  assessed the electricity demand pattern in the relatively temperate climate of the Netherlands. They used daily data over the period from 1970 to 2007 to investigate possible trends in the temperature dependence of electricity demand. Although the Netherlands has the minimum electricity demand in the summer months, however, their empirical results showed significant increases in the temperature dependence of electricity demand in the months of May, June, September, October and during the summer holidays. That is, their alarming result sends a signal to raise future expectations for additional peaks of electricity consumption in summer under the in the influence of climate change.
Moral-Carcedo and Vicéns-Otero  figured out that the relationship between electricity demand and temperature is nonlinear, and the nonlinearity is reflected on the threshold temperatures. They employed the threshold regression model (TR) and the logistic smooth transition regression (LSTR) model to build the relationship between electricity demand and temperature in Spain using daily data from 1995 to 2003.
In their research, they created the variable of working day effect to capture the variations of electricity demand caused by the activities in the industrial and commercial sectors as well as by the behaviors of households during holidays and on working days. Hence, they could eliminate those effects from electricity demand, then focus more on the pure effects of temperature on electricity demand. Their results showed that the threshold temperatures of the TR model are 15.5˚C (59.9˚F) and 18.4˚C (65.1˚F), and the threshold temperature of the STR model is 18˚C (64.4˚F).
Bessec and Fouquau  investigated the relationship between electricity demand and temperature in 15 European countries over the period from 1985 to 2000 using monthly data. They applied a panel smooth transition regression (PSTR) model to describe the relationship between electricity demand and temperature in those countries and find threshold temperatures for those countries. In addition, in order to estimate the pure effects of temperature on electricity demand, they also followed Moral-Carcedo and Vicéns-Otero  , and used dummy variables to represent summer holidays and time trends to filter out other source of electricity consumptions. Their results showed that the nonlinear pattern was more pronounced in the warm countries among the 15 European countries.
Lee and Chiu  used the PSTR model and took into account the potential endogeneity biases to examine the relationship between electricity demand and temperature of 24 OECD countries over the period from 1978 to 2004. They provided evidence of a U-shaped relationship between electricity consumption and temperature of 24 OECD countries, and the threshold temperature is approximately 11.7˚C (53˚F).
In sum, to summarize the literature mentioned above, we can highlight two main findings. First, the relationship between electricity consumption and temperature shows nonlinearity in the past cases, so when establishing an econometric model for cases in Taiwan to estimate the effects of temperature on electricity consumption, we should consider possible nonlinear relationship between electricity consumption and temperature. Secondly, the threshold temperature has some policy implications, such as guidance for the management of electricity demand and supply, strategies for mitigating the impact of climate change on electricity.
To give an example of policy implications on electricity management, the Taiwanese government has introduced a policy since the year 2011 to save energy by asking public sectors to operate air conditioners only if the air temperature is higher than 26˚C (78.8˚F). In addition, once the real threshold temperature is found, it can be applied to computation of the data of CDDs in Taiwan to describe the patterns between temperature and electricity consumption both in the past and in the future. That is, if global warming leads to more temperature degree days, we should consider more power system expansion planning in Taiwan to meet the possible increases in electricity demand in the future. Therefore, we believe that it is worth further discussing how to model the real relationship between electricity consumption and temperature.
The objective of this study is to utilize the nonlinear econometric approach (STR model) to analyze the effects temperature has on electricity consumption in Taiwan. The estimated results of the STR model provide two mainly contributions to this study. First, we show the evidence of positively nonlinear relationship between electricity consumption and temperature in Taiwan. Secondly, we find that the average threshold temperature for Taiwan is about 26.384˚C (79.3˚F) over the period from 1983 to 2012. Furthermore, there are variations of threshold temperatures among different sample periods, and range of threshold temperatures lies between 25.364˚C and 27.156˚C. The contributions of this study could be turned into policy implications for policy makers.
The remainder part of this study is organized as follows: Section 2 describes the data source, data descriptive and data processing. Section 3 presents the econometric methodology and the empirical model. Section 4 provides our empirical results. Section 5 is the conclusion of this study.
2.1. Data Source and Descriptive
In this study, we use monthly time-series data which cover the period from 1983 to 2012. The original data of electricity consumption per capita (kWh) are collected from MOEABOE  , and the gridded dataset of historical climate information from TCCIP  is used to compute the monthly average temperature (˚C) over the period from 1983 to 2012.
Table 1 displays the descriptive statistics on monthly average temperature in the different time period over 1983 to 2012. In Table 1, the mean temperature of the past three decades is between 21.997˚C and 22.383˚C. The coldest month in a year are usually January and February, and the hottest month in a year, July and August. In addition, the mean temperature for summer (June, July and August) is between 27.127˚C and 27.460˚C, and the stand deviation of temperature in
Table 1. Descriptive statistics on temperature.
Notes: () refers to the month of the data value.
summer is between 0.556˚C and 0.840˚C, which indicates that we observe minor variations of temperature in summer in the past three decades in Taiwan.
Figure 1 shows the average electricity consumption and the average temperature, respectively. We can see that the month of the largest electricity consumption per capita in a year is August; however, the month of the highest temperature in a year is July, meaning the non-temperature impacts on electricity consumption should be considered. Therefore, in the Section 2.2, we filter out the effects of non-temperature factors on electricity consumption.
2.2. Filtered Electricity Consumption
In order to examine the pure effects of temperature on electricity consumption, we firstly remove the effects of other factors on electricity consumption   . Especially, Bessec and Fouquau  indicated that three major components must be considered when we filter out the other effects affecting electricity consumption. The first component is the demographic trend, the second component is the technological trend, and the third component is the monthly seasonality related to the activity. However, our data of electricity consumption are divided by population, so we can say that our data of electricity consumption have removed the effects of demographic trends. Then, we follow Moral-Carcedo and Vicéns-Otero  as well as Bessec and Fouquau  , the two last components will be filtered out from electricity consumption by employing Equation (1),
where represents the electricity consumption at time t; t denotes the time trend; D is a dummy variable. We define the value of July and August equal to one, and the values of other months are equal to zero. The dummy variable is used to remove the effects of summer holidays on electricity consumption  . stands for the filtered electricity consumption, and it is also the estimated
Figure 1. The average electricity consumption and average temperature (1983-2012).
residual from Equation (1), and Table 2 shows the estimated result of Equation (1).
Figure 2 is the scatter plot of filtered electricity consumption and temperature over the period from 1983 to 2012. We also repot a regression line between filtered electricity consumption and temperature with a polynomial of order three. However, we cannot see the U-shape relationship between electricity consumption and temperature in Figure 2. In Taiwan, people usually use more gas and oil products for cooking and heating in the cold seasons (from Oct. to Feb.), and air conditioner for cooling in summer, which is the reason for the non U-shape relationship between electricity consumption and temperature in Taiwan. However, although Figure 2 looks like there is a positive linear relationship between electricity consumption and temperature, we still believe that there could be a threshold temperature in the relationship between electricity consumption and
Table 2. The estimated result of filtered electricity consumption.
Notes: 1) () stands for standard deviation. 2) *** refers to the significance level at 1%.
Figure 2. Filtered electricity consumption and temperature (1983-2012).
temperature. At this stage, we will employ smooth transition regression models  to investigate whether the nonlinear relationship exists between electricity consumption and temperature in Taiwan or not.
3. Methodology and Empirical Model
This study employs the STR model to analyze the pure effects of temperature on electricity consumption in Taiwan. However, before the estimation of the STR model, we firstly have to test whether the time series data are stationary or not; that is, we have to conduct a unit root test for each set of time series data. Hence, we will introduce the methodology of unit root test which we have used for this study in the following section.
3.1. Unit Root Test
The stationary of time series data is usually examined by the Augmented and Dickey Fuller (ADF) test and Philips-Perron (P.P.) test in the past literatures, and these two tests were provided by Dickey and Fuller  and Philips and Perron  , respectively.
However, both the ADF test and P.P. test have not considered the possibility of a structural break in the time series data. Therefore, to solve this problem, we employ the unit root test with structural breaks provided by Saikkonen and Lütkepohl  and Lanne et al.  . If there is a shift in the data generating process (DGP) of the level data, it should be taken into account in the unit root testing. The shift function and deterministic trend will be included in the DGP of the time series , such as Equation (2),
where δ and φ are unknown parameters; is the error term generated by AR(p) process with possible unit root. Here, two types of shift function are applied for the unit root test in this study. The first type is a simple shift dummy variable with shift date ; if , the dummy variable is equal to 1; on the contrary, if , the dummy variable is equal to 0. The difference in the shift function makes an impulse dummy, which can be seen as a rational function in the lag operator applied by . Moreover, the actual shift term of impulse dummy is , where δ is a scalar parameter between 0 and 1; and is a two dimensional parameter vector.
Once a break point is fixed, Saikkonen and Lütkepohl  and Lanne et al.  suggested the unit root test (Equation (2)) could be estimated by the generalized least square (GLS) procedure under the null hypothesis of unit root. In addition, Lanne et al.  also provided the critical value for this unit root test.
3.2. Smooth Transition Regression (STR) Model
The STR model is widely used to describe the nonlinear relation of time series data. The univariate form of STR was proposed by Chan and Tong  , subsequently developed by Luukkonen et al.  and Teräsvirta   . Areosa et al.  further showed the estimation of STR models with endogenous variables. Hence, according to Teräsvirta  , the STR model can be specified as Equation (3).
where t represents time dimension; is a dependent variable; represents the transition function with the transition variable ; π and θ represent the linear part of the model and nonlinear part of the model, respectively. γ is a slope parameter which shows the speed of transition from one regime to another regime, and c is also the extreme threshold of the transition variable. In the equation, , indicates the optimal autoregressive term for p lag lengths of dependent variable, and means the independent variable of q lag lengths. π and θ represent the linear part of the model and nonlinear part of the model, respectively. If , the STR model will reduce to the threshold regression (TR) model, meaning that if the transition variable is larger than c, the transition function will be one. However, if the transition variable is smaller than c,the transition function will be zero. In addition, if ,the STR model will change to a linear model.
Generally, the transition function usually can be distinguished into two types of function forms, namely, the logistic function and the exponential function. They are expressed as Equation (4) and Equation (5), respectively.
The first step of the STR estimation starts with examining whether there is a nonlinear relationship between the dependent variable and the transition variable or not. If the nonlinear relationship exists between the variables, then the second step of STR estimation will investigate the number of regime switch. Both of the two steps can be identified via the coefficient test on the following auxiliary regression, as Equation (6) shows:
The null hypothesis of the linearity of the relationship between dependent variables and independent variables versus the alternative hypothesis of the nonlinearity of the relationship between dependent variables and independent variables can be examined by the null hypothesis of coefficient test as follows:
The test statistic of H1 holds F distribution, and F test is suggested for coefficient test because of its better small properties    . If the testing result rejects the null hypothesis of H1, it indicates that the nonlinear model should be selected to describe the relationship between variables. Subsequently, we can select an appropriate nonlinear model via doing three types of coefficient tests, all of which contain F distribution. The three coefficient tests extended from the null hypothesis are as the following:
If the testing result shows that the rejection of H3 is the strongest, the appropriate nonlinear model is either the LSTR2 model or the exponential smooth transition regression model (ESTR), either of which should be chosen by doing the null hypothesis test of . If H0 is not rejected, it means that there should be one regime switch, and the optimal model will be the ESTR model; otherwise, the appropriate model will be the LSTR2 model. However, if the testing result indicates that the rejection of H2 or H4 is the strongest, the appropriate nonlinear model will be the logistic smooth transition regression model with one regime switch (LSTR).
Once the linear model is rejected, the next step is to select the appropriate transition variable and estimate the parameters of transition function. Teräsvirta  suggested the transition function should be standardized to make γ scale-free, which means dividing the exponent in by the standard deviation of , as the expression of Equation (11) shows.
A grid search method is employed here to find the parameter γ and c, and the model with the minimum value of the sum of square residuals (SSR) from the grid search method will be used to provide an initial value of γ and c for the initial estimation of the STR model.
3.3. Empirical Model
Therefore, we establish a smooth transition regression model to describe the nonlinear relationship between electricity consumption and temperature. Considering the assumption of exogeneity for electricity consumption to temperature, Chen et al.  examined the relationship of Granger causality between energy consumption and CO2 emissions using the data of 188 countries, and they only observed the unidirectional causality from energy consumption to CO2 emissions. Chang  also suggested the unidirectional causality from electricity consumption to CO2 emissions. Therefore, we can say that the increasing electricity consumption directly causes the rising of CO2 emissions, indirectly leads to higher average temperature caused by global warming, and electricity consumption further increases again directly due to the higher average temperature. In short, electricity consumption could indirectly affect temperature through climate change in the long term, and in turn, temperature can directly affect electricity consumption  . That is, we can estimate the STR models under the assumption of exogeneity. Hence, our empirical model is written as Equation (12),
where t means time dimension; represents filtered electricity consumption at time t; denotes temperature at time t; stands for the residuals with a mean of zero and constant variance; is the transition function with the transition variable . As mentioned before, the parameters of γ, c, π2, and θ2 are the key parameters of our following estimation.
4. Estimated Results
4.1. Unit Root Testing
Before performing the unit root testing, we divide the samples into six sample groups, and the sample period is four years for each sample group. Table 3 reports the results of unit testing for each data series. We can see that the series is stationary for each sample period from the results of ADF test and P.P. test. On the other hand, the results of ADF test and P.P. test for are similar to the testing results of the series which rejects the null hypothesis of unit root. This means that the series is stationary for each sample period.
When we considers the DGP of series with possible structural breaks, the impulse dummy and shift dummy are used to detect possible structural breaks in this study. The result of unit root testing for structural breaks shows that all the null hypothesis of unit root is rejected at 1% significance level no matter the series are or while the DGP of series includes impulse dummy variables for structural breaks, meaning all series are stationary for each sample period. However, when the DGP of series includes shift dummy variables for structural breaks, some reject the null hypothesis, others do not.
In sum, we conclude that both the series and the series are stationary at level based on our unit root testing; that is, we will estimate the empirical model with level data.
4.2. Estimated Result of the STR Model
Table 4 shows the results of the nonlinear model test for each sample period, F1, F2, F3 and F4 represent the statistic of F test for H1, H2, H3 and H4, respectively. Firstly, we can see that the P-value of F1 for each period all reject the null hypothesis of the linear model at 1% significance level, meaning that we should consider the nonlinear relationship between electricity consumption and temperature in our empirical model for each sample period. Secondly, as the rules of model selection mentioned before indicates, if F2 or F4 has the strongest significant P-value among F2, F3 and F4, the appropriate model will be LSTR. Moreover, the testing results indicate that the appropriate models of each sample period are all LSTR models.
Table 3. Results of unit root test.
Notes: 1) The critical values (1%, 5%, 10%) of UR test with structural break are (−3.48, −2.88, −2.58). 2) ***, **, * refers to the significance level at 1%, 5%, 10%, respectively.
Table 4. Results of appropriate nonlinear model test.
Notes: * denotes the strongest rejection among null hypothesis of H2, H3 and H4.
Table 5. The estimated results of STR models.
Notes: 1) ***, **, * refers to the significance level at 1%, 5%, 10%, respectively. 2) () refers to stand error.
Subsequently, we establish the LSTR model for six sample periods, respectively. Table 5 represents the estimated results of the STR model for each sample periods, the appropriate lag length of the dependent variable is chosen by minimum AIC (Akaike Information Criterion) with maximum lag length of 10 lags. We only report the estimated coefficient of key parameters in Table 5 to focus on threshold temperature of electricity consumption.
We can see that the estimated results of the linear part of all models indicate the positively significant relationships exist between electricity consumption and temperature in Taiwan. Furthermore, the largest coefficient of π2 is estimated during the period from 2003 to 2007 (Model 5), meaning that during these 5 years, when temperature increases by 1˚C, people are the most sensitive to consume excessive electricity compared with other periods.
In addition, for the estimated results of the nonlinear part, Model 2, Model 3 and Model 6 represents the positively significant relationships also exist between electricity consumption and temperature. It implies that when the transition variable (temperature) increases, the relationship between electricity consumption and temperature becomes much more positive.
Turning to the slope parameter γ, the estimated values are from 4.088 to 39.400, and four of six models provide the significant value of parameter γ, in addition, for the estimated results of parameter c, the estimated values are from 25.530˚C to 27.156˚C, which implies that the electricity consumption is continuously transforming with logistic function when the temperature level reaches an inflexion point. For instance, as the result of Model 6, when the temperature reaches 26.884˚C, the relationship between electricity consumption and temperature becomes more sensitive, and increases in temperature cause nonlinear increases in electricity consumption.
At this level, we can make three remarks. Firstly, the estimated value of threshold parameter c (25.530˚C to 27.156˚C) is different from the value of 18.3˚C generally used as the reference temperature of cooling degree days (CDDs) in the past literatures. Secondly, the average of all threshold value is 26.384˚C, and this value is not only close to the official reference temperature of CDDs in Taiwan, but also similar to the threshold temperature of CDDs used by Holtedahl and Joutz  (74˚F and 80˚F). Thirdly, the threshold value of temperature is not always a fixed value; the number could fluctuate in different time periods.
Moreover, we can reasonably infer that El Nino Southern Oscillation (ENSO) is one of the reasons for varying threshold temperatures in different time periods. We believe that El Niño events will lead to warmer winters and hotter summers; that is, there will be more hot days in a year when an El Niño event occurs. Needless to say, more hot days could change how people use electricity. For instance, as an El Niño event makes people feel hot in a warmer winter, people would use air-conditioners to create a comfortable indoor temperature by lowering the temperature, which leads to a lower the threshold temperature. On the contrary, a La Niña event makes colder summers and winters colder, which would also encourage people to use air-conditioner due to a higher temperature; this could lead to a higher threshold temperature.
To add this into consideration, we define ENSO score based on the data of Oceanic Niño Index (ONI) sourced from CPC  . According to the definition of ENSO events, we assign a score of 4 to −4 to distinguish a very strong El Niño year (a score 4), a strong El Niño year (a score 3), a moderate El Niño year (a score 2), a weak El Niño year (a score 1), a neutral year (a score 0), a weak La Niña year (a score −1), a moderate La Niña year (a score −2), and −3 for a strong La Niña year (a score −3). Therefore, we can see the relationship between ENSO scores and threshold temperatures from Figure 3.
4.3. Model Diagnostics
The quality of the estimated nonlinear model would be examined against misspecification like what we conducted on the linear model. Specification tests such as the serial correlation test  , the ARCH-LM test  , a normality test, a parameter constancy test and the no remaining nonlinearity test are employed for
Figure 3. Historical ENSO score and threshold temperature.
Table 6. Testing for the serial autocorrelation of residuals (F-value).
Notes: * refers to the significance level at 10%.
model diagnostic in this section. Table 6 shows the testing results of the serial correlation of residuals for each model with 8 lag terms. We can see that only the fourth lag term of Model 2 rejects the null hypothesis of no serial correlation of residuals at 10% significance level. Overall, we can say that none of six models have the serial correlation.
In addition, Table 7 reports the testing results of the ARCH-LM test, the normality test, the parameter constancy test and the no remaining nonlinearity test, respectively. We can see that all models are no problem of heteroskedasticity up to the lag order of eight as the results of the ARCH-LM test. All results of the Jarque-Bera (J.B.) normality test cannot reject the null hypothesis of normality, meaning that the error terms of all models are normally distributed.
We also check whether there is remaining nonlinearity in the models after all the STR models has been fitted, and the null hypothesis of no remaining nonlinearity test is no additive nonlinearity of the STR model. The testing results are all statistically non-significant, which means there is no additive nonlinearity in any of the STR model. The null hypothesis of parameter constancy test is constant parameters against the alternative hypothesis of smooth continuous changes in parameters. The testing results cannot reject the null hypothesis, and it indicates that the parameters of six models are constant in both regimes.
Figure 4 graphs logistic transition function of six models, and the diagram of Model 1 to Model 6 are listed in Figure 4 from left to right and from up to down. We can see that Model 1 has the steepest slope of logistic transition function among the six models, meaning the speed of its transition between two regimes is the fastest during the period from 1983 to 1988. By contrast, Model 5 has the gentlest slope of logistic transition function among the six models, showing that the speed of its transition between the two regimes is the slowest during the period from 2003 to 2007.
Table 7. Results of specification tests.
Source: Author’s estimation.
Figure 4. Logistic transition function of six models which all are with transition variable of temperature.
This study discusses the relationship between electricity consumption and temperature in Taiwan for the period from 1983 to 2012. In order to get more information from our data, we divide all samples into six groups with sample period of five years for each group before conducting the empirical estimation. Furthermore, we employ the STR model to estimate the nonlinear relationship between electricity consumption and temperature for each sample period. In addition, we also find threshold temperatures on the nonlinear relationship between electricity consumption and temperature for each sample period.
The empirical results show that there are positively significant effects of temperature on electricity consumption in Taiwan for each sample period. When we only focus on the estimated results for the linear part of the model, we can find that the purely linear effects of temperature on electricity consumption keep rising over the whole sample period. However, some estimated results for the nonlinear part of the model are positively significant, while others are not. That is, we cannot conclude that the purely total effects of temperature on electricity consumption also keep magnifying over the whole sample period.
On the other hand, we figure out the threshold temperature estimated using the STR estimation for each sample period. The threshold temperature is 26.892˚C (1983-1987), 25.530˚C (1988-1992), 25.364˚C (1993-1997), 27.156˚C (1998-2002), 26.477˚C (2003-2007), and 26.884˚C (2008-2012). In sum, the average threshold temperature over the period from 1983 to 2012 is 26.384˚C. We can say that the pure effects of temperature on electricity become much more sensitive if temperature reaches the threshold temperature based on our empirical results. In addition, Taiwan has a subtropical climate with higher humidity and a higher yearly average temperature (at about 22˚C), meaning that Taiwan has a higher temperature compared to other countries with a temperate and frigid climate. Hence, estimated temperatures in this study are between 25˚C and 27˚C, which are reasonably comfortable for people living in subtropical climate such as Taiwan. Furthermore, the estimated threshold temperature carries the same meaning with the reference temperature of CDDs. That is, if the air temperature is higher than the threshold temperature, it will lead to increases in electricity consumption. For instance, policy makers could use the threshold temperature to be the reference temperature, and thus they could propose a policy to ask people to reduce electricity consumption when the air temperature is higher than the reference temperature in order to save electricity and to promote efficiency of using electricity.
In addition, increasing use of air conditioning resulted in temperature rise for a comfortable living environment seems to be a serious problem of power supply in Taiwan under the influence of global warming. Santamouris et al.  concluded that a 1% increase in temperature under the threshold temperature of 24˚C in warm countries would lead to a 3.5% increase in peak electricity demand. That is, if the ambient temperature is higher than the threshold temperature, the risk of power shortage sharply rises. Now Taiwan is aiming for a future of a non-nuclear homeland and is actively developing renewable energy (mainly wind and solar power) to meet the estimated losses of nuclear energy. However, wind and solar power are not stable sources of power supply due to current limitations of power storage technology in Taiwan. Therefore, the threshold temperature could send a warning signal to not only promote the idea of energy saving but also suggest more power system expansion planning in preparation for more operating reserve in the future.
Thus, the estimated threshold temperature will have policy implications for policy makers, who can use the threshold temperature in this study as a reference for making electricity management policies in Taiwan.
 Ruth, M. and Lin, A.C. (2006) Regional Energy Demand and Adaptations to Climate Change: Methodology and Application to the State of Maryland, USA. Energy Policy, 34, 2820-2833.
 Parkpoom, S. and Harrison, G.P. (2008) Analyzing the Impact of Climate Change on Future Electricity Demand in Thailand. IEEE Transactions on Power Systems, 23, 1441-1448.
 Ahmed, T., Muttaqi, K.M. and Agalgaonkar, A.P. (2012) Climate Change Impacts on Electricity Demand in the State of New South Wales, Australia. Applied Energy, 98, 376-383.
 Zachariadis, T. and Hadjinicolaou, P. (2014) The Effect of Climate Change on Electricity Needs—A Case Study from Mediterranean Europe. Energy, 76, 899-910.
 Hekkenberg, M., Moll, H.C. and Schoot-Uiterkamp, A.J.M. (2009) Dynamic Temperature Dependence Patterns in Future Energy Demand Models in the Context of Climate Change. Energy, 34, 1797-1806.
 Moral-Carcedo, J. and Vicéns-Otero, J. (2005) Modelling the Non-Linear Response of Spanish Electricity Demand to Temperature Variations. Energy Economics, 27, 477-494.
 Bessec, M. and Fouquau, J. (2008) The Non-linear Link between Electricity Consumption and Temperature in Europe: A Threshold Panel Approach. Energy Economics, 30, 2705-2721.
 Lee, C.C. and Chiu, Y.B. (2011) Electricity Demand Elasticities and Temperature: Evidence from Panel Smooth Transition Regression with Instrumental Variable Approach. Energy Economics, 33, 896-902.
 TCCIP (2014) TCCIP Data Service Platform, Taiwan Climate Change Projection and Information Platform. Ministry of Science and Technology, Executive Yuan, Taiwan.
 Lanne, M., Lütkepohl, H. and Saikkonen, P. (2002) Comparison of Unit Root Tests for Time Series with Level Shifts. Journal of Time Series Analysis, 23, 667-685.
 Terasvirta, T. (1998) Modelling Economic Relationships with Smooth Transition Regression. In: Ullah, A. and Giles, D., Eds., Handbook of Applied Economic Statistics, Marcel Dekker, New York, 229-246.
 Areosa, W.D., McAleer, M. and Medeiros, C.M. (2011) Moment-Based Estimation of Smooth Transition Regression Models with Endogenous Variables. Journal of Econometrics, 165, 100-111.
 Terasvirta, T. (2004) Smooth Transition Regression Modelling. In: Lütkepoh, H. and Kratzig, M., Eds., Applied Time Series Econometrics, Cambridge University Press, Cambridge, 222-242.
 Chen, P.Y., Chen, S.T., Hsu, C.S. and Chen, C.C. (2016) Modeling the Global Relationships among Economic Growth, Energy Consumption and CO2 Emissions. Renewable and Sustainable Energy Reviews, 65, 420-431.
 Chang, C.C. (2010) A Multivariate Causality Test of Carbon Dioxide Emissions, Energy Consumption and Economic Growth in China. Applied Energy, 87, 3533-3537.
 Hao, Y., Chen, H., Wei, Y.M. and Li, Y.M. (2016) The Influence of Climate Change on CO2 (Carbon Dioxide) Emissions: An Empirical Estimation Based on Chinese Provincial Panel Data. Journal of Cleaner Production, 131, 667-677.
 CPC (2015) Climate Prediction Center, National Oceanic and Atomspheric Administration, USA. http://www.cpc.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.shtml
 Santamouris, M., Cartalis, C., Synnefa, A. and Kolokotsa, D. (2015) On the Impact of Urban Heat Island and Global Warming on the Power Demand and Electricity Consumption of Buildings—A Review. Energy and Buildings, 98, 119-124.