ABSTRACT
The paper presents a framework for determining the optimal size of local jurisdictions and whether it varies depending on the geographical heterogeneity of the territory. To that aim, it first develops a theoretical model of cost efficiency that takes into account spatial interactions and spillover effects among neighbouring jurisdictions. The model solution leads to a spatial Durbin panel data specification of local spending as a non-linear function of population size. The model is tested using a large local data set over the period 2003–11 for an aggregate measure of public spending. The empirical findings suggest a ‘U’-shaped relationship between population size and the costs of providing public services. A second step investigates the role of geographical characteristics such as elevation and terrain ruggedness in the determination of the optimal jurisdiction size. The results reveal that optimal city size decreases with elevation and increases with ruggedness.
ACKNOWLEDGEMENTS
The authors gratefully acknowledge the associate editor and three anonymous referees for their numerous insightful comments and constructive suggestions that helped to improve the work significantly. They also thank the participants at the 3rd International Conference on ‘Decentralization after the Great Recession: Fine-tuning or Paradigm Change?’, Santiago de Compostela, Spain, October 2017, for very helpful comments.
DISCLOSURE STATEMENT
No potential conflict of interest was reported by the authors.
Notes
1. See Slack and Bird (Citation2012) and Tavares (Citation2018) for a review of the empirical literature on the existence of economies of scale, the efficiency gains and the democratic outcomes of consolidated local governments.
2. For papers focusing on the optimal size of a coalition and its characteristics, see, for example, Alesina and Spolaore (Citation1997, Citation2005), Blume and Blume (Citation2007), Bolton and Roland (Citation1997), Gordon and Knight (Citation2009), Saarimaa and Tukiainen (Citation2010) and Hanes (Citation2015).
3. Unlike Solé-Ollé and Bosch (Citation2005), where the omission of relevant spatial interaction terms could lead to bias/inconsistent and inefficient estimates (Elhorst, Citation2014; LeSage & Pace, Citation2009), we estimate a panel data version of the SDM that allows one to quantify with accuracy the magnitude of spatial spillovers, thus minimizing the possibilities of overstating or underestimating the optimal municipal size in this context. On the other hand, the analysis displays two key advantages when compared with the cross-sectional analysis based on the general nesting spatial model (GNSM) of Hortas-Rico and Salinas (Citation2014). First, unlike the GNSM, the SDM does not suffer from parameter identification issues (Elhorst, Citation2014).
4. The empirical analysis carried out here presents some similarities with respect to Bastida et al. (Citation2013) and Rios et al. (Citation2017), as they also employ spatial econometric techniques to model local government spending. However, they do not focus on the link between population and spending. Rios et al. (Citation2017) control for population density but do not include population; Bastida et al. (Citation2013) do not consider the possibility of a non-linear effect.
5. This function F is assumed to satisfy the following assumptions: (1) F is continuous, twice differentiable, positive, displays diminishing marginal products and constant returns to scale (i.e., ); and (2) F satisfies the Inada conditions.
6. Note that:
given that
7. is assumed to be the constant across municipalities owing to the perfect mobility of capital.
8. The signs of the first and second partial derivatives of depend on the sign and relative magnitude of the population parameters.
9. For an alternative approach to investigate non-linear effects in panels with large T while accounting for cross-sectional heterogeneity, cross-sectional dependence and feedback effects among the dependent and independent variables, see Chudik, Mohaddes, Pesaran, and Raissi (Citation2017).
10. Accordingly, the recent reform of the local administration (27/2013 Act of Rationalization and Sustainability of Local Administration) establishes measures to encourage the voluntary merger of municipalities and shifting services of municipalities of fewer than 20,000 inhabitants upwards to the provincial councils. Recent examples of territorial reforms in Europe with a similar spirit can be found in Denmark, Greece, Germany and Switzerland. For a more detailed review, see Steiner, Kaiser, and Eythorsson (Citation2016).
11. Given that some of these control variables are not the main objective of the present study, they are discussed here only in brief. See Ladd and Yinger (Citation1989) and Ladd (Citation1992, Citation1994) for a review of arguments that justify their inclusion in the local spending model.
12. Specifically, public provision is compulsory for all municipalities in services such as rubbish collection, street cleaning, water supply, sewerage and street lighting, among others. Municipalities with a population greater than 5000 inhabitants additionally have to provide parks, public libraries and solid waste treatment. Municipalities with a population greater than 20,000 have to provide local police and social services. Finally, municipalities with a population greater than 50,000 inhabitants also have to provide public transport and environmental protection.
13. However, at least in Spain, this does not necessarily translate into different spending levels because local government tends to provide services even without explicit official responsibility. It is ultimately citizens’ demands (and lack of intervention by other layers of government) that help to explain why a service is provided. Thus, the relationship between the level of responsibility and local spending might not be as evident as it might seem at first (Solé-Ollé & Bosch, Citation2005).
14. This implies assuming perfect mobility of labour among the municipalities within the same province.
15. These variables have been provided by Goerlich and Cantarino (Citation2010). According to them, the inclusion of both variables is justified on the grounds that they are almost completely unrelated and both exhibit a very different spatial distribution.
16. These estimated parameters do not refer to the cost function but to the expenditure function. Nonetheless, cost effects can be recovered by dividing the total effects by , where is the estimated total effect of the tax share.
17. The equation implied by our estimates is given by the total effects: .
18. This procedure makes inference slightly more complex than that of linear models with interaction effects or with spatial effects, but without interaction terms. The reason is that in a spatial model with interaction terms, the researcher cannot directly draw from the variance–covariance matrix of the parameters and the estimated parameters but, instead, the draws have to be taken from .
19. To calculate the size of the flypaper effect from our log–log regression and make it comparable with previous studies, elasticities are transformed using:
where and denote the total effects obtained with the log–log specification for transfers and income, respectively; and and are the ratio of spending to grants and income, respectively.
20. According to our data, about 70% of the municipalities with fewer than 1000 inhabitants were included in the high-altitude subsample.
21. Available data show that the number of population clusters is higher in rugged locations, but, at the same time, the average distance between the population clusters within those municipalities is lower.
22. The quality of the results obtained when estimating the pooled model is less than that obtained by the fixed-effects model (R2 = 50.8% versus 89.5%), and, therefore, they should only be taken as indicative. In fact, the likelihood ratio tests on the joint significance of the fixed effects is LR = 76,824.60 (d.f. = 5521), with p = 0.00, suggesting that the simplified model should not be employed to perform inference. Obviously, the fact that in this setting we are not including spatial fixed effects precludes one from performing any comparison with the implied N* reported in .
23. In particular, we estimate alternative specifications of the following augmented SDM version:where N denotes the logarithm of the population; G is the key geographical variable of the regression model interacting with N; and X includes the set of previously defined demand and cost factors (including population density and the other geographical variable). Solving for N in allows one to obtain the optimal population size N* as a function of G. The simulation of the effects for elevation and terrain ruggedness are reported in Tables A6 and A7 in Appendix A in the supplemental data online, whereas the median forecasts of N* stemming from the effects of models (1) to (6) in Tables A6 and A7 are shown in Figures A2 and A3 (online), respectively.