1.Introduction

One of the most discussed concepts in Landscape Archaeology (an “enormous topic” [1]) is the archaeological (inductive) inference, the process of defining the characteristic of the whole context, starting from the analysis of a part of it. Firstly, this is related to the practice of adopting a sampling strategy to analyse the object of the research (the population) and, secondly, to a limit of the archaeological investigation: the sites or artefacts detected on the surface never are the totality existed in the past because several biases affect this result.

The wave of interest in sampling theories and bias effects started in the 1960s [2] and reached its acme in the middle of the 1980s [3], when the post-processual critique started [4]. Particularly in the Americas, trying to avoid the risk of omitting key sites - the so-called “Teotihuacan effect” [5] -, the archaeologists started to adopt a “full-coverage” approach, investigating the entire context, anyhow with a low level of intensity [6]. On the other hand, in the Mediterranean countries, the pedestrian surveys still prefer a high level of intensity and investigate small area of 10-25 sq. km [7]. This was a consequence of the debate on the biases affecting the survey results, a discussion particularly strong in Mediterranean archaeology; in the last two decades of the 20th century, to augment the intensity of surveys seemed the easier solution [8]. Furthermore, the difference between American and Mediterranean archaeology is linked also to a different approach to the landscape. In American archaeology the site remains the fundamental unit of analysis, while European and Mediterranean archaeologists “have moved more clearly toward the landscape approach, with individual features or even artefacts as the basic element of analysis” [9].

The impact of the biases affecting the results of a survey reaches an acme when archaeologists investigate an alluvial floodplain [10]; here, only an integrated approach between different methodologies (geomorphology, topography, artefacts survey, archaeobotany ..) can help to shed a light on past landscapes history and to limit the impact of different kinds of biases [11].

This paper aims to outline the main limits that a landscape archaeological project has to face examining a floodplain with a high rate of vertical growth of ground levels, due mainly to the combined action of subsidence and alluvial sedimentation phenomena. Indeed, the case-study area is a portion of the Po Valley lowland, the Ravenna hinterland, equal to approximately 525 square km and located near the northern coasts of the Adriatic Sea (Fig. 1). This region is currently called “Bassa Romagna›› and since 2009 it has been the research area of the “Bassa Romandiola›› project, which included artefact surveys, geophysics analysis, archaeobotanical and geoarchaeological researches, and written sources analysis. The limits faced here are characteristic of a wide floodplain area and the solutions applied would be helpful for other archaeological research projects.

The first section of the paper is devoted to present the geographical and geomorphological characteristics of the case-study area, the project history, and the research methods; the second section is focused on the characteristics of the archaeological data collected, discussing in particular their limits and the influence of the main biases detected; the third and final section presents some of the solutions applied in order to overcome these limits.

2. The case-study area and the research project

2.1. Geography and Geomorphology of Bassa Romagna region

The Romagna region corresponds to the south-eastern part of the Po Valley; it is today more a cultural and historical region, composed of the provinces of Ravenna, Forlì, and part of that of Bologna. The natural borders of this region are the ancient river Po of Primaro (Reno nowadays), north of Ravenna, the Sillaro river, westward, the Apennine ridge (approximately) and the river Tavollia southward, and, finally, the Adriatic coastline to the east (Fig. 1) [12]. The lowest part of this area is currently characterised by the presence of intensively cultivated lands, extending between several urban centers and towns densely populated (Fig. 2). This lowland was the result of the action of the aggradation of Apennine rivers, interacting with the southern branches of the Po river. From the end of the 13th century this area was identified with the name “Romandiola” and the region started to have a clearer territorial identity [13]. “Bassa Romagna” (namely Low Romagna) is currently the northern and lower part of Romagna, located westward of Ravenna; it corresponds to a historical region called “Territorium Faventino acto Corneliense” during the Middle Ages and, from the 15th-16th centuries, “Romandiola Estensis”, because it was the portion of “Romandiola” that went under the control of the lords of Ferrara, the Estes.

From the geological point of view, the soil surface in Bassa Romagna is exclusively molded on alluvial and swamp sediments mainly represented by the layers of the Subsintema of Ravenna (AES 8) - dated back to last 15.000 years - and in particular the upper part of it, the Unità of Modena (AES 8a), which includes the sediments born from the 4th century CE onwards [14]. The intense activity of alluvial sedimentation of the Apennine rivers is here accentuated by a high rate of subsidence, which reach one of the highest values of the southern part of the Po Basin [15]. The main consequence is that the ancient ground-levels are sometimes buried several metres deep; for instance, sometimes in the investigated region the Late Republican/Imperial stratigraphies are buried 10 metres deep [16]. Moreover, this aspect is characterised by a high rate of variability: this depends on the morphological characteristics of the ancient landscape, which alternated highlands (e.g. inactive fluvial ridges) and depressed area (e.g. wetlands) [17]. These elements also depict a very unstable territory from the geomorphological point of view, already before the Roman Age contrasted by the anthropic actions, which were finalised to manage the landscape to assure water drainage and supply [18]. This makes the archaeologists’ job in this kind of contest more complex. They are obliged since the preliminary phase of the research to have a systematic and wide-range approach, to take into account several possible parallel interpretations, and to consider in a single analysis framework different kinds of data — the only approach that can limit the imponderable.

2.2. Project history and research methods

The landscape archaeology project “Bassa Romadiola” was started in 2009 and, until now, it investigated through intensive artefact surveys approximately the 15% of the territory of Bassa Romagna (i.e. 78 sq. km of the whole context, namely 525 sq. km). Since the beginning of the research several field campaigns took place (Fig. 3) [19]: in 2009, in the area of Lugo [20]; in 2011, in the area of Conselice [21]; in 2012, in the area of Bagnacavallo and Fusignano [22]; in 2016, in the area of Bagnacavallo and Lugo [23]; in 2018, in the area of Cotignola.

The objective of the project was to define the history of the settlement patterns in this portion of the Ravenna hinterland, from a diachronic perspective; in particular, we focused into trying to solve a gap in the historiographic research, which was not able to clearly define the settlements patterns’ transition from the Roman Age to the Middle Ages.

The on-field research method was that of the systematic and intensive survey [24]. In the sample area, every cadastral parcel with a minimal level of visibility (i.e. a “Topographical Unit”) was systematically surveyed by the archaeologists, walking in parallel lines with a distance of 10 metres from each other (the “Intensity” of the survey) [25]. Thus, the “Topographical Unit” was considered the box in which any archaeological find was located, as well as, anyway, the trace on the ground of human and natural activities [26]. During the artefact survey the position of every single find has been recorded with a GPS receiver, like in a siteless survey [27]. Anyhow, any area with a higher density of artefacts, above the background scatter of the Topographical Unit, has been interpreted as an archaeological “Site” [28]. In this case, the intensity of the survey increased and the distance between every archaeologist was set to 1 metre from each other; every archaeological find has been collected, except the architectural ones (i.e. bricks or roof tiles). This double approach (site/siteless) allowed to map the continuum of the artefact distribution [29], preserving anyhow the interpretative category of “Site”, useful for further archaeological analysis, the public bodies policies and the related instruments of protection.

3.Why “Looking through the keyhole”? Limits and characteristics of the collected data set

3.1.First of all: sampling the context

The research started already before fieldwork, with a review of the previously known data, analysing written sources archives, previous archaeological discoveries, and the geomorphological evolution of the context [30]. We were not working in a vacuum, but no archaeological field systematic research was done before in this area. The aim of this phase was to acquire the existing knowledge on the investigated context and to define the sampling strategy of the artefact survey. The data set produced was a geo-database in which each record was associated with a geometric geo-referenced feature (points, lines, and polygons); it included 625 entries, divided into the following categories (Fig. 4 and Tab. 1):

Categories Typologies Geometrical features
Sites and landscape features mentioned in edited written sources Fortified settlements (towers and castles); Isolated houses; Churches (ecclesiae, plebes); Harbors; Hospitals; Productive sites; Land exploitations structures (curtes and massae); Mills; Monasteries and convents; Open villages; Marshes; Channels; Forests Points, lines, and polygons
Archaeological sites Burial areas; Isolated houses; Castles; Churches; Land exploitations sites; Productive sites; Monasteries and convents; Harbors; Nucleated settlements; Channels, bridges, and fords; Road structures; Isolated finds Points
Geomorphological features Fluvial ridges and Fluvial branches, according to edited studies and their hypothesis Polygons and lines

Table 1. The preliminary phase of the Bassa Romandiola project had the objective of summarising the existing knowledge on the investigated area and produced a spatial database; here the entries’ categories:

In this kind of landscape, the main challenge is to investigate the history of fluvial systems; scholars tried to clarify their recent evolution, but the unresolved questions are still numerous [31].

By crossing archaeological, historical, and geomorphological data, this first phase allowed to divide the research area into Landscape Strata, belonging to three main categories (Fig. 4):

To adopt a sampling strategy was the only possibility in order to balance resources (time and costs), with the objectives of the research. Furthermore, the extreme variety of the examined area forced to change the sampling strategy in each stratum, as summarised in the following table (Tab. 2):

Landscape Stratum Geomorphological, archaeological, and historical main characteristics Sample characteristics
i. High floodplain Stable floodplain; Roman Age sites in surface or buried maximum 3-4 metres deep; centuriation mainly conserved Rectangular transect oriented according to the centuriation (approx. N/W-S/E) and a within-stratum sampling fraction of 50/100
ii. Fluvial ridges in the low floodplain Stable floodplain; Roman and Bronze Age buried under plough soil; Roman centuriation retraced in the Late Middle Ages or not conserved and replaced by medieval or modern land divisions Rectangular transects oriented N-S or E-W, with a width of 1-kilometre multiples, which cover the whole fluvial ridges, and a within-stratum sampling fraction of 100/100
iii. Low floodplain Unstable floodplain, generally marshes dried in Modern and Contemporary Age; Roman and Bronze Age buried from 5 to 10 metres deep or more; centuriation mainly absent Area excluded from the sample

Table 2. Landscape Units categories in Bassa Romandiola projects and sampling strategy:

Thus, by choosing different sampling fractions in each stratum - but with a systematic approach -, we applied what we can define a systematic disproportional stratified strategy [32]. Anyhow, we chose high fraction values, from 50% to the full-coverage, because the previous archaeological knowledge was very superficial; moreover, since the designing phase, we left the possibility of doing some minor enlargements of the sample during the fieldwork, with an adaptive strategy [33].

The sampled area covered 150 sq. km, equivalent to 22% of the Bassa Romagna sub-region, and it avoided the “Low floodplain” stratum (Fig. 4) [34]. The transects in the High floodplain stratum are 8, designed according to the main cadastral limits of the landscape and a sampling fraction is 50/100. On other hand, the transect in the “Fluvial ridges in low floodplain” stratum are 15; they are oriented according to the kilometric grid, north-south or east-west, and designed to assure almost a full-coverage of this stratum.

3.2. The “keyhole”: the data set collected

In the ten years of research 71 artefacts concentrations have been detected. After an interpretative phase, they have been classified into different categories, according to their dating and characteristics, as summarised in the following table (tab. 3).

Sites categories Total Number General Chronology
Castle 4 11th/13th-15th
Church 4 11th-15th
Early medieval-High Medieval house 22 8th-12th c.
Late Antique-High Medieval house 4 5th-12th c.
Late Antique house 1 5th-7th
Late Medieval-Modern Age farmhouse 8 14th-18th
Late Medieval-Modern Age productive site 6 14th-18th
Medieval Harbor 1 11th-15th
Medieval not-fortified village 2 11th/12th-14th/15th

Table 3. Archaeological sites detected (2009-2018) [35]:

It is not my aim to discuss in this paper the settlement patterns evolution in this area or the historical phenomena connected to that, e.g. the land reclamation processes [36]; as anticipated, I would like to focus on the methodological aspect of surveying such a kind of landscape. Just for the sake of clarity, the project allowed us to find an answer to our first question: what happens in this part of Romagna after the end of the Roman Age and how the settlement patterns evolved. Summarising, the artefact survey identified a relevant growth in the number of the sites between the 9th-10th centuries, the Carolingian and Ottonian Age, sometimes in continuity with Late Antique sites. This process took place according to two main trends: a dispersed pattern (as recorded in the area of Zagonara, Fig. 5a), as well as a clustered pattern (as recorded in the area of Bagnacavallo, Fig. 5b). The clusterization of the habitat appears to have been a frequent trend in Romagna lowland during the Early Middle Ages, as already documented in other contexts of the Ravenna hinterland itself or in the area of Cesena (40 km south of Ravenna) [37].

The main bias that clearly affected survey results for such a context of a plain with a high rate of vertical growth was a macroscopic factor, acting at site-level [38]: the geomorphological processes. The results of the artefact survey are drastically limited in comparison to the surveyed area. The most ancient sites have been recorded in what we can define as geomorphological windows of less than 10 square km. They are areas in which portions of natural paleo-soils are exposed on the surface, or limitedly buried by alluvial phenomena (so still reachable by plowing). Generally they correspond to paleo-fluvial ridges, higher than the rest of the lowland around. After their deactivation, they were chosen to place settlements and start land-exploitation [39]: they were safe from floods, easily drainable, a suitable environment for agriculture, and natural and secure paths that crossed lowland and marshes and linked localities. Post-depositional alluvial phenomena can bury these landforms, making them impossible to detect simply through a micro-morphological approach [40].

Thus, it is clear that interpreting these data is similar to trying to understand what is happening inside a room looking through the keyhole.

The area around, indeed, was not an empty landscape: sites presence is attested in the medieval written sources; buried sites are detectable in the literature and located under the plough soil level. All that:

It does not mean that sampling is not useful or effective in this kind of context, on the contrary: despite the limits, it can address the research and make it more efficient. Moreover, it can clarify the sustainability of the research questions, which since the beginning move the archaeologists: if your research topic is to define how the transition from the Protohistoric to the Roman landscape occurred, to survey the northern part of Bassa Romagna is not a good option.

If the impact of geomorphological biases is well-known by landscape archaeologists since the 1960s and 1970s [41], anyhow, the existing maps hardly have enough resolution to represent the actual geomorphology of an area with such a local level of variability, as the Bassa Romagna [42]. Not every single window can be predicted and dated a priori.

Furthermore, the systematic survey itself can complete the geomorphological knowledge of the context in a diachronic manner. During this field campaign, we collected all the artefacts, in-site and off-site [43]. The earliest dating of a site, or, if absent, of the artefacts off-site, provides a reliable terminus ante quem to date the geomorphological stability of a Topographic Unit. With a process of spatial interpolation, it is possible to extend the measure to the not sampled portion of the context (Fig. 6) [44]. Obviously, “one swallow does not make a summer”: a single artefact is not enough for that; a consistent set of artefacts is better, and even more an archaeological site or a set of these [45].

Sometimes the data set collected confirms the previous knowledge; sometimes it can also integrate it. Geologists often base their considerations also on morphometry; the fluvial ridges are mapped frequently detecting anomalies in ground elevation. Thanks to that, scholars mapped the fluvial ridges so-called “of San Pietro in Sylvis” (5th-11th centuries) and “of Bagnacavallo”(11th-13th centuries), in the northern part of Bassa Romagna [46]; these hypotheses were confirmed by the survey and this is clearly visible in the sites distribution map (Fig. 5b). Anyhow, the systematic approach of the survey allowed us to detect other geomorphological windows, assessing and improving geologists’ reconstructions. This is the case of Zagonara area: the existence on the surface of Late Antique and Early Medieval sites suggested the presence of an emerging paleo-soil, probably the top-soil of a fluvial ridge; it was not recognizable through the morphometric analysis of the surface because it has been buried and hidden by following alluvial events and currently its surface is - paradoxically - at a lower level than the area around (Fig. 7); the geoarchaeological analysis of the sub-surface confirmed this hypothesis [47]. This paleo-soil can be probably related to the depositional activity of the near so-called paleo-ridge of Barbiano, partially buried here [48].

The second element that affected sites’ visibility after geomorphological factors was the land use/land cover; the proportion between ploughed lands (high level of visibility) and orchards, vineyard or, worst, grassland (low or null visibility), changes radically, according to the geographical context, agrarian traditions and policies, and geopedological characteristics. The stratum of low floodplain has a prominence of extensive ploughed soils, mostly composed of distal alluvial deposits of clay; anyhow, the geomorphological biases do not allow to take any advantage from that (Fig. 8a). The stratum of the middle floodplain is characterized by a proportionate balance between different cultures (Fig. 8b); anyhow, going southwards, towards the high and stable floodplain, the portion of orchards and vineyards becomes bigger and bigger, according to the more and more intense agricultural practices (Fig. 8c). More the burden of the geopedological biases decreases, more the relevance of land use/land cover biases, unfortunately, increases.

Changing the point of view and taking into account the artefact level (the microscopic one [49]), the main bias is related to the action of mechanical agricultural instruments on artefacts’s distribution and dimensions. The Romagna countryside has been intensively cultivated since the beginning of the 20th century and men covered the area with orchards and cereals fields. Since that time the archaeological deposits on the surface have been destroyed and dispersed; artefacts themselves have been fragmented in smaller and smaller sherds. This aspect is clearly recognizable in the field. The fragmentation index of sherds is relevant in any site detected, with a substantial consequence in quantification [50]. Some areas, with a density of artefacts higher than the normal but lower than the sites themselves, do not seem simply a sort of background noise, but sites with a very low level of visibility, because strongly dispersed by mechanical instruments, or partially affected by the “traffic light” effect: a site is not clearly visible every year, but randomly, as a consequence of the ploughing cycle [51].

It is clear that all the mentioned aspects surely deserve future definitions and quantification analysis. The job of surveyors seems pretty hard in a plain with a high rate of vertical growth such as that of Bassa Romagna.

However, is our data set reliable?

Archaeologists have to deal with the paucity of data and frequently they argue inferring from absence, despite the adage “Absence is not evidence of absence” [52]. For instance, that happens when we suggest that a site have been abandoned in a certain period, because we do not find any diagnostic artefacts of that period [53]. Similarly, what has been recorded on the surface is not representative of the sub-soil and, furthermore, it can be also a false-positive: not all the concentrations are actually a site [54]. Anyhow, the frequency of the negative results can have such a degree of likelihood, to “create situations in which inferences from the absence of evidence have a respectable plausibility” [55]. In our case-study all the sites in panes A and B of Fig. 5, do not have any diagnostic artefact more recent than the 12th/13th century; consequence: it is highly plausible they have been abandoned during the 13th century.

At the same time the statistical frequency of the recurrence of a result - so a positive result and not a negative one - suggests a historical trend: in the area of Zagonara between the 9th/10th the 100% of sites are dispersed, while, in the same age, in the area of Bagnacavallo they are clustered.

This is a trend and the answer to my question has to be yes.

4. Filling the gaps

As expected, in this floodplain with a high rate of vertical growth the main bias is obviously the geomorphological one. What was unexpected was the elevated geomorphological variability of this area that makes any inference process harder. The second bias, in order of relevance, was low visibility caused by intensive agriculture, which even more limited the comprehension of the settlement patterns.

These biases - micro and macro - have been compensated in different ways:

Two final considerations to conclude this section: only crossing all this information - archaeological data set, paleo-geomorphology, written sources, environmental factors - we will be able to inference with a probabilistic approach and, thus, apply bayesian statistics [64]. Second, a multiscaled approach seems the right way to investigate this relevant amount of aspects, at different levels of resolution and with a multidisciplinary approach.

5. Conclusion: infinite parallel worlds?

To produce a synthesis on settlement evolution in this kind of landscape, we are forced to reason in a binary way: 0 or 1, presence or absence. Between the two we have to privilege the presences, avoiding the appeal to inference from absence. Any gradient or nuance brings us into a sloping ground. This happens when we try to develop probabilistic previsions or, simply, interpretations of what happens outside our sample or our geomorphological window. What we saw from the keyhole was not enough, but it was something to start with.

It is something like it happens in the quantum physic: a phenomenon can be in different status at the same time, until we have the possibility of observing it; this brings to an infinite number of overlapping possibilities and, potentially, infinite parallel worlds [65]. Luckily in archaeology, we have (at least until now) something of material to touch or see; the number of parallel worlds can be limited with a long process of integration of a set of data, trying to find the right dialogue between different sources and approaches.

Indeed, also this kind of landscape, a plain with a high-rate of vertical growth, has brought surprises and unexpected results (such as the Late Antique sites on the surface in some northern areas). In that, there is a strength and an opportunity: this kind of context allows to investigate the birth of settlements free of any past legacy. There was no Roman villa survived in elevation in the area of Bagnacavallo when a landowner decided to built his farm or his house between the 5th/6th c.; similarly it happened here in the 9th/10th century, when peasants decided or have been forced to build their houses in a clustered manner. In other terms, within the landscape, there were no previous attractor places, which could assure the continuity [66]. After the extreme events that buried the past settlements, in the 5th as well as in the 9th/10th century, the floodplain was an empty canvas in which people could choose to settle as they preferred.

Furthermore, this kind of landscapes is frequently composed of marginal lands: fragile ecosystems exposed first to extreme events and, because of that, obliged to be more resilient to survive [67]; the study of these historical landscapes can help contemporary societies to face future challenges.

Thus, as in any other landscape archaeological project, also in a such difficult context, a systematic and multiscaled approach is always the best one [68].

To conclude, facing a context like this one, we can suggest to:

Acknowledgments

The Bassa Romandiola project is funded by several private and public bodies, in particular: the Centro di Studio sulla Romandiola nord-occidentale; the Comitato per i Beni Culturali del Comune di Lugo; the cities of Conselice, Cotignola, Fusignano, and Lugo; the Flaminia Foundation; the Cassa di Risparmio e Banca del Monte di Lugo Foundation. I would like to thank Michele Abballe, Alice Ferrari, and Celeste Fiorotto for the help in producing some of the figures of this article; Michele Abballe, Andrea Augenti, Federica Boschi, and Stefano Marabini for reading the text and for the precious suggestions. Several topics included in this paper are the result of the discussions I had with dr. Nicola Mancassola, who I have to thank for mentoring.


Bibliography


Notes

1. Orton 2000, 40, with an overview on the related debate.

2. E.g.: Binford 1964; some prodromes reported by E.B. Banning (2020a, 4): Philips, Ford, and Griffin 1951; Vascelius 1960.

3. Bannings 2020a, 5.

4. Shanks and Tilley 1987.

5. B. Mayer-Oakes and R. Nash at the 1964 meetings of the American Anthropological Association performed a theoretical task: they applied a stratified sampling strategy to an existing project, the Teotihuacan Valley survey project of William Sander (1970), using a stratified random sampling strategy, according to L. Binford indications (1964); the sample missed the most important site of the area, the city of Teotihuacan itself (see the discussion reported in Flannery 1976, 133-135).

6. An overview in Banning 2020a.

7. Alcock and Cherry 2004; Opitz et al. 2015.

8. Terrenato 2004, 37; Cremaschi 2000, 220-26; Cowley 2016.

9. Opitz et al. 2015; cfr. also Volpe and Goffredo 2014; Campana 2018, 31-42.

10. Vita-Finzi 1969, 101-02; Brown 1997, 1-5; Dall’Aglio 2000, 237-38; Franceschelli and Marabini 2007, 76-93.

11. Dall’Aglio 2000, 239; Denham 2017.

12. Longhena et al. 1936.

13. Vasina and Mengaldo 1970.

14. Cibin 2014; Severi and Cremonini 2009; Marabini and Vai 2020.

15. Bitelli, Bonsignore, and Unguendoli 2000; Campo, Bruno, and Amorosi 2020.

16. Franceschelli and Marabini 2007, 78 and 194; Cavalazzi et al. 2018; Abballe 2020.

17. For a synthesis cfr. Castiglioni et al. 1997; Franceschelli and Marabini 2007, 76-93; Marabini and Vai 2020.

18. Regarding the analysed area, cfr. for instance the synthesis in: Bottazzi 2000; Dall’Aglio and Franceschelli 2012; Chouquer 2015.

19. For a summary of the results: Cavalazzi et al. 2018.

20. Cavalazzi 2012; Cavalazzi et al. 2018.

21. Cavalazzi et al. 2018.

22. De Felicibus 2017; Cavalazzi et al. 2018.

23. Cavalazzi et al. 2018.

24. Banning 2002, 89-90.

25. Cambi and Terrenato 1994, 161-79.

26. According to Gattiglia and Stagno 2005.

27. Cambi and Terrenato 1994: 256-57. The device was a GPS Garmin 64s.

28. The threshold of site definition on fieldwork was related to every single Topographical Unit and subjectively defined by the field coordinator, dr. Marco Cavalazzi; further statistical analysis did not change substantially this preliminary distribution.

29. Campana 2018: 31-42.

30. The preliminary phase of the research was funded by the Center of Study on the “Bassa Romandiola Nord-Occidentale” of Bagnacavallo and was performed in the winter of 2008 by the author, under the direction of prof. Andrea Augenti and dr. Nicola Mancassola.

31. Fortunately some scholars tried to overcome this limit, as Stefano Marabini (2020, with Gian Battista Vai, or 2007, with Carlotta Franceschelli) and Michele Abballe (2020); see also footnote n. 39.

32. Orton 2000, 30; Banning 2020b, 45. On the contrary, a proportional stratified sample adopts the same proportion between sample and population (the “sampling fraction”) in every strata.

33. Orton 2000, 34-38 and 90-91: an adaptive strategy implies to change the sample facing some predetermined condition, e.g. to find a relevant cluster of sites, as it was in our case. This actually happened during the 2009 field-campaign, when a cluster of archaeological sites was found in the area of Zagonara -Lugo, RA- (Cavalazzi et al. 2018).

34. Banning et al. 2017.

35. The chronology of the sites can be less wide than indicated; the range reported is the widest defined for the sites included in the category.

36. For a synthesis: Cavalazzi et al. 2018.

37. Mancassola 2008; Negrelli 2008.

38. Macroscopic biases are active at site level, while microscopic biases are active at the artefact level (infra, and cfr. Terrenato 2004, 38-39; Given 2004, 19).

39. In this field the studies are numerous; I will mention just some main references related to area analysed and, more in general, to the south-eastern part of the Po Valley: Cremonini 1994; Veggiani 1995; Marabini, Franceschelli 2004; Rucco 2015, 109-14; Mozzi, Piovan, and Corrò 2020.

40. Bondesan and Meneghel 2004, 133-134; Marabini and Franceschelli, 2004, 28.

41. Vita-Finzi 1969; Ammerman and Bonardi 1981; between the last: Terrenato 2004, 39. An archaeological approach to an extreme situation (all the archaeological deposits buried and hidden by a relevant alluvial deposit) in: Corrò, Moine, and Primon 2015.

42. This kind of problem is discussed in: Terrenato 2004, 39.

43. On the use of off-site scattered artefacts: Bintliff and Snodgrass 1988; Whitehead and Stoddart 1991.

44. With the Inverse Distance Weighting (IDW) algorithm, we interpolated the centroids of every Topographical Unit, associated to the value of the dating of the earlier artefact found in it (Collins 2020, 120-21).

45. See below the discussion on inferencing from absence.

46. Franceschelli and Marabini 2007, 31-32; Marabini and Vai 2020.

47. Abballe and Cavalazzi in press and infra.

48. Franceschelli and Marabini 2007, 31;

49. Terrenato 2004, 38-39.

50. Ceci and Santangeli Valenzani 2016, 20-21; Cirelli 2006.

51. Lloyd and Barker 1981, 291.

52. Sagan and Druyan 1997, 213; for a review: Wallach 2019.

53. Banning 2002, 204-05; Wallach 2019.

54. Dall’Aglio 2000, 237-38.

55. Wallach 2019.

56. Tilley 1994.

57. Terrenato 2004; Fahlander 2014.

58. Terrenato 2000; Casarotto 2019.

59. The research is coordinated by Michele Abballe, Ghent University.

60. Abballe 2020; cfr. Franceschelli and Marabini 2007, 78.

61. Fiorotto 2018.

62. Also this part of the research started recently and it is coordinated by Celeste Fiorotto, University of Verona (Fiorotto et al. 2020).

63. It is impossible to refer to all the studies highlighting the importance of remote sensing in landscape archaeology; here just several recent references: Ceraudo and Boschi 2009; Forte and Campana 2016.

64. Otárola-Castillo and Torquato 2018; for an application to this context: Abballe 2017.

65. Rovelli 2020, 40, 68-69.

66. Tosco 2009, 151.

67. Vandam 2019; referring to this area: Abballe and Cavalazzi in press.

68. Dall’Aglio 2000, 235-36; Mills 2007; Cowley 2016, 167.