The Implied Volatility Surface as a Decision-Support Framework for Systematic Cash-Secured Put Strategies

Frédéric Valognes

In this article, Frédéric VALOGNES, lecturer, author and Certified European Financial Analyst (CEFA®), examines whether the dynamics of the implied volatility surface may provide a decision-support framework for systematic cash-secured put strategies.

Abstract

The Black-Scholes-Merton model remains one of the most influential developments in modern financial economics. Whilst its mathematical formulation continues to provide the benchmark for pricing European options, one of its central assumptions — namely that volatility remains constant throughout the life of an option — is persistently contradicted by observed market prices.

Rather than constituting a weakness of the model, these discrepancies reveal valuable information regarding investors’ expectations, market sentiment and the pricing of downside risk. The resulting volatility skews and smiles have therefore become essential components of both academic research and professional option trading.

This paper argues that the implied volatility surface should not be viewed solely as a pricing adjustment. Its geometry and, more importantly, its evolution over time may provide additional information capable of assisting investment decisions. Attention is devoted to cash-secured short put strategies, for which the level of implied volatility alone frequently proves insufficient.

Drawing upon preliminary observations obtained from listed CAC 40 index options across several maturities, the article explores whether the dynamics of the implied volatility surface may constitute a useful decision-support indicator. Rather than proposing a predictive pricing model, the objective is to examine whether changes in the shape, slope and term structure of implied volatility can contribute to a more disciplined framework for identifying favourable market environments in which to initiate systematic cash-secured put strategies.

Introduction

Within option markets, implied volatility occupies a rather singular position. Originally introduced as the unknown parameter required to reconcile observed option prices with the Black-Scholes-Merton valuation model, it has progressively evolved from a purely technical pricing input into one of the most closely monitored indicators in financial markets. Today, implied volatility is commonly interpreted not simply as a pricing parameter, but as a market-based measure of uncertainty, reflecting the aggregate expectations of thousands of market participants.

For investors employing cash-secured short put strategies, that is, selling put options while maintaining sufficient cash reserves to purchase the underlying asset if assignment occurs, implied volatility plays an obvious practical role. Higher implied volatility generally translates into higher option premiums, thereby increasing the potential income associated with selling options. This simple observation has encouraged many practitioners to associate elevated implied volatility with favourable selling opportunities.

Experience, however, suggests that such a conclusion is frequently incomplete. Periods characterised by exceptionally high implied volatility often coincide with episodes of considerable financial stress, during which uncertainty continues to increase and option premiums expand further. Entering short option positions solely because implied volatility appears elevated may therefore expose investors to significant mark-to-market losses before market conditions eventually stabilise.

The question addressed in this article is therefore slightly different.

Rather than asking whether implied volatility is high, it may be more appropriate to ask whether the behaviour of the implied volatility surface itself contains additional information capable of assisting investment decisions.

More specifically, can the dynamics of the implied volatility surface, particularly the evolution of the downside volatility skew, provide useful information regarding changing market conditions? Since deep out-of-the-money put options typically incorporate a substantial premium reflecting institutional demand for portfolio insurance, does a progressive flattening of the skew signal that market stress is easing while option premiums remain comparatively attractive?

If such behaviour can be observed consistently, the volatility surface ceases to be merely an output of an option pricing model. Instead, it becomes a potential decision-support framework, capable of complementing more traditional criteria such as premium level, strike selection or time to maturity.

The purpose of the present article is not to challenge the theoretical foundations of the Black-Scholes-Merton model. On the contrary, the model remains indispensable, since implied volatility itself is extracted from its pricing equation. The objective is rather to investigate whether the systematic departures observed between theoretical assumptions and market prices may themselves convey exploitable information through the dynamics of the implied volatility surface, thereby supporting decisions regarding option selection, strike prices, market conditions and the implementation of systematic cash-secured put strategies.

Figure 1. Transfer of Risk Between Option Buyer and Option Seller

Figure 1. Options transfer market risk between two counterparties with fundamentally different expectations. Whilst the buyer acquires protection against adverse price movements, the seller receives an option premium in exchange for assuming the corresponding contingent obligation. This transfer of risk constitutes the economic foundation upon which option markets operate and explains the central role played by option premiums in systematic short-put strategies.

The following sections revisit the theoretical foundations of implied volatility before examining why market observations systematically depart from the assumptions of constant volatility. Attention is subsequently devoted to the informational content embedded within volatility skews and smiles, leading to the introduction of a practical analytical framework intended to investigate whether changes in the implied volatility surface may contribute to the identification of favourable environments for systematic cash-secured put-selling.

The Black-Scholes-Merton Framework: An Elegant Model Built upon Simplifying Assumptions

Since its publication in 1973, the Black-Scholes-Merton model has become one of the most influential achievements in financial economics. Beyond providing a closed-form solution for the valuation of European options, it established a rigorous mathematical framework linking derivative prices to the stochastic behaviour of the underlying asset. More than half a century later, despite the emergence of increasingly sophisticated numerical models, Black-Scholes remains the common language of option markets.

Its enduring success stems from the remarkable intuition underlying the model. Rather than attempting to forecast future prices directly, Black-Scholes demonstrates that an option may be replicated through a continuously adjusted portfolio combining the underlying asset and a risk-free investment. Under a specific set of assumptions, this replication argument leads to a unique theoretical option value independent of investors’ individual expectations.

These assumptions are well known. Asset prices are assumed to follow a geometric Brownian motion with constant volatility. Markets are perfectly liquid and frictionless, allowing continuous trading without transaction costs or taxes. Interest rates remain constant throughout the life of the contract, whilst European options can only be exercised at maturity. Finally, market participants are assumed to behave rationally and possess homogeneous expectations.

From a practical perspective, few of these assumptions are fully satisfied in real financial markets. Transaction costs exist, volatility varies continuously, liquidity fluctuates and investors frequently react in heterogeneous ways to new information. Nevertheless, the model remains extraordinarily useful because it provides a coherent reference framework from which market observations may subsequently be interpreted.

One of its most significant contributions lies in the concept of implied volatility. Rather than treating volatility as an observable market variable, the Black-Scholes equation can be solved inversely. By inserting the observed option premium together with the remaining market parameters, it becomes possible to determine the level of volatility required for the theoretical model to reproduce the market price exactly. This inferred quantity is known as implied volatility.

Implied volatility therefore represents considerably more than a simple mathematical parameter. It embodies the level of uncertainty collectively embedded within option prices by market participants. Every quoted option premium implicitly reflects the market’s assessment of future price variability, making implied volatility one of the most informative indicators available to option traders.

Yet an important observation immediately follows. If the assumptions of the Black-Scholes model were perfectly satisfied, every option sharing the same maturity would exhibit the same implied volatility, irrespective of its strike price. Reality tells a rather different story.

Figure 2. Call and Put: The Economic Foundations of Option Contracts

Figure 2. A call option grants its holder the right, but not the obligation, to purchase the underlying asset at a predetermined strike price. Conversely, a put option grants the right to sell the underlying asset under identical contractual conditions. In both cases, the buyer acquires a right by paying an option premium, whilst the seller receives that premium in exchange for assuming the corresponding contingent obligation.

Implied volatility: From a Single Parameter to a Market Indicator

The original formulation of Black-Scholes implicitly assumes that volatility constitutes a characteristic of the underlying asset itself. If this were strictly true, every option written on the same asset and sharing an identical maturity would produce the same implied volatility once observed market prices are introduced into the valuation equation.

Empirical evidence has demonstrated otherwise. When implied volatilities are computed across a range of strike prices, they rarely remain constant. Instead, they exhibit systematic patterns whose shape varies according to both the underlying asset and prevailing market conditions. These observations, initially regarded as anomalies, have gradually become recognised as fundamental characteristics of option markets. The discrepancy is not accidental. It reflects the collective behaviour of investors rather than any mathematical imperfection within the pricing equation itself.

Institutional investors, pension funds and asset managers frequently purchase out-of-the-money put options to protect equity portfolios against severe market declines. This persistent demand for downside insurance increases put premiums relative to those predicted under constant volatility assumptions. Consequently, implied volatilities extracted from these option prices become progressively higher as strike prices decrease.

The resulting asymmetry gives rise to what practitioners commonly describe as the volatility skew. Rather than representing a flaw in Black-Scholes, the skew reveals how financial markets collectively price extreme downside events. It therefore provides direct insight into investors’ perception of risk, their appetite for protection and the relative scarcity of option sellers willing to assume such exposure.

Viewed from this perspective, implied volatility ceases to be merely an intermediate calculation. It becomes a market variable, capable of conveying valuable information regarding the balance between fear and confidence prevailing amongst market participants.

From the Volatility smile to the Volatility skew

When implied volatilities are calculated across a range of strike prices for a given maturity, the resulting profile rarely corresponds to the horizontal line predicted by the Black-Scholes-Merton model. Instead, distinct empirical patterns emerge according to both the underlying asset and prevailing market conditions.

The earliest observations concerned currency and commodity options, where implied volatility frequently followed a symmetrical U-shaped profile. Deep in-the-money and deep out-of-the-money options exhibited higher implied volatilities than contracts whose strike prices were close to the prevailing market price. This phenomenon rapidly became known as the volatility smile, reflecting the characteristic curvature obtained when implied volatilities were plotted against strike prices.

The market crash of October 1987 marked a decisive turning point in option pricing. Following the unprecedented decline in global equity markets, practitioners observed that the Black-Scholes-Merton assumption of constant volatility no longer matched market prices. Implied volatilities began to differ substantially across strike prices, particularly for downside put options, reflecting investors’ increased demand for protection against extreme losses. Rather than attempting to force market prices into a single volatility parameter, traders progressively adopted the implied volatility surface itself as the practical input for option valuation. Since then, the smile and, even more prominently, the volatility skew have become standard features of option markets and indispensable tools for pricing, hedging and risk management.

Although initially regarded as an anomaly, the volatility smile gradually became recognised as a natural consequence of market behaviour rather than a failure of financial theory. Financial returns do not follow the perfectly lognormal distribution assumed by the Black-Scholes-Merton framework. Instead, empirical distributions exhibit heavier tails, occasional jumps and varying degrees of asymmetry, all of which contribute to systematic differences in implied volatility across strike prices.

Equity index options, however, generally display a markedly different pattern. Rather than producing a symmetrical smile, implied volatility typically increases as strike prices decrease. Conversely, call options with higher strike prices tend to exhibit progressively lower implied volatilities. The resulting profile no longer resembles a smile but rather a downward-sloping curve commonly referred to as the volatility skew.

This asymmetry is far from accidental. It reflects the structural demand for downside protection that characterises modern equity markets. Pension funds, insurance companies, institutional asset managers and other long-term investors regularly purchase out-of-the-money put options to protect diversified equity portfolios against severe market downturns. Such contracts effectively operate as insurance policies against extreme market events.

As demand for these protective puts increases, their market prices rise beyond the levels predicted by constant-volatility models. Once these prices are translated back into implied volatilities through the Black-Scholes equation, lower strike prices systematically exhibit higher implied volatility. The volatility skew therefore represents considerably more than a graphical curiosity. It provides a direct visual representation of how financial markets collectively price downside risk.

Rather than indicating that the Black-Scholes model has failed, the skew demonstrates that investors attribute different probabilities to upward and downward market movements. In practice, the cost of insuring against a sharp decline is significantly greater than the cost of participating in an equally pronounced upward movement. For option sellers, this distinction is of particular importance.

The additional premium associated with out-of-the-money put options constitutes the primary source of return for many systematic short-put strategies. Yet this additional premium simultaneously reflects the market’s perception of elevated downside risk. The option seller is therefore continuously confronted with a fundamental trade-off: richer premiums are generally accompanied by greater uncertainty.

Understanding this relationship represents the first step towards interpreting implied volatility not merely as a pricing parameter, but as a genuine source of market information.

Figure 3. Black-Scholes-Merton Model with Continuous Dividend Yield

Figure 3. Under the Black-Scholes assumption of constant volatility, implied volatility should remain identical across strike prices. Empirical observations reveal two distinct market structures: the volatility smile, historically observed in several currency option markets, and the downward volatility skew that characterises most equity index options.

The Volatility skew as a Measure of Collective Risk Perception

Traditional option pricing theory treats implied volatility as a parameter required to value derivative contracts. Market practitioners increasingly adopt a rather different perspective. For many traders, implied volatility has progressively become an observable market variable.

Its level reflects the price investors collectively assign to uncertainty, whilst its distribution across strike prices reveals how that uncertainty is allocated between favourable and unfavourable market scenarios. This distinction is fundamental.

If all future price movements were regarded as equally probable, the volatility surface would remain broadly symmetrical. The persistent existence of a downward skew instead demonstrates that investors consistently attribute a greater economic significance to adverse market movements than to equivalent upward fluctuations. In this respect, the volatility skew may be interpreted as a continuously updated measure of collective risk aversion.

Unlike conventional market indicators, which frequently rely upon historical observations, implied volatility incorporates forward-looking expectations embedded directly within option prices. Every transaction reflects the judgement of buyers and sellers regarding future uncertainty. The resulting volatility surface therefore aggregates thousands of independent market assessments into a single observable structure. From the perspective of a systematic put seller, the implications are immediate.

Periods during which the skew becomes exceptionally steep frequently coincide with heightened demand for downside protection. Conversely, a gradual flattening of the skew may indicate that the market is beginning to reassess the likelihood of extreme adverse scenarios.

The central hypothesis explored throughout the remainder of this article is based precisely upon this observation. Rather than considering implied volatility in isolation, greater attention may usefully be devoted to the evolution of the entire volatility surface.

Looking Beyond Implied volatility: Can the Volatility surface Become a Decision-Support Tool?

For most option practitioners, implied volatility is primarily regarded as a pricing variable. Whether calculated directly from market quotations or displayed by professional trading platforms, it is generally interpreted as a measure of the market’s expectation of future uncertainty. Consequently, trading decisions often rely upon a relatively simple observation: higher implied volatility produces higher option premiums.

For investors writing cash-secured puts, this relationship is naturally attractive. Selling options during periods of elevated implied volatility allows the collection of larger premiums whilst maintaining identical contractual obligations. Yet this apparent advantage immediately raises a practical difficulty.

Periods characterised by elevated implied volatility rarely occur in isolation. They are frequently associated with deteriorating market sentiment, increasing downside risk and heightened investor demand for protection. In such circumstances, high option premiums merely compensate sellers for assuming substantially greater uncertainty. The absolute level of implied volatility therefore provides only a partial description of market conditions. A more informative question may instead concern the behaviour of implied volatility itself.

Is the volatility surface continuing to steepen? Has it reached a plateau? Or has it begun to return progressively towards more stable market conditions?

These questions introduce an important distinction between two different approaches to option selling. The first consists simply of identifying expensive options based on their implied volatility. The second seeks to determine whether market conditions themselves have begun to evolve in favour of the option seller. The distinction is subtle but potentially significant.

A market characterised by high implied volatility, and an increasingly steep volatility skew reflects persistent demand for downside protection. Under such circumstances, option premiums may continue to increase despite already appearing historically elevated.

Conversely, if implied volatility remains relatively high whilst the overall structure of the volatility surface begins to normalise, market expectations may be undergoing a gradual transition. Although uncertainty remains elevated, the balance between buyers and sellers of protection may already be changing.

From the perspective of a systematic option seller, such an environment appears fundamentally different. The option premium remains attractive, yet the dynamics of market expectations may already be evolving towards greater stability. This observation forms the central hypothesis explored in the present work.

Rather than evaluating implied volatility solely through its absolute level, the proposed approach investigates whether the progressive normalisation of the implied volatility surface may itself constitute useful information capable of assisting the timing of cash-secured short put strategies.

Importantly, this hypothesis should not be interpreted as an attempt to forecast future market prices. No volatility model can predict future market movements with certainty. Instead, the objective is considerably more modest.

The purpose is to investigate whether the collective information continuously embedded within option prices can be organised into a coherent analytical framework capable of improving the selection of favourable option-selling environments.

Three Market Environments for Systematic Put Selling

Figure 4. The proposed framework focuses less on the absolute level of implied volatility than on the evolution of the volatility surface itself. A gradual normalisation of the skew whilst option premiums remain comparatively elevated may provide a more favourable environment for initiating systematic cash-secured put positions.

Towards a Decision-Support Framework Based on Volatility surface Dynamics

The preceding discussion naturally raises a practical question: if the geometry of the implied volatility surface reflects the collective assessment of market risk, can its evolution also provide useful information regarding the timing of option-selling strategies?

This question forms the starting point of the present investigation. Rather than considering implied volatility as a static variable observed at a single point in time, the proposed framework examines the volatility surface as a dynamic structure whose characteristics evolve continuously in response to changing market expectations. The distinction is important.

Most market participants focus primarily on the absolute level of implied volatility. Elevated implied volatility is generally interpreted as an opportunity to collect richer option premiums, whilst low implied volatility often discourages option-selling strategies. Such reasoning, however, overlooks an essential aspect of market behaviour.

Two market environments may exhibit comparable average implied volatilities whilst reflecting fundamentally different underlying conditions.

In the first case, implied volatility may still be increasing, accompanied by a progressively steeper volatility skew and a persistent demand for downside protection. In the second, implied volatility may remain elevated, but the volatility surface itself may already be beginning to stabilise, suggesting that market participants are gradually reassessing the probability of extreme downside events.

From the perspective of a systematic put seller, these two situations should not necessarily be regarded as equivalent. Although option premiums may appear imilarly attractive, the evolution of collective market expectations differs substantially.

The working hypothesis explored throughout this study is therefore deliberately modest. Rather than attempting to predict future market prices, the objective is to determine whether the progressive normalisation of the implied volatility surface may provide additional information capable of assisting the selection of favourable market environments for initiating cash-secured short put positions.

In this respect, the volatility surface is not viewed as a forecasting instrument. Instead, it is interpreted as a continuously updated representation of market sentiment whose evolution may contribute to a more disciplined investment process.

Decision-Support Framework

Figure 5. General workflow of the proposed analytical framework. Market option prices are first converted into implied volatilities using the Black-Scholes-Merton model. The resulting volatility surface is subsequently analysed through a series of descriptive indicators before being interpreted within a decision-support framework for systematic cash-secured put strategies.

Methodological Approach

The methodology developed in this work follows a sequence of analytical steps intended to transform raw market quotations into interpretable market indicators.

The process begins with the systematic collection of listed option prices for a given underlying asset and maturity. Preference is given to highly liquid option contracts to minimise distortions resulting from wide bid-ask spreads or infrequent trading activity.

Observed market premiums are then converted into implied volatilities through the inverse application of the Black-Scholes-Merton pricing equation. Once computed across the available strike prices, these implied volatilities collectively define the observed volatility surface for the selected maturity.

Rather than analysing each implied volatility independently, several global characteristics of the surface are examined simultaneously.

Attention is devoted to:

  • the overall level of implied volatility;
  • the slope of the volatility skew;
  • the degree of cross-sectional dispersion across strike prices;
  • the temporal evolution of these characteristics between successive market observations.

The purpose of this multidimensional approach is to characterise market conditions more comprehensively than would be possible through the observation of implied volatility alone. Naturally, not all option markets exhibit comparable behaviour.

The preliminary investigations presented in this article suggest that market liquidity and option maturity play a decisive role in determining the regularity of the resulting volatility surface. Highly liquid equity index options with medium- to long-term maturities appear particularly well suited to this type of analysis, whereas shorter maturities or less actively traded underlying assets may generate substantially noisier implied volatility structures.

These observations should not be interpreted as definitive conclusions. Rather, they provide an empirical motivation for the exploratory analyses presented in the following section.

Methodological Approach

Figure 6. Illustrative workflow describing the successive stages of the proposed methodology: market data acquisition, implied volatility computation, volatility surface construction, statistical charac-terisation and decision-support interpretation.

Preliminary Empirical Observations

The analytical framework presented above was subsequently applied to listed option data to examine whether the proposed interpretation of the implied volatility surface could be observed under actual market conditions.

At this stage, the objective was not to perform an exhaustive statistical validation of the methodology. Rather, the purpose was to investigate whether the dynamics of the implied volatility surface exhibited sufficiently regular behaviour to justify further quantitative analysis.

Several option chains were therefore examined, covering different underlying assets and maturities.

Attention was devoted to the CAC 40 index, whose option market offers a high level of liquidity across a broad range of strike prices. Additional observations were conducted on selected individual equities to assess the robustness of the approach under different market conditions.

The first observation concerns the influence of option maturity.

Short-dated options, particularly those approaching expiration, frequently generated irregular implied volatility profiles. Individual quotations occasionally produced local distortions, whilst relatively small pricing discrepancies resulted in disproportionately large variations in calculated implied volatility. Such behaviour appears consistent with the increasing influence of time decay and the reduced amount of remaining time value as maturity approaches.

Consequently, short maturities should be interpreted with caution when constructing continuous volatility surfaces. A markedly different picture emerged for longer maturities.

Options with approximately six months to one year remaining until expiration generally produced substantially smoother implied volatility structures. The resulting volatility skews exhibited the regular downward slope commonly described in the empirical literature, with only limited local distortions across neighbouring strike prices.

These observations proved particularly apparent for the CAC 40 index.

The high liquidity of the option market appeared to facilitate a more stable estimation of implied volatility, thereby providing a significantly more coherent representation of the underlying volatility surface. An equally important observation concerns the distinction between index options and individual equity options.

Whilst the CAC 40 generated relatively stable and interpretable volatility structures, several individual equities produced substantially noisier results. In certain cases, isolated market quotations generated implausibly high or even negative implied volatility estimates, suggesting either temporary pricing inconsistencies or insufficient market liquidity.

Such observations reinforce an important practical consideration.

The proposed methodology appears particularly well suited to highly liquid option markets where quoted premiums reflect continuous interaction between buyers and sellers. Conversely, less liquid markets may introduce local pricing distortions capable of obscuring the global characteristics of the volatility surface.

These preliminary observations do not constitute definitive statistical conclusions.

Nevertheless, they suggest that both liquidity and maturity represent essential prerequisites when analysing implied volatility surfaces for decision-support purposes.

Implied Volatility Curves

Figure 7. Comparison of implied volatility curves obtained for different maturities. Short-dated maturities frequently exhibit irregular local behaviour owing to limited time value and increased pricing sensitivity. Longer maturities generally produce smoother volatility skews, thereby facilitating the interpretation of surface dynamics.

A further observation emerged during the analysis: although several volatility surfaces displayed the expected downward skew, not all of them generated identical decision-support signals.

Certain maturities exhibited a progressive flattening of the skew whilst implied volatility remained at comparatively elevated levels. Others retained a persistent steep slope despite similar average volatility levels. This distinction proved particularly informative. If confirmed through broader empirical investigation, it suggests that the overall geometry of the volatility surface may contain additional information beyond the absolute level of implied volatility alone. From the perspective of systematic option selling, this observation may prove significant.

A market characterised by elevated implied volatility, and a progressively normalising volatility surface appears fundamentally different from one in which both implied volatility and downside protection demand continue to increase simultaneously.

The former may correspond to a market gradually returning towards equilibrium. The latter may still reflect an environment dominated by uncertainty.

Consequently, analysing the dynamics of the volatility surface rather than its static characteristics alone may provide a richer description of prevailing market conditions.

The following section illustrates how these observations may be translated into a practical decision-support framework for systematic cash-secured put strategies.

Discussion

The preliminary observations presented above suggest that the practical usefulness of the implied volatility surface depends upon two essential conditions: the quality of market data and the maturity of the option contracts under consideration.

The first point appears relatively intuitive.

Implied volatility is not directly observable. It is inferred from quoted option prices through the inverse application of the Black-Scholes-Merton model. Consequently, any inconsistency in market quotations is immediately reflected in the calculated implied volatilities.

This phenomenon proved particularly evident during the exploratory analyses conducted on individual equities.

Whilst certain option chains generated coherent volatility structures, others produced isolated implied volatility values that were incompatible with neighbouring strike prices. In a limited number of cases, implausible or unstable implied volatility estimates were obtained despite apparently valid market quotations. Such behaviour most likely reflects temporary liquidity deficiencies, unusually wide bid-ask spreads or isolated transactions executed outside normal market conditions.

These observations underline an important methodological requirement.

The proposed framework should preferably be applied to option markets characterised by sufficient liquidity and a broad distribution of actively traded strike prices. Under such conditions, quoted premiums are more likely to represent the consensus valuation of market participants rather than isolated transactions.

The second observation concerns option maturity.

Short-dated contracts frequently produced irregular volatility profiles whose local fluctuations appeared dominated by pricing noise rather than genuine changes in market expectations. As expiration approaches, the remaining time value becomes progressively smaller, and option prices exhibit increasing sensitivity to relatively minor changes in the underlying asset. Consequently, the resulting implied volatility estimates become substantially less stable.

Conversely, medium- and long-dated maturities generally generated considerably smoother volatility structures.

The downward skew remained clearly identifiable whilst local distortions became significantly less pronounced. This regularity considerably facilitated the interpretation of the surface and its evolution over successive market observations.

Among the datasets examined, listed CAC 40 index options consistently provided the most coherent results. Their combination of high liquidity, narrow bid-ask spreads and broad strike availability produced volatility surfaces whose overall geometry remained remarkably stable. This characteristic makes such instruments particularly well suited to exploratory research concerning the dynamics of implied volatility.

An additional observation deserves particular attention: not every regular volatility surface generated the same analytical conclusion.

Certain maturities displayed a progressive flattening of the volatility skew whilst implied volatility remained comparatively elevated. Others retained a persistent and pronounced downward slope despite exhibiting similar average volatility levels. This distinction appears especially interesting.

If future empirical analyses confirm these preliminary observations, the evolution of the volatility surface may provide information that cannot be obtained from the absolute level of implied volatility alone. Such a conclusion would carry practical implications for systematic option-selling strategies.

Rather than selecting opportunities exclusively according to premium levels or historical volatility, investors may benefit from incorporating the dynamics of the implied volatility surface into their broader decision-making process. Naturally, these findings should be interpreted with appropriate caution.

The present work remains exploratory in nature and does not claim to establish a predictive model. Instead, it proposes an analytical framework intended to organise market information already embedded within option prices into a more coherent decision-support process.

Further empirical investigation involving longer observation periods, multiple market regimes and additional underlying assets will naturally be required before more general conclusions may be drawn.

Evolution of the Implied Volatility Surface

Figure 8. Evolution of the implied volatility surface across successive market observations. The figure illustrates the conceptual distinction between a market in which the volatility skew continues to steepen and one in which the surface progressively normalises whilst implied volatility remains comparatively elevated.

Practical Implications for Systematic Put Sellers

From a practical perspective, the observations discussed throughout this article suggest that implied volatility should perhaps be interpreted less as an isolated numerical indicator and more as one component of a broader analytical framework.

Option sellers have traditionally focused on premium maximisation. Although this objective remains entirely legitimate, premium alone provides only a partial description of prevailing market conditions.

The same premium may arise under markedly different market environments. One may correspond to an increasingly stressed market characterised by rapidly rising demand for downside protection. Another may reflect a market in which uncertainty remains elevated but has already begun to stabilise.

Distinguishing between these situations may prove particularly valuable when implementing systematic cash-secured put strategies. Rather than attempting to forecast market direction, the proposed framework encourages a more disciplined interpretation of the information continuously embedded within option prices.

In this respect, the implied volatility surface becomes considerably more than a graphical representation of option quotations. It evolves into a dynamic indicator describing the collective perception of risk within financial markets.

Conclusion

The Black-Scholes-Merton model remains the fundamental reference upon which modern option pricing is built. Although one of its central assumptions — constant volatility — is systematically contradicted by market observations, these apparent discrepancies have progressively become one of the richest sources of information available to option practitioners.

The implied volatility surface should therefore not merely be regarded as a technical consequence of option pricing theory. It reflects the collective judgement of market participants regarding future uncertainty, the asymmetrical pricing of downside risk and the continuously evolving balance between buyers and sellers of financial protection. The purpose of the present study has been to explore whether this information may be exploited beyond its traditional pricing function.

Rather than concentrating exclusively on the absolute level of implied volatility, this article has proposed a broader analytical perspective based upon the dynamics of the entire volatility surface. Attention has been devoted to the progressive evolution of the volatility skew, whose gradual normalisation may provide additional insight into changing market conditions.

The preliminary empirical observations presented throughout this paper suggest practical conclusions.

First, market liquidity appears to constitute a fundamental prerequisite for obtaining sufficiently stable implied volatility surfaces. Highly liquid option markets, such as listed CAC 40 index options, produce considerably more coherent structures than many individual equity options, whose implied volatilities may occasionally be distorted by isolated transactions or limited trading activity.

Secondly, option maturity also plays a decisive role. Medium- and long-dated contracts generally generate smoother volatility surfaces that appear more suitable for structural analysis than very short-dated maturities, where the increasing influence of time decay frequently introduces substantial local irregularities.

Finally, and perhaps most importantly, the observations suggest that two markets exhibiting comparable average implied volatility levels may nevertheless convey markedly different information through the geometry of their respective volatility surfaces. This distinction may prove particularly relevant for systematic cash-secured put strategies.

Whilst elevated implied volatility undoubtedly increases option premiums, the progressive normalisation of the volatility surface may provide complementary information regarding the evolution of collective market expectations. The proposed framework should therefore not be interpreted as a predictive model. Financial markets remain inherently uncertain, and no analytical methodology can eliminate investment risk.

Instead, the approach presented here seeks to organise information already embedded within option prices into a structured decision-support framework capable of complementing more traditional valuation techniques. Viewed from this perspective, the implied volatility surface ceases to be merely a graphical representation of option prices. It becomes a dynamic description of market behaviour.

Understanding how this structure evolves through time may ultimately prove as informative as measuring its absolute level at any single observation date.

Limitations and Future Research

The present study should be regarded as an exploratory investigation rather than a definitive empirical validation. Several limitations naturally remain.

The observations reported here are based upon a limited number of underlying assets and observation dates. Broader empirical investigations covering multiple market regimes, longer historical periods and additional asset classes will be required before more general conclusions may be established. Future research could also investigate whether quantitative indicators describing the geometry of the implied volatility surface — such as skew slope, local curvature or cross-sectional dispersion — may be systematically incorporated into algorithmic decision-support models for option-selling strategies.

Another promising avenue concerns the comparative behaviour of implied volatility surfaces across different asset classes, including equity indices, individual equities, exchange-traded funds and commodity options.

Finally, machine learning techniques may eventually provide complementary tools capable of identifying recurring patterns within the evolution of volatility surfaces. Such approaches, however, should be viewed as extensions of the present analytical framework rather than substitutes for the economic interpretation of market behaviour. Ultimately, the principal contribution of this work lies less in proposing a new pricing model than in suggesting an alternative way of interpreting information already contained within option markets. If the geometry of the implied volatility surface indeed reflects the collective perception of financial risk, then monitoring its evolution may offer valuable additional insight into the timing of systematic option-selling strategies.

Download the Summary Infographic

Readers wishing to retain a concise visual summary of the concepts presented throughout this article may download the accompanying high-resolution infographic below.

Download the Summary Infographic (High-Resolution PDF)

Related posts on the SimTrade blog

   ▶ Jayati WALIA Brownian Motion in Finance

   ▶ Jayati WALIA Black-Scholes-Merton option pricing model

   ▶ Saral BINDAL Implied Volatility and Option Prices

   ▶ Saral BINDAL Volatility curves: smiles and smirks

   ▶ Saral BINDAL Implied Volatility Surface: Smiles, Smirks and Term Structure

Useful resources

Black, F., & Scholes, M. (1973). The Pricing of Options and Corporate Liabilities. Journal of Political Economy, 81(3), 637-654.

Gatheral, J. (2006). The Volatility Surface: A Practitioner’s Guide. John Wiley & Sons.

Hull, J. C. (2024). Options, Futures and Other Derivatives (11th ed.). Pearson.

Merton, R. C. (1973). Theory of Rational Option Pricing. The Bell Journal of Economics and Management Science, 4(1), 141-183.

Natenberg, S. (2015). Option Volatility and Pricing (2nd ed.). McGraw-Hill Education.

Rebonato, R. (2004). Volatility and Correlation: The Perfect Hedger and the Fox. John Wiley & Sons.

Taleb, N. N. (1997). Dynamic Hedging: Managing Vanilla and Exotic Options. John Wiley & Sons.

About the Author

Tis article was written in July 2026 by Frédéric VALOGNES , who is a lecturer in corporate finance, financial analysis, financial markets and derivatives, with more than twenty-five years of professional experience spanning financial management, higher education, research administration and executive training. He is a Certified European Financial Analyst (CEFA®), a professional designation awarded by the European Federation of Financial Analysts Societies (EFFAS), Frankfurt.

Author’s Note

This article is intended solely for educational and research purposes. It presents the author’s personal reflections on implied volatility, option pricing and systematic option-selling strategies. It should not be construed as investment advice or as a recommendation regarding any financial instrument or trading strategy.

The ideas developed in this article are the result of many years of teaching, professional practice and ongoing research in corporate finance, financial analysis, financial markets and derivatives. They have also been enriched by numerous discussions with academics, finance professionals and market practitioners, whose expertise, critical insights and constructive exchanges have played an important role in shaping the analytical framework presented here.

The author wishes to express his sincere gratitude to all those who have contributed, directly or indirectly, to the development of these ideas. Their encouragement, intellectual generosity and commitment to rigorous financial analysis have been a constant source of inspiration.

Option Implied Risk-Neutral Distribution

Saral BINDAL

In this article, Saral BINDAL (Indian Institute of Technology Kharagpur, Metallurgical and Materials Engineering, 2024-2028 & Research assistant at ESSEC Business School) explains how option prices can be used to build an implied risk-neutral distribution.

Introduction

Derivative markets provide a rich source of information for market expectations. For example, a futures price is the market’s expectation of the future value of an asset. More interestingly, we can derive the moments of the statistical distribution of future asset values from the market prices of options, like the variance (second moment), the skewness (third moment) and the kurtosis (fourth moment). More generally, we can extract the ex-ante risk-neutral probability distribution of future asset prices at a given date from option market prices with the corresponding maturity date.

Physical vs Risk-Neutral Probability Measures

A real-world probability measure represents the statistical distribution of asset returns typically estimated using historical data. These measures incorporate risk premia, market frictions, and investor behaviour, and are primarily used for statistical inference and risk modelling.

In contrast, risk-neutral probability measure is a mathematical pricing measure used in no-arbitrage valuation of financial derivatives. Under this framework, asset prices are evaluated as discounted expected payoffs under an equivalent martingale measure. In this setting, the expected return of any risky asset is adjusted to the risk-free rate within the pricing measure, simplifying valuation by transforming uncertain future payoffs into present values computed via expectation (Hull, 2018; Shreve, 2004).

Historical vs Risk-Neutral Distributions

Historical Distributions are constructed from observed past returns under the physical measure (P-measure). They empirically capture the true statistical behaviour of asset prices, including fat tails, skewness, and volatility clustering driven by real market shocks and investor behaviour. These distributions exhibit higher variance and kurtosis, making them particularly valuable for stress testing, Value-at-Risk estimation, and portfolio risk management where realistic loss scenarios matter.

Risk-Neutral Distributions are derived from option market prices rather than historical data, under the implied measure by no-arbitrage pricing (Q-measure). They reflect market-implied expectations of future payoffs discounted at the risk-free rate resulting in smoother, less skewed densities. While highly effective for pricing derivatives and contingent claims, they tend to underestimate tail risk and do not directly represent the actual probabilities investors assign to future market outcomes.

Risk-neutral distribution: the Black–Scholes–Merton framework

Having distinguished between the physical and risk-neutral probability measures, it is useful to examine the risk-neutral distribution implied by the Black–Scholes–Merton (BSM) model, which is a standard model in quantitative finance. The BSM framework assumes that the underlying asset follows a geometric Brownian motion and provides a simple illustration of how the transition from the physical measure to the risk-neutral measure alters the distribution of future asset prices.

Under the BSM, the standard assumption is that the underlying asset follows a geometric Brownian motion given by the following expressions:


SDE for the geometric Brownian motion (GBM)

where:

  • St = asset price at time t t
  • μ = drift (growth rate of the asset price)
  • r = risk-free rate
  • σ = volatility (standard deviation)
  • dWt/dWtQ = infinitesimal increment of wiener process (N(0,dt)) under respective measures

Solving these stochastic differential equations over the interval [0, T] yields the terminal asset price:


Terminal asset price formulas

Taking logarithms shows that the terminal log-price is normally distributed:


Distributions under the BSM framework

Thus, under the Black–Scholes–Merton framework, the risk-neutral distribution of the terminal asset price is lognormal (as the physical distribution). Relative to the corresponding physical distribution, the volatility remains unchanged, while the drift parameter μ is replaced by the risk-free rate r. This is an important result as the risk-free rate r is known and easily observable while the drift parameter μ has to be estimated and is not directly observable.

Butterfly spread

To extract a continuous risk-neutral probability distribution from the market, we must first understand how to isolate the market’s view on a specific future asset price. The primary tool for this is a classic option trading strategy: the butterfly spread.

A butterfly spread is an options trading strategy designed to achieve limited profit with strictly bounded risk, typically in market environments where relatively small price movements are anticipated. The strategy may be implemented using either call or put options and can be established in either a long or short configuration. For example, a long call butterfly is constructed by purchasing one call option at a lower strike price, selling two call options at an intermediate strike price, and purchasing one call option at a higher strike price. Depending on the relative spacing between the strike prices, a butterfly spread may be either symmetric or asymmetric.

Cost of a Symmetric Butterfly Spread

To understand how option market prices encode the market’s expectations regarding the future distribution of the underlying asset price, we consider a symmetric butterfly. A symmetric butterfly spread is constructed using three European call options with a common maturity T and distinct strike prices. The strategy involves purchasing one call option with strike K – ΔK at a premium of C(K-ΔK,T), selling two call options with strike K at a premium of C(K,T) each, and purchasing one call option with strike K + ΔK at a premium of C(K+ΔK,T).

The price of the resulting butterfly spread is therefore given by


Butterfly spread cost

The net cost of the butterfly spread is obtained by summing the premia paid for the two long call positions and subtracting the premiums received from the two short call positions.

Payoff of a Symmetric Butterfly Spread

The payoff of a symmetric butterfly spread is centred around the strike (K) and can be expressed as


Butterfly spread payoff

Figure 1 illustrates the payoff profile of a symmetric butterfly spread centred at the strike K = 100 with strike spacing ΔK = 5. The payoff reaches its maximum when the terminal asset price ST equals the strike K and declines to zero as ST moves beyond the adjacent strikes K – ΔK and K + ΔK.

Figure 1. Symmetric Butterfly Spread Payoff at Maturity
Symmetric Butterfly Spread Payoff  at Maturity
Source: computation by the author.

As a result, the butterfly spread effectively isolates a narrow range of terminal asset prices, making it a useful instrument for extracting information about the market-implied probability distribution of the underlying asset price at maturity.

Stacked Butterfly Spreads

A stack of butterfly spreads refers to a collection of butterfly spreads constructed across a range of strike prices, such that the central strike of each butterfly is equally spaced from the next. The spacing between successive central strikes is equal to the strike spacing ΔK used in the construction of each individual butterfly spread, as discussed above.

Figure 2 illustrates that a collection of butterfly spreads across strikes at a fixed maturity converges to the market-implied probability density of the underlying asset. Each butterfly corresponds to a discrete approximation of the second derivative of option prices with respect to strike, and aggregating these across strikes recovers the risk-neutral density.

We construct seven butterfly spreads centered at strikes K = 85 to K = 115 in increments of 5, with strike spacing ΔK = 5. The weights are specified using a Gaussian distribution with mean μ = 100 and standard deviation σ = 10, reflecting an assumed market belief about the concentration of terminal prices. The payoff profile is scaled by a factor of 200 to improve visual readability, and it is normalized by ΔK2 to remain consistent with the second-order finite-difference interpretation of butterfly spreads as detailed below.

Figure 2. Approximating the Risk-Neutral Density Using Butterfly Spreads
Approximating the Risk-Neutral Density Using Butterfly Spreads
Source: computation by the author.

As the strike spacing ΔK is reduced, additional butterfly spreads can be constructed between existing butterfly spreads. Consequently, the stacked payoff profile becomes increasingly smooth and, in the limit, approaches a continuous representation of the implied probability distribution.

To better understand this limiting behaviour, it is useful to examine the properties of an individual butterfly spread. As the strike spacing ΔK decreases, the payoff of the butterfly spread becomes increasingly concentrated around its central strike. In the limit as ΔK → 0, the butterfly spread approaches an infinitesimally narrow peak centred at K.

Consequently, the value of the butterfly spread decreases as its payoff becomes increasingly concentrated around its central strike. To obtain a meaningful limiting quantity, the butterfly value must therefore be normalized by (ΔK)2. This normalization is motivated by a well-known result from calculus, central finite-difference approximation of the second derivative.


Normalized Butterfly spread cost

Comparing the two expressions above, reveals that the normalized butterfly value is precisely the finite-difference approximation of the second derivative of the call pricing function with respect to strike.


Second derivative of the call pricing function with respect to strike.

This observation forms the foundation of the Breeden-Litzenberger (1978) result, which establishes that the second derivative of the call pricing function with respect to strike is directly related to the market-implied risk-neutral probability density embedded in option prices, as demonstrated in the derivation below.

You can download the Excel file provided below to generate and visualize the payoff profiles of the butterfly spread and stacked butterfly spread at maturity, as discussed above.

Download the Excel file.

Option implied risk-neutral distribution

This section develops the analytical derivation of the risk-neutral distribution using the seminal Breeden-Litzenberger (1978) result. By exploiting the cross-sectional structure of option prices across strikes, we recover the market-implied risk-neutral density embedded in option market prices.

Analytical derivation

Under the risk-neutral measure, the value of a European call option is given by the present value of its expected payoff at maturity. For a strike price K, continuously compounded risk-free rate r, and time to maturity T, the call pricing function C(K,T) can be expressed as


Call option risk-neutral value.

To obtain a continuous representation of the call price, the expected payoff can be expressed as an integral over the probability density function of the terminal asset price, f(ST).


Call option risk-neutral value PDF.

Note: The integral starts at K because the payoff is zero when St≤K.

Taking the first derivative with respect to K, we get


Call option risk-neutral PDF first derivative

To obtain the risk-neutral probability density function, as shown by Breeden and Litzenberger (1978), we take an additional derivative with respect to the strike


Second derivative of call price with respect to strike.

Rearranging the above formula, we get the risk-neutral distribution


Rearranged Second derivative of call price with respect to strike.

Applying the second-order central difference approximation heuristically developed in the previous section using butterfly spreads, we obtain the following expression:


Implied risk-neutral distribution formula.

This expression shows that the risk-neutral probability density can be recovered directly from the second derivative of the call pricing function with respect to strike. In practice, however, option prices are observed only at a finite set of discrete strike prices, requiring numerical methods to approximate the derivatives and extract the implied risk-neutral distribution.

Numerical methods for extracting the risk-neutral distribution

Methods for extracting the risk-neutral distribution can be broadly classified into non-parametric (data-driven with minimal distributional assumptions), semi-parametric (partial structural assumptions, typically imposed on intermediate quantities such as implied volatility), and parametric or structural (explicit assumptions on the distribution or asset price dynamics) approaches. These methodologies differ in the degree of modelling assumptions imposed on the option pricing function and the terminal asset price distribution, leading to different trade-offs between flexibility, numerical stability, and economic interpretability.

Non-parametric methods

Non-parametric methods aim to recover the risk-neutral distribution directly from observed option prices without imposing any specific parametric structure on either the terminal asset price distribution or the stochastic process governing the evolution of the underlying asset price. Consequently, these methods are highly flexible, but they tend to be sensitive to market microstructure noise, sparse strike coverage, and interpolation error in option quotes.

Risk-neutral histograms: the most direct implementation of the Breeden–Litzenberger result constructs a discrete approximation of the implied risk-neutral density using finite differences across traded strikes (Breeden and Litzenberger, 1978; Neuhaus, 1995). Adjacent butterfly spreads may therefore be interpreted as local estimates of state-contingent probabilities.

Because option contracts are quoted only at discrete strike intervals, the recovered distribution resembles a histogram rather than a smooth continuous density, making the approach highly sensitive to strike spacing and pricing noise.

Kernel regression methods: to mitigate the instability of histogram-based estimates, subsequent research introduced non-parametric smoothing techniques that estimate a continuous option pricing function directly from observed market prices. A prominent example is the kernel regression framework of Aït-Sahalia and Lo (1998).

By reducing the influence of local pricing noise, kernel-based methods generally produce smoother and more stable estimates of the implied risk-neutral density.

Spline-based methods: another widely used class of non-parametric methods employs spline interpolation techniques to construct smooth and arbitrage-consistent call pricing functions across strikes (Bates, 1991). Once a sufficiently smooth pricing function has been obtained, the implied risk-neutral density can be recovered through numerical differentiation.

Spline-based approaches offer substantial flexibility but remain sensitive to data quality and sparse observations in the tails of the distribution.

Semi-parametric approaches

Semi-parametric approaches occupy a middle ground between purely data-driven and fully parametric methodologies. Rather than modelling the risk-neutral density directly, these methods impose structure on intermediate quantities, most commonly the implied volatility smile.

Implied volatility smile methods: in practice, many market participants smooth the implied volatility smile rather than the option prices directly. Observed option prices are first converted into implied volatilities, after which a smooth volatility smile is fitted across strikes using parametric specifications or spline-based interpolation techniques (Shimko, 1993).

The smoothed volatility smile is subsequently mapped back into option prices, allowing the implied risk-neutral density to be recovered through numerical differentiation. These methods generally exhibit greater numerical stability, although tail estimation remains sensitive to extrapolation assumptions in illiquid regions of the smile.

Parametric and structural approaches

Parametric and structural methodologies recover the implied risk-neutral distribution by imposing explicit assumptions on either the terminal distribution of asset prices or the stochastic process governing their evolution.

Parametric density models: a prominent class of methods assumes that the terminal risk-neutral distribution follows a particular parametric specification. One widely used approach models the distribution as a mixture of lognormal densities calibrated to observed option prices (Bahra, 1997; Melick and Thomas, 1997).

Parametric methods are computationally efficient and often yield economically interpretable measures of skewness, kurtosis, and tail risk. Their flexibility, however, is inherently constrained by the assumed functional form.

Dynamic option pricing models: rather than specifying the terminal distribution directly, structural approaches derive the implied density from an assumed stochastic process governing the evolution of the underlying asset price. Examples include stochastic volatility and jump-diffusion frameworks calibrated to observed option prices (Bates, 1995; Malz, 1995).

Within these models, the risk-neutral density emerges endogenously from the dynamics of the underlying asset under the risk-neutral measure. While theoretically appealing, such models are computationally intensive and sensitive to model misspecification.

Application

Implementing the Breeden and Litzenberger (1978) result in practice requires a continuum of European option prices written on the same underlying asset, all sharing a common maturity and spanning a continuous range of strike prices from zero to infinity. Under such idealized conditions, the risk-neutral density can be recovered directly from the cross-section of option prices (at a given maturity date).

In practice, however, listed option markets provide only a sparse and discrete grid of strike prices, typically concentrated around the at-the-money (ATM) region. The absence of a complete continuum of option strikes, particularly in the deep in-the-money and far out-of-the-money regions, necessitates the use of interpolation across observed strikes and extrapolation into the tails in order to recover a smooth and arbitrage-free implied risk-neutral distribution.

Required data

Constructing a risk-neutral distribution requires option chain data (a set of calls and/or puts) for a single maturity, along with the underlying asset price, the prevailing risk-free rate, dividend assumptions, at the exact observation time of the market data.

Such data can be obtained from both free and commercial data providers. One of the most accessible sources is Yahoo! Finance; however, freely available option data is often subject to inconsistencies such as wide bid–ask spreads, stale quotes, and incomplete cross-sectional coverage of strikes, all of which can materially distort empirical estimation of the risk-neutral distribution (RND).

For our application, we employ simulated option data to illustrate the derivation of the implied risk-neutral distribution from an option chain within a controlled and internally consistent setting. This ensures that the resulting distribution remains aligned with the theoretical framework developed above.

Extraction of the implied risk-neutral density

From the collected option chain data, we first apply a series of standard filtering procedures designed to remove illiquid and economically inconsistent observations. In empirical applications, this typically includes liquidity screens, moneyness and maturity filters, implied-volatility sanity checks, and no-arbitrage constraints to mitigate errors arising from stale quotes, asynchronous observations, and market microstructure noise. Since the dataset employed here is simulated and internally consistent by construction, these preprocessing steps can be largely omitted.

Figure 3 below presents the implied volatility smile obtained from the simulated European call option chain after numerical inversion of the Black–Scholes–Merton pricing model. The smile is interpolated using a natural cubic spline over a dense strike grid spanning the filtered strike range of 4,000 to 6,000, under the assumptions of an underlying spot price of $5,300, a continuously compounded risk-free interest rate of 5.2%, and a remaining time-to-maturity of 30 days. The resulting smooth volatility curve serves as the key intermediate input for constructing a continuous and differentiable call pricing function required for subsequent risk-neutral density extraction.

Figure 3. Implied Volatility Smile
Implied Volatility Smile
Source: computation by the author (with python)

The interpolated implied volatility smile is subsequently utilized to reprice European call options across a finely discretized strike grid, thereby constructing a smooth numerical approximation of the cross-sectional call price surface. The option implied risk neutral density is then recovered by applying the Breeden Litzenberger operator, corresponding to the second partial derivative of discounted call prices with respect to strike, to the smoothed pricing function. Figure 4 illustrates the resulting risk neutral density extracted from the simulated European call option chain under an underlying spot level of $5,300, a continuously compounded risk-free interest rate of 5.2%, and a remaining time to maturity of 30 days.

Figure 4. Implied Risk-Neutral Distribution
Implied Risk-Neutral Distribution
Source: computation by the author (with python)

You can download the Python code provided below for generating simulated call option chain data and the option-implied risk-neutral distribution, as discussed above.

Download the Python code.

Alternatively, you can download the R code below with the same functionality as in the Python file.

 Download the R code.

Empirical issues

A primary limitation in empirical recovery of the risk-neutral distribution is the discrete nature of listed option strikes. The Breeden–Litzenberger framework assumes a continuum over strike space, whereas traded options are observed only on a sparse and uneven grid concentrated around the at-the-money region.

A second limitation arises from the unobservability of the distribution tails. Deep in-the-money and far out-of-the-money options are often illiquid or not quoted, implying that tail behaviour of the risk-neutral density must be inferred through extrapolation rather than direct market observation.

A separate issue is asynchronous option quotes. Since option prices across strikes are not necessarily recorded simultaneously, the resulting cross-section may embed timing mismatches, introducing bias in the reconstructed pricing function. This is typically addressed using end-of-day settlement data or synchronized snapshots.

In addition, different levels of market liquidity (due to different levels of bid ask spreads for example) across strikes introduces noise and heterogeneity in observed quotes. Illiquid contracts may exhibit stale or unreliable prices, which can distort the implied volatility surface even after basic filtering.

Finally, the reconstruction procedure does not explicitly impose no-arbitrage conditions or global smoothness constraints across strikes. As a result, when option prices are interpolated to form a continuous surface, the fitted call price function may exhibit local violations of convexity in strike space (e.g., small regions where butterfly spreads imply negative prices or non-monotonic curvature). Such violations are problematic because they imply the possibility of arbitrage and can lead to risk-neutral probability estimates that are not economically consistent.

Despite these limitations, the framework remains a useful reduced-form tool for extracting risk-neutral densities, provided appropriate smoothing and arbitrage constraints are imposed.

Real-life applications

Central Bank Monetary Policy Monitoring

Bahra (1997) and Kim (2009) suggest that policymakers extract ex-ante risk-neutral distributions (RNDs) from interest rate, equity, and currency options to assess market-implied expectations and uncertainty around policy decisions. Unlike futures prices, which only reflect the conditional mean, RNDs incorporate higher-order information such as skewness and kurtosis, allowing for a more complete assessment of perceived tail risks and macro-financial stress. For example, during the February 2007 equity sell-off, the European Central Bank (ECB, 2007) used option-implied probability distributions (“fan charts”) to assess whether the move reflected extreme tail risk and to track the evolution of market expectations after stabilization.

Value-at-Risk (VaR) Forecasting

Risk management units in investment banks use quantiles derived from implied RNDs to forecast extreme portfolio losses in a forward-looking manner. Compared to traditional historical simulation methods, RND-based approaches incorporate market-implied expectations and have been shown to provide improved performance relative to standard volatility-based models such as GARCH(1,1) (Chang, Chang, Huang, & Hsieh, 2011).

Systemic Risk and Stress Testing Indicator

Macroprudential regulators transform option-implied volatility surfaces into arbitrage-consistent risk-neutral distributions to quantify system-wide financial vulnerabilities. By aggregating tail-risk measures across equities, currencies, and interest rates, these distributions can be used to construct time-series indicators of systemic stress and cross-asset fragility (Malz, 2014).

Market Risk Aversion and Investor Sentiment Estimation

By combining option-implied risk-neutral distributions with empirical (physical) distributions, researchers can infer the market’s implicit risk preferences and aggregate degree of risk aversion (Bliss & Panigirtzoglou, 2004). This allows for the identification of time variation in investor sentiment and risk pricing across different investment horizons (Bliss & Panigirtzoglou, 2004; Gemmill & Saflekos, 2000).

Why should you be interested in this post?

The risk-neutral distribution is one of the few tools in finance that reveals how the market prices uncertainty based on the entire distribution of possible future states implied by option prices. It is widely used in practice to understand how the market is pricing downside risk, fat tails, and asymmetry that is directly used in volatility modelling, pricing, and risk management frameworks. From a practical perspective, it is one of the standard tools used to extract forward-looking information from option prices in both research and industry settings.

Related posts on the SimTrade blog

   ▶ Saral BINDAL Historical Volatility

   ▶ Saral BINDAL Implied Volatility and Option Prices

   ▶ Saral BINDAL Volatility curves: smiles and smirks

Useful resources

Academic research on option pricing

Black, F., & Scholes, M. (1973). The pricing of options and corporate liabilities. Journal of Political Economy, 81(3), 637-654.

Hull J.C. (2015) Options, Futures, and Other Derivatives, Eighth Edition, Global Edition, Chapter 14 – The Black-Scholes-Merton model, 299-320.

Merton, R.C. (1973). Theory of rational option pricing. The Bell Journal of Economics and Management Science, 4(1), 141-183.

Academic research on risk neutral distribution

Aït-Sahalia, Y., & Lo, A. W. (1998). Nonparametric estimation of state-price densities implicit in financial asset prices. The Journal of Finance, 53(2), 499-547.

Bahra, B. (1997). Implied risk-neutral probability density functions from option prices: Theory and application. Bank of England Working Paper Series, 66, 1-42.

Bates, D. S. (1991). The crash of ’87: Was it expected? The evidence from options markets. The Journal of Finance, 46(3), 1009-1044.

Bates, D. S. (1995). Testing option pricing models. NBER Working Paper Series, w5135, 1-53.

Bliss, R. R., & Panigirtzoglou, N. (2004). Option-implied risk aversion estimates. The Journal of Finance, 59(1), 407-446.

Breeden, D. T., & Litzenberger, R. H. (1978). Prices of state-contingent claims implicit in option prices. Journal of Business, 51(4), 621-651.

Chang, Y. C., Chang, C. L., Huang, H. T., & Hsieh, T. H. (2011). Value-at-Risk forecasting via option-implied risk-neutral density. Journal of Risk and Financial Management, 4(1), 56-83.

European Central Bank (ECB). (2007). Gauging stock market uncertainty using option-implied distributions. ECB Monthly Bulletin, April, Box 4, 31–32.

Figlewski, S. (2010). Estimating the implied risk neutral density for the U.S. market portfolio. In T. Bollerslev, J. R. Russell, & M. W. Watson (Eds.), Volatility and Time Series Econometrics: Essays in Honor of Robert F. Engle (pp. 43-69). Oxford University Press.

Gemmill, G., & Saflekos, A. (2000). How useful are market-implied probabilities for forecasting sharp changes in asset prices? An application to the UK general election. Market Expectations and the Implications for Monetary Policy, 203-223.

Kim, K. (2009). Monetary policy announcements and market expectations under different monetary policy regimes: An options-based approach. International Finance Discussion Papers (Federal Reserve Board), 977, 1-45.

Malz, A. M. (1996). Using option prices to estimate realignment probabilities in the European Monetary System: the case of sterling-mark. Journal of International Money and Finance, 15(5), 717-748.

Malz, A. M. (2014). A VaR-based systemic risk indicator. Federal Reserve Bank of New York Staff Reports, 668, 1-47.

Melick, W. R., & Thomas, C. P. (1997). Recovering an asset’s pdf from option prices: An application to crude oil during the Gulf crisis. Journal of Financial and Quantitative Analysis, 32(1), 91-115.

Neuhaus, H. (1995). The informational content of derivatives for monetary policy. Deutsche Bundesbank Discussion Paper Series 1: Economic Studies, 1995(03), 1-34.

Shimko, D. (1993). Bounds of probability. Risk, 6(4), 33-37.

Shreve, S. E. (2004). Stochastic calculus for finance II: Continuous-time models. Springer Science & Business Media.

About the author

The article was written in June 2026 by Saral BINDAL (Indian Institute of Technology Kharagpur, Metallurgical and Materials Engineering, 2024-2028 & Research assistant at ESSEC Business School).

   ▶ Discover all articles by Saral BINDAL

Capital safety: μ versus σ

Saral BINDAL

In this article, Saral BINDAL (Indian Institute of Technology Kharagpur, Metallurgical and Materials Engineering, 2024-2028 & Research assistant at ESSEC Business School) explains the expression “Mu versus Sigma” often used in asset management to study the impact of expected performance and risk on the safety of capital invested by individuals.

Introduction

In financial markets, performance over time is governed by three fundamental variables: the drift (μ), volatility (σ), and –perhaps most importantly– time (T). The drift represents the expected growth rate of the price and corresponds to the expected return of assets or portfolios. Volatility measures the uncertainty or risk associated with price fluctuations around this expected growth and corresponds to the standard deviation of returns. Time is not merely the mechanical passage of calendar time; it also embodies the investment horizon chosen by the investor, that is, the period over which risk is borne and returns are realized. The relationship between all variables reflects the trade-off between risk and return. Time, which is related to the investment horizon set by the investor, determines how both performance and risk accumulate. Together, these variables form the foundation of asset pricing to model the behavior of market price over time, and in fine the performance of the investor at their investment horizon.

Modeling asset prices with geometric Brownian motion (GBM)

Financial models commonly usually employ geometric Brownian motion (GBM) because it is a strictly non-negative stochastic process, making it well suited for modeling asset prices (a stock price must be positive under the limited lability rule). GBM provides a mathematically tractable and economically intuitive framework for describing the continuous-time evolution of asset prices under uncertainty.

Under GBM, the proportional change in an asset price is driven by a deterministic drift component and a random shock associated to volatility, represented by the stochastic differential equation:


GBM formula

where:

  • St = asset price at time t
  • μ = drift (growth rate of the asset price)
  • σ = volatility (standard deviation)
  • dWt = infinitesimal increment of wiener process (N(0,dt))

This formulation leads to lognormally distributed prices and allows for closed-form expressions for expected values and confidence bounds, which are particularly useful in Monte Carlo simulations and risk analysis.

In practice, the continuous-time process is discretized (see Maruyama, 1955) to simulate price paths numerically. A more detailed discussion of the discrete GBM model and its implications is available in a separate post, which can be found here.

Figure 1 illustrates three simulated asset price paths generated under a geometric Brownian motion (GBM), sampled at monthly intervals (Δt = 1/12) over a 10-year horizon (T = 10). The expected price trajectory and the associated upper and lower bounds of the 66% confidence interval are shown alongside the simulations. The model assumes an annual drift (μ) of 8% and an annual volatility (σ) of 15%, with the initial asset price set at S0 = €100.

Figure 1. Monte Carlo simulated asset price paths under a Geometric Brownian Motion model.
Monte Carlo–simulated asset price paths for specific cases under a GBM model.
Source: computation by the author.

From Figure 1, the initial capital is set to $100, while the three simulated price paths illustrate alternative realizations of the asset’s stochastic evolution. These scenarios allow us to assess capital preservation over the investment horizon: when the terminal asset price remains above the initial capital level (either above the upper confidence bound or within the confidence band) the invested capital is preserved; conversely, when the price trajectory breaches the lower bound, the initial capital is eroded.

Time Scaling of Drift and Volatility: Mu versus Sigma

Under Brownian-based asset price models, the expected return grows linearly with time, while the dispersion around the mean increases proportionally to the square root of time. Accordingly, the μ and σ functions over time are given by:


Time scaling of mu and sigma

In the formula above, the drift term must be replaced by μ-σ2/2, but we keep the notation μ for simplicity.

The figure below illustrates the linear function μ(t) = t and the square-root function σ(t) = √t. These two functions are represented below for both μ and σ equal to 1.

Figure 2. Comparison of linear (t) and square-root (√t) time scaling.
Linear vs. Square-root growth function in time t.
Source: computation by the author.

From Figure 2 above, we observe that at shorter time horizons (t < 1), the square-root function dominates the linear function, indicating that randomness (risk) has a stronger influence early on. At t = 1, both the functions are equal. For longer horizons (t > 1), the linear growth of the drift increasingly dominates the slower √t growth of volatility.

As discussed previously, the modelled prices also depend on the parameters μ (drift) and σ (volatility). This dependence is illustrated below through the functions μGBM(t) = μt and σGBM(t) = σ√t, with μ = μ-σ2/2 = 6.875% and σ = 15%.

Figure 3. Drift vs. Volatility as a function of time t.
Drift vs. Volatility as a function of time t.
Source: computation by the author.

Capital safety

According to the safety-first principle (see Roy, 1952), the behavior of investors under uncertainty is driven by the tendency to minimize the probability of a disaster, defined as wealth or returns falling below a prespecified critical threshold d, rather than maximizing expected return. Formally, a disaster occurs when R < d, where R denotes expected returns. Typical disaster events include bankruptcy, severe capital loss, or failure to achieve a minimum acceptable outcome. The corresponding disaster probability is therefore given by:


Probability of disaster formula.

When returns are normally distributed, minimizing the disaster probability under the safety-first principle is equivalent to maximizing the Safety-First Ratio (SF Ratio), defined as:


SFRatio formula.

Where μ denotes the mean (expected return, E(R)) and σ the standard deviation (volatility), characterizing the distribution of returns of the asset under consideration.

When modelling asset prices or returns, realized outcomes may deviate from the mean, implying a strictly positive probability that the terminal price falls below any fixed reference level (for example, below $100 when 100 denotes the initial price).

Consequently, the lower bound of a given confidence interval is used to estimate the minimum return, or the disaster level (d), over a specified investment horizon. This lower bound can serve as a capital safety device, as capital safety is achieved by minimising the upper bound on the probability of disaster, or equivalently by maximising the SF Ratio.

In practice, such capital protection objectives are implemented through structured investment strategies such as structured products that combine fixed-income instruments with derivatives to ensure that the terminal payoff remains above a predefined threshold.

“Mu versus Sigma” and capital safety

Under the Geometric Brownian Motion (GBM), the theoretical expression for the lower bound of an asset’s logarithmic returns over time is given by:


Lower bound formula for log-returns under GBM.

For the initial capital to be guaranteed at a given probability level, the lower confidence bound of returns must be positive at the investment horizon. A positive lower confidence bound implies that, with the specified confidence level, the terminal asset value does not fall below the initial investment, thereby providing probabilistic capital protection.

Condition for capital safety
Lower bound inequality.

The evolution of lower bound over time depends on the following two cases:

Case 1: when μ − σ2/2 > 0


Lower bound inequality case 1 t formula.

where t* is the minimum investment horizon for which the lower bound turns positive.

Case 2: when μ − σ2/2 ≤ 0


Lower bound inequality case 2.

However, a non-positive quantity (zero or negative) cannot exceed a strictly positive one. As a result, the inequality can never be satisfied, and the critical time t^*does not exist.

Example: The S&P 500 index

Consider an investment in the S&P 500 index over a 10-year horizon, with an expected return (μ) of 8% and volatility (σ) of 15%. For a 66% confidence level, the critical investment horizon (t*) at which the lower confidence bound of terminal log-returns under a Geometric Brownian Motion (GBM) model becomes non-negative is given by:


Lower bound inequality case 1 example calculation.

The table below shows the evolution of lower bound with time of investment, for returns on the above-mentioned investment,. From the table we can see that for time horizon less than 5, the lower bound is negative (ie no guarantee of initial cap) but after 5 years, the lower bound is positive hence gurantee of cap

The table below reports the evolution of the lower bound of returns as a function of the investment horizon. For short horizons (less than 4.33 years), the lower bound is negative, indicating that at the chosen confidence level there remains a non-zero probability of ending with a value below the initial capital. After 4.33 years the lower bound becomes positive, such that the terminal return is positive with the specified confidence level.

Table 1: Lower bound as function of investment horizon
Lower bound inequality for a given mu and sigma table.

Thus, a minimum investment horizon of 4.33 years is required to guarantee the invested capital at the given probability.

Impact of the parameters on t*

Impact of μ on t*

As discussed above, the condition for t* to exist is (μ − σ2/2 > 0), therefore the value of μ at which t* exist is given by:


Mu* formula.

And the threshold value of μ for which the lower bound of log-returns first turns positive can be derived as follows:


Lower bound inequality for a given time horizon and sigma under GBM.

[F2.1]

For the same S&P 500 investment example discussed above, with volatility (σ) fixed at 15% and the time to maturity set at 10 years, the value of μ* is calculated below.


Mu* value.

The table below reports the minimum investment horizon t* as a function of the drift parameter μ. For drift values below 1.125%, t* does not exist, reflecting the fact that the lower confidence bound of returns never becomes positive. Once the drift exceeds this threshold, t* becomes finite and decreases rapidly as μ increases. In particular, for μ ≥ 5%, the critical horizon falls below 10 years, indicating that capital protection can be achieved over the investment horizons considered, provided the expected growth rate is sufficiently high.

Table 2: t* as a function of μ
Lower bound inequality for a given time horizon and sigma table.

Therefore, under the GBM model, a drift of at least 5.56% is necessary for the lower confidence bound of returns to remain non-negative at a 10-year horizon.

Impact of σ on t*

As discussed above, the condition for t* to exist is (μ − σ2/2 > 0), therefore the value of σ at which t* exist is given by:


Sigma* formula.

And the threshold value of σ for which the lower bound of log-returns first turns positive can be derived as follows:


Lower bound inequality for a given time horizon and mu under GBM.

[F3.1]

For the same S&P 500 investment example discussed above, with drift (μ) fixed at 8% and the time to maturity set at 10 years, the value of σ* is calculated below.


Sigma* value.

The table below reports the minimum investment horizon t* as a function of the volatility parameter &sigm;. For volatility values above 40%, t* does not exist, reflecting the fact that the lower confidence bound of returns never becomes positive. However, volatility values below this threshold, t* becomes finite and increases rapidly as σ increases. In particular, for σ ≥ 20%, the critical horizon exceeds 10 years, indicating that capital protection over a 10-year investment horizon is feasible only if volatility (risk) remains below 20%.

Table 3: t* as a function of σ
Lower bound inequality for a given time horizon and mu table.

Therefore, under the GBM model, capital preservation over a 10-year horizon at the specified confidence level is feasible only if volatility does not exceed approximately 20%.

You can download the Excel file provided below to generate asset prices modeled using geometric Brownian motion and to illustrate the time scaling of drift and volatility, as well as the capital safety calculations discussed above.

Download the Excel file.

Why should I be interested in this post?

Understanding how drift and volatility scale over time is central to the idea of capital safety and guarantee of capital. These concepts form the backbone of many structured products used by investors to protect their wealth, making them powerful tools for managing risk in uncertain markets.

Related posts on the SimTrade blog

   ▶ Saral BINDAL Historical Volatility

   ▶ Saral BINDAL Implied Volatility and Option Prices

   ▶ Jayati WALIA Brownian Motion in Finance

   ▶ Jayati WALIA Monte Carlo simulation method

Useful resources

Academic research

Bachelier L. (1900) Théorie de la spéculation. Annales scientifiques de l’École Normale Supérieure, 3e série, 17, 21–86.

de Haan L., Jansen D. W., Koedijk K. G., & de Vries C. G. (1994) Safety first portfolio selection, extreme value theory and heavy-tailed distributions. In: Extreme Value Theory and Applications, Springer, 471–487.

Kataoka S. (1963) A stochastic programming model. Econometrica, 31, 181–196.

Lawler G.F. (2006) Introduction to Stochastic Processes, 2nd Edition, Chapman & Hall/CRC, Chapter “Brownian Motion”, 201–224.

Maruyama G. (1955) Continuous Markov processes and stochastic equations. Rendiconti del Circolo Matematico di Palermo, 4, 48–90.

Roy A. D. (1952) Safety first and the holding of assets. Econometrica, 20, 431–449.

Samuelson P.A. (1965) Rational theory of warrant pricing. Industrial Management Review, 6(2), 13–39.

Telser L. G. (1955) Safety-first and hedging. Review of Economic Studies, 23, 1–16.

Wiener N. (1923) Differential-space. Journal of Mathematics and Physics, 2, 131–174.

About the author

The article was written in January 2026 by Saral BINDAL (Indian Institute of Technology Kharagpur, Metallurgical and Materials Engineering, 2024-2028 & Research assistant at ESSEC Business School).

Read all posts written by Saral BINDAL.

   ▶ Discover all articles written by Saral BINDAL

Implied Volatility and Option Prices

Saral BINDAL

In this article, Saral BINDAL (Indian Institute of Technology Kharagpur, Metallurgical and Materials Engineering, 2024-2028 & Research assistant at ESSEC Business School) explains how implied volatility is calculated or extracted from option prices using an option pricing model.

Introduction

In financial markets characterized by uncertainty, volatility is a fundamental factor shaping the pricing and dynamics of financial instruments. Implied volatility stands out as a key metric as a forward-looking measure that captures the market’s expectations of future price fluctuations, as reflected in current market prices of options.

The Black-Scholes-Merton model

In the early 1970s, Fischer Black and Myron Scholes jointly developed an option pricing formula, while Robert Merton, working in parallel and in close contact with them, provided an alternative and more general derivation of the same formula.

Together, their work produced what is now called the Black Scholes Merton (BSM) model, which revolutionized investing and led to the award of 1997 Nobel Prize in Economic Sciences in Memory of Alfred Nobel to Myron Scholes and Robert Merton “for a new method to determine the value of derivatives,” developed in close collaboration with the late Fischer Black.

The Black-Scholes-Merton model provides a theoretical framework for options pricing and catalyzed the growth of derivatives markets. It led to development of sophisticated trading strategies (hedging of options) that transformed risk management practices and financial markets.

The model is built on several key assumptions such as, the stock price follows a geometric Brownian motion with constant drift and volatility, no arbitrage opportunities, constant risk-free interest rate and options are European-style (options that can only be exercised at maturity).

Key Parameters

In the BSM model, there are five essential parameters to compute the theoretical value of a European-style option is calculated are:

  • Strike price (K): fixed price specified in an option contract at which the option holder can buy (for a call) or sell (for a put) the underlying asset if the option is exercised.
  • Time to expiration (T): time left until the option expires.
  • Current underlying price (S0): the market price of underlying asset (commodities, precious metals like gold, currencies, bonds, etc.).
  • Risk-free interest rate (r): the theoretical rate of return on an investment that is continuously compounded per annum.
  • Volatility (σ): standard deviation of the returns of the underlying asset.

The strike price (exercise price) and time to expiration (maturity) correspond to characteristics of the option while the current underlying asset price, the risk-free interest rate, and volatility reflect market conditions.

Option payoff

The payoff for a call option gives the value of the option at the moment it expires (T) and is given by the expression below:


Payoff formula for call option

Where CT is the call option value at expiration, ST the price of the underlying asset at expiration, and K is the strike price (exercise price) of the option.

Figure 1 below illustrates the payoff function described above for a European-style call option. The example considers a European call written on the S&P 500 index, with a strike price of $5,000 and a time to maturity of 30 days.

Figure 1. Payoff value as a function of the underlying asset price.
Payoff function
Source: computation by the author.

Call option value

While the value of an option is known at maturity (being determined by its payoff function), its value at any earlier time prior to maturity, and in particular at issuance, is not directly observable. Consequently, a valuation model is required to determine the option’s price at those earlier dates.

The Black–Scholes–Merton model is formulated as a stochastic partial differential equation and the solution to the partial differential equation (PDE) gives the BSM formula for the value of the option.

For a European-style call option, the call option value at issuance is given by the following formula:


Formula for the call option value according to the BSM model

with


Formula for the call option value according to the BSM model

Where the notations are as follows:

  • C0= Call option value at issuance (time 0) based on the Black-Scholes-Merton model
  • K = Strike price (exercise price)
  • T = Time to expiration
  • S0 = Current underlying price (time 0)
  • r = Risk-free interest rate
  • σ = Volatility of the underlying asset returns
  • N(·) = Cumulative distribution function of the standard normal distribution

Figure 2 below illustrates the call option value as a function of the underlying asset price. The example considers a European call written on the S&P 500 index, with a strike price of $5,000 and a time to maturity of 30 days. The current price of the underlying index is $6,000, and the risk-free interest rate is set at 3.79% corresponding to the 1-month U.S. Treasury yield, and the volatility is assumed to be 15%.

Figure 2. Call option value as a function of the underlying asset price.
Call option value as a function of the underlying asset price.
Source: computation by the author (BSM model).

Option and volatility

In the Black–Scholes–Merton model, the value of a European call or put option is a monotonically increasing function of volatility. Higher volatility increases the probability of finishing in-the-money while losses remain limited to the option premium, resulting in a strictly positive vega (the first derivative of the option value with respect to volatility) for both calls and puts.

As volatility approaches zero, the option value converges to its intrinsic value, forming a lower bound. With increasing volatility, option values rise toward a finite upper bound equal to the underlying price for calls (and bounded by the strike for puts). An inflection point occurs where volga (the second derivative of the option value with respect to volatility) changes sign: at this point vega is maximized (at-the-money) and declines as the option becomes deep in- or out-of-the-money or as time to maturity decreases.

The upper limit and the lower limit for the call option value function is given below (Hull, 2015, Chapter 11).


Formula for upper and lower limits of the option price

Figure 3 below illustrates the value of a European call option as a function of the underlying asset’s price volatility. The example considers a European call written on the S&P 500 index, with a strike price of $5,000 and a time to maturity of 30 days. The current price of the underlying index is $6,000, and the risk-free interest rate is set at 3.79% corresponding to the 1-month U.S. Treasury yield. A deliberately wide (and economically unrealistic) range of volatility values is employed in order to highlight the theoretical limits of option prices: as volatility tends to infinity, the option value converges to an upper bound ($6,000 in our example), while as volatility approaches zero, the option value converges to a lower bound $1,015.51).

Figure 3. Call option value as a function of price volatility
 Call option value as a function of price volatility
Source: computation by the author (BSM model).

Volatility: the unobservable parameter of the model

When we think of options, the basic equation to remember is “Option = Volatility”. Unlike stocks or bonds, options are not primarily quoted in monetary units (dollars or euros), but rather in terms of implied volatility, expressed as a percentage.

Volatility is not directly observable in financial markets. It is an unobservable (latent) parameter of the pricing model, inferred endogenously from observed option prices through an inversion of the valuation formula given by the BSM model. As a result, option markets are best interpreted as markets for volatility rather than markets for prices.

Out of the five essential parameters of the Black-Scholes-Merton model listed above, the volatility parameter is the unobservable parameter as it is the future fluctuation in price of the underlying asset over the remaining life of the option from the time of observation. Since future volatility cannot be directly observed, practitioners use the inverse of the BSM model to estimate the market’s expectation of this volatility from option market prices, referred to as implied volatility.

Implied Volatility

In practice, implied volatility is the volatility parameter that when input into the Black-Scholes-Merton formula yields the market price of the option and represents the market’s expectation of future volatility.

Calculating Implied volatility

The BSM model maps five input variables (S, K, r, T, σimplied) to a single output variable uniquely: the call option value (Price), such that it’s a bijective function. When the market call option price (CBSM) is known, we invert this relationship using (S, K, r, T, CBSM) as inputs to solve for the implied volatility, σimplied.


Formula for implied volatility

Newton-Raphson Method

As there is no closed form solution to calculate implied volatility from the market price, we need a numerical method such as the Newton–Raphson method to compute it. This involves finding the volatility for which the Black–Scholes–Merton option value CBSM equals the observed market option price CMarket.

We define the function f as the difference between the call option value given by the BSM model and the observed market price of the call option:


Function for the Newton-Raphson method.

Where x represents the unknown variable (implied volatility) to find and CMarket is considered as a constant in the Newton–Raphson method.

Using the Newton-Raphson method, we can iteratively estimate the root of the function, until the difference between two consecutive estimations is less than the tolerance level (ε).


Formula for the iterations in the Newton-Raphson method

In practice, the inflexion point (Tankov, 2006) is taken as the initial guess, because the function f(x) is monotonic, so for very large or very small initial values, the derivative becomes extremely small (see Figure 3), causing the Newton–Raphson update step to overshoot the root and potentially diverge. Selecting the inflection point also minimizes approximation error, as the second derivative of the function at this point is approximately zero, while the first derivative remains non-zero.


Formula for calculating the volatility at inflexion point.

Where σinflection is the volatility at the inflection point.

Figure 4 below illustrates how implied volatility varies with the call option price for different values of the market price (computed using the Newton–Raphson method). As before, the example considers a European call written on the S&P 500 index, with a strike price of $5,000 and a time to maturity of 30 days. The current level of the underlying index is $6,000, and the risk-free interest rate is set at 3.79% corresponding to the 1-month U.S. Treasury yield.

Figure 4. Implied volatility vs. Call Option value
 Implied volatility as a function of call option price
Source: computation by the author.

You can download the Excel file provided below, which contains the calculations and charts illustrating the payoff function, the option price as a function of the underlying asset’s price, the option price as a function of volatility, and the implied volatility as a function of the option price.

Download the Excel file.

You can download the Python code provided below, to calculate the price of a European-style call or put option and calculate the implied volatility from the option market price (BSM model). The Python code uses several libraries.

Download the Python code to calculate the price of a European option.

Alternatively, you can download the R code below with the same functionality as in the Python file.

 Download the R code to calculate the price of a European option.

Why should I be interested in this post?

The seminal Black–Scholes–Merton model was originally developed to price European options. Over time, it has been extended to accommodate a wide range of derivatives, including those based on currencies, commodities, and dividend-paying stocks. As a result, the model is of fundamental importance for anyone seeking to understand the derivatives market and to compute implied volatility as a measure of risk.

Related posts on the SimTrade blog

   ▶ Akshit GUPTA Options

   ▶ Jayati WALIA Black-Scholes-Merton Option Pricing Model

   ▶ Jayati WALIA Implied Volatility

   ▶ Akshit GUPTA Option Greeks – Vega

Useful resources

Academic research

Black F. and M. Scholes (1973) The pricing of options and corporate liabilities. Journal of Political Economy, 81(3), 637–654.

Merton R.C. (1973) Theory of rational option pricing. The Bell Journal of Economics and Management Science, 4(1), 141–183.

Hull J.C. (2022) Options, Futures, and Other Derivatives, 11th Global Edition, Chapter 15 – The Black–Scholes–Merton model, 338–365.

Cox J.C. and M. Rubinstein (1985) Options Markets, First Edition, Chapter 5 – An Exact Option Pricing Formula, 165-252.

Tankov P. (2006) Calibration de Modèles et Couverture de Produits Dérivés (Model calibration and derivatives hedging), Working Paper, Université Paris-Diderot. Available at https://cel.hal.science/cel-00664993/document.

About the BSM model

The Nobel Prize Sveriges Riksbank Prize in Economic Sciences in Memory of Alfred Nobel 1997

Harvard Business School Option Pricing in Theory & Practice: The Nobel Prize Research of Robert C. Merton

Other

NYU Stern Volatility Lab Volatility analysis documentation.

About the author

The article was written in December 2025 by Saral BINDAL (Indian Institute of Technology Kharagpur, Metallurgical and Materials Engineering, 2024-2028 & Research assistant at ESSEC Business School).

   ▶ Discover all articles written by Saral BINDAL

Historical Volatility

Saral BINDAL

In this article, Saral BINDAL (Indian Institute of Technology Kharagpur, Metallurgical and Materials Engineering, 2024-2028 & Research Assistant at ESSEC Business School) explains the concept of historical volatility used in financial markets to represent and measure the changes in asset prices.

Introduction

Volatility in financial markets refers to the degree of variation in an asset’s price or returns over time. Simply put, an asset is considered highly volatile when its price experiences large upward or downward movements, and less volatile when those movements are relatively small. Volatility plays a central role in finance as an indicator of risk and is widely used in various portfolio and risk management techniques.

In practice, the concept of volatility can be operationalized in different ways: historical volatility and implied volatility. Traders and analysts use historical volatility to understand an asset’s past performance and implied volatility as a forward-looking measure of upcoming uncertainties in the market.

Historical volatility measures the actual variability of an asset’s price over a past period, calculated as the standard deviation of its historical returns. Computed over different periods (say a month), historical volatility allows investors to identify trends in volatility and assess how an asset has reacted to market conditions in the past.

Practical Example: Analysis of the S&P 500 Index

Let us consider the S&P 500 index as an example of the calculation of volatility.

Prices

Figure 1 below illustrates the daily closing price of the S&P 500 index over the period from January 2020 to December 2025.

Figure 1. Daily closing prices of the S&P 500 index (2020-2025).
Daily closing prices of the S&P 500 Index (2020-2025)
Source: computation by the author.

Returns

Returns are the percentage gain or loss on the asset’s investment and are generally calculated using one of two methods: arithmetic (simple) or logarithmic (continuously compounded).


Returns Formulas

Where Ri represents the rate of return, and Pi denotes the asset’s price at a given point in time.

The preference for logarithmic returns stems from their property of time-additivity, which simplifies multi-period calculations (the monthly log return is equal to the sum of the daily log returns of the month, which is not the case for arithmetic return). Furthermore, logarithmic returns align with the geometric mean thereby mathematically capture the effects of compounding, unlike arithmetic return, which can overstate performance in volatile markets.

Distribution of returns

A statistical distribution describes the likelihood of different outcomes for a random variable. It begins with classifying the data as either discrete or continuous.

Figure 2 below illustrates the distribution of daily returns for S&P 500 index over the period from January 2020 to December 2025.

Figure 2. Historical distribution of daily returns of the S&P 500 index (2020-2025).
Historical distribution of daily returns of the S&P 500 index (2020-2025)
Source: computation by the author.

Standard deviation of the distribution of returns

In real life, as we do not know the mean and standard deviation of returns, these parameters have to be estimated with data.

The estimator for the mean μ, denoted by μ̂, and the estimator for the variance σ2, denoted by σ̂2, are given by the following formulas:


Formulas for the mean and variance estimators

With the following notations:

  • Ri = rate of return for the ith day
  • μ̂ = estimated mean of the data
  • σ̂2 = estimated variance of the data
  • n = total number of days for the data

These estimators are unbiased and efficient (note the Bessel’s correction for the standard deviation when we divide by (n–1) instead of n).


Unbiased estimators of the mean and variance

For the distribution of returns in Figure 2, the mean and standard deviation calculated using the formulas above are 0.049% and 1.068%, respectively (in daily units).

Annualized volatility

As the usual time frame for human is the year, volatility is often annualized. In order to obtain annual (or annualized) volatility, we scale the daily volatility by the square root of the number of days in that period (τ), as shown below.


Annual Volatility formula

Where  is the number of trading days during the calendar year.

In the U.S. equity market, the annual number of trading days typically ranges from 250 to 255 (252 tradings days in 2025). This variation reflects the holiday calendar: when a holiday falls on a weekday, the exchange closes ; when it falls on a weekend, trading is unaffected. In contrast, the cryptocurrency market has as many trading days as there are calendar days in a year, since it operates continuously, 24/7.

For the S&P 500 index over the period from January 2020 to December 2025, the annualized volatility is given by


 S&P500index Annual Volatility formula

Annualized mean

The calculated mean for the 5-year S&P 500 logarithmic returns is also the daily average return for the period. The annualized average return is given by the formula below.


Annualized mean formula

Where τ is the number of trading days during the calendar year.

For the S&P 500 index over the period from January 2020 to December 2025, the annualized average return is given by


Annualized mean formula

If the value of daily average return is much less than 1, annual average return can be approximated as


Annualized mean value

Application: Estimating the Future Price Range of the S&P 500 index

To develop an intuitive understanding of these figures, we can estimate the one-standard-deviation price range for the S&P 500 index over the next year. From the above calculations, we know that the annualized mean return is 12.534% and the annualized standard deviation is 16.953%.

Under the assumption of normally distributed logarithmic returns, we can say approximately with 68% confidence that the value of S&P 500 index is likely to be in the range of:


Upper and lower limits

If the current value of the S&P 500 index is $6,830, then converting these return estimates into price levels gives:


Upper and lower price limits

Based on a 68% confidence interval, the S&P 500 index is likely to trade in the range of $6,526 to $8,838 over the next year.

Historical Volatility

Historical volatility represents the variability of an asset’s returns over a chosen lookback period. The annualized historical volatility is estimated using the formula below.


 Historical volatility formula

With the following notations:

  • σ = Standard deviation
  • Ri = Return
  • n = total number of trading days in the period (21 for 1 month, 63 for 3 months, etc.)
  • τ = Number of trading days in a calendar year

Volatility calculated over different periods must be annualized to a common timeframe to ensure comparability, as the standard convention in finance is to express volatility on an annual basis. Therefore, when working with daily returns, we annualize the volatility by multiplying it by the square root of 252.

For example, for the S&P 500 index, the annualized historical volatilities over the last 1 month, 3 months, and 6 months, computed on December 3, 2025, are 14.80%, 12.41%, and 11.03%, respectively. The results suggest, since the short term (1 month) volatility is higher than medium (3 months) and long term (6 months) volatility, the recent market movements have been turbulent as compared to the past few months, and due to volatility clustering, periods of high volatility often persist, suggesting that this elevated turbulence may continue in the near term.

Unconditional Volatility

Unconditional volatility is a single volatility number using all historical data, which in our example is the entire five years data; It does not account for the fact that recent market behavior is more relevant for predicting tomorrow’s risk than events from past years, implying that volatility changes over time. It is frequently observed that after any sudden boom or crash in the market, as the storm passes away the volatility tends to revert to a constant value and that value is given by the unconditional volatility of the entire period. This tendency is referred to as mean reversion.

For instance, using S&P 500 index data from 2020 to 2025, the unconditional volatility (annualized standard deviation) is calculated to be 16.952%.

Rolling historical volatility

A single volatility number often fails to capture changing market regimes. Therefore, a rolling historical volatility is usually generated to track the evolution of market risk. By calculating the standard deviation over a moving window, we can observe how volatility has expanded or contracted historically. This is illustrated in Figure 3 below for the annualized 3-month historical volatility of the S&P 500 index over the period 2020-2025.

Figure 3. 3-month rolling historical volatility of the S&P500 index (2020-2025).
3-month rolling historical volatility of the S&P500 index
Source: computation by the author.

In Figure 3, the 3-month rolling historical volatility is plotted along with the unconditional volatility computed over the entire period, calculated using overlapping windows to generate a continuous series. This provides a clear historical perspective, showcasing how the asset’s volatility has fluctuated relative to its long-term average.

For example, during the start of Russia–Ukraine war (February 2022 – August 2022), a noticeable jump in volatility occurred as energy and food prices surged amid fears of supply chain disruptions, given that Russia and Ukraine are major exporters of oil, natural gas, wheat, and other commodities.

The rolling window can be either overlapping or non-overlapping, resulting in continuous or discrete graphs, respectively. Overlapping windows shift by one day, creating a smooth and continuous volatility series, whereas non-overlapping windows shift by one time period, producing a discrete series.

You can download the Excel file provided below, which contains the computation of returns, their historical distribution, the unconditional historical volatility, and the 3-month rolling historical volatility of the S&P 500 index used in this article.

Download the Excel file for returns and volatility calculation

You can download the Python code provided below, which contains the computation of returns, first four moments of the distribution, and experiment with the x-month rolling historical volatility function to visualize the evolution of historical volatility over time.

Download the Python code for returns and volatility calculation.

Alternatively, you can download the R code below with the same functionality as in the Python file.

Download the R code for returns and volatility calculation.

Alterative measures of volatility

We now mention a few other ways volatility can be measured: Parkinson volatility, Implied volatility, ARCH model, and stochastic volatility model.

Parkinson volatility

The Parkinson model (1980) uses the highest and lowest prices during a given period (say a month) for the purpose of measurement of volatility. This model is a high-low volatility measure, based on the difference between the maximum and minimum prices observed during a certain period.

Parkinson volatility is a range-based variance estimator that replaces squared returns with the squared high–low log price range, scaled to remain unbiased. It assumes a driftless (expected growth rate of the stock price equal to zero) geometric Brownian motion, it is five times more efficient than close-to-close returns because it accounts for fluctuation of stock price within a day.

For a sample of n observations (say days), the Parkinson volatility is given by


Parkinson Volatility formula

where:

  • Ht is the highest price on period t
  • Lt is the lowest price on period t

Implied volatility

Implied Volatility (IV) is the level of volatility for the underlying asset that, when plugged into an option pricing model such as Black–Scholes–Merton, makes the model’s theoretical option price equal to the option’s observed market price.

It is a forward looking measure because it reflects the market’s expectation of how much the underlying asset’s price is likely to fluctuate over the remaining life of the option, rather than how much it has moved in the past.

The Chicago Board Options Exchange (CBOE), a leading global financial exchange operator provides implied volatility indices like the VIX and Implied Correlation Index, measuring 30-day expected volatility from SPX options. These are used by traders to gauge market fear, speculate via futures/options/ETPs, hedge equity portfolios and manage risk during volatility spikes.

ARCH model

Autoregressive Conditional Heteroscedasticity (ARCH) models address time-varying volatility in time series data. Introduced by Engle in 1982, ARCH models look at the size of past shocks to estimate how volatile the next period is likely to be. If recent movements were big, the model expects higher volatility; if they were small, it expects lower volatility justifying the idea of volatility clustering. Originally applied to inflation data, this model has been widely used in to model financial data.

ARCH model capture volatility clustering, which refers to an observation about how volatility behaves in the short term, a large movement is usually followed by another large movement, thus volatility is predictable in the short term. Historical volatility gives a short-term hint of the near future changes in the market because recent noise often continues.

Generalized Autoregressive Conditional Heteroscedasticity (GARCH) extends ARCH by past predicted volatility, not just past shocks, as refined by Bollerslev in 1986 from Engle’s work. Both of these methods are more accurate methods to forecast volatility than what we had discussed as they account for the time varying nature of volatility.

Stochastic volatility models

In practice, volatility is time-varying: it exhibits clustering, persistence, and mean reversion. To capture these empirical features, stochastic volatility (SV) models treat volatility not as a constant parameter but as a stochastic process jointly evolving with the asset price. Among these models, the Heston (1993) specification is one of the most influential.

The Heston model assumes that the asset price follows a diffusion process analogous to geometric Brownian motion, while the instantaneous variance evolves according to a mean-reverting square-root process. Moreover, the innovations to the price and variance processes are correlated, thereby capturing the leverage effect frequently observed in equity markets.

Applications in finance

This section covers key mathematical concepts and fundamental principles of portfolio management, highlighting the role of volatility in assessing risk.

The normal distribution

The normal distribution is one of the most commonly used probability distribution of a random variable with a unimodal, symmetric and bell-shaped curve. The probability distribution function for a random variable X following a normal distribution with mean μ and variance σ2 is given by


Normal distribution function

A random variable X is said to follow standard normal distribution if its mean is zero and variance is one.

The figure below represents the confidence intervals, showing the percentage of data falling within one, two, and three standard deviations from the mean.

Figure 4. Probability density function and confidence intervals for a standard normal varaible.
Standard normal distribution” width=
Source: computation by the author

Brownian motion

Robert Brown first observed Brownian motion was as the erratic and random movement of pollen particles suspended in water due to constant collision with water molecules. It was later formulated mathematically by Norbert Wiener and is also known as the Wiener process.

The random walk theory suggests that it’s impossible to predict future stock prices as they move randomly, and when the timestep of this theory becomes infinitesimally small it becomes, Brownian Motion.

In the context of financial stochastic process, when the market is modeled by the standard Brownian motion, the probability distribution function of the future price is a normal distribution, whereas when modeled by Geometric Brownian Motion, the future prices are said to be lognormally distributed. This is also called the Brownian Motion hypothesis on the movement of stock prices.

The process of a standard Brownian motion is given by:


Standard Brownian motion formula.

The process of a geometric Brownian motion is given by:


Geometric Brownian motion formula.

Where, dSt is the change in asset price in continuous time dt, dXt is a random variable from the normal distribution (N (0, 1)) or Wiener process at a time t, σ represents the price volatility, and μ represents the expected growth rate of the asset price, also known as the ‘drift’.

Modern Portfolio Theory (MPT)

Modern Portfolio Theory (MPT), developed by Nobel Laureate, Harry Markowitz, in the 1950s, is a framework for constructing optimal investment portfolios, derived from the foundational mean-variance model.

The Markowitz mean–variance model suggests that risk can be reduced through diversification. It proposes that risk-averse investors should optimize their portfolios by selecting a combination of assets that balances expected return and risk, thereby achieving the best possible return for the level of risk they are willing to take. The optimal trade-off curve between expected return and risk, commonly known as the efficient frontier, represents the set of portfolios that maximizes expected return for each level of standard deviation (risk).

Capital Asset Pricing Model (CAPM)

The Capital Asset Pricing Model (CAPM) builds on the model of portfolio choice developed by Harry Markowitz (1952), stated above. CAPM states that, assuming full agreement on return distributions and either risk-free borrowing/lending or unrestricted short selling, the value-weighted market portfolio of risky assets is mean-variance efficient, and expected returns are linear in the market beta.

The main result of the CAPM is a simple mathematical formula that links the expected return of an asset to its risk measured by the beta of the asset:


CAPM formula

Where:

  • E(Ri) = expected return of asset i
  • Rf = risk-free rate
  • βi = measure of the risk of asset i
  • E(Rm) = expected return of the market
  • E(Rm) − Rf = market risk premium

CAPM recognizes that an asset’s total risk has two components: systematic risk and specific risk, but only systematic risk is compensated in expected returns.

Returns decomposition fromula.
 Returns decomposition fromula

Where the realized (actual) returns of the market (Rm) and the asset (Ri) exceed their expected values only because of consideration of systematic risk (ε).

Decomposition of risk.
Decompositionion of risk

Systematic risk is a macro-level form of risk that affects a large number of assets to one degree or another, and therefore cannot be eliminated. General economic conditions, such as inflation, interest rates, geopolitical risk or exchange rates are all examples of systematic risk factors.

Specific risk (also called idiosyncratic risk or unsystematic risk), on the other hand, is a micro-level form of risk that specifically affects a single asset or narrow group of assets. It involves special risk that is unconnected to the market and reflects the unique nature of the asset. For example, company specific financial or business decisions which resulted in lower earnings and affected the stock prices negatively. However, it did not impact other asset’s performance in the portfolio. Other examples of specific risk might include a firm’s credit rating, negative press reports about a business, or a strike affecting a particular company.

Why should I be interested in this post?

Understanding different measures of volatility, is a pre-requisite to better assess potential losses, optimize portfolio allocation, and make informed decisions to balance risk and expected return. Volatility is fundamental to risk management and constructing investment strategies.

Related posts on the SimTrade blog

Risk and Volatility

   ▶ Jayati WALIA Brownian Motion in Finance

   ▶ Youssef LOURAOUI Systematic Risk

   ▶ Youssef LOURAOUI Specific Risk

   ▶ Jayati WALIA Implied Volatility

   ▶ Mathias DUMONT Pricing Weather Risk

   ▶ Jayati WALIA Black-Scholes-Merton Option Pricing Model

Portfolio Theory and Models

   ▶ Jayati WALIA Returns

   ▶ Youssef LOURAOUI Portfolio

   ▶ Jayati WALIA Capital Asset Pricing Model (CAPM)

   ▶ Youssef LOURAOUI Optimal Portfolio

Financial Indexes

   ▶ Nithisha CHALLA Financial Indexes

   ▶ Nithisha CHALLA Calculation of Financial Indexes

   ▶ Nithisha CHALLA The S&P 500 Index

Useful Resources

Academic research

Bollerslev, T. (1986). Generalized Autoregressive Conditional Heteroskedasticity, Journal of Econometrics, 31(3), 307–327.

Engle, R. F. (1982). Autoregressive Conditional Heteroscedasticity with Estimates of the Variance of United Kingdom Inflation, Econometrica, 50(4), 987–1007.

Fama, E. F., & French, K. R. (2004). The Capital Asset Pricing Model: Theory and Evidence, Journal of Economic Perspectives, 18(3), 25–46.

Heston, S. L. (1993). A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options, The Journal of Finance, 48(3), 1–24.

Markowitz, H. M. (1952). Portfolio Selection, The Journal of Finance, 7(1), 77–91.

Parkinson, M. (1980). The extreme value method for estimating the variance of the rate of return. Journal of Business, 53(1), 61–65.

Sharpe, W. F. (1964). Capital Asset Prices: A Theory of Market Equilibrium under Conditions of Risk, The Journal of Finance, 19(3), 425–442.

Tsay, R. S. (2010). Analysis of financial time series, John Wiley & Sons.

Other

NYU Stern Volatility Lab Volatility analysis documentation.

Extreme Events in Finance Risk maps: extreme risk, risk and performance.

About the author

The article was written in December 2025 by Saral BINDAL (Indian Institute of Technology Kharagpur, Metallurgical and Materials Engineering, 2024-2028 & Research Assistant at ESSEC Business School).

   ▶ Discover all articles written by Saral BINDAL

Pricing Weather Risk: How to Value Agricultural Derivatives with Climate-Based Volatility Inputs

Mathias DUMONT

In this article, Mathias DUMONT (ESSEC Business School, Global Bachelor in Business Administration (GBBA), 2022-2026) explains how weather risk impacts the pricing of agricultural derivatives like futures and options, and how climate-based data can be integrated into stochastic pricing models. Combining academic insights and practical examples, including a mini-case from the SimTrade Blé de France simulation, the article illustrates adjustments to models such as the Black-Scholes-Merton model for temperature and rainfall variables in valuing agricultural contracts.

Introduction

Extreme weather has always been a critical factor in agriculture, but climate change is amplifying the frequency and severity of these events. From prolonged droughts to unseasonal floods, weather shocks can send crop yields and commodity prices on wild rides. This rising uncertainty has given birth to weather derivatives – financial instruments designed to hedge weather-related risks – and has made volatility forecasting a key challenge in pricing agricultural contracts. In fact, as businesses grapple with climate volatility, trading volume in weather derivatives has surged. CME Group saw a 260% increase last year (CME Group, 2023). The question for traders and risk managers is: how do we quantitatively factor weather risk into the pricing of futures and options on crops like wheat and corn?

Weather Risk and Agricultural Markets

Weather directly affects crop supply. A bumper harvest following ideal weather can flood the market and depress prices, whereas a drought or frost can decimate yields and trigger price spikes. These supply swings translate into volatility for agricultural commodity markets. For example, during the U.S. drought of 2012, corn prices skyrocketed, and the implied volatility of corn futures jumped by over 14 percentage points within a month, reaching ~49% in mid-July. Such surges reflect the market rapidly repricing risk as participants absorb new climate information (in this case, worsening crop prospects). Seasonal patterns are also evident: harvest seasons tend to coincide with higher price volatility because that’s when weather uncertainty is at its peak. Studies show that harvesting cycles create predictable seasonal volatility patterns in crop markets – when a critical growth period is underway, any shift in rainfall or temperature forecasts can send prices swinging.

Beyond affecting supply quantity, weather can influence crop quality (e.g., excessive rain can spoil grain quality) and even logistic costs (flooded transport routes, etc.), further feeding into prices. The interconnected global nature of agriculture means a drought in one region can reverberate worldwide. As noted in the SimTrade Blé de France case, weather conditions in France influence the quantity and quality of wheat the company harvests, while weather conditions around the world influence the international wheat price. In the Blé de France simulation (which models a French wheat producer’s stock), participants see how news of floods or droughts translate into stock price moves. For instance, the company might project a 7-million-ton wheat harvest, but analysts’ forecasts range from 6.5 to 7.2 Mt – with the realized level highly weather-dependent in the final weeks of the season. A poor weather turn not only shrinks the crop but boosts global wheat prices, creating a complex revenue impact on the firm. This mini-case underlines that weather risk entails both volume uncertainty and price uncertainty, a double-whammy for agricultural firms and their investors.

Case Study: Weather Shocks in Wheat Markets

To illustrate the impact of weather risk on commodity pricing, consider three simulated scenarios for an upcoming wheat growing season: (1) **Favorable weather**, (2) **Moderate conditions**, and (3) **Severe weather** such as drought. Each scenario generates a distinct price trajectory in the wheat market. Under favorable weather, prices tend to remain stable or decline slightly, particularly at harvest, due to strong yields and potential oversupply. In moderate conditions, prices may rise modestly as the market adjusts to balanced supply and demand. In contrast, severe weather triggers early price rallies as concerns about yield shortfalls emerge, followed by sharp spikes once crop damage becomes evident. For producers and traders, anticipating these divergent price paths is essential for pricing contracts, managing risk exposure, and structuring hedging strategies effectively.

Figure 1. Simulated commodity price paths under three weather scenarios.
Simulated Price Paths
Source: Author’s simulation.

Figure 1. shows the simulation of commodity price paths under three weather scenarios: severe weather (red), moderate weather (orange), and favorable weather (green). A mid-season weather forecast alert (Day 15) triggers a shift in market expectations, causing price divergence. This simulation illustrates how weather shocks and forecasts impact commodity pricing through volatility and revised yield expectations.

From a risk management perspective, tools exist to handle these contingencies. Farmers or firms concerned about catastrophic weather can turn to weather derivatives for protection. Weather derivatives are financial contracts (often based on indexes like temperature or rainfall levels) that pay out based on specific weather outcomes, allowing businesses to offset losses caused by adverse conditions. They have been used by a wide range of players – from utilities hedging warm winters, to breweries hedging late frosts. These instruments can be customized over-the-counter or traded on exchanges. Notably, CME Group lists standardized weather futures and options tied to indices such as heating degree days (HDD) and cooling degree days (CDD) for various cities. The existence of such contracts means that even when commodity producers cannot fully insure their crop yield, they might hedge certain aspects of weather risk (like an unusually hot summer) via financial markets. In our context, a wheat farmer worried about drought could, say, buy a weather option that pays off if rainfall falls below a threshold, providing funds when their crop output (and thus futures position) suffers.

Climate-Based Volatility in Derivatives Pricing

How can weather uncertainty be incorporated into derivative pricing models? Classic option pricing, such as the Black-Scholes-Merton model, assumes a fixed volatility for the underlying asset’s returns. For agricultural commodities, that volatility is anything but constant – it ebbs and flows with the weather and seasonal progress. Practitioners thus often use stochastic volatility models or at least adjust the volatility input over time. For example, one might use higher volatility estimates during the crop’s growing season and lower volatility post-harvest when output is known. This practice parallels how equity traders anticipate higher volatility in stock prices ahead of major earnings or profit announcements, and lower volatility after the announcement of profits by the firm.

Like companies facing performance surprises, weather shocks inject information asymmetry into the market, which must be priced into the option premiums. This aligns with the observed Samuelson effect, where futures contracts on commodities tend to have higher volatility when they are near maturity (coinciding with harvest uncertainty).

Market prices of options themselves reflect these expectations. When a looming weather event is expected to cause turmoil, options premiums will rise. The metric capturing this is implied volatility – the volatility level implied by current option prices. Implied vol is essentially forward-looking and will jump if traders foresee choppy waters ahead. Empirical evidence shows that extreme weather forecasts translate into higher implied vols for crop options. In 2012, as drought fears intensified, corn option implied volatility spiked (alongside futures prices). Conversely, once a forecasted drought started being relieved by rains, implied volatility eased off, signaling that some uncertainty had been resolved. A recent study also found that integrating meteorological data (like rainfall and temperature anomalies) into volatility modeling significantly improves the ability to hedge risk in agricultural markets. In other words, the more information we feed into our models about the climate, the more accurately we can price and hedge these derivatives.

Figure 2. Implied Volatility of Crop Options Over Time with Weather Events
Line chart showing implied volatility of crop options over 12 months with spikes linked to weather events
Source: Author’s simulation.

This simulation illustrates the evolution of implied volatility over a 12-month crop cycle. Forecasted climate events—drought (Month 3), frost (Month 6), heatwave (Month 8), and rainfall shortage (Month 11)—lead to moderate but distinct volatility spikes. As uncertainty resolves, volatility returns to baseline.

One practical approach to pricing under climate uncertainty is to use scenario-based or simulation-based models. Instead of assuming a single volatility number, an analyst can simulate thousands of possible weather outcomes (perhaps using historical climate data or meteorological forecast models) and the corresponding price paths for the commodity. Each simulated price path yields a payoff for the derivative (e.g. an option’s payoff at expiration), and by averaging those payoffs (and discounting appropriately), one can derive a weather-adjusted theoretical price. This Monte Carlo style approach effectively treats weather as an external random factor influencing the commodity’s drift and volatility. It’s particularly useful for complex derivatives or when the payoff depends explicitly on weather indices (such as a derivative that pays out if rainfall is below X mm).

When the derivative’s underlying is the commodity itself (e.g. a corn futures option), traditional risk-neutral pricing arguments still apply, but the challenge is forecasting volatility. Traders often adjust the volatility smile/skew on agricultural options to account for asymmetric weather risks – for instance, if a drought can cause a much bigger upside move than a rainy season can cause a downside move, call options might embed a higher implied volatility (reflecting that upside risk of price spikes). This is observed in practice as well; extreme weather events can distort the implied volatility “skew” of crop options, as out-of-the-money calls become more sought after as disaster insurance.

In contrast, if the derivative’s underlying is a pure weather index (say an option on cumulative rainfall), then pricing becomes more complex because the underlying (rainfall) is not a tradable asset. In such cases, the Black-Scholes-Merton formula is not directly applicable. Instead, pricing relies on actuarial or risk-neutral methodologies that incorporate a market price of risk for weather. For example, one method is to estimate the probability distribution of the weather index from historical data, then add a risk premium to account for investors’ risk aversion to weather variability, and discount expected payoffs accordingly. Another method uses “burn analysis” – taking historical weather outcomes and the associated financial losses/gains had the derivative been in place, to gauge a fair premium. Academic research has proposed models ranging from modified Black-Scholes-Merton-type formulas for rainfall (with adjustments for the non-tradability) to advanced statistical models (e.g. Ornstein-Uhlenbeck processes with seasonality for temperature indices. The key takeaway is that whether it’s directly in commodity options or in dedicated weather derivatives, climate factors force us to go beyond textbook models and embrace more dynamic, data-driven pricing techniques.

Why should I be interested in this post?

For an ESSEC student or a young finance professional, this topic sits at the intersection of finance and real-world impact. Understanding weather risk in markets is not just about farming – it’s about how big data and climate science are increasingly intertwined with financial strategy. Agricultural commodities remain a cornerstone of the global economy, and volatility in these markets can affect food prices, inflation, and even economic stability in various countries. By grasping how to value derivatives with climate-based volatility inputs, you are gaining insight into a growing niche of finance that deals with sustainability and risk management. Moreover, the skills involved – scenario analysis, simulation modeling, blending of economic and scientific data – are highly transferable to other domains (think energy markets or any sector where uncertainty reigns). In a world facing climate change, expertise in weather-related financial products could open career opportunities in commodity trading desks, insurance/reinsurance firms, or specialized hedge funds. Ultimately, this post encourages you to think creatively and interdisciplinarily: the best hedging or valuation solutions may come from combining financial theory with environmental intelligence.

Related posts on the SimTrade blog

   ▶ Camille KELLER Coffee Futures: The Economic and Environmental Drivers Behind Rising Prices

   ▶ Jayati WALIA Implied Volatility

   ▶ Akshit GUPTA Futures Contract

   ▶ Anant JAIN Understanding Price Elasticity of Demand

Useful resources

Chicago Mercantile Exchange (CME) Weather futures and options product information. (Exchange-traded weather derivative contracts on temperature and other indices)

U.S. Energy Information Administration Drought increases price of corn, reduces profits to ethanol producers (2012). (Article discussing the 2012 drought’s impact on corn prices and volatility)

Nature Communications (2024) Financial markets value skillful forecasts of seasonal climate. (Research showing that seasonal climate outlooks have measurable effects on implied volatility and market uncertainty)

Das, S. et al. (2025) Predicting and Mitigating Agricultural Price Volatility Using Climate Scenarios and Risk Models. (Academic study demonstrating the integration of climate data into volatility models and using Black-Scholes to value a government price support as a put option)

Pai, J. & Zheng, Z. (2013) Pricing Temperature Derivatives with a Filtered Historical Simulation Approach. (Discussion of why Black-Scholes is not directly applicable to weather derivatives and alternative pricing approaches)

About the author

The article was written in May 2025 by Mathias DUMONT (ESSEC Business School, Global Bachelor in Business Administration (GBBA), 2022-2026).

Implied Volatility

Jayati WALIA

In this article, Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022) explains how implied volatility is computed from option market prices and a option pricing model.

Introduction

Volatility is a measure of fluctuations observed in an asset’s returns over a period of time. The standard deviation of historical asset returns is one of the measures of volatility. In option pricing models like the Black-Scholes-Merton model, volatility corresponds to the volatility of the underlying asset’s return. It is a key component of the model because it is not directly observed in the market and cannot be directly computed. Moreover, volatility has a strong impact on the option value.

Mathematically, in a reverse way, implied volatility is the volatility of the underlying asset which gives the theoretical value of an option (as computed by Black-Scholes-Merton model) equal to the market price of that option.

Implied volatility is a forward-looking measure because it is a representation of expected price movements in an underlying asset in the future.

Computation methods for implied volatility

The Black-Scholes-Merton (BSM) model provides an analytical formula for the price of both a call option and a put option.

The value for a call option at time t is given by:

 Call option value

The value for a put option at time t is given by:

Put option value

where the parameters d1 and d2 are given by:,

call option d1 d2

with the following notations:

St : Price of the underlying asset at time t
t: Current date
T: Expiry date of the option
K: Strike price of the option
r: Risk-free interest rate
σ: Volatility of the underlying asset
N(.): Cumulative distribution function for a normal (Gaussian) distribution. It is the probability that a random variable is less or equal to its input (i.e. d₁ and d₂) for a normal distribution. Thus, 0 ≤ N(.) ≤ 1

From the BSM model, both for a call option and a put option, the option price is an increasing function of the volatility of the underlying asset: an increase in volatility will cause an increase in the option price.

Figures 1 and 2 below illustrate the relationship between the value of a call option and a put option and the level of volatility of the underlying asset according to the BSM model.

Figure 1. Call option value as a function of volatility.
Call option value as a function of volatility
Source: computation by the author (BSM model)

Figure 2. Put option value as a function of volatility.
Put option value as a function of volatility
Source: computation by the author (BSM model)

You can download below the Excel file for the computation of the value of a call option and a put option for different levels of volatility of the underlying asset according to the BSM model.

Excel file to compute the option value as a function of volatility

We can observe that the call and put option values are a monotonically increasing function of the volatility of the underlying asset. Then, for a given level of volatility, there is a unique value for the call option and a unique value for the put option. This implies that this function can be reversed; for a given value for the call option, there is a unique level of volatility, and similarly, for a given value for the put option, there is a unique level of volatility.

The BSM formula can be reverse-engineered to compute the implied volatility i.e., if we have the market price of the option, the market price of the underlying asset, the market risk-free rate, and the characteristics of the option (the expiration date and strike price), we can obtain the implied volatility of the underlying asset by inverting the BSM formula.

Example

Consider a call option with a strike price of 50 € and a time to maturity of 0.25 years. The market risk-free interest rate is 2% and the current price of the underlying asset is 50 €. Thus, the call option is ‘at-the-money’. If the market price of the call option is equal to 2 €, then the associated level of volatility (implied volatility) is equal to 18.83%.

You can download below the Excel file below to compute the implied volatility given the market price of a call option. The computation uses the Excel solver.

Excel file to compute implied volatility of an option

Volatility smile

Volatility smile is the name given to the plot of implied volatility against different strikes for options with the same time to maturity. According to the BSM model, it is a horizontal straight line as the model assumes that the volatility is constant (it does not depend on the option strike). However, in practice, we do not observe a horizontal straight line. The curve may be in the shape of the alphabet ‘U’ or a ‘smile’ which is the usual term used to refer to the observed function of implied volatility.

Figure 3 below depicts the volatility smile for call options on the Apple stock on May 13, 2022.

Figure 3. Volatility smile for call options on Apple stock.
Apple volatility smile
Source: Computation by author.

Excel file for implied volatility from Apple stock option

We can also observe that the for a specific time to maturity, the implied volatility is minimum when the option is at-the-money.

Volatility surface

An essential assumption of the BSM model is that the returns of the underlying asset follow geometric Brownian motion (corresponding to log-normal distribution for the price at a given point in time) and the volatility of the underlying asset price remains constant over time until the expiration date. Thus theoretically, for a constant time to maturity, the plot of implied volatility and strike price would be a horizontal straight line corresponding to a constant value for volatility.

Volatility surface is obtained when values for implied volatilities are calculated for options with different strike prices and times to maturity.

CBOE Volatility Index

The Chicago Board Options Exchange publishes the renowned Volatility Index (also known as VIX) which is an index based on the implied volatility of 30-day option contracts on the S&P 500 index. It is also called the ‘fear gauge’ and it is a representation of the market outlook for volatility for the next 30 days.

Related posts on the SimTrade blog

   ▶ All posts about Options

   ▶ Akshit GUPTA Options

   ▶ Jayati WALIA Brownian Motion in Finance

   ▶ Jayati WALIA Brownian Motion in Finance

   ▶ Youssef LOURAOUI Minimum Volatility Factor

   ▶ Youssef LOURAOUI VIX index

Useful resources

Academic articles

Black F. and M. Scholes (1973) “The Pricing of Options and Corporate Liabilities” The Journal of Political Economy, 81, 637-654.

Dupire B. (1994). “Pricing with a Smile” Risk Magazine 7, 18-20.

Merton R.C. (1973) “Theory of Rational Option Pricing” Bell Journal of Economics, 4, 141–183.

Business

CBOE Volatility Index (VIX)

CBOE VIX tradable products

About the author

The article was written in May 2022 by Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022).

Standard deviation

Jayati WALIA

In this article, Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022) presents an overview of standard deviation and its use in financial markets.

Definition

The standard deviation is a measure that indicates how much data scatter around the mean. The idea is to measure how an observation deviates from the mean on average./p>

Mathematical formulae

The first step to compute the standard deviation is to compute the mean. Considering a variable X, the arithmetic mean of a data set with N observations, X1, X2 … XN, is computed as:

img_arithmetic_mean

In the data set analysis, we also consider the dispersion or variability of data values around the central tendency or the mean. The variance of a data set is a measure of dispersion of data set values from the (estimated) mean and can be expressed as:

variance

Note that in the above formula we divide by N-1 because the mean is not known but estimated (usual case in finance). If the mean is known with certainty (when dealing the whole population not a sample), then we divide by N.

A problem with variance, however, is the difficulty of interpreting it due to its squared unit of measurement. This issue is resolved by using the standard deviation, which has the same measurement unit as the observations of the data set (such as percentage, dollar, etc.). The standard deviation is computed as the square root of variance:

standard deviation

A low value standard deviation indicates that the data set values tend to be closer to the mean of the set and thus lower dispersion, while a high standard deviation indicates that the values are spread out over a wider range indication higher dispersion.

Measure of volatility

For financial investments, the X variable in the above formulas would correspond to the return on the investment computed on a given period of time. We usually consider the trade-off between risk and reward. In this context, the reward corresponds to the expected return measured by the mean, and the risk corresponds to the standard deviation of returns.

In financial markets, the standard deviation of asset returns is used as a statistical measure of the risk associated with price fluctuations of any particular security or asset (such as stocks, bonds, etc.) or the risk of a portfolio of assets (such as mutual funds, index mutual funds or ETFs, etc.).

Investors always consider a mathematical basis to make investment decisions known as mean-variance optimization which enables them to make a meaningful comparison between the expected return and risk associated with any security. In other words, investors expect higher future returns on an investment on average if that investment holds a relatively higher level of risk or uncertainty. Standard deviation thus provides a quantified estimate of the risk or volatility of future returns.

In the context of financial securities, the higher the standard deviation, the greater is the dispersion between each return and the mean, which indicates a wider price range and hence greater volatility. Similarly, the lower the standard deviation, the lesser is the dispersion between each return and the mean, which indicates a narrower price range and hence lower volatility for the security.

Example: Apple Stock

To illustrate the concept of volatility in financial markets, we use a data set of Apple stock prices. At each date, we compute the volatility as the standard deviation of daily stock returns over a rolling window corresponding to the past calendar month (about 22 trading days). This daily volatility is then annualized and expressed as a percentage.

Figure 1. Stock price and volatility of Apple stock.

price and volatility for Apple stock
Source: computation by the author (data source: Bloomberg).

You can download below the Excel file for the calculation of the volatility of stock returns. The data used are for Apple for the period 2020-2021.

ownload the Excel file to compute the volatility of stock returns

Related posts on the SimTrade blog

▶ Jayati WALIA Quantitative Risk Management

▶ Jayati WALIA Value at Risk

▶ Jayati WALIA Brownian Motion in Finance

Useful resources

Wikipedia Standard Deviation

About the author

The article was written in November 2021 by Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022).

VIX index

Youssef_Louraoui

In this article, Youssef LOURAOUI (ESSEC Business School, Global Bachelor of Business Administration, 2017-2021) presents the VIX index, which is a financial index that measures the uncertainty in the US equity market.

This article is structured as follows: we begin by defining the grounding notions of the VIX index. We then explain the behavior of this index and its statistical characteristics. We finish by presenting its practical usage in financial markets.

Definition

The CBOE Volatility Index, abbreviated “VIX”, is a measure of the expected S&P 500 index movement calculated by the Chicago Board Options Exchange (CBOE) from the current trading prices of options written on the S&P 500 index.

Known as Wall Street’s “fear index”, the VIX is closely monitored by a broad range of market players, and its level and pattern have become ingrained in market discussion.

Figure 1 illustrates the evolution of the VIX index for the period from 2003 to 2021.
Figure 1 Historical levels of the VIX index from 2003-2021.
VIX_levels_analysis
Source: computation by the author (Data source: Thomson Reuters).

VIX values greater than 20 are regarded to be high by market participants. If the VIX is between 12 and 20, it is considered normal; if it is less than 12, it is considered low. As it is the case with other indices, the VIX is computed using the price of a basket of tradable components (in this case, options expiring within the next month or so). The profit or loss that option buyers and sellers realize during the option’s life will depend, among other things, on how significantly the S&P 500’s actual volatility will differ from the implied volatility given by the VIX at the start of the period (S&P Global Research, 2017).

Behavior of the VIX index

Statistical distribution of the S&P500 index returns and VIX level

Figure 2 displays the statistical distribution of the price variations in the S&P500 index for different levels of the VIX index The higher the VIX index (by convention, greater than 20), the more severe the distribution tends to be, with negative skewness and high kurtosis indicating heightened volatility in the US market, therefore exacerbating both positive and negative swings. An opposite finding may be made for the VIX level at lower levels (often less than 12), when market swings are less evident due to less skewness and lower kurtosis (S&P Global Research, 2017).

Figure 2. The distribution of 30-day return in the S&P500 index for different VIX index levels.
Statistical distribution of the S&P500 index returns
Source: S&P Global Research (2017).

If the VIX is low, market players may benefit by purchasing options; conversely, if the VIX is high, market participants may profit from selling options. The specific utility of anticipated VIX is that it gives us with a more accurate assessment of whether VIX is high, low, or normal at any point in time (S&P Global Research, 2017). Thus, VIX may be regarded of as a crowd-sourced estimate of the S&P 500’s expected volatility. As with interest rates and dividends, one cannot invest directly in them, even though one can guess on their future worth, one cannot invest directly in VIX, and the significance of a specific VIX level is commonly misinterpreted (S&P Global Research, 2017).

Recent volatility in the S&P500 index and VIX level

Figure 3 demonstrates that the VIX index is strongly correlated with recent market volatility. However, there is considerable variance; for example, a recent volatility level of about 20% has been associated with a VIX level of 34 (point B, when VIX was very “high”) and with a VIX level of 12 (point C, when VIX was relatively “low”). Volatility (realized or implied) has a strong propensity to return to its mean. This insight is not especially original, despite its illustrious past. There is an enormous body of data demonstrating that volatility tends to mean revert across markets, and the pioneers of this field were given the Nobel Prize in part for incorporating their results into volatility forecasts and simulations (S&P Global Research, 2017).

Figure 3. Relation between VIX and recent volatility.
VIX_regression_analysis
Source: S&P Global Research (2017).

Realized volatility in the S&P500 index and VIX level

Figure 4 represents the relationship between Realized volatility in the S&P500 index over a period and the VIX level at the begining of the period.

Figure 4. VIX versus next realized volatility.
VIX_realized_graph
Source: S&P Global Research (2017).

Mean reversion

Figure 5 shows how VIX index converge to a certain llong-term level as time passes. This finding is not due to 15% being exceptional in any manner; this figure for M was calculated using historical volatility levels for the S&P 500 and their evolution. It is not implausible that M (else referred to as long-term average volatility in the US equities market) may change over time; changes in the S&P 500’s sector weightings, trade All of these factors have the ability to influence both the pace and the volume and the point at which mean reversion occurs.

Figure 5. Mean-reversion dynamic in recent volatility.
VIX mean reversion
Source: S&P Global Research (2017).

Use of the VIX index in financial markets

There are two methods for determining an asset’s volatility. Either through a statistical calculation of an asset’s realized volatility, also known as historical volatility, which serves as a pointer to the asset’s volatility behavior. This is a limited method that is based on the premise that past volatility tends to replicate itself in the future, without including a forward-looking study of volatility. The second technique is to extract an asset’s volatility from option prices referred to as “implied volatility”.

Why should I be interested in this post?

When investors make investment decisions, they utilize the VIX to gauge the degree of risk, worry, or stress in the market. Additionally, traders can trade the VIX using a range of options and exchange-traded products, or price derivatives using VIX values.

Related posts on the SimTrade blog

All posts about Options

▶ Akshit GUPTA Options

▶ Akshit GUPTA History of Option Markets

▶ Jayati WALIA Implied Volatility

▶ Youssef LOURAOUI Minimum Volatility Factor

Useful resources

Business analysis

CBOE , 2021. VIX

Nasdaq, 2021. Realized Volatility

Nasdaq, 2021. Vix Index Volatility

S&P Global Research, 2017. Reading VIX: Does VIX Predict Future Volatility?

S&P Global Research, 2017. A Practitioner’s Guide to Reading VIX

About the author

The article was written in September 2021 by Youssef LOURAOUI (ESSEC Business School, Global Bachelor of Business Administration, 2017-2021).

Brownian Motion in Finance

Jayati WALIA

In this article, Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022) explains the Brownian motion and its applications in finance to model asset prices like stocks traded in financial markets.

Introduction

Stock price movements form a random pattern. The prices fluctuate everyday resulting from market forces like supply and demand, company valuation and earnings, and economic factors like inflation, liquidity, demographics of country and investors, political developments, etc. Market participants try to anticipate stock prices using all these factors and contribute to make price movements random by their trading activities as the financial and economics worlds are constantly changing.

What is a Brownian Motion?

The Brownian motion was first introduced by botanist Robert Brown who observed the random movement of pollen particles due to water molecules under a microscope. It was in the 1900s that the French mathematician Louis Bachelier applied the concept of Brownian motion to asset price behavior for the first time, and this led to Brownian motion becoming one of the most important fundamental of modern quantitative finance. In Bachelier’s theory, price fluctuations observed over a small time period are independent of the current price along with historical behavior of price movements. Combining his assumptions with the Central Limit Theorem, he also deduces that the random behavior of prices can be said to be represented by a normal distribution (Gaussian distribution).

This led to the development of the Random Walk Hypothesis or Random Walk Theory, as it is known today in modern finance. A random walk is a statistical phenomenon wherein stock prices move randomly.

When the time step of a random walk is made infinitesimally small, the random walk becomes a Brownian motion.

Standard Brownian Motion

In context of financial stochastic processes, the Brownian motion is also described as the Wiener Process that is a continuous stochastic process with normally distributed increments. Using the Wiener process notation, an asset price model in continuous time can be expressed as:

brownian motion equation

with dS being the change in asset price in continuous time dt. dX is the random variable from the normal distribution (N(0, 1) or Wiener process). σ is assumed to be constant and represents the price volatility considering the unexpected changes that can result from external effects. μdt together represents the deterministic return within the time interval with μ representing the growth rate of asset price or the ‘drift’.

When the market is modeled with a standard Brownian Motion, the probability distribution function of the future price is a normal distribution.

Geometric Brownian Motion

weiner notation

with dS being the change in asset price in continuous time dt. dX is the random variable from the normal distribution (N(0, 1) or Wiener process). σ is assumed to be constant and represents the price volatility considering the unexpected changes that can result from external effects. μdt together represents the deterministic return within the time interval with μ representing the growth rate of asset price or the ‘drift’.

When the market is modeled with a geometric Brownian Motion, the probability distribution function of the future price is a log-normal distribution.

Properties of a Brownian Motion

  • Continuity: Brownian motion is the continuous time-limit of the discrete time random walk. It thus, has no discontinuities and is non-differential everywhere.
  • Finite: The time increments are scaled with the square root of the times steps such that the Brownian motion is finite and non-zero always.
  • Normality: Brownian motion is normally distributed with zero mean and non-zero standard deviation.
  • Martingale and Markov Property: Martingale property states that the conditional expectation of the future value of a stochastic process depends on the current value, given information about previous events. The Markov property instead focusses on the ‘no memory’ theory that the expected future value of a stochastic process does not depend on any past values except the current value. Brownian motion follows both these properties.

Simulating Random Walks for Stock Prices

In quantitative finance, a random walk can be simulated programmatically through coding languages. This is essential because these simulations can be used to represent potential future prices of assets and securities and work out problems like derivatives pricing and portfolio risk evaluation.

A very popular mathematical technique of doing this is through the Monte Carlo simulations. In option pricing, the Monte Carlo simulation method is used to generate multiple random walks depicting the price movements of the underlying, each with an associated simulated payoff for the option. These payoffs are discounted back to the present value and the average of these discounted values is set as the option price. Similarly, it can be used for pricing other derivatives, but the Monte Carlo simulation method is more commonly used in portfolio and risk management.

For instance, consider Microsoft stock that has a current price of $258.65 with a growth trend of 55.2% and a volatility of 35.92%.

A plot of daily returns represented as a random normal distribution is:

Normal Distribution

The above figure represents the simulated price path according to the Geometric Brownian motion for the Microsoft stock price. Similarly, a plot of 10 such simulations would be like this:

Microsoft GBM Simulations

Thus, we can see that with just 10 simulations, the prices range from $100 to over $600. We can increase the number of simulations to expand the data set for analysis and use the results for derivatives pricing and many other financial applications.

Brownian motion and the efficient market hypothesis

If the market is efficient in the weak sense (as introduced by Fama (1970)), the current price incorporates all information contained in past prices and the best forecast of the future price is the current price. This is the case when the market price is modelled by a Brownian motion.

Related Posts

▶ Jayati WALIA Black-Scholes-Merton option pricing model

▶ Jayati WALIA Plain Vanilla Options

▶ Jayati WALIA Derivatives Market

▶ Saral BINDAL Modeling Asset Prices in Financial Markets: Arithmetic and Geometric Brownian Motions

Useful Resources

Academic articles

Fama E. (1970) Efficient Capital Markets: A Review of Theory and Empirical Work, Journal of Finance, 25, 383-417.

Fama E. (1991) Efficient Capital Markets: II Journal of Finance, 46, 1575-617.

Books

Malkiel B.G. (2020) A Random Walk Down Wall Street: The Time-tested Strategy for Successful Investing, WW Norton & Co.

Code

Python code for graphs and simulations

Brownian Motion

What is the random walk theory?

About the author

The article was written in August 2021 by Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022).