California earthquake insurance unpopularity: the issue is the price, not the risk perception

. Despite California being a highly seismic prone region, most homeowners are not covered against this risk. This study analyses the reasons for homeowners to purchase insurance to cover earthquake losses, with application in California. A dedicated database is built from 18 different data sources about earthquake insurance, gathering data since 1921. A new model is developed to assess the take-up rate based on the homeowners’ risk awareness and the average annual insurance premium amount. Results suggest that only two extreme situations would lead all owners to cover their home with insurance: (1) a widespread belief that a devastating earthquake is imminent, or alternatively (2) a massive decrease in the average annual premium amount by a factor exceeding 6 (from USD 980 to 160, 2015 US dollars). Considering the low likelihood of each situation, we conclude from this study that new insurance solutions are necessary to ﬁll the protection gap.


Introduction
Since 2002, the rate of homeowners insured against earthquake risk in California has never exceeded 16 %, according to data provided by the California Department of Insurance (CDI, 2017b, California Insurance Market Share Reports).Such a low rate is surprising in a rich area prone to earthquake risk.Consequently, several studies have already investigated homeowners' behaviour regarding earthquake insurance in California, in order to identify people who might have an interest in purchasing earthquake insurance and to understand why they do not do so.Three main variables have been observed as decisive in purchasing an earthquake insurance: the premium amount, the socio-economic background and the risk perception (Kunreuther et al., 1978;Palm and Hodgson, 1992a;Wachtendorf and Sheng, 2002).For the premium amount, both a survey conducted by Meltsner (1978) and the statistics collected from insurance data by Buffinton (1961) and Latourrette et al. (2010) show, as expected, a negative correlation between the premium amount and the insurance adoption.Nevertheless, even uninsured homeowners tend to overestimate the loss they would face in case of a major event and feel vulnerable regardless of their income level (Kunreuther et al., 1978;Palm and Hodgson, 1992a).As they are not expecting much from federal aids, they do not know how they will recover (Kunreuther et al., 1978).
Consequently the decision to purchase insurance to cover earthquake losses is not correlated to the income level (Kunreuther et al., 1978;Wachtendorf and Sheng, 2002).Several other socio-economic factors (e.g. the duration of residence, the neighbours' behaviour, the communication strategies of mass media, real estate agencies and insurance companies) have an impact on the risk perception (Meltsner, 1978;Palm and Hodgson, 1992a;Lin, 2013).
Pricing methods used in the real estate market do not take into account seismic risk.Indeed, based on data from the World Housing Encyclopedia, Pothon et al. (2018) have shown that building construction price is not correlated to the seismic vulnerability rating.Moreover, earthquake insurance is not mandatory in California for residential mortgage (CDI, 2019, Information Guides, Earthquake Insurance).Still, lenders for commercial mortgages are used to requiring earthquake insurance only when the probable maximum loss is high (Porter et al., 2004).However, Porter et al. (2004) have shown that taking into account earthquake risk for calculating the building's net asset value can have a significant impact in earthquake-prone areas like California.
Also, people insured against earthquake risk can receive a compensation lower than the loss incurred because most earthquake insurance policies in California include a deductible amount (i.e. at the charge of the policyholder) and are calibrated based on a total reconstruction cost, declared by the policyholder.As reported by Marquis et al. (2017) after the 2010-2011 Canterbury (NZL) earthquakes, this amount can be inadequate for the actual repair costs.Both high deductible amount (Meltsner, 1978;Palm and Hodgson, 1992b;Burnett and Burnett, 2009) and underestimated total reconstruction cost (Garatt and Marshall, 2003) can make insured people feel unprotected after an earthquake, making this kind of insurance disrepute.
Since homeowners are aware of the destructive potential of a major earthquake, the risk perception reflects their personal estimate of the occurrence (Kunreuther et al., 1978;Wachtendorf and Sheng, 2002).Indeed, most of them disregard this risk because it is seldom, even if destructive earthquake experiences foster insurance underwriting during the following year (Buffinton, 1961;Kunreuther et al., 1978;Meltsner, 1978;Lin, 2015).Consequently, earthquake insurance purchasing behaviour is correlated with personal understanding of seismic risk and is not only related to scientific-based hazard level (Palm and Hodgson, 1992a;Palm, 1995;Wachtendorf and Sheng, 2002).
With limits in earthquake insurance consumption identified, the next step is to assess the contribution of each factor to the take-up rate.Kunreuther et al. (1978) first proposed a sequential model showing that the individual's purchasing behaviour depends on personal risk perception, premium amount and knowledge of insurance solutions.In the same study (Kunreuther et al., 1978), a numerical model was developed, at the same granularity, failing however to accurately reproduce the observed behaviours.According to the authors, the model failed because many surveyed people had a lack of knowledge of existing insurance solutions or were unable to quantify the risk.Models published later (Latourrette et al., 2010;Lin, 2013Lin, , 2015) ) assessed the take-up rate by postal code and used demographic variables to capture the disparity in insurance solution knowledge.Nevertheless because of lack of data, they assumed the premium amount or the risk perception as constant.
The main objective of this study is to introduce a new take-up rate model for homeowners at the scale of California state.Such spatial resolution allows work on most data available in financial reports, which are numerous enough to use both the premium amount and the risk perception as variables.Different from previous studies focusing on risk pricing (Yucemen, 2005;Petseti and Nektarios, 2012;Asprone et al., 2013), this study takes into account the homeowners' risk perception and behaviour.This shift in point of view approach changes the main issue from what the price of earthquake insurance should be considering the risk level to what the acceptable price for consumers to purchase an earthquake insurance cover is.Despite results being at the state level, they bring a new framework to model earthquake insurance consumption which allows us to quantify the gap between premium amount and homeowners' willingness to pay, depending on the homeowners' risk awareness.Last, this study is also innovative by separately modelling the contribution of the risk perception and the premium amount to the level of earthquake insurance consumption.
In the first section, data are collected from several sources, processed and summarized in a new database.A take-up rate model is then developed by introducing the two following explanatory variables: the average premium amount and the subjective annual occurrence probability, defined as the risk perceived by homeowners of being affected by a destructive earthquake in 1 year.The perceived annual probability of occurrence is then studied from 1926 to present, in the next section.Finally, the last section is dedicated to analysing the current low earthquake insurance take-up rate for California homeowners.

Data collection and processing
Developing such a model, even at the state level, faces a first challenge in data collection.Despite the California earthquake insurance market having been widely analysed, data volume before the 1990s remains low and extracted from scientific and mass media publications which refer to original data sources that are no longer available.As a consequence, data description is often sparse and incomplete and values can be subject to errors or bias.Consequently, data quality of each dataset has to be assessed.Here, it is done on the basis of the support (scientific publication or mass media), the data description (quality of information on the data type and the collecting process) and the number of records.The datasets with the highest quality are used and presented in Table 1 (values are available in the Supplement).
This study focuses on the following variables about earthquake insurance policies: the total written premium (W N ), corresponding to the total amount paid by policyholders to insurance companies during the year N; the take-up rate (t N ), defined as the ratio between the number of policies with an earthquake coverage (Nb N ) and those with a fire coverage (Fi N ); the annual average premium (P N ), equal to W N over Nb N ; the average premium rate (p N ) equal to the ratio between the annual average premium amount P N and the average value of the good insured (here a house), later referred to as the average sum insured (ASI N ).
Nat All these variables can be calculated for any set of earthquake insurance policies.Insurance companies usually classify their products as follows: the residential line of business (Res) corresponds to insurance policies covering personal goods (e.g.house, condominium, mobile home, jewellery, furniture).The homeowners line of business corresponds to the insurance policies included in the residential line of business but dedicated to homeowners.As an earthquake insurance cover can be either a guarantee included in a wider policy (e.g. also covering fire or theft risks) or issued in a stand-alone policy, the earthquake line of business (Eq) classifies all insurance policies exclusively covering the earthquake risk.Some insurance policies within the earthquake line of business are dedicated to professional clients and do not belong to the residential line of business.
For clarity, variables hold the title of the line of business when they are not related to the homeowners line of business.Data for all the variables listed are available only since 1996 (Table 1).In order to expand the historical  database, they are also estimated in an indirect way for the following periods : 1926-1930, 1956, 1971-1972, 1976-1978 and 1989-1995, leading us to consider additional data (Table 1) and to use linear regressions (Fig. 1).By definition of the average premium rate (p N ), P N can be approximated by the product of p N and the average sum insured (ASI N ).The latter is estimated from the real building cost index (RBCI N ), which is an economic index (base: 31 December 1979 = 100) capturing the evolution of the cost of building construction works (Shiller, 2018).The following linear regression has been calibrated on data from Table 1 between 2006 and 2015: where ASI N is in 2015 US dollars and R 2 is the coefficient of determination.When the average premium rate corresponds to the residential line of business (p Res N ), it can be converted into p N with the following linear regression built on the CDI database between 2006 and 2016 (Table 1): Again using the CDI database between 1996 and 2016 (Table 1), the following linear regressions have also been developed between P N , t N , and W N and the corresponding metrics for the residential line of business (P Res N , t Res N and W Res N ).Before 1996, the total premium amount for earthquake insurance, mentioned in CDI (2017b) reports, was about the earthquake line of business, and in terms of total earned premium (E N ).While the total written premium (W N ) corresponds to the total amount of premium paid by policyholders to insurance companies during the year N, E N is the amount of premium used to cover the risk during the year N. To illustrate these two definitions, Fig. 2 takes the example of an insurance policy for which the annual premium is paid 1 March every year.As the amount received at year N by the insurance company (W N ) covers the risk until 1 March, N + 1, the insurance company can use only 75 % of W N during the year N (9 months over 12).By adding the 25 % of W N −1 , this makes the total earned premium E N , for the year N.
To estimate W

Eq
N from E Eq N for the earthquake line of business, the following relationship has been defined based on data from the CDI (2017a, b) database and over the period 1992-2016 (Table 1): (7) This strong assumption was used by Garamendi et al. (1992) and was verified for the loss after the 1994 Northridge earthquake reported by Eguchi et al. (1998).
Lastly, when variables cannot be inferred from other variables, they are estimated based on the annual variation for P Res N , and the biennial variation for t Res N , where Y X is equal to the ratio Y over X.Furthermore, the variation in the total written premium for the residential line of business ( ) is linked to as follows: As we have no data regarding the annual variation in the number of fire insurance policies for the residential line of business ( ) before 1996 (Table 1), it is assumed to be equal to the variation in the California population Pop N : Figures 1 and 3 illustrate the goodness of fit of the regressions developed (Eqs. 1 to 6) and the variations in the population compared to the number of fire insurance policies between 1996 and 2016 (Eq.11), respectively.Finally, values of P N , t N and W N collected (Table 1) and estimated in this study (Eqs. 1 to 11) are aggregated into a new database presented in Table 2.
Financial amounts are converted into 2015 US dollars using the Consumer Price Index (US Department of Labor, 2018).The data quality of each variable is also assessed at the weakest data quality used (Table 1), by as many levels as the number of equations used to assess it (listed in the column equation ID).For example, the annual average premium P 1926 has been calculated from the premium rate for the residential line of business p Res 1926 (Table 1) and Eqs. ( 1) and ( 4).The associated data quality is C (Table 2) because the data quality of p Res 1926 is A (Table 1), downgraded by two levels for the two equations used.
The model developed in the next section uses this new database to estimate the take-up rate for the California homeowners (t N ) from the average premium amount (P N ) and another variable capturing the relative earthquake risk awareness.
3 Model development for the period 1997-2016 In order not to presume a linear trend between the consumers' behaviour and explanatory variables (e.g. the average premium amount), this study refers to the expected utility theory (Von Neumann and Morgenstern, 1944) instead of statistical linear models.The expected utility theory is a classical framework in economic science to model the economical choices of a consumer, depending on several variables like the wealth or the risk aversion (Appendix A).Here the application field is the California homeowners instead of a single consumer.Within this framework, homeowners are assumed to be rational and making decisions in order to maximize their utility.Furthermore, their relative risk aversion is considered constant whatever the wealth because the decision to purchase earthquake insurance is independent of the income level (Kunreuther et al., 1978;Wachtendorf and Sheng, 2002).One of the most used utility functions (U ) is (Holt and Laury, 2002) where g N and β are the wealth function at year N and the risk profile controlling the risk aversion level (the larger β, the higher the aversion), respectively.Then, the average household's capital K is assumed to be equal to the average sum insured (ASI 2015 ) since the earthquake insurance consumption is uncorrelated to the wealth   (Kunreuther et al., 1978;Wachtendorf and Sheng, 2002) and equal to USD 604 124 (USD 2015).Consequently, the wealth of uninsured homeowners is g N = K, if no damaging earthquake occurs.Regarding the loss estimation, homeowners are mostly concerned by destructive earthquakes, defined as earthquakes which can potentially damage their home, and tend to overestimate the impact (Buffinton, 1961;Kunreuther et al., 1978;Meltsner, 1978;Lin, 2015).Furthermore, according to Kunreuther et al. (1978), insured people believe that they are fully covered in case of loss.Therefore, we assume that only uninsured homeowners expect to incur a loss equal to K after a damaging earthquake.This leads us to model the fact that uninsured homeowners expect a wealth at g N = 0 after a damaging earthquake.About insured homeowners, we assume that they expect to have a constant wealth at g N = K − P N since they believe to be fully covered in case of losses induced by an earthquake.Table 3 summarizes Table 3.The four different wealth levels of a homeowner, considered in this study, according to the occurrence of a damaging earthquake and the insurance cover.The quantities t N and 1 − t N represent the share of homeowners insured and not insured, respectively.
Occurrence of Insured Not a damaging (t N ) insured earthquake (1 − t N ) these three wealth levels, based on the occurrence of a devastating earthquake and the insurance cover.
Finally, taking into account the share of insured homeowners (t N ) and the results in Table 3, g N can be written for all homeowners as where B EQ (r N ) is a random variable following the Bernoulli distribution (i.e.equal to 1 if a damaging earthquake occurs during the year N and 0 otherwise).The parameter r N , called the subjective annual occurrence probability, controls the homeowners' risk perception through the perceived probability of being affected by a destructive earthquake during the year N (Kunreuther et al., 1978;Wachtendorf and Sheng, 2002).As homeowners want to maximize their utility, t N is solution of where argmax stands for the argument of the maxima function (i.e. which returns the value of t N , which maximizes the quantity E[U (g N (t N ))]) and E is the expected value of U (g N (t N )), which depends on the random variable B EQ (r N ) (Eq. 13).Assuming that homeowners are riskaverse (i.e.β > 0) and noticing that P N K is very small compared to 1, t Estimated N is shown (Appendix B) to be equal to As a first step, the model is calibrated with data for the period 1997-2016 corresponding to the whole activity period of the California Earthquake Authority (CEA) providing high-quality data (Table 2).Created on December 1996, this is a non-profit, state-managed organization selling residential earthquake policies (Marshall, 2017).The stability of the insurance activity and the lack of devastating earthquakes leads us to model parameters β and r N , capturing the homeowners' risk perception, by constant variables.They are estimated at r 1997;...;2016 = 0.027 % and β = 0.93, using the least squared method with the generalized reduced gradient algorithm.The model (Eq.15) fits the observed data with a R 2 = 0.79 (Fig. 4).Furthermore, the values of the two parameters are meaningful: homeowners are modelled as "very risk-averse" (β = 0.93), which is consistent with a high-severity risk like an earthquake (Holt and Laury, 2002).r N is also somehow consistent with the value that can be calculated from a hazard analysis, as presented later.
The model calibration between 1997 and 2016 is not appropriate for other periods corresponding to different seismic activity and insurance economic context.The next section investigates how the homeowners' risk perception has changed since 1926 in order to adapt the model to any period.

Evolution of the homeowners' risk perception since 1926
The homeowners' risk perception is controlled by both the risk materiality (what kind of earthquakes are expected to occur?) and the risk tolerance (how much are homeowners ready to lose?).The associated parameters in the model are r N and β.As we cannot differentiate the one from the other, β is assumed constant in this study and all the variations in the homeowner's risk perception are passed on to r N .According to the model and the relationships developed (10)(11), the variations in the total written premium amount per capita (W N /Pop N ) depend on the variation in the average premium amount (P N ) and the subjective annual occurrence probability (r N ):  Under the assumption that t Estimated N has never reached 100 % in the past and with β = 0.93, Eq. ( 16) can be simplified as follows: The contribution of P N /P N−1 is clearly marginal, leading to consider the regression, built on the variations in W N /Pop N and r N between 1997 and 2016 (Table 2): Equations ( 16) and ( 18) used and the goodness of fit are illustrated in Fig. 5a and b, respectively.Considering that understanding the details of those variations is out of reach given the scarcity of data, this study focuses only on those above 15.5 % in absolute value, qualified hereafter as significant variations (V N ).Other variations are neglected and we assume they cancel each other out.Table 4 lists all of them, together with the most significant event that  : 1933Long Beach, 1952Kern County, 1966Parkfield, 1971San Fernando, 1979Coyote Lake, 1979Imperial Valley and Coyote Lake, 1983Coalinga, 1984Morgan Hill, 1987Whittier Narrows, 1989Loma Prieta, 1992Big Bear and Landers, 1994Northridge, 2003San Simeon, 2004Parkfield and 2014 Napa.occurred the same year, for the earthquake insurance industry.Despite the fact that these events (Table 4) could have led insurance companies to restrict or enlarge the number of earthquake insurance policies (Born and Klimaszewski-Blettner, 2013), contributing in this way to V N , this study focuses only on the homeowners' risk perception variations.
Then, the subjective annual occurrence probability (r N ) is estimated according to the following time series for N ∈ [1926; 2016]: The time series of r N is represented in Fig. 6 and illustrates that some earthquakes (indicated by dots) have increased r N during the year as already published (Buffinton, 1961;Kunreuther et al., 1978;Meltsner, 1978;Lin, 2015), but, more surprisingly, some others had no apparent impact, like the 1952 Kern County and 1989 Loma Prieta earthquakes.The first earthquake damaged over 400 unreinforced masonry buildings (Jones et al., 2012), but none of the buildings were designed under the latest seismic codes (Geschwind, 2001).Consequently, structural engineers claimed that new buildings were earthquake-proof and that no additional prevention measures were needed (Geschwind, 2001).This message was received well among the population despite the warnings of some earthquake researchers like Charles Richter (Geschwind, 2001).

Category
nitude greater than M5 reached the highest peak since 1855 between 1984 and 1998, as illustrated in Fig. 7.The sequence includes in particular the 1983 Coalinga and the 1987 Whittier Narrows earthquakes.This specific temporal burst in seismic activity may have contributed to homeowners' rising risk awareness (Table 4).Some major socio-economic events also had an impact on r N (Table 4), such as the 1929 Great Depression and the post World War II economic expansion.Insurance regulation acts also have an impact.During the period 1984-1985, the California senate assembly voted on the Assembly Bill AB2865 (McAlister, 1984), which required insurance companies to offer earthquake coverage.A second bill in 1995, named AB1366 (Knowles, 1995), amended the Insurance Code Section 10089 to restrict the extent of the mandatory cover and addressed the insurance crisis caused by the 1994 Northridge earthquake.The new "mini-policy" met with limited success mainly because the price was more expensive for a lower guarantee (Reich, 1996).Furthermore, the creation of the CEA was at odds with the trend in insurance market privatization and people perceived it as an insurance industry bailout (Reich, 1996).The combined impact on r N of the two assembly bills is null as (1 4).It means that the efforts to promote earthquake insurance among the population by a first assembly bill were cancelled 20 years later by the other.Last, the discredit of a mistaken prediction can lead to a decrease in r N , as happened in 1928 with the occurrence of a major earthquake in Southern California (Geschwind, 1997;Yeats, 2001).
r N is then compared to a scientifically based earthquake annual occurrence probability to assess the level of the homeowners' risk perception.This requires the assessment of the annual probability for a homeowner of being affected by a destructive earthquake.For that, the ShakeMap footprints, released by the US Geological Survey (Allen et al., 2009), are used.A ShakeMap footprint gives, for a historical earthquake, the modelled ground motions for several metrics, including the macroseismic intensity on the modified Mercalli intensity scale (MMI).In California, a total of 564 ShakeMap footprints are available on the USGS website (https://earthquake.usgs.gov/data/shakemap/)for the period 1952-2018 and the magnitude range M3-M7.3.According to the MMI scale, heavy damage can be observed for an intensity above or equal to VIII.Among the 564 historical earthquakes, the ShakeMap footprints report that only 21 have caused such a high intensity and are labelled as destructive earthquakes.For each of the 21 ShakeMap footprints, the area corresponding to an intensity above or equal to VIII was reported in Table 5 and the population living inside was estimated using the Global Settlement Population Grid (European Commission; Florczyk et al., 2019).
Finally, for each year since 1952, the annual rate of people was calculated by dividing the number of people having experienced an intensity above or equal to VIII by the total population of California according to the Global Human Settlement Population Grid (European Commission, Florczyk et al., 2019).The arithmetic mean of the annual rates during the period 1952-2018 is equal to r Hist = 0.038 %.In this study, we consider r Hist as the true average probability for a Cali-  fornia homeowner to be affected by a destructive earthquake.Therefore, the risk has always been underestimated by California homeowners because AWR N has been lower than r Hist since 1926 (Fig. 6).The severity underestimation is quantified using the earthquake risk awareness ratio (AWR N ) defined as With AWR N being estimated since 1926, the take-up rate model is generalized in the next section and used to understand the current low take-up rate.

Understanding the current low take-up rate
Introducing the earthquake risk awareness ratio (AWR N ), we redefine the take-up rate model (Eq.15) as follows: From the estimated values of AWR N (Eqs.19 and 20) and P N (Table 2), Fig. 8 shows the fit between t Estimated N (Eq.21) and t N (Table 2).The goodness of fit is R 2 = 0.99 and the modelling errors are below ±3 %, as represented by the dashed lines.The record at 100 % corresponds to the 1992 California Residential Earthquake Recovery Fund Program (CRER), which was a mandatory public earthquake insurance, including a cover amount for households up to USD 15 000 (USD 1992), at a very low premium amount (Lagorio et al., 1992).
Figure 9 illustrates the sensitivity of t Estimated N to the average premium P N and the earthquake risk awareness ratio AWR N (Eq.20).
The lines AWR N = 70 %, AWR N = 100 % and AWR N = 425 % stand for the current situation, the true risk level and the target value of AWR N for t Estimated N = 100 %, at the current price of USD 980 (USD 2015), respectively.According to the first value (AWR N = 70 %), half of homeowners (t Estimated N = 50 %) are not willing to pay an average premium amount exceeding P N = USD 310 (USD 2015) per year for an insurance cover.This would represent a 68 % decrease in the price of earthquake insurance coverage.
To reach a 100 % take-up rate with the current premium amount, according to the model, the earthquake risk awareness ratio has to reach 425 %, which is very unlikely.Indeed, the UCERF3 assesses the annual probability of occurrence, in California, of an earthquake with a magnitude above M6.7 at 15 % (Field et al., 2013).Under the assumption that the probability for a homeowner of experiencing a damaging earthquake is proportional to the occurrence of such an earthquake, reaching that awareness ratio would mean that homeowners consider this probability to reach the level of 425 % × 15 % = 64 %.In other words, only the belief that a destructive earthquake is imminent (i.e. a probability of occurrence somewhere in California during the year, for a 1994 Northridge-like earthquake, above 64 %) can bring all homeowners to buy an earthquake insurance plan at the current price.

Conclusions
The model developed in the present study shows that the low earthquake insurance take-up rate, observed until 2016 for the homeowner line of business in California, is due foremost to high premiums.Indeed, no homeowner would prefer to stay uninsured against earthquake risk if the annual average premium were to decrease from USD 980 (as observed in 2016) to USD 160 (USD 2015) or lower.Moreover, Kovacevic and Pflug (2011) have shown that a minimum capital is required to benefit from an earthquake insurance cover, otherwise the cost of the premium would be more than the potential risk.They have also found that the lower the insurance premium, the lower this minimum capital.Consequently, the results of Kovacevic and Pflug (2011) suggest that some uninsured people can be well aware of the earthquake risk but just cannot afford earthquake insurance at the current price.This corroborates the result of the present study because bringing such people to buy earthquake insurance is foremost a matter of price, not of insufficient risk perception.
Nevertheless, as the current average premium amount corresponds to the annualized risk assessment, insurance companies would not have enough reserves to pay future claims following earthquakes if they collect on average only USD 160 per policy (CEA, 2017).Hence, the current insurance mechanism cannot meet the homeowners' demand and be sustainable, as was the case for the 1992 CRER programme (Kunreuther et al., 1998).
Since this initiative, earthquake insurance has not been mandatory in California.Because compulsory insurance is not only hard to implement but also has an unpredictable impact (Chen and Chen, 2013), new insurance solutions have emerged to increase the earthquake insurance market, like parametric insurance.This insurance product stands out from traditional insurance policies by paying a fixed amount when the underlying metrics (e.g. the magnitude and the epicentre location) exceed a threshold, whatever the loss incurred by the policyholder.This new insurance claim process significantly reduces the loss uncertainty, the operating expenses and so the premium rate.In California, parametric insurance has been offered to cover homeowners against earthquake risk since 2016 (Jergler, 2017).
Nevertheless, risk awareness remains important in earthquake insurance consumption, and insurance companies also encourage local policies led by public authorities to improve prevention (Thevenin et al., 2018).
This study supports these initiatives by demonstrating that a new earthquake insurance scheme is necessary to meet the premium expected by homeowners.However, it also shows that the lower the premium amount, the more the earthquake risk awareness contributes in insurance consumption (Fig. 9).Finally, only a global effort involving all stakeholders may fill the adoption gap in earthquake insurance coverage.This also includes the financial stakeholders, involved in the mortgage loan business.A recent study (Laux et al., 2016) analysed the fact that most banks do not require earthquake insurance from their mortgage loan clients, preferring instead to keep the risk and to increase the interest rate by +0.2 % (in 2013).These additional bank fees are, on average, as expensive as the premium rate for the residential earthquake insurance (0.185 % in 2013 according to the CDI, 2017a).Nevertheless, such a market practice puts an estimated USD 50 billion-100 billion loss risk on the US mortgage system (Fuller and Kang, 2018).Consequently, expanding earthquake insurance cover, for homeowners, is also a challenge for the banking industry.In this respect, requiring earthquake insurance for mortgage loans would strongly enhance the insurance and public initiatives.
where β i is the level of relative risk aversion and materializes the decision maker's risk behaviour from highly risk-loving to highly risk-averse.Indeed, one can verify that the RRA measure for the Holt and Laury (2002) utility function is RRA(g 1 ) = β i .Therefore, this function belongs to the constant risk relative aversion utility function family.It means that a decision maker's behaviour is assumed to be uncorrelated to their wealth level (g 1 ).Furthermore, the parameter β i controls the decision maker's behaviour (Table A1).

Figure 1 .
Figure1.Fits of the linear regressions for (a) the average sum insured and the real building cost index, (b) the average premium rate between the residential and the homeowner lines of business, (c) the take-up rate between the residential and the homeowner lines of business, (d) the average premium amount between the residential and the homeowner lines of business, (e) the total written premium amount between the residential and the homeowner lines of business, and (f) the written and the earned premium amounts.Grey points on (a) and (c) represent the extreme values, affecting the linear regression from the solid black line to the dashed grey line, when removed.Financial values are in 2015 US dollars.

Figure 2 .
Figure 2. Illustration of the difference between the earned and the written premium amount.

Figure 3 .
Figure 3.Comparison between the annual variations in the number of fire insurance policies for the residential line of business (solid line) and the California population (dashed line).

Figure 4 .
Figure 4. Fit between the take-up rates modelled and observed for the period 1997-2016.

Figure 5 .
Figure 5. (a) Variations in the average premium amount, the total written premium per capita and the subjective annual occurrence probability between 1997 and 2016.Variations are calculated compared to the previous year.Dashed lines represent the thresholds for a significant variation.(b) Fit between the variation in the total written premium per capita and the subjective annual occurrence probability.

Figure 6 .
Figure 6.Estimated subjective annual occurrence probability between 1927 and 2016.The horizontal dashed line corresponds to the historical average value.Black bars represent the significant evolutions and dots indicate the occurrence of the following big earthquakes in California: 1933 Long Beach, 1952 Kern County,1966 Parkfield, 1971 San Fernando, 1979 Coyote Lake, 1979 Imperial Valley and Coyote Lake, 1983 Coalinga, 1984 Morgan Hill,  1987 Whittier Narrows, 1989 Loma Prieta, 1992 Big Bear and  Landers, 1994 Northridge, 2003 San Simeon, 2004 Parkfield and  2014 Napa.

Figure 7 .
Figure 7. Number of earthquakes with a moment magnitude M5+ occurring during the 5 previous years.For year N , the sum is calculated from N −5 to N −1.Earthquakes occurring the same day count as one.Black bars correspond to the peak of seismicity observed between 1984 and 1988.The dotted line represents the average value from 1855 to 2012.The historical earthquake database used is taken from the UCERF3(Field et al., 2013;Felzer, 2013).

Figure 8 .
Figure 8.Comparison between the take-up rates observed and estimated using the model developed in this study.Grey points correspond to the period 1997-2016, which was used to calibrate the model.Dashed lines represent a buffer of ±0.03.

Figure 9 .
Figure 9. Relationship between the average premium amount, the take-up rate and the earthquake risk awareness.Financial values are in US dollars (2015).

Table 2 .
The California homeowners earthquake insurance database developed for this study.Data quality scale is as follows.A (good): all data from scientific publications with methodology explained; (B) acceptable: at least one data point from a scientific publication without details on the methodology or one processing has been applied; (C) weak: data are uncertain.Financial amounts are in 2015 US dollars.The total written premium is given in millions of US dollars.Column names are Q: data quality; E: equation ID.Abbreviations are E 1 = {3; 6; 7; 10} and E 2 = {5; 6; 7}.

Table 5 .
Florczyk et al., 2019)occurred in California since 1950 with a maximum macroseismic intensity on the modified Mercalli scale (MMI) modelled by the ShakeMap programme (US Geological Survey) above or equal to VIII.The columns "area" and "number of people affected" give the size and the estimated population of the area affected above or equal to VIII, respectively.The population has been assessed using the Global Human Settlement Population Grid (European Commission,Florczyk et al., 2019).