Introduction

The specialty insurance market is a large-scale complex system with many uncertainties, complex business relationships and non-linear dynamics and interactions amongst participants. To quantify the behaviour of this complex system, researchers have turned to various traditional approaches ranging from time-series methods (Boyer et al. 2012; Owadally et al. 2018, 2019; Piotr 2016; Venezian 1985; Weiss & Chung 2004) to differential equation based mathematical models (Cummins 1992; Einav et al. 2010; Norberg 1995). This top-down, aggregate approach of modelling complex systems may fail to capture the interactions which occur at a micro-scale which ultimately lead to large-scale emergent phenomena. Indeed, many researchers attest to this challenge within the specialty insurance modelling literature. For instance, the "Underwriting Cycle", an important phenomena where the market undergoes periods of high profitability, less competition and low profitability, high competition among insurers (in the insurance literature known as the hardening and softening of the market), is seldom accurately captured using the vast number of traditional aforementioned methods. Boyer et al. demonstrates the challenges in modelling this phenomena in their paper (Boyer et al. 2012) where they show time-series methods fail to model this important stylised fact of the specialty insurance market. All of the above points reflect the need for an alternative approach, this is where Individual-Based Modelling (IBM) approaches can help.

Individual-Based Models, which encompass modelling frameworks such as Discrete Event Simulations (DES) to Agent-Based Models (ABMs), can alleviate the above challenges posed by traditional methods utilised in the study of insurance markets. For example, these models allow complex systems to be built from the bottom-up, focusing on interactions at the micro-scale between autonomous "agents" which can represent anything from people, organisations to more granular entities such as molecules. Since the late 90’s, IBM applications have grown in popularity, starting from computational social science to ecological studies, biology and environmental sciences. Given its influence in various research fields, the IBM method continues to advance across disciplines, to more relevant areas of research including the modelling of economic markets such as housing, insurance and energy (Arthur 2006; Baptista et al. 2016; Farmer & Foley 2009; Ge 2017; Groff et al. 2019; Hamill & Gilbert 2015; Heinrich et al. 2021; Heppenstall et al. 2007; Johnson & Groff 2014; McLane et al. 2011; Owadally et al. 2018, 2019; Yun & Moon 2020).

The IBM methodology, can enable researchers to observe the emergent phenomena at the aggregate level which arise from the micro level interactions of autonomous heterogeneous agents. Each type of agent can embed behaviours and actions reflective of real-world concepts. In fact, these powerful IBM methodologies have reached industry practitioners within the wider retail insurance market, evidencing the appeal of such approaches in modelling market dynamics. The Institute and Faculty of Actuaries (IFOA), the professional body and regulator of actuaries in the UK, have gone so far as to start an "Agent-Based Modelling Working Party" where workshops have taken place, and research such as (Zhou 2013) presented. Moreover, the Bank of England has published working papers utilising the approach (Baptista et al. 2016). On the other hand, academics from the Bayes Business School (Owadally et al. 2018, 2019) and The Institute for New Economic Thinking (Farmer & Axtell 2022; Farmer & Foley 2009; Heinrich et al. 2021) have all applied ABMs to the insurance market with significant success.

Given the spotlight on the modelling of insurance markets and a lack thereof in specialty insurance markets, this paper will utilise the findings from the existing ABM insurance literature. We combine these findings with the expert knowledge gained through workshops and interviews with underwriters, actuaries, capital modellers, portfolio analysts, brokers and algorithm engineers at Ki and Brit syndicates (insurers) within Lloyd’s of London (the oldest specialty insurance marketplace). These findings allow us to create a novel, discrete event simulation (DES) to study the emergent properties and drivers of the specialty insurance market at various spatio-temporal resolutions. The model proposed in this paper is novel in several aspects:

  • The model utilises a DES framework open-sourced by Ki insurance (HADES). HADES consists of two main components:
    • Processes - these components are responsible for performing some actions and subsequently emitting events as a response to input events, i.e., Broker process responding to events which occur Daily, Monthly, Yearly.
    • Events - a piece of information that can be altered given its interactions with processes.
  • The model incorporates functionality for different event types, such as catastrophe losses (typically low-frequency, high-severity events) and attritional losses (high-frequency, low-severity events). While previous research endeavours have only sought a single type of loss.
  • Interviews and workshops conducted with experts within Ki and Brit insurance have shaped the conceptualisation of the agents represented in the model.
  • The model simulates the dynamics of the Lloyd’s of London specialty insurance market by incorporating unique features such as lead and follow insurers.

This paper will conduct four experiments with varying complexity to demonstrate the emergent properties of the specialty insurance market. The first experiment, explores the interactions between syndicates and brokers with simplistic pricing methods to demonstrate a profitable market with low volatility. The second experiment incorporates catastrophe events that increase volatility in the market which can lead to insolvencies and pronounced cyclicality. Thirdly, the decision-making of syndicates is advanced with the introduction of improved exposure management, where syndicates become better equipped in dealing with catastrophe losses. Finally, we introduce the syndication of risk in the marketplace, via lead-follow mechanics, which significantly reduces volatility and more tightly couples syndicates’ loss experiences.

Some initial findings from the model have demonstrated stylised facts observed in the specialty insurance market, Lloyd’s of London, thus providing some validation for the model. For example, when insurers price, based on past loss history with uniform risks and attritional losses only, we find that premium prices tend towards the fair analytical price but with substantial variance. When catastrophe events occur which affect multiple risks simultaneously, this leads to large industry losses, which cause steep drops in syndicate capitals and increases in premium prices immediately afterwards. When syndicates are able to use exposure management processes to manage tail risks, this leads to smaller and better placed syndicate portfolios. Lastly, syndication enhances the accuracy of actuarial price estimates and narrows the divergence among syndicates.

The following sections of the paper are Literature Review, where pre-existing insurance based literature is disseminated and strengths/weaknesses highlighted. Model Description section describes every process, event and underlying functionality of the model. Results & Discussion introduces the specialty insurance market in the real-world, then delves into the aforementioned four experiment scenarios and discusses the observed quantitative and qualitative results. Lastly, the Conclusion highlights the findings from the model, how it contributes to specialty insurance research and future avenues to be explored.

Literature Review

The specialty insurance market differs from the conventional retail insurance market. The main difference is that the risks are intrinsically complex, offering cover for perils such as kidnap and ransom, cybersecurity breaches, terrorist attacks and political violence, commercial property damage and personal accident. These risks typically require specialist assessment by insurers and/or researchers De Mot & Faure (2014). Simplifying the process extensively, a client reaches out to a broker with a risk, the broker brings the risk to the Lloyd’s of London marketplace. Underwriters that operate within the Lloyd’s market are reached out to, the underwriters utilise their subject matter expertise, actuarial pricing strategies and portfolio management approaches to decide if they should offer a lead, follow line or no line at all. If an underwriter offers a lead line, they underwrite a portion of the risk and tender a policy as a lead insurer, usually the lead covers a bigger portion of the risk and thus is paid a bigger share of the premiums, subsequently taking on more risk. Every risk is insured by a single lead insurer who usually shapes the policy and multiple follow insurers who agree to the terms and offer a line size. Once the risk has been covered proportionally across many syndicates, premiums are distributed depending on the line size agreed by all parties. Conversely, if a claim comes through, then everyone insurer pays out proportionally, see Thoyts (2010), James (2007).

The hardening and softening of the insurance market has been an area of academic enquiry over the years. The specialty insurance market can be a very costly, dynamic and highly volatile market to operate within, one characteristic of the market is the underwriting cycle. In short, when capacity is high, rates decrease and profits decrease, conversely, losses occur, insurers go insolvent or withdraw and capacity decreases, with subsequent upwards pressure on premium levels to generate acceptable return on investment. This is a typical soft to hard market transition, respectively (Boyer et al. 2012), which attracts new capacity. The cycle then repeats.

In the literature, insurance market phenomena such as the underwriting cycle, have been investigated through a variety of approaches ranging from agent-based to time-series models in order to understand the emergence and drivers of these phenomena (Boyer et al. 2012; Harrington & Danzon 1994; Heinrich et al. 2021; Owadally et al. 2018, 2019; Venezian 1985; Weiss & Chung 2004). Surprisingly, there seems to be contradictory views with regard to the causes of market cycles, for example, Cummins & Danzon (1997) argues that institutional lags are the main reason behind the liability insurance market cycle in the mid-1980s in North America when analysing insurance company level data. However, when analysing market-level data, they claim the opposite. Most in-depth literature reviews (Boyer et al. 2012; Boyer & Owadally 2015; Harrington & Danzon 1994; Weiss & Chung 2004) agree that research into underwriting cycles intensified in the mid-1980s due to the "liability crisis". During this time, an interesting discovery made by Venezian (1985) was that different lines of insurance have quasi-cyclical behaviour but others can not be discerned from white noise fluctuations, i.e., random walk. Due to these discoveries, Venezian et al. proposed auto-regressive (AR) models which replicated these cyclical behaviours. Following the findings from (Venezian 1985), most researchers primarily utilised AR models to investigate cyclicality in the insurance market.

There are several hypothesised causes of underwriting cycles discerned by researchers, these are:

One of the prevailing reasons for why the analysis of market cycles has been inadequate is due to the cycle periods being longer than periods over which data is sampled (Boyer et al. 2012; Boyer & Owadally 2015). Furthermore, traditional time-series approaches work well in analysing deterministic systems, however, given the volatile nature of financial systems, this approach can be limited (Ding & Granger 1996). Moreover, if a modelled system can be affected by external factors then these factors should explicitly be included in the time-series model (Hyndman et al. 2013). For this and many other aforementioned reasons, this paper is well-timed and necessary.

Given the interest in modelling the insurance market, researchers have adopted individual-based modelling approaches, to capture new insights and overcome some of the earlier described issues (Dubbelboer et al. 2017; England et al. 2022; Heinrich et al. 2021; Owadally et al. 2018, 2019; Zhou 2013). Some researchers have utilised individual-based models to discover stylised facts (Giardina & Bouchaud 2003) such as the emergence of market crashes and bubbles. Owadally et al. (2019), Boyer et al. (2012) argue that traditional methods such as AR models, fall short in quantifying the existence of a cycle in time-series data and questions the ability for these methods to forecast accurately. The paper by Owadally et al. (2019) suggests that linear models are insufficient and other approaches used to detect periodicity should be considered. Individual-based modelling approaches are more appealing as they can model the individual behaviours and social interactions at the micro level and generate complex aggregate behaviours such as cycles at the macro level (Farmer & Axtell 2022; Owadally et al. 2018). Some researchers agree that the insurance market is complex, heterogeneous and hard to represent with simple rational models, evidenced by the interactions of individuals such as actuaries, underwriters, claims adjusters, brokers and organisations such as syndicates, regulators, managing agents and capital providers (Feldblum 2001; Fitzpatrick et al. 2004; Ingram & Bush 2013; Ingram & Underwood 2010). The model developed by 1 tests the hypothesis that agents with simple decision-making rules interacting at the individual-level can produce aggregate complex behaviours depicting empirically valid market dynamics. 1 also tests the theory of plural rationality (Ingram & Bush 2013; Ingram & Underwood 2010), where they can parameterise the insurers mark-up calculation where a policy is either priced aggressively or priced using the actuarial fair price (market value). They found that as the pricing aggression increased, market volatility also increased due to large price fluctuations.

Research conducted in Heinrich et al. (2021) proposes an agent-based simulation of the speciality insurance market. An important aim of the article is to measure the effect of bankruptcy under different risk model regiments i.e., all syndicates using the same or a diverse range of risk models. The model proposed agents and processes which represented re-insurers, insurers, customers, shareholders and catastrophe bonds that interact with each other and share information. The paper finds that when the number of risk models increases for an insurer, i.e., an insurer utilises several catastrophe models (these are processes used to evaluate and manage natural and man-made catastrophe risk from perils, for example, hurricanes and wildfires) more risks are insured and conversely fewer risk models lead to lower profitability and higher risk of defaulting. The paper provides many interesting results with regards to catastrophe modelling approaches, to advanced exposure management features, whereby each insurer follows real-world regulations such as Solvency II and the "Solvency Capital Requirement". This requirement is described as Insurers are required to have 99.5% confidence they could cope with the worst expected losses over a year. That is, they should be able to survive any year-long interval of catastrophes with a recurrence frequency equal to or less than 200 years. However, in the model of Heinrich et al. (2021) some important aspects of the specialty insurance market have been neglected. These include attritional losses (England 2002; Owadally et al. 2018) as well as the syndication of risk, a unique feature of Lloyd’s of London (Herschaft 2004), and finally the temporal granularity is fixed. We understand that these features may not have fallen within the scope of the article, but believe them to be important to modelling the specialty insurance market. Our proposed solution addresses these limitations. We believe our model is the first discrete event simulation of the specialty insurance market inspired by real-world actors within Lloyd’s of London.

A key insight and novelty of our paper is in the use of a novel DES framework as opposed to the traditional ABM approach adopted by the aforementioned papers (Dubbelboer et al. 2017; England et al. 2022; Farmer & Foley 2009; Heinrich et al. 2021; Owadally et al. 2018, 2019; Zhou 2013). Unlike the classical ABM approach, the DES approach captures the time-irregular and asynchronous nature of events in the specialty insurance market whereby events (such as insurance claims) occur at irregular time intervals from as little to as many events occurring simultaneously. Handling such unpredictable, irregular and asynchronous events in a classic ABM approach can be challenging, but is seamlessly handled using the DES approach. The benefits of DES have thoroughly been discussed in the following papers (Hill et al. 2001; Jacobson et al. 2006; Nevins et al. 1998). Our model architecture is highly modular compared to the other approaches (Farmer & Foley 2009; Heinrich et al. 2021; Owadally et al. 2018, 2019), as evidenced in the workflow diagram in Figure 1, this allows the user to "plug-in" with ease different pricing models (actuarial pricing, underwriter markup, exposure management), relationship networks (the relationships between syndicates and brokers), loss generators (attritional, catastrophe or both) and many output metrics such as yearly, monthly syndicate capitals, premiums offered and accepted, industry statistics (market health indicators), syndicate performance (syndicate health indicators), syndicate insolvency, catastrophe event (the risks affected and extent of losses), claims per risk and syndicates assigned to risk. These configuration processes and metrics will be described in detail in the Model Description and Results & Discussion sections.

There are many reasons why insurers and/or researchers would utilise our proposed model, or an evolution of this model, to simulate the specialty insurance market dynamics. Practitioners may want a deeper understanding of the market dynamics, such as underwriting cycles, which drive the stability and profitability of the market. Moreover, the model provides both a quantitative and qualitative picture of the impact and drivers of these dynamics. Lastly, the model could be used to develop heuristics which enable better estimation of the current and future market conditions.

Model Description

The purpose of the Lloyd’s of London model is to investigate the emergent, complex characteristics which arise from individual-level interactions of the specialty insurance market. In this pursuit, the model conceptualises the behaviours of syndicates and brokers with which complex processes are undertaken as shown in Figure 1. These individual-level interactions from the bottom-up lead to the emergence of stylised facts at the aggregate level which correspond to empirical market trends. In this section, we describe the main agents/processes which the model encapsulates, the events they respond to and emit (i.e. the actions they undertake) and how these features all fit together to produce the DES. We start this section by discussing the notion of time in the model and its importance.

Time

Discrete Event Simulations do not necessarily require the notion of time compared to other methods such as ABMs. Instead, events trigger processes which the process in turn responds to by emitting additional events and/or performing actions, i.e., updating an internal state variable. In theory, there is no need for the concept of time. However, in reality, there is a notion of time associated with events. For simulations of complex financial systems, such as the specialty insurance market, the challenge is that events occur asynchronously while the timescale is irregular and varying. Capturing these features of the specialty insurance market is both paramount (in order to ensure the model reflects empirical trends) and difficult to achieve with a typical ABM as mentioned in the following (Heinrich et al. 2021; Owadally et al. 2019; Zhou 2013). Conversely, the flexibility of the DES framework, allows events to have an associated timestep, which only triggers once the timestep in question occurs, e.g., \(t = 3\). In our DES, a regular timer is incorporated with a standard timestep of days, months and years. This allows us to capture events which occur in the simulation at different levels of granularity. Moreover, we are able to investigate different emergent phenomena independent of the timescale across which they occur.

The regular timer in our model generates the events shown in Table 1. We also detail the processes which respond to the events and the description of the events.

Table 1: Events which the time process generates.
Event Processes which respond to event Description
Day Broker Process The smallest unit of time in the model. Trigger processes that must respond on a daily basis.
Month Syndicate Process, Actuarial Sub-Process, Central Risk Repository Process Trigger processes that must respond on a monthly basis.
Year Syndicate Process, Actuarial Sub-Process, Underwriting Sub-Process, Premium Exposure Management Sub-Process, Value at Risk (VaR) Exposure Management Sub-Process Trigger processes that must respond on a yearly basis.

Broker process

The primary purpose of the Broker process is to bring new risks to the market. This occurs in response to a Day event as noted in Table 1. Upon the triggering of a Day event, the broker generates \(n_{r}\) risks according to a Poisson distribution where the \(\lambda\) variable for the distribution is given by the risks per day (RPD) input parameter as described in Table 13. The Broker process emits a number of events when the risk is generated in order to broadcast the risk to the insurers/syndicates and to set deadlines for the quotes from the insurers to be finalised.

The broker responds to and generates the events shown in Table 2.

Table 2: Events which the broker process responds to and generates.
Event Processes which respond to event Processes which generate event Description
Day Broker Process Time Process Triggers the broker process to generate new risks.
Risk Broadcasted Broker-Syndicate Network Process, Central Risk Repository Process, Attritional Loss Generator Process, Industry Statistics Process Broker Process Broadcasts a new risk has entered the market. The responding processes perform certain actions in response.
Lead Quote Consolidation Deadline Reached Syndicate Process Broker Process Sets a deadline for any lead syndicates who plan on sending an offer to consolidate and send the quote offer.
Lead Quote Selection Deadline Reached Central Risk Repository Process Broker Process Sets a deadline for a lead syndicate to be chosen based on the minimum quote offered by the syndicates.
Follow Quote Consolidation Deadline Reached Syndicate Process Broker Process Sets a deadline for any follow syndicates who plan on sending an offer to consolidate and send the quote offer
Follow Quote Selection Deadline Reached Central Risk Repository Process Broker Process Sets a deadline for follow syndicates to be chosen based on the line sizes offered by the syndicates.

Broker-syndicate network process

The broker-syndicate network process is responsible for requesting quotes for new risks entering the market from the registered syndicates. This is facilitated by a network topology which for a given risk, selects a number of syndicates to which a quote is requested. Therefore, a risk does not necessarily have to be broadcasted to all syndicates depending on the parameter choices. The options for the topologies are the following:

  • Circular topology: inspired by the model presented in Owadally et al. (2018)
  • Network topology: inspired by the interviews and workshops with stakeholders at Ki & Brit Syndicates.
  • Random topology: is a base case feature.

Prevalent in all the topologies is a lead and follow \(top\_k\), parameter as mentioned in Table 13. This parameter selects the best \(k\) syndicates based on the chosen topology methods. For instance, in the circular topology (Owadally et al. 2018), the distance between the brokers and syndicates is used as a measure of the "ease of doing business" where a small distance implies a low cost of doing business conversely a large distance implies a larger cost. The syndicates are then ordered based on the \(top\_k\) parameter which selects the \(k\) lowest cost syndicates. The network topology represents a connected network/graph between the brokers and syndicates where the edge weightings represent the ease of doing business between the brokers and syndicates. The larger the edge weighting, the more likely for brokers and syndicates to do business and vice-versa. Once again, the syndicates are ordered based on the strongest relationship with the \(top\_k\) parameter filtering this down to the strongest \(k\) relationships. Lastly, in the random topology (the adopted network for the experiments conducted below, inspired by the Erdős–Rényi model), as the name suggests, the syndicates are randomly ordered with the \(top\_k\) parameter selecting the first \(k\) syndicates in the list.

Given the above, the broker-syndicate network responds to and generates the events shown in Table 3.

Table 3: Events which the broker-syndicate network process responds to and generates.
Event Processes which respond to event Processes which generate event Description
Risk Broadcasted Broker-Syndicate Network Process Broker Process Triggers broker-syndicate network to select \(k\) syndicates to request a Leader Quote and \(k\) syndicates for a Follower Quote.
Lead Quote Requested Syndicate Process Broker-Syndicate Network Process Triggers syndicate process to calculate a price and lead line size for the risk if possible.
Follow Quote Requested Syndicate Process Broker-Syndicate Network Process Triggers syndicate process to calculate a follow line size based on leaders terms if possible.

Central risk repository process

The central risk repository process is responsible for tracking all of the risks, quotes and underwritten policies in the market. This includes all of the quotes offered for a risk, the policy information if the risk has been underwritten such as who the leader and follower syndicates are. Finally, the central risk repository is responsible for applying any losses, whether attritional or catastrophic, to the underlying syndicates on the policy in question.

Given the above, the central risk repository responds to and generates the events shown in Table 4.

Table 4: Events which the central risk repository process responds to and generates.
Event Processes which respond to event Processes which generate event Description
Month Central Risk Repository Process Time Process Triggers central risk repository process to calculate monthly statistics about the risks, quotes and policies.
Risk Broadcasted Central Risk Repository Process Broker Process Registers new risk entering the market.
Lead Quote Offered Central Risk Repository Process Syndicate Process Registers a lead quote offered by a syndicate for a given risk.
Lead Quote Selection Deadline Reached Central Risk Repository Process Broker Process Triggers the deadline by which no more lead quote offers are accepted. The central risk repository process selects the leader based on the cheapest quote.
Follow Quote Offered Central Risk Repository Process Syndicate Process Registers a follow quote offered by a syndicate for a given risk.
Follow Quote Selection Deadline Reached Central Risk Repository Process Broker Process Triggers the deadline by which no more follow quote offers are accepted. The central risk repository checks to ensure that a leader has already been selected, how much line size remains and assigns the requested line size to the followers as long as it does not exceed 100% whereby they are proportionally signed down.
Attritional Loss Occurred Central Risk Repository Process Attritional Loss Generator Process Triggers the attritional loss to be broadcasted to the syndicates on the policy.
Catastrophe Loss Occurred Central Risk Repository Process Catastrophe Loss Generator Process Triggers the catastrophe loss to be broadcasted to the syndicates on the policy.
Lead Quote Accepted Syndicate Process Central Risk Repository Process Notifies syndicate that their lead quote was accepted.
Follow Quote Accepted Syndicate Process Central Risk Repository Process Notifies syndicate that their follow quote was accepted.
Claim Received Syndicate Process Central Risk Repository Process Notifies syndicate that a claim was received for a risk they underwrite.
Industry Pricing Statistics Reported Industry Statistics Process Central Risk Repository Process Passes industry statistics for the risks, quotes and policies in the central risk repository.

Syndicate process

The syndicate process represents one of the most detailed agents within the model, as it is responsible for a number of important functions. Primarily, the syndicate is responsible for pricing any risks which come to market and decide which line size to give. This is all done in the context of a capital management framework whereby syndicates must ensure they are appropriately capitalised even in the case of tail loss events occurring, in order to avoid insolvency. Lastly, the syndicate must also provide any dividend back to capital providers in the case of profitable performance. As the syndicate is responsible for a number of functions, this section has been split into a number of sub-sections according to the sub-processes which compose the main syndicate process.

Before moving on to the sub-processes, we note that the main syndicate process is responsible for coordinating and organising the sub-processes which compose it. In this regard, the syndicate process responds to and generates the events shown in Table 5.

Table 5: Events which the syndicate process responds to and generates.
Event Processes which respond to event Processes which generate event Description
Month Syndicate Process Time Process Triggers the syndicate process to share state variables, i.e., capital with its sub-processes.
Year Syndicate Process Time Process Triggers the dividend sub-module to evaluate the profitability of the syndicate and pay a dividend if profitable.
Lead Quote Requested Syndicate Process Broker-Syndicate Network Process Triggers the syndicate process to note that a lead quote has been requested so that the actuarial, underwriting, line size and exposure management sub-processes can prepare or refuse a quote.
Follow Quote Requested Syndicate Process Broker-Syndicate Network Process Triggers the syndicate process to note that a follow quote has been requested so that the actuarial, underwriting, line size and exposure management sub-processes can prepare or refuse a quote.
Quote Component Computed Syndicate Process Actuarial Sub-Process, Underwriting Sub-Process, Premium Exposure Management Sub-Process, VaR Exposure Management Sub-Process Triggers the syndicate process to take the different components of the quote and store all of these components until a quote consolidation deadline is reached.
Lead Quote Consolidation Deadline Reached Syndicate Process Broker Process Triggers the syndicate to take all the stored quote components of a lead quote and to either emit a LeadQuoteOffered event or if one of the components rejects the quote (i.e., exposure management deems it risky) then nothing is emitted.
Lead Quote Accepted Syndicate Process Central Risk Repository Process Notifies the syndicate that their lead quote has been accepted by the broker and they now underwrite the policy (i.e., responsible for claims and collect premiums).
Follow Quote Consolidation Deadline Reached Syndicate Process Broker Process Triggers the syndicate to take all the stored quote components of a follow quote and to either emit a FollowQuoteOffered event or if one of the components rejects the quote (i.e., exposure management deems it risky) then nothing is emitted.
Follow Quote Accepted Syndicate Process Central Risk Repository Process Notifies the syndicate that their follow quote has been accepted by the broker and they now underwrite the policy (i.e. responsible for claims and collect premiums).
Claim Received Syndicate Process Central Risk Repository Process Notifies the syndicate that a claim has been recieved for one of the policies they underwrite. The claim is paid and in case of insolvency, the syndicate emits a SyndicateBankrupted event.
Lead Quote Offered Central Risk Repository Process Syndicate Process Triggers the central risk repository to note that a lead quote is offered by the syndicate for the risk in question.
Follow Quote Offered Central Risk Repository Process Syndicate Process Triggers the central risk repository to note that a follow quote is offered by the syndicate for the risk in question.
Syndicate Bankrupted Actuarial Sub-Process, Underwriting Sub-Process, Premium Exposure Management Sub-Process VaR Exposure Management Sub-Process, Line Size Sub-Module, Dividend Sub-Module Syndicate Process Triggers the DES framework to remove the syndicate and all of its sub-processes from the simulation.
Syndicate Capital Reported Actuarial Sub-Process, Underwriting Sub-Process, Premium Exposure Management Sub-Process VaR Exposure Management Sub-Process, Line Size Sub-Module, Dividend Sub-Module Syndicate Process Updates the shared state of capital among the syndicate sub-processes.
Actuarial sub-process

Pricing of risks in the insurance market is a complex process; risks are heterogeneous, extremely variable, and fundamentally stochastic and ambiguous, which are commonly employed assumptions in various papers (Barrett 1999; Dionne 2000; Heinrich et al. 2021; Owadally et al. 2019). Our model captures all of these features but makes the simplifying assumption that the type of risks are homogeneous and that they belong to a catastrophe-exposed class such as property insurance. The price of a risk will depend on several factors, primarily, the nature of the risk itself (i.e., the frequency and severity with which the risk can occur) and the market’s view of the risk (i.e., the supply and demand for insurance services) (Kunreuther 1989; Owadally et al. 2019; Zhou 2013).

These two factors differ from one another. The first is related to the quantification of the risk itself without considering any external market influences. This process is often carried out by an actuary (Owadally et al. 2018; Zhou 2013). The actuary will assess the risk based on experience-based metrics (historical data on the risk or similar risks) and/or exposure-based metrics (quantification of the risk in the absence of data but based on a risk profile). Fundamentally, the actuary attempts to come up with a "fair price" for the risk based on the expectation of losses. That is to say, a price which would cover the expected losses of the risk over its lifetime. The second factor considers the price of the risk based on the market’s view, i.e., the supply and demand for insurance services. This process is typically carried out by an underwriter. The underwriter will take the fair price guideline from the actuary and, in very simplistic terms, scale this price up or down based on the supply and demand in the market. For instance, although the actuary may propose a much higher fair price, market trends may force the underwriter to reduce this price in order to remain competitive.

In this section, we will focus on the actuarial sub-process, inspired by the work of Owadally et al. (2019), Zhou (2013). In the next section, we will discuss the underwriting sub-process.

Based on the above, the actuarial sub-process price is given by two main components. The first component is the insurer’s expected claim cost:

\[ \tilde{P}_t = z\bar{X}_t + (1-z)\lambda'_t\mu'_t\] \[(1)\]
where \(\tilde{P}_t\) is the insurers expected claim cost at timestep \(t\), \(\bar{X}_t\) is the insurers past weighted average claims, \(\lambda'_t\) is the industry-wide average claim frequency, \(\mu'_t\) is the industry-wide expected claim cost and finally \(z\) is the internal experience weight input parameter Table 13 which decides whether a syndicate weighs their own loss or the industry loss experience as more important. As per Owadally et al. (2019), Zhou (2013), the insurers expected claim cost, \(\bar{X}_t\), is calculated as a simple exponentially weighted moving average where the weight is an input parameter to the model Table 13 called the loss experience recency weight.

The final actuarial price, which we denote as \(P_{at}\), is the sum of the insurers expected claim cost, \(\tilde{P}_t\), and a "risk loading term", \(\alpha F_t\), where \(F_t\), is the standard deviation of the insurer’s claims while \(\alpha\) is a input parameter to the model Table 13, called the volatility weight:

\[ P_{at} = \tilde{P}_t + \alpha F_t\] \[(2)\]

As can be seen, the actuarial price in Equation 2, captures the main idea of syndicates pricing a risk to cover their expected claim losses as well as to allow for some volatility in the losses.

The actuarial sub-process responds to and generates the events shown in Table 6.

Table 6: Events which the actuarial sub-process responds to and generate.
Event Processes which respond to event Processes which generate event Description
Quote Requested Actuarial Sub-Process Syndicate Process Triggers actuarial sub-process to calculate the actuarial price as per Equations 1 and 2.
Industry Loss Statistics Reported Actuarial Sub-Process Industry Statistics Process Distributes the latest industry statistics so that the actuarial sub-process can use this in Equation 1.
Industry Pricing Statistics Reported Actuarial Sub-Process Industry Statistics Process Distributes the latest industry pricing statistics so that the actuarial sub-process can use this in Equation 1.
Claim Received Actuarial Sub-Process Central Risk Repository Process Triggers the actuarial sub-process to update the syndicate’s claim experience so it can be used in Equation 1.
Syndicate Bankrupted Actuarial Sub-Process Syndicate Process Triggers the actuarial sub-process to be removed from the simulation after its associated syndicate becomes insolvent.
Quote Component Computed Syndicate Process Actuarial Sub-Process Sends the associated syndicate the actuarial pricing component i.e., Equation 2.
Underwriting sub-process

As explained in the previous section, in this simple model the objective of the underwriting sub-process is to "scale" the actuarial price in order to match market supply and demand. We again, employ the equations used by Owadally et al. (2019), Zhou (2013), where they apply neoclassic price theory for the price-elasticity of demand. The details of the derivation can be found in Owadally et al. (2019), Zhou (2013). The underwriter scaling is given by:

\[ P_{t} = P_{at}e^{m_t}\] \[(3)\]
where \(P_{t}\) is the final price offered by the syndicate after the underwriters scaling, \(P_{at}\) is the actuarial price from Equation 2 and \(m_t\) is the underwriter log markup which attempts to model the price-elasticity of demand in the market.

The underwriter markup, \(m_t\), is also calculated as a simple exponentially weighted moving average where the weight is an input parameter to the model Table 13 called the underwriter markup recency weighting. Details on how the underwriter markup is calculated can be found in Owadally et al. (2019), Zhou (2013).

To illustrate the definition of \(m_t\) within a discrete event model, it is imperative to perceive \(m_t\) not simply as a fixed parameter but as a dynamic entity that evolves in response to market fluctuations. Conceptually, \(m_t\) is the underwriter’s adjustment to the complex balance of market supply, demand, and competitive pressures. This dynamic is essential for simulating the cyclical nature observed in insurance markets, particularly the influence of competition on premium adjustments.

In a discrete event framework, \(m_t\) is re-calibrated at each event occurrence, for instance, upon a quote request or market update, employing the following conceptual formula:

\[m_{t+1} = \alpha m_t + (1 - \alpha)\left(\ln\left(\frac{Demand_{t}}{Supply_{t}}\right) + Competitive\_Factor\right)\] \[(4)\]

Here, \(\alpha\) signifies the underwriter markup recency weighting, reflecting the balance between recent and historical market conditions. \(Demand_{t}\) and \(Supply_{t}\) represent the market demand and supply at time \(t\), respectively, whilst \(Competitive\_Factor\) is an adjustment reflecting the competitive dynamics within the market, influencing the underwriting pricing strategy.

This conceptualisation underscores the adaptability of \(m_t\) to market dynamics, facilitating the simulation of cycles induced by competitive premium adjustments. The inclusion of a competitive factor emphasises the model’s capacity to replicate market behaviours arising solely from competitive interactions, consistent with the cycles observed in empirical studies.

The underwriting sub-process responds to and generates the events shown in Table 7.

Table 7: Events which the underwriting sub-process responds to and generate.
Event Processes which respond event Processes which generate event Description
Quote Requested Syndicate Process Underwriting Sub-Process Triggers underwriting sub-process to calculate the final price as per Equations 3.
Syndicate Bankrupted Syndicate Process Underwriting Sub-Process Triggers the underwriting sub-process to be removed from the simulation after its associated syndicate becomes insolvent.
Quote Component Computed Underwriting Sub-Process Syndicate Process Sends the associated syndicate the underwriting pricing component i.e., Equation 3.
Value at Risk (VaR) exposure management sub-process

Exposure management is a critical function for all insurance firms, to such that appropriate levels of exposure management are regulated by law for insurers. Exposure management allows insurers to have a quantitative understanding of the tail risk associated with the policies that they underwrite. By doing so, insurers can quantify the impact of worst case scenarios on their portfolio. This informs insurer’s underwriting strategy. For instance, given the state in Figure 2, an over-exposed strategy would be for lead syndicate B to underwrite all risks within region 1. Therefore, exposure management can alert syndicates if they are over-exposed to a given catastrophe peril region.

The Value at Risk (VaR) of an insurer’s portfolio is a common measure used to quantify the level of risk taken by an insurer. The VaR with some exceedance probability, which we denote as \(\alpha\) Table 14, identifies the amount of syndicate capital which is at risk in case of any tail events occurring e.g. exceedingly large catastrophe events. The exposure management sub-process is therefore responsible for ensuring that the syndicate capital remains above this threshold value. Given its importance, our model employs the VaR exposure management detailed in Heinrich et al. (2021).

The VaR Exposure Management sub-process responds to and generates the events shown in Table 8.

Table 8: Events which the VaR EM sub-process responds to and generates.
Event Processes which respond to event Processes which generate event Description
Simulation Started VaR Exposure Management Sub-Process Predefined Event\(^1\) Triggers the VaR EM sub-process to perform the VaR simulations for the various peril regions.
Quote Requested VaR Exposure Management Sub-Process Syndicate Process Triggers the VaR EM sub-process to check its exposure management to the peril region in question. If the risk can not be underwritten then the QuoteComponentComputed event reflects this.
Syndicate Capital Reported VaR Exposure Management Sub-Process Syndicate Process Updates the common state of capital between the syndicate process and the EM sub-process.
Syndicate Bankrupted VaR Exposure Management Sub-Process Syndicate Process Triggers the VaR EM sub-process to be removed from the simulation after it’s associated syndicate becomes insolvent.
Quote Component Computed Syndicate Process VaR Exposure Management Sub-Process Sends the associated syndicate the VaR EM quote component. This either allows the quote to go ahead or vetoes the quote in case the capital requirements of the syndicate do not satisfy the minimum exposure management standards.
\(^{1}\) A special event initialised by the DES framework
Premium exposure management sub-process

In the previous section, one approach to exposure management was explored via the VaR measure. However, as extensively discussed by Heinrich et al. (2021), VaR exposure management is a complicated process, often relying on computationally expensive monte carlo simulations and in many cases, difficult-to-measure variables. Often insurers seek approximations or proxy approaches which capture the essence of the VaR exposure management. In this section, we detail such a methodology which we refer to as Premium Exposure Management.

The premium an insurer collects for underwriting a risk can be thought of as a proxy for exposure. In particular, when the total premium written is compared to the total capital available, this gives a proxy measure of how exposed an insurer is to the potential risk of insolvency. For example, if an insurer has underwritten a large number of risks, and therefore is collecting a large premium, but their total capital is comparatively smaller. This indicates that the insurer might have insufficient capital available to cover the full range of potential losses. On the other hand, if premiums written were large but the capital available was also large then this would indicate that the insurer has suffered minimal losses and that their current underwriting strategy was profitable. This is the essence of the premium exposure management sub-process.

The Premium Exposure Management sub-process responds to and generates the events shown in Table 9.

Table 9: Events which the Premium EM sub-process responds to and generates.
Event Processes which respond to event Processes which generate event Description
Quote Requested Premium Exposure Management Sub-Process Syndicate Process Triggers the Premium EM sub-process to calculate the premium EM scaling factor which is used to scale the quote price based on the current exposure.
Syndicate Capital Reported Premium Exposure Management Sub-Process Syndicate Process Updates the common state of capital between the syndicate process and the EM sub-process.
Syndicate Bankrupted Premium Exposure Management Sub-Process Syndicate Process Triggers the Premium EM sub-process to be removed from the simulation after it’s associated syndicate becomes insolvent.
Quote Component Computed Syndicate Process Premium Exposure Management Sub-Process Sends the associated syndicate the Premium EM quote component. This component is the scaling factor applied to the quote price based on the capital requirements of the syndicate.
Line size sub-module

A important feature of the Lloyd’s of London specialty insurance market is the syndication of risk i.e., there will typically be several syndicates on a policy which includes one lead syndicate and several follow syndicates. Typically the lead syndicate negotiates the price, as well as the contractual terms of the policy, with a broker and take on some fraction of the risk. Follow syndicates will accept similar terms to the lead syndicate and also take on a fraction of the risk. The fraction of a risk offered is referred to as a line size. In our simple model, the lead syndicates offer a default line size and set the price of the policy based on the actuarial and underwriting sub-processes. The candidate follow syndicates, also price the risk based on their actuarial and underwriting sub-processes which they then compare to the actual price from the lead syndicate. This allows them to assess the "pricing strength" which in turn is used to decide the line size to offer. The pricing strength is defined as the ratio of the follower’s proposed price to the lead price. That is, if the pricing strength is above one, this implies the price of the risk is good and a larger line size is offered and vice versa.

Dividend sub-module

As explored in Heinrich et al. (2021), our DES model also includes, the ability for syndicates to pay a dividend to capital holders on a yearly basis provided they have made a profit. The reason for why this is a class object and not a DES process, is because the only process that relies on its outputs is the syndicate process.

The dividend sub-module only becomes active when a Year event triggers the syndicate process. When this occurs, the syndicate uses the dividend sub-module to check whether a profit has been made and if so calculates the dividend as follows:

\[ D = \gamma Pr_t\] \[(5)\]
where \(D\) is the dividend paid, \(Pr_t\) is the profit made by the syndicate, and \(\gamma\), profit fraction, is an input parameter to the model Table 13, which represents the fraction of the profit to pay out as a dividend.

Attritional loss generator process

Attritional losses (as opposed to catastrophe losses) are defined as those losses which are generally uncorrelated with each other in both space and time, have high frequency, low severity and are fairly predictable (England 2002). On the other hand, catastrophe losses, which will be discussed in the next sub-section, tend to be spatially correlated (affecting a number of policies concentrated by a given peril region, class or industry), low frequency, high severity and difficult to predict.

In our model, we develop the attritional loss generator process inspired by the work of 1, Owadally et al. (2019), Zhou (2013). The attritional loss generator process pre-generates a number of AttritionalLossOccurred events when the risk is first brought to the market. The number of claim events is given by the Poisson distribution with the \(\lambda\) value set as the yearly claim frequency, an input parameter as defined in Table 13. The severity of the loss is defined by the gamma distribution. The shape parameter of the distribution is given by \(\frac{1}{COV^2}\), where \(COV\) is the gamma distribution’s coefficient of variation as defined in Table 13. The scale parameter of the distribution is given by \(\mu COV^2\) where \(\mu\) is the mean of the gamma distribution. These AttritionalLossOccurred events are then placed within the event queue within the expiration date of the risk.

Given the above, the attritional loss generator responds to and generates the events shown in Table 10.

Table 10: Events which the attritional loss generator process responds to and generates.
Event Processes which respond to event Processes which generate event Description
Risk Broadcasted Attritional Loss Generator Process Broker Process Triggers the attritional loss generator process to pre-generate a number of AttritionalLossOccurred events as specified above.
Attritional Loss Occurred Central Risk Repository Attritional Loss Generator Process Triggers the central risk repository to broadcast the attritional loss to the syndicates on the policy in question.

Catastrophe loss generator process

In our model, we develop the catastrophe loss generator process inspired by Heinrich et al. (2021). The catastrophe loss generator attempts to capture the phenomena that catastrophe losses are both correlated spatially (e.g. a number of policies concentrated by a given peril region) and temporally. For this reason, all risks in the model have an associated peril region, as observed in Figure 2. Unlike the attritional loss generator, which generates attritional events and losses on a risk by risk basis, the catastrophe loss generator generates catastrophe events and losses on a peril region basis. The total loss affecting a given peril region is then cascaded down to the affected risks within the peril region, and subsequently the lead and follow insurers.

When the simulation starts, the catastrophe loss generator pre-generates a number of catastrophe events over the length of the entire simulation. This is done via the Poisson distribution with the \(\lambda\) value set as the product of the mean number of catastrophe events per year Table 13 and the number of years in the simulation. A peril region is randomly assigned to the CatastropheLossOccurred event. The total loss affecting the peril region is given by a truncated Pareto distribution with the minimum value set as the minimum catastrophe damage Table 13. The CatastropheLossOccurred events are then added within the event queue up until the end of the simulation.

Given the above, the catastrophe loss generator responds to and generates the events shown in Table 11.

Table 11: Events which the catastrophe loss generator process responds to and generates.
Event Processes which respond to event Processes which generate event Description
SimulationStarted Catastrophe Loss Generator Process Pre-defined Event Triggers the catastrophe loss generator process to pre-generate CatastropheLossOccurred events for the different peril regions as described above and shown in Figure 2..
Catastrophe Loss Occurred Central Risk Repository Catastrophe Loss Generator Process Triggers the central risk repository to cascade the total loss imposed on the peril region to the specific policies within the region.

Industry statistics process

The industry statistics process, as the name suggests, simply keeps a track of all the relevant industry statistics and metrics. This is necessary, as many of the syndicate sub-processes rely on market-wide metrics. Note, that this process does not attempt to mimic any real market agents/processes. It is simply just an aggregator of data.

The industry statistics process responds to and generates the events shown in Table 12.

Table 12: Events which the industry statistics process responds to and generates.
Event Processes which respond to event Processes which generate event Description
Year Industry Statistics Process Time Triggers the industry statistics process to send out yearly market-wide statistics.
Risk Broadcasted Industry Statistics Process Broker Triggers the industry statistics process to log new risks entering the market.
Lead Quote Accepted Industry Statistics Process Central Risk Repository Triggers the industry statistics process to log newly underwritten risk.
Claim Received Industry Statistics Process Central Risk Repository Triggers the industry statistics process to log new claim events.
Industry Loss Statistics Reported Actuarial Sub-Process Industry Statistics Process Sends industry loss statistics to the actuarial sub-process in order to perform relevant calculations.
Table 13: Model input parameter descriptions
Parameter Process Description Default Value
Risks per day (RPD) Broker \(\lambda\) value for the Poisson distribution used by the broker process to generate new risks. 0.06
Number of peril regions Broker Number of peril regions for the catastrophes. 10
Risk limit Broker The maximum value of a risk. $10,000,000
Lead \(top\_k\) Broker-Syndicate Network The number of lead syndicates a broker reaches out to. 2
Follow \(top\_k\) Broker-Syndicate Network The number of follow syndicates a broker reaches out to. 5
Yearly Claim Frequency Attritional Loss Generator \(\lambda\) value for the Poisson distribution for the number of attritional claims generated per year. 0.1
\(COV\) Attritional Loss Generator Coefficient of Variation for the gamma distribution which generates the severity of attritional claim events. 1
\(\mu\) Attritional Loss Generator Mean of the gamma distribution which generates the severity of attritional claim events. $3,000,000
Mean catastrophe events per year Catastrophe Loss Generator \(\lambda\) value for the Poisson distribution for the number of catastrophe claims generated per year. 0.05
Pareto shape Catastrophe Loss Generator The shape parameter of the Pareto distribution which generates the severity of catastrophe claim events. 5
Minimum catastrophe damage Catastrophe Loss Generator The minimum value for an event to be considered a catastrophe (fraction of the risk limit). 0.25
Capital Syndicate Process The initial syndicate capital. $10,000,000
Default lead quote line size Syndicate Process Default value for the lead line size offered. 0.5
Default follow quote line size Syndicate Process Default value for the follow line size offered. 0.1
Internal experience weight Actuarial Sub-Process Weighting which dictates whether the actuarial pricing is based primarily on syndicate history or industry-wide history. 0.5
Loss experience recency weight Actuarial Sub-Process Weighting which dictates whether actuarial pricing weighs the past losses more than recent losses. 0.2
Volatility weight (\(\alpha\)) Actuarial Sub-Process Scaling factor which dictates how much actuarial pricing considers the standard deviation of losses i.e. volatility of losses. 0
Underwriter markup recency weighting Underwriter Sub-Process Weighting which dictates whether the underwriter markup weights past mark-ups more than recent ones. 0.2
Profit fraction (\(\gamma\)) Dividend Sub-Module Fraction of profit vented as dividend. 0.4
VaR EM exceedance probability (\(\alpha\)) VaR EM Sub-Process The tail probability used in the VaR calculations. 0.05
VaR EM safety factor VaR EM Sub-Process Scaling safety factor applied to the VaR value - larger values employ more conservative EM. 1
Premium reserve ratio Premium EM Sub-Process Premium to capital reserve ratio. 0.5
Minimum capital reserving ratio Premium EM Sub-Process Reserved capital to working capital ratio. 1
Maximum scaling factor Premium EM Sub-Process Minimum scaling factor applied to premium based on EM. 1

Results & Discussion

The Lloyd’s of London specialty insurance market is a complex system. Many in academia and the insurance industries have tried to identify the underlying processes that lead to the emergence of phenomena at the aggregate level. One particular phenomena that has been an area of focus, is the underwriting cycle (as described in the Literature review section). As a proof-of-concept we propose several experiments that aim to simulate conditions in the market with varying complexity to reproduce existing industry phenomena. Whereby practitioners and researchers can investigate the underlying features of the market that lead to these trends, better preparing those involved in the market to deal with exogenous shocks.

This section first presents empirical market dynamics, then delves into the four experiments, starting with a base case, introduction of catastrophe events, syndicates adapting to these events and finally the introduction of lead and follow dynamics (described in the previous section).

In Figure 3 the market experiences catastrophe events which exacerbates periods of hardening and softening of the market. This phenomena is conveyed as the underwriting cycle, which as a benchmark could be a phenomena we reproduce using our DES model. In the following sub-section, we describe each model scenario (experiment) and subsequently present the results while discussing the findings in relation to Lloyd’s of London.

Model scenarios

To ensure the model is able to tell a story from simplistic behaviours to more advanced features and outcomes, we start with a base case experiment, this and subsequent experiments are described below (see Model Description Section for descriptions of each model feature):

  1. Base case actuarial pricing, all syndicates use actuarial pricing models, premium exposure management and attritional losses occur (Scenario 1).
  2. Catastrophe event, all syndicates use actuarial pricing, premium exposure management, attritional and catastrophe losses occur (Scenario 2).
  3. Syndicates adapting to catastrophe events, all syndicates use actuarial pricing, VaR exposure management, attritional and catastrophe losses occur (Scenario 3).
  4. Leaders and Followers, syndicates can either be leaders or followers and use actuarial pricing, premium exposure management and only attritional losses occur (Scenario 4).
Table 14: Input parameter values for the four experiment scenarios
Parameters Scenario 1 Scenario 2 Scenario 3 Scenario 4
Simulation end time (years) 50 50 50 50
RPD 0.06 0.06 0.06 0.06
Number of peril regions 10 10 10 10
Risk limit 10,000,000 10,000,000 10,000,000 10,000,000
Lead \(top\_k\) 2 2 2 2
Follow \(top\_k\) - - - 5
Yearly Claim Frequency \(\lambda\) 0.1 0.1 0.1 0.1
\(COV\) 1 1 1 1
\(\mu\) 3,000,000 3,000,000 3,000,000 3,000,000
Mean catastrophe events per year - 0.05 0.05 -
Pareto shape - 5 5 -
Minimum catastrophe damage - 0.25 0.25 -
Capital 10,000,000 10,000,000 10,000,000 10,000,000
Default lead quote line size - - - 0.5
Default follow quote line size - - - 0.1
Internal experience weight 0.5 0.5 0.5 0.5
Loss experience recency weight 0.2 0.2 0.2 0.2
Volatility weight (\(\alpha\)) 0 0 0 0
VaR EM exceedance probability (\(\alpha\)) - - 0.05 -
VaR EM safety factor - - 1 -
Premium reserve ratio 0.5 0.5 - 0.5
Minimum capital reserving ratio 1 1 - 1
Maximum scaling factor 1 1 - 1

Base case actuarial pricing (scenario 1)

This experiment utilises the most basic components of the model, in order to demonstrate the ability for core features to represent important market phenomena. In this experiment, we include five syndicates and twenty-five brokers. This maintains the 1/5 ratio of syndicates to brokers as is the case in the Lloyd’s of London market pocket guide (Neal 2019). As we are only utilising the actuarial pricing model, we expect the premium prices to converge around the fair price, which given the model parameters Table 14 is $300,000. Secondly, we expect the loss ratio of syndicates to fluctuate between periods of profitability and losses, i.e., early signs of cyclicality. Note, due to computational complexity, the model was only run for a limited number of times.

In Figure 4a, the syndicate response over time, is indicative of capital fluctuations in reality, as observed in (Heinrich et al. 2021) where the authors discuss capital profiles over time. For instance, some syndicates, i.e., syndicate 0, are bankrupted as a result of their underwriting strategy, while others ride the boom and bust period better. As we hypothesised, the premiums offered have converged to the actuarial fair price of approximately $300,000 Figure 4b, the reason for this (refer to the Sub-Section 4.12) is because the syndicate’s actuarial pricing, offers a price which attempts to cover its prior losses. Given the current model setup, the mean of the prior losses is $300,000.

For context, the loss ratio when \(\geq 1\) indicates an unprofitable syndicate as their losses are greater than their income, i.e., premiums. Conversely, \(< 1\) leads to profitable syndicates. As mentioned, periods of profitability and losses is a crucial trait of all models that simulate the insurance market as described in 1, Owadally et al. (2019), Zhou (2013), Heinrich et al. (2021). Clearly, our model outputs also demonstrate this behaviour as observed in Figure 5.

World shock events, where catastrophes meet the insurance market (scenario 2)

In this experiment, we repeat the input parameters earlier, however, now we introduce the catastrophe loss generator. When catastrophe events occur, this should lead to major losses and the severity should vary across syndicates depending on their underwriting strategies. We expect to see some syndicates go insolvent, while others may not. Furthermore, insurance insiders and academics claim that the primary cause of cyclicality in the market, is due to unexpected catastrophe events (Boyer et al. 2012; Manikowski & Weiss 2012; Venezian 1985). Our model allows us to test this hypothesis, which we attempt in this section. These experiments should show exaggerated cyclical trends compared to the previous experiment.

The most intriguing result from this experiment, is shown in Figure 6, which demonstrates not only the presence of cyclicality in the premiums offered, but provides an explanation for why this phenomena occurs. Put simply, the cyclicality occurs, in the initial phase, when premiums begin to converge to a fair price, however, catastrophe events which result in large losses, forces syndicates to price premiums higher resulting in an increase in the prices offered. Eventually, once the effect of the catastrophe wears off, the syndicates once again try to converge towards the fair price and the cycle repeats itself.

As mentioned previously, cyclicality is pronounced and in some cases, due to large catastrophe losses, loss ratios go above one Figure 7.

Adapting to market shocks, the utility of VaR exposure management (scenario 3)

The purpose of this experiment, is to specifically showcase the effects of advanced exposure management methods. Due to brevity, we only showcase results indicative of this objective. To this end, syndicates have utilised simplistic premium exposure management. Now, we introduce the VaR exposure management, which allows syndicates to manage their exposure in a more sophisticated manner representing behaviours of real-world syndicates more closely.

VaR exposure management (Heinrich et al. 2021) should enable syndicates to adopt underwriting strategies which avoid over-subscribing to a given peril-region. A good exposure management strategy, is to distribute risk underwritten uniformly across all peril-regions (Heinrich et al. 2021). For this reason, we present the uniform deviation (measures how much the real distribution of risks in the peril regions varies from a perfectly uniform peril-region distribution, where 0 is perfectly uniform). For comparison, we include no exposure management and premium exposure management outputs. As can be seen from Figure 8: From top-left to bottom-right, we observe that as the exposure management becomes more sophisticated and stringent, it is clear to see that the uniform deviation moves closer to zero.

Leaders and followers among syndicates (scenario 4)

The Lloyd’s of London specialty insurance market as described in the Literature review section differs from other insurance markets. The most prominent difference is the syndication of risk, i.e., the ability for multiple syndicates to underwrite a risk either as leaders or followers. Lloyd’s of London claims that the syndication of complex risks across multiple syndicates, allows these unconventional risks and the potential losses to be distributed among several insurers as opposed to one. This should result in less volatility, and lower likelihoods of insolvency as described in the documentation published by Lloyd’s (Lloyd’s London 2019; Neal 2019). However, as far as the authors are aware, these claims have not been quantitatively tested in academic research, until now where our DES model allows us to investigate this particular claim.

As observed in Figure 9a compared to Figure 4 b the volatility of the premium offered is significantly lower, with the premiums tightly converging towards the fair price. Due to the syndication of risks, this allows more syndicates to participate significantly in the market. This means that any losses, are shared among all syndicates, implying their loss experience is similar to each other. As a result, they all offer similar prices. This behaviour can also be observed in Figure 9b where the loss ratios of the syndicates are tightly coupled/correlated, indicating a similar loss experience.

In this scenario, we observe that no insolvencies occurred, compared to the same scenario minus the lead and follow dynamics in Sub-section 5.2. These findings, provide strong quantitative justification for the market structure imposed by the regulator Lloyd’s of London.

In summary, the proof-of-concept DES can reproduce quantitatively and qualitatively market phenomena, such as the underwriting cycle, merits of lead - follow market structures, and the importance of exposure management.

Conclusion

This paper proposes a novel DES of the Lloyd’s of London specialty insurance market. The model captures granular interactions of significant actors within the marketplace, i.e., syndicates and brokers. The model has demonstrated significant results with regards to market dynamics observed in the real-world, this includes signs of cyclical behaviours congruent of the underwriting cycle, which becomes more severe with catastrophe events. Given the unique model architecture, users can swap components such as pricing models, broker-syndicate relationship networks and other features with ease that may not have been possible in past attempts such as 1, Owadally et al. (2019), Zhou (2013), Heinrich et al. (2021). Furthermore, the paper has proposed a conceptualisation of the lead and follow syndicate dynamics which is unique to Lloyd’s of London, the results from these specific experiments have reinforced the innovative market structure, employed by Lloyd’s which reduces volatility in their market.

The proposed DES has many strengths as discussed previously, for example, from a technical perspective, we adopt first-hand knowledge from experts in Lloyd’s of London syndicates Ki and Brit and propose a novel, modular DES with changeable components. From a market dynamics and results perspective, we incorporate the lead and follow mechanics unique to Lloyd’s of London, integrate both attritional and catastrophe loss events, and quantitatively study many important market phenomena. However, several areas of improvement can also be highlighted, for example, given the nature of the complexity which drives the market, it is difficult to abstract and quantify all aspects of the market appropriately. An example is quantifying the broker-syndicate relationship which is the cornerstone of the market. Capturing these individual-level human relationships is difficult with any modelling framework, however, individual-based modelling has allowed us to acquire new insights regarding this relationship. Furthermore, reinsurance is an important function of the insurance market. Reinsurance has a big effect on the capacity a syndicate can underwrite as they manage the tail risk of a portfolio. Modelling the availability and cost of reinsurance policies is therefore a significant driver of market dynamics. Conventional calibration and validation has not been pursued in this iteration of the proof-of-concept model, i.e., global sensitivity analysis. However, the behaviours of the model have been verified by Ki and Brit. In future there are plans to incorporate proprietary Ki and Brit data in order to calibrate the model so that these validated results can be disseminated.

Given the wealth of engineered features of the model, many future avenues can be explored. The underwriting markup (pricing model) (Owadally et al. 2018, 2019) was not utilised in the experiments conducted in this paper, however, comparing the different pricing models, i.e., actuarial and underwriting may lead to new insights with regards to market competition. Additionally, as discussed and utilised by Heinrich et al. (2021), another major component of a syndicate’s activities is to pay out dividends in case of profitable performance. Given that the dividend feature is available in the model as described in the Model description section this can easily be explored as a path for future work.

The simplification to homogeneous risks, while facilitating the initial proof-of-concept, acknowledges the complexity and diversity of risks within Lloyd’s specialty markets. This assumption is a limitation for those looking to apply the model directly to real-world data, where risk heterogeneity plays a crucial role. Future enhancements to the model could include mechanisms for incorporating risk heterogeneity, thus allowing for a more nuanced simulation of the insurance market.

One potential avenue for such enhancement is the integration of the lead/follow process more explicitly into the model’s architecture. This process, a hallmark of Lloyd’s marketplace, helps mitigate the impact of the winner’s curse, particularly prevalent in markets dealing with heterogeneous risks. By allowing insurers to follow the pricing and underwriting lead of others, the model could capture the dynamic interplay of competition and risk assessment accuracy, thereby reducing the adverse effects of underestimating risk. This modification would not only enhance the realism of the model but also provide insights into the strategic interactions within the insurance market. We hope this model provides researchers and specialty insurance practitioners with new insights that enable future R&D projects and provide the means to reduce market volatility and enhance insurance businesses.


Acknowledgements

A special thanks to the Alan Turing Institute’s Internship Network (TIN) in supporting the project. Accenture for their partnership with academia. Ki and Brit syndicates for bridging the relationship between academia and industry, providing the resources to undertake the research project. A special thanks to Reuben Thomas-Davis for developing and providing access to the HADES framework and Sachin Suri for their time devoted to improving the quality of the article.

Supplementary Materials

The Hades framework is made open-source by Ki Insurance and can be found at the following link: https://pypi.org/project/hades-framework/

References

ARTHUR, W. B. (2006). Out-of-equilibrium economics and agent-based modeling. In L. Tesfatsion & K. L. Judd (Eds.), Handbook of Computational Economics (Vol. 2, pp. 1551–1564). Amsterdam: Elsevier.

BAPTISTA, R., Farmer, J. D., Hinterschweiger, M., Low, K., Tang, D., & Uluc, A. (2016). Macroprudential policy in an agent-Based model of the UK housing market. Available at: https://www.bankofengland.co.uk/working-paper/2016/macroprudential-policy-in-an-agent-based-model-of-the-uk-housing-market.

BARRETT, M. I. (1999). Challenges of EDI adoption for electronic trading in the London Insurance Market. European Journal of Information Systems, 8, 1–15.

BOYER, M. M., Jacquier, E., & Van Norden, S. (2012). Are underwriting cycles real and forecastable? Journal of Risk and Insurance, 79(4), 995–1015.

BOYER, M. M., & Owadally, I. (2015). Underwriting apophenia and cryptids: Are cycles statistical figments of our imagination? The Geneva Papers on Risk and Insurance - Issues and Practice, 40, 232–255.

CHOI, S., Hardigree, D., & Thistle, P. D. (2002). The property/liability insurance cycle: A comparison of alternative models. Southern Economic Journal, 68(3), 530–548.

CUMMINS, J. D. (1992). Financial Pricing of Property and Liability Insurance. Berlin Heidelberg: Springer.

CUMMINS, J. D., & Danzon, P. M. (1997). Price, Financial Quality, and Capital Flows in Insurance Markets. Journal of Financial Intermediation, 6(1), 3–38.

DE Mot, J., & Faure, M. G. (2014). Special insurance systems for motor vehicle liability. The Geneva Papers on Risk and Insurance - Issues and Practice, 39, 569–584.

DING, Z., & Granger, C. W. J. (1996). Modeling volatility persistence of speculative returns: A new approach. Journal of Econometrics, 73(1), 185–215.

DIONNE, G. (2000). Handbook of Insurance. Berlin Heidelberg: Springer.

DUBBELBOER, J., Nikolic, I., Jenkins, K., & Hall, J. (2017). An agent-based model of flood risk and insurance. Journal of Artificial Societies and Social Simulation, 20(1), 6.

EINAV, L., Finkelstein, A., & Levin, J. (2010). Beyond testing: Empirical models of insurance markets. Annual Review of Economics, 2, 311–336.

ENGLAND, P. D. (2002). Financial Simulation Models in General Insurance. Available at: https://www.actuaries.org/EVENTS/Seminars/New_Delhi/chapters/page-73to89.pdf.

ENGLAND, R., Owadally, I., & Wright, D. (2022). An Agent-Based Model of Motor Insurance Customer Behaviour in the UK with Word of Mouth. Journal of Artificial Societies and Social Simulation, 25(2), 2.

FARMER, J. D., & Axtell, R. L. (2022). Agent-Based modeling in economics and finance: Past, present, and future. Journal of Economic Literature, Forthcoming.

FARMER, J. D., & Foley, D. (2009). The economy needs agent-based modelling. Nature, 460, 685–686.

FELDBLUM, S. (2001). Underwriting cycles and business strategies.

FITZPATRICK, A. L., Powe, N. R., Cooper, L. S., Ives, D. G., & Robbins, J. A. (2004). Barriers to health care access among the elderly and who perceives them. American Journal of Public Health, 94(10), 1788–1794.

GE, J. (2017). Endogenous rise and collapse of housing price: An agent-based model of the housing market. Computers, Environment and Urban Systems, 62, 182–198.

GIARDINA, I., & Bouchaud, J.-P. (2003). Bubbles, crashes and intermittency in agent based market models. The European Physical Journal B - Condensed Matter and Complex Systems, 31, 421–437.

GRACE, M. F., & Hotchkiss, J. L. (1995). External impacts on the property-liability insurance cycle. Journal of Risk and Insurance, 62(4), 738–754.

GROFF, E. R., Johnson, S. D., & Thornton, A. (2019). State of the art in agent-based modeling of urban crime: An overview. Journal of Quantitative Criminology, 35, 155–193.

HALEY, J. D. (1995). A by-line cointegration analysis of underwriting margins and interest rates in the property-liability insurance industry. Journal of Risk and Insurance, 62(4), 755–763.

HAMILL, L., & Gilbert, N. (2015). Agent-Based Modelling in Economics. Hoboken, NJ: John Wiley & Sons.

HARRINGTON, S. E., & Danzon, P. M. (1994). Price cutting in liability insurance markets. Journal of Business, 67(4), 511–538.

HEINRICH, T., Sabuco, J., & Farmer, J. D. (2021). A simulation of the insurance industry: The problem of risk model homogeneity. Journal of Economic Interaction and Coordination, 17, 535–576.

HARRINGTON, S. E., & Yu, T. (2003). Do property-casualty insurance underwriting margins have unit roots? Journal of Risk and Insurance, 70(4), 715–733.

HEPPENSTALL, A. J., Evans, A. J., & Birkin, M. H. (2007). Genetic algorithm optimisation of an agent-based model for simulating a retail market. Environment and Planning B: Planning and Design, 34(6).

HERSCHAFT, J. A. (2004). Not your average coffee shop: Lloyd’s of London - A twenty-first-century primer on the history, structure, and future of the backbone of marine insurance. Tulane Maritime Law Journal, 29, 169.

HILL, R. R., Miller, J. O., & McIntyre, G. A. (2001). Applications of discrete event simulation modeling to military problems. Proceeding of the 2001 Winter Simulation Conference (Cat. No. 01CH37304).

HYNDMAN, R. J., Booth, H., & Yasmeen, F. (2013). Coherent mortality forecasting: The product-ratio method with functional time series models. Demography, 50(1), 261–283.

INGRAM, D., & Bush, E. (2013). Collective approaches to risk in business: An introduction to Plural Rationality Theory. North American Actuarial Journal, 17(4), 297–305.

INGRAM, D., & Underwood, A. (2010). The human dynamics of the insurance cycle and implications for insurers: An introduction to the theory of plural rationalities. ERM Symposium/SOA Monograph.

JACOBSON, S. H., Hall, S. N., & Swisher, J. R. (2006). Discrete-event simulation of health care systems. In R. W. Hall (Ed.), Patient Flow: Reducing Delay in Healthcare Delivery (pp. 211–252). Berlin Heidelberg: Springer.

JAMES, J. (2007). Lloyd’s and the London insurance market: An overview. Handbook of International Insurance: Between Global Dynamics and Local Contingencies, 903–924.

JOHNSON, S. D., & Groff, E. R. (2014). Strengthening theoretical testing in criminology using agent-based modeling. Journal of Research in Crime and Delinquency, 51(4), 509–525.

KARL, J. B., & Nyce, C. (2019). How cellphone bans affect automobile insurance markets. Journal of Risk and Insurance, 86(3), 567–593.

KUNREUTHER, H. (1989). The role of actuaries and underwriters in insuring ambiguous risks. Risk Analysis, 9(3), 319–328.

LENG, C.-C., & Meier, U. B. (2006). Analysis of multinational underwriting cycles in property-liability insurance. The Journal of Risk Finance, 7(2), 146–159.

LLOYD’S London. (2019). Business at Lloyd’s is still conducted face-to-face, and the bustling underwriting room is central to the smooth running of the market. Available at: https://www.lloyds.com/about-lloyds/our-market/lloyds-market.

MANIKOWSKI, P., & Weiss, M. A. (2012). Cyclicality or volatility? The satellite insurance market. Space Policy, 28(3), 192–198.

MCLANE, A. J., Semeniuk, C., McDermid, G. J., & Marceau, D. J. (2011). The role of agent-based models in wildlife ecology and management. Ecological Modelling, 222(8), 1544–1556.

NEAL, J. (2019). Lloyd’s of London, Pocket Guide. Lloyd’s of London.

NEVINS, M. R., Macal, C. M., & Joines, J. C. (1998). A discrete-event simulation model for seaport operations. Simulation, 70(4), 213–223.

NORBERG, R. (1995). A time-continuous markov chain interest model with applications to insurance. Applied Stochastic Models and Data Analysis, 11(3), 245–256.

OWADALLY, I., Zhou, F., & Wright, D. (2018). The insurance industry as a complex social system: Competition, cycles and crises. Journal of Artificial Societies and Social Simulation, 21(4), 2.

OWADALLY, M. I., Zhou, F., Otunba, R., Lin, J., & Wright, I. D. (2019). An agent-based system with temporal data mining for monitoring financial stability on insurance markets. Expert Systems with Applications, 123, 270–282.

PIOTR, M. (2016). Underwriting cycles and crises. The World of the New Economy, 4, 76–81.

THOYTS, R. (2010). Insurance Theory and Practice. London: Routledge.

VENEZIAN, E. C. (1985). Ratemaking methods and profit cycles in property and liability insurance. The Journal of Risk and Insurance, 52(3).

WEISS, M. A., & Chung, J.-H. (2004). U.S. reinsurance prices, financial quality, and global capacity. The Journal of Risk and Insurance, 71(3), 437–467.

WINTER, R. A. (1994). The dynamics of competitive insurance markets. Journal of Financial Intermediation, 3(4), 379–415.

YUN, T.-S., & Moon, I.-C. (2020). Housing market agent-based simulation with loan-to-value and debt-to-income. Journal of Artificial Societies and Social Simulation, 23(4), 5.

ZHOU, F. (2013). Application of agent based modeling to insurance cycles. PhD Thesis, City University London.