Next Article in Journal
Identification and Risk Assessment of Priority Control Organic Pollutants in Groundwater in the Junggar Basin in Xinjiang, P.R. China
Next Article in Special Issue
Blackland Conservation and Utilization, Carbon Storage and Ecological Risk in Green Space: A Case Study from Heilongjiang Province in China
Previous Article in Journal
Potential Clinical Applications of Ozone Therapy in Dental Specialties—A Literature Review, Supported by Own Observations
Previous Article in Special Issue
Does Increasing Farm Plot Size Influence the Visual Quality of Everyday Agricultural Landscapes?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Analysis of the Spatial Development of European Cities Based on Their Geometry and the CORINE Land Cover (CLC) Database

1
Department of Geoinformation and Cartography, Faculty of Geoengineering, University of Warmia and Mazury in Olsztyn, Prawocheńskiego 15, 10-720 Olsztyn, Poland
2
Department of Socio-Economic Geography, Faculty of Geoengineering, University of Warmia and Mazury in Olsztyn, Prawocheńskiego 15, 10-720 Olsztyn, Poland
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(3), 2049; https://doi.org/10.3390/ijerph20032049
Submission received: 18 December 2022 / Revised: 13 January 2023 / Accepted: 18 January 2023 / Published: 22 January 2023
(This article belongs to the Special Issue Land Management for Territorial Spatial Planning)

Abstract

:
The study demonstrated that the rate of spatial development is correlated with its fractal dimension. The presented results indicate that the fractal dimension can be a useful tool for describing different phases of urban development. Therefore, the formulated research hypothesis states that the fractal dimension of cities’ external boundaries is correlated with the rate of spatial development in urban areas. The above implies that the higher the rate of spatial development, the smoother the external boundaries of urban investment. Rapidly developing cities contribute to considerable changes in land management, in particular in municipalities surrounding the urban core. Urban development processes often induce negative changes in land management and contribute to chaotic and unplanned development. To address these problems, new methods are being developed for modeling and predicting the rate of changes in transitional zones between urban and rural areas. These processes are particularly pronounced in urban space, whose expansion proceeds at an uneven pace. The aim of this study was to propose a method for describing urbanization processes that are based on the dependence between the urban growth rate, the fractal dimension, and basic geometric parameters, such as city area and the length of city boundaries. Based on the calculated changes in the values of these parameters, a classification system was proposed to identify distinctive phases of urban development. The study revealed that land cover databases are highly useful for such analyses. The study was conducted on 58 medium-size European cities with a population of up to 300,000, including France, Germany, Italy, Poland, and Croatia. The study demonstrated that the fractal dimension and the basic geometric parameters of urban boundaries are significantly correlated with the rate of the spatial development of cities. The proposed indicators can be used to describe the spatial development of urban areas and the rate of urban growth. The development of the analyzed cities was modeled with the use of CORINE Land Cover (CLC) data for 2000–2006–2012–2018, acquired under the EU Copernicus program.

1. Introduction

1.1. Spatial Development of Cities

Cities are the key areas in the development of every civilization, and they attract considerable attention in various fields of scientific inquiry [1,2,3]. However, urban development in Europe and the world is accompanied by negative phenomena. The accumulation of uncontrolled social, economic, and spatial processes exacerbates the scale of the encountered problems and affects a growing number of people residing in cities and the surrounding areas [4,5,6,7]. Increased mobility and unlimited access to mobile devices and the Internet contributes to the growth of cities and the urban population [8].
Urban development analyses indicate that urban expansion contributes to adverse changes in land use and management [9]. The rapid development of cities leads to urban crowding as well as land, water, and air pollution, which poses a public health risk [10]. Cities generally differ in structure and design, but upon closer inspection, certain similarities can be observed in urban forms [11,12].
Wilson, Ware, and Ware described three types of processes in the development of urban areas: infill development, expansion of the existing areas, and the creation of new areas that are located remotely from the urban core. The new settlements can be isolated, can be linear, or can form clusters [13].
One approach to identifying new urban areas postulates that infill development takes place when at least 40% of the transformed land is surrounded by urbanized areas [14]. Most of these areas are equipped with basic technical infrastructure, such as roads and public utilities, and with social infrastructure [13]. According to Ellman, the existing infrastructure is the main prerequisite for infill development [15]. This is a rational approach that accounts for the importance of infrastructure in the process of building cities. The emergence of new settlements puts pressure on the local authorities to deliver the required infrastructure [15,16]. From the point of view of landscape preservation and environmental protection, infill development contributes to the loss of open areas and urban ventilation corridors [17]. This approach confirms that urban development is strictly linked with a city’s geometry; therefore, it can be described with the use of geometric parameters such as surface area, a perimeter, and the fractal dimension.
Urban expansion takes place when the percentage of urbanized land surrounding a given area does not exceed 50% [18]. This approach focuses mainly on changes in land use, and urban expansion occurs when urban land-use types become predominant in each territory [19]. The conversion of undeveloped land to urban land can involve two processes [20]. The first process includes changes in the transitional zone, which is characterized by low availability of infrastructure. The second process involves changes outside the transitional zone, where new settlements are highly dispersed, and they are referred to as urban villages. Therefore, the existing land-use types are indicative of the stage of urban development, which implies that land cover data can also be used to analyze the spatial expansion of cities [21,22].
Urban expansion has far-reaching consequences for the functioning of entire ecosystems as well as the residents of urban, transitional, and rural areas [23]. This type of urban development is referred to as boundary or border development, where urban areas expand in parallel bands with an outer edge [1]. The above phenomenon is linked with the concept of the urban growth boundary (UGB), which has been extensively explored in the literature [24,25,26,27]. New settlements can be also remotely established from the existing urban centers [20]. The resulting changes in land use take place outside the transitional zone, and they are termed as leapfrog development [9]. The establishment of new settlements at a certain distance from the urban core is also known as isolated development, and it is characteristic of rural areas where new additions to the existing clusters of buildings match the local architectural style.

1.2. Fractal Dimension in Analyses of Urban Development

Urbanization is one of the most dynamic global processes. Urbanization drives social and economic growth, which is why cities continue to expand their area and population. Urbanization induces profound changes in space, and it gradually modifies land-use structures. Urban spaces and changes in areas that are directly subjected to urbanization pressure should be monitored. New analytical methods and techniques are needed to reliably assess the nature and rate of urbanization, in particular the spatial expansion of cities [28].
An analysis of the literature indicates that the boundaries of a city rarely coincide with its administrative boundaries and are difficult to define. The course and geometry of city boundaries play a key role in urban development, which is why they attract considerable research interest. Urban areas have a complex and extensive spatial structure, and the distance from the city center to city boundaries is difficult to determine. The boundaries of urban areas play a key role in urban research. However, the length of urban boundaries and the area enclosed by these boundaries cannot be objectively defined [29,30,31,32]. Therefore, the length of the urban boundary has to be accurately measured to determine the size of the urban population.
Increasingly accurate models of urban development are being developed to address progressing urbanization, which poses one of the greatest civilizational problems in the world. Research into local factors that promote urban development indicates that cities are self-organizing structures [33]. One of the most interesting approaches to analyzing urban development and changes in urban boundaries posits that cities should be regarded as fractal units [34]. A fractal is a geometrical figure whose individual parts are similar to the whole [35,36]. By definition, a fractal is a set for which the Hausdorff–Besicovitch dimension strictly exceeds the topological dimension [37,38]. A fractal is scaleless, and it cannot be described with traditional units such as length, surface area, volume, or density.
The fractal dimensions of cities are generally defined in two-dimensional space on the basis of digital maps and remote-sensing images [9,10,20]. If an urban fractal is defined in two-dimensional space, the urban area and the urban boundary can be described with the use of fractal dimensions. The urban boundary can be regarded as a fractal line [36,39,40,41]. A closed curve representing the urban boundary is defined as the urban periphery, where the Euclidean pattern of urban space can be described [37,42]. In practice, the upper boundary and the lower boundary of the fractal dimension of urban space are influenced by the method applied to define the study area. Two approaches can be used to derive the fractal time series of urban growth [43,44]. The first approach relies on a fixed study area [45] and the second approach on a variable study area [46,47].
In the literature, all fractal images represent prefractals, rather than true fractals in a mathematical sense. A true fractal has an infinite number of iterations that can be revealed only in the world of mathematics. A prefractal has a limited number of hierarchical levels with a fractal geometry. Cities are accidental prefractals rather than true fractals because the urban form cannot be described on a characteristic scale [48].
Two approaches have been proposed for determining the fractal dimensions of urban areas. In the first approach, fractal dimensions are calculated for the shape of a city’s boundary, whereas in the second approach, fractal dimensions are linked with the density of urban development [33]. These approaches imply that urban development resembles the growth of two-dimensional particle aggregates [49]. Several concepts have been proposed for modeling the spatial development of cities. A recent approach makes a reference to cluster expansion in terms of statistical physics, and it posits that new objects are not added to the cluster but rather are linked with it [50]. A model describing the development of urban settlements as miniature cities (subclusters) relies on the assumption that urban growth fuels further growth. According to the various approaches to urban development in the literature, urban growth should be regarded as a process of organic growth that begins in the urban core and spreads to miniature cities outside the urban core [28,51].
Spatial databases containing information about land cover as well as advanced techniques for processing and modeling spatial data are vast sources of knowledge, and they can be deployed to develop new tools for identifying and monitoring urban sprawl [52,53,54]. The CORINE Land Cover (CLC) database supports broad spatial analyses because the data describing land cover in Europe are characterized by spatial continuity and enable the nonambiguous identification of various land-use types. Most importantly, CLC databases are regularly developed, which facilitates analyses of the dynamics and rate of changes and supports forecasting. Other sources of data, such as the Urban Atlas and the Global Human Settlement Layer, are less versatile in this respect [15]. According to the literature, the CLC is a far more useful resource for small-scale studies, but it is a less reliable tool for analyses conducted on a larger scale [17]. The CLC project, initiated in 1985 and updated in 2000, 2006, 2012, and 2018, provides information on the EEA member countries (39 countries, EEA39). The minimum mapping unit (MMU) in the CLC is 25 hectares (ha) for aerial phenomena and a minimum width of 100 m for linear phenomena. Land cover is mapped mainly on the basis of the visual interpretation of high-resolution satellite images. The CLC data set consists of 44 land cover classes grouped into five main categories: artificial surfaces, agricultural areas, forest and seminatural areas, wetlands, and water bodies [55,56,57]. Studies that rely on CLC data have certain limitations, such as the detailed nature of input data and interpretation methods and, consequently, a high degree of generalization. In areas characterized by considerable land fragmentation, the results can be generalized into dominant land-use types, which leads to a certain loss of information.
In this study, the development of urbanized areas was described and modeled with the use of the fractal dimension calculated on the basis of CLC data. The main aim of this study was to propose a method for describing the spatial development of a city. Basic geometric parameters were computed in selected European cities, and changes in these parameters were analyzed over time to propose a formula for describing and classifying cities on the basis of their stage of evolution. This research is important because of the importance of monitoring urbanization processes. The adopted method makes it possible to identify the trends and dynamics of the processes taking place, providing an important indicator for optimizing the spatial policy of the cities under analysis. The proposed procedure may contribute to extending the scope of research in spatial terms (selection of analyzed cities), in temporal terms (successive time intervals), and in terms of indicators (the analysis may be extended by adding subsequent indicators, e.g., number of inhabitants). The application of new elements in the analysis may contribute to the enrichment of empirical research in the identification of urban development processes. The spatial expansion of cities, the fractal dimension, and CLC data are described in the Introduction. The cities selected for the study and the procedure for calculating the fractal dimensions of city boundaries and the rate of urban development are described in the Section 2. The indicators calculated for the analyzed cities, the increase in the urban area, and the fractal dimension of the analyzed cities are presented, and the examined cities are classified on the basis of the stages of urban expansion in different time intervals in the Section 3. The results are discussed in the Section 4, and the recommendations for future research are formulated in the Section 5.

2. Materials and Methods

2.1. Study Area

The rate of spatial development was analyzed in 58 European cities, which were selected for the study on the basis of the availability of CLC data for 2000, 2006, 2012, and 2018 (Figure 1). The evaluated medium-size cities had a population of up to 300,000 at the beginning of the analyzed period. Owing to considerable differences in the social and economic development of European cities, only cities in the European Union, including France, Germany, Italy, Poland, United Kingdom, and Croatia, were selected for the study. In line with the New Urban Agenda [58], the study was conducted in European cities with the greatest potential to make Europe a global reference point for identifying, experimenting, and applying solutions to future urban challenges [59]. Land cover data for the examined cities were acquired from the CLC database for 2000, 2006, 2012, and 2018. CORINE Land Cover data are publicly available, and they are stored and managed by the respective national authorities in the EU [60]. In addition, the use of data from the CLC database, owing to the different levels of accuracy in the data sets, where older data sets have less accuracy than newer ones, required visual verification of the different land-use types occurring for each urban area. The authors in this paper wanted to accurately determine geometric indicators for urban areas, so they focused on 58 European cities, and they bore in mind that historical CLC data sets have some limitations, such as the detailed nature of the input data and methods of interpretation.
The cities located in different European countries and characterized by various spatial attributes were compared to verify the fractal dimension as an indicator for evaluating the rate of urban development. The analyzed set included coastal cities (Aberdeen, Rijeka, Trieste) as well as inland cities (Cheltenham, Strasbourg, Crawley). The limiting effect of rivers’ intersecting urban areas was also considered. In the past, rivers played an important role in the establishment of urban settlements, and they presently pose a significant barrier to the expansion of local transportation systems. Therefore, cities intersected by rivers (Verona, Norwich, Nantes) as well as cities that were not built around water courses (Oviedo, Białystok, Erfurt) were included in the analysis. The social, economic, legal, and spatial aspects of urban development differ across European countries; therefore, the evaluated cities most probably differ according to other distinguishing features. The fractal dimension was determined to objectivize the process of describing urban development.
The available sources of land cover data were analyzed for the needs of this study. The CLC inventory was selected thanks to its availability, broad coverage, and confirmed usefulness in evaluations of urban development [61,62,63,64]. The boundaries of urban areas had to be determined in detail; therefore, land use was described on the basis of level 3 CLC data. The land cover map for selected cities is presented in Figure 2. In the studied group of 48 cities, land cover was determined at two points in time on the basis of CLC data for 2000, 2006, 2012, and 2018. This approach supported the identification of significant changes in the fractal dimension of city boundaries. The boundaries of the studied cities were identified in detail, which facilitated analyses of the rate of urbanization.

2.2. Procedure for Calculating the Fractal Dimension of Boundaries and the Rate of Urban Development

This study was inspired by the referenced works and by research into the dependence between the fractal dimension of the external boundaries of various land-use types and the degree of anthropogenic pressure [64].
For the needs of the study, urban areas were selected by determining which of the 58 land cover classes in the CLC inventory should be used for further analysis. The main land cover classes characteristic of urban areas were described on the basis of the principles for identifying urban morphological zones (UMZs). In the CLC inventory, level 1 landscape patterns represented artificial surfaces. In some cities, level 3 patterns denoting forests and seminatural areas were also classified as urban areas. This was the case when forests and seminatural ecosystems were surrounded in their entirety by artificial surfaces. The selected areas were then aggregated into a single group. In the following step, urban boundaries were verified on the basis of high-resolution satellite images for the analyzed time periods [65]. The procedure was repeated for CLC 2000, 2006, 2012, and 2018 data. The urban boundaries identified on the basis of CLC data for 2000 2006, 2012, and 2018 were compared to identify changes in urbanization processes in the evaluated cities (Figure 3).
In this study, the fractal dimension of urban development was analyzed on the basis of research into land cover and land-use types [51,66]. The fractal dimension of urban areas was determined with the box-counting method. The hierarchical grid method proposed by De Cola was used, where a grid of boxes is superimposed on a fractal, and the objects intersected by the fractal are counted [67]. In this approach, the number of boxes represents the fractal’s surface area, and the combined length of box edges is the fractal’s perimeter. The size of the boxes in the grid overlay is decreased in successive iterations. The box-counting method is used to identify changes in the fractal’s perimeter when the length of the box edges is decreased [68]. In urban structures, the number of elements in successive iterations is not constant; therefore, the box-counting dimension is defined as a boundary value where the length of the box sides decreases toward zero. According to Equation (1),
O = c S D f
where
  • O is the perimeter;
  • c is the shape constant;
  • S is the surface area;
  • Df is the box-counting fractal dimension.
Therefore, the box-counting dimension is defined by Equation (2):
D f = log O log S
The box-counting dimension is determined as the slope of a regression line in a graph presenting the correlations between log values. In practice, real-world objects are usually modeled with the box-counting method and the compass-walking method. The relevant technique has to be selected when planning the experiment because it will influence the preparation of input data and the choice of calculation methods.
The obtained results were compared with CLC data, which were acquired under the Copernicus program. The CLC inventory was used because not all urbanized areas are contained within the administrative boundaries of cities. As a result, these areas are not taken into account in analyses of the spatial development of cities. There is no single harmonized definition of a city and a functional urban area, which considerably impedes analyses of European cities [19]. For this reason, the boundaries of urban areas were identified by introducing the urban morphological zone (UMZ) concept, which denotes continuous urban fabric with different population density. An urban morphological zone can be defined as “a set of urban areas laying less than 200 m apart”. These urban areas are defined on the basis of land cover classes describing the urban tissue and function [69]. The UMZ database was created by the European Environment Agency (EEA) on the basis of CLC data and automated methods for defining the boundaries of urban agglomerations. Maximum distance is the key criterion. Urban morphological zones comprise areas of dense urban development (separated by a maximum distance of 200 m), and they consist of the following land cover classes: continuous or discontinuous urban fabric, industrial or commercial units, green urban areas, selected forest areas, port areas, airports, sport and leisure facilities, and road and railway networks [70]. The adopted solution for identifying urban areas relied on a set of spatial data developed on the basis of the CLC inventory and the assumptions for defining UMZs [71].
The calculation of the fractal dimension and the determination of the rate of urban development in the analyzed cities between 2000–2006, 2006–2012, and 2012–2018 were important steps in the adopted procedure. A standard urban development model was adopted to describe population growth and human behaviors relating to optimal decision-making within the allocation of time and according to required effort and resources. There are many examples linking the applied model with economic, social, and health-related behaviors, including descriptions of the exponential growth of urban populations [72,73].
The fractal dimension (Df) was calculated with the box-counting method, also known as the area-perimeter method. In the adopted approach, the variables were the surface area of the fragments of the grids covering the studied object (S), which is described by the number of boxes, and the object’s perimeter (O), which is expressed by an equal number of box edges in a given fragment of the grid. All calculations were performed in the QGIS program. A geographic database containing information about land-use types in the analyzed years and classes of objects representing grid boxes with specific dimensions was generated. An algorithm was created for calculating the surface area and the perimeter of grid boxes. The fractal dimension was calculated on the basis of the values computed for boxes with a minimum surface area of 625 m2 and a maximum surface area of 625 km2. The box area was determined in 30 steps. The algorithm was tested and applied to boxes in each size category on the basis of CLC data for 2000, 2006, 2012, and 2018. The increase in the fractal dimension (ΔDf) was calculated with the use of Equation (3):
Δ D f = ( D 2 i D 1 i D 2 i )
where
  • ΔDf is the increase in the fractal dimension;
  • D1 is the fractal dimension in time t;
  • D2 is the fractal dimension in time t + 1;
  • i is the number for the analyzed city.
The increase in the fractal dimension of cities’ external boundaries is correlated with the rate of spatial development in urban areas. Higher values of increase in the fractal dimension imply a higher rate of spatial development, creating a ragged and dendritic urban fabric. On the other hand, the lower values of the increase indicator are related mostly to filling the areas between the main traffic routes and areas vulnerable to the pressures of urbanization processes while avoiding natural and anthropogenic obstacles. The study was expanded to include an analysis of changes in the area and perimeter of geometric figures that represent urban land-use types and were used to calculate the fractal dimension. The rate of changes in figure geometry and the fractal dimension was determined in successive time intervals: 2000–2006, 2006–2012, and 2012–2018.
The spatial development of cities (increase in urban area) in each time interval (2000–2006, 2006–2012, and 2012–2018) was calculated with the use of Equation (4):
Δ A i = ( A 2 i A 1 i A 2 i )
where
  • ΔA i is the urban growth rate i city;
  • A1 is the urban area in time t;
  • A2 is the urban area in time t + 1,
  • i is the number for the analyzed city.
The values of the urban growth Δ A i rate indicator can range from 0 to several points. Values close to 0 indicate that a city’s area did not increase in the analyzed period, which implies that its development was inhibited. This is often the case when urban development is constrained by natural or artificial barriers. Values higher than 1 are indicative of dynamic urban development.
The increase in the perimeter of the geometric figure, i.e., the increase in the length of city boundaries, was calculated with the use of Equation (5):
Δ P i = ( P 2 i P 1 i P 2 i )
where
  • ΔAPi is the increase in the length of i city boundaries;
  • P1 is the boundary length in time t;
  • P2 is the boundary length in time t + 1;
  • i is the number for the analyzed city.
The increase in boundary length denotes the degree of figure filling. Ragged boundaries are longer, and they could imply that a rapidly developing city has annexed the surrounding areas. In turn, smoother boundaries are shorter, which suggests that the urban fabric is more compact and urban growth relies mainly on infill development (for example, in vacant areas between transport routes).
In order to organize the steps taken in the procedure for determining the phases of urban development, we can identify three stages. The first involves the determination of the study area. Within the framework of the activities undertaken, it is necessary to make a selection of cities and characterize the data sets that allow the identification of urban processes. The second stage involves data analysis and the determination of boundaries and classes characteristic of urbanized areas. The next step assumes the calculation of basic geometrical parameters (area, perimeter) and the fractal dimension. According to the acquired data, an increase in the above indicators is determined for the available time intervals. The last step involves the verification of the obtained results and the determination of a pattern for classifying the stages of urban development for every city. The respective increases in the fractal dimension (ΔDf), area (ΔA), and boundary (perimeter) length (ΔP) were calculated in the next step. The relationships between the calculated values were analyzed in all time intervals, but significant correlations were not observed. Therefore, it was assumed that the calculated increases can provide additional information on spatial expansion. An urban expansion matrix was developed on the basis of the calculated increase in the analyzed parameters and literature data (Table 1). The examined cities were assigned to different urban expansion classes on the basis of the calculated changes in perimeter, area, and the fractal dimension. All possible change combinations in 2000–2006, 2006–2012, and 2012–2018 were determined, and the increase (↑) and decrease (↓) in the analyzed indicators were calculated. The resulting urban expansion matrix was composed of eight classes, and it was used to determine the rate of spatial changes in each city. The classification process based on the urban expansion matrix is presented in Table 1.
On the basis of the calculated changes in perimeter, area, and the fractal dimension, the analyzed cities were divided into eight classes, denoting different phases (stages) of spatial expansion. The proposed classes are described below:
Class 8—
in this rapid urban expansion, the city expands in an uncontrolled manner in all directions. New urban fabric is ragged and dendritic.
Class 7—
urban expansion involves mainly infill development between major transport routes, and obstacles are bypassed. Areas that are relatively sensitive to urbanization pressure are annexed by the city. City boundaries are more compact and regular.
Class 6—
urban expansion takes place only in the vicinity of urban infrastructure. Urban boundaries are smoother, but the city has a dendritic shape.
Class 5—
urban expansion takes place in areas that are more resistant to anthropogenic pressure. Urban boundaries are smoother and more compact.
Class 4—
urbanized areas disappear. City boundaries are more ragged and dendritic in shape. Examples of the above include reclaimed areas with urban infrastructure, in particular linear infrastructure.
Class 3—
in urban regression, urbanized areas disappear from the urban periphery. Urban boundaries are compact but ragged.
Class 2—
in the regression phase, urbanization regresses, and urbanized areas are found mainly in the vicinity of linear infrastructure. Urban boundaries evolve into a dendritic shape.
Class 1—
full regression occurs in all directions.
The developed classification system was used to describe the phases of urban development in the examined cities in 2000–2006, 2006–2012, and 2012–2018. The proposed method was applied to examine the intensity and rate of changes in urban development.

3. Results

The study was conducted in 58 European cities. Their surface area, perimeter, and fractal dimension (Df) calculated on the basis of CLC data for 2000, 2006, 2012, and 2018 are presented in Table A1. The urban growth rate indicator was calculated for all cities on the basis of the increase in their area, perimeter, and fractal dimension in 2000–2006, 2006–2012, and 2012–2018 (Table 2).
In terms of the increase in city area (ΔA), the values of the urban growth rate indicator considerably differed across the examined time intervals. In 2000–2006, the greatest increase in urban area was noted in Verona (0.673), Alicante (0.609), and Rijeka (0.480). In 2006–2012, the greatest increase in urban area was observed in Monchengladbach (0.578), Chemnitz (0.346), and Almeria (0.2410). These cities are popular tourist destinations with similar locations and geographical features, which suggests that these factors played a role in their development in the analyzed periods. Similar observations were made in the last time interval. In 2012–2018, the greatest increase in urban area was noted in Västerås (0.353), Burgos (0.323), and Tampere (0.050). An analysis of changes in the fractal dimension (ΔDf) over time produced equally interesting results in 2000–2006: the greatest increase in the fractal dimension was observed in Burgos (0.102), Vigo (0.056), and Trieste (0.041). In 2006–2012, the greatest increase in the fractal dimension was noted in Munster (0.020), Santander (0.012), and Osijek (0.010). In 2012–2018, the greatest increase in the fractal dimension was observed in Västerås (0.061), Burgos (0.030), and Messina (0.009). The greatest increase in city perimeter (ΔP) was noted in Vigo (1.338), Verona (1.216), and Burgos (0.946) in 2000–2006; in Chemnitz (0.369), Monchengladbach (0.362), and Santander (0.242) in 2006–2012; and in Västerås (2.557), Burgos (0.274), and Verona (0.106) in 2012–2018. Moreover, the greatest increase in the analyzed parameters was observed in cities that are popular tourist destinations and, in two cases (Chemnitz and Munster), in rapidly growing industrial centers. In some cases, the analyzed parameters increased at a steady rate, which confirms the hypothesis that urban expansion proceeds at an uneven rate.
In line with the adopted procedure, the analyzed cities were divided into eight classes describing different phases (stages) of spatial development in the examined time intervals. The examined cities were allocated to different classes on the basis of the previously calculated increase (↑) or decrease (↓) in city area, perimeter, and the fractal dimension in 2000–2006, 2006–2012, and 2012–2018 (Table 1 and Table 2). The resulting classification of the analyzed cities is presented in Table 3.
The proposed classification system presents the stages of urban expansion in 58 European cities in 2000–2006, 2006–2012, and 2012–2018. Rapid changes in the rate of spatial development can be observed across the examined time intervals. Some cities were assigned to radically different classes in each period, including Alicante (8-1-8), whereas a rapid decrease in the rate of spatial development was observed in Split, Częstochowa, and Toruń (7-8-1). Between 2000 and 2018, a high and steady rate of urban development was observed in Białystok, Derby, Freiburg, Rijeka, Northampton, and Norwich, where the values of all examined parameters increased. Aberdeen, Almeria, Erfurt, Geneva, Trieste, Trondheim, and Verona were characterized by a relatively steady rate of urban expansion, and none of these cities moved up or down by more than one or two classes in the analyzed period. In the remaining cities, spatial expansion proceeded in rapid spurts. The phases of spatial expansion in the studied cities are presented in maps in Figure 4.
Figure 4. Phases of spatial expansion in the analyzed cities, labeled as in Table 3 and Table 1.
Figure 4. Phases of spatial expansion in the analyzed cities, labeled as in Table 3 and Table 1.
Ijerph 20 02049 g004

4. Discussion

Geospatial data, including CLC data, are a rich source of knowledge, and they can be used to develop new tools for identifying and monitoring urbanization processes. In areas that are most exposed to urbanization pressure, the degree and rate of urban expansion can be most effectively monitored by analyzing land cover data. However, this approach also has many limitations, which became apparent at the stage of calculating the boundaries and the fractal dimension of the cities.
The maximum city area for calculating the fractal dimension in the box-counting method was set at a maximum side length of 25,000 m, and the minimum side length was set at 25 m. In the calculations, 30 intermediate values of box side length (min-max) were adopted. Box sizes were compared in all evaluated periods. The comparison for 2018 was based on Sentinel data, which are more accurate (to the nearest 10 m) than other data sets. However, the acquired data were analyzed only in the vector format on the assumption that CLC data for 2000, 2006, 2012, and 2018 can be used to analyze changes in the surface area of the studied cities. Urban morphological zones were established for the needs of the analysis. However, a detailed analysis of the data set prepared on the basis of UMZ guidelines revealed that despite the absence of changes in land use, the same areas were differently classified in the CLC database in the studied periods. For example, the classification of undeveloped areas changed in Alicante (black circle in Figure 5). Undeveloped areas in the CLC database were marked with code 3.2.3 (sclerophyllous vegetation; 3.2 scrub and/or herbaceous vegetation associations; and 3 forest and seminatural areas) (Figure 5a) in 2000 and with code 3.2.1 (natural grassland) in 2006 (Figure 5b).
According to the research conducted, it can be concluded that there are specific cases such as Alicante that require in-depth analysis. The background of the changes that have taken place is related to the development of transport infrastructure and the temporary identification of adjacent areas as built-up areas. The consequence of the processes that have occurred is the infilling of land between the main traffic routes. Such processes affect the dynamics of change in terms of urban development phases. In view of the above, the authors believe that it makes sense to continue the research and extend it to include further time scales that would allow a more accurate identification of the urbanization processes.
In Lubeck, the area marked with code 1.1.2 (discontinuous urban fabric; 1.1 urban fabric; 1 artificial surfaces) in 2000 (black circle in Figure 6a) was marked with code 3.1.1 (broad-leaved forest; 3.1 forest; and 3 forest and seminatural areas) in 2006 (black circle in Figure 6b).
In Milton Keynes, the area marked with code 2.3.1 (pastures) in 2000 (black circle in Figure 7a) was classified as an urban area in 2006 (code 1.4.1 green urban areas) (black circle in Figure 7b) despite the absence of any changes in land use in this period.
These observations suggest that the fractal dimension, area, and perimeter of objects classified as urban areas could have been considerably influenced by the image interpretation method adopted in the selected periods.
The update frequency of land cover data in the vector format is yet another limitation in analyses on urbanization processes. The update cycle for CLC is 6 years. This update frequency appears reasonable for urbanization studies, but it significantly affects the time scale of the analysis. On the basis of the quality and specificity of the available information, land cover data can be reliably compared beginning from 2000, which implies that only four data sets can be used for analysis. This is a relatively small sample, in particular in studies that explore dynamic urbanization processes. The quality of the results is also affected by rapid technological advancements in data acquisition methods, which affects the accuracy of the available information. The above could pose a considerable problem, in particular in comparative analyses.
The rate of urban growth is determined mainly by the urban morphology (urban form) and its relationship with the size of a city, its functions, and the economic, social, technical, and environmental determinants of development. Observations of long-term, continuous urbanization processes generate valuable and sometimes-surprising results. In many cities, the urban growth rate indicator assumed negative values in the studied period. This could imply that many areas that are influenced and transformed by urbanization processes ultimately evolve into nonurban forms. The above suggests that urban development is a more dynamic and bidirectional process than was previously thought. Urban expansion is not always represented by positive values of the urban growth rate indicator, and negative values can be observed in certain stages of the urban life cycle [74]. Urban development or expansion begins in land that is most accessible. Therefore, urban land-use types can be compared to a predator in Lotka–Volterra equations, and they aggressively dominate over the remaining (surrounding) land-use forms [54,55,56]. In this analogy, “weaker” land-use types are more rapidly transformed into the dominant land-use forms. The results of the analysis presented in the above tables produce interesting conclusions concerning the rate of urbanization in the evaluated cities. Cities with the highest rate of urban growth in the analyzed period were in Germany (Freiburg im Breisgau), Croatia (Grad Rijeka), the United Kingdom (Derby, Northampton, and Norwich), and Poland (Białystok), a finding that is largely consistent with the economic trends reported in these countries in the studied period.
Most of the analyzed 58 cities were characterized by positive values of the urban growth rate indicator in each of the three studied time intervals. Negative indicator values were observed in only five cases (in total, 174 indicator values were calculated) in the examined period. This result confirms that urbanization exerts considerable anthropogenic pressure, and areas affected by urbanization rarely regress to other land-use types. The above could also imply that area expansion is not a highly reliable indicator for analyzing the development of areas that have already undergone urbanization. The territorial expansion of cities could be more effectively examined on the basis of the increase in the length of their boundaries and the fractal dimension, which, as demonstrated by this study, are characterized by much greater variation.
The proposed classification method effectively identifies trends in the spatial development of cities [75]. The number of classes is directly associated with the number of parameters adopted for the analysis. The phases (stages) of urbanization and changes in the spatial development of cities were described by assigning cities to different classes on the basis of changes in their area, perimeter, and fractal dimension. High and sudden increases in perimeter and the fractal dimension denote the rapid and often-uncoordinated spatial development of urbanized areas.

5. Conclusions

Analyses of the rate of urban development play key roles in land management in the context of monitoring the growth of cities and the adjacent areas and undertaking preventive measures. The rate of urban development can be accurately evaluated with the use of spatial data. The proposed procedure for evaluating the spatial development of cities requires access to regularly updated land cover data. Land cover data sets and advanced techniques for processing and modeling spatial data are rich sources of knowledge, and they can be used to develop new tools for identifying and monitoring the spatial development of cities. The CLC inventory is characterized by spatial continuity, and it can be reliably used to identify various land-use types and analyze the geometry of cities. CORINE Land Cover data are acquired at regular time intervals; therefore, changes in land use can be reliably monitored. Several key conclusions can be formulated on the basis of the analysis of the external boundaries of urban areas.
The proposed classification of the development phases of medium-size European cities indicates development trends and can be used as a tool for monitoring urban urbanization. On the basis of the analyzed cases, the stability of development of cities in the central part of Europe can be observed. Additionally, we can find individual cases, such as Munster, Vigo, Kiel, and Augsburg, where there is a steady decline in development phases. The authors point out that the usefulness of the tool will be even greater in the future, when it will be possible to carry out analyses for further time intervals. The method adopted also makes it possible to identify specific cases of changes in urban development phases, using Alicante as an example (grades 8-1-8). Such cases require in-depth analyses. It should be noted that they may be due to the processes of intensive development of infrastructure investments taking place. Consequently, in their surroundings, many areas act as servicing areas for the investments, thus being classified as urbanized areas. On the other hand, after the investment has completed, they return to their original uses, being identified as nonurbanized areas. However, in the long term, the completed investment over time induces urbanization processes in the neighboring areas. Consequently, neighboring areas are again identified as urbanized land-use classes.
The results of the conducted analysis indicate that the CLC inventory is a useful resource for describing the rate of urban growth. A high rate of urbanization leads to rapid development regardless of the attributes of built-up land, which can speed up or slow down this process to a varied degree. These findings could suggest that slow urbanization promotes the development of areas with optimal spatial attributes.
The results of this study and the derived conclusions can be used to formulate a preliminary rule concerning the evolution of the external boundaries of cities: cities where urban development proceeds at a faster rate are characterized by smoother external boundaries. Therefore, cities (urban morphological zones) whose spatial development proceeds at a faster rate tend to have ragged external boundaries. The reverse also applies: the fractal dimension of cities’ external boundaries is higher in urban areas characterized by a lower rate of spatial development.
The study demonstrated that urbanization not only increases the area of urbanized areas but also leads to changes in the length and shape of city boundaries. These geometric parameters could play key roles in describing the intensity of urbanization processes. The proposed procedure can be applied in the preliminary stage of identifying the phase of urban development. The results can be used to determine whether cities evolve in a planned manner and whether spatial transformation processes lead to the creation of optimal, compact urban forms.
The identification of the directions and rate of the urban development in Europe and in other parts of the world poses a considerable challenge. Information on potential land use is key to enhancing inclusive and sustainable urbanization. Urban sprawl and inefficient use of land continue to pose a problem, with varying impacts in different contexts. Modern technologies, such as satellite data, support the continuous monitoring of the changes in, standardization of, and protection of citizens’ privacy. Urban planning solutions that make optimal use of the available space will maximize social benefits, support the identification of areas that require careful regulation, and promote forward-looking urban growth strategies.

Author Contributions

Conceptualization, S.C., K.S., T.B., I.C. and A.B.; methodology, S.C., K.S., I.C., T.B. and A.B.; software, S.C. and K.S.; validation, I.C. and A.B.; formal analysis, T.B. and I.C.; investigation, K.S., I.C. and A.B.; resources, S.C. and K.S.; data curation, S.C. and K.S.; writing—original draft preparation, S.C., A.B. and I.C.; writing—review and editing, S.C., I.C. and A.B.; visualization, K.S. and I.C.; project administration, S.C. and A.B.; supervision, T.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Indicators calculated for the analyzed cities.
Table A1. Indicators calculated for the analyzed cities.
No.CityArea
2000
Perimeter 2000Df 2000Area
2006
Perimeter 2006Df 2006Area
2012
Perimeter 2012Df 2012Area
2018
Perimeter 2018Df 2018
1Aberdeen67,263,637102,2371.33869,519,042108,2491.34070,710,440103,8931.33472,249,601104,3081.338
2Alicante50,496,580192,3601.42281,238,955262,1961.44677,068,189202,3821.41279,707,784212,7261.420
3Almeria11,640,21247,3081.14716,880,30051,0131.16220,947,70244,3281.14620,964,78744,3581.146
4Augsburg89,180,879191,8701.32792,362,516193,8661.330107,547,113201,4741.313107,909,640200,8131.313
5Basel157,934,812387,5081.366170,283,024421,1031.382182,923,668401,6531.376184,934,808406,4821.367
6Białystok64,455,729100,1121.32875,273,396156,9221.37680,130,273159,0101.37683,448,907168,1361.378
7Bordeaux308,775,959507,0421.368322,066,733523,4541.373358,002,287524,8261.374361,462,345516,8011.372
8Bournemouth142,992,121217,1211.323147,159,663207,1521.317150,338,702213,8051.324150,513,812214,6681.324
9Brighton and Hove127,602,753252,3041.269133,183,489244,0201.265135,963,558236,7621.263137,381,993243,3051.266
10Brunswick68,843,109138,6181.38172,661,920151,1601.39388,771,563130,1071.36688,957,609130,0871.366
11Burgos20,668,59944,8401.15623,074,20787,2691.27427,158,45583,5901.27335,929,603106,5241.311
12Cambridge37,317,910111,8821.33449,680,391139,5531.34951,974,652142,7521.35353,166,626140,3071.351
13Cheltenham29,106,66040,7981.15230,191,30845,3561.16730,684,75944,6741.16430,918,81245,8001.167
14Chemnitz95,025,243237,1921.400100,236,312227,2381.394134,918,653311,0871.388132,202,434323,7191.393
15Colchester31,736,11573,8221.26036,473,41763,2031.23337,425,12360,9711.22938,124,25261,6651.232
16Crawley45,806,99495,7861.29664,077,035128,9681.32064,920,112128,8771.32065,106,310129,0901.320
17Częstochowa80,208,849251,2421.419109,027,693360,2771.399110,550,809366,0661.402110,727,554366,1741.402
18Derby70,670,12097,4001.32671,285,720108,9301.35171,452,501117,9701.36172,826,329123,7821.370
19Erfurt57,566,592129,6521.36763,353,156151,4061.38567,939,309158,1821.39568,624,471161,2751.390
20Freiburg im
Breisgau
39,406,61673,8351.27341,128,14181,0331.29043,809,28581,5681.29244,000,78482,5421.301
21Geneva145,904,944437,9141.339168,884,570459,8231.351172,782,425470,7741.35417,506,1836473,5321.353
22Gloucester42,890,39671,1971.27445,536,29476,2611.28847,975,86870,3191.27950,161,41269,6181.276
23Osijek24,599,96960,6701.22024,698,58160,6021.22026,061,04065,7731.23526,064,76865,8041.235
24Rijeka47,180,036163,6761.38269,847,365254,6351.38770,362,184256,8541.38970,832,830258,8901.393
25Split46,084,157132,0891.30861,027,336197,3201.28861,440,402199,4551.29061,440,405199,4551.290
26Graz128,847,606330,5781.379159,050,472462,7021.389159,756,599459,5611.388160,518,832461,3991.389
27Karlsruhe86,889,721222,5991.37089,062,547225,8481.373102,202,011211,5311.365104,230,470221,6981.364
28Kiel89,888,134180,7041.38292,140,510180,3461.384100,675,617177,5161.371100,617,931174,3451.373
29Linz98,715,706272,9481.431106,758,145291,9371.424110,305,359297,4181.427110,912,553296,4691.427
30Lubeck78,665,421162,2101.39181,268,488174,3631.40190,610,074172,1191.39693,274,423174,8761.398
31Luton54,020,40076,8521.29754,562,31774,5941.29254,756,66376,4241.29655,096,48776,9311.297
32Magdeburg90,037,58497,1321.29792,257,93297,8281.30195,883,524106,2591.30896,013,646106,2291.308
33Messina36,905,747144,3251.28739,970,198135,6241.23640,630,620137,2921.23841,305,243146,5691.249
34Milton Keynes83,929,32987,0821.29294,139,52487,8311.29494,081,15894,0751.30496,267,59497,6781.310
35Monchengladbach77,042,656230,7231.37377,549,121231,0281.373122,412,255314,5481.345123,498,839313,5071.339
36Munster9,176,26627,8900.9989,544,80831,0761.01410,192,35430,8581.03410,192,35530,8581.034
37Nantes196,300,406497,2021.402210,214,281508,6591.408220,020,040489,2201.401222,787,045484,6511.400
38Newport92,690,567196,9171.337104,052,121226,1911.315107,181,083218,3451.310107,432,449220,1981.311
39Northampton63,236,55983,8861.31272,886,65391,6331.32374,017,54494,5041.32775,220,30194,7751.328
40Norwich58,350,202105,1021.34168,785,504106,6011.34871,163,865113,7341.35572,476,976116,8941.361
41Oviedo30,270,394130,1941.32640,454,565159,6811.35146,695,897159,5491.34246,981,231162,1241.344
42Padua113,177,353480,6541.367122,373,858522,5871.382130,302,472520,8171.380132,809,463523,2601.378
43Radom63,372,848227,4871.39490,990,809300,4461.409102,731,711309,3101.407105,125,141316,7931.410
44Rennes84,397,742147,3011.35297,395,636188,5831.40199,621,628186,8421.399103,629,305180,5971.393
45Rostock80,216,869203,4821.42884,941,850215,4001.44295,138,942214,7341.43795,422,326214,2701.437
46Salzburg60,770,398159,0211.40863,451,964189,8341.41862,272,285197,4661.42362,282,265198,6021.425
47Santander46,175,380206,7761.41258,128,853190,5151.41463,353,448236,6781.43163,459,655236,0631.431
48Stavanger72,917,640188,0311.34174,362,384193,4891.33277,402,931193,7731.33479,872,570190,6961.335
49Strasbourg140,282,522265,1111.341147,162,531264,4281.345150,239,157272,1741.346150,481,183271,2741.345
50Tampere167,864,760501,0691.337173,900,228524,0511.344193,083,347527,9091.348202,683,538508,2791.344
51Taranto37,633,35175,2721.27438,159,63478,1131.27839,231,71974,7041.27640,245,20072,3571.271
52Toruń51,195,714117,3971.36573,447,782181,0931.36380,300,568198,1991.36580,519,719198,5451.365
53Trieste35,666,462119,9881.32945,250,584157,6161.38450,507,992166,3661.37350,950,231168,8361.375
54Trondheim60,497,462118,2991.32661,437,964126,2461.34066,239,648128,9831.33966,917,558129,2001.340
55Uppsala45,758,47996,3041.29946,587,69995,8121.29948,250,041100,9581.31248,504,97399,3961.310
56Västerås54,461,02575,1941.28856,785,93686,0651.30757,440,33485,5261.30677,730,254304,2101.386
57Verona45,782,281137,6321.37076,586,414305,0481.38677,356,561304,8731.38679,615,225337,2881.383
58Vigo63,774,393144,6901.31078,190,764338,2341.38379,615,225337,2881.38360,542,61897,2471.310

References

  1. Fang, C.; Yu, D. Urban Agglomeration: An Evolving Concept of an Emerging Phenomenon. Landsc. Urban Plan. 2017, 162, 126–136. [Google Scholar] [CrossRef]
  2. Samuelsson, K.; Giusti, M.; Peterson, G.D.; Legeby, A.; Brandt, S.A.; Barthel, S. Impact of Environment on People’s Everyday Experiences in Stockholm. Landsc. Urban Plan. 2018, 171, 7–17. [Google Scholar] [CrossRef]
  3. Malamis, S.; Katsou, E.; Inglezakis, V.J.; Kershaw, S.; Venetis, D.; Folini, S. Urban Environment. In Environment and Development: Basic Principles, Human Activities, and Environmental Implications; Elsevier: Amsterdam, The Netherlands, 2016. [Google Scholar]
  4. Cieślak, I.; Biłozor, A.; Szuniewicz, K. The Use of the CORINE Land Cover (CLC) Database for Analyzing Urban Sprawl. Remote Sens. 2020, 12, 282. [Google Scholar] [CrossRef] [Green Version]
  5. Maktav, D.; Erbek, F.S. Analysis of Urban Growth Using Multi-Temporal Satellite Data in Istanbul, Turkey. Int. J. Remote Sens. 2005, 26, 797–810. [Google Scholar] [CrossRef]
  6. Batty, M.; Besussi, E.; Chin, N. Traffic, Urban Growth and Suburban Sprawl; Centre for Advanced Spatial Analysis (UCL): London, UK, 2003; Volume 44. [Google Scholar]
  7. Power, A. Social Exclusion and Urban Sprawl: Is the Rescue of Cities Possible? Reg. Stud. 2001, 35, 731–742. [Google Scholar] [CrossRef]
  8. Wong, A. The Impact of Mobile Phones on the New Urban Poor: Leaving an Urban Footprint? J. Urban Technol. 2008, 15, 25–38. [Google Scholar] [CrossRef]
  9. Viana, C.M.; Oliveira, S.; Oliveira, S.C.; Rocha, J. Land Use/Land Cover Change Detection and Urban Sprawl Analysis. In Spatial Modeling in GIS and R for Earth and Environmental Sciences; Elsevier: Amsterdam, The Netherlands, 2019. [Google Scholar]
  10. Gill, S.E.; Handley, J.F.; Ennos, A.R.; Pauleit, S.; Theuray, N.; Lindley, S.J. Characterising the Urban Environment of UK Cities and Towns: A Template for Landscape Planning. Landsc. Urban Plan. 2008, 87, 210–222. [Google Scholar] [CrossRef]
  11. Batty, M.; Fotheringham, A.S.; Longley, P. Fractal Geometry and Urban Morphology. Fractals in Geography; Prentice Hall: Kent, OH, USA, 1993. [Google Scholar]
  12. Biłozor, A.; Cieślak, I.; Czyża, S. An Analysis of Urbanisation Dynamics with the Use of the Fuzzy Set Theory-A Case Study of the City of Olsztyn. Remote Sens. 2020, 12, 1784. [Google Scholar] [CrossRef]
  13. Wilson, I.D.; Ware, J.M.; Ware, J.A. A Genetic Algorithm Approach to Cartographic Map Generalisation. Comput. Ind. 2003, 52, 291–304. [Google Scholar] [CrossRef]
  14. Wilson, A. Entropy in Urban and Regional Modelling: Retrospect and Prospect. Geogr. Anal. 2010, 42, 364–394. [Google Scholar] [CrossRef]
  15. Ellman, T. Infill: The Cure for Sprawl? Ariz. Issue Anal. 1997, 146, 7–9. [Google Scholar]
  16. Bruegmann, R. Urban Sprawl. In International Encyclopedia of the Social & Behavioral Sciences, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2015; pp. 778–783. ISBN 9780080970875. [Google Scholar]
  17. Forman, R.T.T. Some General Principles of Landscape and Regional Ecology. Landsc. Ecol. 1995, 10, 133–142. [Google Scholar] [CrossRef]
  18. Biłozor, A.; Cieślak, I. Review of Experience in Recent Studies on the Dynamics of Land Urbanisation. Land 2021, 10, 1117. [Google Scholar] [CrossRef]
  19. UN-HABITAT. Urbanization and Development: Emerging Futures. World Cities Report 2016; UN-Habitat: Nairobi, Kenia, 2016. [Google Scholar]
  20. Heimlich, R.E.; Anderson, W.D. Development at the Urban Fringe and Beyond: Impacts on Agriculture and Rural Land; Economic Research Service/USDA: Washington, DC, USA, 2001.
  21. Rana, M.S.; Sarkar, S. Prediction of Urban Expansion by Using Land Cover Change Detection Approach. Heliyon 2021, 7, e08437. [Google Scholar] [CrossRef]
  22. Hens, L.; Thinh, N.A.; Hanh, T.H.; Cuong, N.S.; Lan, T.D.; van Thanh, N.; Le, D.T. Sea-Level Rise and Resilience in Vietnam and the AsiaPacific: A Synthesis. Vietnam. J. Earth Sci. 2018, 40, 127–153. [Google Scholar] [CrossRef] [Green Version]
  23. Chen, J.; Gao, J.; Yuan, F. Growth Type and Functional Trajectories: An Empirical Study of Urban Expansion in Nanjing, China. PLoS ONE 2016, 11, e0148389. [Google Scholar] [CrossRef]
  24. He, Q.; Tan, R.; Gao, Y.; Zhang, M.; Xie, P.; Liu, Y. Modeling Urban Growth Boundary Based on the Evaluation of the Extension Potential: A Case Study of Wuhan City in China. Habitat. Int. 2018, 72, 57–65. [Google Scholar] [CrossRef]
  25. Liang, X.; Liu, X.; Li, X.; Chen, Y.; Tian, Y.; Yao, Y. Delineating Multi-Scenario Urban Growth Boundaries with a CA-Based FLUS Model and Morphological Method. Landsc. Urban Plan. 2018, 177, 47–63. [Google Scholar] [CrossRef]
  26. Zheng, X.; Lv, L. A WOE Method for Urban Growth Boundary Delineation and Its Applications to Land Use Planning. Int. J. Geogr. Inf. Sci. 2015, 30, 691–707. [Google Scholar] [CrossRef]
  27. Tayyebi, A.; Pijanowski, B.C.; Tayyebi, A.H. An Urban Growth Boundary Model Using Neural Networks, GIS and Radial Parameterization: An Application to Tehran, Iran. Landsc. Urban Plan. 2011, 100, 35–44. [Google Scholar] [CrossRef]
  28. Cieślak, I.; Biłozor, A.; Salvati, L. Land as a Basis for Recent Progress in the Study of Urbanization Dynamics. Land 2022, 11, 118. [Google Scholar] [CrossRef]
  29. Arbesman, S. New Ways to Measure Science. Wired, 9 January 2012; 1. [Google Scholar]
  30. Frankhauser, P. Aspects Fractals Des Structures Urbaines. Espace Geogr. 1990, 19, 45–69. [Google Scholar] [CrossRef]
  31. Sonis, M. Book Review: The Dynamics of Cities: Ecological Determinism, Dualism and Chaos. Urban Stud. 1993, 30, 1613–1616. [Google Scholar] [CrossRef]
  32. Rosser, J.B. The Dynamics of Cities: Ecological Determinism, Dualism and Chaos. J. Econ. Behav. Organ. 1994, 24, 115–118. [Google Scholar] [CrossRef]
  33. Batty, M. New Ways of Looking at Cities. Nature 1995, 377, 574. [Google Scholar] [CrossRef]
  34. Tannier, C.; Thomas, I.; Vuidel, G.; Frankhauser, P. A Fractal Approach to Identifying Urban Boundaries. Geogr. Anal. 2011, 43, 211–227. [Google Scholar] [CrossRef]
  35. Feder, J. The Fractal Dimension. In Fractals; Springer: Berlin/Heidelberg, Germany, 1988. [Google Scholar]
  36. de Castro, L.N. Fractal geometry of nature. In Fundamentals of Natural Computing; Chapman and Hall/CRC: Boca Raton, FL, USA, 2020. [Google Scholar]
  37. Mandelbrot, B.B.; Wheeler, J.A. The Fractal Geometry of Nature. Am. J. Phys. 1983, 51, 286–287. [Google Scholar] [CrossRef]
  38. Sparrow, C.; Mandelbrot, B. The Fractal Geometry of Nature. J. R. Stat. Soc. Ser. A 1984, 147, 616. [Google Scholar] [CrossRef]
  39. de Keersmaecker, M.L.; Frankhauser, P.; Thomas, I. Using Fractal Dimensions for Characterizing Intra-Urban Diversity: The Example of Brussels. Geogr. Anal. 2003, 35, 310–328. [Google Scholar] [CrossRef]
  40. Chen, Y.G. Logistic Models of Fractal Dimension Growth of Urban Morphology. Fractals 2018, 26, 1850033. [Google Scholar] [CrossRef] [Green Version]
  41. Sreelekha, M.G.; Krishnamurthy, K.; Anjaneyulu, M.V.L.R. Fractal Assessment of Road Transport System. Eur. Transpor—Trasp. Eur. 2017, 5, 1–13. [Google Scholar]
  42. Thomas, I.; Frankhauser, P.; Biernacki, C. The Morphology of Built-up Landscapes in Wallonia (Belgium): A Classification Using Fractal Indices. Landsc. Urban Plan. 2008, 84, 99–115. [Google Scholar] [CrossRef]
  43. Batty, M. Space, Scale, and Scaling in Entropy Maximizing. Geogr. Anal. 2010, 42, 395–421. [Google Scholar] [CrossRef]
  44. Whitehand, J.W.R.; Batty, M.; Longley, P. Fractal Cities: A Geometry of Form and Function. Geogr. J. 1996, 162, 113. [Google Scholar] [CrossRef]
  45. Purevtseren, M.; Tsegmid, B.; Indra, M.; Sugar, M. The Fractal Geometry of Urban Land Use: The Case of Ulaanbaatar City, Mongolia. Land 2018, 7, 67. [Google Scholar] [CrossRef] [Green Version]
  46. Molinero, C. A Fractal Theory of Urban Growth. Front. Phys. 2022, 10, 418. [Google Scholar] [CrossRef]
  47. Lerner, A.M.; Eakin, H. An Obsolete Dichotomy? Rethinking the Rural-Urban Interface in Terms of Food Security and Production in the Global South. Geogr. J. 2010, 177, 311–320. [Google Scholar] [CrossRef]
  48. Son, N.T.; Chen, C.F.; Chen, C.R.; Thanh, B.X.; Vuong, T.H. Assessment of Urbanization and Urban Heat Islands in Ho Chi Minh City, Vietnam Using Landsat Data. Sustain. Cities Soc. 2017, 30, 150–161. [Google Scholar] [CrossRef]
  49. Mulligan, G.F. Fractal Cities: A Geometry of Form and Function. Cities 1997, 14, 54–55. [Google Scholar] [CrossRef]
  50. Acevedo, W.; Masuoka, P. Time-Series Animation Techniques for Visualizing Urban Growth. Comput. Geosci. 1997, 23, 423–435. [Google Scholar] [CrossRef]
  51. Makse, H.A.; Havlin, S.; Stanley, H.E. Modelling Urban Growth Patterns. Nature 1995, 377, 608–612. [Google Scholar] [CrossRef]
  52. Ustaoglu, E.; Aydinoglu, A.C. Regional Variations of Land-Use Development and Land-Use/Cover Change Dynamics: A Case Study of Turkey. Remote Sens. 2019, 11, 885. [Google Scholar] [CrossRef]
  53. Pelorosso, R.; Leone, A.; Boccia, L. Land Cover and Land Use Change in the Italian Central Apennines: A Comparison of Assessment Methods. Appl. Geogr. 2009, 29, 35–48. [Google Scholar] [CrossRef]
  54. Cieslak, I.; Szuniewicz, K.; Czyża, S. Analysis of the Variation of the Areas under Urbanization Pressure Using Entropy Index. Procedia Eng. 2016, 161, 2001–2005. [Google Scholar] [CrossRef] [Green Version]
  55. Perdigão, V.; Annoni, A. Technical and Methodological Guide for Updating CORINE Land Cover Database; EARSeL: Paris, France, 1997. [Google Scholar]
  56. Biłozor, A.; Czyża, S.; Bajerowski, T. Identification and Location of a Transitional Zone between an Urban and a Rural Area Using Fuzzy Set Theory, CLC, and HRL Data. Sustainability 2019, 11, 7014. [Google Scholar] [CrossRef] [Green Version]
  57. Cieślak, I.; Szuniewicz, K.; Pawlewicz, K.; Czyża, S. Land Use Changes Monitoring with CORINE Land Cover Data. IOP Conf. Ser. Mater. Sci. Eng. 2017, 245, 052049. [Google Scholar] [CrossRef] [Green Version]
  58. United Nations. New Urban Agenda; United Nations: New York, NY, USA, 2017. [Google Scholar]
  59. Valentina, A.; Maria, A.R.; Carmelo, A.; Davide, A.; Ricardo, R.B.; Filipe, B.E.S.; Peter, B.; Paolo, B.; Flavio, B.; Ioris, B.; et al. The Future of Cities; Publications Office of the European Union: Luxembourg, 2019. [Google Scholar]
  60. European Environment Agency; Copernicus Land Monitoring Service. Corine Land Cover; European Environment Agency: Copenhagen, Denmark, 1995.
  61. Vaz, E.; Nijkamp, P.; Painho, M.; Caetano, M. A Multi-Scenario Forecast of Urban Change: A Study on Urban Growth in the Algarve. Landsc. Urban Plan. 2012, 104, 201–211. [Google Scholar] [CrossRef] [Green Version]
  62. Petrişor, A.-I. Using CORINE Data to Look at Deforestation in Romania: Distribution & Possible Consequences. Urban Arhit. Constr. 2015, 6, 83–90. [Google Scholar]
  63. Feranec, J.; Hazeu, G.; Christensen, S.; Jaffrain, G. Corine Land Cover Change Detection in Europe (Case Studies of the Netherlands and Slovakia). Land Use Policy 2007, 24, 234–247. [Google Scholar] [CrossRef]
  64. Bajerowski, T. Metodyka Wyboru Optymalnego Użytkowania Ziemi Na Obszarach Wiejskich; Wydawnictwo ART: Olsztyn, Poland, 1996. [Google Scholar]
  65. Schulz, D.; Yin, H.; Tischbein, B.; Verleysdonk, S.; Adamou, R.; Kumar, N. Land Use Mapping Using Sentinel-1 and Sentinel-2 Time Series in a Heterogeneous Landscape in Niger, Sahel. ISPRS J. Photogramm. Remote Sens. 2021, 178, 97–111. [Google Scholar] [CrossRef]
  66. de Jong, S.M.; Burrough, P.A. A Fractal Approach to the Classification of Mediterranean Vegetation Types in Remotely Sensed Images. Photogramm. Eng. Remote Sens. 1995, 61, 1041–1053. [Google Scholar]
  67. de Cola, L. Fractal Analysis of a Classified Landsat Scene. Photogramm. Eng. Remote Sens. 1989, 55, 601–610. [Google Scholar]
  68. Peitgen, H.-O.; Jürgens, H.; Saupe, D.; Feigenbaum, M.J. Chaos and Fractals: New Frontiers of Science; Springer: Berlin/Heidelberg, Germany, 1992; Volume 7. [Google Scholar]
  69. European Environment Agency. Urban Morphological Zones; European Environment Agency: Copenhagen, Denmark, 2006.
  70. Śleszyński, P.; Gibas, P.; Sudra, P. The Problem of Mismatch between the CORINE Land Cover Data Classification and the Development of Settlement in Poland. Remote Sens. 2020, 12, 2253. [Google Scholar] [CrossRef]
  71. Xie, Z.; Ye, X.; Zheng, Z.; Li, D.; Sun, L.; Li, R.; Benya, S. Modeling Polycentric Urbanization Using Multisource Big Geospatial Data. Remote Sens. 2019, 11, 310. [Google Scholar] [CrossRef] [Green Version]
  72. Krebs., C.J. Ecology: Pearson New International Edition: The Experimental Analysis of Distribution and Abundance; Pearson: London, UK, 2013. [Google Scholar]
  73. Jobe, J.M.; Crow, E.L.; Shimizu, K. Lognormal Distributions: Theory and Applications. Technometrics 1989, 31, 392. [Google Scholar] [CrossRef]
  74. Bretagnolle, A.; Paulus, F.; Pumain, D. Time and Space Scales for Measuring Urban Growth. Cybergeo Eur. J. Geogr. 2002, 2002. [Google Scholar] [CrossRef]
  75. van den Berg, L.; Drewett, R.; Klaassen, L.H.; Rossi, A.; Vijverberg, C.H.T. Urban Europe Vol. 1: A Study of Growth and Decline. Ann. Assoc. Am. Geogr. 1981, 73, 630–632. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions, and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions, or products referred to in the content.
Figure 1. Location of the analyzed cities.
Figure 1. Location of the analyzed cities.
Ijerph 20 02049 g001
Figure 2. Land-use types identified on the basis of level 3 CLC data for selected cities.
Figure 2. Land-use types identified on the basis of level 3 CLC data for selected cities.
Ijerph 20 02049 g002
Figure 3. Changes in the urban boundaries of selected cities between 2000, 2006, 2012, and 2018.
Figure 3. Changes in the urban boundaries of selected cities between 2000, 2006, 2012, and 2018.
Ijerph 20 02049 g003
Figure 5. Classification of land-use types in Alicante: (a) 2000, (b) 2006.
Figure 5. Classification of land-use types in Alicante: (a) 2000, (b) 2006.
Ijerph 20 02049 g005
Figure 6. Classification of land-use types in Lubeck: (a) 2000, (b) 2006.
Figure 6. Classification of land-use types in Lubeck: (a) 2000, (b) 2006.
Ijerph 20 02049 g006
Figure 7. Classification of land-use types in Milton Keynes: (a) 2000, (b) 2006.
Figure 7. Classification of land-use types in Milton Keynes: (a) 2000, (b) 2006.
Ijerph 20 02049 g007
Table 1. Urban expansion matrix for dividing the analyzed cities into eight classes. ↑ represents an increase in perimeter (P), area (A), or fractal dimension (Df) in the analyzed period; ↓ represents a decrease in perimeter, area, or fractal dimension in the analyzed period.
Table 1. Urban expansion matrix for dividing the analyzed cities into eight classes. ↑ represents an increase in perimeter (P), area (A), or fractal dimension (Df) in the analyzed period; ↓ represents a decrease in perimeter, area, or fractal dimension in the analyzed period.
ClassPADf
8
7
6
5
4
3
2
1
Note: The colours used in the table symbolize the different phases of spatial expansions, labelled as in Table 3 and Figure 4.
Table 2. Increase in the area, perimeter, and fractal dimension of the analyzed cities.
Table 2. Increase in the area, perimeter, and fractal dimension of the analyzed cities.
No.CityΔA
2000–2006
ΔP
2000–2006
ΔDf
2000–2006
ΔA
2006–2012
ΔP
2006–2012
ΔDf
2006–2012
ΔA
2012–2018
ΔP
2012–2018
ΔDf
2012–2018
1Aberdeen0.0340.0590.0020.017−0.040−0.0040.0220.0040.003
2Alicante0.6090.3630.017−0.051−0.228−0.0230.0340.0510.006
3Almeria0.4500.0780.0140.241−0.131−0.0140.0010.0010.000
4Augsburg0.0360.0100.0020.1640.039−0.0130.003−0.0030.000
5Basel0.0780.0870.0120.074−0.046−0.0050.0110.012−0.007
6Białystok0.1680.5670.0370.0650.0130.0000.0410.0570.002
7Bordeaux0.0430.0320.0040.1120.0030.0000.010−0.015−0.001
8Bournemouth0.029−0.046−0.0050.0220.0320.0060.0010.0040.000
9Brighton and Hove0.044−0.033−0.0030.021−0.030−0.0020.0100.0280.003
10Brunswick0.0560.0900.0090.222−0.139−0.0190.0020.0000.000
11Burgos0.1160.9460.1020.177−0.042−0.0010.3230.2740.030
12Cambridge0.3310.2470.0110.0460.0230.0030.023−0.017−0.002
13Cheltenham0.0370.1120.0130.016−0.015−0.0030.0080.0250.003
14Chemnitz0.055−0.042−0.0040.3460.369−0.004−0.0200.0410.003
15Colchester0.149−0.144−0.0210.026−0.035−0.0030.0190.0110.002
16Crawley0.3990.3460.0180.013−0.0010.0000.0030.0020.000
17Częstochowa0.3590.434−0.0140.0140.0160.0020.0020.0000.000
18Derby0.0090.1180.0190.0020.0830.0080.0190.0490.006
19Erfurt0.1010.1680.0140.0720.0450.0070.0100.020−0.004
20Freiburg im Breisgau0.0440.0970.0140.0650.0070.0010.0040.0120.007
21Geneva0.1580.0500.0090.0230.0240.0030.0130.006−0.001
22Gloucester0.0620.0710.0110.054−0.078−0.0070.046−0.010−0.002
23Osijek0.004−0.0010.0000.0550.0850.0120.0000.0000.000
24Rijeka0.4800.5560.0040.0070.0090.0010.0070.0080.003
25Split0.3240.494−0.0150.0070.0110.0010.0000.0000.000
26Graz0.2340.4000.0070.004−0.007−0.0010.0050.0040.000
27Karlsruhe0.0250.0150.0030.148−0.063−0.0060.0200.048−0.001
28Kiel0.025−0.0020.0010.093−0.016−0.010−0.001−0.0180.002
29Linz0.0820.070−0.0050.0330.0190.0020.006−0.0030.000
30Lubeck0.0330.0750.0070.115−0.013−0.0040.0290.0160.001
31Luton0.010−0.029−0.0040.0040.0250.0040.0060.0070.000
32Magdeburg0.0250.0070.0030.0390.0860.0050.0010.0000.000
33Messina0.083−0.060−0.0400.0170.0120.0010.0170.0680.009
34Milton Keynes0.1220.0090.002−0.0010.0710.0080.0230.0380.005
35Monchengladbach0.0070.0010.0000.5790.362−0.0200.009−0.003−0.005
36Munster0.0400.1140.0160.068−0.0070.0200.0000.0000.000
37Nantes0.0710.0230.0040.047−0.038−0.0050.013−0.009−0.001
38Newport0.1230.149−0.0160.030−0.035−0.0040.0020.0080.001
39Northampton0.1530.0920.0080.0160.0310.0030.0160.0030.001
40Norwich0.1790.0140.0050.0350.0670.0050.0190.0280.004
41Oviedo0.3360.2260.0190.154−0.001−0.0060.0060.0160.001
42Padua0.0810.0870.0110.065−0.003−0.0010.0190.005−0.002
43Radom0.4360.3210.0110.1290.030−0.0010.0230.0240.002
44Rennes0.1540.2800.0360.023−0.009−0.0010.040−0.033−0.005
45Rostock0.0590.0590.0100.120−0.003−0.0040.003−0.0020.000
46Salzburg0.0440.1940.007−0.0190.0400.0040.0000.0060.001
47Santander0.259−0.0790.0010.0900.2420.0120.002−0.0030.000
48Stavanger0.0200.029−0.0060.0410.0010.0020.032−0.0160.001
49Strasbourg0.049−0.0030.0030.0210.0290.0010.002−0.003−0.001
50Tampere0.0360.0460.0060.1100.0070.0030.050−0.037−0.003
51Taranto0.0140.0380.0030.028−0.044−0.0020.026−0.031−0.004
52Toruń0.4350.543−0.0020.0930.0940.0020.0030.0020.000
53Trieste0.2690.3140.0410.1160.056−0.0080.0090.0150.001
54Trondheim0.0160.0670.0100.0780.0220.0000.0100.0020.001
55Uppsala 0.018−0.005−0.0010.0360.0540.0100.005−0.015−0.001
56Västerås0.0430.1450.0150.012−0.006−0.0010.3532.5570.061
57Verona0.6731.2160.0120.010−0.0010.0000.0290.106−0.002
58Vigo0.2261.3380.0560.018−0.0030.000−0.240−0.712−0.052
Table 3. Classification of cities in the analyzed time intervals.
Table 3. Classification of cities in the analyzed time intervals.
No.City200620122018
1Aberdeen858
2Alicante818
3Almeria858
4Augsburg875
5Basel857
6Białystok888
7Bordeaux885
8Bournemouth588
9Brighton and Hove558
10Brunswick856
11Burgos858
12Cambridge885
13Cheltenham858
14Chemnitz574
15Colchester558
16Crawley858
17Częstochowa781
18Derby888
19Erfurt887
20Freiburg im Breisgau 888
21Geneva887
22Gloucester855
23Osijek588
24Rijeka888
25Split781
26Graz858
27Karlsruhe857
28Kiel652
29Linz785
30Lubeck858
31Luton588
32Magdeburg885
33Messina588
34Milton Keynes848
35Monchengladbach875
36Munster861
37Nantes855
38Newport758
39Northampton888
40Norwich888
41Oviedo858
42Padua857
43Radom878
44Rennes855
45Rostock855
46Salzburg848
47Santander685
48Stavanger786
49Strasbourg685
50Tampere885
51Taranto855
52Toruń781
53Trieste878
54Trondheim878
55Uppsala585
56Västerås858
57Verona867
58Vigo851
Note: The colours used in the table symbolize the different phases of spatial expansions, labelled as in Table 1 and Figure 4.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Czyża, S.; Szuniewicz, K.; Cieślak, I.; Biłozor, A.; Bajerowski, T. An Analysis of the Spatial Development of European Cities Based on Their Geometry and the CORINE Land Cover (CLC) Database. Int. J. Environ. Res. Public Health 2023, 20, 2049. https://doi.org/10.3390/ijerph20032049

AMA Style

Czyża S, Szuniewicz K, Cieślak I, Biłozor A, Bajerowski T. An Analysis of the Spatial Development of European Cities Based on Their Geometry and the CORINE Land Cover (CLC) Database. International Journal of Environmental Research and Public Health. 2023; 20(3):2049. https://doi.org/10.3390/ijerph20032049

Chicago/Turabian Style

Czyża, Szymon, Karol Szuniewicz, Iwona Cieślak, Andrzej Biłozor, and Tomasz Bajerowski. 2023. "An Analysis of the Spatial Development of European Cities Based on Their Geometry and the CORINE Land Cover (CLC) Database" International Journal of Environmental Research and Public Health 20, no. 3: 2049. https://doi.org/10.3390/ijerph20032049

APA Style

Czyża, S., Szuniewicz, K., Cieślak, I., Biłozor, A., & Bajerowski, T. (2023). An Analysis of the Spatial Development of European Cities Based on Their Geometry and the CORINE Land Cover (CLC) Database. International Journal of Environmental Research and Public Health, 20(3), 2049. https://doi.org/10.3390/ijerph20032049

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