arXiv:1901.05086v1 [physics.soc-ph] 15 Jan 2019 aos[7.Rcado’ oko ete forecasting weather on work Richardson’s turbu- world [27]. make of later famous would equations Lorenz Norton the Edward that in very lence same Richardson lurking the way, chaos earlier, century this mathematical half In a un- almost discovered, developed later. nearly be decades not several would maketil to predictions power weather computing sufficient useful [38], that weather fact the the forecast the despite to pioneered approaches numerical also of Richardson use called number. is Richardson systems turbulent the related in flows quantity the shear and dimensionless develop buoyancy recognition to a to in contributions, aimed and those here forecasting, of many work weather for his of focus of mathematics primary Much his was decades. which meteorology, in wars. of statistics the a in of pattern discovery “scale-free” his foreshadowed also of [30]. insight systems studies Richardson’s technological numerous and biological, informed social, fractal has complex on which work the [28], celebrated paved geometry Mandelbrot’s effect, Benoit Richardson length. for the total way called it—the coastline’s now measure the effect, to longer This used the ruler ruler, the the of paradoxi- shorter length depends, coastline the Richardson a on of cally, coastline? length British the the that simple showed is deceptively long a how by captured question: is which paradox,” many line its in humans of behavior collective forms. galax- the of and formation ies, the biology, as developmental systems understand as to varied con- efforts which scientific science, inform complexity to in of tinue pieces roles major critical deadly played other other two fre- also and Richardson the wars of [39–41]. interstate analysis conflicts of One his severity was and how contingency. works quency celebrated and of most rules chaos his simple the of of from interaction complex- emerges how the structure sci- both from complexity understand arises of to ity aims field which the modern [30], of the ence one of as fathers stands founding (1881–1953) Richardson Fry Lewis h eodaoefo ihrsnspoern work pioneering Richardson’s from arose second The “coast- the on work his in arose these of first The ehnssta rdc niiulwr n eid fpeac ou of on conflict. periods constraints in and important trends wars about place individual debate facts produce These that mechanisms ori caveats. Richardson’s en that few mid-20th have show we a the conflicts data, violent in and methods data about statistical and insights the methods his and his comprehensive, process Although arrival Poisson distribution. fol simple between, a in followed everything frequency and wars interstate to homicides ei r ihrsnage httefeunyadseverity and frequency the that argued Richardson Fry Lewis 1 eateto optrSine nvriyo ooao Bo Colorado, of University Science, Computer of Department 2 iFotesIsiue nvriyo ooao Boulder, Colorado, of University Institute, BioFrontiers ntefeunyadsvrt fitrtt wars interstate of severity and frequency the On 3 at eIsiue at e M851 USA 87501, NM Fe, Santa Institute, Fe Santa ao Clauset Aaron Pr( htgnrt asadproso ec 1 3 3 5 35, 25, 23, 23, [1, peace of periods mechanisms and dy- wars political it generate and long-term that social Richardson, the underlying the on by of constraints namics claimed strong patterns indicate would simple the follow new time. over generating change of not inter- rules do is of the wars model statistics that Richardson’s meaning the Crucially, “stationary,” for worldwide. onsets—represent model wars their testable state for and process war simple Poisson for a a distribution and power-law comprehen- sizes patterns—a less far these were data sive, his and standards modern between years for distribution war [39]. geometric new wars a simple for a probability im- annual and process, constant Poisson a that simple both a argued plying followed he wars of Second, timing the parameter. “scaling” the called in nwihtepoaiiyta a kills war a distribu- wars that power-law these probability a the of which called in “severities,” pattern, tion, or precise sizes, a the followed First, events, that pattern. these argued statistical of their he basis about claims the two On made he [40]. 1820–1945 from wars this In statis- 34]. wars. the 22, interstate about of 11, ideas civil tics Richardson’s [9, [4], on wars focus insurgency interstate we chapter, and [16], 24], terrorism includ- [9, [3], conflict, wars political unrest of civil like study have ing events the events small-scale to larger central about about become claims ideas his Richardson’s to homicides, attention paid uni- little now Although followed is wars, [41]. patterns interstate homicides statistical like like versal events events small-scale large-scale “quar- from deadly to kinds, of all severity of and he rels” frequency which the in that conflict, violent argued of studies foundational son’s statistics long-term the wars. in of interest his foreshadowed also ,2 3, 2, 1, fteeprclsaitc fitrtt asral do really wars interstate of statistics empirical the If by rigorous not were methods statistical his Although ihrsnsoiia nlssol oee interstate covered only analysis original Richardson’s Richard- in together came threads intellectual These x ) oe nvra ttsia atrs their patterns: statistical universal lowed ∝ ia lisapa agl orc,with correct, largely appear claims ginal ∗ rsvrt olwdasml power-law simple a followed severity ir ue.I hscatr sn modern using chapter, this In dured. x − ,adse ih ntepersistent the on light shed and e, fdal oflcso l id,from kinds, all of conflicts deadly of nesadn fteunderlying the of understanding r α etr eenihrrgru nor rigorous neither were century o all for , le,C 00,USA 80309, CO ulder, O833 USA 80303, CO x ≥ x min > ,adwhere and 0, x epeis people > α is 1 2

107

106

105

104

3

Severity (battle deaths) 10

102 1823 1838 1853 1868 1883 1898 1913 1928 1943 1958 1973 1988 2003 Year

FIG. 1: The Correlates of War (CoW) interstate war data [44] as a conflict time series, showing both severity (battle deaths) and onset year for the 95 conflicts in the period 1823–2003. The absolute sizes of wars range from 1000, a minimum by definition, to 16,634,907, the recorded battle deaths of the Second World War. Delays between consecutive war onsets range from 0 to 18 years, and average 1.91 years. Most wars (79%) ended no more than 2 years after their onset. Originally published in Ref. [13].

47]. A long-running debate within the study of conflict and longest fluctuations occur for a scaling parameter has focused on whether or not such conflicts are char- of α < 3, when one or both the mean and variance are acterized by genuine trends. (See Ref. [19] for a recent mathematically infinite, i.e., they never converge, even review of this debate.) for infinite-sized samples. For interstate wars, this prop- If the underlying mechanisms that produce wars are erty could produce long transient patterns of low-severity stationary, then any “trend” is inherently illusory. How- or the absence of wars, making it difficult to distinguish ever, deciding whether trends exist has proved difficult a genuine trend toward peace from a mere fluctuation in to resolve, in part because there are multiple ways to an- a stationary process. swer this question, depending on what type of conflict To illustrate the counter-intuitive nature of power- is chosen, what variable is analyzed, and how the notion law distributions, consider a world where the heights of of change is formalized. Different choices can lead to Americans are power-law distributed, but with the same opposite conclusions about the existence or direction of average as reality (about 1.7 meters), and we line them change in the statistics of conflict [6, 12, 13, 19, 22, 33]. up in a random order. In this world, nearly 60,000 Amer- Here, we consider a more straightforward question: icans would be as tall as the tallest adult male on record given modern statistical tools and interstate war data, (2.72 meters), 10,000 individuals would be as tall as an do Richardson’s claims about statistical patterns hold adult male giraffe, one would be as tall as the Empire up, and if so, what does that imply about the long peace State Building (381 meters), and 180 million diminutive of the post-war period? For this investigation, we ap- individuals would stand only 17 cm tall. As we run down ply state-of-the-art methods [14, 15] to the set of inter- the line of people, we would repeatedly observe long runs state wars 1823–2003 given in the Correlates of War data of relatively short heights, one after another, and then, set [44] (Fig. 1). This data set provides comprehensive rarely, we would encounter a person so astoundingly tall coverage in this period, with few artifacts and relatively that their singular presence would dramatically shift our low measurement bias, and allows us to focus on a period estimate of the average or variance of all heights. This is during which Richardson’s model is plausible [11]. the kind of pattern that we see in the sizes of wars. Identifying a power-law distribution within an empiri- cal quantity can suggest the presence of exotic underly- PRELIMINARIES ing mechanisms, including nonlinearities, feedback loops, and network effects [31], although not always [36], and Before analyzing any data, we must clarify several epis- power laws are believed to occur broadly in complex so- temological issues and the impact of different assump- cial, technological, and biological systems [14]. For in- tions on the accuracy and interpretation of the analysis. stance, the intensities or sizes of many natural disas- Power-law distributions have unusual mathematical ters, such as earthquakes, forest fires, and floods [31], properties [31], which can require specialized statistical as well as many social disasters, like riots and terrorist tools to analyze. (For a primer on power-law distribu- attacks [3, 16], are well-described by power laws. tions in conflict, see Refs. [9, 16].) For instance, when Testing if some quantity does or does not follow a observations are generated by a power law, time series power law requires specialized statistical tools [14, 37], of summary statistics like the mean or variance can ex- because uncertainty tends to be greatest in the upper hibit long fluctuations resembling a trend. The largest tail, which governs the frequency of the largest and rarest 3 events, i.e., the frequency of large wars. Modern sta- to modeling the underlying processes that generate in- tistical tools provide rigorous methods for estimating terstate wars. This represents an important direction for and testing power-law models, distinguishing them from future work. other “heavy-tailed” distributions, and even using them to make statistical forecasts of future events [15]. Poisson processes pose fewer statistical issues than THE SIZE AND TIMING OF WARS power-law distributions. However, for consistency, our analysis applies similar methods to both war timing and In Richardson’s view, the size or severity of an inter- size data. Specifically, we estimate an ensemble of mod- state war (number of battle deaths) follows a power-law els, each fitted to a bootstrap of the empirical data, which distribution of the form Pr(x) ∝ x−α for some α> 1 and better represents our statistical uncertainty than would a for all x ≥ xmin > 0. single model. Technical details are described in Ref. [13]. Using standard techniques to estimate α and xmin, and Richardson’s models are defined in terms of absolute to test the fitted distribution [14], we find that the set numbers, i.e., the number of interstate wars per year and of observed interstate wars sizes, from 1823–2003, are the number of battle deaths per war. Hence, we consider statistically indistinguishable from an iid draw from a war variables in their unnormalized forms and consider power-law distribution (Fig. 2) [13]. all recorded interstate wars, meaning that our analysis Similarly, Richardson posits that the onset of a new takes the entire world as a system. interstate war follows a Poisson process, meaning that Analyses of interstate war statistics sometimes normal- wars arrive at a constant rate q, and the time t between ize either the number of wars or their sizes by some kind onsets of new wars follows a geometric distribution of the − of reference population, e.g., dividing a war’s size by the form Pr(t) ∝ e qt for some q> 0 and for all t ≥ 1. global population at the time. Such normalizations rep- Using the same techniques as above to estimate q resent additional theoretical assumptions about the un- and test the fitted distribution, we find that the set of derlying data generating process. observed delays between onsets are statistically indis- For instance, normalizing the number of wars per year tinguishable from an iid draw from a Poisson process by the number of pairs of nations that could be at war as- (Fig. 3) [13]. sumes that war is a dyadic event and that dyads indepen- That is, both the size and timing of interstate dently generate conflicts with equal probability [47]. This wars 1823–2003 are statistically indistinguishable from choice of normalizer grows quadratically with the num- Richardson’s simple model of a power-law distribution for ber of nations, and will create the appearance of a trend toward fewer wars, even if Richardson’s stationary model of wars is correct in an absolute sense. Considerable ev- 100 idence indicates that dyads do not independently gener- ate conflicts, and that dyadic likelihoods vary across time x and space, and by national covariates [1, 23, 25, 35, 47]. 0.25

Similarly, normalizing a war’s size by the nations’ or Density x world’s total population, producing a per-capita rate [12, 0.50 34], assumes that individuals contribute independently 1.2 1.4 1.6 1.8 -1 x Power-law exponent, α and with equal probability to potential or actual vio- 10 0.75 lence, regardless of who or where they are. Normalizing by world population is thus equivalent to assuming that doubling Canada’s population would linearly increase the level of violence in a war in Yemen. In general, human populations have increased so dramatically over the past Fraction of wars with at least x deaths

200 years that this normalizer nearly always produces the 10-2 appearance of a decline in violence, even if wars were, in 103 104 105 106 107 108 an absolute sense, stationary. However, there is little ev- Battle deaths, x idence that real conflict sizes or rates increase linearly with population [5, 17, 32]. FIG. 2: Interstate wars sizes, 1823–2003. The maxi- mum likelihood power-law model of the largest-severity wars That said, per-capita variables can be useful for other (solid line,α ˆ = 1.53 ± 0.07 for x ≥ xˆmin = 7061) is a statisti- reasons [34], and population surely does play some role in cally plausible data-generating process of these 51 empirical the sizes of wars, albeit probably not a simple one [17, 32]. severities (Monte Carlo, pKS = 0.78 ± 0.03). For reference, A realistic per-capita normalizer should instead account distribution quartiles are marked by vertical dashed lines. In- for the effects of alliances, geographic proximity, geopo- set: bootstrap distribution of maximum likelihood parameters litical stability, democratic governance, economic ties, Pr(ˆα), with the empirical value (black line). Originally pub- etc. [11], in addition to population, and would be akin lished in Ref. [13]. 4

100 Wars, x 90 Wars, x≥ x 0.25 Wars, x≥ x 0.50 Wars, x≥ x 0.75 Density 70

-1 0.3 0.4 0.5 0.6 10 Geometric exponent, q 50

Cum. number of wars 30 Fraction lasting at least T years

-2 10 10

0 5 10 15 20 1823 1853 1883 1913 1943 1973 2003 Years between war onsets, t Year

FIG. 3: Times between interstate war onsets, 1823– FIG. 4: Historical and simulated accumulation curves 2003. The maximum likelihood geometric model (solid line, of interstate wars. Empirical cumulative counts of wars qˆ = 0.428 ± 0.002 for t ≥ 1) is a statistically plausible of different sizes (dark lines) over time, alongside ensembles data-generating process of the empirical delays (Monte Carlo, of simulated counts from a stationary model (light lines), in pKS = 0.13 ± 0.01), implying that the apparent discontinuity which empirical severities are replaced iid with a bootstrap at t = 5 is a statistical artifact. Inset: the bootstrap distribu- draw from the empirical severity distribution. Dashed lines tion of maximum likelihood parameters Pr(ˆq), with the em- mark the end of the Second World War and the end of the pirical estimate (black line). Originally published in Ref. [13]. Cold War. Originally published in Ref. [13].

lation during the long peace is statistically unusual under war sizes and a Poisson process for their arrival [39, 41]. the reference distribution, it would indicate support for This agreement is remarkable considering the overall sim- an underlying change in the generating processes. plicity of the model compared to the complexity and con- We define “large” wars as those in the upper quartile of tingency of international relations over this period, and the historical war size distribution, meaning x ≥ x0.75 = the fact that this time period includes nearly 60 years of 26, 625 battle deaths, but similar results are obtained for additional data over Richardson’s original analysis. other large thresholds. The initial 1823–1939 period con- tains 19 such large wars, for an arrival rate of one per 6.2 years, on average. The “great violence” pattern of 1914– ARE LARGE WARS DECLINING? 1939, which spans the onsets of the First and Second World Wars, includes 10 large wars, or about one every Combining the two parts of Richardson’s model allows 2.7 years. The long peace of the 1940–2003 post-war pe- us to generate simulated interstate war data sets, drawn riod contains only 5 large wars, or about one every 12.8 from a stationary process, with similar onset times and years. Figure 4 shows the historical accumulation curves war sizes as the historical record. The statistics of these of these events, and for smaller wars, as a function of simulated histories define a reference distribution against time. which we can compare aspects of the historical record. Our combined model takes the historical war onset We now apply this model to address a long-running years as given, and then for each of these 95 conflicts, debate in international relations: did the underlying pro- draws a synthetic war size iid from the empirical size dis- cesses that generate interstate wars change after the Sec- tribution (with replacement), as in a simple bootstrap of ond World War? This point in history is commonly the data. Ref. [13] considers two additional models of proposed as the beginning of a “long peace” pattern this flavor, which produce larger variances in the accu- in interstate conflict, meaning a pronounced decrease mulation curves for large wars but yield similar results in the frequency and severity of wars, especially large and conclusions. ones [6, 18, 26, 34, 47]. To test the long peace hypothesis from Richardson’s perspective, we consider the accumulation of large inter- Evaluating the past state wars over time, and assess whether or when that ac- cumulation in the post-war period represents a low prob- Within the historical accumulation curve for large ability event under the stationary model. If the accumu- wars, the long peace is a visible pattern, in which the 5

Empirical pattern Formalization Model great violence Pr(V ≡ n ≥ 10 large wars over t ≤ 27 years) 0.107(1) long peace Pr(P ≡ n ≤ 5 large wars over t ≥ 64 years) 0.622(2)

TABLE I: Stationary likelihood of empirical conflict patterns. Under a simple stationary model of conflict generation (see text), the estimated likelihoods of observing two particular large-war patterns over the period 1823–2003: a “great violence,” meaning 10 or more large war onsets (x ≥ x0.75) over a 27 year period (the empirical count of such onsets, 1914–1939); or, a “long peace,” meaning 5 or fewer large war onsets over a 64 year period (the empirical count of such onsets, 1940–2003). Probabilities estimated by Monte Carlo. Parenthetical values indicate the standard error of the least significant digit. Originally published in Ref. [13].

arrival rate (the curve’s slope) is substantially more flat past future than in the preceding great violence period (Fig. 4). How- 170 Wars, x Wars, x≥ x ever, under the stationary model, this pattern is well 0.25 150 Wars, x≥ x within the envelope of simulated curves, and the observed 0.50 Wars, x≥ x pattern is statistically indistinguishable from a typical 130 0.75 excursion, given the heavy-tailed nature of historical war sizes. 110 In fact, most simulated war sequences contain a period 90 of “peace” at least as long in years and at least a peace- 70 ful in large-war counts as the long peace (Table I). Fifty years or more of relatively few large wars is thus entirely 50 typical, given the empirical distribution of war sizes, and Cumulative number of wars 30 observing a long period of peace is not necessarily ev- idence of a changing likelihood for large wars [12, 13]. 10 Even periods comparable to the great violence of the 1823 1863 1903 1943 1983 2023 2063 2103 World Wars are not statistically rare under Richardson’s Year model (Table I). Taking these findings at face value implies that the FIG. 5: How long must the peace last? Simulated ac- probability of a very large war is constant. Under the cumulation curves for wars of different sizes under a simple stationary model, overlaid by the empirical curves up to 2003 model, the 100-year probability of at least one war with (dark lines) and linear extrapolations of the empirical post- 16, 634, 907 or more battle deaths (the size of the Second war trends (the long peace) for the next 100 years (dashed World War) is 0.43 ± 0.01, implying about one such war lines). Originally published in Ref. [13]. per 161 years, on average.

In this extrapolated future, the post-war pattern of Peering into the future relatively few large wars becomes progressively more un- likely under a stationary hypothesis (Fig. 5). However, We can also use Richardson’s model to simulate future it is not until 100 years into the future that the long sequences of interstate wars, and thereby evaluate how peace becomes statistically distinguishable from a large long the long peace must last before it becomes com- but random fluctuation in an otherwise stationary pro- pelling evidence that the underlying processes did, in- cess. Even if there were no large wars anywhere in the deed, change after the Second World War. To extend world after 2003, the year of significance would arrive the simulated war sequences beyond 2003, for each year, only a few decades earlier. we create a new war onset according to a simple Bernoulli The consistency of the historical record of interstate process, with the historical production rate (on average, wars with Richardson’s stationary model places an im- a new war every 1.91 years). We then linearly extrap- plicit upper bound on the magnitude of change in the olate the long peace pattern, in which a new large war underlying conflict generating process since the end of occurs on average every 12.8 years, out into the future the Second World War [8]. Our modeling effort here can- until 95% of the simulated accumulation curves exceed not rule out the existence of a change in the rules that the extrapolated pattern’s curve. At that moment in generate interstate conflicts, but if it occurred, it cannot time, the long peace pattern will have become statisti- have been a dramatic shift. The results here are entirely cally significant, by conventional standards, relative to consistent with other evidence of genuine changes in the the stationary model, and we could say with confidence international system, but they constrain the extent to that the time since the Second World War was governed which such changes could have genuinely impacted the by a different, more peaceful underlying process. global production of interstate wars. 6

DISCUSSION also include unknown mechanisms. In the long run, processes that promote interstate war The agreement between the historical record of inter- may be consequences of those that reduce it over the state wars and Richardson’s simple model of their fre- shorter term, through feedback loops, tradeoffs, or back- quency and severity is truly remarkable, and it stands as lash effects. For example, the persistent appeal of na- a testament to Richardson’s lasting contribution to the tionalism, whose spread can increase the risk of interstate study of violent political conflict. war [45], is not independent of deepening economic ties via globalization [46]. Investigating such interactions is There are, however, a number of caveats, insights, and a vital direction of future research, and will facilitate a questions that come out of our analysis. For instance, more complete understanding of the processes that gov- Richardson’s Law—a power-law distribution in conflict ern the likelihood of patterns like the long peace. event sizes—appears to hold only for sufficiently large “deadly quarrels,” specifically those with 7061 or more More concretely, our results here indicate that the battle deaths. The lower portion of the distribution is post-war efforts to reduce the likelihood of large inter- slightly more curved than expected for a simple power state wars have not yet changed the observed statistics law, which suggests potential differences in the processes enough to tell if they are working. This fact does not that generate wars above and below this threshold. lessen the face-value achievement of the long peace, as a large war today between major powers could be very With only 95 conflicts and a heavy-tailed distribution large indeed, and there are real benefits beyond lives of war sizes, there are relatively few large wars to con- saved [42] that have come from increased economic ties, sider. This modest sample size surely lowers the statis- peace-time alliances, and the spread of democracy. How- tical power of any test and is likely partly to blame for ever, it does highlight the continued relevance of Richard- needing nearly 100 more years to know whether the long son’s foundational ideas as the appropriate null hypoth- peace pattern is more than a run of good luck under a esis for patterns in interstate war. stationary process. One explanation for the apparent stationarity of wars One could imagine increasing the sample size by in- since 1823 is the existence of compensatory trends in re- cluding civil wars, which are about three times more lated conflict variables that mask a genuine change in numerous than interstate wars over 1823–2003. Includ- the conflict generating processes. Patterns across mul- ing these, however, would confound the resulting inter- tiple conflict variables do seem to indicate a broad shift pretation, because civil wars have different underlying toward less violence [20, 21, 33, 34]. But, not all conflict causes [10, 43, 48], and because the distribution of civil variables support this conclusion, and some, such as mil- war sizes is shifted toward smaller conflicts and exhibits itary disputes and the frequency of terrorism, seem to be relatively fewer large ones [24]. increasing, instead [15, 22]. Untangling conflict variables’ Putting aside these technical issues, the larger ques- interactions and characterizing their trends and differ- tion our analysis presents is this: how can it be possible ences across groups of nations will be a valuable line of that the frequency and severity of interstate wars are future work. so consistent with a stationary model, despite the enor- An alternative explanation is that the mechanisms mous changes and obviously non-stationary dynamics in that govern the likelihood of war have unfolded heteroge- human population, in the number of recognized states, neously across time and geographic regions over the past in commerce, communication, public health, and tech- 200 years, thereby creating an illusion of global station- nology, and even in the modes of war itself? The fact ary by coincidence. The long peace pattern is sometimes that the absolute number and sizes of wars are plausibly described only in terms of peace among largely European stable in the face of these changes is a profound mystery powers, who fell into a peaceful configuration after the for which we have no explanation. great violence for well understood reasons. In parallel, There is, of course, substantial evidence for a genuine however, conflicts in other parts of the world, most no- post-war trend toward peace, based on mechanisms that tably Africa, the Middle East, and Southeast Asia, have reduce the likelihood of war [23, 25, 35] and on statistical became more common, and these may have statistically signatures of a broad and centuries-long decline in gen- balanced the books globally against the decrease in fre- eral violence [20, 21, 33, 34] or the improvement of other quency in the West, and may even be causally dependent aspects of human welfare [42]. on the drivers of European war and then peace. But, a full accounting of the likelihood that the long Developing a more mechanistic understanding of how peace will endure must also consider mechanisms that changes in the likelihood of conflict in one part of the increase the likelihood of war (for example, Refs. [2, 7, world may induce compensatory changes in the likeli- 29]). War-promoting mechanisms certainly include the hood of conflict in other parts of the world would have reverse of established peace-promoting mechanisms, e.g., enormous value if it could explain how some regions can the unraveling of alliances, the slide of democracies into fall into more peaceful patterns as a group, while other autocracy, and the fraying of economic ties, but they may regions go the opposite direction. The evident stability 7 of Richardson’s Law may, in the end, be an artifact of [11] Cederman, L.-E., T. C. Warren, and D. Sornette (2011). these kinds of complex, macro-scale dynamics, playing Testing Clausewitz: Nationalism, mass mobilization, and out across the global stage. the severity of war. International Organization 65, 605– 638. Finally, it is worth reiterating how remarkable and [12] Cirillo, P. and N. N. Taleb (2015). On the statistical counter-intuitive it is that Richardson’s original mod- properties and tail risk of violent conflicts. Physica A 429, els of the frequency and severity of wars, first proposed 252–260. more than half a century ago, successfully hold up un- [13] Clauset, A. (2018). Trends and fluctuations in the sever- der modern, more rigorous statistical methods of evalua- ity of interstate wars. Science Advances 4, eaao3580. tion applied to far more comprehensive data. A sobering [14] Clauset, A., C. R. Shalizi, and M. E. J. Newman (2009). implication of this success is that the probability of a Power-law distributions in empirical data. SIAM Re- view 51, 661–703. large interstate appears to have remained constant, de- [15] Clauset, A. and R. Woodard (2013). Estimating the his- spite profound collective efforts to lower it. torical and future probabilities of large terrorist events. More importantly for the general study of conflict, Annals of Applied Statistics 7, 1838–1865. Richardson’s work presents a simple and enduring mys- [16] Clauset, A., M. Young, and K. S. Gleditsch (2007). On tery: how can the frequency and severity of interstate the frequency of severe terrorist events. Journal of Con- flict Resolution 51 wars be so consistent with a stationary model, despite (1), 58–87. [17] Falk, D. and C. Hildebolt (2017). Annual war deaths the dramatic changes and obviously non-stationary dy- in small-scale versus state societies scale with population namics in so many other aspects of human civilization? size rather than violence. Current Anthropology 58 (6), Answering this question will shed new light on the under- 805–813. lying causes of war, and greatly inform efforts to promote [18] Gaddis, J. L. (1986). The long peace: elements of sta- peace. Richardson, who was an avowed pacifist and who bility in the postwar international system. International worked as an ambulance driver during the First World Security 10 (4), 99–142. War, would surely be pleased if his work on the statistics [19] Gleditsch, K. S. and A. Clauset (2018). Trends in con- flict: What do we know and what can we know? In W. C. of war ultimately helped devise better policies to promote Wohlforth and A. Gheciu (Eds.), Handbook of Interna- peace. tional Security, pp. 227–244. Oxford University Press. [20] Goldstein, J. S. (2011). Winning the War on War. Hialeah, FL: Dutton. [21] Gurr, T. R. (2000). Ethnic warfare on the wane. Foreign Affairs 79, 52–64. ∗ Electronic address: [email protected] [22] Harrison, M. and N. Wolf (2012). The frequency of wars. [1] Alesina, A. and E. Spolaore (1997). On the number and Economic History Review 65 (3), 1055–1076. size of nations. Quarterly Journal of Economics 112, [23] Jackson, M. O. and S. M. Nei (2015). Networks of mili- 1027–1056. tary alliances, wars, and international trade. Proc. Natl. [2] Barbieri, K. (1996). Economic interdependence: A path Acad. Sci. USA 112 (50), 15277–15284. to peace or a source of interstate conflict? Journal of [24] Lacina, B. (2006). Explaining the severity of civil wars. Peace Research 33 (1), 29–49. Journal of Conflict Resolution 50 (2), 276–289. [3] Biggs, M. (2005). Strikes as forest fires: Chicago and [25] Leeds, B. A. (2003). Do alliances deter aggression? The Paris in the late 19th century. American Journal of So- influence of military alliances on the initiation of milita- ciology 111 (1), 1684–1714. rized interstate disputes. American Journal of Political [4] Bohorquez, J. C., S. Gourley, A. R. Dixon, M. Spagat, Science 47, 427–439. and N. F. Johnson (2009). Common ecology quantifies [26] Levy, J. S. and W. R. Thompson (2010). Causes of War. human insurgency. Nature 462, 911–914. West Sussex, UK: Wiley-Blackwell. [5] Bowles, S. (2009). Did warfare among ancestral hunter- [27] Lorenz, E. N. (1963). Deterministic nonperiodic flow. gatherers affect the evolution of human social behaviors? Journal of the Atmospheric Sciences 20 (2), 130–141. Science 324, 1293–1298. [28] Mandelbrot, B. (1967). How long is the coast of Britain? [6] Braumoeller, B. F. (2013). Is war disappearing? Statistical self-similarity and fractional dimension. Sci- APSA Chicago 2013 Meeting. Available at SSRN: ence 156, 3775. https://ssrn.com/abstract=2317269. [29] Mansfield, E. D. and J. Snyder (1995). Democratization [7] Bremer, S. A. (1992). Dangerous dyads: Conditions af- and the danger of war. International Security 20 (1), 5– fecting the likelihood of interstate war, 1816-1965. Jour- 38. nal of Conflict Research 36 (2), 309—341. [30] Mitchell, M. (2011). Complexity: A Guided Tour. Oxford [8] Cederman, L.-E. (2001). Back to Kant: Reinterpreting University Press. the democratic peace as a macrohistorical learning pro- [31] Newman, M. E. J. (2005). Power laws, Pareto distri- cess. American Political Science Review 95 (1), 15–31. butions and Zipf’s law. Contemporary Physics 46 (5), [9] Cederman, L.-E. (2003). Modeling the size of wars: From 323–351. billiard balls to sandpiles. American Political Science [32] Oka, R. C., M. Kissel, M. Golitko, S. G. Sheridan, N. C. Review 97 (1), 135–150. Kim, and A. Fuentes (2017). Population is the main [10] Cederman, L.-E., K. S. Gleditsch, and H. Buhaug (2013). driver of war group size and conflict casualties. Proc. Inequality, grievances, and civil war. Cambridge: Cam- Natl. Acad. Sci. USA 114 (52), E11101–E11110. bridge University Press. [33] Payne, J. L. (2004). A History of Force: Exploring the 8

Worldwide Movement against Habits of Coercion, Blood- [42] Roser, M., E. O. Ospina, and J. Mispy (2017). Our World shed, and Mayhem. Sandpoint, ID: Lytton. in Data. https://ourworldindata.org. [34] Pinker, S. (2012). Better Angels of Our Nature. Penguin [43] Salehyan, I. and K. S. Gleditsch (2006). Refugees and Books. the spread of civil war. International Organization 60, [35] Ray, J. L. (1998). Does democracy cause peace? Annual 335–366. Review of Political Science 1, 27–46. [44] Sarkees, M. R. and F. Wayman (2010). Resort to War: [36] Reed, W. J. and B. D. Hughes (2002). From gene families 1816 – 2007. Washington, DC: CQ Press. and genera to incomes and internet file sizes: Why power [45] Schrock-Jacobson, G. (2012). The violent consequences laws are so common in nature. Physical Review E 66, of the nation: Nationalism and the initiation of interstate 067103. war. Journal of Conflict Resolution 56 (5), 825–852. [37] Resnick, S. I. (2006). Heavy-Tail Phenomena: Probabilis- [46] Smith, A. D. (1992). National identity and the idea of tic and Statistical Modeling. New York: Springer-Verlag. European unity. International Affairs 68 (1), 55–76. [38] Richardson, L. F. (1922). Weather Prediction by Numer- [47] Ward, M. D., R. M. Siverson, and X. Cau (2007). Dis- ical Process. Cambridge: Cambridge University Press. putes, democracies, and dependencies: A reexamination [39] Richardson, L. F. (1944). The distribution of wars in of the Kantian peace. American Journal of Political Sci- time. Journal of the Royal Statistical Society 57, 242– ence 51 (3), 583–601. 250. [48] Wucherpfennig, J., P. Hunziker, and L.-E. Cederman [40] Richardson, L. F. (1948). Variation of the frequency of (2016). Who inherits the state? Colonial rule and fatal quarrels with magnitude. Journal of the American postcolonial conflict. American Journal of Political Sci- Statistical Association 43, 523–546. ence. 60 (4), 882–898. [41] Richardson, L. F. (1960). Statistics of Deadly Quarrels. Pittsburgh: The Boxwood Press.