A Global Analysis of the Relationship Between Urbanization and Fatalities in Earthquake-Prone Areas

Urbanization can be a challenge and an opportunity for earthquake risk mitigation. However, little is known about the changes in exposure (for example, population and urban land) to earthquakes in the context of global urbanization, and their impacts on fatalities in earthquake-prone areas. We present a global analysis of the changes in population size and urban land area in earthquake-prone areas from 1990 to 2015, and their impacts on earthquake-related fatalities. We found that more than two thirds of population growth (or 70% of total population in 2015) and nearly three quarters of earthquake-related deaths (or 307,918 deaths) in global earthquake-prone areas occurred in developing countries with an urbanization ratio (percentage of urban population to total population) between 20 and 60%. Holding other factors constant, population size was significantly and positively associated with earthquake fatalities, while the area of urban land was negatively related. The results suggest that fatalities increase for areas where the urbanization ratio is low, but after a ratio between 40 and 50% occurs, earthquake fatalities decline. This finding suggests that the resistance of building and infrastructure is greater in countries with higher urbanization ratios and highlights the need for further investigation. Our quantitative analysis is extended into the future using Shared Socioeconomic Pathways to reveal that by 2050, more than 50% of the population increase in global earthquake-prone areas will take place in a few developing countries (Pakistan, India, Afghanistan, and Bangladesh) that are particularly vulnerable to earthquakes. To reduce earthquake-induced fatalities, enhanced resilience of buildings and urban infrastructure generally in these few countries should be a priority.


Introduction
The process of urbanization remains an ongoing global phenomenon. Migration to and concentration of global population in urban areas has reached 55% with an estimated 68% level anticipated by 2050 (UN 2019). The outcome of this concentration of people, buildings, and other assets attributed to urbanization (Bloom et al. 2008) is an increase in the exposure of populations to natural catastrophes, for example, earthquakes (UN 2015), floods (Pesaresi et al. 2017), and disease spread (Alirol et al. 2011). Due to the high cost of living in proximity to urban area (Ravallion and van de Walle 1991;Shucksmith et al. 2009), the most vulnerable proportion of the population often lives in substandard residential building (Lankao and Qin 2011;UNHCR 2014;Nagendra et al. 2018).
While the growth of urban population is substantial (58% increase in global urban population from 1990 to 2015), it has been outpaced by the rate of urban expansion (74%) UNDP 2018). Urban expansion often involves the conversion of agriculture and forest lands to low-density, single-story, residential building land uses (Spence 2007;Güneralp et al. 2020). The speed of building creation on the urban periphery can lead to construction that is often of questionable quality (Green 2008;Oteng-Ababio 2012), since the focus of development is often on gaining affordable access to urban occupational opportunities and amenities (Bilham 2009) rather than the quality of living produced by development.
For developing countries, urbanization is usually accompanied by a growth in vulnerable buildings and infrastructure (Jackson 2006;Spence 2007;Ambraseys and Bilham 2011;Wyss 2018), because proper construction codes may be lacking or poorly implemented (Bilham 2004;Aitsi-Selmi et al. 2015). Developing countries also experience the fastest rate of urbanization (Jiang and O'Neill 2017;Chen et al. 2020a), with many of those fastgrowing urban areas located in earthquake-prone zonesup to 31% of global built-up areas (Pesaresi et al. 2017). While global changes in population (Leyk et al. 2019) and land-use (Potere and Schneider 2007) have been recorded and their association with urbanization quantified (Montgomery 2008;IPCC 2018), the relationship between urbanization in earthquake-prone areas and its impacts on earthquake fatalities is focused at the regional and city scales (Sarris et al. 2010;Abramson and Qi 2011;Martins et al. 2012;Huang et al. 2018).
Little is known about the relationship between earthquake fatalities and global urbanization (Peduzzi et al. 2009;Daniell et al. 2011;UN 2015). While population exposure to earthquakes has increased (Djordjević et al. 2016;He et al. 2016;Dou et al. 2018) and built-up surfaces in earthquake-prone areas have increased globally from 1975 to 2015 (Pesaresi et al. 2017), an examination of the relationship between population size or urban land area and earthquake fatalities is lacking. Exploring this relationship is critical to predicting seismic risks and mitigating earthquake fatalities.
To overcome this gap in knowledge, we quantify the changes in global population and urban land area in earthquake-prone places and their impacts on associated fatalities in the context of urbanization. We also discuss the implications of this study in the context of expected global urbanization under different Shared Socioeconomic Pathways (SSPs) in the future.

Materials and Methods
We used multiple sources of data to examine the changes in population and urban land in global earthquake-prone areas and their association with earthquake fatalities. First, we quantified the change in population and urban land in earthquake-prone settings from 1990 to 2015. Then, we examined the relationship between urbanization (that is, population, urban land, and the urbanization ratio) and earthquake fatalities in the earthquake-prone areas during 1990-2018, using a series of regression models and the STIRPAT (STochastic Impacts by Regression on Population, Affluence, and Technology) approach (York et al. 2003) in which fatalities are estimated as a function of hazard, exposure, and vulnerability (Lavell et al. 2012;UNISDR 2017).

Data
To quantify the changes in global population and urban land in the earthquake-prone areas we utilized four types of data. First, population and urban land data were acquired from the HYDE 3.2 dataset at 10 km resolution (Goldewijk 2016). The History Database of the Global Environment (HYDE) 3.2 dataset is the only global dataset generated using a consistent approach to estimate historical and future population and urban land dynamics. Gridded HYDE 3.2 population and urban land data for 1990 and 2015 as well as 2050 under five Shared Socioeconomic Pathways (SSPs) were used. We report the mean, maximum, and minimum values of the five SSP scenarios to depict uncertainties and variation in the future van Vuuren et al. 2014).
Second, global peak ground acceleration (PGA) data were acquired to map earthquake-prone areas. Areas with a PGA equal to or greater than 0.2 g (gravitational acceleration), when the collapse probability of buildings starts to increase remarkedly, were defined as the earthquake-prone areas (Jaiswal et al. 2011a;Holzer and Savage 2013). The most updated PGA data and resulting map were used, which collated more than 30 national and regional maps from various institutions by the Global Earthquake Model (GEM) foundation (Pagani et al. 2018). The PGA data were generated using the OpenQuake engine (Pagani et al. 2014) and expert input by hundreds of global scientists and engineers over the last 20 years. Although the PGA map provided by the GEM may under or overestimate the PGA (Ishibe and Shimazaki 2012;Wyss 2015) and alternatives were proposed (for example, Panza and Bela 2020), it still remains the most up-to-date, only freely accessible, and widely used global map available (Gerstenberger et al. 2020).
Third, fatalities from earthquakes were derived from the global CATDAT (CATastrophy DATabase) database (Daniell et al. 2011), a catalog that emphasizes earthquakes, their triggering impacts, and downstream costs. Although alternative datasets on earthquake damage, such as the National Geophysical Data Center/World Data Service (NGDC/WDS 2018), Emergency Events Database (EM-DAT 2018), and PAGER-CAT (Prompt Assessment of Global Earthquakes for Response) catalog (Allen et al. 2009) datasets are available, CATDAT was used due to its accessibility, information completeness, and wide use. Furthermore, the CATDAT database is the largest known catalogue that has been validated against four major earthquake databases (PAGER-CAT, NGDC, UTSU (Utsu 2004), and MRNATHAN (Munich Re 2009)) earthquakeby-earthquake at the global scale including many other original data sources, which has resulted in substantial error reduction (Daniell et al. 2011). The CATDAT catalog records the impacts of damaging earthquakes, including direct and secondary effects since 1900 for over 9,900 earthquakes, from over 35,000 individual literature sources in 90 languages. Since it is difficult to quantify accurately the exact amount of earthquake damage, the database reports the preferred estimates and lower/upper bounds of earthquake damages after a careful analysis of all damage reports.
Fourth, auxiliary data on global administrative boundaries were obtained from the Resource and Environment Data Cloud Platform, Chinese Academy of Science. 1 Socioeconomic data were obtained from the World Bank. 2 All gridded data were reprojected to the Albers projection with a spatial resolution of 10 km.

Examining the Impacts of Urbanization on Fatalities in Earthquake-Prone Areas
Fatalities caused by earthquakes are typically estimated as a function of the earthquake hazard (mainly from ground shaking, but also from landslides and tsunami), the size of the population impacted, number of buildings exposed to earthquakes, and the vulnerability of exposed population and buildings (Lavell et al. 2012;UNISDR 2017). Similar to other approaches (for example, Peduzzi et al. 2009;Li et al. 2015), nationally reported data were used to represent these three components (hazard, exposure, and vulnerability). In terms of the seismic hazard, we chose the extent of the earthquake-prone areas, and the average depth and magnitude of earthquakes in the earthquake-prone areas during 1990-2015 for each country. For exposure to earthquakes, we used the average population and the average urbanized land in the earthquake-prone areas between 1990 and 2015.
Vulnerability is difficult to quantify globally in the absence of extensive building inventories and substantial socioeconomic data (Jaiswal and Wald 2010). When mapped globally, vulnerability is typically aggregated from national scale variables that act as proxies for specific variables of interest (Ruiter et al. 2017;Ward et al. 2020). Quantifying and mapping global vulnerability to earthquakes require the same approach as local factors, such as building quality and code, which are not always available at the national scale (Jaiswal and Wald 2010). To overcome this data gap, indicators such as the human development index (HDI) (Jaiswal et al. 2011b;Li et al. 2015) and vulnerability code (Daniell et al. 2011) are typically used.
We opted to use a five-class vulnerability code developed by the USGS (United States Geological Survey) for inclusion in their Prompt Assessment of Global Earthquakes for Response (PAGER) system (Jaiswal et al. 2011a). The vulnerability code is strongly correlated (r = -0.55, p\0.01) with HDI and considers the materials and methods used for building construction as well as the efforts and capacity to respond to earthquake disasters.
To estimate the impacts of urbanization in earthquakeprone areas on fatalities, fatalities were estimated, as in previous studies (Peduzzi et al. 2009;Li et al. 2015), by using a simple multiplicative function combining hazard, exposure, and vulnerability. Similar to the widely used environmental impact = f(population, affluence, technology) or IPAT model (York et al. 2003), we used the STIRPAT approach (York et al. 2003) to develop four models to examine the relative influence of each of the three components and their relationship with urbanization on earthquake fatalities.
The four models were used to examine the impacts of hazard, exposure, vulnerability, and urbanization on earthquake fatalities. First, we used a parsimonious model (Model I-hazard and vulnerability) to account for the impacts of hazard and vulnerability. The equation for this model can be expressed as follows: where the hazard was represented by the extent of the earthquake-prone areas and the average depth and magnitude of earthquakes in the earthquake-prone areas for each country, and vulnerability was denoted by the USGS vulnerability code. All the variables were used in their logarithmic form. Then, we incorporated the effects of population and urban land exposure as follows: in Model II-population exposure and Model III-urban land exposure. Based on previous findings (Jackson 2006;Spence 2007;Ambraseys and Bilham 2011), we expected an inverted Ushaped relationship between urbanization ratio and earthquake fatalities. To test this conjecture, we added the urbanization ratio and its quadratic term (both in logarithmic form) to generate Model IV-urbanization ratio, following the well-known environmental Kuznets curve hypothesis, expressed as follows: When a 4 is significantly positive and a 5 is significantly negative, an inverted U-shaped relationship between fatalities and urbanization ratio can be confirmed.
Changing the values of a 4 and a 5 can affect the shape of the inverted U-shaped curve between urbanization ratio and deaths.
In examining the global relationship between urbanization and earthquake-induced fatalities, our analysis is constrained to a relatively short time period of 25 years with data comprising a large amount of uncertainty. In terms of data uncertainty, we report the upper and lower bounds of fatalities of our regression results based on the CATDAT catalogue. In addition to the analysis using all the available data, we conducted an extra analysis by excluding fatalities from catastrophic earthquakes (for example, the 2010 Haiti earthquake) as the recorded fatalities are highly controversial (Doocy et al. 2013a), and the results are similar to our original analysis.

Results
We report the characteristics of earthquake-prone areas, the changes in population and urban land in earthquake-prone areas, and the relationship between exposure change and earthquake fatalities in this section.

Earthquake-Prone Areas
Total global earthquake-prone area was estimated to cover 1.36 9 10 7 km 2 , representing 9.2% of global terrestrial area (Figs. 1, 2). These areas were unevenly distributed globally with the largest proportion of earthquake-prone areas concentrated in several countries along the Circum-Pacific Belt and the Mediterranean-Asiatic Belt, corroborating others (England and Jackson 2011;Bilham 2014;Liu and Stein 2016). A total of 94 countries have earthquake-prone areas. Among these countries, the 10 countries (Iran, China, Peru, Indonesia, India, Pakistan, United States, Argentina, Chile, and Turkey) with the largest amount of earthquake-prone area have a total area of 7,730,000 km 2 , accounting for approximately 57% of the globe's earthquake-prone areas.
People living in global earthquake-prone areas were subject to disproportionally high levels of earthquake risk (Fig. 2). Globally, only a small proportion of population and urban land (20.8% and 15.5%, respectively) were located in earthquake-prone areas. Yet more than 56% of the earthquakes that caused casualties, nearly half of deaths caused by earthquakes, and more than 91% of deaths resulting from building collapse worldwide from 1900 to 2018 were attributed to earthquakes occurring in these earthquake-prone areas.
While only 9.2% of the global areas affected by earthquakes experience these higher levels of ground acceleration, 47.8% of all earthquake deaths and 91.7% of all deaths due to shaking reside within these earthquake-prone areas. Average population and urban land data between 1990 and 2015 are used to depict population and urban land exposure. Earthquakes, total deaths, and shaking deaths refer to the number of earthquakes, the total deaths and shaking deaths of earthquakes, respectively, during 1900-2018 recorded by CATDAT.

Population Change in Global Earthquake-Prone
Areas from 1990 to 2015 The population in global earthquake-prone areas increased at a faster rate than the global population (Table 1). The rate of population increase was 40.2% in global earthquake-prone areas, rising from 1075.94 million in 1990 to 1508.79 million in 2015. During the same period, global population increased by 36.6%. The result of a 9% faster growth rate of population in earthquake-prone areas is that, in 2015, those areas held more than 1/5 of the global population.
The population residing in urban areas relative to the total population of a country, or other spatial unit, defines a measurement known as the urbanization ratio (UN 2019). Population increases in global earthquake-prone areas were concentrated in countries with a relatively low urbanization ratio (Table 1, Fig. 3). Among the 94 countries that have earthquake-prone areas, more than 2/3 of the population increase from 1990 to 2015 occurred in these areas within 43 countries that have a low urbanization ratio (between 20 and 60%). The average earthquake-prone area population increase of the 43 countries (47.8%) was 1.2 times the average global population increase for all earthquake-prone areas. While rates of population increase in earthquakeprone areas were high, they were also highly variable (ranging from 37 to 72%). Furthermore, countries with the highest urbanization ratio had the lowest rate of population increase in earthquake-prone areas. Only those countries with the highest urbanization ratios-80-100%, for example, the United States, New Zealand, and Japan-had a rate of population increase (16.0%) substantially less than the global rate (36.6%).
For the five countries with the largest growth in earthquake-prone-area population (that is, India, Pakistan, Indonesia, the Philippines, and Bangladesh), the rate of Bars represent totals associated with all earthquakes while the proportion delineated in orange represents those in earthquakeprone areas (that is, with a peak ground acceleration equal to or greater than 0.2 gravitational acceleration; see the Data section (Sect. 2.1) for details). Data source CATDAT population increase ranged from 37% to 72% and the total growth accounted for more than half of the global growth in earthquake-prone area population (Figs. 3a,4a). It is no coincidence that all five of these countries reside in Asia. Seismic faults interface with drinking water sources along the Mediterranean-Asiatic Belt. These faults act as aquifers for many semiarid and arid countries in the Middle East and Central Asia, India, Pakistan, and Bangladesh (Jackson 2006). The combination of drinking water sources with the rapid population growth of Southeast Asian countries (for example, Indonesia and the Philippines) has led to the location of vulnerable populations along the Circum-Pacific Seismic Belt (UNHCR 2014).
Countries with earthquake-prone areas and low or high urbanization ratios (0-20% and 80-100%) in 2015 experienced low levels of population growth (Fig. 3). While the absolute population increase among the five countries with a low urbanization ratio only accounted for just over 17 million, their average rate of change was over 74% between 1990 and 2015. In contrast, the 14 countries with the highest urbanization ratio held an average population growth rate of 16%, substantially lower than the 37% global average (Fig. 3, Table 1).  3 Relationship between seismic exposure change and urbanization ratio in global earthquake-prone areas. a Population; b urban land. Scatterplots represent changes in total population or urban land for each country, while the bars represent the aggregated change for countries within each interval of urbanization ratio. According to the Northam Curve (Northam 1979;Mulligan 2013), the urbanization ratio of a nation can be divided into at least three stages, that is, the initial stage with an urbanization ratio between 0 and 20%, the growth stage with an urbanization ratio between 20% and 80%, and the terminal stage with a ratio above 80%

Urban Land Change in Global Earthquake-Prone Areas from 1990 to 2015
Land area classified as urban increased faster than population in global earthquake-prone areas from 1990 to 2015 (Table 1). The urban lands in global earthquake-prone areas expanded from 55,466 km 2 in 1990 to 93,867 km 2 in 2015, representing an increase of 69.2%. This increase was 1.7 times greater than the increase in population in those same areas. In addition, the expansion of urban lands in global earthquake-prone areas was 1.2 times greater than the global urban expansion. Consequently, the percentage of global urban land located in the global earthquake-prone areas increased from 15.1% to 16.8% from 1990 to 2015.
In contrast to the rate of population increase being greatest in earthquake-prone areas in countries with 20-60% urbanization ratio, the increase in global urban earthquake-prone land area was concentrated in a small number of countries with high urbanization ratios (60-100%) (Fig. 3b). Among the 94 countries with earthquake-prone areas, nearly 69% of the increased urban land in global earthquake-prone areas occurred in 46 countries with an urbanization ratio greater than 60% in 2015. The total growth of urban land in these countries' earthquakeprone areas reached 26,473 km 2 and the corresponding rate of increase was 63.6%, which is similar to the average rate in the earthquake-prone areas throughout the world. Among the 46 countries, the top five countries that experienced the largest increases of urban land in their earthquake-prone areas were the United States (48.4%), Iran (742.3%), Turkey (340.7%), Mexico (118.7%), and Japan (15.7%). The cumulative outcome of urban expansion in these five countries amounted to 17,872 km 2 , more than 12 times the earthquake-prone area urban expansion in Japan.
Large increases of urban land area in earthquake-prone areas also occur in some countries with an urbanization level lower than 60%, such as Indonesia, the Philippines, and India (Figs. 3b, 4b). The corresponding increases in urban land in the three countries' earthquake-prone areas were 2821, 1673, and 1524 km 2 , respectively, and all of them are developing countries in Asia. The growth of urban land in the earthquake-prone areas of the five countries with a low level of urbanization (urbanization ratio between 0 and 20%) was marginal, with a total amount of 232 km 2 , or 0.4% of the total growth in global earthquakeprone areas.
Earthquake-prone area urban land expanded slower in countries with a high urbanization ratio (80-100%) than countries with a lower level of urbanization ratio (Table 1). This difference suggests that to fulfill increasing residential demands in the earthquake-prone areas of developing countries more multistory buildings and crowded informal settlements are frequently constructed (Frolking et al. 2013;Henderson et al. 2016). The lack of building codes and guidelines for coping with earthquake strikes and postdisaster management put these vulnerable populations in danger (Doberstein and Stager 2013).

Earthquake Fatalities in Global Earthquake-Prone Areas during 1990-2018
Earthquake fatalities in global earthquake-prone areas during 1990-2018 for the 94 countries at risk were estimated at 413,576, with lower and upper bounds of 378,152 and 516,722, respectively (see Sect. 2.1 Data for the explanation of the lower/upper bounds). The top five countries that had the largest numbers of earthquake fatalities in their earthquake-prone areas were China, Pakistan, Haiti, India, and Afghanistan. Their earthquake fatalities in the earthquake-prone areas during this period were 88,710, 87,698, 80,009, 28,912, and 11,561, respectively, which together accounted for 71% of total deaths in global earthquake-prone areas. These countries had an urbanization ratio below 60% in 2015 and belong to the group of countries most vulnerable to earthquakes (Fig. 5).
As one would expect, earthquake fatalities were concentrated in the most vulnerable countries (Fig. 5). However, the proportion of all earthquake fatalities in these countries is striking. Approximately 96% (96.1-97.0%) of all deaths in global earthquake-prone areas occurred in countries (33) classified with the highest vulnerability score (a vulnerability code of 4 or 5). Seven of these countries (Pakistan, China, Iran, India, Afghanistan, Nepal, and Bangladesh) had the highest level of vulnerability (vulnerability code of 5) with total earthquake fatalities in earthquake-prone areas reaching 305,286 (280,531-321,280), accounting for approximately 74% (62.2-74.2%) of all fatalities in global earthquake-prone areas. The average number of deaths per thousand residents in countries with the highest level of vulnerability was 5.77 (5.30-6.07), which is 11.5 times (11.3-12.1) greater than the three countries with the lowest classification of vulnerability (New Zealand, Japan, and the United States).
The death rates per thousand population were disproportionally higher in countries with a low urbanization ratio and a high level of vulnerability (Fig. 6). The average number of deaths per thousand people in global earthquake-prone areas was 4.11 (3.79-5.37) in countries with an urbanization ratio below 60%. The average number of deaths per thousand people was nearly half (2.11) in countries with an urbanization ratio above 60% (1.88-2.36). Similarly, average deaths per km 2 urban land in global earthquake-prone areas was 10 times higher in countries with an urbanization ratio below 60%, 17.8 per km 2 (16.4-23.2), than those above 60% (1.9 per km 2 , or 1.7-2.0). The average deaths per km 2 of urban land in countries with the highest vulnerability (vulnerability codes of 4 and 5) was 140 times (137-147) greater than what was found in the three countries with the lowest vulnerability.

The Relationship Between Urbanization and Earthquake Fatalities
Overall, the average population living in earthquake-prone areas during 1990-2015 and the extent of the earthquakeprone areas were found to be significantly and positively associated with earthquake fatalities ( Table 2). The correlation coefficient between these two variables and earthquake induced deaths were 0.66 and 0.58, which were the most influential among eight tested variables (Table 2). A regression analysis further confirmed that the amount of population in the earthquake-prone areas and the extent of the earthquake-prone areas were positively associated with fatalities caused by earthquakes. Standardized regression coefficients for the average population in earthquake-prone areas varied from 0.53 to 0.92 from Model II to Model IV, whereas the coefficients for the extent of the earthquake-prone area were 0.31 and 0.60 in Model IV and Model I, respectively. In contrast, an increase in the average depth of earthquakes was significant and negatively associated with fatalities (Model I-IV), with standardized coefficients between -0.25 and -0.32. An increase in average amount of urban land in the earthquake-prone areas was significantly and negatively associated with earthquake deaths, whereby standardized of urbanization ratio (for example, 10-20% for the first box). Vulnerability was attributed using the vulnerability code in the PAGER system by USGS (see the Methods section (Sect. 2.2) for details) Fig. 5 Distribution of the fatalities caused by earthquakes in the earthquake-prone areas of countries during 1990-2018. Vulnerability was attributed using the vulnerability code in the USGS' PAGER (https:// earthquake.usgs.gov/data/pager/ ) system (see the Methods section (Sect. 2.2) for details) coefficients of -0.36 and -0.58 were acquired by Model III and IV, respectively. These findings corroborate existing literature, where population exposure and seismic hazard were strongly associated with earthquake fatalities (Jackson 2006;Peduzzi et al. 2009;Jaiswal et al. 2011b).
The introduction of the urbanization ratio and its quadratic term (Model IV) improved model fit. The eight variables in Model IV explained 58.2% of the cross-national variance in earthquake fatalities in the earthquakeprone areas, while the most parsimonious model (Model I), which only considered the impacts of the hazard and vulnerability, explained 40.9% of the variance. These results confirm that urbanization is critically related to earthquake fatalities. More importantly, the coefficient between the number of fatalities caused by earthquakes and the urbanization ratio is significant and positive, whereas the coefficient for the quadratic term of urbanization ratio is significant and negative, which suggests that an inverted U relationship similar to the environmental Kuznets curve (York et al. 2003) may exist. When the three components of risk (hazard, exposure, and vulnerability) are held constant, an increase in low to moderate levels of the urbanization ratio produces an increase in earthquake fatalities. When the urbanization ratio is moderate to high, however, a further increase in the urbanization ratio produces a decrease in earthquake fatalities (Fig. 7). The transition or point of inflection is estimated to occur when the urbanization ratio reaches around 44%. Twenty-four countries have an urbanization ratio B 44%; all are developing countries and all are located in Africa or Asia, with the exception of Papua New Guinea. These countries are still experiencing increasing seismic risk and will continue to do so if earthquake fatality mitigation strategies are not improved upon and operationalized.

Discussion
The potential causes for the nonlinear relationship between urbanization and earthquake fatalities, and its implications for mitigating earthquake risk in the future, are discussed in this section.

The Relationship Between Urbanization and Earthquake Fatalities
The nonlinear relationship between urbanization and earthquake fatalities may partly be a function of population   The number of samples used for the regression was 50 because some countries did not have earthquake with fatalities during the studied period. The values in bracket are the p values **p \ 0.01; *p \ 0.05; ? p \ 0.1 density. For developing countries with low levels of urbanization, a large portion of the population is rural and resides in single-story buildings that lack adequate protection from earthquake-induced shaking. When a large earthquake event occurs, the collapse of these buildings led to a large number of deaths (Jackson 2006). This is especially true for rural households in semiarid and arid climates (for example, Iran and Pakistan), where the main source of building materials are mud bricks and shelters have heavy clay roofing (Jackson 2006;Spence 2007;Motamedi et al. 2012). This may also explain why the estimated proportion of earthquake fatalities in rural areas to total fatalities were over 80% from 2003 to 2018; such was the case resulting from the 2005 Kashmir and 2015 Nepal earthquakes (Wyss 2018). In addition, countries with low levels of urbanization typically struggle with poverty and corruption, and the earthquake resistance of buildings is largely overlooked, especially with the expansion of informal settlements (Green 2008;Ambraseys and Bilham 2011;Bilham and Gaur 2013).
With the movement of population from rural areas to urban areas with a high population density (that is, urbanization), the exposure rate to earthquakes changes (Jackson 2006;Doocy et al. 2013b;Bai et al. 2017). Furthermore, most developing countries with earthquakeprone areas lack or have weak enforcement of building codes for earthquakes (Daniell et al. 2014), which affects population vulnerability. Among the 94 countries having earthquake-prone areas over the last 28 years, more than 70% of total deaths occurred in countries with a vulnerability code of 4 or 5 in the PAGER system and an urbanization ratio between 20 and 60% (Fig. 5). If an earthquake with a large magnitude strikes a rapidly urbanizing area with a weak building resistance, the potential exists for a large number of fatalities . For example, substantially more deaths occurred in urban areas than rural areas from both the 1976 Tangshan Earthquake in China (242 thousand deaths) and the 31thousand deaths of the 2003 Bam Earthquake in Iran (EM-DAT 2018). In other words, among countries with a middle level of urbanization, both the collapse of fragile single-floor rural settlements as well as unreinforced multifloor urban apartments contribute to earthquake fatalities.
A third conceptual urbanization stage occurs when earthquake fatalities substantially decline, likely due to the enactment and strong enforcement of building codes (Crowley and Elliott 2012) and increased economic prosperity that situates more of the population in higher quality structures (Jackson 2006). This is especially true for highly developed countries, such as the United States, New Zealand, and Japan. The data support this transition with only a small portion of total deaths (1.3-1.7%) in global earthquake-prone areas during 1990-2018 occurring in the United States, New Zealand, and Japan (Fig. 5). Previous studies also supported this finding (Ambraseys and Bilham 2011;He et al. 2018). For example, Spence (2007) found that the annual fatalities caused by earthquakes to the total population in the United States and Japan decreased from 41 and 2670 persons per million to 3.4 and 230 victims, respectively. In contrast, the fatality number in Iran increased from 308 persons per million pre-1960 to 2970 since the 1980s because of the weak enforcement of building codes. Other researchers further support the argument that urbanization provides a reduction in earthquake mortality risk by improving: (1) the technology for sending earthquake alerts (Suárez et al. 2009); (2) knowledge about mitigating seismic risk (Takeuchi and Shaw 2014); and (3) access to shelter during earthquakes (Alçada-Almeida et al. 2009).

Implications of Future Urbanization in Global
Earthquake-Prone Areas Population and urban land area will likely continue to increase in global earthquake-prone areas in the future (Table 3). According to the projections of population and urban land area under the five Shared Socioeconomic Pathways (SSPs) examined by O'Neill et al. (2014O'Neill et al. ( , 2017, the population in global earthquake-prone areas will increase by 376.23 (215.05-593.87) million from 2015 to ð Þ ¼ 44% (42-46%, 95% confidence interval). The shaded area represents the 95% confidence interval 2050, with a rate of growth of 24.9% (14-39.4%). Meanwhile, the increase in urban land area in global earthquakeprone areas will likely be 60,991 (47,761-73,522) km 2 , with a rate of growth of 65.0% (50.9-78.3%). Population density in urban areas in global earthquake-prone areas will decrease along with the increase in urbanization.
More than half of the total increase in population (51.5%) in global earthquake-prone areas will occur in nations that will have an urbanization ratio between 40 and 60% (Table 3). In particular, Pakistan, India, Afghanistan, and Bangladesh will rank in the top four. The total population growth in these four countries' earthquake-prone areas will reach 181.62 (92.62-300.05) million, accounting for 48.3% (43.1-50.5%) of total population growth in global earthquake-prone areas. This is consistent with previous findings that most of the urban agglomeration growth at risk will be in Asian developing countries (Bilham 1999(Bilham , 2009. Our study reveals that these countries are mostly expected to pass the 44% turning point in the urbanization ratio by 2050, which implies that adequate strategies to improve building resistance to earthquakes should be adopted more rapidly than has been done in the past. For the next three decades in these countries, a window exists in which to provide a safer living environment for local residents (Ambraseys and Bilham 2011;Uitto and Shaw 2016). Among these developing countries, increased urbanization and population growth is typically accompanied by the expansion of informal settlements and multistory buildings, which are especially vulnerable to catastrophic natural hazards and disasters (Jaiswal et al. 2011a;Henderson et al. 2016). Thus, mitigating the seismic risk of the urban poor and enhancing the seismic resilience *The proportion refers to the change of population and urban land in each group of countries to the total change in global earthquake-prone areas.
UR stands for the urbanization ratio, percentage or urban population to total population. The values in bracket are the minimum and maximum values among the five Shared Socioeconomic Pathways of high-rise buildings should be prioritized to decrease seismic risk in developing countries. Approximately 45% of urban land growth in global earthquake-prone areas will take place in countries with an urbanization ratio above 80% in 2050 (Table 3). This could be good news for global seismic risk mitigation, as historical urban expansion in the developed countries is characterized by reinforced buildings for earthquake resistance and sprawling development that disperses impact (Spence 2007). Several developed countries have successful historical experiences dealing with seismic risk and mitigating fatalities during large earthquakes (Spence 2007;Crowley and Elliott 2012). If careful and rigorous disaster planning and operationalization of policy to ensure code enforcement for seismic resistance and retrofit of buildings are upheld (Bilham 2009;Samaddar et al. 2017), policy initiative will lead to further success stories. Enhancing the resilience of buildings, infrastructure, and urban systems in rapidly urbanizing cities would increase the capacity to cope with catastrophic earthquakes by medium and small urban centers.

Limitations and Future Perspectives
We analyzed the spatiotemporal dynamics of population and urban land exposed to earthquake hazard globally, and investigated the association between urbanization and earthquake fatalities. The results suggest that as the urbanization ratio increases from low to moderate values, earthquake-induced fatalities also increase. However, as the urbanization ratio increases from moderate to higher levels (that is, 40-60%), fatalities decreased. Under data limitations for our global analysis, the range of uncertainty is highest at the ends of the urbanization ratio (that is, low and high; Fig. 7). But with inclusion of variation in earthquake fatality estimates the presented evidence supports a Kuznet's curve type of response in fatalities due to an increasing urbanization ratio. While we present an analysis of the relationship between urbanization and earthquake fatalities globally, the limited time span of 25 years and coarse spatial resolution of population and urban land data may not be sufficient to form a definitive conclusion. We hope it will provide a point of departure for further studies that use newly available population and urban land datasets on various scales (Merkens et al. 2016;Boke-Olén et al. 2017;Chen et al. 2020a;Chen et al. 2020b;Chen et al. 2020c).
Furthermore, some factors that can affect earthquake fatalities, such as age and gender distribution were omitted due to a lack of data availability at the global scale. Given the disparate effect of catastrophes on different demographics (Vink and Takeuchi 2013), the geographic heterogeneity of demography could provide additional insight into at-risk populations or guide mitigations efforts targeted at specific demographic segments of the population (Jackson 2006).
Lastly, while it is likely that poor building resistance to earthquakes in developing countries is a dominant driver of the observed inverted U-shaped relationship between urbanization and earthquake fatalities at the global scale, new national or global data on the spatial distribution of informal settlements and poorly constructed buildings are required to understand the contribution of informal settlements relative to other factors impacting earthquake fatalities.

Conclusions
Over the past 25 years, approximately 67.8% of population growth and 74.5% of the total deaths, or 307,918 (282,938-405,431) deaths in global earthquake-prone areas took place in developing countries with an urbanization ratio between 20 and 60% in 2015. In contrast, only 30.5% of urban land expansion in global earthquake-prone areas occurred in these countries. The relation between the urbanization ratio and fatalities exhibited an inverted U-shape, with a turning point at approximately 44% (42-46%, 95% confidence interval), after considering other factors such as hazard, exposure, and vulnerability.
By 2050, more than 50% of the total increase in population in global earthquake-prone areas will still take place in developing countries with an urbanization ratio between 40 and 60% in 2050. Among the 94 countries having earthquake-prone areas, 16 of them will pass the urbanization turning point (44%) before 2050. Special attention to building resistance to earthquake fatalities should be initiated in developing countries (Haiti, Pakistan, India, Afghanistan, and Bangladesh) in the future to mitigate earthquake risk.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.