Unraveling population trends in Italy (1921–2021) with spatial econometrics

Testing density-dependence and path-dependence in long-term population dynamics under differentiated local contexts contributes to delineate the changing role of socioeconomic forces at the base of regional disparities. Despite a millenary settlement history, such issue has been rarely investigated in Europe, and especially in highly divided countries such as those in the Mediterranean region. Using econometric modeling to manage spatial heterogeneity, our study verifies the role of selected drivers of population growth at ten times between 1921 and 2021 in more than 8000 Italian municipalities verifying density-dependent and path-dependent dynamics. Results of global and quantile (spatial) regressions highlight a differential impact of density and (lagged) population growth on demographic dynamics along the urban cycle in Italy. Being weakly significant in the inter-war period (1921–1951), econometric models totalized a high goodness-of-fit in correspondence with compact urbanization (1951–1981). Model’s fit declined in the following decades (1981–2021) reflecting suburbanization and counter-urbanization. Density-dependence and path-dependence were found significant and, respectively, positive or negative, with compact urbanization, and much less intense with suburbanization and counter-urbanization. A spatial econometric investigation of density-dependent and path-dependent mechanisms of population dynamics provided an original explanation of metropolitan cycles, delineating the evolution of socioeconomic (local) systems along the urban-rural gradient.


Econometric analysis
We assumed different spatial regimes of population dynamics associated with each stage of the cycle 64 , modelling the variability in population growth rates ( Pop.Growth (t=2,1) ) as a function of (i) population growth rate ( Pop.Growth (t=2,1) ) in the previous (i.e.lag (−1) ) time interval, (ii) population density ( Dem.Density (t=0) ), and (iii) the overall size of the resident population as a proxy of agglomeration ( Pop.Size (t=0) ), both measured at the beginning of the related observation time.The analysis has also considered (iv) an average measure of Elevation for each municipality, as derived from the official source of Istat municipal atlas, (v) a dummy of closeness to the sea coastline (Sea.Prox) classifying each municipality as 'coastal' (code 1) or 'inland' (code 0) and (vi) a dummy separating municipalities acting as the 'head town' (Cap.city) of a given province (with code 1) from the remaining municipalities classified with code 0. All variables were standardised prior to analysis 50 .Use of these variables in econometric models testing density-dependent and path-dependent mechanisms of population dynamics was discussed in 42,50,63 .Model specification was summarised as follows: where α is the regression constant (model's intercept), β 1 , β 2 . . . ., β 6 are the regression coefficients (slope), and ǫ is the stochastic error of the model.Models were run controlling for time (i.e.distinguishing the impact of the four stages of the cycle mentioned above) and space (i.e. using spatial econometrics approaches whose results were compared with those from standard approaches).The adopted specification allow (i) discriminating the impact of economic from non-economic forces of population growth, (ii) distinguishing the role of densitydependent mechanisms of population growth and decline from the more general path-dependency of local population dynamics, (iii) highlighting the importance of direct spatial effects, and (iv) separating them from the indirect ones (i.e.spillovers).The individual stages of the metropolitan cycle in Italy were defined as follows 63 : urbanisation (1951-1981), suburbanisation (1981-2001), counter-urbanisation mixed with early re-urbanisation (2001-2021); population dynamics during inter-war decades were more mixed and prepared the system to the sudden shift toward compact urbanisation.Cross-section regressions were run assuming each observation decade as a separate stage of the cycle representative of specific socioeconomic contexts and population dynamics at the local scale.A Variance Inflation Factor (VIF) was finally calculated for each time interval.Values systematically below 5 for all variables delineate a non-redundant structure of predictors' matrix, in line with the basic assumption of non-collinearity typical of most econometric models.

Standard models
Assuming linear changes over time in population distribution over space, Eq. (1) was preliminary tested with a linear specification adopting global Ordinary Least Squares (OLS) regressions.The models' goodness of fit was checked by way of adjusted R 2 coefficients and the Akaike Information Criterion (AIC).Inference on regression results (i.e.Fisher-Snedecor F tests and Student t tests respectively on the overall regression fit and on individual coefficients, testing against the null hypothesis of zero coefficients with p < 0.001 ) provided an additional cri- terion for model's evaluation 74 .To verify the violations of the basic assumptions of a general linear model, a Durbin-Watson (DW) statistic checking for serial correlation, a Breusch-Pagan (BP) index for heteroscedasticity, and a Moran (M) spatial autocorrelation coefficient for spatial dependence of residuals were run for each model, testing for significance at p < 0.05 against the null hypothesis of no serial correlation, no heteroscedasticity, and no spatial autocorrelation structure, respectively.

Spatially explicit models
Equation (1) was additionally estimated comparing the results of global models that make spatial relations explicit using spatial weights among municipalities calculated as (i) a contiguity ( 0 − 1 ) Queen matrix (Q) and (ii) a linear distance matrix (W).While presenting a variable goodness-of-fit, consistent regression outputs (i.e. the same significant predictors with comparable intensity and sign) may identify a statistically stable (and conceptually relevant) relationship between population growth rates and the selected predictors 42 .By investigating the dependence of a given variable's values on the values of the same variable recorded at neighbouring locations, spatial autocorrelation assumes outcome in one area to be affected by outcomes, covariates or errors in nearby areas, meaning that models may contain spatial lags of the outcome variable, spatial lag of covariates, and autoregressive errors, respectively 56 .Regressions run in this study include a Spatial Autoregressive Model (SAR), a spatial autoregressive error term (SDE), and a Spatial Durbin Model (SDM).A flowchart of the adopted spatial models has been shown in Figure 1.
Both direct and indirect (spillover) effects between municipalities were detected.Best-fit estimation of the proposed models using empirical data was evaluated using pseudo R 2 .As in the case above, non-parametric (quantile) regressions with spatial weights separately from Q and W matrices were run to estimate Eq. 1 at four percentiles of the dependent variable assuming (i) deviation from normality, (ii) non-linear dependence between predictors and the dependent variable, and (iii) spatial relations among input variables.Model's outcomes include estimates of intercept and slope coefficients and the associated significance level testing for the null-hypothesis of non-significant regression coefficient at p < 0.05.

The interwar period
Ordinary Least Square (OLS) regressions performed on standardized input variables show positive and significant values of lagged population growth rates and proximity to the sea coast despite a relatively small adjusted-R 2 (1) Pop.Growth (t=2,1) = α + β 1 Pop.Growth (t=1,0) + β 2 Dem.Density (t=0) + β 3 Pop.Size (t=0) + β 4 Elevation+ + β 5 Sea.Prox + β 6 Cap.city + ǫ (0.05) typical of the time period (1921-1951) preceding the sharp wave of urbanisation in the aftermath of World War II.The regression coefficient for head towns is also positive.Negative and significant coefficients are observed for elevation.The variance inflation factor (VIF) is systematically low and largely below 5, suggesting a nonredundant structure of the predictor matrix (Table 1).The results obtained from quantile regressions confirm the results of the OLS regression with reference to lagged population growth, population size, elevation, proximity to the sea coast and head towns.In contrast, the differences with the OLS regression emerge when referring to population density.In general, the slope coefficients of lagged population growth rates increase from the first ( τ = 0.25 ) to the fourth ( τ = 0.99 ) quartile.Starting from the coefficients of the variable linked to population density, the sign-in the first three quartiles-is consistent with the one performed in the OLS, while in the fourth quartile is negative.A significant and positive coefficient is also shown by the variable linked to the size of the population from the first to the third quartile, to then become negative-even if significant-if referred to the fourth quartile.With reference to the last two variables, i.e. proximity to the sea and head town, the coefficient maintains the sign of OLS model with highly significant values in the quartiles from τ = 0.25 to τ = 0.75 for the first variable, and no significance in τ = 0.99 for the first predictor.With reference to the second predictor, no coefficient appears significant except that of the third and forth quartiles, maintaining a coherent coefficient sign with OLS model.Compared to OLS models, lower AIC values have been shown in the first three quartiles of the standard quantile regression; a higher value is shown when τ = 0.99.
The results of the econometric tests indicate OLS estimations as partly biased, since tests for serial correlation, heteroscedasticity, and spatial dependence are all significant.Thus, we adopted spatial econometric techniques starting from a Spatial Durbin model (hence SDM) for the global regressions by explicitly considering the spatial structure of the input data and we compared the results with reference techniques such as Spatially Autoregressive (SAR), and Spatial Error (SDE) models.In direct and indirect SDM regressions we find positive coefficients related to the variable describing lagged population growth rates, coherently with the OLS regression.This predictor displays significant and positive coefficients for both direct and indirect effects.With reference to population size and density, elevation, as well as head town, we observe opposite signs considering direct and indirect effects.For instance, the indirect effect of population size is negative and statistically significant, while being positive but statistically insignificant when considering direct effects.Proximity to the sea coast has a positive impact (both direct and indirect) on population growth rates.Moving on to spatially quantile regressions, significant and negative coefficients were found for elevation -even if the intensity of this significance is strongest in the first three quartiles and then decreases in the fourth quartile.Positive coefficients, on average, are associated with lagged population growth and demographic density, proximity to the sea coast, and head town.Finally, the positive role of head town on population growth is positive and significant at least in the second, third and forth quartile.Positive coefficients are associated with population size in the first and second quartile, reverting to a negative coefficient in the forth quartile.Similar results were found comparing the outcomes of quantile regressions run with spatial matrices based on contiguity and linear distances among municipalities.Models based on the distance-weighted matrix tend to have lower AIC values than models based on contiguity.www.nature.com/scientificreports/These results are in line with R 2 , which, in this case, shows greater values than contiguity-based spatial models.
In particular, SDM performs the lowest value of AIC whereas SAR performs the greatest value of AIC in distance weighted models.When considering contiguity, SDE and SDM produce the lowest AIC values, while SAR performs the greatest AIC value.Generally speaking, the goodness of fit of spatial models is consistently higher than that of standard models.

The breakdown of intense urbanisation
Moving to the subsequent time interval, Table 2 illustrates the results of standard econometrics and spatial models using both contiguity and linear distance spatial weights.Here the results of OLS estimates and standard quantile regressions appear rather coherent in assigning positive and significant coefficients to lagged population growth rates, demographic density, proximity to the sea coast and head town, with increasing values of the adjusted-R 2 .Econometric diagnostics also indicate that the OLS estimate is (moderately) biased.The tests of serial correlation, heteroscedasticity and spatial dependence are all significant, suggesting the appropriateness of using spatial models.When comparing the econometric results from different spatial weighting schemes, models such as SAR and SDE give similar results with OLS as far as sign and significance of the regression coefficients.A general comparison of OLS values of AIC with that of spatial weighted models highlight that SAR and SDM seem to be more effective than OLS in explaining the variance of the phenomenon, while SDE has an     1951), elevation, proximity to the sea coast and a dummy indicating municipalities that act as provincial head town as predictors (*p < 0.05 ; **0.001 < p < 0.05 ; ***p < 0.001).
In general, regressions based on a spatially weighted distance matrix seem to perform better than those based on contiguity.With distance, the three AIC values are consistent with the corresponding R 2 (lower in SAR and greater in SDM).With contiguity, SDM displays a lowest AIC value that corresponds to the greater R 2 value.Considering quantile regressions, lagged population growth, elevation, and proximity to the sea coast maintain the OLS coefficients' sign and significance almost for all quartiles.Population size coefficients were found positive for the first three quartiles, becoming negative in the fourth quartile.
Considering a subsequent time interval with intense expansion of human settlements in Italy, lagged population growth rates, demographic density, elevation, proximity to the sea coast, and head town display almost positive and significant regression coefficients both for OLS and quantile regressions, with markedly improved R 2 (Table 3).The positive coefficients estimated for demographic density and population size are antithetical to those observed for earlier time windows.Variance inflation factors (VIFs) for each predictor are systematically less than 5, suggesting a non-redundant structure of the regressors matrix.Test diagnostics for serial correlation, heteroscedasticity, and spatial dependence suggest the appropriateness of using spatial modelling.Considering both contiguity and distance weighting, SAR and SDE results are rather consistent with those of OLS with respect to the sign of the regression coefficients.The direct effects of SDM mostly resemble the outcome of SAR and SDE models.Direct and indirect coefficients for demographic density and population size assume opposite signs.Spatial models, in general, have higher values of AIC than OLS.Quantile regressions provided similar outcomes irrespective of the spatial weighting scheme adopted.

From compact urbanisation to suburbanisation
Table 4 illustrates the results of econometric modeling investigating a distinctive time interval as far as population dynamics in Italy is concerned.Both OLS and quantile regressions improved their goodness-of-fit (adjusted-R 2 > 0.4 ), corresponding to contained values of AIC, at least in OLS and in the first three quartiles of standard quartile regressions.Lagged population growth rates and demographic density have a positive impact of current population growth rates; population size, elevation, proximity to the sea coast and head town have a negative-while less intense and more mixed-impact on current population growth rates.Following the results of econometric diagnostics, spatial models provide outcomes adjusted to the spatial structure of both predictors and the dependent variable.However, in the essence, global spatial models give results well aligned with those of the OLS regression; moreover, spatial quantile regressions confirm the outcomes of standard quantile regressions.Coherent results are also displayed comparing adjusted-R 2 with AIC, i.e., higher values of the first indicator correspond to lower values of the second one.Adjusted-R 2 confirms the best performance of SDM in both cases, while low performances (low adjusted-R 2 and high AIC), are found in the two specifications of SAR.
Results of both OLS and quantile regressions for the subsequent time interval (Table 5) are also aligned with the models' outcomes described for the previous decade (see above), with model's goodness-of-fit maintaining generally high.Lagged population growth rates, demographic density, and proximity to the sea coast reveal a positive and significant impact on the dependent variable.Elevation and head town show the reverse effect, while population size display more mixed results, possible associated with a higher spatial heterogeneity characteristic of this predictor.Although the non-redundant structure of predictors (VIF systematically below 2), econometric diagnostics suggested the appropriateness of spatial modelling also in this case.Improving slightly the goodness-of-fit in respect with non-spatial modeling, results of global regressions (SAR, SDE and SDM, direct effects) are mostly aligned with the OLS model.Moreover, spatial quantile regressions estimate a comparable structure of coefficients for the selected predictors, irrespective of the spatial weighting scheme.All these models show appreciable goodness-of-fit that reflects convergent results, as far as sign and significance of regression coefficients.Based on the values of AIC, the best performance is attributed to SDE and SDM weighted by distance In Table 5.Conversely, SAR model exhibits a lower value of AIC with contiguity.All outcomes are in line with adjusted-R 2 .

Moving toward counter-urbanisation
Similarly to what reported in Tables 4 and 5, modeling population growth rates between 1981 and 1991 (Table 6) with standard econometrics delineate the role of head town (negative impact) as well as lagged population growth rates and demographic density (positive impacts).In partial disagreement with what has been observed in earlier decades, population size, elevation, and proximity to the sea coast, show more heterogeneous results, with significant and non-significant coefficients and contrasting signs when moving from OLS to quantile regressions.Despite a satisfactory goodness-of-fit, econometric diagnostics suggested the use of spatial models that reached a systematically higher fit than standard models.Results of global models, irrespective of the spatial weighting scheme, confirm the role of lagged population growth rates and demographic density as predictors of current population growth rates.Outcomes from AIC, confirmed by adjusted-R 2 , follow the same tendency.Elevation and population size give more mixed, while significant, results.Quantile regressions provide useful insights when describing the latent heterogeneity associated with predictors such as population size, elevation, proximity to the sea coast and head town.For instance, the negative impact of population size is spatially polarized, i.e. higher for the forth quartile and lower moving from the third to the first quartile.A reverse pattern is observed for elevation and proximity to the sea coast.
Standard and spatial modeling explaining the variability in population growth rates (1991-2001) gain significance, with a considerable goodness-of-fit increasing further when the spatial structure of predictors is considered (Table 7).Lagged population growth rates and demographic density show a positive and highly significant coefficient irrespective of the model used.The impact of population size (negative and significant coefficients) is also homogeneous across model's specifications.The impact of elevation is also negative, although minor differences were found comparing the results of quantile regressions run with the two spatial weighting schemes.
Vol:.( 1234567890 Global models (both standard and spatial) also suggest how proximity to the sea coast exerted a negative and significant impact on the dependent variable.Quantile regressions provide similar results, with the only exception of the forth quartile.Head town as a predictor of population growth rates is associated to mostly negative coefficients (global models).Outcomes of adjusted-R 2 and AIC, jointly confirm the better performance of distance weighted models compared to contiguity.Negative coefficients were also found in quantile regressions (both standard and spatial) for the first and the second quartiles, but not for the third and the forth quartiles.

Latent trends toward re-urbanisation
Table 8 shows the results of econometric models estimating the spatial variability of population growth rates between 2001 and 2011.All models, starting from OLS, significantly improved their goodness-of-fit compared with what has been estimated in earlier decades.Econometric diagnostics, goodness-of-fit (adjusted-R 2 or pseudo-R 2 ) and the Akaike Information Criterion (AIC) coherently document how spatial models performed better than standard models.Concerning spatial models, in line with adjusted-R 2 , the distance weighted scenario seems to better perform than contiguity.Lagged population growth rates and demographic density have a positive and highly significant impact on the dependent variable, irrespective of the model's specification and the spatial weighting scheme.Global and quantile models provide comparable results for these two predictors.
The impact of elevation is, in turn, negative and highly significant, for all models.However, the indirect impact of this variable estimated via SDM results slightly significant or completely insignificant.The impact of the proximity to the sea coast and head town as predictors of population growth rates is negative -as depicted in the outcomes of global (both standard and spatial) econometric models -and more mixed in quantile regressions, being moderately significant for the first and second quartiles (proximity to the sea coast) and for the first and forth quartiles (head town).The outcomes of regressions modeling the spatial variability of population growth rates between 2011 and 2021 indicate substantially different demographic dynamics, as far as intensity and spatial structure are concerned (Table 9).The overall estimate of model's goodness-of-fit is moderate or rather low for all econometric specifications, and improved slightly moving from standard to spatial models, despite test diagnostics were convergent in suggesting the appropriateness of using spatially explicit approaches.From the results of AIC and adjusted-R 2 in spatially weighted models, it emerges the higher fit of distance models in comparison with contiguity models.Quantile regressions performed better than global models; more specifically, results for the first and second quartiles had a satisfactory goodness-of-fit, declining for the third and forth quartiles.Among predictors, the impact of lagged population growth rates on the dependent variable is confirmed, in line with the results of estimates for previous decades, although (global) regression coefficients are less intense and significant, irrespective of the model's specification.The same (positive) impact was found in quantile regressions (first to third quartile), with the exception of the forth quartile.The impact of demographic density is really mixed and heterogeneous moving from global to quantile regressions.Population size seems to have a slightly negative impact on the dependent variable, although econometric estimates were rather mixed, as in the case of demographic density.Elevation, proximity to the sea coast and head town display almost insignificant and close-to-zero regression coefficients, apart from few exceptions (mainly in quantile regressions).

Discussion
The present study illustrates a diachronic analysis of demographic dynamics at the municipal scale in Italy, verifying density-dependence, path-dependence, agglomeration/scale impacts, and spatial effects over a complete (demographic-urban) cycle from urbanisation to re-urbanisation.The approach developed in this study was based on the comparison of different statistical models, both parametric (global econometrics) and nonparametric (quantile regressions).This approach allowed a precise identification of the factors underlying processes of demographic growth and decline, in turn consolidating complex urban-rural hierarchies in advanced economies 75 .While global models provide a gross assessment of the impact of various predictors of population growth rate at a sufficiently detailed spatial scale 76 , quantile regressions allow an even more accurate inspection of the trends characteristic of specific parts of the statistical distribution of the dependent variable 77 .These statistical loci correspond to specific demographic behaviours, which reflect diversified but internally homogeneous territorial contexts.Examples include demographically dynamic contexts with positive and sustained growth rates, and demographically shrinking contexts having systematically negative growth rates.While non-spatial models (both global and quantile) provided the baseline knowledge to a refined understanding of population Table 5. Results of standard (OLS, Ordinary Least Square, and quantile) regressions as well as global (SAR: Spatial Autoregressive model; SDE: Spatial Error model; SDM: Spatial Durbin model) and quantile spatial models run with both contiguity and linear distance spatial weighting matrices; population growth rate (% annual) in 1971-1981 as dependent variable; population growth rate (1961-1971), demographic density (1961), population size (1961), elevation, proximity to the sea coast and a dummy indicating municipalities that act as provincial head town as predictors (*p < 0.05 ; **0.001 < p < 0.05 ; ***p < 0.001).www.nature.com/scientificreports/dynamics and the underlying factors and contexts, spatial models proved to be innovative tools analysing regional variability in the dependent variable.Moreover, the comparison between the results of global and quantile models specifying the geographical structure of the elementary units adopted in this study allows a more accurate examination of the role of spatial heterogeneity in population dynamics 78 .As far as the case study, while global models have satisfactorily explained population growth rates at a sufficiently detailed spatial level, the outcomes of quantile models often were in line with the results of global models.These outcomes were characteristic of the 'urbanisation' phase (1951-1961, 1961-1971, 1971-1981) with medium-high growth rates in Italy 79 .On the contrary, during the ' counter-urbanisation' and 're-urbanisation' waves (1991-2001, 2001-2011, 2011-2021), global models fitted the dependent variable less effectively.At such times, quantile regressions provided likely more accurate indications of demographic behaviours in conditions of spatial heterogeneity 80 .This may highlight latent spatial patterns that are characteristic of territories with systematically high or low population growth rates 64 , in turn corresponding with specific quartiles of the statistical distribution of the dependent variable.Taken together, the results of the econometric models document how density-dependence has been observed in correspondence with urbanisation, suggesting a role for economic agglomeration and immigration 24 .Densitydependence was less significant over both suburbanisation and counter-urbanisation, when population tends to be more dispersed across regions 44 .In these contexts, the role of agglomeration and scale reduced proportionally, and path-dependent factors regulating population growth took the lead.With re-urbanisation, the positive rate of population growth observed in rural areas counterbalanced the stable (or negative) pattern observed in urban areas, indicating a relationship with population growth that reflects congestion externalities and subtle processes of peri-urbanisation intensifying in recent decades 5 .The outcomes of quantile regressions have more specifically delineated the existence of a non-linear relationship between population growth and density for all time intervals, although with important differences as far as the impact of individual factors is concerned.Results of quantile regressions document a positive effect of density on population growth rates, being stronger at higher levels of urban concentration, while declining slightly over time.Such findings are in line with the documented outcomes of sequential waves of urbanization, suburbanization and re-urbanization typical of postwar Italy 11,49,53 .In other words, the density-growth relationship is indicative of sequential stages characteristic of the metropolitan cycle in Mediterranean Europe 73 .All in all, our study demonstrates how sequential waves of concentration and de-concentration of urban and rural locations were associated with density-dependent mechanisms of population growth and decline.This process, shaping the expansion of rural/accessible districts, and the abandonment of marginal districts, accentuated the divide in high-density and low-density areas 63,81,82 .The flexibility of this approach justifies an extended use of global econometric models and quantile regressions analysing the drivers of population growth in socioeconomic contexts distinct from the one studied in this application, and at vastly differentiated (spatial and temporal) scales.An extensive use of spatial panel techniques applied to both global and quantile models is also recommended when sufficiently long time series of predictors are available at homogeneous and stable spatial units 77 .In addition, comparing the results of econometric models specifying multiple spatial weighting schemes (e.g. based on different contiguity and distance metrics) seems to be an appropriate tool 83 when inferring about the stability of regression coefficients (sign and significance) across models 84 and when taking decisions about the best performing models 85 .In this last case, diagnostics such as Akaike Information Criterion can easily and effectively complement such approaches 86 .Future studies should also investigate the appropriateness of local techniques (e.g. the Geographically Weighted Regression, GWR) to the investigation of regional variability and local heterogeneity in population dynamics 87 , considering both cross-section approaches and panel extensions, when input data allow such improvements.Although the process of demographic growth and decline at a local scale seem to follow similar underlying logics in various regions of Europe (for instance, in the Mediterranean countries 88 ), a more accurate examination of the latent factors at the base of formation (and/or consolidation) of regional disparities in the geographical distribution and density of resident population may support a spatial planning aimed at a balanced, polycentric, and sustainable development of territories and local communities 89 .While identifying distinctive (demographic) regimes at the local scale, the empirical results of our study outline the intrinsic characteristics of local contexts and the differences in the relationship between population growth and density over time 15 , corroborating the assumption that density-dependent regulation was intrinsically associated with exogenous dynamics depending on urban cycles.In this perspective, long-term demographic processes in Mediterranean Europe 59 can be seen as representative of more general dynamics at the continental scale.Based on a comparative approach 76 , our study definitely offers an exploratory econometric perspective to regional studies of population dynamics that can be easily adapted to different spatial scales (from local to regional levels), temporal schedule (from decadal to annual windows), and variables, e.g.moving from strictly demographic indicators to economic predictors of leading and lagging contexts in the old continent and beyond.

Conclusions
A spatial econometric investigation of density-dependent and path-dependent mechanisms of population dynamics provided an original explanation of metropolitan cycles, delineating the evolution of socioeconomic (local) systems along the urban-rural gradient.With accelerated population dynamics, empirical results delineate compact urbanisation (1951-1981) as the main factor consolidating spatial disparities in Italy.As a matter of fact, econometric models-being only weakly significant in the inter-war period (1921-1951) -showed a high goodness-of-fit in correspondence with compact urbanisation that declined moderately with suburbanisation and counter-urbanisation (1981-2021).Density-dependence and path-dependence were found significant and, respectively, positive and negative, with compact urbanisation, and much less intense with suburbanisation and counter-urbanisation.The results of our study justify a renewed (diachronic and spatially explicit) analysis of

Figure 1 .
Figure 1.Graphical description of the spatial models adopted in the work.

Table 4 .
Results of standard (OLS, Ordinary Least Square, and quantile) regressions as well as global (SAR: Spatial Autoregressive model; SDE: Spatial Error model; SDM: Spatial Durbin model) and quantile spatial models run with both contiguity and linear distance spatial weighting matrices; population growth rate (% annual) in 1961-1971 as dependent variable; population growth rate (1951-1961), demographic density (1951), population size (