Whose Urban Green? Mapping and Classifying Public and Private Green Spaces in Padua for Spatial Planning Policies

The rising environmental issues on contemporary cities urgently calls for sustainable planning policies. Implementation of nature-based solutions, ecosystem services, and green infrastructures associated to green spaces management is at present of paramount importance. In contrast to policies mainly focused on public greenery, the inclusion of private green in planning strategies might be a promising pathway. The general aim is mapping and classifying urban green spaces in Padua, a city of 93.3 km2 (Northeast Italy). Specific aims are (i) testing an NDVI-derived extraction from very high-resolution orthophotos; (ii) classifying property status; (iii) highlighting multilevel relationships and strategies for urban green spaces implementation and management; (iv) assessing greenery in relation to per capita population. By performing remote sensing and GIS analyses, a first detailed global map of urban green spaces in Padua was created; then, binary classification and thematic maps for rural/non-rural, public/private, municipal/non-municipal greenery were produced for all urban units. Results show that, among total green spaces (52.23 km2), more than half are rural. Moreover, private green spaces represent 80%, while within public areas (20%) less than 10% are municipal (5 km2). We therefore highlight scenarios for planning policies in Padua by providing tools to policymakers for an integrated management of green spaces, where private greenery might also contribute to ecosystem services implementation for common urban well-being.


Green Spaces Management in Contemporary Urban Planning
In contemporary territorial and urban planning, nature-based policies and practices to improve the health and social life of citizens are growing and becoming even more relevant. Two terms largely used to encompass planning principles and initiatives related to this framework are ecosystem services and nature-based solutions (NBS) [1][2][3][4][5].
Ecosystem services are inclusive approaches for natural resources management [6][7][8], orienting environmental protection measures and biodiversity conservation to preserve nature and recognise social and economic benefits provided by ecosystems to human societies [9,10]. In particular, urban ecosystem services (UES) are the benefits provided to people living in cities by the environmental and nature-related urban elements [11][12][13].
NBS is an umbrella term including all planning and management strategies trying to cope with climatic, environmental, and socioeconomic issues, including adaptation dynamics [14,15]. In this conceptual framework, the development of policies and initiatives is driven and inspired by ecological processes and functions. In fact, NBS are well recognised within the innovation policies by the European Union (EU) as a supporting tool to foster the goals of the European Green Deal, EU Biodiversity Strategy, and EU Adaptation Strategy [16].
A concept strictly related to UES and NBS is green infrastructure (GI): it can be defined as a network of natural or semi-natural areas preserving biodiversity and providing several types of benefits and services to people [17]. GI typically presents a multiscale dimension, connecting urban and non-urban landscapes as well as supporting sustainable mobility networks; since water bodies such as rivers and streams are often part of the core of such infrastructures, they are also referred to with the expression 'green-and-blue infrastructure' [18].
Including UES and NBS in urban planning strategies is a relevant topic in current planning studies, especially those related to performance-based planning [19][20][21][22]. This research area aims to cope with the complexity of contemporary cities and territories by integrating land use zoning and metric standards for public services with a broader consideration for spatial and environmental resources and their performance in achieving collective health and equality goals. Adopting such a framework, a quantitative assessment of UES may be used to drive qualitative maps of territorial values in support of a local GI strategy [23].
In all the above-mentioned frameworks, urban green spaces (UGS) play an important role in addressing and supporting urban planning policies and actions. They provide relevant UES, especially by the contribution of vegetation systems, permeable soils, and city farmlands. Therefore, they may represent a supporting basis to implement integrated NBS able to mitigate air and water pollution, hydrogeological/hydraulic risk, and urban heat islands; improving their spatial continuity and connections leads to the creation of an effective urban GI [24]. Moreover, UGS play a widely recognised role in improving public health [25].
Traditionally, urban planners and policymakers focused on the development and management of public green spaces [26]. In recent years, the rise of new planning strategies involving a wider role for private actors, together with new approaches such as UES, led to take into consideration the role of private green spaces as a crucial issue of sustainable urban management [26][27][28][29].
Current literature highlights the important role of private gardens for network connectivity and biodiversity conservation, as well as their role as hubs in the urban GI and their contribution in providing regulating and supporting UES such as carbon storage and sequestration, urban heat island mitigation, stormwater attenuation, and noise reduction [26,30,31]. For these reasons, low-density housing with private green spaces is considered by some scholars as a form of sustainable urban planning [32]. However, other studies warn about economic and social issues related to the proliferation and prevalence of private green spaces, especially in situations of lacking or inadequate public UGS [33]. In addition, it is still debated whether urban and peri-urban agricultural plots should be considered as a part of urban GI [34,35]. Indeed, from both GI and UES perspectives, they appear to play a remarkable role [36]: rural greenery is responsible for provisioning regulatory ecosystem services to citizens and may be included, together with other green-and-blue urban features, in urban GI supplied with sustainable mobility networks connecting cities to their natural or semi-natural surroundings.
In this context, accurate localisation and quantification of UGS and mapping their property status become paramount for GI/NBS implementation and sustainable spatial planning.
Detecting locations, extension, and use of green spaces within a city, together with their spatial relationships with urban fabric and infrastructures, is indeed pivotal to the development of GI and to the implementation of site-specific NBS [31]; moreover, data about land property applied to UGS are necessary to urban planners in order to elaborate ISPRS Int. J. Geo-Inf. 2021, 10, 538 3 of 28 feasible, effective strategies to enhance benefits provided by the urban environment. For example, it is possible to integrate the role of private or public soils and vegetation in contributing to the provision of some regulation UES by creating environmental corridors or buffer zones [30,32], while focusing on soft mobility connections between public and accessible UGS in implementing GI for public fruition.

Spatial Planning in Padua: Between Soil Sealing and Urban Greening
Padua is a medium-sized city located in Veneto, northeast of Italy, approximately 35 km west of Venice ( Figure 1). The municipal territory spans 93.3 km 2 with a population of about 211,000 inhabitants (2020). The city lies on an almost completely plain surface, south of the first pre-alpine reliefs. It is an ancient city, playing an important cultural and economic role in Italy since the Middle Ages and characterised by a wall system realised in the 16th century, when the city was under Venetian rule [37]. Lamentably, detailed estimation of the spatial dimension of UGS at the urban scale is disaggregated and incomplete, making spatial planning at present critical and inconsistent. Moreover, identification of suitable areas for soil sealing mitigation/compensation and integration of GI into UES management require urban professionals to also map the property status by classifying public and private UGS.

Aims of the Study
The general aim of this study is to map and classify all UGS within the municipal territory of Padua by integrating NDVI-based analyses with ancillary and ground data.
The specific aims are (i) to test the integration of NDVI extraction from very highresolution orthophotos, together with a municipal geodatabase, to accurately estimate the spatial dimension of UGS; (ii) to classify property status and land use of UGS; (iii) to highlight multilevel relationships and possible strategies for UGS implementation and management; (iv) to assess the distribution and the amount of U in relation to per capita pop- As in other European cities, in recent decades, Padua experimented with urban development consisting of medium-density residential districts around the city centre and sparse new buildings or complexes spreading on the urban fringes between the core city and the surrounding countryside [38]. This phenomenon and the resulting growth of grey mobility infrastructures led to a relevant fragmentation of peri-urban natural and rural areas in Padua [39].
Regarding administrative units, Padua is currently divided into 6 neighbourhoods and 40 sub-neighbourhood urban units [40].
In this urban development context, the latest National Report from the Italian Institute for Environmental Protection and Research ISPRA (2020) ranked Veneto as the second region in Italy affected by soil sealing phenomenon, with 217,619 ha completely covered by sealed surfaces (11.87%). In addition, with 45-50% of sealed surfaces Padua is one of the most affected cities in Italy, ranked as fifth by ISPRA assessment and showing a remarkable rate increase of about 20-25 ha per year [41][42][43]. To limit soil sealing and its impacts on UES, different regional rules and regulations are currently issued for a more sustainable territory planning agenda, involving public institutions, policymakers, and stakeholders. In fact, within the 'no net land take by 2050,' European Commission's Roadmap to a Resource Efficient Europe [44], recent revisions of the Veneto regional planning law issued strict limitations to new buildings.
The regional planning law of Veneto, known as 'Regulations for territorial and landscape management', approved in 2004 [45] and lastly revised in 2017 and 2019, rules that every municipality must adopt an urban development plan consisting of two parts: (i) a 'Territorial Management Plan' ('Piano di Assetto del Territorio', PAT) which includes strategies and policies and (ii) a 'Plan of Interventions' ('Piano degli Interventi', PI) which is based on operational provisions for the implementation of planning decisions [46,47].
By revisions and integrations of the regional planning law [48,49], new regulations and incentives to limit and mitigate soil sealing are issued for all the municipalities of Veneto Region, by promoting regeneration and renaturalisation practices. They designate and allocate the total amount of soil sealing for all the municipalities of the Veneto Region and by implementing GI, renaturalisation and urban regeneration processes [48]. In 2019, a regional decree reallocated for the whole municipal territory of Padua a maximum amount of 262.48 ha of soil sealing [50], fuelling the debate about the protection of existing green spaces, implementation of GIs, and increase of UES.
Lamentably, detailed estimation of the spatial dimension of UGS at the urban scale is disaggregated and incomplete, making spatial planning at present critical and inconsistent. Moreover, identification of suitable areas for soil sealing mitigation/compensation and integration of GI into UES management require urban professionals to also map the property status by classifying public and private UGS.

Aims of the Study
The general aim of this study is to map and classify all UGS within the municipal territory of Padua by integrating NDVI-based analyses with ancillary and ground data.
The specific aims are (i) to test the integration of NDVI extraction from very highresolution orthophotos, together with a municipal geodatabase, to accurately estimate the spatial dimension of UGS; (ii) to classify property status and land use of UGS; (iii) to highlight multilevel relationships and possible strategies for UGS implementation and management; (iv) to assess the distribution and the amount of U in relation to per capita population.
Coherently with our purposes, in this paper, we use a definition of UGS which includes all the pervious surfaces able to host some kind of vegetation and to play a role in green management policies by local authorities. In particular, given the conformation of the municipality territory and the importance of rural land plots in a site-specific analysis focused on Padua, we chose to include peri-urban agricultural areas in our UGS mapping, as proposed in the scientific literature focusing on UGS assessment for urban planning purposes [31,34,35].

Data Sources
For the first screening and detection of UGS, we adopted the Normalised Difference Vegetation Index (NDVI) to perform a multispectral analysis on very high-resolution aerial images. This vegetation index, based on spectral signature [51] is largely used in vegetation mapping [52,53].
NDVI analysis was performed on public orthophotos dating to the early summer of 2015, provided by AGEA/Veneto Region. Orthophotos were structured in two multiband, image datasets at very high spatial resolution (20 cm pixel −1 ): one with red, green, and blue bands in the visible range and the other one with near-infrared (NIR), red, and blue bands.
In order to integrate the results, we used ancillary data from the Topographic Database (DB) of Padua Municipality, part of a series of digital land use/land cover maps implemented by Veneto Region in collaboration with municipalities involved.
For data validation, we used UGS maps of four sample urban units in Padua that had been carried out in the framework of a study aiming to calculate the Biotope Area Factor ecological index [43]. This study was developed by adopting a methodology based on highly detailed photointerpretation performed by skilled researchers, supervised and checked by the authors of the study. The source dataset was AGEA/Regione Veneto 2015 orthoimagery, and the same was used for NDVI extraction. Therefore, in our work, we adopted results of the abovementioned paper [43] as macro-area data validation for our results.
The classification of UGS by property was realised by means of a query on cadastral codes of public institutions performed on the cadastral database of Padua Municipality.
Finally, tables and analyses concerning population in the urban units in which the city is subdivided were based on data from Padua 2019 statistical yearbook, released by the municipality [40].
All spatial and spectral analyses were performed by using the open source GIS software QGIS (version 3.x). Table 1 shows the main input spatial data used in the workflow, schematising their use and sources.

Urban Green Spaces Detection
Firstly, orthophotos were merged and clipped on the boundary of the Padua municipality. Then, the NDVI value was calculated for each image pixel, by performing the standard normalised band ratio as follows: Afterwards, different tests on NDVI threshold values were performed in order to extract green areas. As reported in scientific literature, due to different factors, NDVI threshold values to detect vegetation are not univocal, although a value below 0.1 generally represents the absence of greenness [54]. As the NDVI analysis is able to extract only vegetated surfaces we performed different tests in the range of 0.1-0.3 to identify a suitable threshold value according to the aim of the study.
Visual analysis was performed with the help of updated aerial imagery on different spots in the city and carried out for the different thresholds tested, the results of which confirmed that the value of 0.15 was the most suitable for UGS extraction in this case study. In fact, values more than 0.15 allowed the extraction of very low vegetated areas, avoiding the presence of artificially sealed soils from surfaces detected as UGS.
The NDVI-derived raster image was then vectorised into a shapefile. The adoption of multispectral indexes derived from remote sensing to map urban vegetation, indeed, may produce flaws consisting in the underrepresentation of ploughed or fallow land but also in the overrepresentation of tree canopies spreading on streets or other artificial features.
Basing on topographic DB classes, we performed an integration function and a subtraction on the NDVI-derived vector map. In the first case, we gathered topographic DB classes representing potentially vegetated land, such as 'green area', 'agricultural area', 'pasture and fallow', 'vegetation-free area', 'river banks' and added to the map the areas not already detected by NDVI. In the second case, we gathered DB classes related to sealed or permanently non-vegetated land, such as 'street area', 'pedestrian paths area', 'railway', 'building', 'water body' and subtracted the features overlapping our vegetation map. In this way, we fixed the most common misrepresentations. It should be said that such corrections are useful to improve quantitative mapping results but may cause little flaws in detailed visualisation.
At the end of the process, we obtained a polygonal shapefile representing actually or potentially vegetated land (UGS) within the borders of Padua Municipality.
Later, in order to verify correlations between the size of green areas and their distance from the centre of the city, we extracted centroids of polygons representing UGS and created a scatterplot matching UGS extent and distance from the city centre. Since our goal was to analyse UGS dimensions with respect to Padua spatial development, we calculated distances starting from 'Piazza delle Erbe', the main square in Padua and centre of the historical city since the Middle Ages.

Data Validation
To validate data derived from NDVI extraction and topographic DB integrations, we performed a spatial linear regression with the results of detailed mapping of green surfaces in four macro-sample areas of the city (Figure 2), derived from the study described in Section 2.1. The regression was calculated using the GRASS GIS function 'r.regression.line' in QGIS.
The validation process was strengthened by calculating a confusion matrix for each sample area with QGIS semi-automatic classification plugin, to determine how many pixels resulted as 'green' or 'not green' according to the two datasets. Then, we calculated user's and producer's accuracy and kappa coefficients to measure global accordance.
The macro-areas used for our comparison are four urban units such as Brentelle (macro-area 1, 2.60 km 2 ), Forcellini (macro-area 2, 2.66 km 2 ), Sacra Famiglia (macro-area 3, 2.77 km 2 ), and San Lazzaro (macro-area 4, 3.44 km 2 ). In their diversity, they represent a good sample of the relationships between green spaces and urban fabric in Padua.

Urban Green Spaces Definition and Classification
Once obtained a map of the distribution of green areas in Padua, we chose to organise and classify our data by adopting three binary categories, with the purpose of achieving more structured information and coping with the needs of spatial planners and local policymakers. The binary categories of UGS adopted in our classification are the following: Rural-non-rural. It allows defining urban and peri-urban UGS associated with agricultural land use, which may serve as a basis for community gardens projects and provide both regulating and provisioning UES; Public-private. It is a key categorisation, as it identifies areas that, thanks to their public property status, may be included in mobility or leisure urban GI and, in opposition, private UGS suitable to contribute to the creation of urban ecological corridors and the enhancement of regulating UES; Municipal-non-municipal. It defines UGS owned by the Municipality, that is, the primary features of urban green spatial planning and management oriented to improve green accessibility and reduce inequalities in UGS fruition by citizens.
Firstly, we classified rural and non-rural green areas. We used once again the topographic DB, gathering the class representing agricultural areas with a selection of farmlands belonging to other classes of the DB (e.g., 'green areas', 'fallow', 'non-vegetated areas') detected by visual analysis on orthoimagery and then intersecting the aggregated features with the overall UGS map derived from NDVI. The sum of the resulting features corresponds to the share of rural green areas in Padua. The remaining percentage of total green areas is therefore non-rural.

Urban Green Spaces Definition and Classification
Once obtained a map of the distribution of green areas in Padua, we chose to organise and classify our data by adopting three binary categories, with the purpose of achieving more structured information and coping with the needs of spatial planners and local policymakers. The binary categories of UGS adopted in our classification are the following: Rural-non-rural. It allows defining urban and peri-urban UGS associated with agricultural land use, which may serve as a basis for community gardens projects and provide both regulating and provisioning UES; Public-private. It is a key categorisation, as it identifies areas that, thanks to their public property status, may be included in mobility or leisure urban GI and, in opposition, private UGS suitable to contribute to the creation of urban ecological corridors and the enhancement of regulating UES; Municipal-non-municipal. It defines UGS owned by the Municipality, that is, the primary features of urban green spatial planning and management oriented to improve green accessibility and reduce inequalities in UGS fruition by citizens.
Firstly, we classified rural and non-rural green areas. We used once again the topographic DB, gathering the class representing agricultural areas with a selection of farmlands belonging to other classes of the DB (e.g., 'green areas', 'fallow', 'non-vegetated areas') detected by visual analysis on orthoimagery and then intersecting the aggregated features with the overall UGS map derived from NDVI. The sum of the resulting features corresponds to the share of rural green areas in Padua. The remaining percentage of total green areas is therefore non-rural.
Then, we classified UGS according to their property status, in order to obtain a map of public and private green areas. For this purpose, we had to perform an overlay between cadastral maps of land parcels owned by public institutions in Padua and green surfaces derived by our mapping. This task was conducted in cooperation with Padua Municipality, which provided a table of Italian public institutions that we filtered to select the ones that were likely to own land plots in Padua (for instance, public parks and gardens are mostly the property of Padua Municipality; the University of Padua holds several studies and research facilities equipped with green spaces; Province of Padua is the owner of some secondary schools building complexes surrounded by green areas; Veneto Region owns branch offices of regional agencies provided with inner gardens; vegetated embankments of streams and rivers are usually state property; the Italian Ministry of Defence still owns several unused barracks and military complexes often equipped with vegetation). This selection was then used to query the cadastral database of Padua in order to extract publicly owned land parcels.
The following step was to overlay such parcels with our map of green areas: in this way, we obtained a map of the spatial distribution of public green areas in Padua. The remaining share of total green represents private green areas.
Starting from the public-private green classification, since our spatial data on public green areas were equipped with information about the owner institutions, with a further selection we were able to map municipality-owned green areas. This classification allows detecting the green spaces that can be directly managed by the local city planning office and the ones that require a change of ownership or an agreement to undergo changes such as the development of urban gardens or the planting of new trees.
Finally, in order to examine the relationships between the location of green areas and the distribution of population in the city, by clipping our maps, we calculated amounts and per capita values of each UGS category in the 40 urban units in which the municipality is officially divided.

Data Quality and Update Assessment
The main dataset used for green areas detection is the 2015 series of orthophotos released by AGEA/Veneto Region. Such orthophotos are updated every three years but at the time this work was performed the 2018 series was still not available. The overall image quality is very good; it was not possible to determine the exact data of the survey, but the amount of vegetation in the images suggests placing the photographic campaign in the spring-summer period.
For future developments, a good suggestion may be to repeat the workflow of the study on the subsequent series of orthophotos in order to detect changes in overall green areas and in their classification.
The dataset used to integrate NDVI results is the Topographic DataBase of Padua: it is gradually updated to cope with changes in land cover, but its official date of update is 2013 [55]. This potential lack of update was balanced out with supervised checks and corrections on surfaces that underwent recent changes in land cover.
Green spaces mapping of the four sample areas used for data validation are based on the same orthoimagery used for NDVI detection. These data were generated from a strongly supervised and detailed visual analysis, and therefore, they are highly reliable.
For the classification based on property, we operated a selection starting from an overall table of Italian public administrations and institutions provided by Padua Municipality, which can be considered thorough enough for the purposes of this study. The matching with Padua cadastral database was based upon codes designing land-owning institutions. Since Italian cadastral data suffer from flaws regarding the update of codes, the query and overlay that resulted in the map of public green areas may not have identified all the required surfaces. After performing sample checks of the location and structure of known public and private green spaces in Padua, nevertheless, we can state that our results appear to fit well with the actual property layout.

Overall Workflow
In order to facilitate the comprehension of the workflow carried out for the investigation, the main steps presented in the previous paragraphs are schematised in the flowchart in Figure 3. Data input, as well as intermediate and final output, are inserted in parallelepipeds with yellow, blue, and red borders, respectively. Processes are enclosed in rectangles. A green rectangle comprises the two result validation processes carried out to validate the final overall map of Padua UGS.

Overall Workflow
In order to facilitate the comprehension of the workflow carried out for the investigation, the main steps presented in the previous paragraphs are schematised in the flowchart in Figure 3. Data input, as well as intermediate and final output, are inserted in parallelepipeds with yellow, blue, and red borders, respectively. Processes are enclosed in rectangles. A green rectangle comprises the two result validation processes carried out to validate the final overall map of Padua UGS.

Urban Green Space Extraction (NDVI), Integration, and Cross Validation
UGS location after NDVI extraction and before integration based on the topographic DB is shown in Figure 4. The chromatic value scale adopted for representation, based on a linear interpolation of NDVI values, allows locating areas where vegetation is most flourishing, often coincident with the tree canopy. Moreover, the pattern of peri-urban agricultural fields and the state of cultivations can be distinguished in a yellow-light-green range of colours.

Urban Green Space Extraction (NDVI), Integration, and Cross Validation
UGS location after NDVI extraction and before integration based on the topo DB is shown in Figure 4. The chromatic value scale adopted for representation, b a linear interpolation of NDVI values, allows locating areas where vegetation flourishing, often coincident with the tree canopy. Moreover, the pattern of pe agricultural fields and the state of cultivations can be distinguished in a yellow-lig range of colours.  Tables 2-8. Table 2 shows the number of green areas in the four macro-sample Brentelle (macro-area 1), Forcellini (macro-area 2), Sacra Famiglia (macro-area 3), Lazzaro (macro-area 4), according to NDVI and visual analysis (see also Figure 5 parative spatial analyses between NDVI and classification by expert photo interp are very similar for each neighbourhood. Results of validation performed in the four sample areas after extracting UGS by a proper NDVI threshold and integrating them with topographic DB classes are presented in Tables 2-8. Table 2 shows the number of green areas in the four macro-sample units of Brentelle (macro-area 1), Forcellini (macro-area 2), Sacra Famiglia (macro-area 3), and San Lazzaro (macro-area 4), according to NDVI and visual analysis (see also Figure 5). Comparative spatial analyses between NDVI and classification by expert photo interpretation are very similar for each neighbourhood.   Results of the spatial linear regressions calculated on the sample areas between NDVIderived data, combined with the topographic DB, and data from visual analysis [42,43] ( Table 3), show that R values range from 0.80 and 0.84, indicating a high efficiency of the adopted methodology. Results of the error matrices calculated on the sample macro-areas (Tables 4-7) show that both user's and producer's accuracy values range from 84% to 95%, indicating a high correspondence between the two datasets. Overall accuracy and kappa coefficient for the sample macro-areas are shown in Table 8. Overall accuracy ranks above 90% in all the four areas; kappa coefficients, which measures overall agreement in classification, range from 0.80 to 0.84 with a pattern similar to R values.
In conclusion, we can state that results of comparative spatial analyses in the four macro-areas allow us to validate the UGS detection operated at the city level.

Urban Green Spaces: Estimation and Binary Classification
According to NDVI-based classification and validation, the total area of UGS is 52.23 km 2 , which represents 56% of the municipal territory of Padua. Rural green areas are 28.8 km 2 (55.14% of total UGS), while non-rural green areas are 23.43 km 2 (44.86% of total UGS).
Results about property-based UGS classification show a large preponderance of private green spaces (41.98 km 2 , 80.38% of total UGS); on the contrary, public UGS sums up to 10.25 km 2 which is less than 20% of total UGS (Table 9). Finally, Padua Municipality owns about 50% of public UGS, for an amount of about 5 km 2 (less than 10% of the total green amount). UGS not owned by the Municipality represent about 47.21 km 2 which represents about 90% of the overall mapped green spaces ( Table 9).
The visual analysis highlights a general pattern in the spatial distribution of UGS, from the city centre towards the borders: the urban core consists mostly of built-up surfaces, with a minority of sparse, small, or medium-sized green plots. Moreover, outward from the historical city, UGS become denser and broader so that their relationship with built-up areas changes ( Figure 6). One exception is represented by the industrial area of Padua, which spreads on the eastern sector of the city.  Table 11 for corresponding urban unit names). In yellow are the distances at 1 km, 3 km, 5 km, and 7 km from the city centre (yellow dot). Figure 7 (see Section 2.2) shows the distribution of UGS feature dimensions and their distance from the city centre. The scatter diagram confirms that wider green areas are mostly located in a 3-7 km distance range from the urban centre, with a peak near a 5 km distance.
Visual analysis of the map of rural and non-rural UGS (Figure 8) points out a centreborders pattern similar to the one highlighted in Figure 6, with a clear prevalence of crops in the transition belt between the city and the surrounding countryside. This map also displays a size opposition, because non-rural green appears to be located in multiple, small land parcels, while agricultural areas mainly cover wide, compact land plots.  Table 11 for corresponding urban unit names). In yellow are the distances at 1 km, 3 km, 5 km, and 7 km from the city centre (yellow dot). Figure 7 (see Section 2.2) shows the distribution of UGS feature dimensions and their distance from the city centre. The scatter diagram confirms that wider green areas are mostly located in a 3-7 km distance range from the urban centre, with a peak near a 5 km distance.
Visual analysis of the map of rural and non-rural UGS (Figure 8) points out a centreborders pattern similar to the one highlighted in Figure 6, with a clear prevalence of crops in the transition belt between the city and the surrounding countryside. This map also displays a size opposition, because non-rural green appears to be located in multiple, small land parcels, while agricultural areas mainly cover wide, compact land plots.     Classification of public and private UGS ( Figure 9) shows a predominance of privately owned greenery: it includes both a vast majority of crops and cultivated areas towards the city edges and a lot of small-sized courts, gardens, and backyards located in the old town or beside low-rise houses spreading beyond the city walls. Most of large public green areas correspond to major city public parks. Classification of public and private UGS ( Figure 9) shows a predominance of privately owned greenery: it includes both a vast majority of crops and cultivated areas towards the city edges and a lot of small-sized courts, gardens, and backyards located in the old town or beside low-rise houses spreading beyond the city walls. Most of large public green areas correspond to major city public parks. Finally, municipal UGS appear to be discontinuous, scattered subsets of municipal greenery; most recognisable areas are located outside of the core of the city ( Figure 10).
As UES are provided by all ecosystems, it is interesting to intersect the individual contribution of public (municipal-non-municipal) and private greening with rural-nonrural UGS. This matrix suggests the possibility to focus on policies aiming to enhance ecosystem services provided by private rural areas in Padua, which represent more than 60% of total private green (Table 10). Finally, municipal UGS appear to be discontinuous, scattered subsets of municipal greenery; most recognisable areas are located outside of the core of the city ( Figure 10).
As UES are provided by all ecosystems, it is interesting to intersect the individual contribution of public (municipal-non-municipal) and private greening with rural-nonrural UGS. This matrix suggests the possibility to focus on policies aiming to enhance ecosystem services provided by private rural areas in Padua, which represent more than 60% of total private green (Table 10).

Urban Green Spaces and Population in Padua Urban Units
In the current literature, UGS per capita is widely used as an indicator to assess environmental quality and availability of green spaces by citizens [56][57][58][59]. On average, European citizens have access to about 18 m 2 of public green space within the boundary of their city, with a benchmarking of 20 m 2 per person [12]. In Italy, standards fixed at a national level in 1968 and still in force set a minimum of 9 m 2 of accessible public green areas per capita [60].
Several local case studies have been performed on this topic [57,61,62]; however, comparative analyses of results might be misleading due to the differences in the definition of green spaces, extraction methodology, and classification criteria. Moreover, the results of municipality-level research are widely influenced by the extent of administrative boundaries and land use patterns of study areas.

Urban Green Spaces and Population in Padua Urban Units
In the current literature, UGS per capita is widely used as an indicator to assess environmental quality and availability of green spaces by citizens [56][57][58][59]. On average, European citizens have access to about 18 m 2 of public green space within the boundary of their city, with a benchmarking of 20 m 2 per person [12]. In Italy, standards fixed at a national level in 1968 and still in force set a minimum of 9 m 2 of accessible public green areas per capita [60].
Several local case studies have been performed on this topic [57,61,62]; however, comparative analyses of results might be misleading due to the differences in the defini-tion of green spaces, extraction methodology, and classification criteria. Moreover, the results of municipality-level research are widely influenced by the extent of administrative boundaries and land use patterns of study areas.
Therefore, in our study, we focused on a comparative assessment of UGS categories in Padua urban units in order to highlight variability and imbalances in their spatial distribution.
As illustrated in Section 2.5, we calculated extents and per capita rates of the different types of classified UGS for each one of the urban units in Padua (Table 11). Population in urban units is represented in Figure 11. Therefore, in our study, we focused on a comparative assessment of UGS categories in Padua urban units in order to highlight variability and imbalances in their spatial distribution.
As illustrated in Section 2.5, we calculated extents and per capita rates of the different types of classified UGS for each one of the urban units in Padua (Table 11). Population in urban units is represented in Figure 11.  Table  11 for the corresponding urban units).
With respect to Table 11, we focus on overall, public, and municipal UGS, since these categories are, respectively, relevant for their contribution to ES provisioning, for their potential public accessibility, and for their actual public fruition managed by Padua Municipality. Values for the aforementioned UGS are graphically represented in Figure 12; density maps concerning public and municipal UGS are displayed in Figure 13.  Table 11 for the corresponding urban units).
With respect to Table 11, we focus on overall, public, and municipal UGS, since these categories are, respectively, relevant for their contribution to ES provisioning, for their potential public accessibility, and for their actual public fruition managed by Padua Municipality. Values for the aforementioned UGS are graphically represented in Figure 12; density maps concerning public and municipal UGS are displayed in Figure 13. Analysis of values in urban units shows that those included in the core city or in the densely populated northern district show the lowest rates of total UGS (e.g., 13.88 m 2 per capita in unit 1.1), due to the remarkable share of built-up spaces. The highest rates belong to scarcely populated, mostly agricultural peripheral urban units (e.g., 1503.08 m 2 per capita in unit 12); unit 30.1 (4096.68 m 2 per capita) is a special case because it corresponds to the industrial area of the city with a very low population density. A similar pattern applies to public UGS, for example, urban parks, where, nevertheless, there is a difference between unit 1.1 (where the city Hall, the city cathedral, and some of the main squares belong), which shows the lowest rate with 2.44 m 2 per capita, and the remaining urban units in the historical centre, showing slightly higher rates (up to 31.76 m 2 per capita in unit 1.3); other critical urban units are the dense residential districts located immediately outside the city wall system (e.g., 6.24 m 2 per capita in unit 3, 8.71 m 2 per capita in unit 25.1), implying that, in some cases, the city development outside its ancient core did not adequately consider the need for public green spaces.  As regards municipal UGS areas, their distribution in the compact city is comparable to public UGS, with the lowest rate in unit 1.1 (1.11 m 2 per capita) and remarkably low values for units located north and west of the historical centre and adjacent to the wall system (e.g., 3.83 m 2 per capita in unit 3). A special case is unit 15 on the western municipality border (6.49 m 2 per capita), where a compact residential suburb developed in a formerly rural area crossed by two water streams: the area is therefore quite valuable from a landscape and environmental perspective but currently lacks municipal green facilities.
Then, we briefly focus on the five most densely populated urban units ( Padua). These values show how issues related to UGS availability are particularly relevant in the compact city and require an overall vision to be tackled.
Finally, we gathered together results on different UGS categories in a general scheme (Figure 14), where values of every UGS category assessed are divided into five equally numerous subsets or quintiles.
The scheme in Figure 14 enables an aggregate assessment for every urban unit and may help to direct policies and strategies at the municipality level, as discussed in Section 3.4. Predictably, urban units with higher population density prove less adequate in UGS per capita equipment. Moreover, a comparative analysis of density variations for different categories in critical urban units may serve as an indicator of strengths and weaknesses to start building green areas management and development strategies, thus reinforcing the link between quantitative assessments and a planning perspective.

Urban Green Spaces: From Mapping and Classifying to Spatial Planning
The study presented in this paper and, in particular, the results of green areas classifications, endows urban planning with useful geospatial data and applicative tools for UGS management and implementation.
Performance-based urban planning focuses on the environmental performance of natural and artificial features in planned areas. Although studies and publications concerning the inclusion of UES in urban plans are more and more increasing [63,64], to this day, they still find limited applications in planning activities by local administrators [23,65,66].
In our target area, the current urban development plan is the Territorial Management Plan of Padua [46], which outlines a thematic spatial framework for the implementation of a green-and-blue infrastructure leaning upon recognisable urban features such as the wall system, the network of streams, and rivers, the major public parks and the peri-urban rural areas. However, it does not embrace UES as an analytical or planning approach. Accurate data and thematic maps about the distribution and property status of UGS, together with an assessment of related UES, could help planners in designing a site-specific, citizen-oriented urban GI.
More generally, issues concerning inequalities in access and fruition of UGS caused by built-up density in some districts would require multiple solutions, ranging from the increase of vegetation and equipped areas in municipal UGS to a focus on designing public green areas within reuse or regeneration projects for neglected sites. As an example from our case study, from Figure 14, it can be inferred that more densely populated districts show a general lack of green spaces; among them, we chose four urban units (marked with letters a-d in Figure 14) as an example of different possible general UGS management approaches: Unit 24 (a) is a very densely populated district with low levels of per capita UGS; however, it relies on a certain amount of municipal green spaces. UGS management may lean on improving accessibility, connectivity, and vegetation equipment of municipal green areas;  Table 11 and subdivided into quintiles. Units a-d are quoted as samples for outlining UGS management strategies in Section 3.4.
Unit 9 (b) is also densely populated, has a higher overall UGS equipment than unit 24 but inadequate public and municipal green areas; at the same time, it is endowed with some private green areas. UGS management will be primarily based on maintaining adequate vegetation amounts in private green areas and fostering public accessibility, especially by walking on cycling, among these areas, also considering the proximity to a water body; Unit 1.3 (c) is a quite densely populated district where municipal UGS equipment stands out on the other UGS categories, due to the presence of city parks. It may serve as a green hub for the city centre and as a soft mobility gateway to the stream network, in the framework of implementing an urban GI; Unit 15 (d) is an outer district with medium population density, and good overall public UGS equipment but lack municipal UGS (see also Section 3.3): it would be strategic to develop green facilities by agreements among local administrators and public institutions which own the green areas, enhancing their accessibility and soft mobility connections.
A further action, to be undertaken especially in core cities and surrounding districts, is improving connections between built-up areas and urban GI or peri-urban open spaces. In Padua, for example, a potentially useful planning tool may be the project of the new city cycling network [67], which could be integrated with the design of green-and-blue infrastructure, thus encouraging and easing soft mobility for leisure and health.
Furthermore, the location of non-municipal public green areas may help municipalities in developing urban GI strategies in collaboration with other public institutions. Additionally, we noted how domestic gardens and other private UGS are able to provide a fair amount of regulating and supporting UES [68] and should be considered, together with public and accessible UGS, as components of urban GI networks. For these reasons, spatial planners should include private UGS in ecosystem services mapping in support of urban design scenarios. It would also be useful for municipalities to involve private citizens that own significant land plots in some form of sharing of green spaces with citizens. Finally, to increase continuity of public, accessible UGS, trade-ins between private and municipal green areas may be promoted by local authorities, within multilevel governance of urban land use and land cover.
To sum up, it can be stated that spatial distribution of categorised green spaces suggests some possible guidelines for managing these spaces, partially derived from our case study but provided with general validity.
In Table 12, we present a simple layout of planning strategies related to the UGS categories discussed in our case study and to the main UES they can provide.

Private
-Undertake initiatives to prevent the abandonment of crop fields and allow accessibility of their borders -Incentive peri-urban social farming, agroecology, and short chains to supply city consumers [69] -Protect and incentive private vegetation development [27] -Consider the location and distribution of relevant private green spaces for environmentally centred plans and initiatives [31] Public, non-municipal -Foster connectivity and accessibility by pedestrian and cycle paths connected to urban GI -Encourage teaching and educational activities involving children and students [70] -Promote public accessibility through agreements between Municipalities and other public institutions -Create and enhance urban GI leaning on existing city network systems Municipal -Promote their use as allotments or community gardens [71] -Consider the reuse of uncultivated fields as semi-natural peri-urban buffer zones connected to urban GI -Increase accessibility to public uses with pocket parks or small public urban green spaces [72] -Promote trade-ins of land plots with other public or private actors to enhance green spaces continuity

Conclusions
Urban green space management is a crucial topic in current urban planning, and it will presumably gain further importance in the near future, especially by considering relationships with climate change and other urgent environmental issues. Classification of green spaces by type and property can help planners and authorities to develop site-specific policies and initiatives in order to enhance the environmental and social benefits provided by vegetated and pervious areas to the urban population.
In this framework, we analysed the case study of Padua, a middle-sized historical city in the northeast of Italy. We mapped urban green spaces in Padua municipality by adopting the NDVI vegetation index; then, we classified them into three binary categories: ruralnon-rural, public-private, municipal-non-municipal. Finally, we calculated statistics on the distribution of green areas by urban units and on their relationship with the city population. Results indicate that rural green areas are predominant over non-rural, especially in the peripheral neighbourhoods, by a prevalence of private green areas over public. Finally, municipally owned green spaces represent less than 10% of total UGS. Then, we focused on imbalances in green space distribution within the municipality, finding that most of the densely populated urban units display values for public and municipal greenery that are definitely below the city average.
Using our results, we attempted to link quantitative assessments with a planning perspective and outline for local spatial planners and administrators aiming to foster environmental benefits for each green space category and reduce spatial inequalities, in a performance-based planning framework which considers the contribution of private green spaces for UES provision, also considering that in the city of Padua private UGS are by far predominant over public ones.
Our work could serve as a basis for further analysis on urban green connectivity and accessibility for citizens, involving an assessment of green spaces by their reachable range by walking or bike; other possible analyses may include mapping selected ecosystem services at the municipality level, with focuses on strategic districts. Future uses and developments of this study would require periodic updates of results and thematic maps to cope with changes in land cover, land use, or property.