Forest Transition and Metropolitan Transformations in Developed Countries: Interpreting Apparent and Latent Dynamics with Local Regression Models

: Metropolitan fringes in Southern Europe preserve, under different territorial contexts, natural habitats, relict woodlands, and mixed agro-forest systems acting as a sink of biodiversity and ecosystem services in ecologically vulnerable landscapes. Clarifying territorial and socioeconomic processes that underlie land-use change in metropolitan regions is relevant for forest conservation policies. At the same time, long-term dynamics of fringe forests in the northern Mediterranean basin have been demonstrated to be rather mixed, with deforestation up to the 1950s and a subsequent recovery more evident in recent decades. The present study makes use of Forest Transition Theory (FTT) to examine spatial processes of forest loss and expansion in metropolitan Rome, Central Italy, through local regressions elaborating two diachronic land-use maps that span more than 80 years representative of different socioeconomic and ecological conditions. Our study evaluates the turnaround from net forest area loss to net forest area gain, considering together the predictions of the FTT and those of the City Life Cycle (CLC) theory that provides a classical description of the functioning of metropolitan cycles. The empirical ﬁndings of our study document a moderate increase in forest cover depending on the forestation of previously abandoned cropland as a consequence of tighter levels of land protection. Natural and human-driven expansion of small and isolated forest nuclei along fringe land was demonstrated to fuel a polycentric expansion of woodlands. The results of a Geographically Weighted Regression (GWR) reveal the importance of metropolitan growth in long-term forest expansion. Forest–urban dynamics reﬂect together settlement sprawl and increased forest disturbance. The contemporary expansion of fringe residential settlements and peri-urban forests into relict agricultural landscapes claims for a renewed land management that may reconnect town planning, reducing the intrinsic risks associated with fringe woodlands (e.g., wildﬁres) with environmental policies preserving the ecological functionality of diversiﬁed agro-forest systems.


Introduction
Providing a global interpretation of socio-environmental changes across a development gradient, Forest Transition Theory (FTT) predicts the inherent shift from net forest area loss, typical of emerging economies, to net forest area expansion [1]. This shift is characteristic of advanced economies and is coupled with a generalized recovery of mixed agricultural-wildland ecosystems [2][3][4]. This process has occurred in recent times via natural regeneration [5], active planting [6][7][8], or a combination of the two [9]. From this perspective, the notion of 'forest transition' is associated intimately with underlying socioeconomic forces [10], such as urbanization, late industrialization, tourism growth, The predictions of the FTT have been tested empirically at various spatial scales [17][18][19][20][21]. Although the main forces of deforestation act locally, forest transitions in developed countries are far less investigated at that scale, especially in metropolitan regions [22][23][24]. As a result of urbanization, these regions represent a challenging opportunity to investigate the response of agro-forest landscapes to ecological disturbance and socioeconomic shocks [15,[25][26][27]. Continuous urbanization has reflected the long-term interactions between nature and humans characteristic of the Mediterranean region [7,15,20,[25][26][27][28][29]. Following urbanization, the abandonment of cropland, clearcutting, but also wood re-colonization in some cases, may result in sequential stages of forest decline and recovery, depending on the local context [30][31][32]. Wildland-Urban Interfaces (WUIs) in Southern Europe frequently mix low-density settlements with natural habitats [17]. As a result, fringe landscapes feature irregular boundaries with discontinuous, low-density settlements that reflect suburbanization (i.e., population relocation to suburbs) and the delocalization of economic activities in peripheral areas [33]. These fragmented and mixed landscapes derive from a latent (and quite recent) process of peri-urban forest expansion into cropland, shrubland, and pastures [34], which has followed a significant decline in fringe forest cover following compact urban expansion [35]. As a matter of fact, economic development was demonstrated to fuel, in Italy as in Spain, and in Portugal as in Greece, compact urban growth in the 1950s and the 1960s and the uneven expansion of residential, low-density settlements in the 1970s and the 1980s [36][37][38][39]. From this perspective, FTT The predictions of the FTT have been tested empirically at various spatial scales [17][18][19][20][21]. Although the main forces of deforestation act locally, forest transitions in developed countries are far less investigated at that scale, especially in metropolitan regions [22][23][24]. As a result of urbanization, these regions represent a challenging opportunity to investigate the response of agro-forest landscapes to ecological disturbance and socioeconomic shocks [15,[25][26][27]. Continuous urbanization has reflected the long-term interactions between nature and humans characteristic of the Mediterranean region [7,15,20,[25][26][27][28][29]. Following urbanization, the abandonment of cropland, clearcutting, but also wood recolonization in some cases, may result in sequential stages of forest decline and recovery, depending on the local context [30][31][32]. Wildland-Urban Interfaces (WUIs) in Southern Europe frequently mix low-density settlements with natural habitats [17]. As a result, fringe landscapes feature irregular boundaries with discontinuous, low-density settlements that reflect suburbanization (i.e., population relocation to suburbs) and the delocalization of economic activities in peripheral areas [33]. These fragmented and mixed landscapes derive from a latent (and quite recent) process of peri-urban forest expansion into cropland, shrubland, and pastures [34], which has followed a significant decline in fringe forest cover following compact urban expansion [35]. As a matter of fact, economic development was demonstrated to fuel, in Italy as in Spain, and in Portugal as in Greece, compact urban growth in the 1950s and the 1960s and the uneven expansion of residential, low-density settlements in the 1970s and the 1980s [36][37][38][39]. From this perspective, FTT may appropriately predict cyclical trends in forest decline and expansion in complex peri-urban landscapes including relict agro-forest systems [40]. At the same time, predictions of the FTT were compared with those delineated by the City Life Cycle (CLC) theory. This theory provided a classical description of metropolitan cycles in advanced economies [41], explaining-at least indirectly-the socioeconomic dynamics at the base of peri-urban forest decline and Land 2022, 11, 12 3 of 14 expansion [38]. CLC indicates different, sequential stages of urban growth forming a cycle of urbanization, suburbanization, counter-urbanization, and re-urbanization [33]. The first two stages, urbanization and suburbanization, have been particularly well studied in Mediterranean Europe [16], and they respectively reflect a compact-dense/mono-centric growth and a low-density, spatially heterogeneous expansion of residential settlements [42].
For the first time, to the best of our knowledge, the present study mixes these two dynamic theories (CLC and FTT) with the aims of (i) improving forest assessment in Mediterranean fringe districts and (ii) informing dedicated planning strategies for the sustainable management of WUIs [2,18,43,44]. Mediterranean forests are widely recognized as an invaluable sink of biodiversity [45], forming green infrastructures that may (indirectly) contain urbanization and the negative externalities of economic activity [34]. However, earlier studies have demonstrated that built-up areas around forests expanded in parallel with a decrease in cropland. In other words, residential settlements and forests act as competing land uses with negative impacts on agricultural systems [5].
Going beyond a model of forest expansion from small nuclei of pristine, high-quality woodlands (concentrated in peripheral districts), recent land-use dynamics reflect a more complex path with the growth of scattered and fragmented forest nuclei embedded in a consolidated network of relict peri-urban and rural woodlands [31,46,47]. The shift from a traditional model of forestation irradiating from a few peripheral locations (hereafter defined as 'mono-centric') to a 'polycentric' structure of forest growth (irradiating from multiple locations even close to cities), may be seen at the base of FTT in Mediterranean contexts [48]. This pattern may inform a refined interpretation of forest transitions in metropolitan regions [39]. This 'polycentric' network of forests is also at the base of green infrastructures, assuring the efficient protection of natural habitats in fringe districts [48].
The present study introduces a refined analysis of long-term changes in forest cover as a result of urban expansion in a large metropolitan area (Rome, Central Italy) representative of compact, mono-centric cities in Southern Europe [33]. Reflecting socio-ecological dynamics that impact ecologically fragile Mediterranean environments, land-use complexity warrants further investigation of forest trends  in light of the FTT and the parallel urban cycle. A particularly long time interval was investigated in this study covering a productive cycle of forests and a period that encompassed various stages of metropolitan growth including urbanization and suburbanization [35]. Peri-urban landscapes surrounding Rome (Central Italy) have been identified as the most dynamic across the country because of city size and human pressure [34], resulting in sequential waves of forest decline and recovery [47]. While 55 Natura 2000 sites involving 22% of the study area have been recently established in Rome [48], natural land converted to settlements grew by 0.4% per year during the last half-century, endangering fragile coastal ecosystems more than mountain systems [49].
To check the validity of the polycentric model applied to forest expansion, local econometric estimations based on a Geographically Weighted Regression (GWR) contributed to ascertain the relationship between settlement expansion and changes over time in the composition and morphology of fringe forest landscapes [45]. Being deliberately simplified, our approach tests the shift from a mono-centric model of forest distribution over space toward a more polycentric structure, as hypothesized by Colantoni et al. (2015) [35]. The empirical verification of this hypothesis benefits from a spatially explicit analysis of a dependent variable (percent share of forest cover in total landscape area) and of a predictor, the distance from a central location, which is assumed to represent a metropolitan gradient from central to peripheral locations typical of mono-centric regions [41]. To account for the residual variability that the predictor is unable to explain [16], space was explicitly considered in the econometric specification. The empirical results of this study may inform policy and planning strategies preserving relict agro-forest systems in ecologically sensitive contexts typical of Mediterranean metropolises.

Study Area
The investigated area is the metropolitan region of Rome (5355 km 2 ), which is an administrative partition of the Latium region (Central Italy) encompassing the previously established province of Rome (NUTS-3 level of the European Nomenclature of Territorial Statistics). The area includes 15% flat land (<100 m at the sea level) and 20% mountainous land (>600 m at the sea level). The Simbruini mountains (highest elevation: 1820 m above sea level), which belong to the Apennine district, are the most important relief in the area (Figure 2, left); the alluvial flat area of the Tiber river is made up of lowlands (the so-called "Agro Romano" district; Figure 2, right). The climate is typically Mediterranean, with rainfalls concentrated in autumn and spring, and mild winter temperatures [35]. The average annual rainfall and mean daily temperature in Rome (1971Rome ( -2000 were nearly 700 mm and 17 • C, respectively [47].
centric regions [41]. To account for the residual variability that the predictor is unable to explain [16], space was explicitly considered in the econometric specification. The empirical results of this study may inform policy and planning strategies preserving relict agro-forest systems in ecologically sensitive contexts typical of Mediterranean metropolises.

Study Area
The investigated area is the metropolitan region of Rome (5355 km 2 ), which is an administrative partition of the Latium region (Central Italy) encompassing the previously established province of Rome (NUTS-3 level of the European Nomenclature of Territorial Statistics). The area includes 15% flat land (<100 m at the sea level) and 20% mountainous land (>600 m at the sea level). The Simbruini mountains (highest elevation: 1820 m above sea level), which belong to the Apennine district, are the most important relief in the area (Figure 2, left); the alluvial flat area of the Tiber river is made up of lowlands (the so-called "Agro Romano" district; Figure 2, right). The climate is typically Mediterranean, with rainfalls concentrated in autumn and spring, and mild winter temperatures [35]. The average annual rainfall and mean daily temperature in Rome (1971Rome ( -2000 were nearly 700 mm and 17 °C, respectively [47]. The total population almost doubled during the study period from 2 million to 4 million inhabitants [5]. The percent share of the population living in inner Rome in the total population of the study area increased from 20% to 40% between 1936 and 2020 [39]. Even though urban settlements are becoming prevalent in the study area, the majority of the province's land is still comprised of forests, pastures, and cultivated fields [45]. Summer wildfires and human pressure, because of Rome expansion, have degraded the pristine forests, but relict woodlands have been sometimes preserved in both flat and hilly areas [48]. Based on official statistics [49], forest composition was changed significantly in the study area. Chestnuts and conifers increased considerably between the mid-1930s and late 2010s (respectively from 13.8% to 18.5% and from 1.7% to 3.0% of total forest stock); beech experienced a moderate decline (from 10.1% to 9.8% of total forest stock). Cropland abandonment around the central city was the base of natural afforestation uplands [47]; urbanization determined, at the same time, a progressive fragmentation of relict forests around the city, especially along the coastal rim [45]. Industrial areas were traditionally located in the eastern part of the 'Agro Romano' district, while traditional rural areas were (and still are) relatively abundant in the Western side of the study area [49]. The total population almost doubled during the study period from 2 million to 4 million inhabitants [5]. The percent share of the population living in inner Rome in the total population of the study area increased from 20% to 40% between 1936 and 2020 [39]. Even though urban settlements are becoming prevalent in the study area, the majority of the province's land is still comprised of forests, pastures, and cultivated fields [45]. Summer wildfires and human pressure, because of Rome expansion, have degraded the pristine forests, but relict woodlands have been sometimes preserved in both flat and hilly areas [48]. Based on official statistics [49], forest composition was changed significantly in the study area. Chestnuts and conifers increased considerably between the mid-1930s and late 2010s (respectively from 13.8% to 18.5% and from 1.7% to 3.0% of total forest stock); beech experienced a moderate decline (from 10.1% to 9.8% of total forest stock). Cropland abandonment around the central city was the base of natural afforestation uplands [47]; urbanization determined, at the same time, a progressive fragmentation of relict forests around the city, especially along the coastal rim [45]. Industrial areas were traditionally located in the eastern part of the 'Agro Romano' district, while traditional rural areas were (and still are) relatively abundant in the Western side of the study area [49].
Land-use maps were reproduced in Figure 3 to illustrate the mono-centric organization of Rome's metropolitan region, with compact settlements located downtown and a moderate spread of low-density settlements around central locations [34]. The shift toward dispersed urbanization in 2018 was not altering the polarization in central places (downtown Rome) and peripheral districts West and East of Rome [16]. Overall, the study area is administered by 122 municipalities, including the municipality of Rome that encompasses the historical city ('within the Aurelian Walls') and the surrounding area (1285 km 2 ) within and outside the main Ring Road ('Grande Raccordo Anulare', GRA). In this study, the central municipality of Rome was further partitioned in 115 local districts (a sort of formal neighborhood known as 'suddivisioni toponomastiche' in Italian, and it is mainly used for statistical reporting at a sub-municipal scale). Following earlier studies [33], the present study made use of a specific (territorial) partition of the area into 236 elementary spatial units that derive from the union (i.e., spatial merging) of 121 municipalities (excluding Rome) and 115 local districts in Rome.
Land-use maps were reproduced in Figure 3 to illustrate the mono-centric organization of Rome's metropolitan region, with compact settlements located downtown and a moderate spread of low-density settlements around central locations [34]. The shift toward dispersed urbanization in 2018 was not altering the polarization in central places (downtown Rome) and peripheral districts West and East of Rome [16]. Overall, the study area is administered by 122 municipalities, including the municipality of Rome that encompasses the historical city ('within the Aurelian Walls') and the surrounding area (1285 km 2 ) within and outside the main Ring Road ('Grande Raccordo Anulare', GRA). In this study, the central municipality of Rome was further partitioned in 115 local districts (a sort of formal neighborhood known as 'suddivisioni toponomastiche' in Italian, and it is mainly used for statistical reporting at a sub-municipal scale). Following earlier studies [33], the present study made use of a specific (territorial) partition of the area into 236 elementary spatial units that derive from the union (i.e., spatial merging) of 121 municipalities (excluding Rome) and 115 local districts in Rome.

Logical Framework
To test the inherent shift toward a 'polycentric' model of forest expansion (Colantoni et al. 2015), the percent of forest land area in the total landscape (i.e., municipal/district area) was derived, for each elementary analysis' unit, from the spatial overlay of a shapefile map illustrating the boundaries of the selected domains (n = 236 municipalities/local districts, see above) separately with each forest map (1936 and 2018). The average distance of each spatial unit (polygon) from downtown Rome was computed as the linear distance of the polygon centroid from a central place in the historical cityconsidered the heart of political, social, economic, touristic, and cultural life [16]. Adopting administrative boundaries as the elementary unit of investigation gave room to reliable comparisons with official statistics and constitutes a relevant domain for planning purposes [34]. In this study, we assumed Rome's metropolitan region as a laboratory of forest transitions [47], e.g., hypothesizing a classical ('mono-centric') model of forest decline around settlements and forest expansion in peripheral districts, and testing a possible shift toward a more 'polycentric' model of forest expansion around suburban settlements [48]. To verify such assumptions, changes over time in the spatial distribution of forest land along the urban gradient were analyzed comparing the results of spatially implicit and spatially explicit regressions that test the predictions of the FTT at a local scale [5]. More specifically, global and local regressions assessed the empirical relationship between percent forest area (dependent variable) and the distance from downtown Rome (predictor) at two times (1936 and 2018). We assumed a linear and spatially homogeneous relationship between the two variables in a 'mono-centric' model of forest expansion (i.e., a low percentage of forest area close to the inner city that increases with the distance from the city itself). With 'poly-centric' forest expansion [35], the relationship between the two variables is expected to be less polarized along the metropolitan gradient, with the predictor's impact varying largely across space [33]. The highest impact can be reflective of a local context (directly or indirectly) supporting forest expansion [5].

Statistical Analysis
Following the logical framework exposed above, both spatially implicit and explicit models, estimated separately for two time points (1936 and 2018) have been adopted to

Logical Framework
To test the inherent shift toward a 'polycentric' model of forest expansion (Colantoni et al. 2015), the percent of forest land area in the total landscape (i.e., municipal/district area) was derived, for each elementary analysis' unit, from the spatial overlay of a shapefile map illustrating the boundaries of the selected domains (n = 236 municipalities/local districts, see above) separately with each forest map (1936 and 2018). The average distance of each spatial unit (polygon) from downtown Rome was computed as the linear distance of the polygon centroid from a central place in the historical city-considered the heart of political, social, economic, touristic, and cultural life [16]. Adopting administrative boundaries as the elementary unit of investigation gave room to reliable comparisons with official statistics and constitutes a relevant domain for planning purposes [34]. In this study, we assumed Rome's metropolitan region as a laboratory of forest transitions [47], e.g., hypothesizing a classical ('mono-centric') model of forest decline around settlements and forest expansion in peripheral districts, and testing a possible shift toward a more 'polycentric' model of forest expansion around suburban settlements [48]. To verify such assumptions, changes over time in the spatial distribution of forest land along the urban gradient were analyzed comparing the results of spatially implicit and spatially explicit regressions that test the predictions of the FTT at a local scale [5]. More specifically, global and local regressions assessed the empirical relationship between percent forest area (dependent variable) and the distance from downtown Rome (predictor) at two times (1936 and 2018). We assumed a linear and spatially homogeneous relationship between the two variables in a 'mono-centric' model of forest expansion (i.e., a low percentage of forest area close to the inner city that increases with the distance from the city itself). With 'poly-centric' forest expansion [35], the relationship between the two variables is expected to be less polarized along the metropolitan gradient, with the predictor's impact varying largely across space [33]. The highest impact can be reflective of a local context (directly or indirectly) supporting forest expansion [5].

Statistical Analysis
Following the logical framework exposed above, both spatially implicit and explicit models, estimated separately for two time points (1936 and 2018) have been adopted to interpret long-term forest expansion/decline and the underlying territorial dynamics along the metropolitan gradient in Rome. A preliminary check of the linear or non-linear relationship between forest cover and distance from downtown Rome was performed using parametric (Pearson) and non-parametric (Spearman) pair-wise correlation coefficients, testing for significance at p < 0.05. A detailed description of the regression techniques Land 2022, 11, 12 7 of 14 subsequently adopted to test the polycentric expansion of forests has been provided in the following two paragraphs, distinguishing global from local approaches.

Global Approaches
Spatially implicit global approaches adopted in this study included Ordinary Least Square (OLS) regressions using linear and non-linear specifications, namely first (linear), second-order (squared), and third-order (cubic) forms. More specifically, the mono-centric structure was first tested adopting a linear specification [33] estimated through logarithmic transformations of both dependent and predictor variables as follows: where F is the percent share of forest cover in total landscape, D is the distance of each spatial unit from downtown Rome (km, log-transformed), a and b are the regression coefficients, and e is the regression error. In this case, forest cover increases (or decreases) linearly along the metropolitan gradient, satisfying the assumption of the mono-centric model [33].
To verify the possible impact of the varying size (i.e., surface area) of each spatial partition adopted as the elementary analysis' unit in this study (i.e., municipalities/urban districts), we run OLS linear estimations of Equation (1) both un-weighted and weighted with the surface area of each elementary unit (km 2 ). District area was the result of the application of an ArcGIS 'calculate area' tool on the shapefile depicting municipalities/urban districts provided by the National Institute of Statistics (ISTAT) and statistical office of Rome municipality (Figure 2, left).
A squared specification was also tested with the aim at defining a quadratic relationship of the metropolitan gradient on forest cover: where F is the percent share of forest cover in the total landscape, D is the distance of each spatial unit from downtown Rome (km, log-transformed), a, b, and c are the regression coefficients, and e is the regression error. More complex patterns based on a cubic specification were tested to highlight place-specific responses of forest cover as the distance from downtown Rome increases [48], as follows: where F is the percent share of forest cover in the total landscape, D is the distance of each spatial unit from downtown Rome (km, log-transformed), a, b, c, and d are the regression coefficients, and e is the regression error. Significant specifications were assessed using standard diagnostics, namely the adjusted R 2 .

Local Approaches
We adopted a spatially explicit approach referring to the Geographically Weighted Regression (GWR) introduced by Fotheringham et al. (2002) to interpret complex spatial patterns of environmental and socioeconomic processes along a given geographical gradient. A refined understanding of the spatial variability related to the mono-centric expansion of metropolitan regions contributes to identify significant, local deviations from global patterns [48], informing spatial policies better adapted to specific territorial conditions. The use of administrative districts, such as municipalities, taken as policy-relevant spatial domains, allows a comprehensive investigation of context-based interactions that are frequently demised in 'centralized' (spatially implicit) interpretative models [51]. Based on these premises, GWRs estimate local regression models for each spatial domain in a given area, accounting for spatial dependence and heterogeneity [38], i.e., controlling for spatial structures that characterize the relationship between the dependent variable and the  (Ali et al., 2007). We adopted a standard GWR specification for a given location s = 1 to n, as follows: F(s) = B(s)D(s) + e(s) (4) where F(s) and D(s) are, respectively, the dependent variable (percent forest cover) and the predictor (distance from downtown Rome) illustrated in Equation (1) above, which were both measured at each location s; B(s) means the column vector of the regression coefficients at location s, and e(s) is the random error at location s. Regression parameters were estimated separately for 1936 and 2018 at each location by weighted least squares and are thus spatially explicit, i.e., a function of s. Spatial coordinates of each elementary unit considered in this study were determined based on a shapefile map of administrative boundaries (see Figure 2, left) managed in ArcGIS (ESRI Inc., Redwoods, CA, USA) environment. A Kernel spline function was adopted to calculate weights for the estimation of local regressions [38]. Limitations regarding the use of GWR only arise when drawing conclusions based on a reduced number of sample observations, which is not the case of this study. GWRs were run separately for each time interval. For each regression, the model's results included (i) a global measure of goodness-of-fit (adjusted R 2 ) that was compared with the same index obtained from the respective OLS regression (both weighted and un-weighted for the surface area of each spatial partition, see above) and (ii) a spatial distribution of local coefficients (R 2 , slope coefficient) illustrated through maps (Salvati et al. 2016). The GWR results clarify the role of spatial heterogeneity as a novel contribution to the empirical test of the mono-centric model [33]. The intrinsic shift from a classical 'mono-centric' model of forest expansion toward a more 'poly-centric' structure was reflected in a different spatial distribution of both slope and R 2 regression coefficients [16]. In fact, we expect (i) a more fragmented/heterogeneous distribution of R 2 coefficients (possibly indicating the increase in sparse nuclei of forest expansion) and (ii) an expansion of the spatial domains with high slope coefficients, reflecting a latent, but generalized, increase in forest cover [35]. With this perspective in mind, comparing the intensity and direction of the relationship between percent forest cover and the distance from downtown Rome for different econometric models and specifications (e.g., global vs. local, spatially implicit vs. spatially explicit) provides an intrinsic measure of the role of territorial factors [48]. Such contextual aspects should be more explicitly considered in any environmental policy targeting the mitigation of (and adaptation to) urbanization processes [49].

Results
Forest cover in Rome increased from 18.3% in 1936 to 19.9% in 2018. The spatial distribution of forest area ( Figure 5) reflects the intense landscape transformations between 1936 and 2018, with forest expansion in relict croplands East of Rome and moderate decline in flat districts surrounding the 'Agro Romano' and along the sea coast. These districts included, up to now, the main remnants of the forest cover extending along the Tyrrhenian Sea coast from Tuscany to Campania in ancient times. Accelerated rates of forest recovery have been observed in rural, upland districts west and north of Rome, which were originally devoted to shrublands and olive groves (e.g., 'Sabina' district). Forest cover was relatively stable in the 'Castelli Romani' district, a peri-urban, upland area with high human pressure (population density around 500 inhabitants/km 2 ) and mostly devoted to chestnut production.

Global Approaches
All statistical models testing changes in the spatial relationship between forest cover and the distance from downtown Rome (Table 1)  and third-order (i.e., squared and cubic) OLS regressions performed better for both years (adjusted R 2 = 0.58 for 2018). With global R 2 coefficients increasing from 0.65 in 1936 to 0.75 in 2018, local regressions outperformed OLS (linear and non-linear) models, documenting the importance of the spatial dimension in long-term forest transitions characteristic of the study area. Moran's spatial autocorrelation indexes demonstrate that the percent share of forest cover in the total municipal/district area was not randomly distributed over space (1936: I(z) = 4.29, p < 0.05; 2018: I(z) = 3.69, p < 0.05). Spatial autocorrelation justifies the use of a local regression approach when assessing the relationship between percent share of forest cover in the total municipal/district area as the dependent variable and the distance from downtown Rome as an intrinsic test of the mono-centric (or poly-centric) forest structure.

Global Approaches
All statistical models testing changes in the spatial relationship between forest cover and the distance from downtown Rome (Table 1) indicated how the correlation between the two variables intensified over time. Non-parametric Spearman rank coefficients delineate a positive and strong correlation between forest cover and distance from a central place in Rome (rs = 0.76 and 0.86, respectively for 1936 and 2018). However, while the first-order (linear) Ordinary Least Square (OLS) model performed rather low for 1936 and 2018, second-and third-order (i.e., squared and cubic) OLS regressions performed better for both years (adjusted R 2 = 0.58 for 2018). With global R 2 coefficients increasing from 0.65 in 1936 to 0.75 in 2018, local regressions outperformed OLS (linear and nonlinear) models, documenting the importance of the spatial dimension in long-term forest transitions characteristic of the study area. Moran's spatial autocorrelation indexes demonstrate that the percent share of forest cover in the total municipal/district area was not randomly distributed over space (1936: I(z) = 4.29, p < 0.05; 2018: I(z) = 3.69, p < 0.05). Spatial autocorrelation justifies the use of a local regression approach when assessing the relationship between percent share of forest cover in the total municipal/district area as the dependent variable and the distance from downtown Rome as an intrinsic test of the mono-centric (or poly-centric) forest structure.

Local Approaches
The satisfactory result of local regressions suggests how territorial contexts distinctive of municipalities and urban districts in the study area shape forest dynamics at different distances from downtown Rome. Local regressions modeled the spatial dimension of both 1936 and 2018 data by incorporating a location distance matrix that allowed the local estimation of both slope (measuring the predictor's impact) and R 2 (measuring the model's goodness-of-fit) regression coefficients for each spatial domain. Vector maps were used Land 2022, 11, 12 10 of 14 to illustrate the spatial distribution of local regression slope and local R 2 coefficients. The results of local regression models for forest cover as the dependent variable and the distance from downtown Rome as the predictor were illustrated in maps considering separately the spatial distribution of regression slopes and local R 2 goodness-of-fit coefficients ( Figure 6). Slope coefficients depicted the (more or less intense) contribution of specific areas to a mono-centric spatial structure.

Local Approaches
The satisfactory result of local regressions suggests how territorial contexts distinctive of municipalities and urban districts in the study area shape forest dynamics at different distances from downtown Rome. Local regressions modeled the spatial dimension of both 1936 and 2018 data by incorporating a location distance matrix that allowed the local estimation of both slope (measuring the predictor's impact) and R 2 (measuring the model's goodness-of-fit) regression coefficients for each spatial domain. Vector maps were used to illustrate the spatial distribution of local regression slope and local R 2 coefficients. The results of local regression models for forest cover as the dependent variable and the distance from downtown Rome as the predictor were illustrated in maps considering separately the spatial distribution of regression slopes and local R 2 goodness-of-fit coefficients ( Figure 6). Slope coefficients depicted the (more or less intense) contribution of specific areas to a mono-centric spatial structure.  The empirical results of local regressions indicate a distinctive spatial distribution of forest cover for 1936 and 2018, in turn reflecting a substantially different relationship with the distance from a central place in Rome. Comparing the models' goodness-of-fit, the regression for the last observation year (2018) was performing slightly better than the model for the first observation year (1936). Spatial units (municipalities/districts) with local R 2 > 0.5 (regarded as a satisfactory model's fit at that spatial scale) increased during the study period. The spatial distribution of local R 2 coefficients in 1936 delineated particularly high values in economically peripheral municipalities of the Apennine mountain districts ('Lucretili' and 'Lepini' districts), as well as in rural, remote municipalities in the northwest of Rome ('Tolfa-Bracciano' district). In 2018, the highest R 2 coefficients were found either in upland/mountainous municipalities East of Rome ('Lucretili', 'Simbruini' and 'Lepini' districts); the same pattern was observed in sparse coastal/lowland municipalities featuring human disturbance and low ecological quality ('Anzio-Nettuno' coastal district South of Rome; 'Bracciano Lake' district). In both cases, spatial units with high R 2 were found closer to Rome in 2018 than in 1936. As far as the impact of the distance from downtown, the highest regression slopes for 1936 were recorded in remote districts East of Rome and in more sparse flat areas close to the city. High slopes were also recorded in a rural, remote district West of Rome, and they decreased in intensity along the sea coast and in flat districts around the city. Eighty years later, while peripheral municipalities in the Apennine district still had the highest slope coefficients, more heterogeneous, peri-urban districts both west and east of Rome displayed high (or very high) slope coefficients. Taken together, local slope and R 2 coefficients reflect the polarization in districts with forest recovery and decline, which is likely because of wildfires, clearcutting, or other human/ecological disturbances. The gradual disappearance of coastal woodlands because of progressive clearcutting made the spatial distribution of forests more heterogeneous, in turn contributing to a less polarized landscape in urban and rural areas. At the same time, urban sprawl and the rising human pressure driven by sprawled urbanization has consolidated some smaller nuclei of forest expansion, acting as an 'endogenous source' of forest transitions (i.e., fueling the shift from net decline to net recovery) at the local scale.

Discussion
During the last century, a significant decrease in cropland and a contemporary increase in residential settlements and forests has been observed in Italy [50] and, more generally, in Southern Europe [23]. Land-use change was responsible for the expansion of 'contact areas' between developed (i.e., residential and/or industrial settlements) and natural/seminatural landscapes in Mediterranean Europe [24,25,52,53]. By confirming the predictions of the Forest Transition Theory, our study documents the relevance of a comprehensive landscape analysis along urban-forest interfaces to inform policies containing urbanization and mitigating the negative impact of human disturbance on the surrounding natural habitats [7]. Barbati et al. (2013) found that the spatial distribution of natural landscapes in recent times became more polarized along the metropolitan continuum [5]. A thorough understanding of the intrinsic mechanisms of urban growth and the consequent dynamics of natural landscapes-mainly forests-along this gradient provides an integral vision and future perspective to landscape research [48]. From this perspective, our study contributes to regional science suggesting how the notion of mono-centric (or poly-centric) urban expansion can be (more or less directly) extended to dynamics related to other land use, namely forests, underlying socioeconomic mechanisms and territorial processes that are governed (and are in turn influenced by) the intrinsic linkage between humans and nature at the base of regional landscape dynamics [8]. A global analysis failed in assessing these complex dynamics, suggesting how spatially explicit approaches are appropriate tools to investigate the importance of heterogeneous local contexts [33], especially when they assume a stronger role than regional drivers of change and local socioeconomic pressures [9].
The results of local regressions presented in our study provided evidence on the importance of a mono-centric structure of urban growth in Rome, which is in line with the empirical findings of earlier studies [5,16,42]-as plastically illustrated in Figure 2 earlier in this paper with the parallel consolidation of a moderately polycentric model of forest expansion [35]. This complex spatial pattern is assumed to be the joint impact of two processes: (i) intense (natural) forestation in mountainous and economically depressed areas and (ii) moderate (and mostly human-induced) forest recovery in flat and coastal areas, depending on the intensity of conflicts over alternative land-use, i.e., on the basis of the 'economic strength' of the agricultural sector [31,38,45,47]. As a matter of fact, cropland abandonment typical of economically marginal districts was observed also in some periurban areas originally devoted to agriculture because of particularly fertile soils [51]. In such areas, forest expansion occurred more likely than in other socioeconomic contexts of the same geographical region [5,39,48].
Being representative of more general dynamics at the regional scale, distinctive transformations of peri-urban landscapes were described in the study area considering together the empirical results of this study and earlier literature: (i) economic decline of remote districts resulting in re-forestation [28,50,54], (ii) crop intensification in more accessible rural districts [5,29,34], and (iii) urban sprawl, agricultural land consumption, and forest fragmentation in flat districts [35,55,56]. As a counter-intuitive result, forest increase was intrinsically associated with reduced proximity between settlements and natural land [51], as a result of forestation of abandoned cropland in accessible districts indirectly fueled by urban expansion [37][38][39]. From a socioeconomic standpoint, such transformations fueled the heterogeneous expansion of forest land cover [25], with ecosystem dynamics remaining essentially under-investigated [30].
From the spatial planning perspective, the results of our study highlight the need to adopt a multi-target strategy for the conservation of natural and semi-natural areas in peri-urban contexts [35]. This strategy should recognize forest expansion in peri-urban areas together as a positive process of greening [31] and a potentially negative phenomenon, as the increase in forest cover as well as in the length of the margin directly in contact with residential settlements causes an increased risk of wildfires in already vulnerable areas [47]. Future studies should better interpret long-term changes in metropolitan landscapes [37], clearly distinguishing the impact of socioeconomic drivers and planning/biophysical constraints to building, in turn clarifying the possible role of different forest and agroforest types in that change [5,45,48]. More generally, the analysis of changes in the forest landscape in peri-urban areas [46] cannot be separated from an accurate investigation of the causes of expansion and decline of agriculture in the facing rural contexts.

Conclusions
The present study refers to long-term changes in forest area that encompass a sufficiently long time period covering the productive cycle of a high forest. Long-term changes in forest cover include both natural and human-induced transformations and may account for a truly comprehensive picture of landscape changes. In this perspective, we tried to estimate the synergic impact of both drivers. The complex results of this study (i.e., the emergence of a polycentric model of forest expansion in a context of mono-centric expansion of urban settlements) outlines the urgent need for policy strategies managing urban-forest interactions along the fringe. Effective strategies governing peri-urban landscape dynamics should include (i) urban containment policies reducing the environmental impact of sprawl on fringe land and (ii) measures controlling the unwanted expansion of peri-urban woodlands (starting from the multiplication of a sort of 'hotspots' indirectly acting as nuclei of forest spread) at the expense of abandoned cropland.
Based on the empirical findings of this study and a pertinent review of the recent literature, these actions should target the integral preservation of the traditional diversity of agro-forest mosaics in Southern European fringe lands. Our study highlights the importance of monitoring techniques for the quantification of urban-forest ecosystem services, starting from a recent field experience in the Mediterranean context. New logical approaches and information technologies for the 'holistic' study of peri-urban forests are requested to integrate operational tools, methodologies, and analysis' techniques, with increasing, comparative evidence from case studies. Assessing the quality and health of peri-urban forests implies a reduction of monitoring costs, the improvement of sampling efficiency and accuracy, and a more effective integration of the existing measurement networks in a long-term vision that refers to the paradigms of ecological resilience and socioeconomic sustainability.