Urban Geomorphology Methods and Applications as a Guideline for Understanding the City Environment

: Cities all over the world have developed on different geological-geomorphological substrates. Different kinds of human activities have operated for millennia as geomorphic agents, generating numerous and various erosion landforms and huge anthropogenic deposits. Considering the increasing demand for land and the expansion of the built-up areas involving and disturbing any kind of natural system inside and surrounding the actual urban areas, it is not negligible how important the dynamics of the urban environment and its physical evolution are. In this context, this manuscript addresses insights into eight case studies of urban geomorphological analyses of cities in Italy, Greece, and Brazil. The studies are based on surveying and mapping geomorphological processes and landforms in urban areas, supporting both geo-hazard assessment, historical evolution, and paleomorphologies, as well as disseminating knowledge of urban geoheritage and educating about the anthropogenic impact on urban sustainability. We hypothesize that urban geomor-phological analysis of several case studies addresses the physical environment of modern cities in a multi-temporal, multidisciplinary, and critical way concerning global changes. Thus, this study aims to illustrate and propose a novel approach to urban geomorphological investigation as a model for the understanding and planning of the physical urban environment on a European and global scale.


Introduction
"The urban geomorphology is the study of man as a physical process of change, whereby he metamorphoses a more natural terrain to an anthropogene cityscape" [1].In the framework of anthropogenic geomorphology [2], urban geomorphology focuses on the changes to nature's physical systems caused by human activities in cities (urban residential, economic, and traffic functions) [3,4].Considering the contemporary world challenge of coping with the high rate of change of the urban population (i.e., population density), the consequent urbanization growth rate [5][6][7], and their management, the discipline is of increasing interest.Growing urbanization represents a threat to the natural environment due to land consumption and to the urban environment itself, which is very sensitive from an ecological point of view and to the natural risks deriving from global changes.
An urban environment is defined as an ecosystem of an urban area in which the urban residents interact with biotic and abiotic factors [8,9].The residents' interaction with the abiotic factors of the urban environment implies their knowledge and, in this context, this paper deals with urban physical environment analysis aimed at knowledge of physical processes, and geomorphological ones in particular, acting in cities. Cities all over the world have different geological-geomorphological substrates, and urban soils, and several kinds of human activities have operated since different times.Humans, as a geomorphic agent, have always generated anthropogenic erosion and accumulation, moving "tremendous amounts of soil and rock" since the Paleolithic period [10].The urbanization of historical European cities is characterized by the century-old stratification of urban expansion phases [11], from their ancient foundation, to the Middle Ages and major development during the Industrial Revolution, until the intense and progressive growth that started after the Second World War and is still in action [12].
In the event that by 2050, 70% of the world's population will live in cities, the increasing demand for land will favor uncontrolled expansion of built-up areas, involving and disturbing any kind of natural system surrounding the actual urban areas.This condition favors the increase of natural risk in present-day densely populated urban areas and will affect cities globally in the future [13,14].The current trend towards urbanization affects Africa and Asia in particular, where there is a city model greatly different from what was conceived in Europe.The African and Asian megalopolises and megacities arising are a little-known but very interesting field for understanding the differences in terms of hazards and resources for such large areas [15][16][17][18][19].At the same time, the European cities are growing too, with evident problems which are challenging the idea of sustainable development.In this changing world, urban geomorphology can aid better understanding of the landscape's natural, historical, and anthropogenic evolution, the settlement's choices in relation to natural morphologies, and the effects of city development on natural geomorphological processes.
Urban geomorphological studies in European and Italian cities have strongly contributed to the building of a methodological approach.The world's scientific literature is enriched by those studies, as shown by the number of Italian papers in the scientific literature collected for this article's background (see Supplemental Material), which is summarized in Table 1.
In this context, the Working Group on Urban Geomorphology of the Italian Association of Physical Geography and Geomorphology (AIGeo) strongly contributed to the discipline with its activity, considering the importance of urban geomorphological analysis in Europe and in Italy especially, where the cities were founded in historical times, expanded in the Middle Ages and progressively entered a period of uncontrolled urban sprawl in the 20th century.Several examples of urban geomorphological analysis are available for Italian cities all along the Italian peninsula.Moreover, the AIGeo WG, collaborates with international colleagues to share and improve methodologies applied to foreign cities.The eight case studies presented in this paper (Figure 1) address the following: • surveying and mapping geomorphological processes and landforms in urban areas • geomorphological analysis supporting geo-hazard assessment • analyzing landscape historical evolution and paleomorphologies • disseminating knowledge of urban geoheritage and educating about anthropogenic impacts on urban sustainability.The hypothesis at the core of this paper is that European and Italian cities' urban geomorphological analyses address the physical environment of modern cities in a multitemporal, multidisciplinary and critical manner with respect to global changes.Thus, the objective of this study is to illustrate the Italian approach to urban geomorphological investigation and propose it and its results as a model and example for physical urban environment knowledge and planning on the European and global scales.

Description of the Study Areas
The selected case studies shown in the location map (Figure 1) are six Italian cities selected all along the Italian peninsula, three along the western coast or a few km from it (from north to south: Genoa, Rome, and Pozzuoli) and three along the central-southern Apennine axis (Perugia, Benevento, and Potenza); furthermore, two other non-Italian cases were considered: a city in the Greek mediterranean area-Patras-and a city located in southern Brazil-São João da Barra.
Table 2 summarizes the main features of the cities: the city typology, in terms of size [20] and degree of urbanization [21][22][23][24]; the geographical context, in terms of landscape; the main geomorphological processes acting; the study topic; and the methodological approach used.
A brief overview of the cities, including their geological and geomorphological frameworks, precedes the discussion of each case study in the specific paragraph.The description will facilitate understanding of the structure of each city to compare it with other cities of interest and evaluate the investigation method applicability.
Table 2. Main features of the 8 case studies.In the column typology, each city's size [20] and degree of urbanization [21][22][23][24] are described.The column headed 'geographical context' describes the landscape character.The column headed 'main morphogenetic processes' is compiled considering the cities geomorphological analysis.The study topics are introduced and the methodological approaches summarized in the last two columns.

Methods and Materials for Urban Geomorphology Analysis and Applications
Geomorphological surveying is the basis of other kinds of analysis and relative applications for each case study; it consists of the following: a. scientific literature and technical materials collection; b. a multitemporal and multidisciplinary approach; c. field surveying.The multitemporal and multidisciplinary approach is based on aerial photograph interpretation, historical cartographic documents analysis and comparison to actual topography, nature of bedrock and deposits analysis through geognostic surveys, qualitative and quantitative assessment, reference to historical documents and archaeological reports, and ancient iconographic materials such as maps, paintings, or photographs [25][26][27].Field surveys are fundamental for observations and direct data checking.The analysis method exploits the principles of data overlay and spatial analysis in the GIS environment.
In the urban environment, what emerges is the importance of the multidisciplinary approach, consisting of cultural (stricto sensu), archaeological, geoarchaeological, and geohistorical information about the areas.In fact, in the urban context, landscape evolution is strongly influenced by anthropogenic action, which buries, transforms, deletes, or hides the geological and geomorphological features.Multidisciplinary data integration allows the recognition and reconstruction of geological and geomorphological characteristics; in particular, anthropogenic erosion and accumulation.
The above-described approach required the availability of the multitemporal and multidisciplinary materials listed below:

•
Base maps: Recent Regional Technical Maps (C.T.R.), at detailed scales, from the 1990s to the present (several Italian case studies).See references therein [26].The different case studies widen the urban environmental analysis, focusing on the topics listed in the introduction and in Table 2.They employed the methods summarized in the last column in Table 2 and described in detail in the case studies paragraph.
To discuss the proposed methodological approaches and applications as a guideline for knowledge of the physical urban environment, a brief review is implemented in Table 1 aimed at exploring how the scientific literature has addressed the following topics: 1. a-surveying and mapping urban geomorphology 2. b-geomorphology and urban hazards 3. c-historical evolution and paleomorphologies 4. d-disseminating education and knowledge of urban geoheritage.
The core collections of Scopus and Google Scholar are used as the search databases for papers, based on the following keywords: Anthropogenic Landforms; Urban Geomorphology Method; Urban Environment; Urban Geomorphology; Anthropogenic Geomorphology; and Anthropocene Geomorphology.

Surveying and Mapping Urban Geomorphological Processes and Landforms
Geomorphological maps are widely used in land planning and territorial management, environmental impact assessment, visualization, and communication of scientific information [27,29].They are the starting point for many applications and for realizing thematic, derived, and applied maps (e.g., hazard maps).Despite researchers from several countries' efforts to set up a standardized mapping procedure and legend, currently, a great variety of geomorphological legends are still in use, which differ in their contents and symbols.
In Italy, the Italian National Institute for Environmental Protection and Research (IS-PRA) proposed national common and standard geomorphological mapping criteria at 1:50,000 in scale [30].Following these criteria, landforms and associated deposits are depicted by using symbols of different colors, and grouped according to genetic criteria, whereas their state of activity is marked by means of color shades.The lithological characteristics of bedrock are emphasized using solid colors by grouping lithotypes according to their similar mechanical behavior.Recently the AIGeo, in collaboration with ISPRA, reviewed the official Italian geomorphological legend, revising and updating graphical symbols due to the great progress made in geomorphology regarding new knowledge on geomorphological processes and landforms, anthropogenic landforms in particular [31].
Using this framework, thanks to the case studies described in this section, the mapping and understanding of anthropogenic landforms was achieved by proposing innovative strategies for the recognition of urban environments, erased natural morphology, and the evolution of the landscape under the influence of human activities.

Rome
The Rome urban area analysis has been one of the first case studies of urban geomorphological surveying and mapping at the Italian and European levels.The study area is 600 km 2 wide in the municipality of Rome, partially covering the territory of the Italian capital city.It is in the central part of the Italian peninsula, on the western side of the Lazio region.Rome is a large metropolis, with about three million inhabitants living in the municipal area.The city arose just south of the confluence where the Aniene River meets the Tiber about 3000 years ago; thus, the area has been densely urbanized for millennia.Until the 19th century, the urbanized area was limited to the historical city walls (Aurelian and Leonine) which delimit the historical center.Between 1924 and 1936, the city's expansion reached an area that is now enclosed in a railway ring surrounding the city.After 1949, and during the 1960s in particular, the city's expansion followed a south-southeast direction and in recent years it has also involved the northeastern area and the western part.The relief is divided between the minimum altitude of 10-20 m in the floodplains of the Tiber and Aniene and the maximum altitude of 100-200 m on the top of the volcanic plateau.A Plio-Pleistocene marine sedimentary basement characterizes the area [32].On this basement, volcanic, tectonic, and glacio-eustatic processes have interplayed since Middle Pleistocene.The volcanic plateau was generated by the Sabatini and Alban Volcanos [33], and is made of effusive rocks which alternate with continental sediments.Late Pleistocene fluvial erosion shaped the volcanic plateau, resulting in the historically hilly, fluvial landscape of Rome, which has been deeply modified by 3 millennia of urbanization.
An urban geomorphological survey method was proposed in 2016 and improved in 2020 [26,34] based on multitemporal analysis, multidisciplinary investigations, and innovative strategies for anthropogenic landform mapping.The method's strength is the multidisciplinary approach: it integrates historical, archaeological, iconographical, and artistic data with multitemporal topographical, geological, geomorphological, and geognostic data.The GIS environment allows for the management of data overlay and the morphometric processing of historical data.In particular, the anthropogenic erosion and accumulation survey is supported by the processing of the DEM of Differences (DoD), a map algebra raster difference between the present Digital Elevation Model (DEM) and historical ones obtained from early maps [35,36].
Due to 3 millennia of human activities modifying the city landscape, the case study deeply contributed to anthropogenic landform classification and new cartographic solutions for their representation.Historically, European cities like Rome have experienced drainage network modifications to obtain wide surfaces for building activities and in order to manage hydrogeological risks.The valleys in Rome have been partially or completely filled by anthropogenic deposits.So, natural modified landforms are present, such as valleys, whose shapes are still recognizable, but anthropogenic deposits modified the section (Figure 2a).There are also anthropogenic landforms, such as valleys which are no longer readable in the landscape because they are completely filled (Figure 2b).Artificial reliefs are also present because of anthropogenic accumulation, such as man-made hills and tabular reliefs due to retaining walls (Figure 2c,d).An example of anthropogenic erosion is excavation associated with mining or obtaining space for construction activities (Figure 2e).Lastly, terraced slopes are mixed examples of both anthropogenic erosion and accumulation (Figure 2f).As shown in the map legend (Figure 2), natural landforms modified by man, but still recognizable, were grouped into the original morphogenetic process category (in the example, fluvial landforms are depicted in green), adding a special column named 'modified', whereas anthropogenic landforms were included in a specific category, depicted with black symbols, when they completely erased the natural morphology.

Genoa
Genoa is a metropolis of 583,000 inhabitants spread over an area of 240 km 2 .The urban area is developed along a narrow coastal strip of the Liguria region for almost 30 km, mainly on the terminal section of the floodplains that cross the shoreline.The main drainage axes include those of the Leira, Cerusa, Varenna, Polcevera, Bisagno, Sturla, and Nervi streams.The current municipal area, however, is recent and dates back to 1926 following the unification of 19 municipalities that were added to the six incorporated in 1874.The historic Genoese part, known since the Middle Ages as the Maritime Republic and remaining virtually unchanged until the first half of the 19th century, lies in the morphological amphitheater between the Lanterna promontory to the west and the Carignano hill to the east and culminates with the summit of Fort Sperone at 500 m a.s.l.[37].This area of about 8.5 km 2 is characterized by seven small catchments, largely modified by anthropogenic activities through culverts and urbanization, ranging in size from 0.5 to 2.5 km 2 , with slopes of up to 50 percent and drainage density between 2.5 and 4 km −1 [38].
The geological profile is rather simple because it is almost entirely characterized by a flysch of marly limestones with shales, siltstones and marl interlayers [39].Only at the historic center, in the areas featuring low steepness, are Pliocene clays mostly inserted in graben-type tectonic structures.
To recognize the geomorphology in the urban landscape of the morphological amphitheater on which the historic center of Genoa lies, a thematic map and bibliography were first studied.Then, by adopting a GIS platform, the stratigraphic data derived from the information layers of drill surveys were inserted, which allowed for an assessment of the anthropogenic filling thickness.Multi-temporal cartographic comparison aimed at identifying landforms, especially before the 19th century urbanization, has been crucial in identifying both excavations and fillings, changes to the hydrographic network, and sea embankments [40].Field work and observations confirmed and improved the survey of urban landforms previously identified by remote surveying (Figure 3).Nearly 70 percent of Genoa's historic area is characterized by anthropogenic landforms, of which the most visible involve the slopes and coastline, as well as changes in the hydrographic network.Approximately 17 km of the riverbed, out of a total length of about 24 km, is culverted, with changes in both the planimetric and longitudinal profile.The end section of the watercourses has been culverted since the Middle Ages, and the modifications include artificial diverting of the watercourse, as in the case of the San Gerolamo in the Carbonara stream.Many buildings of great cultural value are built over the culverts, such as the Palazzo San Giorgio and the Commenda di Pre complex.Anthropogenic landforms on slopes include fills, especially at valley bottoms and for road construction; many historic buildings, such as the Albergo dei Poveri, are built on top of filling materials.They also include dams for water reservoirs, later filled in for development purposes.Among the slope cuts, related to road opening and quarrying, is described the case of the Lanterna Promontory, once the subject of stone extraction activity and dismantled in the 1930s for city roads.The remains of this promontory, represented by an isolated rock outcrop, are still observable within the port area, where the "Lanterna of Genova" stands.The construction of coastal embankments for land reclamation represents one of the major morphological modifications which interested this area following the development of the Genoa harbor.The sea embankments started in the Middle Ages with the growing power of the Republic of Genoa.The extension of the coastline became significant in the 19th century; meanwhile, the port reached its current dimensions after the building of basins between the end of the same century and the beginning of the 20th.

Geomorphological Analysis Supporting Geo-Hazards Assessment
Geomorphological analysis through mapping and field monitoring can be an important tool for diagnosing and supporting the management of hazards and risk events.In this sense, using different geotechnologies in the field and the office allows us to obtain significant results to guide public policies in the context of anthropogenic and natural relations in urban areas.Therefore, this section includes case studies in four cities in three different countries, whose urbanization processes are completely different, but who have similar environmental problems and challenges.The cities of Pozzuoli, Patras, São João da Barra, and Potenza present problems that can be supported through analysis, monitoring, and geomorphological support with different methodological approaches.In these cases, procedures such as bathymetric surveying, sediment collection and analysis, multiscale mapping with orbital and non-orbital images to understand coastline displacements due to coastal erosion, systemic analysis of the processes that resulted in the current stand out natural and anthropogenic forms of the landscape, geological models of subsidence, and geostatistics should be employed.

Pozzuoli
The city of Pozzuoli faces the homonymous large bay, with a mean depth of about -60 m and a maximum depth of −110 m [41], and into the wider Gulf of Naples (Figure 4, upper).Situated on the Margin of the Tyrrhenian Sea in the Campania region, it originated during a Pliocene extensional tectonic phase [42].Some fifty craters have been recognized both inland and offshore.Two main active volcanic districts, the Phlegraean Fields to the northwest and Mt.Somma-Vesuvius to the southeast, developed along the coastland.Pozzuoli represents the emerged-submerged sector of the Phlegraean Fields as characterized by the marine terrace of La Starza between 30 and 80 m a.s.l.[43], bordered by the pyroclastic monogenic volcanoes of Monte Nuovo (1538 CE), Solfatara (4 ka), and Averno (5.2-4.2 ka), and a seabed central plain, with epiclastic sedimentation, delimited to the south by several submerged volcanic banks.The tectonics of the Phlegraean Fields are very dynamic [44]: historical records of bradyseismic vertical movements date back to Greek colonization over 2000 years BP, and the area mainly features a general downlift of the bay and its urban structures that today are partially submerged down to 4-6 m in depth (Figure 4, lower).In 1538, a new mountain, Monte Nuovo, beside the town of Pozzuoli, was formed by a pyroclastic eruption within eight days of an uplift phase of 4.5 m.In 1970In -1972In and 1982In -1984, , the ground uplifted at a rate of 3 mm per day due to positive bradyseismic movement, with a maximum uplift of about 3.8 m followed by a gradual downlift [45,46].From 2005 until today, the uplift started again with a total rise of over 1.15 m, an acceleration since January 2023, and many seismic shocks with 1 < Md <4.2.Gas emissions (fumaroles), both on land (Solfatara Volcano) and underwater (Pozzuoli Bay), are also associated with volcanic activity.The presence of submarine hydrothermal springs located along the northwest-southeast axis in the eastern sector of the gulf [47] are related to the secondary volcanic activity of the Phlegraean Fields.
Since 2015 [48,49], morphological and stratigraphic surveys along the emerged and submerged urbanized area, based on the interpretation of 45 boreholes carried out since the 80s, bathymetric and underwater surveys, and sediment sample analyses, highlighted the main geomorphic elements of Pozzuoli (Figure 4, upper).The volcanic and marine sediment distributions in the town are strictly correlated with climatic and tectonic events.The urban geomorphology of Pozzuoli is characterized by vertical soil motions due to bradyseismic phenomena and volcano-tectonics [50], is rich in displaced Roman-age archaeological ruins [51,52], and hit by littoral erosion and flash floods.The underwater area facing the city shows terraced surfaces down to −40 m of depth, a receding edge of the narrow continental shelf, and the two relic-submerged valleys of Epitaffio (west) and Bagnoli (east) [53], starting from shallow water and with alluvial-detritic fans at their end in the gulf depocenter at −90 m of depth.Depth is in meters a.s.l. and the geographic coordinate system is WGS84.Lower: Geothematic map of the lowering curves of the Phlegraean Fields coastland related to vertical movements between the Greco-Roman period and the present (after [51]).Subsidence is in mm/year and the geographic coordinate system is WGS84; DTM Lidar from MATTM-Environmental Remote Sensing Plan (PT-A).

Patras
The city of Patras is the third-largest urban area in Greece, home to about 171,000 people [54] in the northwest Peloponnese.It lies in a belt of low land between the Panahaikon mountains and the Patraikos Gulf.The wider area of the city of Patras is founded on Quaternary deposits, Plio-Pleistocene sediment, and unconsolidated soils [55], such as swamp, coastal and fluvial sands and gravels, Holocene deposits, and the weathering mantle of Plio-Pleistocene sediments, as well as alluvial deposits.The Patras area falls within the Corinth-Patras rift zone, which is about 100 km long and trends west northwest-east southeast across the Hellenic Mountain range.The rift zone is considered one of the most rapidly developed inter-continental rifts in the Mediterranean region, with normal faults trending west northwest-east southeast and with a maximum surface trace exceeding 25 km [56][57][58][59][60].The various catastrophic geological events, earthquakes, and landslides, that have affected the area are closely related to fault tectonics and to the geodynamically active grabens.The city's drainage network mainly consists of three rivers, all flowing in the Gulf of Patras, the Diakoniaris, Glafkos, and Milichos.Patras city is classified as a region with a Mediterranean climate.It has an average annual rainfall of around 830 mm, with seasonal peaks concentrated in autumn and spring and minimum values in summer.City development typically occurs without preventive measures for natural hazards, despite the scientific knowledge of areas where geological or geomorphological hazards are present.However, there is scientific knowledge of the areas where geological or geomorphological hazards are visible.So, urban centers are often founded on steep slopes, or in poor soil conditions, even in areas with high seismicity.
The Patras case study explores the interconnection between the geomorphology and urbanization of a coastal area, locating and analyzing the behavior and interactions of the landforms of the city.Generally, the city underwent fast and continuous expansion in the mid-1800s, which was further accelerated in 1960s.The city's expansion resulted in the diversion of the river channels (Figure 5, left).Additionally, it led to the rivers and coastal areas being filled to create the city's infrastructure.Since 1960, the erosional status of our coastal study area resulted in coastal retreat of up to 7.5 m (Figure 5, right).The coastal belt is an area in which the interaction between the marine and land areas is caused by biotic and abiotic components.Moreover, climate change alters the basic conditions under which shores evolve.In both cases, estimates of long-term or short-term displacement of the coastline are crucial to mitigate the negative impacts and protect the coastal structures from corrosion and other processes that alter the coastal stability.In our case, the coastal belt was exposed to storms from strongly bimodal winds from the northern to western directions.

São João da Barra
The city of São João da Barra (Rio de Janeiro, Brazil) has a historical context and space-time dynamics that are different from the other case studies in this work.During the medieval period, for example, all the other cities mentioned here already had complex spatial organizations and considerable urban infrastructure, while São João da Barra was just a town characterized by the presence of indigenous peoples.In 1676 [61] the place came to be categorized as a Village or "Vila" (in Portuguese terminology), with more established socio-spatial relations and economic activity based on fishing, cattle raising, and the rise of sugar cane cultivation.During this period, and until the 18th century, several access roads were created following the village expansion; moreover, a river port arose due to increased sugar cane production and the need for its transport.It is noteworthy that São João da Barra is located south of the mouth of the Paraíba do Sul River, an important watershed of approximately 55.5 km 2 [62]; nowadays, it is a rapidly growing small urban area hosting about 33,000 inhabitants.The city's geological-geomorphological framework is characterized by Quaternary sediments related to coastal plains, lagoons, mangrove deposits and sandy sediments, originated by different evolutionary phases of the deltaic system during periods of transgressions and regressions of the coastline, with a fundamental role of fluvial sedimentation due to the tropical climate [63].Morphologically, beach ridges, dunes and sandbanks are found, in addition to river and fluvio-marine plains with a predominance of old beach deposits [64][65][66].
In the current context, paleoenvironmental studies have been carried out to understand the genesis and evolution of this deltaic complex; however, the approaches are diverse and still far from reaching a scientific consensus [65,67,68].In this sense, we highlight the efforts made by [65] through a geochronological approach; even without presenting the complexities associated with the deltaic evolutionary environment, they made it possible to carry out more detailed studies, such as that of [67].In recent decades, the city of São João da Barra has stood out in the Brazilian and international scene because of the intense process of coastal erosion that has been progressively consuming the coast, mainly in the south of the delta.Regarding this aspect, [69,70] carried out a multi-scale analysis to understand, using a systemic approach, the possible causes of the coastal erosion that has been plaguing the local community.In this sense, the main causes are inherent to anthropogenic pressures in the hydrographic basin (dams, water capture, and other types of use), and to a lesser extent due to local factors such as the construction of the Açú port complex 7 km south of the delta, whose construction greatly altered the sedimentary dynamics in the region.As a result, an area of intense coastal erosion is delimited (Figure 6), consisting of 458,124.65 m 2 , in the northern portion (first north-south transect in Figure 6).In this sense, it is possible to infer that in the period of 37 years from 1985 to 2022, the sea advanced an average of 12,381.74m 2 over the continent in the locality, which is a significant loss when compared to other areas of the Brazilian coast.In the aerial photograph in Figure 6 (lower image), it is possible to visualize one of the coastal sectors devastated by erosion.
São João da Barra has great touristic potential thanks to the diversity of biotic and abiotic elements present in the landscape.However, what actually draws the most attention is the erosion and the related destruction evident due to the sea ingression.The intense coastal erosion has a negative local socio-environmental impact: the city is visited by groups interested in observing the devastating scenario induced by man.

Potenza
Potenza municipality is the capitol of the Basilicata region (southern Italy), classified as a small urban area of some 67,000 inhabitants.It is in the axial, active seismic belt (30 to 50 km wide) of the southern Apennines, characterized by high seismic activity and where strong earthquakes have occurred.From a geomorphological point of view, the study area is located on a long and narrow asymmetrical ridge southwest-northeast oriented, delimited along the northern sector by steep escarpments.This ridge hosted the city's historic center, partly still delimited by the medieval walls, until the post-war period.Subsequently, urban expansion forcefully invaded the rural and peri-urban territories, obliterating much of the original forms, fundamentally due to river terracing, landslides, and erosive incisions on the sides of the hill ridge.Geologically, it is characterized by a sedimentary sequence of Pliocene deposits included in the Potenza Basin [71], characterized by two sedimentary cycles known as the Altavilla and Ariano units (Figure 7 left).The older Altavilla Unit (Early-Middle Pliocene) consists of conglomerates and sands and occurs predominantly in the northern part of the basin.The younger Ariano Unit, of Middle to Late Pliocene origin, occurs in the southern part of the basin.This unit mainly comprises clinostratified conglomerates, marine sandstones, and calcarenites with an over-consolidated clayey substratum.The urban area of Potenza town is characterized by a clayey substratum and sandy conglomerate deposit layers, which varies in thickness along both west-east and north-south directions (Figure 7).The geodatabase implemented for this study case consists of a remarkable search of unpublished data, such as a topographic map (1:4000), a geological and geomorphological map, borehole logs, geotechnical parameters, geophysical data, historical photographs, and historical macroseismic data at the building scale.The stratigraphic data, deriving from the recovery of dozens of continuous core drillings, made it possible to reconstruct a general geological model of the site, including the trend of the clayey and overconsolidated Pliocene substrate, as well as the thickness of the sandy conglomerate cover.These data have made it possible to exhume a geomorphology buried and modified by urbanization (valleys filled with fills, areas with stagnant water and peat in probable counterslope landslide areas, river terraces located at multiple levels, and ancient landslides).The stratigraphic and geomorphological structure of the sub-flat hill on which the town stands, were also responsible for the historical local seismic response, the effects of which were highlighted by the reconstruction of macroseismic damage scenarios analyzed with geostatistical techniques [72,73].
Applying autocorrelation techniques [74] it was possible to evaluate several relationships, comparing between them the data relating to the damage to the buildings, the degree of clustering and the DEM of the clay substrate, reconstructed using the historical borehole data (Figure 7, left).In particular, it was possible to observe the level of autocorrelation between the level of seismic damage (from the 1857 earthquake) and the geomorphological and geological-stratigraphic conditions of the urbanized hill (understood as the variation in the depth of the overconsolidated clayey substrate, from −3 to −18m, and the presence, above it, of a saturated sandy-conglomerate layer).
Three different areas were thus distinguished (Figure 7, right), characterized as follows: 5.The simultaneous presence of autocorrelation (and therefore damage, shown by white points) and high values in the thickness of the conglomerate cover (14-18 m); in this area the most important structural damage is recorded linked to the amplification of the horizontal seismic acceleration due to morphological (ridge effect) and stratigraphic causes.6.The high thickness of the conglomerate cover and poor presence of positive autocorrelation; this is probably due to the limited influence of the geomorphological factor (top of the relief wide compared to the base), which reduces the possibility of seismic amplification (Pagano square and Largo Pignatari).7. The presence of high localized damage and low values of the thickness of the covering (Porta Salza and Torre Guevara).Here the geomorphological factor (elongated and narrow ridge with steep slopes) seems to play an important role in the amplification of horizontal seismic acceleration, determining high levels of damage although concentrated.
Figure 7. On the left: Geological-stratigraphic structure of the historic center of Potenza (modified from [71]).Below left is an old photo from 1857 (with evidence of the earthquake damage), which shows the asymmetric ridge on which Potenza stands and the entire non-urbanized southern side of which the original morphology can be observed.On the right: Geological schematic section and model along the main axis of the Potenza hilltop town (section A-A', west-east direction), with evidence of the 3 main areas analyzed and the digital model of the top of the clay substrate.

Analysis of Landscape Historical Evolution and Paleomorphology
Anthropogenic geomorphological processes create landforms that are completely artificial or natural features modified by human activities, so it is fundamental to build up strategies for geomorphological surveys of landscapes that are completely transformed.Reconstructing the natural morphology is a fundamental step.The methodological approach in [26] promotes multitemporal topographic data processing and the multidisciplinary analysis of anthropogenic erosion and accumulation landforms.This qualitative approach represents a strength of natural morphology reconstruction.The case studies described in this section show examples of a quantitative approach to natural morphology reconstruction, using the creation of multitemporal digital terrain models and DEM of Differences to estimate anthropogenic excavation/accumulation volumes, together with multitemporal and multidisciplinary qualitative analysis of the urban expansion.

Perugia
Perugia, with a population of 167,000, is the capital city of the Umbria region.The historic city center is situated atop a hill with a maximum altitude of 493 m a.s.l.From its peak, five ridges branch out towards the south, overlooking the Tiber River valley, which is around 170 m a.s.l.Therefore, Perugia's two main morphological elements are the ridgelines at the top and the fluvial valleys at the bottom.The hill's eastern and southeastern area is characterized by three well-shaped ridges directed towards the north, northeast, and southeast.In contrast, the last two ridges, which point towards the southwest and west, are shorter and less defined.In general, the eastern side of the hill is steeper than the western side, where the morphology is gentler.The rivers that flow down also exhibit different behaviors.The rivers on the eastern and southeastern sectors are shorter, steeper, and far from the equilibrium condition.The headwaters are affected by gravitational and fluvial erosion, resulting in landslides and significant linear erosion along the thalwegs.The rivers that originally flowed along the western sector have been culverted over the centuries due to better morphological attitudes and are no longer visible.After the Second World War, this city area experienced the greatest expansion of its building fabric and infrastructure.Although the fluvial valleys are no longer present, the morphology of the areas close to the historic city center resembles that of a drainage basin.During the most important past civilizations that modified the urbanization of the downtown area (Etruscan: 6th-3rd century BC and in the Middle Ages) the headwaters of the main rivers and the thalwegs have been the sites where the morphology has been deeply changed [75,76].The two primary changes made by humans were filling in the depressions where erosion occurred and leveling the tops of the reliefs.Some of the most iconic places in the historic center of Perugia are therefore forms of "morphological false".Geomorphological surveying can recreate the initial topographic assessment and reconstruct the most significant phases of the morphological evolution of the hill.
To understand the volumes and the perimeters of the areas filled with anthropogenic deposits or levelled to a lower altitude, it is fundamental to have a DEM of the paleomorphology for comparison with the present layout.To do this, old paintings and maps were put together.Where it is possible, georeferencing in a GIS environment, the ancient documents and contour lines have been derived.Moreover, the geognostic surveys have been added as vector points where, in the attribute table, the depth and the thickness of the anthropogenic deposits are collected.Polylines (contour lines) and points (the geognostic surveys) are the mandatory input data for the geostatistical analysis that allow us to obtain a final DEM where the paleomorphology is defined [36].The horizontal resolution and vertical accuracy may vary, depending on the spatial distribution of input data and the interpolation method.The next step is to obtain a DEM for the present topographic surface with a similar horizontal resolution.By subtracting the two models (current minus ancient one) with the algebraic function "minus" in map algebra, it is possible to obtain a DEM of Differences raster [35,77] with three possible scenarios: pixels with values close to zero (unmodified surface), pixels with positive values corresponding to infilling areas, and pixels with negative values for the leveled zones.The absolute value for each pixel measures the thickness of the anthropogenic deposits or the volume of removed soil.
At least eight areas have been identified in Perugia where the paleomorphology differs significantly from the current arrangement (Figure 8).One meaningful example is the Fortebraccio Square, located in the northern part of Perugia close to the headwater of two rivers (Lat 43.11483, Long.12.38990, point 6 in the map on the left in Figure 8); the first one is the Bulagaio river, flowing from the square toward northeast and the second one is the Conca River, moving towards the west and nowadays ducted and covered by a street.Since the Etruscan period, strong headwater erosion has affected the saddle between the two rivers.However, it has remained a crucial point for the urban network.The most important Etruscan arch rises in this area, which was also the only access point to the city for those entering from the north.In 1536, the cardinal Marino Grimani ordered the deep depression to be filled, creating the current Piazza Grimana and connecting the Etruscan walled city with the ancient northern district.In addition, the point where the Etruscan arch was built has been lowered over time to reduce erosion and landslide movements.The result is a paleomorphology very different from the current topography.The applied methodology allows for the creation of a detailed final map (Figure 8, on the right) derived from the map algebra operations applied on two DEMs with a 10 × 10 m horizontal resolution, where the modified areas and the range of the volumes have been mapped.

Benevento
The post-war urbanization of Benevento, an inner city of the Apennines in the Campania region, followed territorial development guidelines aimed at replacing the primordial rural economy and its georesources with modern residential and infrastructural conditions.The geomorphological result of the change has enhanced ancient risk conditions and generated new criticality and hydrological instability.The urban territory of Benevento occupies about a third of the municipal area, mostly consisting of small rural centers connected to each other by a modest extra-urban road network.From a geological and geomorphological point of view (Figure 9, upper box) it developed on the bottom and on the heights of a large tectonic depression deeply incised by the Calore and Sabato rivers and their tributaries.The bedrock, like the surrounding hills, is made up of arenaceousmarly-clayey Miocene deep-sea sediments deformed by the compressive tectonics also responsible for the fold-thrust structuring of the Campanian Apennine Orogen.The depression fill is mostly made up of Pliocene clastic, sandy-clayey, and conglomerate marinecoastal sediments, also deformed by tectonics, on which lie discordant Pleistocene and Holocene alluvial (gravels) and fluvial-lacustrine (sands and clays) terraced sediments.
Geological-geomorphological surveying, historiographic research, and archaeological reports analysis [78] allowed us to find evidence of human occupation of the Benevento area dating back at least to the Neolithic period.However, it was in the Samnite Period (VII-IV century BC) that Benevento established itself as a rural village scattered along the valley floors of the Sabato and Calore rivers.From the fourth century BC, with Roman predominance, a proto-urban settlement developed in an elevated position with infrastructures that was more easily defensible and safe from recurring river floods (Figure 9, lower box, panel A).In the late Middle Ages, the Lombards and later the Vatican domain further encouraged the development of urbanization towards the high hills and towards the east, starting the progressive abandonment of the insecure alluvial plains.This intramoenia urban development paradigm lasted until the post-unification period (1870).From this moment and until 1936, the depopulation of the countryside, the advent of industrialization, the new administrative structure and demographic growth imposed a decisive exit from the city walls, with re-occupation of the valley floors following new development guidelines northward and southward (Figure 9, lower box, panel B).The systematic occupation of the valley floors produced a profound alteration in the river system, whose regimentation allowed the affirmation of factories for the milling industry.The construction of new road axes made it possible to connect the nascent urban districts to the historic center, distorting the oro-hydrographic structure of the hill-valley system, with consequent widespread instability of the slopes which was not slow in manifesting itself in calamitous events.Since 1945, the urban center has undergone a progressive and disorderly radial expansion, with the occupation of all the available space from the alluvial valley floor to the sandy-clayey hill (Figure 9, lower box, panel C).This senseless development has been carried out to the detriment of an already lithologically vulnerable territorywhich has aggravated the erosion (stream erosion, mass movements, etc.) of the steep sandy-clayey slopes of the Benevento hill-and has favored fatal floods due to the stiffening of the river system [79][80][81].The study area is subject to high seismicity, with earthquakes that have repeatedly damaged urban settlements; therefore, the increase in geomorphological hazards induced by anthropization has further increased existing risk conditions.

Dissemination of Urban Geoheritage and Education about Anthropogenic Impact on Urban Sustainability
Knowledge of urban geomorphology can contribute to education about anthropogenic impact on urban sustainability.Urban geoheritage consisting of sites testifying to the anthropogenic geomorphological evolution of the Earth is a quite new field of research and a promising resource for touristic and didactic issues aimed at the comprehension of anthropogenic impacts on natural landscapes.Several of the case studies above also worked on the urban geomorphoheritage inventory, assessment, and valorization, based on the geomorphological investigation's results.
The case studies contributed to the urban geomorphoheritage definitions and methods.The Rome and Perugia inventories contributed to the "urban geomorphosite" definition and distinction in lato sensu, sites of any geological-geomorphological interest within the city, and stricto sensu sites specifically testifying to the anthropogenic morpho-evolution of cities [82].The Rome case study highlighted the importance of cultural and geohistorical values in urban geomorphoheritage assessment methods, proposing criteria [83] that are based on the urban environment's peculiarities and the geo-scientific value of anthropogenic geomorphosites as testimonies to the local Earth's history, especially related to the Anthropocene.
The implementation of geotourist products is a tool for geoheritage valorization, Some examples are the Rome GeoGuide app [84], a digital geotourist itinerary across the city historical center, and the Perugia Upside-Down exhibition [75], a tool for the popularization of the city geological context.
The AIGeo Working Group on Urban Geomorphology has also proposed a series of educational books on the urban geomorphology of the main Italian cities: Rome [85], Palermo [86], and Naples [87].

Discussion
At the beginning of this manuscript, we underlined the significance of the urban environment as an ecosystem in which the urban residents interact with biotic and abiotic factors [8,9].Considering the ongoing, increasing demand for land and the expansion of the built-up areas involving and disturbing any kind of natural system inside and surrounding the actual urban areas, it is not negligible how important are the urban environment dynamics and knowledge of its physical evolution.In this context the authors presented their methodological and applicative contribution to the urban environment investigation aimed at the knowledge of the physical processes; geomorphological, historical, and anthropogenic ones in particular.
The eight case studies presented in this paper contributed to the urban geomorphological knowledge of four topics.The hypothesis is that the AIGeo WG on Urban geomorphology methodological approaches and applications can result in a guideline for physical urban environment knowledge and planning on a European and global scale.With this aim, a discussion is open about the usefulness and applicability of the approaches in other cities in the world.
The surveying and mapping of geomorphological processes and landforms in urban areas highlights some differences with respect to the analog in the natural environment.The anthropogenic processes in particular act over very short timescales and with invasive transformation of the landscape and its natural dynamics due to anthropogenic erosion and accumulation.This is important as a base for knowledge of the landscape evolution of an area and its implications.In Europe, population growth and human activities have gradually transformed natural geomorphological processes over the centuries; meanwhile, the actual population growth is more rapid and intense, especially in world regions different to Europe, where the rapidity of construction activities may have dramatically different impacts on natural topography.The strength of the multitemporal and multidisciplinary approach shown by the case studies in Section 3.1 lies in recognizing a very large variety of human activities and related anthropogenic processes and landforms in historical cities.The same human activities have had the same impact on the geomorphological evolution of urban areas since time immemorial, and it happens the same way everywhere.The wide anthropogenic landform variety is similar worldwide, and all the cities around the world can be classified and mapped thanks to the proposed multitemporal and multidisciplinary analysis.The scientific literature in Table 1 supports our opinion (see also Supplemental Material-column topic code, topic code a).It is evident in [88,89], describing and mapping anthropogenic landforms in the municipality of João Pessoa (Brazil), similar to the Italian ones, and in [90] (Vienna, VA , USA), using multitemporal topographic data (stereo aerial and ortho-imagery) for anthropogenic erosion and accumulation detection.Numerous other examples are available from Poland [91][92][93][94][95][96].The proposal again, and in almost all cities around the world of the same landforms, offers the possibility of applying the Italian approach to cities of recent evolution or sprawling ones.
The geomorphological analysis supporting hazard assessment is important considering that human activities hardly transform natural dynamics in urban areas, inducing natural hazards and considering the high risk related to the huge amount of population, infrastructure and heritage exposed.The strength of the examples of urban geomorphology in support of hazard assessment presented in Section 3.2 consists in the natural landscape geomorphological evolution reconstruction, which is transversal to the different hazard analysis methodologies based on geotechnologies.The reconstructions often offer the key for correct geological interpretations, for urban and emergency planning in particular.The proposed methods can be replicated in multiple contexts and allow the creation of urban zoning in terms of measurable risks and consequently also have an impact on the properties' value due to their position and typological construction quality.Furthermore, the multi-temporal investigation and the integration with stratigraphic data allow us to better define the past impact of morpho-evolutionary processes, also as a function of different climatic conditions.The international scientific literature in Table 1 offers a few examples (see also Supplemental Material-topic code b) using the multidisciplinary and multitemporal approach: only in [19] (Khartoum-Omdurman, Sudan) did the archaeological value of the study area require this kind of analysis.Considering what we stated above about human activities and natural dynamics, in our opinion the interaction between anthropogenic processes and natural ones over time should be considered in natural hazard assessment.
The analysis of landscape historical evolution and paleomorphologies in urban areas provides useful knowledge for a better comprehension of the interactions between natural and human-driven geomorphological processes, and for a more refined understanding of Holocene landscape evolution.The multitemporal and multidisciplinary analysis of topographical data and the DoD methodology may be successfully applied in further historical urban centers, where historical topographic maps with unfailing planoaltimetric data have been surveyed worldwide during the 19th and 20th centuries.The DoD, in particular, provides useful elevation information for more refined geological mapping, since trends in natural stratigraphic unit boundaries are required to be drawn consistently with the paleotopography, which in an urban setting is very often nullified by anthropogenic erosion or buried by anthropogenic deposits.In addition, the thickness evaluation seems to be valuable for urban planners and archaeologists.Compared to the previous studies in Table 1 (see also Supplemental Material-topic code c), the strengths mentioned above have not been fully applied to the analysis of interactions between natural and humandriven geomorphological processes.Several studies have used historical relief models and historical topographic maps for cultural heritage documentation [97][98][99][100]; meanwhile, many more studies have used this kind of data for the cities' landscape historical evolution but not aimed at analyzing the natural and anthropogenic morphogenesis interactions [101][102][103][104][105][106].
The case studies regarding the dissemination of urban geoheritage value presented in Section 3.4 deal with urban geotourism and geoeducation.The strength of this kind of dissemination consists in storytelling about the anthropogenic impact on urban sustainability, which is incisive thanks to the attractiveness of historical cities' heritage and mix of natural and cultural values.In Table 1 (see also Supplemental Material-topic code d) mainly European examples are available on the topic [107][108][109][110][111][112][113][114][115][116][117], even though every city in the world has this opportunity.The relationship between nature and culture is certainly very strong in the European cities and it is important for successful geotourism.Thus, from the geoheritage and geo-education point of view, the Italian and European case studies should form a guideline on the topic to be replicable everywhere.

Conclusions
The main aim of this paper was to illustrate an effective approach to urban geomorphological investigation and propose it and the results as a model and example for physical urban environment knowledge and planning on a European and global scale.The geomorphological research into the urban environment has potential for enhancing basic knowledge of natural and anthropogenic morphoevolution (surveying and mapping urban geomorphology) and for making cities and human settlements safe, resilient, and sustainable, as planned in Sustainable Development Goal 11 (geo-hazards assessment, historical evolution and paleomorphologies analysis, disseminating human activities' impact on the landscape through geoeducation).The six Italian case studies and the two international ones (Greece and Brazil) describe the same approach and several applications in similar but different urban contexts, addressing the physical environment of modern cities in a multi-temporal, multidisciplinary and critical manner with respect to global changes.In the discussion, the following strengths of the approach are highlighted: the multi temporal analysis of multidisciplinary data allows the analysis of anthropogenic erosion and accumulation of types and rates and their impact on natural morphodynamics, and to plan and manage the urban sprawl adequately: historical topographic information was integrated with recent findings, and comparison by overlay highlights the variety of human activities and the resulting landscape transformations; integration with historical and archaeological data allows us to recognize the geomorphological constraints to settlements, the solutions adopted, and to hypothesize regarding the suitable directions of the urban sprawl in the past and currently; the investigation of these factors enhances the knowledge of hardly transformed natural dynamics and the interaction between anthropogenic processes and natural ones in triggering natural hazards; the above mentioned analysis and related results are useful for inventorying urban geomorphosites, testifying to the anthropogenic geomorphological evolution of the Earth, and it allows the planning of geoducation and geotourism activities.
In conclusion, the actualism principle applies for anthropogenic processes too.In fact, the human activities acting nowadays have acted since the beginning of human settlements and will act similarly also in the future.It requires analysis of historical city expansion over the centuries and learning from best practices and, moreover, from bad ones.For this reason, the multidisciplinary and multitemporal approach should be a basic guideline for urban environment geomorphological analysis, useful for several applications such as the ones presented in this work.
Furthermore, urban geomorphology is directly related to the historical context of the land use and occupation of these cities, which occurred with different dynamics in time and space.This fact allows the knowledge and material produced in works like this to be used as one of the foundations for urban planning and the prevention of negative environmental impacts.However, it is important to highlight that this must be a continuous effort, as the complexity inherent to the theme, given the breadth of urban contexts and associated global changes, demands special attention to improve the quality of life and environmental challenges arising from anthropogenic pressures.

Figure 1 .
Figure 1.Location map of the selected case studies.On the top, a global view of the countries involved; below left, the 6 Italian cities selected all along the Italian peninsula: Genoa, Perugia, Rome, Pozzuoli, Benevento, and Potenza, and the Greek city-Patras; below right, the city located in southern Brazil-São João da Barra.

Table 1 .
Number of papers per country per topic, collected from the core collections of Scopus and Google Scholar.The full collected literature is available as Supplementary Material.The following four topic codes summarize the main topics in the papers: a-surveying and mapping urban geomorphology; b-geomorphology and urban hazards; c-historical evolution and paleomorphologies; d-disseminating urban geoheritage and educating.

Figure 2 .
Figure 2. Schematic drawings of the landform diversity in Rome, focusing on still recognizable natural landforms modified by human activities (a) and anthropogenic landforms of accumulation (b-d), erosion (e), and mixed examples of both (f).In the top right corner, a geomorphological map excerpt and related legend are depicted showing the mapping of the above-mentioned landforms.

Figure 3 .
Figure 3. On the left, drainage network modifications and flood prone areas of the morphological amphitheater on which the historic center of Genoa lies.Legend: (1) poorly modified and/or natural riverbed; (2) culverted stream; (3) concrete channel; (4) eaves channel; (5) abandoned channel; (6) exposed buildings; (7) flooding area with returned period > 200 years; and (8) historical flooded area.On the right, the figures show the upper sector of the Lagaccio stream valley affected by relevant man-made morphological modifications: (a) the present-day situation, with the sports facilities on the fills along the stream; (b) the Lagaccio dam lake in the 1960s; and (c) geomorphological section (BH = boreholes).

Figure 5 .
Figure 5. Patras (urban area), damage recovery on road in coastal zone affected by powerful sea erosion in November 2021.

Figure 6 .
Figure 6.Satellite image of 1985 coastline overlapped with the current one, showing the areal gap due to erosion (modified from [69]).Below: coastal sector of the municipality of São João da Barra (Rio de Janeiro) consumed by coastal erosion.The yellow dotted line represents the past position of an avenue (N-S view) (Photo: Laboratory of Geomorphology UNICAMP collection, 2023).

Figure 8 .
Figure 8.The slope angle map of Perugia city is on the right of the figure.The map clearly shows the areas with the highest slope angle values (in red), which are all located in the eastern sector of the city or on the eastern facing slope of the fluvial valleys.Additionally, the map indicates areas with a paleomorphological order (black numbered dots) that differs from the current one.All the areas are located close to the downtown area, where the highest and oldest settlement is limited from the ancient Etruscan Wall.On the left side of the image is Grimana square (point 6 in left figure); (a) the initial topographic layout pre-urbanization; (b) the actual topographic layout; (c) volumes of filled material (in red) and eroded areas (in blue).

Figure 9 .
Figure 9. Upper box-Schematic geomorphological map of the Benevento urban territory.Lower box-Simplified maps of the urban expansion of Benevento town from Roman times to the presentday.