Introduction
Major demographic changes have taken place in Europe since the 19th century. This social evolution has resulted on a relevant urban growth and on significant transformations on land use and land cover (Meirich, 2008). These variations can be approached from several disciplines since they are fostered by complex processes. However, in this specific study these changes have been analysed from a cartographic and geographic perspective.
Soil sealing is one of the main consequences of urban growth. This is a matter of grave concern as it entails ground waterproofing and leads to almost irreversible damages on soil ecosystems (Nizeyimana, et al, 2001; Montanarella, 2007; Scalenghe & Ajmone-Marsan, 2009; EEA, 2012; Siedentop & Fina, 2012). Spatio-temporal dynamics of soil sealing have been regularly studied with the use of 20th century maps, aerial photographs, orthophotographs and satellite images (Ridd, 1995; Zhang, 2001; Wu, 2004; Moeller, 2005; Kampouraki et al., 2006, etc.). However, the use of information predating the mid-20th century is rather unusual since only a few organisations have digitised ancient documents on a large scale and made them accessible to the general public. Instead, this study has included 19th century maps provided by the National Geographic Institute of Spain (IGN). They have been used to analyse patterns of urban growth and the extent to which land has been lost to sealing over the last 150 years in a borough located in the outskirts of Madrid. This is a rather innovative approach that enables this retrospective analysis to be extended in time an additional century to what is usually considered. The effects of soil sealing resulting from urban growth are intensely felt in Spain on coastal regions (Ojeda & Villar, 2006; Valera et al., 2011a) and on major metropolitan areas, like Madrid (García & Pérez, 2007, 2011 and 2016; García et al, 2014).
This study has focused on a specific town called Colmenar Viejo, which is located in the northern outskirts of Madrid (40° 39' 31'' N, 3° 45' 58'' W). This borough has been chosen because it has got both common characteristics to other Spanish towns affected by metropolitan processes, but also distinctive attributes which are only to be observed in this territory, the most outstanding being:
- An environmentally significant and well preserved extension (18,256 ha).
- Proximity to the city of Madrid (30 km) and to several Royal Estates, i.e.
El Pardo y El Soto de Viñuelas, La Granja de San Ildefonso, San Lorenzo de El Escorial and Aranjuez (Fig. 1). Colmenar Viejo is well connected to the main city and to the surrounding towns by both road and rail.
- An urban character right from the beginning of the period under study as Colmenar Viejo has been a medium-size town since the 15th century. Nevertheless, it is still nowadays physically apart from other built-up areas within the metropolitan area of Madrid. This eases its separate study.
- Various attempts to accomplish a scientific Map of Spain based on a geodetic network were carried out in the surroundings of Colmenar Viejo since the 1860's. Therefore, many maps of this area from the mid and late 19th century can be found at the National Geographic Institute of Spain (IGN), which is the only organisation that stores this type of scientific ancient cartography. This large amount of available ancient maps is most uncommon in Spain, where some territories had not been scientifically surveyed even a century later.
One of the main disadvantages of analysing soil sealing is that it is rarely confined to one administrative unit (borough, province, region or country) (Seifert, 2010). As a result, analysis is more about how it has evolved in a specific territory than about comparing its evolution with other areas. Instead, this study has aimed to show the spatio-temporal results of sealing at a municipal level so that it can be compared with other boroughs or territorial units. All information has been integrated in a GIS so it can be measured and updated.
Alarge variety of sources of information has been used to measure and analyse soil sealing. This variability results from the different techniques used for mapping over the last 150 years, i.e. classic topography (18601960 approx), photogrammetry (1960-2005 approx) and remote sensing (2005-2017 approx). This enables the use and integration of various cartographic techniques and its use in combination with GIS programs, what makes this study rather innovative.
In the first place, a careful selection of 365 ancient maps of Colmenar Viejo has been provided by the Topographic Archive at the National Geographic Institute of Spain (IGN). From 1860, when the Spanish geodetic network started being calculated, to 1960, when photogrammetry started being widely used, several mapping projects using classic topographic methods were carried out in Colmenar Viejo. The material gathered for this lengthy period ranges from highly detailed land registry records at a scale of 1:500, to sheets from the National Topographic Map at a scale of 1:50,000. The characteristics of these old maps are summarised in Table 1 and some examples of this ancient cartography are shown in Fig. 2. Should anyone be interested in going in depth in these sources, see Muro et al., 1996; Arístegui, 2014; Urteaga & Nadal, 2001; Moreno, 2008.
Table 1 Ancient maps of Colmenar Viejo created using classic topography [period 1860-1960] (Source: National Geographic Institute of Spain (IGN))
![](/img/revistas/esrj/v21n3//1794-6190-esrj-21-03-00111-gt1.jpg)
![](/img/revistas/esrj/v21n3//1794-6190-esrj-21-03-00111-gf2.jpg)
Figure 2 Examples of ancient maps of Colmenar Viejo created using classic topography [period 1860-1960] (Source: National Geographic Institute of Spain (IGN))
In the second place, the period from 1960 to 2005 has been analysed by means of visually interpreting some aerial photographs. These images have been chosen from various photogrammetric flights, one per decade. They were captured at different scales, varying from 1:18,000 to 1:40,000, on analogue or digital formats, and on a colour or on a black & white scale, as shown on Table 2.
Table 2 Aerial photographs [period 1960-2005] (Source: National Geographic Institute of Spain (IGN))
![](/img/revistas/esrj/v21n3//1794-6190-esrj-21-03-00111-gt2.jpg)
In the third place, the most recent period has been examined by means of remote sensing images from the year 2011. A panchromatic image (with a spatial resolution of 2.5 m) and a multispectral image (with a spatial resolution of 10 m) from Spot 5 satellite (path 69, row 033) were analysed. These images have been supplied by the National Geographic Institute of Spain (IGN) over the National Remote Sensing Plan (PNT).
In the fourth place, some auxiliary sources have also been used. On the hand, an orthophoto from the year 2013 has been selected. It has been provided by the National Geographic Institute of Spain (IGN) over the Aerial Orthophoto National Plan (PNOA). This layer has been used for the first step of the method, which has consisted of georeferencing all maps. On the other hand, an up-to-date cadastral map of the built-up area has been chosen. This map has been provided by the Spanish Cadastral Service. It is in a polygonal vector format (.shp) and represents urban land plots in the year 2013 (Fig. 3). This map has been used as a base layer for the second stage of the method, which has consisted of editing this layer to help calculating soil loss by means of sealing.
![](/img/revistas/esrj/v21n3//1794-6190-esrj-21-03-00111-gf3.jpg)
Figure 3 Base layer: each polygon represents an urban land plot within the built-up area as it stood in the year 2013 (Source: Spanish Cadastral Service (Catastro))
Finally, the different population censuses carried out in Spain since the mid-19th century have been acquired at the National Institute for Statistics of Spain (INE) (Table 3).
Method
The first step of the method has included changing all maps to a raster format and georeferencing them in datum ETRS89, projection UTM, zone 30 T, which is the official reference system in central Spain at present. This process has been carried out with Global Mapper 14 GIS program, using the 2013 orthophoto as a base layer. Between 30 and 50 control points have been assigned to each map. These control points have been located in places that seem to have remained invariable over the last decades. Afterwards a rubber-sheeting conversion with triangulation has been used and rotations, transformations and changes in local and directional scales have been applied. Mean RMS at control points has been 3.4 m. All control points with a RMS above 6 m have not been taken into account. Finally, georeferenced maps have been saved in .ecw format.
The second stage of the method has included part of the method proposed by Valera et al. (2011 b). First, the urban area has been analysed using maps dating from 1861, 1918 and 1952, photos from 1968, 1979, 1983, 1991 and 2002, and satellite images from 2011. Second, the 2013 cadastral base layer, downloaded on Global Mapper 14 GIS application, has been progressively edited in inverse chronological order by means of visually interpreting maps, photos and satellite images for each date. As the process has gone back in time, plot erasing has always been needed (it has never been necessary to add them) since the built-up area has always shown an increase with time, never a decrease.
Due to the differing nature of mapping techniques, certain decisions had to be taken. On the one hand, spectral analysis using Erdas Imagine 2011 remote sensing program has been used to edit the 2013 cadastral base layer with respect to the 2011 satellite image. On the other hand, to estimate the urban area in 1952, for which only rural maps were available, it has been necessary to erase those plots that appeared on the 1968 edited base layer, but not on the 1952 map. However, the territory comprising the urban area could not be modified because no information from 1952 was available. Therefore they have been left the same as they were in 1968. Also the variability in the different urban units represented on ancient maps had to be overcome. For instance, when it came to adjusting the base layer to the urban area in 1918, maps referred to city blocks and not to land plots. In other words, the available information was less detailed on the 1918 map than on the base layer, since each block includes several plots. Less spatial information meant that blocks had to be edited as a whole as it was not possible to discriminate information related to each plot.
The third level of the method has included analysing patterns and determining urban growth by means of using ArcGIS10.1 GIS program. In this regard, Pinho & Oliveira, 2009 and Tombolini et al., 2015 have been consulted. Two aspects have been assessed:
- Two-dimensional growth considering only the x and y axis. These calculations help examining the extent of soil sealing in quantitative terms. The sealed soil area (in hectares), the sealing increase between consecutive dates (in percentages), the mean annual sealing growth (in hectares per year) and the sealed soil among the total municipal territory (in percentages) have been calculated for each date.
- Three-dimensional growth considering also the z axis. These figures help analysing the patterns followed by urban extension from a qualitative perspective. Since building heights were not represented on maps and there was therefore a lack of information, this analysis has been based on auxiliary information related directly to vertical growth, i.e. population and population density. Population of the town for each date from 1861 to 2011 has been obtained from a linear function between the two population censuses closest in time. Population density for each date has been determined by dividing the estimated population by the urban area. Three-dimensional urban growth has been calculated by means of the following parameters: population (in inhabitants), population increase between consecutive dates (in percentages), mean annual population growth (in inhabitants per year), population density (in inhabitants per sq km) and population density mean annual increase (in percentages).
The fourth step of the method has consisted of examining the territory in situ and verifying the obtained figures and calculations.
Results
The main calculations related to soil sealing and the characteristics of urban growth have been summarised in Table 4 and Fig. 4.
With respect to the two-dimensional growth of the urban area, which refers to the extent of the process from a quantitative perspective, the following should be noted:
- The old town was small in relation to the present-day built-up area. In 1861 it comprised 25.3 hectares, i.e. only 0.14% of the municipal territory.
- From 1861 to 1918 urban growth was subdued, with a total increase in the built-up area of 40%, which implies a mean annual growth of less than 0.20 hectares. This means that the loss of soil due to sealing in Colmenar Viejo affected hardly 0.19% of the municipal territory by 1918.
- From 1952 to 1991 a considerable urban extension took place and thus the built-up area of the town easily tripled. Mean annual growth showed ever greater increases, going from 3 hectares/year at the initial time to nearly 6 hectares/year at the final stage. At the beginning of the 1990's, 1.31% of the municipal territory was covered by built-up area.
- The greatest urban growth has taken place over the last two decades. From 1991 to 2011 the built-up area has doubled, with a mean annual growth of between 9 and 16 hectares. At the end of the period, the percentage of sealed soil has already reached 2.64%.
Therefore, the two-dimensional growth of the built-up area throughout the period under study has been both progressive, i.e. all periods have recorded positive increases, and exponential, i.e. the mean annual growth has increased over an ever larger area. The urban area has doubled in size during the period of fastest growth which comprises only twenty years, i.e. from 1991 to 2011. As a result, 2.64% of the municipal territory is now sealed as a consequence of urban processes. This percentage is lower than in other European cities, where soil sealed from 1990 to 2006 within a range of 6% to 9% (EEA, 2012), and it is also lower than in other boroughs within the Region of Madrid, where these figures vary from 8.5% to 13.5% (García et al., 2014). These discrete figures are probably related to the proximity of Colmenar Viejo to some Royal Estates and the influence of their environmental protection. It could therefore be concluded that the data recorded in Colmenar Viejo should not give rise to much concern. However, it should also be borne in mind that since the borough under study covers the second largest municipal territory within the Region of Madrid, small quantities in relative terms are highly significant in absolute terms.
With respect to the three-dimensional growth, which shows the patterns followed by urban processes from a qualitative perspective, the main features are:
- Population growth has been steady throughout the whole period under study. The greatest mean annual growths were recorded in the 1970's, with around 700 inhabitants per year, and during the final stage, with over 1,000 inhabitants per year.
- Population density has varied constantly. It registered a plus sign over some periods, which shows a rise on population density, and a minus sign during some others, which denotes a drop on population density.
- Variations in population density have remained modest, with mean annual changes within a range of ±1.5%.
If only population statistics were to be taken into account, it could therefore be presumed that while population has risen steadily in numerical terms, it has not followed the same patterns in terms of population density. This infers that there have been some periods in which the built-up area has spread loosely and population density has dropped, and others in which population density has risen as a result of high building construction. However, if, in addition, these data are broken down into periods corresponding to available maps, and statistics are complemented with the ground truth, it can be concluded that the dominant processes are consistent with the general model of urban growth followed by metropolitan outskirts of major Spanish cities. These can be summarised as follows:
- Initial period from 1861 to 1918: modest yet steady extension, following the same patterns observed historically up to that time, i.e. narrow streets and one-storey buildings
- Periods of vertical growth:
1968-1979: tower blocks with very few gardens and leisure areas
1991-2002: residential estates of 3-4 storey buildings in landscaped surroundings
1918-1952: similar houses to those in the old town but including bigger courtyards and vegetable patches
1952-1968: detached houses for holiday home use
1979-1991: detached and semi-detached houses with ample gardens
2002-2011: vast urban extension and exorbitant soil sealing; extensive residential fabric with a predominance of houses over blocks
Conclusions
This study has analysed the spatio-temporal dynamics and has enabled the measuring of the soil sealing process due to urban growth over the last 150 years in a borough located in the northern outskirts of Madrid, i.e. Colmenar Viejo. This town is one of the boroughs of highest environmental value around the capital city of Spain and has got both common characteristics to other Spanish towns influenced by metropolitan processes and specific attributes that are only to be observed in this territory.
The sources and the method used in this study are rather innovative. On the one hand, most soil sealing analysis only consider maps since the mid-20th century. Instead, this study has enabled this retrospective study to go back in time an additional century with the use of 19th century maps obtained from the National Geographic Institute of Spain (IGN). On the other hand, integrating in a GIS program a large volume of maps made by different cartographic techniques, i.e. classic topography, photogrammetry and remote sensing, has guaranteed the capture of a detailed and reliable urban cartography since the mid-19th century and has enabled the use and integration of various mapping techniques. Furthermore, this study has aimed to show the spatio-temporal results of sealing at a municipal level so that it can be measured, updated and compared with other boroughs or territorial units.
Calculating the extension of urban growth from a quantitative perspective has been fulfilled by georeferencing ancient maps, digital maps, aerial photographs and satellite images, integrating them into a GIS program and analysing some parameters regards the two-dimensional urban extension. Observing the patterns followed by urban growth from a qualitative point of view has been achieved by using population censuses accomplished since the mid-19th century and using an indicator system to analyse the three-dimensional urban growth.
Part of Valera's et al. (2011b) method for editing and adjusting maps to the urban area according to different dates has been used. On the one hand, always using the same layer as the base map avoids slivers. On the other hand, using a vector layer for a plot detail provides accurate and precise information for quantifying and analysing urban growth.
Results show that, in terms of quantity, urban growth has been both progressive, i.e. all periods have recorded positive increases, and exponential, i.e. the mean annual growth has increased over an ever larger area. Along the whole period, urban extension has increased from 25.3 hectares to 482.5 hectares, i.e. from 0.14% to 2.64% of the municipal territory, what means multiplying 19-fold its size. The trend followed by urban growth has been therefore similar to other Spanish towns influenced by metropolitan processes. However, the intensity of this process has remained slightly under the European and Madrilian average, which is probably related to its proximity to some Royal Estates, what has conferred this territory a stronger environmental protection.
In terms of quality, the patterns observed on urban extension comply with the general urban growth model recorded during the 20th and 21st centuries in the outskirts of major Spanish cities.