I am an archaeologist specialising in the 1st millennium CE Italian countryside, with a particular focus on quantitative approaches and Bayesian modelling of environmental archaeological data.
As part of my doctoral project, I have just released an open access dataset which includes archaeobotanical macroremains and zooarchaeological NISP data from mainland Italy. I also co-coordinate the Northern Apulia Coastal Survey, a long-term project investigating settlement dynamics in the wetland and coastal landscapes of northern Apulia, where I focus on fieldwork methodology and the modelling of survey data.
I am currently working on the Marie Skłodowska-Curie Actions project BTChron which aims to better integrate chronological uncertainty into typo-chronologically dated archaeological datasets and to develop stronger tools for large-scale temporal analysis, ultimately hoping to improve our understanding of agricultural change during the Roman period!
-
Integrating Historical Sources with Climate Simulations: Temperature and Precipitation Anomalies in Southern Italy (1000–1800 CE)
Roberto Ragno, and Roberto Goffredo
Quaternary Science Reviews, Nov 2026
This paper presents the first systematic palaeoclimatic reconstruction of southern Italy, primarily focusing on the Puglia region, from 1000 to 1800 CE. Although documentary sources provide a wealth of climatic information for this area, they primarily record anomalies rather than normal climatic conditions, offering no baseline for comparison. To address this issue, we compiled a relational dataset comprising 526 documentary climate events, 285 of which were selected for quantitative analysis using the Pfister seven-point scale. We developed a Bayesian state-space model integrating these observations with CHELSA-TraCE21k palaeoclimatic simulations, which provide informative priors linked by a first-order autoregressive process. Observations enter the model through an ordinal logistic likelihood. The resulting posterior should be interpreted as a simulation-anchored latent climate state, informed primarily by documented extremes, rather than as a direct reconstruction of centennial annual-mean conditions. The temperature posterior indicates predominantly negative anomalies from the 15th century onward, broadly consistent with a Little Ice Age cooling phase, although the magnitude is influenced by the strong representation of extreme cold events and the 16th-century estimate is largely determined by the prior and temporal persistence. The precipitation reconstruction reveals a more complex picture: the CHELSA-TraCE21k prior exhibits a progressive drying trend not reproduced by an alternative forced simulation (ModE-Sim), while the documentary posterior remains closely aligned with the forced simulation. This agreement suggests that the prior–posterior divergence may reflect a regional feature of CHELSA-TraCE21k, although alternative explanations, including documentary reporting bias and model assumptions, cannot be excluded. Posterior uncertainty reflects data availability, with narrower credible intervals in well-documented periods. We compare both posterior reconstructions with two products from the Modern Era Reanalysis framework (ModE-RA and ModE-Sim), using anomalies re-baselined to a common reference period, and find broad agreement for temperature but a more informative divergence for precipitation. This comparison assesses consistency across climate products and highlights the sensitivity of the inferred precipitation trajectory to the simulation product used as prior, rather than providing an independent validation of either reconstruction. This framework can be applied to other regions where discontinuous historical evidence must be integrated with simulation-derived climate baselines.
-
Classifying Functional Areas at the Roman Villa of Mascherone (Manfredonia, Italy): Modeling Sparse Surface Survey Count Data Using a Bayesian Zero-Inflated Negative Binomial Model
Roberto Ragno, Roberto Goffredo , and Luciano Piepoli
Nov 2025
Statistical approaches to the study of archaeological surface survey datasets are problematic, as they are often characterized by excessive zero counts and overdispersion. This paper introduces a Bayesian Zero-Inflated Negative Binomial (ZINB) model to classify survey grid units into user-defined functional areas based on artifact distributions. This approach was applied to the Roman maritime villa of Mascherone, located near the city of Siponto (northern Apulia, Italy), surveyed using a total sampling strategy. After filtering, 52 squares (20 × 20 m) containing 21 distinct artifact types were analyzed to identify three hypothesized functional areas: residential, storage, and craft. The model explicitly accounts for structural and sampling-derived zeros in the dataset while also handling overdispersion. Furthermore, it provides probabilistic classifications with quantified uncertainty for each square unit. Results indicate a residential core consistent with legacy aerial evidence, while storage and craft zones remain less certain due to limited indicators. This approach effectively addresses zero-inflation in survey datasets and offers a scalable framework for broader archaeological and landscape analyses.
-
Peasants, Agriculture, and Environment in the 1st Millennium CE Italian Countryside: A Bayesian Approach
Roberto Ragno
Nov 2025
This dissertation explores the subsistence methods, economic systems and environmental adaptations of Italian peasant communities in the 1st millennium CE, with a particular focus on the transitional period from the Roman Empire to the early medieval era. Existing work on agricultural production in this period has been based on literary sources and field surveys, or has focused on individual sites or regional collec-tions, while a multi-source archaeological study is absent from the discourse. This work addresses this gap by using environmental proxies to reconstruct the historical agricultural landscape through plant and animal occurrence patterns in legacy data. To this end, 190 botanical and 466 faunal assemblages from 309 sites are quantitatively analysed within a Bayesian framework, revealing a strong trend towards re-gionalisation in agricultural strategies during the early medieval period. In addition, these findings expose variations in agricultural techniques and dietary patterns across Roman settlements, shedding light on the extent to which Roman agricultural and economic frameworks persisted or changed during the early medieval transition, and the adaptive agricultural strategies adopted by farmers. The quantitative analyt-ical findings are also contextualised alongside wider historical sources, archaeological evidence, and cur-rent debate, allowing for a bottom-up understanding of the agricultural regimes in question. This work represents the first attempt to use temporally and geographically diverse bioarchaeological data to visualise the Italian agricultural landscape across the longue durée.
-
Archaeobotanical Data from the Italian Peninsula in the 1st Millennium CE
Roberto Ragno
Nov 2025
This dataset contains raw counts of archaeobotanical (macro-)remains from archaeological sites located in mainland Italy, dating from the 1st century BCE to the 11th century CE. The 195 carpological assemblages have been collected from a variety of published sources and reports and have been stored in a database with the relevant contextual and geographical metadata. The collection focuses on 40 plant taxa that were most commonly represented across the reviewed archaeobotanical reports, prioritising species with clear associations with agricultural activities and subsistence practices. This data may be used for statistical assessment of spatial and temporal differences in dietary patterns, agrarian strategies and crop distribution in the Italian peninsula between the Roman imperial period and the early medieval age.
-
Sheep and Goats Taxonomic Abundance Trends in 1st Millennium CE Southern Italy: Multilevel Bayesian Modelling of NISP Datasets
Roberto Ragno
Nov 2024
The 1st millennium CE represents a period of significant change in the agricultural landscape of southern Italy. Sheep and goats are among the most common faunal remains recovered from archaeological excavations of this period, but the contribution of these animals to the agricultural economy (particularly wool production) is often discussed through textual sources. This paper synthesises caprine taxonomic abundance trends using a Bayesian multilevel modelling approach that employs a beta-binomial distribution to address the problems of overdispersion and unequal assemblage/group sizes. Our models contribute directly to the problem by suggesting a period of change in livestock management practices around the 4th and 6th centuries CE, when the region’s shift to cereal farming appears to be accompanied by an increase in sheep and goat numbers.
-
Cereal Farming Practices in Italy during the 1st Millennium CE: An Integrated Approach to Regionality
Roberto Ragno
Nov 2023
This study investigates the cereal farming practices of the Italian peninsula during the first millennium CE, with a particular focus on the Early Middle Ages. Using non-parametric multivariate statistics and a dimensionality reduction algorithm (PERMANOVA; nMDS), this research presents and compares 177 archaeobotanical caryopses assemblages from three areas of Italy. The results showed that differences in cereal farming practices between Northern and Southern Italy were not statistically significant during the Roman period, but became significant during the early medieval period. The research suggests that after the collapse of the Roman Empire, northern peasants had more autonomy in selecting their crops, while southern farmers were more resistant to change. These findings challenge the prevailing assumption that early medieval peasants across the entire peninsula uniformly adopted cereal crops like millets and rye. Overall, the study sheds new light on the diversity of cereal farming practices in early medieval Italy and presents a comprehensive collection of cereal macroremains from mainland Italy.