Modeling Temperature and Pricing Weather Derivatives Based on Temperature

This study first proposes a temperature model to calculate the temperature indices upon which temperature-based derivatives are written. The model is designed as a mean-reverting process driven by a Levy process to represent jumps and other features of temperature. Temperature indices are mainly measured as deviations from a base temperature, and, hence, the proposed model includes jumps because they may constitute an important part of this deviation for some locations. The estimated value of a temperature index and its distribution in this model apply an inversion formula to the temperature model. Second, this study develops a pricing process over calculated index values, which returns a customized price for temperature-based derivatives considering that temperature has unique effects on every economic entity. This personalized price is also used to reveal the trading behavior of a hypothesized entity in a temperature-based derivative trade with profitmaximization as the objective.Thus, this study presents a new method that does not need to evaluate the risk-aversion behavior of any economic entity.


Introduction
Temperature-based derivatives represent a new financial tool to buy and sell a natural phenomenon.Doing so requires two things: a unit of measurement for the natural phenomenon that everyone agrees upon and a price that may facilitate a transaction.This study is designed to evaluate these two requirements.
Some preexisting measures already appear in the form of indices to meet the first requirement.To find values for these indices, the literature contains several temperature models using mean-reverting processes as the main tool.The most cited study develops an Ornstein-Uhlenbeck (OU) process to model temperature [1].Using the equivalent martingale measures approach, the authors determine the price of an option.Benth and Šaltytė-Benth [2] model temperature as a continuous time autoregressive process for Stockholm and report a clear seasonal variation in regression residuals.They propose a model using a higher-order continuous time autoregressive process, driven by a Wiener process with seasonal standard deviation.While pricing futures and options, they consider a Gaussian structure in the temperature dynamics.In another study [3], they model temperature with an OU process driven by a generalized Levy process.The model contains seasonal mean and volatility.Instead of dynamic models, some authors offer time-series models to represent temperature.Campbell and Diebold [4] apply a time series approach to model temperature, including trend seasonality represented by a low-ordered Fourier series and cyclical patterns represented by autoregressive lags.The contributions to conditional variance dynamics are coming from seasonal and cyclical components.The authors used Fourier series and GARCH processes to represent seasonal volatility components and cyclical volatility components, respectively.Jewson and Caballero [5] discuss the use of weather forecasts in pricing weather derivatives, presenting two methods for strong seasonality in probability distributions and the autocorrelation structure of temperature anomalies.Elias et al.
Advances in Meteorology [6] develop four regime-switching models of temperature for pricing temperature based derivatives and find that a twostate model governed by a mean-reverting process as the first state and by a Brownian motion as the second state was superior to the others.Schiller et al. [7] and Oetomo and Stevenson [8] provide a comparison of different models.
To comply with the first requirement, the current study offers a temperature model based on Alaton et al. [1], which was defined after analyzing temperature data from different locations.The temperature model in this study is a meanreverting Levy process.The Levy part contains a Brownian motion and two mean reverting jump processes driven by compound Poisson processes.For some flexibility, the jumps are designed as slow and fast mean-reverting processes, which are independent.The main difference with the model proposed here is its inclusion of jumps.Because temperature indices are mainly calculated as deviation of temperature from a base temperature, the model assumes that jumps are inevitable, at least for certain locations.The numerical estimates in this study contain test results related to this issue.The solution to the proposed temperature model is applied inversion formula to obtain approximated expected value of a specific index type and to obtain the approximated distribution of the same index.
Notably, temperature has unique behavior for any location in which it is measured.Therefore, it is not possible to develop a single model that explains every temperature behavior in every location.In addition, more than 100,000 weather stations worldwide measure temperature for different periods.It may even be difficult to develop a temperature model that is valid for all time at a single location.Thus, this study aims to cover more locations and periods by simply using a flexible model that can include or exclude jumps.
The second requirement, temperature-based derivatives pricing, is more complicated.Because the underlying commodity is not a traded asset, weather derivatives based on temperature have an incomplete market [9].Carr et al. [10], Magill and Quinzii [11], and El Karoui and Quenez [12] provide a general discussion of incomplete markets.Pricing temperature-based derivatives is mainly based on two approaches: dynamic valuation and equilibrium asset pricing.The dynamic pricing approaches [1,2] were discussed above.For equilibrium pricing, Cao and Wei [13] use a generalization of the Jr. Lucas model [14], which considers weather as another source of uncertainty.Richards et al. [15] suggest another equilibrium model.Davis [16] uses the marginal substitution value approach for pricing in incomplete markets.In addition, Xu et al. [17] use another classification for pricing temperature-based derivatives and add actuarial pricing and extended risk-neutral valuation in addition to equilibrium asset pricing, where the former is based on Jewson and Brix [18] and the latter on Hull [19] and Turvey [20].In addition, some researchers used Monte-Carlo simulations in pricing temperature-based derivatives [21].
This study bases its pricing on the monetary effect of the natural phenomenon on economic entities.Further, this study shows that temperature has different effects on different entities.The same temperature may have a positive effect on one entity and a negative effect on another and is therefore personal, ceteris paribus.Thus, the study develops a personal price, which may require a determination of the entity's risk aversion behavior.To address this problem, this study focuses on entity-specific trading behavior rather than the entity's risk aversion behavior in order to develop a more realistic approach by avoiding an inconclusive debate over the risk premiums and utility functions used to calculate risk premiums.Moreover, a benefit of the proposed pricing model is that it is independent from how researchers measure temperature.
Critics may object to the move from a stochastic temperature model to some form of actuarial pricing model.There are several reasons for this move: first, this study demonstrates that risk-neutral pricing ends up with superhedging; second, the discussion about risk premiums in the literature is unclear; and, finally, the calculations of jump processes needed approximations to obtain certain results.These considerations led to this study's development of a more appropriate and practical method.
The second section of the paper provides the approximated index calculation and distribution of temperature after presenting a temperature model.Third section develops individualized prices and discusses the trading behavior of a hypothetical entity.The paper then presents the study's conclusions.

Model
Some basic terminology is defined in the following: where  is the daily average temperature,  represents a certain day, and  max  and  min  are the maximum and minimum temperatures of the given day, respectively.
For the Heating Degree Day (HDD) temperature index used in temperature-based derivatives, for a given day, where Base is a predetermined temperature level and   is the average temperature calculated as in (1) for a given day .
Cumulative HDD (CHDD): where HDD  is calculated as in (2) and  is the time horizon, which is generally a month or a season.
2.1.The Temperature Model.Based on Alaton et al. [1], the temperature model is an OU process driven by a Levy process that contains independent processes as Brownian motions and two mean-reverting compound Poisson processes.The model is represented as follows: where    is a cyclical process of temperature and represented in (5).Additionally,  is the mean-reversion parameter, and subscript  represents time.
=  +  +  sin ( + ) where  = 2 365 ,  is the phase angle. ( The differential of the driving Levy process   is defined as follows: The Brownian component of   will be approximated by the ARCH (1) model.To represent the different jump structures in temperature in the form of a single jump and a series of jumps,   and   are defined as fast and slow mean-reverting OU processes driven by compound Poisson processes with intensities of   and   , and  and  being mean-reversion parameters, respectively.Hayfavi and Talasli [22] use a somewhat similar mean-reverting jump process combination in their model of spot electricity prices.
The solutions to these non-Gaussian processes are [23] the following: The solution to (4) is given as To find the value of a temperature-based derivative, one needs the distribution of the underlying temperature given in (9).However, this does not have a closed-form solution.One way to address this problem is to use a characteristic function of the temperature and apply inversion techniques to find the value of an HDD, an approximated distribution of CHDD, and an approximated distribution for temperature itself.

Characteristic Function of
Temperature.This study follows Cont and Tankov [23] to find the characteristic function given in (9).First, using  1 , the characteristic exponent of (6) will be determined, where characteristic exponent () is defined as   =  () .The solution to  1 is Then, the characteristic exponents of the Levy components will be where  =  (∫ and for jump processes () can be written similarly.It is not possible to evaluate the integral in (12).Consequently, the following approximation method was developed.
Then, by using linear approximation, Again,   () can be written similarly.Finally, the characteristic function of the temperature model can be written explicitly.Referring to Cont and Tankov [23], In explicit form,

HDD and Distribution
Function.This part of the study focuses on measuring HDDs.It is easy to apply the calculations into other types of indices.In the current case, inversion techniques will be used to find the value of an HDD and its distribution and hence the CHDD values.

Approximating Density Function of Temperature.
To find an approximating density function of temperature, inversion formula will be applied to the characteristic function of the temperature defined in (15).Before applying the inversion formula, the following shortcuts are derived from (15).Let () and 0() be the density function and characteristic function of temperature, respectively.
Then, by inversion formula Because weather derivatives are defined on CHDDs, their distribution will be defined.Clearly, HDD values may show autocorrelation.In addition, due to the nature of the proposed temperature model in terms of the independence of the included processes and motivation to keep the process simple, the model assumes the independence of the HDDs.With this assumption, the approximated distribution of CHDD can be found using CHDD ∼  ( −  −  * , ) . (21)

Measuring HDD.
HDDs are clearly contingent claims on how temperature deviates from a base temperature.As one way to find the expected value of an HDD, this study will first find its Fourier transform.Then, the inverse Fourier transform will be applied to both the HDD's Fourier transform and the characteristic function of temperature [24].
Let  =   , () is HDD's payoff function given in ( 2), Base = , and ŵ( Now, the inversion will be applied to ŵ()0  (−), where ŵ() is defined in (22) and 0  is the characteristic function defined in (15).Let temperature in (9) be defined in shorthand notation as   =  * + Λ  , where  * is defined as in (16) and The characteristic function of Λ  can be obtained from (15) and written as 0 where  represents expectations However, it was not possible to evaluate this integral analytically, and thus the elliptic package of the  statistical software package [25] was applied to evaluate the integral numerically.
The results indicated that the integral is equal to  −  −  * ; therefore, 2.3.3.Numerical Estimates.The success of the proposed temperature model and (24) were tested in terms of forecasting Cooling Degree Day (CDD), which is another index based on temperature, and HDD values for the 12 cities listed in Tables 1 and 2. CDD is calculated as CDD  = max (0,   − Base), where Base is a predetermined temperature level and   is the average temperature calculated as in (1).Cumulative CDD (CCDD) is calculated using ∑  =1 CDD  , where CDD  is calculated as in the previous sentence and  is the  Design.A simulation was designed to run in two dimensions: the first on different cities and the second to capture changes in the parameters through time for each city.In this respect, the initial simulations used the last 5 years of data.They were then continued by including one more year of data to the existing data in each turn where HDD and CDD estimates were calculated.A turn consisted of 10,000 simulation runs.Then, the simulation results were compared with actual HDD and CDD values.Finally, the parameters of the year that offered the best estimates of the HDDs and CDDs were chosen for use in the one-year-ahead predictions.
Discretization.Discretization was done using Euler approximation [26][27][28] as follows: where  +1 = √ 0 +  1     ,  0 and  1 are being ARCH parameters and   ∼  (0, 1) . ( Parameter Estimation.Parameters were estimated as defined in [1,27,29]. Simulation of Jumps.Because they have a different structure, jumps were simulated separately and the results added to the discretized model.For this aim, the jumps were detected by first removing the mean from the actual data and selecting the values above two standard deviations.These jumps were then separated into two categories: single and multiple jumps to constitute fast and slow mean reverting jumps, respectively.The sample means and sample standard deviations of these two jump groups were found to simulate jump sizes.In addition, the intensities were found by λ = no. of jumps/no. of observations.The simulation consisted of 10,000 runs, during which the jump times were first found by using intensities.Then, for each run, random draws were taken using means and standard deviations obtained from data.Finally, the discretized jumps were added to the discretized jump model.After these phases, one-year-ahead predictions were conducted along with three other models: the Campbell and Diebold [4] model (Campbell Model), the Benth and Šaltytė-Benth [2] model (Benth Model), and the Historical Burn Analysis (HBA) model that calculates historical averages.In addition, HDD predictions were calculated based on (24).This yielded the results in Tables 1 and 2. In the Appendix, the parameter estimates for HDD and CDD values are presented, in Tables 3 and 4, respectively, including the time period used to obtain these parameters.In addition, statistical test values of these parameters are shown in Tables 5 and 6.

Analysis of Numerical
Estimates.The best estimates of HDDs and CDDs were obtained for different periods as shown in the Appendix.This is mainly a characteristic of the temperature since it changes in its long-term behavior.However, this may not be a good way to use all existing data for a city.Instead, every location must be scanned and evaluated for different time periods to obtain the best prediction results.In addition, the current model is equally successful in HDD and CDD predictions.
Having a good estimate of HDD and CDD values does not necessarily correspond to the best fit of the model to temperature data.This may be because including the jumps may result in a better estimation of index values while deteriorating the fit of the model to the data.
The current model demonstrates its capacity related to the changing conditions of the temperature data.For example, in Chicago, both jump types were statistically significant and the model predicted HDDs accurately.On the other hand, Tokyo did not have any jumps during entire period, and the current model was still able to make accurate predictions for HDDs.Interestingly, the Historical Burn Analysis containing long-term HDD and CDD averages provided successful predictions.These results suggest that temperatures do not change significantly for certain locations.
As expected, the approximated HDD calculations obtained from (24) were less accurate than the simulations.Nevertheless, the predictions based on the equation were still successful.The estimated HDDs of Los Angeles and Washington were better than any other model.Finally, while this study conducted comparisons for only 12 stations, there are 125,000 weather stations worldwide.It is therefore impossible to say that one model is superior to the others, though the current model and (24) were successful in certain locations and periods, which merits evaluation.

Pricing
This section addresses temperature risk and its differences from classical asset risks before providing a fair price for a temperature-based derivative written of HDD with real probabilities.Next, this section shows that the price of the derivative will be super-hedging when using risk-neutral probabilities.Finally, a personal price will be developed based on personal temperature risk.3.1.Temperature Risk.In simple terms, temperature risk is volume risk that affects sales [13].This study will first focus on a single company with an obvious exposure to temperature to reveal the relationship between temperature index and the company's sales.To maintain the focus on the relationship between sales and temperature, a very simple linear model that omits other possible factors that may affect sales will be introduced.As a candidate company, consider a retail gas seller concerned about sales and profit in the following January.The company's risk will be measured through the effect of sales on profits.Consider the following: ES =  + CHDD where ES is expected sales.
A positive value for  in (26) indicates that the company can sell a certain amount of its products, even in the case of zero CHDD.Similarly, a positive value for  in the same equation represents a positive relationship between sales and CHDDs.
The relationship between CHDDs and profit can be established using cost and revenues.
To guarantee a positive profit after a certain amount of CHDDs is realized, assume that Price > 1 and constant.
The aim is to construct a relationship between CHDD and profit functions to ascertain the magnitude of the effect of temperature on the company.Figure 1  relationship, with the assumption that  >  * Price initially provides negative profits.
Figure 1 constructs a relationship between CHDD and profit such that every CHDD value now represents a monetary value in terms of positive and negative profits. Figure 1 shows a deterministic relationship between CHDD level and profit.The CHDD value for the proposed period is unknown.Now, the company's risk can be defined as having a low value of CHDD for a certain period that leads to a loss.In other words, the company's risk will fall left of CHDD 1 and will be covered by the  Line.This possible loss will be referred to as temperature risk (TR) and is equal to the area covered by 0 CHDD 1 (Price − 1) −  in Figure 1.This is actually the total TR (TTR) and can be realized if CHDD for the next January becomes zero.In reality, the CHDD for the next January is not known because it is a stochastic value.Therefore, TR can be written as where chdd is unknown value of CHDD for a certain period.(30) It is clear that the magnitude of TR will depend on business type and size.For example, for a gas company, a decrease in the index value means a lower value in sales.However, at the same time, this decrease may result in an increase in the sales of a beverage company.The magnitude of the decrease or the increase in profits, on the other hand, will be directly related to the size of the business.

An Approximated Fair Price of a Temperature Based Put
Option.Under linear approximation, the temperature has a normal distribution with mean  * and variance , as defined in (18) and (19) The distribution function is given in (21).Let  be the strike value,  = CHDD, let () be the probability density function of CHDD, and tick the monetary value for each Equation ( 32) provides an actuarial price since it is based on expected values.Derivative pricing using real probabilities may require an evaluation of risk premiums.The literature offers only inconclusive discussions.For example, Hull [19] states that it is possible to calculate the payoffs of weather derivatives with real probabilities because these derivatives have no systemic risk.Turvey [30] supports this idea.On the other hand, Chincarini [31] examines the efficiency of weather futures in CME in HDD and CDD futures assuming an efficient market and risk premiums varying from negative to positive values across cities.In addition, Cao and Wei [13] highlighted the importance of the market price of risk for weather derivatives.This study will follow a different path by evaluating the trading behavior of a candidate company instead of the company's risk behavior.The next section will show the outcome of using risk-neutral probabilities in (32).

Risk-Neutral Pricing.
In this part, the put option of the previous section will be priced using risk-neutral probabilities in a binomial model.CHDD is clearly composed of HDDs.A closer look at HDDs reveals that realizations of HDDs can be represented by a binomial model such that if HDD is realized, there will be an upward movement as in the binomial model; otherwise, in case of a downward movement, there will be 0. The calculations were done according to Björk [32].For the up movement, the best candidate for HDD will be the mean value of an HDD calculated using (24).Let  = (24).Figure 2 represents the binomial model.
The probability of an upward movement,   , is found by ∫  −∞ (), while the probability of downward movement is equal to   = 1 −   .Let risk-free rate  = 0.The model satisfies the condition of being arbitrage-free in the form of  ≤ 1 +  ≤  by definition.The martingale measure for the Definition 2. When [max( − , 0)]( −  +  * ) = ,  is the general price, which is valid for any economic entity.
The above setup can be extended by the concept of a shadow price, which describes the effects of the resources in a production process on profit.In the current case, the profit function is a deterministic function that is a payoff function of the index value CHDD.Therefore, CHDD can be seen as the resource that produces the profit.It is possible to measure the effect of one unit of change in CHDD on profit and use it as the price of one unit of CHDD.Again, assuming CHDD = , consider the following: The value given in (39) is a good candidate for the tick value mentioned in (32).Now, using ( 32) and ( 39), a new definition can be given.From here, a connection between the current approach and the utility approach can be established.Since the current setup is based on profit maximization, it coincides with the utility approach based on wealth maximization.The gain is that this statement is true for any utility function choices.
A numerical one-day-ahead estimate of temperature for the price of an HDD for Ankara was developed using (33) and (32).The mean and standard deviation of the approximated distribution were calculated according to (18) and (19).The value of  in (19) was approximated by conditional variance of the ARCH model.The tick value was taken as $1, and the strike value  was taken as an interval from 65 to 100.The estimated values are shown in Figure 4.
As mentioned earlier, (33) defines a general price and trading behaviors for any company since the equation does not include the profit function.This generality does not give much insight into what a put option with a strike value of  actually means for a specific company.This deficiency was corrected by replacing the tick value with (40) because, unlike ordinary assets, temperature affects economic entities on different scales.Thus, a personalized price must apply for each economic entity.Moreover, (40) and Definition 3 state that the hypothetical company will enter a trade for the option if there is a possibility for arbitrage.If the fair price is available, the company will be indifferent to entering a trade or doing nothing.In the new pricing setup, the expected profit will always be the maximum, as will be the company's utility.Risk aversion will negatively affect the maximum profit and utility.Therefore, having a maximum profit and the resulting maximum utility will direct the company to follow the presented approach rather than the one of suboptimal risk aversion.

Conclusions
Derivatives written on temperature are based on index values obtained from temperature data, which are essentially measured as deviations of temperature from a threshold value.This makes measuring deviations from a base temperature in the form of jumps important for any temperature model for some locations.This study proposed and demonstrated a temperature model that included different kinds of jumps that was then handled using different techniques.In addition, unlike existing models that consider temperature risk as the result of the temperature itself, like in stocks, the proposed model shows that financial risk caused by temperature differs from classical asset risk, though this risk depends on the type of business.This study demonstrated a method to measure this temperature risk.Moreover, almost all of the existing pricing methods are based on risk-neutral valuations.The results from this study showed that risk-neutral valuation in temperature-based derivatives ends with super-hedging.Finally, the current study offers a pricing scheme that differs from classical pricing approaches that are based on risk-neutrality or risk-aversion concepts.Instead of utility functions, this study employs a more realistic and practical approach in terms of objective functions set by the firm itself.
In return, the model provides a personalized price based on company-specific temperature risk to realize an objective in terms of profit.

Table 1 :
One-year-ahead prediction and error values for HDDs.Number of years of data included in estimations.6 Advances in MeteorologyTable 2: One-year-ahead prediction and error values for CDDs.
Number of years of data included in estimation.Advances in Meteorology 7 time horizon, which is generally a month or season.The test assumes a base temperature of 18 degrees Celsius and proceeds in the following manner:Data.Temperature data is available for 12 cities covering 38 years from 1974 to 2011.The first part uses 37 years of data to estimate the parameters.The temperature data for 2011 were used to compare with one-year-ahead predictions.Temperature data were obtained from the National Climatic Data Center.

Figure 2 :
Figure 2: Construction of a CHDD tree using a binomial model.

Figure 4 :
Figure 4: Estimated value of an HDD for Ankara.

Table 3 :
Parameter estimation for HDD calculations.Year represents amount of data in years that produced best estimation results. *

Table 4 :
Parameter estimation for CDD calculations.Year represents amount of data in years that produced best estimation results.* * *

Table 5 :
values of parameters for HDD calculations.

Table 6 :
values of parameters for CDD calculations.
NS: not statistically significant.NA: not available.