Copernicus Land Monitoring Service Research Articles

Integrating Copernicus LMS with ground measurements data for leaf area index and biomass assessment for grasslands in Poland and Norway

ABSTRACT The integration of satellite data from the Copernicus Land Monitoring Service (CLMS) with ground-based measurements represents a pioneering approach to the assessment of leaf area index (LAI) and in situ biomass at grasslands in Poland and Norway. The aim of this study was to develop a method for assessing grass growth conditions and predicting biomass yield based on Sentinel-2 data and CLMS products. LAI values derived from CLMS were compared with in-situ measurements on grassland plots in the two regions of Podlaskie (PL84) and Wielkopolskie (PL41). Small random statistical errors observed in the differences represent a significant opportunity to use Copernicus data to develop a relationship model for biomass prediction. The NDII index calculated using Sentinel2 was considered important for biomass assessment in the humid areas of Norway. This relational biomass prediction model could provide valuable information to farmers, improving their ability to manage grasslands effectively. As a result of the research, relational models were developed has been developed to predict grassland fresh biomass yield with an R 2 accuracy of 0.72 for the first cut, 0.81 for the second cut and 0.91 for the third cut. This allows farmers to effectively manage and monitor grassland throughout the growing season.

Twenty Years of Land Accounts in Europe

Land use and its change impact food security, carbon cycling, biodiversity, and, hence, the condition of ecosystems to mitigate and adapt to climate change, support economic prosperity, and human well-being. To support and guide policy actions between the economy and the environment, harmonized time series datasets, transparent methodologies, and easily interpretable statistics are needed. Therefore, monitoring of the function and condition of lands and their change, along with properly agreed methodologies and freely accessible data, are essential. The Copernicus Land Monitoring Service has produced over 20 years of Corine Land Cover datasets for 39 countries in Europe, which allows continental-wide harmonized and comparable monitoring and accounting of land cover and land use change at a high thematic resolution and in a long time series (2000–2018). With the upcoming 2024 update, the time series will reach a unique product worldwide in terms of time series length, spatial resolution, extent, and thematic detail, enabling policymakers and the scientific community to address the main anthropogenic drivers of land and ecosystem degradation. This paper describes a unified approach for producing continental-wide land accounts that aligns with internationally agreed-upon standards for measuring the environment and its relationship with the economy. Furthermore, the study provides a harmonized time series of geospatial data for deriving land accounts and provides statistics of land cover and land use status and changes for a twenty-year period. All geospatial data and statistics presented in this paper are freely accessible and downloadable to serve other studies.

Conceptual Model for Integrated Meso-Scale Fire Risk Assessment in the Coastal Catchments in Croatia

Various factors influence wildfire probability, including land use/land cover (LULC), fuel types, and their moisture content, meteorological conditions, and terrain characteristics. The Adriatic Sea coastal area in Croatia has a long record of devastating wildfires that have caused severe ecological and economic damages as well as the loss of human lives. Assessing the conditions favorable for wildfires and the possible damages are crucial in fire risk management. Adriatic settlements and ecosystems are highly vulnerable, especially during summer, when the pressure from tourist migration is the highest. However, available fire risk models designed to fit the macro-scale level of assessment cannot provide information detailed enough to meet the decision-making conditions at the local level. This paper describes a model designed to assess wildfire risks at the meso-scale, focusing on environmental and anthropogenic descriptors derived from moderate- to high-resolution remote sensing data (Sentinel-2), Copernicus Land Monitoring Service datasets, and other open sources. Risk indices were integrated using the multi-criteria decision analysis method, the analytic hierarchy process (AHP), in a GIS environment. The model was tested in three coastal catchments, each having recently experienced severe fire events. The approach successfully identified zones at risk and the level of risk, depending on the various environmental and anthropogenic conditions.

Analysis of the temporal and spatial evolution of turbidity in Tonle Sap Lake and its influencing factors

Turbidity is a crucial indicator of water quality. The European Commission's Copernicus Land Monitoring Service Platform provides free turbidity data for large lakes to monitor the global water quality of lakes. However, the data were missing from April 2012 to April 2016, severely limiting long-term analysis. Based on MODIS and turbidity data, Random Forest and XGBoost models are used to invert Tonle Sap Lake's turbidity. Random Forest outperformed the XGBoost model. Based on Random Forest model, missing data were filled in to construct long-term series data of Tonle Sap Lake turbidity (2004–2021). Trend, persistence and correlation analyses were conducted to reveal spatiotemporal characteristics and driving mechanism of turbidity. The results showed that: (1) spatially, the average annual, monthly, and seasonal turbidity was higher in the north but lower in the south, with regions of higher turbidity exhibiting more significant changes; (2) temporally, the annual turbidity mean was 53.99 NTU and showed an increasing trend. Monthly, turbidity values were higher from March to August and lower from September to February, with the highest and lowest recorded in June and November at 110.06 and 5.82 NTU, respectively. Seasonally, turbidity was higher in spring and summer compared to autumn and winter, with mean turbidity values of 84.16, 93.47, 15.33 and 23.21 NTU, respectively; (3) In terms of sustainability, the Hurst exponent for annual turbidity was 0.23, indicating a reverse trend in the near future; (4) Dam construction's impact on turbidity was not significant. Compared with natural factors (permanent wetlands, grasslands, lake surface water temperature, and remote sensing ecological index), human activities (barren, urban and built-up lands, croplands and population density) had a more significant impact on turbidity. Turbidity was highly correlated with croplands (r = 0.76), followed by population density (r = 0.71), and urban and built-up lands (r = 0.69).

Estimation of Spatiotemporal Changes of Land Cover in Llapi River Catchment (Northeastern Part of Kosova)

Aim of the studyIn recent decades, Kosova in general and the Llapi River catchment have experienced landscape transformation which resulted in Land Use and Land Cover (LULC) changes, mainly in flat terrains where most of country’s population is concentrated. Other landscapes have experienced changes too, but on a small scale. The aim of the study was to identify landscape transformation and driving forces behind it.Material and methodsSpatiotemporal changes of Land Use and Land Cover (LULC) were made using GIS techniques. Land Cover data from Copernicus Land Monitoring Service (CLMS) with spatial resolution 100 m for three different years: 2000, 2006 and 2018 were analyzed. DEM with 10 m spatial resolution was used to estimate slope steepness and catchments’ altitude. Population data for settlements inside the Llapi River catchment represented in features were downloaded from the Statistical Agency of Kosova, and those data include population numbers for official censuses between 1948 until 2011, while ArcMap 10.8 environment was used for analysis.Results and conclusionsWith GIS technique in analyzing Land Use and Land Cover data (LULC) and other open source data, we found that agricultural land have decreased in area by 9.53%, while the increase of artificial surface is 160.87%, or expressing by absolute values, total agricultural areas in 2000 were 37,827 ha, while in 2018 were 34,334 ha, meaning 0.3ha per capita, less than European level. In order to protect agricultural land, zoning and land use plans should be implemented.

A novel remote sensing-based approach to determine loss of agricultural soils due to soil sealing — a case study in Germany

Soils provide habitat, regulation and utilization functions. Therefore, Germany aims to reduce soil sealing to 30 ha day-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{-1}$$\\end{document} by 2030 and to eliminate it by 2050. About 55 ha day-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{-1}$$\\end{document} of soil are damaged (average 2018–2021), but detailed information on its soil quality is lacking. This study proposes a new approach using geo-information and remote sensing data to assess agricultural soil loss in Lower Saxony and Brandenburg. Soil quality is assessed based on erosion resistance, runoff regulation, filter functions, yield potential and the Müncheberg Soil Quality Rating from 2006 to 2015. Data from the German Soil Map at a scale of 1:200,000 (BÜK 200), climate, topography, CORINE Land Cover (CLC) and Imperviousness Layer (IMCC), both provided by the Copernicus Land Monitoring Service (CLMS), are used to generate information on soil functions, potentials and agricultural soil loss due to sealing. For the first time, soil losses under arable land are assessed spatially, quantitatively and qualitatively. An estimate of the qualitative loss of agricultural soil in Germany between 2006 and 2015 is obtained by intersecting the soil evaluation results with the quantitative soil loss according to IMCC. Between 2006 and 2015, about 73,300 ha of land were sealed in Germany, affecting about 37,000 ha of agricultural soils. This corresponds to a sealing rate of 11 ha per day for Germany. In Lower Saxony and Brandenburg, agricultural soils were sealed at a rate of 1.9 ha day-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{-1}$$\\end{document} and 0.8 ha day-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{-1}$$\\end{document} respectively, removing these soils from primary land use. In Lower Saxony, 75% of soils with moderate or better biotic yield potential have been removed from primary land use, while in Brandenburg this figure is as high as 88%. Implementing this approach can help decision-makers reassess sealed land and support Germany’s sustainable development strategy.

German and European Ground Motion Service: a Comparison

Since the end of 2022, two ground motion services that cover the complete area of Germany are available as web services: the German Ground Motion Service (Bodenbewegungsdienst Deutschland, BBD) provided by the Federal Institute for Geosciences and Natural Resources (BGR), and the first release of the European Ground Motion Service (EGMS) as part of the Copernicus Land Monitoring Service. Both services are based on InSAR displacement estimations generated from Sentinel‑1 data. It would seem relevant to compare the products of the two services against one another, assess the data coverage they provide, and investigate how well they perform compared to other geodetic techniques. For a study commissioned by the surveying authority of the state of Baden-Württemberg (Landesamt für Geoinformation und Landentwicklung Baden-Württemberg, LGL), BBD and EGMS data from different locations in Baden-Württemberg, Saarland, and North Rhine-Westphalia (NRW) were investigated and validated against levelling and GNSS data. We found that both services provide good data quality. BBD shows slightly better calibration precision than EGMS. The coverage provided by EGMS is better than that of BBD on motorways, federal roads, and train tracks of the Deutsche Bahn. As an example, where both services have difficulties in determining the correct displacements, as they cannot be described well by the displacement models used for processing, we present the test case of the cavern field at Epe (NRW). Finally, we discuss the implications of our findings for the use of the products of BBD and EGMS for monitoring tasks.

Employing Copernicus Land Service and Sentinel-2 Satellite Mission Data to Assess the Spatial Dynamics and Distribution of the Extreme Forest Fires of 2023 in Greece

In 2023, Greece faced its worst wildfire season, with nine major fires causing unprecedented environmental damage of 1470.31 km2. This article uses Copernicus Land Monitoring Service and Sentinel-2 data, employing advanced remote sensing and GIS techniques to analyze spatial dynamics, map burn severity, assess fire extent, and highlight pre-fire tree density and land cover. The study focuses on the catastrophic fire in the Evros region and the damage to the National Forest Park of Dadia–Lefkimmi–Soufli. It also analyzes significant fires in Rhodes, Attica, Thessaly, Evia, Corfu, and Magnesia, emphasizing the compounded challenges posed by terrain, climate, and human factors in those areas. Additionally, the climate data for each affected area were compared with the weather conditions prevailing at the time of the fires. Copernicus Land Cover and Tree Density data are integrated to aid future management, assessment, and restoration. The analysis of maps and fire statistics underscores a notable pattern: areas with higher pre-fire tree density experienced correspondingly higher burn severity. This research underscores the crucial role of such data in assessing wildfire impact. In addition, compared with Copernicus Emergency Management Service, the burned area maps validate the accuracy and reliability of the utilized satellite data. The total burned area was assessed with a high accuracy rate of 96.28%.

Links between Land Cover and In-Water Optical Properties in Four Optically Contrasting Swedish Bays

The optical complexity of coastal waters is mostly caused by the water discharged from land carrying optical components (such as dissolved and particulate matter) into coastal bays and estuaries, and increasing the attenuation of light. This paper aims to investigate the links between in-water optical properties in four Swedish bays (from the northern Baltic proper up to the Bothnian bay) and the land use and land cover (LULC) in the respective catchment of each bay. The optical properties were measured in situ over the last decade by various research and monitoring groups while the LULC in each bay was classified using the Copernicus Land Monitoring Service based on Landsat 8/OLI data. The absorption coefficient of colored dissolve organic matter (CDOM) at 440 nm, aCDOM (440), was significantly correlated to Wetlands which may act as sources of CDOM, while Developed areas (Agricultural and Urban classes) were negatively correlated. The Agriculture class was also negatively related to suspended particulate organic matter (SPOM), whilst Coniferous Forests and Mixed Forests as well as Meadows were positively correlated. SPOM seems thus to mostly originate from Natural classes, possibly due to the release of pollen and other organic matter. Overall, the methods applied here allow for a better understanding of effects of land use and land cover on the bio-optical properties, and thus coastal water quality, on a macroscopic scale.

ADVANCED ANALYSIS TOOLS FOR THE EUROPEAN GROUND MOTION SERVICE DATA

Abstract. This paper is focused on Advanced Differential Interferometric SAR (A-DInSAR). In the first part, the paper describes the European Ground Motion Service (EGMS), a new service of the Copernicus Land Monitoring Service. The EGMS provides consistent, regular, standardized, harmonized and reliable information regarding natural and anthropogenic ground motion over Europe. It is derived from A-DInSAR analysis of Copernicus Sentinel-1 data. The EGMS includes three main products: Basic, i.e. basic deformation maps delivered for individual and consistent frames of image stacks; Calibrated, i.e. deformation map integrated into a reference Global Navigation Satellite System network; and Ortho, which includes the East-West and Up-Down deformation, obtained by combining the ascending and descending results. The second part of the paper describes examples of data analysis. The first one is a clustering of PSI results, which highlights the so-called Active Deformation Areas (ADAs). The second one is an automatic classification of the ADAs. The third analysis concerns the differential deformation maps, which are useful to highlight areas where high deformation gradients occur. Such gradients are often associated to damages in buildings and infrastructures.

Monitoring and Mapping Urban Sprawl Over Heritage Hotspots Using Copernicus Land Monitoring Services: The case of periurban large-scale, wind-powered water extraction mills in Palma (Mallorca)

Among the varied group of human constructions of heritage interest are the old wind-powered water extraction mills that sometimes form large-scale sets in rural areas with high wind potential. On the Mediterranean island of Mallorca (Spain), up to 2,400 windmills dating from the 19th and early 20th centuries have been counted. The areas in which they are located are considered heritage ‘hotspots’, i.e., areas prone to specific problems, such as the progressive abandonment of ethnological heritage resulting from urban sprawl over areas with an agricultural orientation. This article aims to monitor urban sprawl in the municipality of Palma, to quantify and map its impact on a set of windmills located mainly in the plain of Sant Jordi, to the east of the city. The study has been carried out using methodologies and analysis techniques from the Copernicus Land Monitoring Service’s Urban Atlas and Imperviousness Density products. The study shows that areas with agricultural land uses have been progressively transformed into urbanised ones. This transformation has impacted, above all, windmills located in peri-urban areas adjacent to the city. The analysis aims to show the analytical possibilities of Copernicus services and products, and their applicability in the planning and management of peri-urban agro-industrial heritage.

Dealing with the uncertainty of technical changes in the CORINE Land Cover dataset: The Portuguese approach

Land Use and Land Cover (LULC) datasets are widely used across disciplines, with many users demanding more and better information. Understanding the uncertainties and errors associated to the main LULC datasets is a required step to facilitate their correct use, as well as to identify what could be improved in the future production of these products. CORINE Land Cover is probably the most well-known and used LULC dataset in Europe, especially valuable for the rich time-series that it provides. Despite being produced through a change mapping first approach, which tries to avoid technical errors and uncertainties in the temporal analysis of LULC changes, the Copernicus Land Monitoring Service distributes status layers of CORINE (CSL), which are not valid for change analysis because of their associated errors and uncertainties. The CORINE layers of changes (CHA) remove a lot of these issues, but do not meet the needs of many users. In Portugal, the national authority in charge of producing CORINE, the DGT, has implemented a backdating approach to produce consistent CSL layers that allow change analysis with low levels of uncertainty. Throughout this paper, we evaluate the changes that can be analyzed through all available CORINE layers in Portugal: Copernicus CSL layers; the national DGT CSL layers; and CHA layers. To this end, we aim to assess what type of changes can be studied through each type of layer, their associated sources of uncertainty and the relevance and utility of the Portuguese backdating approach to produce a consistent time-series of LULC maps. The results prove how the Portuguese CORINE layers distributed by the Copernicus Land Monitoring Service contain important sources of uncertainty, which however have been removed through the national backdating methodology. This methodology can be therefore exported for the production of CORINE in other European countries.

N°368 – Non-neurologic diseases mimicking stiff person syndrome

In climate modelling, a critical issue is the accurate representation of urban areas in space and variation over time. Urban features may be parameterised through urban canopy parameters (UCPs), including the impervious surface area (ISA). One of the aims of the Consortium for Small-scale Modelling (COSMO) and the CLM-Community is to improve urban numerical weather prediction (NWP) and climate simulations by implementing the urban canopy scheme TERRA_URB in the ICOsahedral Nonhydrostatic (ICON) modelling framework. To harmonise and update the datasets adopted for the UCPs, the Consortium has launched the Priority Project “City Induced Temperature change Through A'dvanced modelling” (PP CITTÀ). Within this framework, the study explores recent and detailed land cover datasets to update ISA and harmonise it with datasets adopted for the UCPs. We identify and compare some relevant datasets at the pan-European level. The results show the benefits obtained by correlating the UCPs from urban Local Climate Zones taken from the ECOCLIMAP-Second Generation with the ISA depending on the imperviousness degree provided by the Copernicus Land Monitoring Service. These results will foster an updating of NWP and climate models with urban surface parameterisations to match urban features while avoiding inconsistencies consistently.

EGMStream, a Desktop App for EGMS Data Downstream

The recent release of European Ground Motion Service (EGMS) products implemented under the responsibility of the Copernicus Land Monitoring Service (CLMS) guarantees free and accessible Europe-wide ground motion data for ground deformation analysis at the local and regional scales. The need for value-adding services and tools for optimal dissemination of radar data from the Copernicus Sentinel-1 satellite mission urges the scientific community to find efficient solutions. A desktop R-based application with a user-friendly interface capable of automatically downloading and transforming EGMS products delivered as large .csv tiles, equivalent to a radar burst into geospatial databases, is presented here. EGMStream is a self-contained desktop app that enables users to systematically store, customize, and convert ground movement data into geospatial databases, burst per burst or for an area of interest directly selectable on the app interface.

Survey of land cover datasets for updating the imperviousness field in urban parameterisation scheme TERRA_URB for climate and weather applications

In climate modelling, a critical issue is the accurate representation of urban areas in space and variation over time. Urban features may be parameterised through urban canopy parameters (UCPs), including the impervious surface area (ISA). One of the aims of the Consortium for Small-scale Modelling (COSMO) and the CLM-Community is to improve urban numerical weather prediction (NWP) and climate simulations by implementing the urban canopy scheme TERRA_URB in the ICOsahedral Nonhydrostatic (ICON) modelling framework. To harmonise and update the datasets adopted for the UCPs, the Consortium has launched the Priority Project “City Induced Temperature change Through A'dvanced modelling” (PP CITTÀ). Within this framework, the study explores recent and detailed land cover datasets to update ISA and harmonise it with datasets adopted for the UCPs. We identify and compare some relevant datasets at the pan-European level. The results show the benefits obtained by correlating the UCPs from urban Local Climate Zones taken from the ECOCLIMAP-Second Generation with the ISA depending on the imperviousness degree provided by the Copernicus Land Monitoring Service. These results will foster an updating of NWP and climate models with urban surface parameterisations to match urban features while avoiding inconsistencies consistently.

DIGITAL ENVIRONMENTALLY SPECIFIC NEEDS OF A SMALL OPEN PIT EXPLOITATION IN A PROTECTED AREA (PART II)

This paper is Part II of specific digital need of a small open pit exploitation. Using several intersections extracted from Copernicus of European Union's Copernicus Land Monitoring Service information, and also the provided GIS Protected Area shapes of the Romanian Ministry of Environment [1], an open pit/quarry project is analyzed concerning the Protected Area Measurements in order to minimize the environmental impact and the investments costs. An ArcGIS Pro License of the Technical University of Cluj Napoca, North University Center of Baia Mare owned in the Faculty of Engineering was used combined with the Commercial License of Rhinoceros 3D.

Influence of urbanisation on the phenology of evergreen coniferous and deciduous broadleaf trees in Madrid (Spain)

Plant phenology plays a crucial role in terrestrial ecosystems and is a key indicator of ecological changes. It is influenced by alterations in meteorological conditions at micro- and mesoscale in cities as a result of urbanisation and the intensification of anthropogenic activities. One example is the urban heat island (UHI), i.e. the increase in temperature in city cores compared to their surroundings. In this study, we used plant phenological data from the High-Resolution Vegetation Phenology and Productivity dataset of the Copernicus Land Monitoring Service to explore the relationship between urbanisation intensity (i.e., proportion of built-up area) and the phenology and productivity of broadleaf deciduous and evergreen coniferous trees in the city of Madrid and its surroundings, one of the largest metropolitan areas in Europe. We found the average start of season date advanced by 0.79–11.79 days in broadleaf deciduous trees in highly urbanised areas, and was delayed by 5.77–9.16 days in evergreen coniferous trees. Although differences were also observed for the end and length of the growing season, these patterns showed great interannual variability. Both broadleaf deciduous and evergreen coniferous trees were most productive in highly urbanised areas. Although previous studies have explored the response of different land-cover types to urbanisation, to our knowledge, none have looked at this relationship at a fine spatial scale in a Mediterranean city with different vegetation functional groups. The analysis of more specific classes of vegetation enables a better understanding of the phenological response of vegetation to urbanisation intensity.

Artificial neural networks for assessing forest fire susceptibility in Türkiye

Wildfires often threaten natural and economic resources and human lives. Wildfire susceptibility assessments have become essential for efficient disaster management and increasing resilience. In this study, we assessed the forest fire susceptibility in Istanbul Province and Thrace Region, Türkiye using a well-known machine learning technique, Artificial Neural Networks (ANN). Benefiting from freely available Earth Observation datasets such as Sentinel-2 images, Tree Cover Density from European Union (EU) European Environment Agency (EEA) Copernicus Land Monitoring Service, Shuttle Radar Topography Mission (SRTM) data, etc., and a forest inventory with ignition locations recorded over a period of eight years, we utilized a total of 16 independent and one dependent variables. The variables can be categorized as anthropogenic, topographic, vegetation, and hydrological factors. A ratio of 1:2 was preferred for the fire/non-fire location samples. The results show that the ANN exhibited high prediction performance with Area Under the Receiver Operating Characteristic Curve (AUC) value and F-1 score of 0.94 and 0.80, respectively. Based on feature importance analyses, we found that a human-related factor, proximity to forest roads, was the most predictive input variable. The ANN model trained with openly available data (i.e., without forest database) also yielded a high F-1 score, but produced maps with fewer details. Our results confirm that data-driven machine learning methods are promising for regional forest fire susceptibility assessments and can be extended further for other regions by deriving similar parameters from freely available Earth Observation datasets.

Land cover data provided by the Copernicus Land Monitoring Service as important sources for thematic cartography

Land cover data provided by the Copernicus Land Monitoring Service as important sources for thematic cartography

"Evolution of Built Surfaces Based on Copernicus High Resolution Layers. The Case of Growth Poles-Based Metropolitan Areas, Romania "

The evolution of built areas is a widely debated subject mostly because their high dynamics and expansion can be related to the economic development. Technically, the availability of European datasets that include the built areas extracted from satellite imagery is an advantage that could be used for further research. The study aims to analyse the spatial and temporal dynamics of the built areas in six metropolitan areas of Romania, in the period 2006-2018. Functional areas were delimited and the new cover of built areas was calculated for each administrative unit. The Imperviousness Classified Change datasets for the periods 2006-2009, 2009-2012, 2012-2015 and 2015-2018 were downloaded in raster format from Copernicus Land Monitoring Service website. The new surfaces, added from one period to another, were outlined in cartographic representations and Google Earth was used for validation. The results showed that, for the entire reference period, the largest increase values of the newly built area were registered in the metropolitan areas of Cluj, followed by Braşov and Timişoara, whilst lower increase values were recorded for Constanţa, Craiova and Iaşi. Also, by analysing the extension of the built areas, we were able to illustrate some spatial development directions within each of the metropolitan areas. The major differences can be explained by the economic development of the examined areas, which included the expansion of both residential and transport infrastructure. Moreover, the pace of development was more intense in the case of the administrative units located nearby the urban centre as opposed to those situated at a considerable distance. This conclusion emphasises the importance of proximity within the metropolitan area studies.