Next Article in Journal
Urban Overheating Assessment through Prediction of Surface Temperatures: A Case Study of Karachi, Pakistan
Next Article in Special Issue
Spatially Varying Effects of Street Greenery on Walking Time of Older Adults
Previous Article in Journal
Fast Geo-Location Method Based on Panoramic Skyline in Hilly Area
Previous Article in Special Issue
Hydrological Modeling of Green Infrastructure to Quantify Its Effect on Flood Mitigation and Water Availability in the High School Watershed in Tucson, AZ
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Guglielmo Pristeri
1,
Francesca Peroni
2,*,
Salvatore Eugenio Pappalardo
3,
Daniele Codato
3,
Antonio Masi
4 and
Massimo De Marchi
3
1
GIScience and Unmanned Systems for the Integrated Management of the Territory and the Natural Resources, University of Padua, 35100 Padua, Italy
2
Department of Historical and Geographic Sciences and the Ancient World (DiSSGeA), University of Padua, 35100 Padua, Italy
3
Department of Civil, Environmental and Architectural Engineering (ICEA), University of Padua, 35100 Padua, Italy
4
Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padua, 35100 Padua, Italy
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2021, 10(8), 538; https://doi.org/10.3390/ijgi10080538
Submission received: 31 May 2021 / Revised: 29 July 2021 / Accepted: 5 August 2021 / Published: 10 August 2021

Abstract

:
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.

1. Introduction

1.1. 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 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.

1.2. 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 km2 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].
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.

1.3. 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 high-resolution 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].

2. Data and Methods

2.1. 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.

2.2. 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:
N D V I = N I R R E D N I R + R E D
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.

2.3. 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 km2), Forcellini (macro-area 2, 2.66 km2), Sacra Famiglia (macro-area 3, 2.77 km2), and San Lazzaro (macro-area 4, 3.44 km2). In their diversity, they represent a good sample of the relationships between green spaces and urban fabric in Padua.

2.4. 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.

2.5. 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.

2.6. 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.

3. Results and Discussion

3.1. 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.
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 Table 2, Table 3, Table 4, Table 5, Table 6, Table 7 and Table 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.
Table 2. Green surfaces from NDVI and from visual analysis in the four sample areas.
Table 2. Green surfaces from NDVI and from visual analysis in the four sample areas.
Macro-Sample AreaGreen Areas from NDVI/Topo DB (km2)Green Areas from Visual Analysis (km2)
Brentelle (1)1.841.80
Forcellini (2)1.591.64
Sacra Famiglia (3)1.731.71
San Lazzaro (4)1.231.24
Results of the spatial linear regressions calculated on the sample areas between NDVI-derived 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.
Table 3. Results of spatial linear regressions in the four sample areas.
Table 3. Results of spatial linear regressions in the four sample areas.
Sample AreaR Value
Brentelle0.80
Forcellini0.81
Sacra Famiglia0.84
San Lazzaro0.84
Results of the error matrices calculated on the sample macro-areas (Table 4, Table 5, Table 6 and Table 7) show that both user’s and producer’s accuracy values range from 84% to 95%, indicating a high correspondence between the two datasets.
Table 4. Error matrix for Brentelle. 0 = not green; 1 = green. Rows = visual analysis; columns = NDVI + Topographic DB. UA = user’s accuracy; PA = producer’s accuracy.
Table 4. Error matrix for Brentelle. 0 = not green; 1 = green. Rows = visual analysis; columns = NDVI + Topographic DB. UA = user’s accuracy; PA = producer’s accuracy.
01TOTALPA (%)
02,785,201351,4553,136,65684.15
1524,7806,870,4277,395,20795.13
TOTAL3,309,9817,221,88210,531,863
UA (%)88.8092.90
Table 5. Error matrix for Forcellini. 0 = not green; 1 = green. Rows = visual analysis; columns = NDVI + Topographic DB. UA = user’s accuracy; PA = producer’s accuracy.
Table 5. Error matrix for Forcellini. 0 = not green; 1 = green. Rows = visual analysis; columns = NDVI + Topographic DB. UA = user’s accuracy; PA = producer’s accuracy.
01TOTALPA (%)
04,899,198646,4425,545,64091.01
1484,2675,864,3636,348,63090.07
TOTAL5,383,4656,510,80511,894,270
UA (%)88.3492.37
Table 6. Error matrix for Sacra Famiglia. 0 = not green; 1 = green. Rows = visual analysis; columns = NDVI + Topographic DB. UA = user’s accuracy; PA = producer’s accuracy.
Table 6. Error matrix for Sacra Famiglia. 0 = not green; 1 = green. Rows = visual analysis; columns = NDVI + Topographic DB. UA = user’s accuracy; PA = producer’s accuracy.
01TOTALPA (%)
04,382,170390,4914,772,66189.31
1524,7176,323,1816,847,89894.18
TOTAL4,906,8876,713,67211,620,559
UA (%)91.8292.34
Table 7. Error matrix for San Lazzaro. 0 = not green; 1 = green. Rows = visual analysis; columns = NDVI + Topographic DB. UA = user’s accuracy; PA = producer’s accuracy.
Table 7. Error matrix for San Lazzaro. 0 = not green; 1 = green. Rows = visual analysis; columns = NDVI + Topographic DB. UA = user’s accuracy; PA = producer’s accuracy.
01TOTALPA (%)
08,320,253403,7488,724,00193.54
1574,4174,359,7104,934,12791.52
TOTAL8,894,6704,763,45813,658,128
UA (%)95.3788.36
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.
Table 8. Overall accuracy and kappa coefficient in the four sample areas.
Table 8. Overall accuracy and kappa coefficient in the four sample areas.
Sample AreaOverall Accuracy (%)Kappa Coefficient
Brentelle91.680.80
Forcellini90.490.81
Sacra Famiglia92.120.84
San Lazzaro92.840.84

3.2. Urban Green Spaces: Estimation and Binary Classification

According to NDVI-based classification and validation, the total area of UGS is 52.23 km2, which represents 56% of the municipal territory of Padua. Rural green areas are 28.8 km2 (55.14% of total UGS), while non-rural green areas are 23.43 km2 (44.86% of total UGS).
Results about property-based UGS classification show a large preponderance of private green spaces (41.98 km2, 80.38% of total UGS); on the contrary, public UGS sums up to 10.25 km2 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 km2 (less than 10% of the total green amount). UGS not owned by the Municipality represent about 47.21 km2 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.
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 centre–borders 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.
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–non-rural 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).

3.3. 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 m2 of public green space within the boundary of their city, with a benchmarking of 20 m2 per person [12]. In Italy, standards fixed at a national level in 1968 and still in force set a minimum of 9 m2 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.
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.
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.
Regarding municipality-scale UGS, in Padua, there is quite a remarkable total of 247.35 m2 of per capita UGS; among them, public UGS areas amount to 48.54 m2 per capita, and municipal UGS areas are 23.77 m2 per capita.
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 m2 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 m2 per capita in unit 12); unit 30.1 (4096.68 m2 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 m2 per capita, and the remaining urban units in the historical centre, showing slightly higher rates (up to 31.76 m2 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 m2 per capita in unit 3, 8.71 m2 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 m2 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 m2 per capita in unit 3). A special case is unit 15 on the western municipality border (6.49 m2 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 (1.1; 9; 24; 25.1; 25.2), hosting almost 23% of the city population. In these units, UGS values rank far below municipality average values: total UGS areas range from 13.88 to 68.36 m2 per capita (247.35 m2 per capita in Padua); public UGS areas range from 2.44 to 16.75 m2 per capita (48.54 m2 per capita in Padua); municipal UGS areas range from 1.11 to 15.11 m2 per capita (23.77 m2 per capita in 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.

3.4. 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;
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.

4. 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: rural–non-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.

Author Contributions

Conceptualisation, Guglielmo Pristeri, Francesca Peroni, and Salvatore Eugenio Pappalardo; Methodology, Guglielmo Pristeri, Daniele Codato, and Salvatore Eugenio Pappalardo; Analysis, Guglielmo Pristeri and Francesca Peroni; Writing—original draft preparation, Guglielmo Pristeri, Francesca Peroni, and Salvatore Eugenio Pappalardo; Writing—review and editing, Salvatore Eugenio Pappalardo, Antonio Masi, and Massimo De Marchi; Validation, Daniele Codato; Visualisation, Francesca Peroni and Guglielmo Pristeri; Project administration, Massimo De Marchi; Funding acquisition, Massimo De Marchi; Supervision, Massimo De Marchi and Antonio Masi. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Master in GIScience and Unmanned Systems for the integrated management of the territory and the natural resources (ICEA Department, University of Padua (Italy) and the research scholarship ‘First stage of the implementation of a Geographic Information System of non-municipal urban green spaces in Padua municipality territory’, funded by ICEA Department, University of Padua (Italy) and Department of Environment and Territory, Municipality of Padua.

Data Availability Statement

No ethical approval was required due to the use of publically available and/or aggregate secondary data.

Acknowledgments

The authors would like to thank members of the Department of Environment and Territory and the Department of Informatics Services, Municipality of Padua, especially Fiorenza Colombo and Daniela Luise.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Borgström, S.; Similä, J. Integration of ecosystem services into decision-making. In Towards a Sustainable and Genuinely Green Economy. The Value and Social Significance of Ecosystem Services in Finland (TEEB for Finland) Synthesis and Roadmap; Jäppinen, J.-P., Heliölä, J., Eds.; Finland Ministry of the Environment: Helsinki, Finland, 2015; pp. 73–80. [Google Scholar]
  2. Colding, J. The role of ecosystem services in contemporary urban planning. In Urban Ecology: Patterns, Processes and Applications; Niemelä, J., Breuste, J.H., Elmqvist, T., Guntenspergen, G., James, P., McIntyre, N.E., Eds.; Oxford University Press: Oxford, UK, 2011; pp. 228–237. [Google Scholar]
  3. Cortinovis, C.; Geneletti, D. A framework to explore the effects of urban planning decisions on regulating ecosystem services in cities. Ecosyst. Serv. 2019, 38, 100946. [Google Scholar] [CrossRef]
  4. Davies, C.; Lafortezza, R. Transitional path to the adoption of nature-based solutions. Land Use Policy 2019, 80, 406–409. [Google Scholar] [CrossRef]
  5. Faivre, N.; Fritz, M.; Freitas, T.; de Boissezon, B.; Vandewoestijne, S. Nature-based solutions in the EU: Innovating with nature to address social, economic and environmental challenges. Environ. Res. 2017, 159, 509–518. [Google Scholar] [CrossRef]
  6. Conceptual Framework Working Group of the Millennium Ecosystem Assessment. Ecosystems and Human Well-Being. A Framework for Assessment; Island Press: Washington, DC, USA, 2003; pp. 7–19. [Google Scholar]
  7. Haines-Young, R.; Potschin, M. Common International Classification of Ecosystem Services (CICES): Consultation on Version 4, August–December 2012; European Environment Agency: Copenhagen, Denmark, 2013.
  8. Torres, A.V.; Tiwari, C.; Atkinson, S.F. Progress in ecosystem services research: A guide for scholars and practitioners. Ecosyst. Serv. 2021, 49, 101267. [Google Scholar] [CrossRef]
  9. Karimi, A.; Yazdandad, H.; Fagerholm, N. Evaluating social perceptions of ecosystem services, biodiversity and land management: Trade-offs, synergies and implications for landscape planning and management. Ecosyst. Serv. 2020, 45, 101188. [Google Scholar] [CrossRef]
  10. Hatan, S.; Fleischer, A.; Tchetchik, A. Economic valuation of cultural ecosystem services: The case of landscape aesthetics in the agritourism market. Ecol. Econ. 2021, 184, 107005. [Google Scholar] [CrossRef]
  11. European Commission. Mapping and Assessment of Ecosystems and Their Services. Urban Ecosystems, 4th Report–Final May 2016; Office for Official Publications of the European Communities: Luxembourg, 2016. [Google Scholar]
  12. Maes, J.; Zulian, G.; Günther, S.; Thijssen, M.; Raynal, J. Enhancing Resilience of Urban Ecosystems through Green Infrastructure (EnRoute) Final Report; Publications Office of the European Union: Luxembourg, 2019; ISBN 978-92-79-98984-1. [CrossRef]
  13. Andersson-Sköld, Y.; Klingberg, J.; Gunnarsson, B.; Cullinane, K.; Gustafsson, I.; Hedblom, M.; Knez, I.; Lindberg, F.; Sang, Å.O.; Pleijel, H.; et al. A framework for assessing urban greenery’s effects and valuing its ecosystem services. J. Environ. Manag. 2018, 205, 274–285. [Google Scholar] [CrossRef]
  14. Bush, J.; Doyon, A. Building urban resilience with nature-based solutions: How can urban planning contribute? Cities 2019, 95, 102483. [Google Scholar] [CrossRef]
  15. Xie, L.; Bulkeley, H. Nature-based solutions for urban biodiversity governance. Environ. Sci. Policy 2020, 110, 77–87. [Google Scholar] [CrossRef]
  16. Giachino, C.; Pattanaro, G.; Bertoldi, B.; Bollani, L.; Bonadonna, A. Nature-based solutions and their potential to attract the young generations. Land Use Policy 2021, 101, 105176. [Google Scholar] [CrossRef]
  17. Seiwert, A.; Rößler, S. Understanding the term green infrastructure: Origins, rationales, semantic content and purposes as well as its relevance for application in spatial planning. Land Use Policy 2020, 97, 104785. [Google Scholar] [CrossRef]
  18. Li, J.; Wang, Y.; Ni, Z.; Chen, S.; Xia, B. An integrated strategy to improve the microclimate regulation of green-blue-grey infrastructures in specific urban forms. J. Clean. Prod. 2020, 271, 122555. [Google Scholar] [CrossRef]
  19. Almenar, B.J.; Elliot, T.; Rugani, B.; Philippe, B.; Gutierrez, N.T.; Sonnemann, G.; Geneletti, D. Nexus between nature-based solutions, ecosystem services and urban challenges. Land Use Policy 2021, 100, 104898. [Google Scholar] [CrossRef]
  20. Kendig, L.; Connor, S. Performance Zoning; APA Planners Press: Chicago, IL, USA, 1980. [Google Scholar]
  21. Frew, T.G. The Implementation of Performance-Based Planning in Queensland under the Integrated Planning Act 1997: An Evaluation of Perceptions and Planning Schemes; Queensland University of Technology, Faculty of Built Environment and Engineering: Brisbane, Australia, 2011. [Google Scholar]
  22. Nin, M.; Soutullo, A.; Rodríguez-Gallego, L.; Minin, D.E. Ecosystem services-based land planning for environmental impact avoidance. Ecosyst. Serv. 2016, 17, 172–184. [Google Scholar] [CrossRef]
  23. Ronchi, S.; Arcidiacono, A.; Pogliani, L. Integrating green infrastructure into spatial planning regulations to improve the performance of urban ecosystems. Insights from an Italian case study. Sustain. Cities Soc. 2020, 53, 101907. [Google Scholar] [CrossRef]
  24. van Oijstaeijen, W.; van Passel, S.; Cools, J. Urban green infrastructure: A review on valuation toolkits from an urban planning perspective. J. Environ. Manag. 2020, 267, 110603. [Google Scholar] [CrossRef] [PubMed]
  25. World Health Organization. Urban Green Spaces: A Brief for Action; WHO Regional Office for Europe: Copenhagen, Denmark, 2017. [Google Scholar]
  26. Cameron, R.W.F.; Blanuša, T.; Taylor, J.E.; Salisbury, A.; Halstead, A.J.; Henricot, B.; Thompson, K. The domestic garden-its contribution to urban green infrastructure. Urban For. Urban Green. 2012, 11, 129–137. [Google Scholar] [CrossRef]
  27. Clark, C.; Ordóñez, C.; Livesley, S.J. Private tree removal, public loss: Valuing and enforcing existing tree protection mechanisms is the key to retaining urban trees on private land. Landsc. Urban Plan. 2020, 203, 103899. [Google Scholar] [CrossRef]
  28. Cervinka, R.; Schwab, M.; Schönbauer, R.; Hämmerle, I.; Pirgie, L.; Sudkamp, J. My garden-my mate? Perceived restorativeness of private gardens and its predictors. Urban For. Urban Green. 2016, 16, 182–187. [Google Scholar] [CrossRef]
  29. Glavan, M.; Schmutz, U.; Williams, S.; Corsi, S.; Monaco, F.; Kneafsey, M.; Rodriguez, G.P.A.; Čenič-Istenič, M.; Pintar, M. The economic performance of urban gardening in three European cities–examples from Ljubljana, Milan and London. Urban For. Urban Green. 2018, 36, 100–122. [Google Scholar] [CrossRef]
  30. Mimet, A.; Kerbiriou, C.; Simon, L.; Julien, J.F.; Raymond, R. Contribution of private gardens to habitat availability, connectivity and conservation of the common pipistrelle in Paris. Landsc. Urban Plan. 2020, 193, 103671. [Google Scholar] [CrossRef]
  31. Liu, O.Y.; Russo, A. Assessing the contribution of urban green spaces in green infrastructure strategy planning for urban ecosystem conditions and services. Sustain. Cities Soc. 2021, 68, 102772. [Google Scholar] [CrossRef]
  32. Tahvonen, O.; Airaksinen, M. Low-density housing in sustainable urban planning–scaling down to private gardens by using the green infrastructure concept. Land Use Policy 2018, 75, 478–485. [Google Scholar] [CrossRef]
  33. Pearsall, H.; Eller, J.K. Locating the green space paradox: A study of gentrification and public green space accessibility in Philadelphia, Pennsylvania. Landsc. Urban Plan. 2020, 195, 103708. [Google Scholar] [CrossRef]
  34. Contesse, M.; van Vliet, B.J.M.; Lenhart, J. Is urban agriculture urban green space? A comparison of policy arrangements for urban green space and urban agriculture in Santiago de Chile. Land Use Policy 2018, 71, 566–577. [Google Scholar] [CrossRef]
  35. Feltynowski, M.; Kronenberg, J. Urban green spaces—An underestimated resource in third-tier towns in Poland. Land 2020, 9, 453. [Google Scholar] [CrossRef]
  36. Caputo, P.; Zagarella, F.; Cusenza, M.A.; Mistretta, M.; Cellura, M. Energy-environmental assessment of the UIA-OpenAgri case study as urban regeneration project through agriculture. Sci. Total Environ. 2020, 729, 138819. [Google Scholar] [CrossRef]
  37. Walls of Padua. A guide to the Renaissance Fortification System (in Italian: Mura di Padova. Guida al Sistema Bastionato Rinascimentale; Fadini, U. (Ed.) Edibus: Vicenza, Italy, 2013. [Google Scholar]
  38. Padua and Metropolitan City. 1807-2007 Metamorphosis of the Urban Landscape (in Italian: Padova e la Città Metropolitana. 1807–2007 la Metamorfosi del Paesaggio Urbano; Marzari, S. (Ed.) I Antichi Editori Venezia: Venezia, Italy, 2008. [Google Scholar]
  39. The Image of Metropolitan Territory. The Metropolitan City of Padua (in Italian: L’immagine del Territorio Metropolitano. La Città Metropolitana di Padova; Boschetto, P.; Schiavon, A. (Eds.) Cleup: Padova, Italy, 2011. [Google Scholar]
  40. Padova, C.D. 2019 Statistical Yearbook (in Italian: Annuario Statistico 2019); Comune di Padova: Padua, Italy, 2020. [Google Scholar]
  41. Munafò, M. (Ed.) Land Take, Territorial Dynamics and Ecosystem services (in Italian: Consumo di Suolo, Dinamiche Territoriali e Servizi Ecosistemici. Edizione 2020; Sistema Nazionale Protezione Ambiente: Rome, Italy, 2020. [Google Scholar]
  42. Pristeri, G.; Peroni, F.; Pappalardo, S.E.; Codato, D.; Castaldo, A.G.; Masi, A.; de Marchi, M. Mapping and assessing soil sealing in Padua municipality through biotope area factor index. Sustainability 2020, 12, 5167. [Google Scholar] [CrossRef]
  43. Peroni, F.; Pristeri, G.; Codato, D.; Pappalardo, S.E.; de Marchi, M. Biotope area factor: An ecological urban index to geovisualize soil sealing in Padua, Italy. Sustainability 2020, 12, 150. [Google Scholar] [CrossRef] [Green Version]
  44. European Commission. Roadmap to a Resource Efficient Europe; COM(2011) 571; European Commission: Brussels, Belgium, 2011. [Google Scholar]
  45. Veneto, R.D. Regional Law April 23, 2004, n. 11 (in Italian: Legge Regionale 23 Aprile 2004, n. 11; Regione del Veneto: Venice, Italy, 2004. [Google Scholar]
  46. Padova, C.D. Territorial Management Plan. Technical Implementation Standards (in Italian: Piano di Assetto del Territorio. Norme Tecniche di Attuazione; Comune di Padova: Padua, Italy, 2014. [Google Scholar]
  47. Padova, C.D. Plan of Interventions. Technical Implementation Standards (in Italian: Piano Degli Interventi. Norme Tecniche di Attuazione; Comune di Padova: Padua, Italy, 2018. [Google Scholar]
  48. Veneto, R.D. Regional Law June 6, 2017, n. 14 (in Italian: Legge Regionale 6 Giugno 2017, n. 14; Regione del Veneto: Venice, Italy, 2017. [Google Scholar]
  49. Veneto, R.D. Regional Law April 4, 2019, n. 14 (in Italian: Legge Regionale 4 Aprile 2019, n. 14; Regione del Veneto: Venice, Italy, 2017. [Google Scholar]
  50. Veneto, R.D. Regional Government Decree April 29, 2019, n. 64 (in Italian: Decreto Della Giunta Regionale 29 Aprile 2019, n. 64; Regione del Veneto: Venice, Italy, 2019. [Google Scholar]
  51. Xue, J.; Su, B. Significant remote sensing vegetation indices: A review of developments and applications. J. Sens. 2017, 2017, 1353691. [Google Scholar] [CrossRef] [Green Version]
  52. Rondeaux, G.; Steven, M.; Baret, F. Optimization of soil-adjusted vegetation indices. Remote Sens. Environ. 1996, 55, 95–107. [Google Scholar] [CrossRef]
  53. Spadoni, G.L.; Cavalli, A.; Congedo, L.; Munafò, M. Analysis of normalized difference vegetation index (NDVI) multi-temporal series for the production of forest cartography. Remote Sens. Appl. Soc. Environ. 2020, 20, 100419. [Google Scholar] [CrossRef]
  54. Atasoy, M. Monitoring the urban green spaces and landscape fragmentation using remote sensing: A case study in Osmaniye, Turkey. Environ. Monit. Assess. 2018, 190, 713. [Google Scholar] [CrossRef]
  55. Metadata for Topographic DataBase of Padua Municipality. Available online: https://geodati.gov.it/geoportale/ (accessed on 19 May 2021).
  56. Wüstemann, H.; Kalisch, D. Towards a National Indicator for Urban Green Space Provision and Environmental Inequalities in Germany: Method and Findings. SFB 649 Discussion Paper, No. 2016-022; Humboldt University of Berlin, Collaborative Research Center 649-Economic Risk: Berlin, Germany, 2016. [Google Scholar]
  57. Kabisch, N.; Strohbach, M.; Haase, D.; Kronenberg, J. Urban green space availability in European cities. Ecol. Indic. 2016, 70, 586–596. [Google Scholar] [CrossRef]
  58. Russo, A.; Cirella, G.T. Modern compact cities: How much greenery do we need? Int. J. Environ. Res. Public Health 2018, 15, 2180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Pafi, M.; Siragusa, A.; Ferri, S.; Halkia, S. Measuring the Accessibility of Urban Green Areas: A Comparison of the Green ESM with Other Datasets in Four European Cities; Publications Office of the European Union: Luxembourg, 2016; ISBN 9789279612862.
  60. Repubblica Italiana. Interministerial Decree, April 2, 1968, n. 1444 (in Italian: Decreto Interministeriale 2 Aprile 1968, n. 1444; Repubblica Italiana: Roma, Italy, 1968. [Google Scholar]
  61. Badiu, D.L.; IojǍ, C.I.; PǍtroescu, M.; Breuste, J.; Artmann, M.; NiţǍ, M.R.; GrǍdinaru, S.R.; Hossu, C.A.; Onose, D.A. Is urban green space per capita a valuable target to achieve cities’ sustainability goals? Romania as a case study. Ecol. Indic. 2016, 70, 53–66. [Google Scholar] [CrossRef]
  62. Shi, L.; Halik, Ü.; Abliz, A.; Mamat, Z.; Welp, M. Urban green space accessibility and distribution equity in an arid oasis city: Urumqi, China. Forests 2020, 11, 690. [Google Scholar] [CrossRef]
  63. Gómez-Baggethun, E.; Barton, D.N. Classifying and valuing ecosystem services for urban planning. Ecol. Econ. 2013, 86, 235–245. [Google Scholar] [CrossRef]
  64. Cortinovis, C.; Geneletti, D. Mapping and assessing ecosystem services to support urban planning: A case study on brownfield regeneration in Trento, Italy. One Ecosyst. 2018, 3, e25477. [Google Scholar] [CrossRef]
  65. Frew, T.; Baker, D.; Donehue, P. Performance based planning in Queensland: A case of unintended plan-making outcomes. Land Use Policy 2016, 50, 239–251. [Google Scholar] [CrossRef]
  66. Arcidiacono, A.; Ronchi, S.; Salata, S. Managing multiple ecosystem services for landscape conservation: A green infrastructure in Lombardy region. Procedia Eng. 2016, 161, 2297–2303. [Google Scholar] [CrossRef] [Green Version]
  67. Padova, C.D. Padua Bici Masterplan 2018-2022. Relation (in Italian: Bici Masterplan di Padova 2018–2022. Relazione; Comune di Padova: Padua, Italy, 2019. [Google Scholar]
  68. Dewaelheyns, V.; Rogge, E.; Gulinck, H. Putting domestic gardens on the agenda using empirical spatial data: The case of Flanders. Appl. Geogr. 2014, 50, 132–143. [Google Scholar] [CrossRef]
  69. Diehl, J.A.; Sweeney, E.; Wong, B.; Sia, C.S.; Yao, H.; Prabhudesai, M. Feeding cities: Singapore’s approach to land use planning for urban agriculture. Glob. Food Secur. 2020, 26, 100377. [Google Scholar] [CrossRef]
  70. Gómez-Villarino, M.T.; Ruiz-Garcia, L. Adaptive design model for the integration of urban agriculture in the sustainable development of cities. A case study in northern Spain. Sustain. Cities Soc. 2021, 65, 102595. [Google Scholar] [CrossRef]
  71. Cabral, I.; Keim, J.; Engelmann, R.; Kraemer, R.; Siebert, J. Urban forestry & urban greening ecosystem services of allotment and community gardens: A Leipzig, Germany case study. Urban For. Urban Green. 2017, 23, 44–53. [Google Scholar] [CrossRef]
  72. Peschardt, K.K.; Schipperijn, J.; Stigsdotter, U.K. Use of small public urban green spaces (SPUGS). Urban For. Urban Green. 2012, 11, 235–244. [Google Scholar] [CrossRef]
Figure 1. Geographical framework of Padua (Italy) and a zoomed-in image on different types of urban fabric of the city.
Figure 1. Geographical framework of Padua (Italy) and a zoomed-in image on different types of urban fabric of the city.
Ijgi 10 00538 g001
Figure 2. Location of the macro-sample areas within the municipality borders of Padua adopted as data validation.
Figure 2. Location of the macro-sample areas within the municipality borders of Padua adopted as data validation.
Ijgi 10 00538 g002
Figure 3. Overall flowchart showing the main methodologies (black boxes), data inputs (yellow boxes), urban green spaces validation (green box), and outputs (red boxes).
Figure 3. Overall flowchart showing the main methodologies (black boxes), data inputs (yellow boxes), urban green spaces validation (green box), and outputs (red boxes).
Ijgi 10 00538 g003
Figure 4. Urban green spaces in Padua by NDVI extraction, before integration with municipal topographic database.
Figure 4. Urban green spaces in Padua by NDVI extraction, before integration with municipal topographic database.
Ijgi 10 00538 g004
Figure 5. Comparison of green surfaces obtained from NDVI and from previous results in the four sample macro-areas [42,43]: (a) Brentelle neighbourhood resulted by NDVI analysis (macro-area 1); (b) Brentelle neighbourhood resulted by visual analysis (macro-area 1); (c) Forcellini neighbourhood resulted by NDVI analysis (macro-area 2); (d) Forcellini neighbourhood resulted by visual analysis (macro-area 2); (e) San Lazzaro neighbourhood resulted by NDVI analysis (macro-area 3); (f) San Lazzaro neighbourhood resulted by visual analysis (macro-area 3); (g) Basso Isonzo neighbourhood resulted by NDVI analysis (macro-area 4); (h) Basso Isonzo neighbourhood resulted by visual analysis (macro-area 1).
Figure 5. Comparison of green surfaces obtained from NDVI and from previous results in the four sample macro-areas [42,43]: (a) Brentelle neighbourhood resulted by NDVI analysis (macro-area 1); (b) Brentelle neighbourhood resulted by visual analysis (macro-area 1); (c) Forcellini neighbourhood resulted by NDVI analysis (macro-area 2); (d) Forcellini neighbourhood resulted by visual analysis (macro-area 2); (e) San Lazzaro neighbourhood resulted by NDVI analysis (macro-area 3); (f) San Lazzaro neighbourhood resulted by visual analysis (macro-area 3); (g) Basso Isonzo neighbourhood resulted by NDVI analysis (macro-area 4); (h) Basso Isonzo neighbourhood resulted by visual analysis (macro-area 1).
Ijgi 10 00538 g005
Figure 6. Spatial distribution of all urban green spaces (for definition UGS see paragraph 1.3) in the city of Padua, divided by urban units (municipal codes; see 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 6. Spatial distribution of all urban green spaces (for definition UGS see paragraph 1.3) in the city of Padua, divided by urban units (municipal codes; see 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).
Ijgi 10 00538 g006
Figure 7. Scatterplot of the distribution of urban green spaces feature dimensions (y axis, m2) and distance from city centre (x axis, m).
Figure 7. Scatterplot of the distribution of urban green spaces feature dimensions (y axis, m2) and distance from city centre (x axis, m).
Ijgi 10 00538 g007
Figure 8. Map of rural and non-rural urban green spaces in Padua.
Figure 8. Map of rural and non-rural urban green spaces in Padua.
Ijgi 10 00538 g008
Figure 9. Map of public and private urban green spaces in Padua.
Figure 9. Map of public and private urban green spaces in Padua.
Ijgi 10 00538 g009
Figure 10. Map of public urban green spaces: municipal green spaces (light green), non-municipal green spaces (green).
Figure 10. Map of public urban green spaces: municipal green spaces (light green), non-municipal green spaces (green).
Ijgi 10 00538 g010
Figure 11. Population in Padua urban units. Urban units are represented in the map by their municipal codes (see Table 11 for the corresponding urban units).
Figure 11. Population in Padua urban units. Urban units are represented in the map by their municipal codes (see Table 11 for the corresponding urban units).
Ijgi 10 00538 g011
Figure 12. Total, public, and municipal UGS (x axis per capita m2) in Padua urban units (y axis).
Figure 12. Total, public, and municipal UGS (x axis per capita m2) in Padua urban units (y axis).
Ijgi 10 00538 g012
Figure 13. Spatial analyses of urban green space and population at sub-neighbourhood urban unit levels: (a) map of total urban green spaces; (b) map of urban green spaces per capita; (c) map of municipal urban green spaces; (d) map of municipal urban green spaces per capita.
Figure 13. Spatial analyses of urban green space and population at sub-neighbourhood urban unit levels: (a) map of total urban green spaces; (b) map of urban green spaces per capita; (c) map of municipal urban green spaces; (d) map of municipal urban green spaces per capita.
Ijgi 10 00538 g013
Figure 14. Density classes for population and urban green space categories in the urban units of Padua, based on per capita (PC) values displayed in Table 11 and subdivided into quintiles. Units a–d are quoted as samples for outlining UGS management strategies in Section 3.4.
Figure 14. Density classes for population and urban green space categories in the urban units of Padua, based on per capita (PC) values displayed in Table 11 and subdivided into quintiles. Units a–d are quoted as samples for outlining UGS management strategies in Section 3.4.
Ijgi 10 00538 g014
Table 1. Main input spatial data used and their sources.
Table 1. Main input spatial data used and their sources.
NameDescriptionUseSource
Orthophotos20 cm/pixel resolution, multiband (RGB-NIR), June–July 2015NDVI calculationAGEA/Veneto Region (2015)
Topographic database (DB)Digital land use/land cover mapIntegration and refinement of UGS from NDVI, rural–urban classificationPadua Municipality
UGS map of four sample areasUGS mapped by expert photointerpretation on orthophotos (ground truth) of four areas of PaduaUGS from NDVI + Topographic DB result validationPeroni et al., 2020
Cadastral database (DB)Cadastral DB of public and municipal propertiesProperty (public, municipal, private) classificationPadua Municipality
Table 9. Classification of green areas in Padua.
Table 9. Classification of green areas in Padua.
Green Area CategoryAreaPercentage of Total Green Areas
Total52.23 km2100%
Rural28.80 km255.14%
Non-rural23.43 km244.86%
Public10.25 km219.62%
Private41.98 km280.38%
Municipal5.02 km29.61%
Non-municipal47.21 km290.39%
Table 10. Matrix of public–private and rural–non-rural green areas.
Table 10. Matrix of public–private and rural–non-rural green areas.
PropertyTotal
(km2)
Rural Green
(km2)
Non-Rural Green (km2)Percentage of Rural Green (%)Percentage of Non-Rural Green (%)
Municipal UGS 5.031.223.8124.2575.75
Public, non-municipal UGS 5.221.713.5132.7667.24
Private UGS41.9825.8916.0961.6738.33
Table 11. Areas and per capita (PC) values of green space categories in the urban units of Padua.
Table 11. Areas and per capita (PC) values of green space categories in the urban units of Padua.
Urban Unit CodeUrban UnitPop.Area
(km2)
Pop. Density
per km2
Total
UGS
(m2 PC)
Rural
UGS
(m2 PC)
Public
UGS
(m2 PC)
Private
UGS (m2 PC)
Municipal
UGS
(m2 PC)
1.1 Piazze68930.808616.2513.880.192.4411.451.11
1.2 Savonarola65601.195535.8663.281.0324.8038.4813.14
1.3 Santo - Portello69751.634292.3165.110.6931.7633.3517.05
1.4 Prato della Valle32720.764305.2668.011.3027.8540.1615.51
1.5 Città Giardino41680.785371.1364.862.9916.6348.2212.89
2 Sacra Famiglia73272.772640.36231.8383.7490.61141.2238.14
3 San Giuseppe74631.246023.4140.760.306.2434.523.83
4.1 Porta Trento Sud24680.633936.2059.340.006.0253.324.77
4.2 Porta Trento Nord6250.461346.98457.83184.84123.97333.86123.24
5.1 Fiera20941.022061.0296.8613.7746.3150.5543.55
5.2 Stazione22740.832749.7077.2411.1237.0940.159.98
6 Stanga37511.402671.65117.3919.1358.0859.3148.55
7 Forcellini98362.663697.74132.6331.7041.3091.3323.10
8 Madonna Pellegrina109642.234909.9991.0222.4526.5264.5118.50
9 Sant’Osvaldo66461.086148.0168.3610.779.7558.617.03
10 Voltabarozzo52722.072549.32227.63130.8820.90206.7314.85
11 SS. Crocefisso46122.441887.07365.24221.5946.17319.0713.80
12 Salboro25964.71551.401503.081274.3294.641408.4432.94
13 Guizza127704.253007.54203.04135.6023.96179.0723.02
14 Mandria102488.921149.27662.69461.4548.45614.2425.39
15 Brusegana72923.572041.43319.30176.9995.66223.756.49
16 Cave42102.072030.87255.52126.3975.43180.0972.55
17 Brentelle42742.601645.11425.47236.5745.74379.8425.76
18 Sant’ignazio37911.372775.26242.66156.4633.00209.6631.21
19 Montà11960.911309.97613.94416.8175.28538.6668.88
20 Ponterotto27682.82982.95801.02607.42119.33681.7053.08
21 Sacro Cuore49034.96988.71767.44520.22133.60633.8483.24
22 Altichiero41113.541162.94625.81388.6673.02552.8052.34
23 Pontevigodarzere53021.912771.56205.12108.1131.42173.7013.98
24 San Carlo15,0442.236761.3553.6910.2916.7536.9415.11
25.1Arcella15,9442.257092.5340.992.568.7132.287.60
25.2San Bellino34600.3410267.0631.790.303.1828.612.74
26 Mortise65031.893438.92153.6372.9045.68107.9530.95
27 Torre44963.071466.41467.90288.9840.02427.8827.17
28 San Lazzaro18383.44534.61669.18270.77311.46357.72108.09
29 Ponte di Brenta34701.272727.99150.8821.2317.74133.1413.22
30.1Zona Industriale5758.0771.294096.681464.421888.522208.1755.48
30.2Terranegra2421.18205.613541.261941.741749.261792.00472.61
31 Camin40522.251800.09377.63261.1644.56333.0731.92
32 Granze8771.73505.771023.89647.74311.21712.6950.32
PADUA211,16293.312263.09247.35136.3948.54198.8023.77
Table 12. Possible planning strategies for the UGS categories classified in the study.
Table 12. Possible planning strategies for the UGS categories classified in the study.
PropertyStrategies for Rural UGSStrategies for Non-Rural UGS
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
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pristeri, G.; Peroni, F.; Pappalardo, S.E.; Codato, D.; Masi, A.; De Marchi, M. Whose Urban Green? Mapping and Classifying Public and Private Green Spaces in Padua for Spatial Planning Policies. ISPRS Int. J. Geo-Inf. 2021, 10, 538. https://doi.org/10.3390/ijgi10080538

AMA Style

Pristeri G, Peroni F, Pappalardo SE, Codato D, Masi A, De Marchi M. Whose Urban Green? Mapping and Classifying Public and Private Green Spaces in Padua for Spatial Planning Policies. ISPRS International Journal of Geo-Information. 2021; 10(8):538. https://doi.org/10.3390/ijgi10080538

Chicago/Turabian Style

Pristeri, Guglielmo, Francesca Peroni, Salvatore Eugenio Pappalardo, Daniele Codato, Antonio Masi, and Massimo De Marchi. 2021. "Whose Urban Green? Mapping and Classifying Public and Private Green Spaces in Padua for Spatial Planning Policies" ISPRS International Journal of Geo-Information 10, no. 8: 538. https://doi.org/10.3390/ijgi10080538

APA Style

Pristeri, G., Peroni, F., Pappalardo, S. E., Codato, D., Masi, A., & De Marchi, M. (2021). Whose Urban Green? Mapping and Classifying Public and Private Green Spaces in Padua for Spatial Planning Policies. ISPRS International Journal of Geo-Information, 10(8), 538. https://doi.org/10.3390/ijgi10080538

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop