1. Introduction
Archaeological evidence suggests that Homo sapiens emerged in Africa as early as about 300,000 years ago.Footnote 1 For most of its history, humans were hunter-gatherers. Then, the Neolithic Revolution (the transition from hunting-gathering to agriculture) occurred in the Fertile Crescent about 12,000 years ago and then in other parts of the world.Footnote 2 In the late 18th and early 19th century, the Industrial Revolution (the transition from agriculture to the manufacturing of goods) took place in Britain and then in continental Europe and the USA.Footnote 3 Are these transitions in the economic evolution of human society inevitable? If not, what are the different conditions that could have potentially made the transitions more or less likely to occur?
This study develops a Malthusian model that captures the economic evolution of human society from hunting-gathering to agriculture and then from agriculture to an industrial economy. In our model, human society evolves across these stages as the size of the population grows. However, under endogenous population growth determined by the fertility decisions of optimizing agents in our microfounded model, the population may stop growing at any stage and never reach the next threshold. If it fails to reach the first threshold, then the population remains in a hunting-gathering Malthusian trap. If the population size reaches the first threshold, then an agricultural society emerges; therefore, both the Boserupian and Malthusian forces are present in our model.Footnote 4 The Neolithic Revolution occurs under the following conditions: a high level of agricultural productivity, a low cost of fertility, and a strong preference for fertility. In the main text, we discuss the intuition of these results and their relation to existing hypotheses. We also provide empirical evidence that high agricultural productivity leads to an earlier transition to agriculture.
After an agricultural society emerges, the economy eventually becomes completely agricultural until it reaches the next threshold. If it fails to reach the next threshold, the economy remains in an agricultural Malthusian trap and does not experience industrialization. Industrialization is influenced by the same conditions as the Neolithic Revolution (namely, a high level of agricultural productivity, a low cost of fertility, and a strong preference for fertility) and also other conditions: a high level of industrial productivity and a low fixed cost of industrial production. Therefore, the conditions (e.g., a high level of agricultural productivity) that trigger the Neolithic Revolution also trigger the subsequent industrialization, but not necessarily vice versa. Here, the importance of the population size on industrialization is due to its increasing returns to scale (i.e., having a large enough market to cover the fixed costs associated with industrial production) as in Murphy et al. (Reference Murphy, Shleifer and Vishny1989), whereas the importance of the population size on the Neolithic Revolution is due to the decreasing returns to scale in hunting-gathering as in North and Thomas (Reference North and Thomas1977). If the population reaches the industrial threshold, then a modern economy emerges and exhibits positive steady-state population growth.
Empirically, we use cross-country data to examine our key theoretical result that high agricultural productivity triggers not only the Neolithic Revolution but also the subsequent industrialization. We follow previous empirical studies, such as Olsson and Hibbs (Reference Olsson and Hibbs2005), Ashraf and Galor (Reference Ashraf and Galor2011), and Ang (Reference Ang2015), to consider an index of prehistoric biogeographic conditions that affect agricultural productivity and explore how it affects the transition to agriculture. Specifically, we use the index of biogeographic conditions as an instrument for the timing of transitions to agriculture and find that an earlier transition to agriculture has a positive effect on the degree of industrialization in the modern era.
This study relates to the literature on the economic modeling of the transition from hunting-gathering to agriculture,Footnote 5 see Smith (Reference Smith1975) and North and Thomas (Reference North and Thomas1977) for early studies and Weisdorf (Reference Weisdorf2005) for an excellent review of this literature.Footnote 6 A subsequent study by Locay (Reference Locay1989) develops a dynamic general equilibrium model with endogenous fertility to explore the transition of the human population from nomadic hunter-gatherers to a sedentary agricultural society, see also the interesting studies by Olsson (Reference Olsson2001) and Weisdorf (Reference Weisdorf2003).Footnote 7 Baker (Reference Baker2008) estimates an extended version of the Locay model using historical data on the incidence of agriculture and finds empirical support for the model. Weisdorf (Reference Weisdorf2011) explores the case in which the agricultural transition is caused by an exogenous discovery of agricultural technology.Footnote 8 Recent studies by Dow et al. (Reference Dow, Reed and Olewiler2009) and Dow and Reed (Reference Dow and Reed2015, Reference Dow and Reed2022) consider climate change as a cause of the transition to agriculture, whereas Bowles and Choi (Reference Bowles and Choi2019) explore the origins of private property as a reason for adopting agriculture.
Our model is based on Locay (Reference Locay1989) and Baker (Reference Baker2008) with the introduction of an industrial economy as the third stage of the economic evolutionary process, without which the population remains either in a hunting-gathering or an agricultural Malthusian trap in the long run. An important finding is that the transition from hunting-gathering to agriculture and the transition from agriculture to industrial production are both endogenous and share a similar set of determinants. For example, high agricultural productivity triggers not only the Neolithic Revolution but also the subsequent industrialization.Footnote 9 Ashraf and Galor (Reference Ashraf and Galor2011) provide empirical evidence that agricultural productivity has a significant positive effect on population density in the preindustrial era. Olsson and Hibbs (Reference Olsson and Hibbs2005) also provide evidence to show that prehistoric biogeographic conditions that are favorable to agriculture can trigger the Neolithic Revolution and the subsequent development in the industrial era. Using an index of biogeographic conditions as instruments, Ang (Reference Ang2015) finds that the timing of transitions to agriculture has a significant effect on technology adoption from 1000 BC to 1500 AD. Our empirical analysis follows the footsteps of these studies but instead examines the effects of agricultural productivity on the timing of transitions to agriculture and the degree of industrialization in the modern era.
This study also relates to the literature on unified growth theory, see Galor and Weil (Reference Galor and Weil2000) for the seminal study and Galor (Reference Galor2005, Reference Galor2011) for a comprehensive review. Studies in this literature explore the endogenous transition of an agricultural economy in a Malthusian trap to a modern industrial economy with technological progress and long-run economic growth. This study complements the interesting studies in this literature by developing a simple unified model that captures both the first transition from hunting-gathering to agriculture and the second transition from agriculture to a modern industrial economy. In the spirit of Diamond (Reference Diamond1997), Olsson and Hibbs (Reference Olsson and Hibbs2005) also model both of these important transitions in human history using a theoretical framework that focuses on the causal relationship between initial biogeographic conditions and the subsequent development of the economy. Specifically, they assume that a better biogeographic endowment causes a higher growth rate of productive knowledge, which in turn triggers the transitions once productive knowledge reaches certain exogenous thresholds. We take a more microfounded approach in which population growth is endogenously determined by optimizing agents and the transitions occur only when population size crosses thresholds that are also endogenously determined within the model. As a result, the transition from hunting-gathering to agriculture and the transition from agriculture to industry may not always occur depending on parameter conditions in our model.
The rest of this study is organized as follows. Section 2 presents the static model with exogenous population. Section 3 develops the dynamic model with endogenous population. Section 4 presents the empirical evidence. Section 5 concludes.
2. A static model of economic evolution
 Our model is based on Locay (Reference Locay1989) and Baker (Reference Baker2008). We extend the Locay model to introduce an industrial economy as the third stage of economic evolution. In the first stage, the population engages in hunting-gathering. In the second stage, an agricultural society emerges. In the third stage, an industrial economy emerges. The population consists of 
 $N$
 identical agents. Each agent is endowed with
$N$
 identical agents. Each agent is endowed with 
 $l$
 units of labor, which can be allocated to hunting-gathering
$l$
 units of labor, which can be allocated to hunting-gathering 
 $l_{H}$
, farming
$l_{H}$
, farming 
 $l_{F}$
, or industrial production
$l_{F}$
, or industrial production 
 $l_{Y}$
. Therefore, the labor constraint faced by each agent is
$l_{Y}$
. Therefore, the labor constraint faced by each agent is
 \begin{equation} l_{H}+l_{F}+l_{Y}=l\text{.} \end{equation}
\begin{equation} l_{H}+l_{F}+l_{Y}=l\text{.} \end{equation}
In the preindustrial era, industrial production does not yet exist, and hence, the constraint simplifies to 
 $l_{H}+l_{F}=l$
. There is also a fixed amount of land denoted as
$l_{H}+l_{F}=l$
. There is also a fixed amount of land denoted as 
 $Z$
, which can be used for hunting-gathering or farming.
$Z$
, which can be used for hunting-gathering or farming.
2.1. Hunting-gathering
 Hunting-gathering takes place in available land that is not occupied for farming. We use 
 $\overline{l}_{H}$
 to denote the average amount of labor endowment devoted to hunting-gathering. Then, total food production from hunting-gathering is given by
$\overline{l}_{H}$
 to denote the average amount of labor endowment devoted to hunting-gathering. Then, total food production from hunting-gathering is given by
 \begin{equation} H=\theta (\overline{l}_{H}N)^{\gamma }(Z_{H})^{1-\gamma }\text{,} \end{equation}
\begin{equation} H=\theta (\overline{l}_{H}N)^{\gamma }(Z_{H})^{1-\gamma }\text{,} \end{equation}
where 
 $\overline{l}_{H}N$
 and
$\overline{l}_{H}N$
 and 
 $Z_{H}\leq Z$
 are, respectively, the total amount of labor and land devoted to hunting-gathering. The parameters
$Z_{H}\leq Z$
 are, respectively, the total amount of labor and land devoted to hunting-gathering. The parameters 
 $\theta \gt 0$
 and
$\theta \gt 0$
 and 
 $\gamma \in (0,1)$
 measure, respectively, the productivity and labor intensity of the hunting-gathering process. An agent, who contributes
$\gamma \in (0,1)$
 measure, respectively, the productivity and labor intensity of the hunting-gathering process. An agent, who contributes 
 $l_{H}$
 units of labor to hunting-gathering, receives
$l_{H}$
 units of labor to hunting-gathering, receives 
 $h$
 units of food production given by
$h$
 units of food production given by
 \begin{equation} h=\frac{l_{H}}{\overline{l}_{H}N}\theta (\overline{l}_{H}N)^{\gamma }(Z_{H})^{1-\gamma }\text{,} \end{equation}
\begin{equation} h=\frac{l_{H}}{\overline{l}_{H}N}\theta (\overline{l}_{H}N)^{\gamma }(Z_{H})^{1-\gamma }\text{,} \end{equation}
in which the agent takes 
 $\overline{l}_{H}$
 and
$\overline{l}_{H}$
 and 
 $Z_{H}$
 as given.
$Z_{H}$
 as given.
2.2. Agriculture
 Farming also requires both labor and land. The farming production of an agent, who devotes 
 $l_{F}$
 units of labor to farming, is
$l_{F}$
 units of labor to farming, is
 \begin{equation} f=\varphi (l_{F})^{\alpha }z^{1-\alpha }\text{,} \end{equation}
\begin{equation} f=\varphi (l_{F})^{\alpha }z^{1-\alpha }\text{,} \end{equation}
where the parameters 
 $\varphi \gt 0$
 and
$\varphi \gt 0$
 and 
 $\alpha \in (0,1)$
 measure, respectively, the productivity and labor intensity in agriculture.
$\alpha \in (0,1)$
 measure, respectively, the productivity and labor intensity in agriculture. 
 $z$
 is the amount of land used by the agent. We follow Baker (Reference Baker2008) to assume a fixed ratio
$z$
 is the amount of land used by the agent. We follow Baker (Reference Baker2008) to assume a fixed ratio 
 $\rho$
 of land to farming labor given by
$\rho$
 of land to farming labor given by
 \begin{equation} z=\rho l_{F} \end{equation}
\begin{equation} z=\rho l_{F} \end{equation}
when agricultural land is not scarce (i.e., 
 $\rho \overline{l}_{F}N\lt Z$
); in this case,
$\rho \overline{l}_{F}N\lt Z$
); in this case, 
 $f=\varphi \rho ^{1-\alpha }l_{F}$
. Weisdorf (Reference Weisdorf2005) argues that the temporary constant returns to farming labor, which is also present in the analysis of North and Thomas (Reference North and Thomas1977), are a reasonable assumption when there is abundant agricultural land. When agricultural land becomes scarce, it is equally divided between agents, that is,
$f=\varphi \rho ^{1-\alpha }l_{F}$
. Weisdorf (Reference Weisdorf2005) argues that the temporary constant returns to farming labor, which is also present in the analysis of North and Thomas (Reference North and Thomas1977), are a reasonable assumption when there is abundant agricultural land. When agricultural land becomes scarce, it is equally divided between agents, that is,
 \begin{equation} z=Z/N\text{.} \end{equation}
\begin{equation} z=Z/N\text{.} \end{equation}
In this case, there is no more land available for hunting-gathering (i.e., 
 $Z_{H}=0$
), see North and Thomas (Reference North and Thomas1977) for a discussion that with communal property rights on agricultural land, farmers have better access to land than hunter-gatherers.
$Z_{H}=0$
), see North and Thomas (Reference North and Thomas1977) for a discussion that with communal property rights on agricultural land, farmers have better access to land than hunter-gatherers.
2.3. Industrial production
 As in Murphy et al. (Reference Murphy, Shleifer and Vishny1989), the operation of modern industrial production requires a fixed cost 
 $\delta \gt 0$
 under which total industrial output is given byFootnote 
10
$\delta \gt 0$
 under which total industrial output is given byFootnote 
10
 \begin{equation} Y=A(\overline{l}_{Y}N-\delta )\text{,} \end{equation}
\begin{equation} Y=A(\overline{l}_{Y}N-\delta )\text{,} \end{equation}
where 
 $\overline{l}_{Y}N$
 is the total amount of labor devoted to industrial production, and the parameter
$\overline{l}_{Y}N$
 is the total amount of labor devoted to industrial production, and the parameter 
 $A\gt 0$
 determines the level of industrial productivity. The fixed cost is shared by all agents when the industrial economy operates. Then, the output of industrial production received by an agent, who devotes
$A\gt 0$
 determines the level of industrial productivity. The fixed cost is shared by all agents when the industrial economy operates. Then, the output of industrial production received by an agent, who devotes 
 $l_{Y}$
 units of labor, is
$l_{Y}$
 units of labor, is
 \begin{equation} y=A(l_{Y}-\delta/N)\text{.} \end{equation}
\begin{equation} y=A(l_{Y}-\delta/N)\text{.} \end{equation}
Due to the fixed cost 
 $\delta$
, the industrial market would not operate unless the population size
$\delta$
, the industrial market would not operate unless the population size 
 $N$
 is sufficiently large.
$N$
 is sufficiently large.
2.4. From Neolithic Revolution to industrialization
 In this section, we explore the evolution of the economy and impose the following parameter assumption: 
 $A\gt \varphi \rho ^{1-\alpha }\gt \theta \rho ^{1-\gamma }$
. The population begins as hunter-gatherers and evolves into an agricultural society before an industrial economy emerges. We will impose parameter restrictions to ensure the realistic scenario in which industrialization takes place only after the complete transition from hunting-gathering to agriculture.
$A\gt \varphi \rho ^{1-\alpha }\gt \theta \rho ^{1-\gamma }$
. The population begins as hunter-gatherers and evolves into an agricultural society before an industrial economy emerges. We will impose parameter restrictions to ensure the realistic scenario in which industrialization takes place only after the complete transition from hunting-gathering to agriculture.
We begin by assuming that each agent maximizes consumption given by
 \begin{equation} c=x+y=h+f+y\text{.} \end{equation}
\begin{equation} c=x+y=h+f+y\text{.} \end{equation}
Here, we make a simplifying assumption that there is perfect substitutability between food production 
 $x$
 and industrial production
$x$
 and industrial production 
 $y$
 in the consumption of agents. This assumption helps to keep our analysis tractable and is not entirely unrealistic because industrial production includes modern methods of food production that requires fixed investment.
$y$
 in the consumption of agents. This assumption helps to keep our analysis tractable and is not entirely unrealistic because industrial production includes modern methods of food production that requires fixed investment.
 In the initial stage, there is no industrial production, so we have 
 $l_{Y}=0$
. An agent’s decision is to choose labor allocation between hunting-gathering
$l_{Y}=0$
. An agent’s decision is to choose labor allocation between hunting-gathering 
 $l_{H}$
 and farming
$l_{H}$
 and farming 
 $l_{F}$
 to maximize food production
$l_{F}$
 to maximize food production 
 $x$
 given by
$x$
 given by
 \begin{equation} x=h+f=\frac{l_{H}}{\overline{l}_{H}N}\theta (\overline{l}_{H}N)^{\gamma }(Z_{H})^{1-\gamma }+\varphi (l_{F})^{\alpha }z^{1-\alpha }=\left ( l-l_{F}\right ) \theta \left ( \frac{Z_{H}}{\overline{l}_{H}N}\right ) ^{1-\gamma }+\varphi \rho ^{1-\alpha }l_{F}\text{,} \end{equation}
\begin{equation} x=h+f=\frac{l_{H}}{\overline{l}_{H}N}\theta (\overline{l}_{H}N)^{\gamma }(Z_{H})^{1-\gamma }+\varphi (l_{F})^{\alpha }z^{1-\alpha }=\left ( l-l_{F}\right ) \theta \left ( \frac{Z_{H}}{\overline{l}_{H}N}\right ) ^{1-\gamma }+\varphi \rho ^{1-\alpha }l_{F}\text{,} \end{equation}
where we have used the resource constraint on labor 
 $l_{H}+l_{F}=l$
 and the fixed ratio of land to farming labor
$l_{H}+l_{F}=l$
 and the fixed ratio of land to farming labor 
 $z=\rho l_{F}$
. The first-order condition is given by
$z=\rho l_{F}$
. The first-order condition is given by
 \begin{equation} \frac{\partial x}{\partial l_{F}}=-\theta \left ( \frac{Z_{H}}{\overline{l}_{H}N}\right ) ^{1-\gamma }+\varphi \rho ^{1-\alpha }=-\theta \left [ \frac{Z-\rho l_{F}N}{(l-l_{F})N}\right ] ^{1-\gamma }+\varphi \rho ^{1-\alpha }\text{,} \end{equation}
\begin{equation} \frac{\partial x}{\partial l_{F}}=-\theta \left ( \frac{Z_{H}}{\overline{l}_{H}N}\right ) ^{1-\gamma }+\varphi \rho ^{1-\alpha }=-\theta \left [ \frac{Z-\rho l_{F}N}{(l-l_{F})N}\right ] ^{1-\gamma }+\varphi \rho ^{1-\alpha }\text{,} \end{equation}
where we have invoked symmetry 
 $\{l_{H},l_{F}\}=\{\overline{l}_{H},\overline{l}_{F}\}$
 and also used the resource constraint on land
$\{l_{H},l_{F}\}=\{\overline{l}_{H},\overline{l}_{F}\}$
 and also used the resource constraint on land 
 $Z_{H}=Z-\rho l_{F}N$
. In (11),
$Z_{H}=Z-\rho l_{F}N$
. In (11), 
 $\varphi \rho ^{1-\alpha }$
 is the marginal product of farming labor
$\varphi \rho ^{1-\alpha }$
 is the marginal product of farming labor 
 $l_{F}$
, whereas
$l_{F}$
, whereas 
 $\theta \left [ \frac{Z-\rho l_{F}N}{(l-l_{F})N}\right ] ^{1-\gamma }$
 is the average product of hunting labor
$\theta \left [ \frac{Z-\rho l_{F}N}{(l-l_{F})N}\right ] ^{1-\gamma }$
 is the average product of hunting labor 
 $l_{H}=l-l_{F}$
. In the following subsections, we first compare these two objects under different population levels.
$l_{H}=l-l_{F}$
. In the following subsections, we first compare these two objects under different population levels.
2.4.1. Stage 1: Hunting-gathering
Equation (11) implies that if the following inequality holds:
 \begin{equation} N\lt \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\text{,} \end{equation}
\begin{equation} N\lt \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\text{,} \end{equation}
then 
 $\partial x/\partial l_{F}\lt 0$
 even at
$\partial x/\partial l_{F}\lt 0$
 even at 
 $l_{F}=0$
. In this case, all labor is allocated to hunting-gathering
$l_{F}=0$
. In this case, all labor is allocated to hunting-gathering 
 $l_{H}=l$
 and the per capita output of food production is given by
$l_{H}=l$
 and the per capita output of food production is given by
 \begin{equation} x=h=\theta l^{\gamma }\left ( \frac{Z}{N}\right ) ^{1-\gamma }\text{,} \end{equation}
\begin{equation} x=h=\theta l^{\gamma }\left ( \frac{Z}{N}\right ) ^{1-\gamma }\text{,} \end{equation}
which is increasing in hunting productivity 
 $\theta$
, labor supply
$\theta$
, labor supply 
 $l$
, and the amount of land
$l$
, and the amount of land 
 $Z$
 but decreasing in the population size
$Z$
 but decreasing in the population size 
 $N$
 due to the decreasing returns to labor in hunting-gathering.
$N$
 due to the decreasing returns to labor in hunting-gathering.
2.4.2. Stage 2: From hunting-gathering to agriculture
 Equation (11) and 
 $\rho l_{F}N\lt Z$
 imply that if the following inequalities hold:
$\rho l_{F}N\lt Z$
 imply that if the following inequalities hold:
 \begin{equation} \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\lt N\lt \frac{Z}{\rho l}\text{,} \end{equation}
\begin{equation} \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\lt N\lt \frac{Z}{\rho l}\text{,} \end{equation}
then 
 $\partial x/\partial l_{F}=0$
 at some interior values of
$\partial x/\partial l_{F}=0$
 at some interior values of 
 $\{l_{F},l_{H}\}\in (0,l)$
. In this case, the transition from hunting-gathering to agriculture begins. The first inequality shows that a reduction in hunting productivity
$\{l_{F},l_{H}\}\in (0,l)$
. In this case, the transition from hunting-gathering to agriculture begins. The first inequality shows that a reduction in hunting productivity 
 $\theta$
 or an increase in population size
$\theta$
 or an increase in population size 
 $N$
 could trigger this transition. In our static model, the reduction in hunting productivity
$N$
 could trigger this transition. In our static model, the reduction in hunting productivity 
 $\theta$
 can capture the extinction of large herding animals analyzed in Smith (Reference Smith1975),Footnote 
11
 whereas an exogenous increase in population size
$\theta$
 can capture the extinction of large herding animals analyzed in Smith (Reference Smith1975),Footnote 
11
 whereas an exogenous increase in population size 
 $N$
 can capture the population pressure theory discussed in Cohen (Reference Cohen1977). However, as we will show, these results would be quite different in our dynamic model with endogenous population growth.
$N$
 can capture the population pressure theory discussed in Cohen (Reference Cohen1977). However, as we will show, these results would be quite different in our dynamic model with endogenous population growth.
During the gradual transition from hunting-gathering to agriculture, the per capita output of food production is given by
 \begin{equation} x=h+f=\left ( l-l_{F}\right ) \theta \left ( \frac{Z_{H}}{\overline{l}_{H}N}\right ) ^{1-\gamma }+\varphi \rho ^{1-\alpha }l_{F}=\varphi \rho ^{1-\alpha }l\text{,} \end{equation}
\begin{equation} x=h+f=\left ( l-l_{F}\right ) \theta \left ( \frac{Z_{H}}{\overline{l}_{H}N}\right ) ^{1-\gamma }+\varphi \rho ^{1-\alpha }l_{F}=\varphi \rho ^{1-\alpha }l\text{,} \end{equation}
which uses 
 $\theta \left [ Z_{H}/(\overline{l}_{H}N)\right ] ^{1-\gamma }=\varphi \rho ^{1-\alpha }$
 from (11). Equation (15) shows that
$\theta \left [ Z_{H}/(\overline{l}_{H}N)\right ] ^{1-\gamma }=\varphi \rho ^{1-\alpha }$
 from (11). Equation (15) shows that 
 $x$
 is increasing in labor supply
$x$
 is increasing in labor supply 
 $l$
 and agricultural productivity
$l$
 and agricultural productivity 
 $\varphi \rho ^{1-\alpha }$
.
$\varphi \rho ^{1-\alpha }$
.
2.4.3. Stage 3: Complete transition to agriculture
 When 
 $N\gt Z/(\rho l)$
, the transition from hunting-gathering to agriculture is complete (i.e.,
$N\gt Z/(\rho l)$
, the transition from hunting-gathering to agriculture is complete (i.e., 
 $l_{F}=l$
) because
$l_{F}=l$
) because 
 $Z_{H}=0$
. At this stage of the economy, an industrial market still does not emerge if the population size is insufficient to cover the fixed cost
$Z_{H}=0$
. At this stage of the economy, an industrial market still does not emerge if the population size is insufficient to cover the fixed cost 
 $\delta$
. This threshold value of
$\delta$
. This threshold value of 
 $N$
 is implicitly determined by the following equality:
$N$
 is implicitly determined by the following equality:
 \begin{equation} \varphi l^{\alpha }\left ( \frac{Z}{N}\right ) ^{1-\alpha }=A\left ( l-\frac{\delta }{N}\right ) \text{,} \end{equation}
\begin{equation} \varphi l^{\alpha }\left ( \frac{Z}{N}\right ) ^{1-\alpha }=A\left ( l-\frac{\delta }{N}\right ) \text{,} \end{equation}
in which the left-hand side is farming output per capita when 
 $l_{F}=l$
 and decreasing in
$l_{F}=l$
 and decreasing in 
 $N$
, whereas the right-hand side is the industrial output per capita when
$N$
, whereas the right-hand side is the industrial output per capita when 
 $l_{Y}=l$
 and increasing in
$l_{Y}=l$
 and increasing in 
 $N$
. A simple graphical analysis would confirm that there exists a unique cutoff value of
$N$
. A simple graphical analysis would confirm that there exists a unique cutoff value of 
 $N$
 for the emergence of an industrial economy, which is denoted as
$N$
 for the emergence of an industrial economy, which is denoted as 
 $N_{I}$
 and has the following comparative statics:
$N_{I}$
 and has the following comparative statics:
 \begin{equation} N_{I}(\underset{+}{\varphi },\underset{+}{Z},\underset{+}{\delta },\underset{-}{A},\underset{-}{l})\text{.} \end{equation}
\begin{equation} N_{I}(\underset{+}{\varphi },\underset{+}{Z},\underset{+}{\delta },\underset{-}{A},\underset{-}{l})\text{.} \end{equation}
This implies that by making agriculture more productive, higher agricultural productivity 
 $\varphi$
 delays industrialization, which contradicts the evidence discussed in Nurkse (Reference Nurkse1953).Footnote 
12
 As we will show, this counterfactual result will be overturned under endogenous population growth.
$\varphi$
 delays industrialization, which contradicts the evidence discussed in Nurkse (Reference Nurkse1953).Footnote 
12
 As we will show, this counterfactual result will be overturned under endogenous population growth.
In summary, if the following inequality holds:Footnote 13
 \begin{equation} \frac{Z}{\rho l}\lt N\lt N_{I}\text{,} \end{equation}
\begin{equation} \frac{Z}{\rho l}\lt N\lt N_{I}\text{,} \end{equation}
then the agents would be better off allocating all their labor to farming (i.e., 
 $l_{F}=l$
). In this case, the level of output per capita is given by
$l_{F}=l$
). In this case, the level of output per capita is given by
 \begin{equation} x=f=\varphi l^{\alpha }\left ( \frac{Z}{N}\right ) ^{1-\alpha }\text{,} \end{equation}
\begin{equation} x=f=\varphi l^{\alpha }\left ( \frac{Z}{N}\right ) ^{1-\alpha }\text{,} \end{equation}
which is increasing in agricultural productivity 
 $\varphi$
, labor supply
$\varphi$
, labor supply 
 $l$
, and the amount of land
$l$
, and the amount of land 
 $Z$
 but decreasing in the population size
$Z$
 but decreasing in the population size 
 $N$
 due to the decreasing returns to labor in farming when agricultural land is scarce.
$N$
 due to the decreasing returns to labor in farming when agricultural land is scarce.
2.4.4. Stage 4: Industrial economy
 If 
 $N\gt N_{I}$
, then the transition from agriculture to an industrial economy occurs. In this case, the level of output per capita is given by
$N\gt N_{I}$
, then the transition from agriculture to an industrial economy occurs. In this case, the level of output per capita is given by
 \begin{equation} y=A\left ( l-\frac{\delta }{N}\right ) \text{,} \end{equation}
\begin{equation} y=A\left ( l-\frac{\delta }{N}\right ) \text{,} \end{equation}
which is increasing in industrial productivity 
 $A$
, labor supply
$A$
, labor supply 
 $l$
, and population size
$l$
, and population size 
 $N$
 but decreasing in the fixed cost
$N$
 but decreasing in the fixed cost 
 $\delta$
 of industrial production. Equation (20) is obtained by setting
$\delta$
 of industrial production. Equation (20) is obtained by setting 
 $l_{Y}=l$
 in (8). When the population size is sufficiently large, the agents would immediately allocate all their labor to industrial production because the marginal product of industrial labor is greater than the marginal product of agricultural labor,Footnote 
14
 that is,
$l_{Y}=l$
 in (8). When the population size is sufficiently large, the agents would immediately allocate all their labor to industrial production because the marginal product of industrial labor is greater than the marginal product of agricultural labor,Footnote 
14
 that is,
 \begin{equation*} A\gt \varphi \rho ^{1-\alpha }\gt \varphi (l_{F})^{\alpha -1}\left ( \frac {Z}{N}\right ) ^{1-\alpha }\gt \alpha \varphi (l_{F})^{\alpha -1}\left ( \frac {Z}{N}\right ) ^{1-\alpha } \end{equation*}
\begin{equation*} A\gt \varphi \rho ^{1-\alpha }\gt \varphi (l_{F})^{\alpha -1}\left ( \frac {Z}{N}\right ) ^{1-\alpha }\gt \alpha \varphi (l_{F})^{\alpha -1}\left ( \frac {Z}{N}\right ) ^{1-\alpha } \end{equation*}
for 
 $l_{F}\gt Z/(\rho N)$
.Footnote 
15
 Naturally, we assume that it is infeasible for humans to return to hunting-gathering at this stage.Footnote 
16
$l_{F}\gt Z/(\rho N)$
.Footnote 
15
 Naturally, we assume that it is infeasible for humans to return to hunting-gathering at this stage.Footnote 
16
2.4.5. Summary
In this section, we summarize the level of consumption per capita at different levels of population as follows:
 \begin{equation} c=x+y=\left \{ \begin{array}{l@{\quad}l} h=\theta l^{\gamma }\left ( \frac{Z}{N}\right ) ^{1-\gamma } & \text{ for } N\lt \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l} \\[4pt] h+f=\varphi \rho ^{1-\alpha }l & \text{ for } \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\lt N\lt \frac{Z}{\rho l} \\[4pt] f=\varphi l^{\alpha }\left ( \frac{Z}{N}\right ) ^{1-\alpha } & \text{ for } \frac{Z}{\rho l}\lt N\lt N_{I} \\[4pt] y=A\left ( l-\frac{\delta }{N}\right ) & \text{ for } N\gt N_{I} \end{array}\right. \text{.} \end{equation}
\begin{equation} c=x+y=\left \{ \begin{array}{l@{\quad}l} h=\theta l^{\gamma }\left ( \frac{Z}{N}\right ) ^{1-\gamma } & \text{ for } N\lt \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l} \\[4pt] h+f=\varphi \rho ^{1-\alpha }l & \text{ for } \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\lt N\lt \frac{Z}{\rho l} \\[4pt] f=\varphi l^{\alpha }\left ( \frac{Z}{N}\right ) ^{1-\alpha } & \text{ for } \frac{Z}{\rho l}\lt N\lt N_{I} \\[4pt] y=A\left ( l-\frac{\delta }{N}\right ) & \text{ for } N\gt N_{I} \end{array}\right. \text{.} \end{equation}
Equation (21) presents the level of per capita consumption 
 $c$
 as population
$c$
 as population 
 $N$
 increases. In summary,
$N$
 increases. In summary, 
 $c$
 is initially falling due to the decreasing returns to labor in hunting-gathering. Then,
$c$
 is initially falling due to the decreasing returns to labor in hunting-gathering. Then, 
 $c$
 reaches to a stationary level (from above) when the gradual transition from hunting-gathering to agriculture begins. Therefore, before the transition to agriculture, hunter-gatherers enjoy a higher level of consumption than the later farmers, which is consistent with archaeological evidence, see, for example, Cohen and Armelagos (Reference Cohen, Armelagos, Cohen and Armelagos1984). However, our model implies that the hunter-gatherers would have experienced a subsequent fall in consumption if they did not adopt farming due to the decreasing returns to labor in hunting-gathering. When the transition from hunting-gathering to agriculture is complete,
$c$
 reaches to a stationary level (from above) when the gradual transition from hunting-gathering to agriculture begins. Therefore, before the transition to agriculture, hunter-gatherers enjoy a higher level of consumption than the later farmers, which is consistent with archaeological evidence, see, for example, Cohen and Armelagos (Reference Cohen, Armelagos, Cohen and Armelagos1984). However, our model implies that the hunter-gatherers would have experienced a subsequent fall in consumption if they did not adopt farming due to the decreasing returns to labor in hunting-gathering. When the transition from hunting-gathering to agriculture is complete, 
 $c$
 becomes falling again due to the decreasing returns to labor in farming when agricultural land is scarce. When the industrial economy emerges,
$c$
 becomes falling again due to the decreasing returns to labor in farming when agricultural land is scarce. When the industrial economy emerges, 
 $c$
 becomes rising due to the increasing returns to scale in the presence of a fixed cost of industrial production and converges toward a steady-state level given by
$c$
 becomes rising due to the increasing returns to scale in the presence of a fixed cost of industrial production and converges toward a steady-state level given by 
 $y^{\ast }=Al$
 as
$y^{\ast }=Al$
 as 
 $N\rightarrow \infty$
.
$N\rightarrow \infty$
.
3. A dynamic model with endogenous population growth
 The previous section presents a static model with an exogenous level of population. This section extends the model into a dynamic setting with endogenous population growth. We follow Locay (Reference Locay1989) and Baker (Reference Baker2008) to consider overlapping generations of agents. Each agent lives for two periods. Each adult agent at time 
 $t$
 has the following utility function:
$t$
 has the following utility function:
 \begin{equation} u_{t}=(1-\sigma )\ln c_{t}+\sigma \ln n_{t+1}\text{,} \end{equation}
\begin{equation} u_{t}=(1-\sigma )\ln c_{t}+\sigma \ln n_{t+1}\text{,} \end{equation}
where the parameter 
 $\sigma \in (0,1)$
 measures the preference for fertility and
$\sigma \in (0,1)$
 measures the preference for fertility and 
 $n_{t+1}$
 is the agent’s number of children, who then become adults at time
$n_{t+1}$
 is the agent’s number of children, who then become adults at time 
 $t+1$
. Raising children is costly, and the level of consumption net of the fertility cost is given by
$t+1$
. Raising children is costly, and the level of consumption net of the fertility cost is given by
 \begin{equation} c_{t}=x_{t}+y_{t}-\beta n_{t+1}\text{,} \end{equation}
\begin{equation} c_{t}=x_{t}+y_{t}-\beta n_{t+1}\text{,} \end{equation}
where the parameter 
 $\beta \gt 0$
 determines the cost of fertility. Substituting (23) into (22), we derive the utility-maximizing level of fertility
$\beta \gt 0$
 determines the cost of fertility. Substituting (23) into (22), we derive the utility-maximizing level of fertility 
 $n_{t+1}$
 as
$n_{t+1}$
 as
 \begin{equation} n_{t+1}=\frac{\sigma }{\beta }(x_{t}+y_{t}) \end{equation}
\begin{equation} n_{t+1}=\frac{\sigma }{\beta }(x_{t}+y_{t}) \end{equation}
and 
 $c_{t}=(1-\sigma )(x_{t}+y_{t})$
 in which the agent maximizes
$c_{t}=(1-\sigma )(x_{t}+y_{t})$
 in which the agent maximizes 
 $x_{t}+y_{t}$
 as in Section 2. Each adult agent has
$x_{t}+y_{t}$
 as in Section 2. Each adult agent has 
 $n_{t+1}$
 children, and the number of adult agents at time
$n_{t+1}$
 children, and the number of adult agents at time 
 $t$
 is
$t$
 is 
 $N_{t}$
. Therefore, the law of motion for the adult population size (i.e., the labor force) is given by
$N_{t}$
. Therefore, the law of motion for the adult population size (i.e., the labor force) is given by
 \begin{equation} N_{t+1}=n_{t+1}N_{t}=\frac{\sigma }{\beta }(x_{t}+y_{t})N_{t}\text{,} \end{equation}
\begin{equation} N_{t+1}=n_{t+1}N_{t}=\frac{\sigma }{\beta }(x_{t}+y_{t})N_{t}\text{,} \end{equation}
and the adult population growth rate at time 
 $t$
 is
$t$
 is
 \begin{equation} \frac{\Delta N_{t}}{N_{t}}\equiv \frac{N_{t+1}-N_{t}}{N_{t}}=\frac{\sigma }{\beta }(x_{t}+y_{t})-1\text{,} \end{equation}
\begin{equation} \frac{\Delta N_{t}}{N_{t}}\equiv \frac{N_{t+1}-N_{t}}{N_{t}}=\frac{\sigma }{\beta }(x_{t}+y_{t})-1\text{,} \end{equation}
which will be simply referred to as the population growth rate. In the following subsection, we will use the information from Section 2 to derive the population dynamics.
3.1. Stage 1: Hunting-gathering
Given an initial level of population:
 \begin{equation} N_{0}\lt \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\text{,} \end{equation}
\begin{equation} N_{0}\lt \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\text{,} \end{equation}
the human population engages in hunting-gathering only. Substituting (13) into (26) yields the growth rate of population as
 \begin{equation} \frac{\Delta N_{t}}{N_{t}}=\frac{\sigma }{\beta }\theta l^{\gamma }\left ( \frac{Z}{N_{t}}\right ) ^{1-\gamma }-1\text{,} \end{equation}
\begin{equation} \frac{\Delta N_{t}}{N_{t}}=\frac{\sigma }{\beta }\theta l^{\gamma }\left ( \frac{Z}{N_{t}}\right ) ^{1-\gamma }-1\text{,} \end{equation}
which yields the following steady-state level of population in the hunting-gathering era:
 \begin{equation} N_{H}^{\ast }=\left ( \frac{\sigma }{\beta }\theta l^{\gamma }\right ) ^{1/(1-\gamma )}Z\text{.} \end{equation}
\begin{equation} N_{H}^{\ast }=\left ( \frac{\sigma }{\beta }\theta l^{\gamma }\right ) ^{1/(1-\gamma )}Z\text{.} \end{equation}
If the following inequality holds:
 \begin{equation} N_{H}^{\ast }\lt \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\Leftrightarrow \frac{\sigma }{\beta }\varphi \rho ^{1-\alpha }l\lt 1\text{,} \end{equation}
\begin{equation} N_{H}^{\ast }\lt \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\Leftrightarrow \frac{\sigma }{\beta }\varphi \rho ^{1-\alpha }l\lt 1\text{,} \end{equation}
then the human population would remain as hunter-gatherers indefinitely. Substituting (29) into (13) yields 
 $x^{\ast }=\beta/\sigma$
, which is increasing in fertility cost
$x^{\ast }=\beta/\sigma$
, which is increasing in fertility cost 
 $\beta$
 and decreasing in the degree
$\beta$
 and decreasing in the degree 
 $\sigma$
 of fertility preference but independent of hunting productivity
$\sigma$
 of fertility preference but independent of hunting productivity 
 $\theta$
 and land
$\theta$
 and land 
 $Z$
. In other words, the population is in a hunting-gathering Malthusian trap, in which higher hunting productivity
$Z$
. In other words, the population is in a hunting-gathering Malthusian trap, in which higher hunting productivity 
 $\theta$
 and more land
$\theta$
 and more land 
 $Z$
 increase the level of population
$Z$
 increase the level of population 
 $N_{H}^{\ast }$
 but not the level of income
$N_{H}^{\ast }$
 but not the level of income 
 $x^{\ast }$
.
$x^{\ast }$
.
 Alternatively, if 
 $\sigma \varphi \rho ^{1-\alpha }l\gt \beta$
, then an agricultural society would emerge. Therefore, the transition from hunting-gathering to agriculture occurs under the following conditions: a low fertility cost
$\sigma \varphi \rho ^{1-\alpha }l\gt \beta$
, then an agricultural society would emerge. Therefore, the transition from hunting-gathering to agriculture occurs under the following conditions: a low fertility cost 
 $\beta$
, a strong fertility preference
$\beta$
, a strong fertility preference 
 $\sigma$
, a high level of agricultural productivity
$\sigma$
, a high level of agricultural productivity 
 $\varphi \rho ^{1-\alpha }$
, and a high level of labor supply
$\varphi \rho ^{1-\alpha }$
, and a high level of labor supply 
 $l$
. A strong fertility preference
$l$
. A strong fertility preference 
 $\sigma$
 and a low fertility cost
$\sigma$
 and a low fertility cost 
 $\beta$
 give rise to a higher level of population and make it more likely to cross the population threshold for the emergence of agriculture in a Boserupian manner, but they also reduce income
$\beta$
 give rise to a higher level of population and make it more likely to cross the population threshold for the emergence of agriculture in a Boserupian manner, but they also reduce income 
 $x^{\ast }=\beta/\sigma$
 in case the population remains in a hunting-gathering Malthusian trap. Although a higher level of hunting productivity
$x^{\ast }=\beta/\sigma$
 in case the population remains in a hunting-gathering Malthusian trap. Although a higher level of hunting productivity 
 $\theta$
 and a larger amount of land
$\theta$
 and a larger amount of land 
 $Z$
 also increase population, they increase the endogenous threshold for agriculture as well by making hunting-gathering more attractive. These opposite effects cancel each other, and hence, hunting productivity
$Z$
 also increase population, they increase the endogenous threshold for agriculture as well by making hunting-gathering more attractive. These opposite effects cancel each other, and hence, hunting productivity 
 $\theta$
 and the amount of land
$\theta$
 and the amount of land 
 $Z$
 do not affect the transition to agriculture, which stands in stark contrast to the case of exogenous population.
$Z$
 do not affect the transition to agriculture, which stands in stark contrast to the case of exogenous population.
 Finally, high agricultural productivity 
 $\varphi \rho ^{1-\alpha }$
 reduces the endogenous threshold by making agriculture more attractive, and hence, a higher level of agricultural productivity
$\varphi \rho ^{1-\alpha }$
 reduces the endogenous threshold by making agriculture more attractive, and hence, a higher level of agricultural productivity 
 $\varphi \rho ^{1-\alpha }$
 can trigger the Neolithic Revolution. This finding is consistent with the empirical evidence in Olsson and Hibbs (Reference Olsson and Hibbs2005), who find that favorable biogeographic conditions can trigger the transition to agriculture. Olsson (Reference Olsson2001) examines the archeological evidence in the Jordan Valley and concludes that the abundance of species suitable for agriculture was one of the key reasons for the transition to agriculture. This abundance of agricultural species corresponds to a high level of agricultural productivity in our model. Furthermore, our analysis implies that climate change that affects agricultural productivity would also affect the transition to agriculture, for example, Richerson et al. (Reference Richerson, Boyd and Bettinger2001) argue that the climatic conditions during the most recent Ice Age were hostile to agriculture and made the transition to agriculture impossible at that time.
$\varphi \rho ^{1-\alpha }$
 can trigger the Neolithic Revolution. This finding is consistent with the empirical evidence in Olsson and Hibbs (Reference Olsson and Hibbs2005), who find that favorable biogeographic conditions can trigger the transition to agriculture. Olsson (Reference Olsson2001) examines the archeological evidence in the Jordan Valley and concludes that the abundance of species suitable for agriculture was one of the key reasons for the transition to agriculture. This abundance of agricultural species corresponds to a high level of agricultural productivity in our model. Furthermore, our analysis implies that climate change that affects agricultural productivity would also affect the transition to agriculture, for example, Richerson et al. (Reference Richerson, Boyd and Bettinger2001) argue that the climatic conditions during the most recent Ice Age were hostile to agriculture and made the transition to agriculture impossible at that time.
3.2. Stage 2: From hunting-gathering to agriculture
 Suppose the population size 
 $N_{t}$
 crosses the first threshold, that is,
$N_{t}$
 crosses the first threshold, that is,
 \begin{equation} \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\lt N_{t}\lt \frac{Z}{\rho l}\text{.} \end{equation}
\begin{equation} \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\lt N_{t}\lt \frac{Z}{\rho l}\text{.} \end{equation}
Then, the transition from hunting-gathering to agriculture begins. We can substitute (15) into (26) to derive the population growth rate as
 \begin{equation} \frac{\Delta N_{t}}{N_{t}}=\frac{\sigma }{\beta }\varphi \rho ^{1-\alpha }l-1\gt 0\text{,} \end{equation}
\begin{equation} \frac{\Delta N_{t}}{N_{t}}=\frac{\sigma }{\beta }\varphi \rho ^{1-\alpha }l-1\gt 0\text{,} \end{equation}
which is positive if and only if the transition to agriculture occurs (i.e., 
 $\sigma \varphi \rho ^{1-\alpha }l\gt \beta$
) and implies that population
$\sigma \varphi \rho ^{1-\alpha }l\gt \beta$
) and implies that population 
 $N_{t}$
 increases over time during the gradual transition from hunting-gathering to agriculture.
$N_{t}$
 increases over time during the gradual transition from hunting-gathering to agriculture.
3.3. Stage 3: Complete transition to agriculture
 Given (32), the level of population 
 $N_{t}$
 eventually crosses the second threshold, that is,
$N_{t}$
 eventually crosses the second threshold, that is,
 \begin{equation} \frac{Z}{\rho l}\lt N_{t}\lt N_{I}\text{,} \end{equation}
\begin{equation} \frac{Z}{\rho l}\lt N_{t}\lt N_{I}\text{,} \end{equation}
where 
 $N_{I}$
 is implicitly given in (16) and (17). At this stage, we can substitute (19) into (26) to derive the growth rate of population as
$N_{I}$
 is implicitly given in (16) and (17). At this stage, we can substitute (19) into (26) to derive the growth rate of population as
 \begin{equation} \frac{\Delta N_{t}}{N_{t}}=\frac{\sigma }{\beta }\varphi l^{\alpha }\left ( \frac{Z}{N_{t}}\right ) ^{1-\alpha }-1\text{,} \end{equation}
\begin{equation} \frac{\Delta N_{t}}{N_{t}}=\frac{\sigma }{\beta }\varphi l^{\alpha }\left ( \frac{Z}{N_{t}}\right ) ^{1-\alpha }-1\text{,} \end{equation}
which yields a steady-state level of population in agriculture as
 \begin{equation} N_{A}^{\ast }=\left ( \frac{\sigma }{\beta }\varphi l^{\alpha }\right ) ^{1/(1-\alpha )}Z\text{.} \end{equation}
\begin{equation} N_{A}^{\ast }=\left ( \frac{\sigma }{\beta }\varphi l^{\alpha }\right ) ^{1/(1-\alpha )}Z\text{.} \end{equation}
If 
 $N_{t}$
 reaches
$N_{t}$
 reaches 
 $N_{A}^{\ast }$
 before reaching
$N_{A}^{\ast }$
 before reaching 
 $N_{I}$
, then the economy would remain as an agricultural society indefinitely. Substituting (35) into (19) yields
$N_{I}$
, then the economy would remain as an agricultural society indefinitely. Substituting (35) into (19) yields 
 $x^{\ast }=\beta/\sigma$
, which is once again increasing in fertility cost
$x^{\ast }=\beta/\sigma$
, which is once again increasing in fertility cost 
 $\beta$
 and decreasing in the degree
$\beta$
 and decreasing in the degree 
 $\sigma$
 of fertility preference but independent of agricultural productivity
$\sigma$
 of fertility preference but independent of agricultural productivity 
 $\varphi$
 and land
$\varphi$
 and land 
 $Z$
. In other words, the population is now in an agricultural Malthusian trap, in which higher agricultural productivity
$Z$
. In other words, the population is now in an agricultural Malthusian trap, in which higher agricultural productivity 
 $\varphi$
 and more land
$\varphi$
 and more land 
 $Z$
 increase the level of population
$Z$
 increase the level of population 
 $N_{A}^{\ast }$
 but not the level of income
$N_{A}^{\ast }$
 but not the level of income 
 $x^{\ast }$
.
$x^{\ast }$
.
3.4. Stage 4: Industrial economy
 If the level of population 
 $N_{t}$
 manages to cross the third threshold
$N_{t}$
 manages to cross the third threshold 
 $N_{I}$
, then an industrial economy emerges. In this case, we can substitute (20) into (26) to derive the population growth rate as
$N_{I}$
, then an industrial economy emerges. In this case, we can substitute (20) into (26) to derive the population growth rate as
 \begin{equation} \frac{\Delta N_{t}}{N_{t}}=\frac{\sigma A}{\beta }\left ( l-\frac{\delta }{N_{t}}\right ) -1\text{,} \end{equation}
\begin{equation} \frac{\Delta N_{t}}{N_{t}}=\frac{\sigma A}{\beta }\left ( l-\frac{\delta }{N_{t}}\right ) -1\text{,} \end{equation}
which is increasing in 
 $N_{t}$
. Setting
$N_{t}$
. Setting 
 $\Delta N_{t}/N_{t}=0$
 yields the following level:
$\Delta N_{t}/N_{t}=0$
 yields the following level:
 \begin{equation} N_{I}^{\ast }=\frac{\delta }{l-\beta/(\sigma A)}\text{,} \end{equation}
\begin{equation} N_{I}^{\ast }=\frac{\delta }{l-\beta/(\sigma A)}\text{,} \end{equation}
above which the population grows over time during the industrial era.
 Figure 1 shows that if and only if 
 $N_{A}^{\ast }\gt N_{I}^{\ast }$
, then
$N_{A}^{\ast }\gt N_{I}^{\ast }$
, then 
 $N_{t}$
 would reach the third threshold
$N_{t}$
 would reach the third threshold 
 $N_{I}$
 and trigger the emergence of an industrial economy. When
$N_{I}$
 and trigger the emergence of an industrial economy. When 
 $N_{t}\gt N_{I}$
, the output level
$N_{t}\gt N_{I}$
, the output level 
 $x_{t}+y_{t}$
 is higher under industrial production than under agricultural production. From (35) and (37), the inequality
$x_{t}+y_{t}$
 is higher under industrial production than under agricultural production. From (35) and (37), the inequality 
 $N_{A}^{\ast }\gt N_{I}^{\ast }$
 is equivalent to
$N_{A}^{\ast }\gt N_{I}^{\ast }$
 is equivalent to
 \begin{equation} \left ( l-\frac{\beta }{\sigma A}\right ) \left ( \frac{\sigma }{\beta }\varphi l^{\alpha }\right ) ^{1/(1-\alpha )}\frac{Z}{\delta }\gt 1\text{.} \end{equation}
\begin{equation} \left ( l-\frac{\beta }{\sigma A}\right ) \left ( \frac{\sigma }{\beta }\varphi l^{\alpha }\right ) ^{1/(1-\alpha )}\frac{Z}{\delta }\gt 1\text{.} \end{equation}
Therefore, the transition from an agricultural economy to an industrial economy occurs under the following conditions: a low fertility cost 
 $\beta$
, a strong fertility preference
$\beta$
, a strong fertility preference 
 $\sigma$
, a high level of agricultural productivity
$\sigma$
, a high level of agricultural productivity 
 $\varphi$
, a high level of labor supply
$\varphi$
, a high level of labor supply 
 $l$
, a large amount of land
$l$
, a large amount of land 
 $Z$
, a high level of industrial productivity
$Z$
, a high level of industrial productivity 
 $A$
, and a low fixed cost
$A$
, and a low fixed cost 
 $\delta$
 for operating industrial firms.
$\delta$
 for operating industrial firms.

Figure 1. Industrial threshold. Figure plots the population growth rate in (34) and (36) and shows that the economy switches from agriculture to industrial production when population crosses the threshold 
 ${\rm N}_{I}$
.Footnote 
17
${\rm N}_{I}$
.Footnote 
17
 As before, a strong fertility preference 
 $\sigma$
 and a low fertility cost
$\sigma$
 and a low fertility cost 
 $\beta$
 give rise to a higher level of population and make it more likely to cross the population threshold
$\beta$
 give rise to a higher level of population and make it more likely to cross the population threshold 
 $N_{I}$
 for the emergence of an industrial economy, but they also reduce income
$N_{I}$
 for the emergence of an industrial economy, but they also reduce income 
 $x^{\ast }=\beta/\sigma$
 in case the population remains in an agricultural Malthusian trap. Interestingly, unlike the case of exogenous population, a high level of agricultural productivity
$x^{\ast }=\beta/\sigma$
 in case the population remains in an agricultural Malthusian trap. Interestingly, unlike the case of exogenous population, a high level of agricultural productivity 
 $\varphi$
 can now trigger industrialization by raising the level of population. This result is consistent with the early work of Nurkse (Reference Nurkse1953) and Murphy et al. (Reference Murphy, Shleifer and Vishny1989) and also supported by the empirical evidence in Olsson and Hibbs (Reference Olsson and Hibbs2005) and Ang (Reference Ang2015), who find that favorable initial biogeographic conditions can contribute to the subsequent development in the industrial era and technology adoption in as late as 1500 AD, in addition to the Neolithic Revolution.
$\varphi$
 can now trigger industrialization by raising the level of population. This result is consistent with the early work of Nurkse (Reference Nurkse1953) and Murphy et al. (Reference Murphy, Shleifer and Vishny1989) and also supported by the empirical evidence in Olsson and Hibbs (Reference Olsson and Hibbs2005) and Ang (Reference Ang2015), who find that favorable initial biogeographic conditions can contribute to the subsequent development in the industrial era and technology adoption in as late as 1500 AD, in addition to the Neolithic Revolution.
 Furthermore, a high level of industrial productivity 
 $A$
 and a low fixed cost
$A$
 and a low fixed cost 
 $\delta$
 of industrial production reduce the endogenous threshold by making industrial production more attractive and can also trigger industrialization. Finally, if the population size reaches the industrial threshold, then a modern economy emerges and the population growth rate rises toward a steady-state value given by
$\delta$
 of industrial production reduce the endogenous threshold by making industrial production more attractive and can also trigger industrialization. Finally, if the population size reaches the industrial threshold, then a modern economy emerges and the population growth rate rises toward a steady-state value given by 
 $\Delta N/N=\frac{\sigma }{\beta }Al-1$
 in the long run.Footnote 
18
$\Delta N/N=\frac{\sigma }{\beta }Al-1$
 in the long run.Footnote 
18
3.5. Summary
We summarize all the above results in the following proposition:
Proposition 1. Under exogenous population growth, human society evolves from hunting-gathering to agriculture and then an industrial economy. Under endogenous population growth, the population may stop growing in a hunting-gathering society; in this case, the population remains as hunter-gatherers. The Neolithic Revolution occurs under a low fertility cost, strong fertility preference, high agricultural productivity, and high labor supply. The population may also stop growing in an agricultural society; in this case, the economy remains in an agricultural Malthusian trap. Industrialization occurs under a low fertility cost, strong fertility preference, high agricultural productivity, high labor supply, a large amount of agricultural land, high industrial productivity, and a low fixed cost of industrial production.
 
Proof. The population growth rate is summarized in (39). From (30), if 
 $\sigma \varphi \rho ^{1-\alpha }l\gt \beta$
, then
$\sigma \varphi \rho ^{1-\alpha }l\gt \beta$
, then 
 $N_{t}$
 reaches the agricultural threshold before the hunting-gathering steady-state
$N_{t}$
 reaches the agricultural threshold before the hunting-gathering steady-state 
 $N_{H}^{\ast }$
. If (38) holds, then
$N_{H}^{\ast }$
. If (38) holds, then 
 $N_{t}$
 reaches the industrial threshold
$N_{t}$
 reaches the industrial threshold 
 $N_{I}$
 before the agricultural steady state
$N_{I}$
 before the agricultural steady state 
 $N_{A}^{\ast }$
.
$N_{A}^{\ast }$
.
If the population manages to evolve from hunting-gathering to agriculture and then activate the emergence of an industrial economy, the dynamics of the population growth rate can be summarized as follows:
 \begin{equation} \frac{\Delta N_{t}}{N_{t}}=\frac{\sigma }{\beta }(x_{t}+y_{t})-1=\left \{ \begin{array}{l@{\quad}l} \frac{\sigma }{\beta }\theta l^{\gamma }\left ( \frac{Z}{N_{t}}\right ) ^{1-\gamma }-1 &\text{ for } N_{t}\lt \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l} \\[5pt] \frac{\sigma }{\beta }\varphi \rho ^{1-\alpha }l-1 & \text{ for } \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\lt N_{t}\lt \frac{Z}{\rho l} \\[5pt] \frac{\sigma }{\beta }\varphi l^{\alpha }\left ( \frac{Z}{N_{t}}\right ) ^{1-\alpha }-1 & \text{ for } \frac{Z}{\rho l}\lt N_{t}\lt N_{I} \\[5pt] \frac{\sigma }{\beta }A\left ( l-\frac{\delta }{N_{t}}\right ) -1 & \text{ for }s N_{t}\gt N_{I}\end{array}\right. \text{.} \end{equation}
\begin{equation} \frac{\Delta N_{t}}{N_{t}}=\frac{\sigma }{\beta }(x_{t}+y_{t})-1=\left \{ \begin{array}{l@{\quad}l} \frac{\sigma }{\beta }\theta l^{\gamma }\left ( \frac{Z}{N_{t}}\right ) ^{1-\gamma }-1 &\text{ for } N_{t}\lt \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l} \\[5pt] \frac{\sigma }{\beta }\varphi \rho ^{1-\alpha }l-1 & \text{ for } \left ( \frac{\theta }{\varphi \rho ^{1-\alpha }}\right ) ^{1/(1-\gamma )}\frac{Z}{l}\lt N_{t}\lt \frac{Z}{\rho l} \\[5pt] \frac{\sigma }{\beta }\varphi l^{\alpha }\left ( \frac{Z}{N_{t}}\right ) ^{1-\alpha }-1 & \text{ for } \frac{Z}{\rho l}\lt N_{t}\lt N_{I} \\[5pt] \frac{\sigma }{\beta }A\left ( l-\frac{\delta }{N_{t}}\right ) -1 & \text{ for }s N_{t}\gt N_{I}\end{array}\right. \text{.} \end{equation}
Figure 2 plots the population growth rate 
 $\Delta N_{t}/N_{t}$
 for the following three scenarios: (a) the population converges to a hunting-gathering Malthusian trap as discussed in Section 3.1; (b) the population converges to an agricultural Malthusian trap as discussed in Section 3.3; and (c) the population achieves long-run growth as discussed in Section 3.4.
$\Delta N_{t}/N_{t}$
 for the following three scenarios: (a) the population converges to a hunting-gathering Malthusian trap as discussed in Section 3.1; (b) the population converges to an agricultural Malthusian trap as discussed in Section 3.3; and (c) the population achieves long-run growth as discussed in Section 3.4.

Figure 2. Dynamics of population growth. Figure plots the following three scenarios: case (a) plots the population growth rate given by (28); case (b) plots the population growth rate given by (28), (32), and then (34); and case (c) plots the population growth rate given by (28), (32), (34), and then (36).
4. Empirical evidence
In this section, we use cross-country data to evaluate the effects of agricultural productivity on the transition to agriculture and the subsequent industrialization. Specifically, we follow Ashraf and Galor (Reference Ashraf and Galor2011) and Ang (Reference Ang2015) to employ a two-stage least squares regression:
 \begin{align}{\tau }_{j} ={\kappa _{1}\varphi }_{j}+\Phi _{j}+\epsilon _{1,j} \end{align}
\begin{align}{\tau }_{j} ={\kappa _{1}\varphi }_{j}+\Phi _{j}+\epsilon _{1,j} \end{align} 
 \begin{align} y_{j} ={\kappa _{2}}\widehat{{\tau }}_{j}+\Phi _{j}+\epsilon _{2,j} \end{align}
\begin{align} y_{j} ={\kappa _{2}}\widehat{{\tau }}_{j}+\Phi _{j}+\epsilon _{2,j} \end{align} 
where (40a) describes the first-stage regression and (40b) describes the second-stage regression.Footnote 
19
 In (40a), 
 ${\varphi }_{j}$
 denotes the prehistoric level of agricultural productivity in present-day country
${\varphi }_{j}$
 denotes the prehistoric level of agricultural productivity in present-day country 
 $j$
, for which we use an index of prehistoric biogeographic endowments as a proxy. We use data from Olsson and Hibbs (Reference Olsson and Hibbs2005) on domesticable wild animals and plants known to exist in prehistory (in location that corresponds to present-day country
$j$
, for which we use an index of prehistoric biogeographic endowments as a proxy. We use data from Olsson and Hibbs (Reference Olsson and Hibbs2005) on domesticable wild animals and plants known to exist in prehistory (in location that corresponds to present-day country 
 $j$
).Footnote 
20
 Then, we follow Ang (Reference Ang2015) to combine the two dimensions (animals and plants) into a single index by computing their first principal component.
$j$
).Footnote 
20
 Then, we follow Ang (Reference Ang2015) to combine the two dimensions (animals and plants) into a single index by computing their first principal component.
 In (40a), 
 ${\tau }_{j}$
 denotes the timing of agricultural transition, measured by the number of years before 2000 CE, whereas
${\tau }_{j}$
 denotes the timing of agricultural transition, measured by the number of years before 2000 CE, whereas 
 $\widehat{{\tau }}_{j}$
 in (40b) denotes the predicted value of
$\widehat{{\tau }}_{j}$
 in (40b) denotes the predicted value of 
 ${\tau }_{j}$
 from the first-stage regression.
${\tau }_{j}$
 from the first-stage regression. 
 $y_{j}$
 denotes the degree of industrialization, measured by the share of non-agricultural employment in 1991. We use data in 1991 because data on employment shares are only available from 1991 for many countries, so using earlier data would lead to a sharp reduction in the sample size. Our underlying assumption is that countries experiencing earlier industrialization should have higher shares of non-agricultural employment in the modern era.
$y_{j}$
 denotes the degree of industrialization, measured by the share of non-agricultural employment in 1991. We use data in 1991 because data on employment shares are only available from 1991 for many countries, so using earlier data would lead to a sharp reduction in the sample size. Our underlying assumption is that countries experiencing earlier industrialization should have higher shares of non-agricultural employment in the modern era.
 Our theory predicts that 
 $\kappa _{1}$
 and
$\kappa _{1}$
 and 
 $\kappa _{2}$
 are significantly positive.
$\kappa _{2}$
 are significantly positive. 
 $\kappa _{1}$
 being significantly positive implies that a higher level of agricultural productivity triggers an earlier transition to agriculture.
$\kappa _{1}$
 being significantly positive implies that a higher level of agricultural productivity triggers an earlier transition to agriculture. 
 $\kappa _{2}$
 being significantly positive implies that an earlier transition to agriculture (triggered by a higher level of agricultural productivity) also causes a higher degree of industrialization in the modern era.
$\kappa _{2}$
 being significantly positive implies that an earlier transition to agriculture (triggered by a higher level of agricultural productivity) also causes a higher degree of industrialization in the modern era.
 
 $\Phi$
 is a set of control variables, including the constant term. We follow Nunn and Puga (Reference Nunn and Puga2012) to control for a terrain ruggedness index, total land area, the ratio of fertile soil to land area, and the percentage of land area that is within 100 km of the nearest ice-free coast. We also consider continent fixed effects. Finally,
$\Phi$
 is a set of control variables, including the constant term. We follow Nunn and Puga (Reference Nunn and Puga2012) to control for a terrain ruggedness index, total land area, the ratio of fertile soil to land area, and the percentage of land area that is within 100 km of the nearest ice-free coast. We also consider continent fixed effects. Finally, 
 $\epsilon _{1,j}$
 and
$\epsilon _{1,j}$
 and 
 $\epsilon _{2,j}$
 are the error terms. Table 1 presents the summary statistics and data sources.
$\epsilon _{2,j}$
 are the error terms. Table 1 presents the summary statistics and data sources.
Table 1. Summary statistics

Data sources: World Bank Data for the share of non-agricultural employment; Putterman and Trainor (Reference Putterman and Trainor2006) for the timing of agricultural transition; Nunn and Puga (Reference Nunn and Puga2012) for other variables.
 Table 2 presents our regression results. Columns (1)–(2) do not include control variables except for the constant term and continent fixed effects, whereas columns (3)–(4) include control variables. The odd columns show the first-stage regression, whereas the even columns show the second-stage regression. The estimation results show that 
 $\kappa _{1}$
 and
$\kappa _{1}$
 and 
 $\kappa _{2}$
 are positive and statistically significant at least at the 5% level. Using the estimates in columns (3)–(4), we find that increasing the index of prehistoric biogeographic endowments by one (recall that this index ranges from –1.3 to 1.9) causes an earlier agricultural transition by about one millennium (1040 years) and a higher degree of industrialization reflected by a larger share of non-agricultural employment of 8.3% (
$\kappa _{2}$
 are positive and statistically significant at least at the 5% level. Using the estimates in columns (3)–(4), we find that increasing the index of prehistoric biogeographic endowments by one (recall that this index ranges from –1.3 to 1.9) causes an earlier agricultural transition by about one millennium (1040 years) and a higher degree of industrialization reflected by a larger share of non-agricultural employment of 8.3% (
 $=$
 1040
$=$
 1040 
 $\times$
 0.008%) in 1991.Footnote 
21
$\times$
 0.008%) in 1991.Footnote 
21
Table 2. Effects of agricultural productivity on agricultural transition and industrialization

 Notes: Robust standard errors in parentheses. 
 ${}^{\ast}$
${}^{\ast}$
 
 $p\lt 0.1,^{\ast\ast}$
$p\lt 0.1,^{\ast\ast}$
 
 $p\lt 0.05,^{\ast\ast\ast}$
$p\lt 0.05,^{\ast\ast\ast}$
 
 $p\lt 0.01$
. In Table 2, productivity refers to the index of prehistoric biogeographic endowments; transition refers to the timing of agricultural transition; and industrialization refers to the non-agricultural employment share in 1991. The IV regressions employ the index of prehistoric biogeographic endowments as an instrument for the timing of agricultural transition.
$p\lt 0.01$
. In Table 2, productivity refers to the index of prehistoric biogeographic endowments; transition refers to the timing of agricultural transition; and industrialization refers to the non-agricultural employment share in 1991. The IV regressions employ the index of prehistoric biogeographic endowments as an instrument for the timing of agricultural transition.
 To ensure the robustness of our results, we explore other proxies for industrialization in the modern era. First, we consider the non-agricultural share of GDP as a dependent variable in the second-stage regression.Footnote 
22
 In this case, we can use earlier data in 1980 and still retain a sample size of 62.Footnote 
23
 Second, we consider the log of GDP per capita as a dependent variable in the second-stage regression.Footnote 
24
 In this case, we can use even earlier data in 1950 and retain a sample size of 88.Footnote 
25
 Table 3 presents the regression results. The estimated coefficients of 
 $\kappa _{1}$
 and
$\kappa _{1}$
 and 
 $\kappa _{2}$
 remain positive and statistically significant at the 1% level. These estimates imply that increasing the index of prehistoric biogeographic endowments by one causes an earlier agricultural transition by about one millennium as before and gives rise to a larger share of non-agricultural GDP of 7.2% (
$\kappa _{2}$
 remain positive and statistically significant at the 1% level. These estimates imply that increasing the index of prehistoric biogeographic endowments by one causes an earlier agricultural transition by about one millennium as before and gives rise to a larger share of non-agricultural GDP of 7.2% (
 $=$
1031
$=$
1031 
 $\times$
 0.007%) in 1980 and also an increase in GDP per capita by 21.5% (
$\times$
 0.007%) in 1980 and also an increase in GDP per capita by 21.5% (
 $=$
 1077
$=$
 1077 
 $\times$
 0.0002
$\times$
 0.0002 
 $\times$
 100%) in 1950.Footnote 
26
$\times$
 100%) in 1950.Footnote 
26
Table 3. Robustness tests

 Notes: Robust standard errors in parentheses. 
 ${}^{*}\,p\lt 0.1, {}^{**}\,p\lt 0.05, {}^{***}\,p\lt 0.01$
. In Table 3, productivity refers to the index of prehistoric biogeographic endowments; transition refers to the timing of agricultural transition; GDP share refers to the non-agricultural share of GDP in 1980; and GDP per capita refers to the log of GDP per capita in 1950. The IV regressions employ the index of prehistoric biogeographic endowments as an instrument for the timing of agricultural transition.
${}^{*}\,p\lt 0.1, {}^{**}\,p\lt 0.05, {}^{***}\,p\lt 0.01$
. In Table 3, productivity refers to the index of prehistoric biogeographic endowments; transition refers to the timing of agricultural transition; GDP share refers to the non-agricultural share of GDP in 1980; and GDP per capita refers to the log of GDP per capita in 1950. The IV regressions employ the index of prehistoric biogeographic endowments as an instrument for the timing of agricultural transition.
5. Conclusion
In this study, we have developed a simple Malthusian model that captures the economic evolution of human society across the three stages of hunting-gathering, agriculture, and industrial production. We find that under endogenous population growth, the evolution to the next stage is not inevitable. If the population fails to reach the agricultural threshold, then the human population remains as hunter-gatherers. If the population fails to reach the industrial threshold, then the human population remains as agriculturalists. Our model identifies several potential causes for the Neolithic Revolution: a high level of agricultural productivity, a low cost of fertility, a strong preference for fertility, and a high level of labor supply. An implication is that the transitions to agriculture in different parts of the world (such as Central Mexico, China, the Middle East, and Sub-Saharan Africa) at different time periods could have been triggered by different reasons. Furthermore, the above conditions that trigger the Neolithic Revolution can also trigger the subsequent industrialization, but not necessarily vice versa because other conditions (such as a high level of industrial productivity and a low fixed cost of industrial production) may also trigger industrialization. Although our simple model is unlikely to capture all possible causes for the Neolithic Revolution and the subsequent industrialization, we find empirical support for agricultural productivity as a determinant for the timing of transitions to agriculture and the degree of industrialization in the modern era.
Acknowledgements
We would like to thank Francesco Zanetti for helpful comments. Chu gratefully acknowledges financial support from the Asia-Pacific Academy of Economics and Management at the University of Macau. The usual disclaimer applies.
Appendix A: Monopolistic market in the industrial era
In this appendix, we replace the reduced-form industrial production function in (7) by a modern monopolistic market with a standard CES aggregator:Footnote 27
 \begin{equation} Y=\left \{ \int _{0}^{1}[Y(i)]^{\varepsilon }di\right \} ^{1/\varepsilon }\text{,} \end{equation}
\begin{equation} Y=\left \{ \int _{0}^{1}[Y(i)]^{\varepsilon }di\right \} ^{1/\varepsilon }\text{,} \end{equation}
where 
 $\varepsilon \in (0,1)$
 determines the elasticity of substitution
$\varepsilon \in (0,1)$
 determines the elasticity of substitution 
 $1/(1-\varepsilon )$
 between differentiated products
$1/(1-\varepsilon )$
 between differentiated products 
 $Y(i)$
 for
$Y(i)$
 for 
 $i\in \lbrack 0,1]$
. Profit maximization yields the conditional demand function:
$i\in \lbrack 0,1]$
. Profit maximization yields the conditional demand function:
 \begin{equation} Y(i)=\left [ \frac{p}{p(i)}\right ] ^{1/(1-\varepsilon )}Y\Leftrightarrow p(i)=p\left [ \frac{Y}{Y(i)}\right ] ^{1-\varepsilon }\text{,} \end{equation}
\begin{equation} Y(i)=\left [ \frac{p}{p(i)}\right ] ^{1/(1-\varepsilon )}Y\Leftrightarrow p(i)=p\left [ \frac{Y}{Y(i)}\right ] ^{1-\varepsilon }\text{,} \end{equation}
where 
 $p$
 and
$p$
 and 
 $p(i)$
 are, respectively, the prices of
$p(i)$
 are, respectively, the prices of 
 $Y$
 and
$Y$
 and 
 $Y(i)$
 for
$Y(i)$
 for 
 $i\in \lbrack 0,1]$
.
$i\in \lbrack 0,1]$
.
 As in Krugman (Reference Krugman1979), operating an industrial firm requires a fixed cost 
 $\delta \gt 0$
 under which the output of
$\delta \gt 0$
 under which the output of 
 $Y(i)$
 is
$Y(i)$
 is
 \begin{equation} Y(i)=A[l_{Y}(i)-\delta ]\text{,} \end{equation}
\begin{equation} Y(i)=A[l_{Y}(i)-\delta ]\text{,} \end{equation}
where 
 $l_{Y}(i)$
 is the labor devoted to the production of
$l_{Y}(i)$
 is the labor devoted to the production of 
 $Y(i)$
. The profit function for firm
$Y(i)$
. The profit function for firm 
 $i$
 is
$i$
 is
 \begin{equation} \pi (i)=p(i)Y(i)-wl_{Y}(i)=pY^{1-\varepsilon }[Y(i)]^{\varepsilon }-w\left [ \frac{Y(i)}{A}+\delta \right ] \text{,} \end{equation}
\begin{equation} \pi (i)=p(i)Y(i)-wl_{Y}(i)=pY^{1-\varepsilon }[Y(i)]^{\varepsilon }-w\left [ \frac{Y(i)}{A}+\delta \right ] \text{,} \end{equation}
where 
 $w$
 is the wage rate of industrial labor. Profit maximization yields markup pricing:
$w$
 is the wage rate of industrial labor. Profit maximization yields markup pricing:
 \begin{equation} p(i)=\frac{1}{\varepsilon }\frac{w}{A}\gt \frac{w}{A}\text{,} \end{equation}
\begin{equation} p(i)=\frac{1}{\varepsilon }\frac{w}{A}\gt \frac{w}{A}\text{,} \end{equation}
where 
 $w/A$
 is the marginal cost of producing
$w/A$
 is the marginal cost of producing 
 $Y(i)$
. The amount of monopolistic profit is
$Y(i)$
. The amount of monopolistic profit is
 \begin{equation} \pi (i)=p(i)A[l_{Y}(i)-\delta ]-wl_{Y}(i)=\frac{1-\varepsilon }{\varepsilon }w\left [ l_{Y}(i)-\frac{\delta }{1-\varepsilon }\right ] \text{,} \end{equation}
\begin{equation} \pi (i)=p(i)A[l_{Y}(i)-\delta ]-wl_{Y}(i)=\frac{1-\varepsilon }{\varepsilon }w\left [ l_{Y}(i)-\frac{\delta }{1-\varepsilon }\right ] \text{,} \end{equation}
which is positive if and only if 
 $l_{Y}(i)=\overline{l}_{Y}N\gt \delta/(1-\varepsilon )$
 for all
$l_{Y}(i)=\overline{l}_{Y}N\gt \delta/(1-\varepsilon )$
 for all 
 $i\in \lbrack 0,1]$
. As before, due to the fixed cost
$i\in \lbrack 0,1]$
. As before, due to the fixed cost 
 $\delta$
, the industrial market would not operate unless population
$\delta$
, the industrial market would not operate unless population 
 $N$
 is sufficiently large.
$N$
 is sufficiently large.
 To be consistent with our baseline model, we assume that agents produce their own food output 
 $x$
 for their own consumption and raising children. However, they need to purchase industrial output
$x$
 for their own consumption and raising children. However, they need to purchase industrial output 
 $y$
 (when available) using their industrial labor income
$y$
 (when available) using their industrial labor income 
 $wl_{Y}$
. Therefore, in the industrial era, each agent maximizes
$wl_{Y}$
. Therefore, in the industrial era, each agent maximizes 
 $x+y=f+y$
 subject to farming production in (4), labor constraint
$x+y=f+y$
 subject to farming production in (4), labor constraint 
 $l_{F}+l_{Y}=l$
, and the following budget constraint:
$l_{F}+l_{Y}=l$
, and the following budget constraint:
 \begin{equation} py=wl_{Y}+\frac{1}{N}\int _{0}^{1}\pi (i)di\text{,} \end{equation}
\begin{equation} py=wl_{Y}+\frac{1}{N}\int _{0}^{1}\pi (i)di\text{,} \end{equation}
where profits 
 $\pi (i)\geq 0$
 are redistributed to all
$\pi (i)\geq 0$
 are redistributed to all 
 $N$
 agents equally. The first-order condition is
$N$
 agents equally. The first-order condition is
 \begin{equation} \frac{\partial (x+y)}{\partial l_{F}}=\underset{\equiv MPL_{F}}{\underbrace{\alpha \varphi (l_{F})^{\alpha -1}\left ( \frac{Z}{N}\right ) ^{1-\alpha }}}-\frac{w}{p}\text{,} \end{equation}
\begin{equation} \frac{\partial (x+y)}{\partial l_{F}}=\underset{\equiv MPL_{F}}{\underbrace{\alpha \varphi (l_{F})^{\alpha -1}\left ( \frac{Z}{N}\right ) ^{1-\alpha }}}-\frac{w}{p}\text{,} \end{equation}
where 
 $w/p=w/p(i)=\varepsilon A$
 from symmetry and markup pricing in (A5). Figure A1 plots (A8) and shows that there are two scenarios: (a) interior solution (i.e.,
$w/p=w/p(i)=\varepsilon A$
 from symmetry and markup pricing in (A5). Figure A1 plots (A8) and shows that there are two scenarios: (a) interior solution (i.e., 
 $\varphi \rho ^{1-\alpha }\gt \varepsilon A$
) and (b) corner solution (i.e.,
$\varphi \rho ^{1-\alpha }\gt \varepsilon A$
) and (b) corner solution (i.e., 
 $\varphi \rho ^{1-\alpha }\lt \varepsilon A$
). Recall that we have only assumed
$\varphi \rho ^{1-\alpha }\lt \varepsilon A$
). Recall that we have only assumed 
 $\varphi \rho ^{1-\alpha }\lt A$
 but
$\varphi \rho ^{1-\alpha }\lt A$
 but 
 $\varepsilon \lt 1$
.
$\varepsilon \lt 1$
.

Figure A1. Labor market. Figure plots (A8) and shows that there can be a gradual transition (i.e., 
 $l_{F}>0$
) or an immediate transition (i.e.,
$l_{F}>0$
) or an immediate transition (i.e., 
 $l_{F}>0$
) from agriculture to industrial production.
$l_{F}>0$
) from agriculture to industrial production.
 
Interior solution: If 
 $\varphi \rho ^{1-\alpha }\gt \varepsilon A$
, then the equilibrium level of agricultural labor
$\varphi \rho ^{1-\alpha }\gt \varepsilon A$
, then the equilibrium level of agricultural labor 
 $l_{F}$
 from (A8) is
$l_{F}$
 from (A8) is
 \begin{equation} l_{F}=\left ( \frac{\alpha \varphi }{\varepsilon A}\right ) ^{1/(1-\alpha )}\frac{Z}{N}\text{,} \end{equation}
\begin{equation} l_{F}=\left ( \frac{\alpha \varphi }{\varepsilon A}\right ) ^{1/(1-\alpha )}\frac{Z}{N}\text{,} \end{equation}
which implies that the equilibrium level of industrial labor is
 \begin{equation} l_{Y}=l-l_{F}=l-\left ( \frac{\alpha \varphi }{\varepsilon A}\right ) ^{1/(1-\alpha )}\frac{Z}{N}\text{.} \end{equation}
\begin{equation} l_{Y}=l-l_{F}=l-\left ( \frac{\alpha \varphi }{\varepsilon A}\right ) ^{1/(1-\alpha )}\frac{Z}{N}\text{.} \end{equation}
An industrial market would only emerge if 
 $N$
 is sufficiently large to cover the fixed cost
$N$
 is sufficiently large to cover the fixed cost 
 $\delta$
 such that
$\delta$
 such that 
 $l_{Y}N\geq \delta/(1-\varepsilon )$
, which is required for nonnegative profit
$l_{Y}N\geq \delta/(1-\varepsilon )$
, which is required for nonnegative profit 
 $\pi (i)\geq 0$
. Then, (A10) yields
$\pi (i)\geq 0$
. Then, (A10) yields
 \begin{equation} N\geq \frac{1}{l}\left [ \left ( \frac{\alpha \varphi }{\varepsilon A}\right ) ^{1/(1-\alpha )}Z+\frac{\delta }{1-\varepsilon }\right ] \equiv N_{I}(\underset{+}{\varphi },\underset{+}{Z},\underset{+}{\delta },\underset{-}{A},\underset{-}{l})\text{,} \end{equation}
\begin{equation} N\geq \frac{1}{l}\left [ \left ( \frac{\alpha \varphi }{\varepsilon A}\right ) ^{1/(1-\alpha )}Z+\frac{\delta }{1-\varepsilon }\right ] \equiv N_{I}(\underset{+}{\varphi },\underset{+}{Z},\underset{+}{\delta },\underset{-}{A},\underset{-}{l})\text{,} \end{equation}
which is now given by a closed-form solution and has the same comparative statics as (17).
 Before the emergence of industrial production, the population growth rate 
 $\Delta N_{t}/N_{t}$
 and the steady-state population level
$\Delta N_{t}/N_{t}$
 and the steady-state population level 
 $N_{A}^{\ast }$
 in the agricultural era are given by (34) and (35) in Section 3.3. If
$N_{A}^{\ast }$
 in the agricultural era are given by (34) and (35) in Section 3.3. If 
 $N_{t}$
 reaches
$N_{t}$
 reaches 
 $N_{A}^{\ast }$
 before reaching
$N_{A}^{\ast }$
 before reaching 
 $N_{I}$
, then the economy would remain as an agricultural society indefinitely. From (35) and (A11), the inequality
$N_{I}$
, then the economy would remain as an agricultural society indefinitely. From (35) and (A11), the inequality 
 $N_{A}^{\ast }\gt N_{I}$
 is equivalent to
$N_{A}^{\ast }\gt N_{I}$
 is equivalent to
 \begin{equation} \left ( \frac{\sigma l}{\beta }\right ) ^{1/(1-\alpha )}\gt \left ( \frac{\alpha }{\varepsilon A}\right ) ^{1/(1-\alpha )}+\frac{\delta }{(1-\varepsilon )\varphi ^{1/(1-\alpha )}Z}\text{,} \end{equation}
\begin{equation} \left ( \frac{\sigma l}{\beta }\right ) ^{1/(1-\alpha )}\gt \left ( \frac{\alpha }{\varepsilon A}\right ) ^{1/(1-\alpha )}+\frac{\delta }{(1-\varepsilon )\varphi ^{1/(1-\alpha )}Z}\text{,} \end{equation}
which shows that the gradual transition from an agricultural economy to an industrial economy begins under the following conditions: a low fertility cost 
 $\beta$
, a strong fertility preference
$\beta$
, a strong fertility preference 
 $\sigma$
, a high level of agricultural productivity
$\sigma$
, a high level of agricultural productivity 
 $\varphi$
, a high level of labor supply
$\varphi$
, a high level of labor supply 
 $l$
, a large amount of land
$l$
, a large amount of land 
 $Z$
, a high level of industrial productivity
$Z$
, a high level of industrial productivity 
 $A$
, and a low fixed cost
$A$
, and a low fixed cost 
 $\delta$
 for operating industrial firms. These conditions are the same as in Section 3.4, except that the transition in this case is gradual (i.e.,
$\delta$
 for operating industrial firms. These conditions are the same as in Section 3.4, except that the transition in this case is gradual (i.e., 
 $l_{F}\gt 0$
) until
$l_{F}\gt 0$
) until 
 $N_{t}\rightarrow$
$N_{t}\rightarrow$
 
 $\infty$
.
$\infty$
.
Under the interior solution, the level of output per capita in the industrial era is given by
 \begin{equation} x+y=f+y=\varphi (l_{F})^{\alpha }\left ( \frac{Z}{N_{t}}\right ) ^{1-\alpha }+A\left ( l_{Y}-\frac{\delta }{N_{t}}\right ) =\varphi (l_{F})^{\alpha }\left ( \frac{Z}{N_{t}}\right ) ^{1-\alpha }-Al_{F}+A\left ( l-\frac{\delta }{N_{t}}\right ) \text{,} \end{equation}
\begin{equation} x+y=f+y=\varphi (l_{F})^{\alpha }\left ( \frac{Z}{N_{t}}\right ) ^{1-\alpha }+A\left ( l_{Y}-\frac{\delta }{N_{t}}\right ) =\varphi (l_{F})^{\alpha }\left ( \frac{Z}{N_{t}}\right ) ^{1-\alpha }-Al_{F}+A\left ( l-\frac{\delta }{N_{t}}\right ) \text{,} \end{equation}
which is decreasing in 
 $l_{F}$
 because
$l_{F}$
 because 
 $\alpha \varphi (l_{F})^{\alpha -1}\left ( Z/N\right ) ^{1-\alpha }=w/p=\varepsilon A\lt A$
. Then, (A9) shows that
$\alpha \varphi (l_{F})^{\alpha -1}\left ( Z/N\right ) ^{1-\alpha }=w/p=\varepsilon A\lt A$
. Then, (A9) shows that 
 $l_{F}$
 is decreasing in
$l_{F}$
 is decreasing in 
 $N$
. Substituting (A9) and (A13) into (26) yields the population growth rate, which as before converges toward the same steady state
$N$
. Substituting (A9) and (A13) into (26) yields the population growth rate, which as before converges toward the same steady state 
 $\Delta N/N=\frac{\sigma }{\beta }Al-1$
 as
$\Delta N/N=\frac{\sigma }{\beta }Al-1$
 as 
 $N_{t}\rightarrow \infty$
.
$N_{t}\rightarrow \infty$
.
 
Corner solution: If 
 $\varphi \rho ^{1-\alpha }\lt \varepsilon A$
, then the level of industrial labor
$\varphi \rho ^{1-\alpha }\lt \varepsilon A$
, then the level of industrial labor 
 $l_{Y}$
 increases sharply from 0 to
$l_{Y}$
 increases sharply from 0 to 
 $l$
 when
$l$
 when 
 $N_{t}$
 crosses the threshold
$N_{t}$
 crosses the threshold 
 $N_{I}\equiv \delta/[(1-\varepsilon )l]$
. In this case, the inequality
$N_{I}\equiv \delta/[(1-\varepsilon )l]$
. In this case, the inequality 
 $N_{A}^{\ast }\gt N_{I}$
 is equivalent to
$N_{A}^{\ast }\gt N_{I}$
 is equivalent to
 \begin{equation} (1-\varepsilon )\left ( \frac{\sigma }{\beta }\varphi l\right ) ^{1/(1-\alpha )}\frac{Z}{\delta }\gt 1\text{,} \end{equation}
\begin{equation} (1-\varepsilon )\left ( \frac{\sigma }{\beta }\varphi l\right ) ^{1/(1-\alpha )}\frac{Z}{\delta }\gt 1\text{,} \end{equation}
which uses (35) and has the same comparative statics for 
 $\{\beta,\sigma,\varphi,l,Z,\delta \}$
 as in Section 3.4. The only exception is industrial productivity
$\{\beta,\sigma,\varphi,l,Z,\delta \}$
 as in Section 3.4. The only exception is industrial productivity 
 $A$
; however, a larger
$A$
; however, a larger 
 $A$
 makes the corner solution more likely to apply in which case industrialization could be triggered as a result because the threshold
$A$
 makes the corner solution more likely to apply in which case industrialization could be triggered as a result because the threshold 
 $N_{I}$
 decreases from (A11) to
$N_{I}$
 decreases from (A11) to 
 $N_{I}\equiv \delta/[(1-\varepsilon )l]$
.
$N_{I}\equiv \delta/[(1-\varepsilon )l]$
.
 It is useful to note that although the industrial transition is immediate in this case, 
 $N_{I}\equiv \delta/[(1-\varepsilon )l]$
 is not the same as
$N_{I}\equiv \delta/[(1-\varepsilon )l]$
 is not the same as 
 $N_{I}$
 in (16)–(17) and that there exists a unique interior value of
$N_{I}$
 in (16)–(17) and that there exists a unique interior value of 
 $\varepsilon \in (0,1)$
 above which
$\varepsilon \in (0,1)$
 above which 
 $\delta/[(1-\varepsilon )l]$
 is greater than
$\delta/[(1-\varepsilon )l]$
 is greater than 
 $N_{I}$
 in (16)–(17) in which case industrialization occurs later because the markup ratio
$N_{I}$
 in (16)–(17) in which case industrialization occurs later because the markup ratio 
 $1/\varepsilon$
 is too small to cover the fixed cost
$1/\varepsilon$
 is too small to cover the fixed cost 
 $\delta$
. Finally, under the corner solution, the level of output per capita and the population growth rate in the industrial era are the same as (20) and (36), respectively. In the long run, the population growth rate rises toward the same steady state
$\delta$
. Finally, under the corner solution, the level of output per capita and the population growth rate in the industrial era are the same as (20) and (36), respectively. In the long run, the population growth rate rises toward the same steady state 
 $\Delta N/N=\frac{\sigma }{\beta }Al-1$
 as
$\Delta N/N=\frac{\sigma }{\beta }Al-1$
 as 
 $N_{t}\rightarrow \infty$
.
$N_{t}\rightarrow \infty$
.
Appendix B: Robustness check
Table B1. Separate regression 
 $y_{j}={\kappa _{3}\varphi }_{j}+\Phi _{j}+\epsilon _{3,j}$
$y_{j}={\kappa _{3}\varphi }_{j}+\Phi _{j}+\epsilon _{3,j}$

 Notes: Robust standard errors in parentheses. 
 $^{\ast}$
$^{\ast}$
 
 $p\lt 0.1,^{\ast \ast}$
$p\lt 0.1,^{\ast \ast}$
 
 $p\lt 0.05,^{\ast \ast \ast}$
$p\lt 0.05,^{\ast \ast \ast}$
 
 $p\lt 0.01$
. Here, we estimate the effects of agricultural productivity on the share of non-agricultural employment in 1991 in column 1, the non-agricultural share of GDP in 1980 in column 2, and the log of GDP per capita in 1950 in column 3.
$p\lt 0.01$
. Here, we estimate the effects of agricultural productivity on the share of non-agricultural employment in 1991 in column 1, the non-agricultural share of GDP in 1980 in column 2, and the log of GDP per capita in 1950 in column 3.
 
  
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 











