Detecting Classic Maya Settlements with Lidar-Derived Relief Visualizations

Identification of Maya ruins covered by jungle using Sentinel-1

Archaeologists commonly use airborne LIDAR technology to produce 3D models of structures, even when obscured by a forest canopy. However, this technology has a high cost, both from the plane itself and from the processing of the LIDAR point cloud. Furthermore, this technique can only be used over limited regions. This paper proposes a technique that uses SAR satellite imagery to identify man-made structures hidden by a forest canopy. To do so, we exploit the Ascending and Descending passes of Sentinel-1 so that we obtain two images of the candidate site but from different sight directions. Because of cardinal effects, a large enough building will sign differently from the comparatively isotropic forest canopy it is obscured by. Practically, the technique is based on the ratio of backscattered intensity from these two illumination angles and is well adapted for large areas. The advantages and shortcomings are discussed for the specific case of Sentinel-1 SAR images over two Maya archaeological sites in Central America. Our analysis shows that SAR satellite imagery might provide a free, global-scale way of preselecting sites with large or tall structures to complement LIDAR technology.

Decline in seasonal predictability potentially destabilized Classic Maya societies

Classic Maya populations living in peri-urban states were highly dependent on seasonally distributed rainfall for reliable surplus crop yields. Despite intense study of the potential impact of decadal to centennial-scale climatic changes on the demise of Classic Maya sociopolitical institutions (750-950 CE), its direct importance remains debated. We provide a detailed analysis of a precisely dated speleothem record from Yok Balum cave, Belize, that reflects local hydroclimatic changes at seasonal scale over the past 1600 years. We find that the initial disintegration of Maya sociopolitical institutions and population decline occurred in the context of a pronounced decrease in the predictability of seasonal rainfall and severe drought between 700 and 800 CE. The failure of Classic Maya societies to successfully adapt to volatile seasonal rainfall dynamics likely contributed to gradual but widespread processes of sociopolitical disintegration. We propose that the complex abandonment of Classic Maya population centres was not solely driven by protracted drought but also aggravated by year-to-year decreases in rainfall predictability, potentially caused by a regional reduction in coherent Intertropical Convergence Zone-driven rainfall.

Ancient Lowland Maya neighborhoods: Average Nearest Neighbor analysis and kernel density models, environments, and urban scale

Many humans live in large, complex political centers, composed of multi-scalar communities including neighborhoods and districts. Both today and in the past, neighborhoods form a fundamental part of cities and are defined by their spatial, architectural, and material elements. Neighborhoods existed in ancient centers of various scales, and multiple methods have been employed to identify ancient neighborhoods in archaeological contexts. However, the use of different methods for neighborhood identification within the same spatiotemporal setting results in challenges for comparisons within and between ancient societies. Here, we focus on using a single method-combining Average Nearest Neighbor (ANN) and Kernel Density (KD) analyses of household groups-to identify potential neighborhoods based on clusters of households at 23 ancient centers across the Maya Lowlands. While a one-size-fits all model does not work for neighborhood identification everywhere, the ANN/KD method provides quantifiable data on the clustering of ancient households, which can be linked to environmental zones and urban scale. We found that centers in river valleys exhibited greater household clustering compared to centers in upland and escarpment environments. Settlement patterns on flat plains were more dispersed, with little discrete spatial clustering of households. Furthermore, we categorized the ancient Maya centers into discrete urban scales, finding that larger centers had greater variation in household spacing compared to medium-sized and smaller centers. Many larger political centers possess heterogeneity in household clustering between their civic-ceremonial cores, immediate hinterlands, and far peripheries. Smaller centers exhibit greater household clustering compared to larger ones. This paper quantitatively assesses household clustering among nearly two dozen centers across the Maya Lowlands, linking environment and urban scale to settlement patterns. The findings are applicable to ancient societies and modern cities alike; understanding how humans form multi-scalar social groupings, such as neighborhoods, is fundamental to human experience and social organization.

Towards Better Visualisation of Alpine Quaternary Landform Features on High-Resolution Digital Elevation Models

Alpine topography is formed by a complex series of geomorphological processes that result in a vast number of different landforms. The youngest and most diverse landforms are various Quaternary sedimentary bodies, each characterised by its unique landform features. The formation of Quaternary sedimentary bodies and their features derive from the dominant building sedimentary processes. In recent years, studies of Quaternary sedimentary bodies and processes have been greatly aided by the use of digital elevation models (DEMs) derived by airborne laser scanning (ALS). High-resolution DEMs allow detailed mapping of sedimentary bodies, detection of surface changes, and recognition of the building sedimentary processes. DEMs are often displayed as hillshaded reliefs, the most common visualisation technique, which suffers from the limitation of a single illumination source. As a result, features can be barely visible or even invisible to the viewer if they are parallel to the light source or hidden in the shadow. These limitations become challenging when representing landforms and subtle landscape features in a diverse alpine topography. In this study, we focus on eleven visualisations of Quaternary sedimentary bodies and their sedimentary and morphological features on a 0.5 m resolution DEM. We qualitatively compare analytical hillshading with a set of visualisation techniques contained in the Raster Visualisation Toolbox software, primarily hillshading from multiple directions RGB, 8-bit sky view factor and 8-bit slope. The aim is to determine which visualisation technique is best suited for visual recognition of sedimentary bodies and sedimentation processes in complex alpine landscapes. Detailed visual examination of previously documented Pleistocene moraine and lacustrine deposits, Holocene alluvial fans, scree deposits, debris flow and fluvial deposits on the created visualisations revealed several small-scale morphological and sedimentary features that were previously difficult or impossible to detect on analytical hillshading and aerial photographs. Hillshading from multiple directions resulted in a visualisation that could be universally applied across the mountainous and hilly terrains. In contrast, 8-bit sky view factor and 8-bit slope visualisations created better visibility and facilitated interpretation of subtle and small-scale (less than ten metres) sedimentary and morphological features.

Airborne Lidar Survey, Density-Based Clustering, and Ancient Maya Settlement in the Upper Usumacinta River Region of Mexico and Guatemala

We present results from the archaeological analysis of 331 km2 of high-resolution airborne lidar data collected in the Upper Usumacinta River basin of Mexico and Guatemala. Multiple visualizations of the DEM and multi-spectral data from four lidar transects crossing the Classic period (AD 350–900) Maya kingdoms centered on the sites of Piedras Negras, La Mar, and Lacanja Tzeltal permitted the identification of ancient settlement and associated features of agricultural infrastructure. HDBSCAN (hierarchical density-based clustering of applications with noise) cluster analysis was applied to the distribution of ancient structures to define urban, peri-urban, sub-urban, and rural settlement zones. Interpretations of these remotely sensed data are informed by decades of ground-based archaeological survey and excavations, as well as a rich historical record drawn from inscribed stone monuments. Our results demonstrate that these neighboring kingdoms in three adjacent valleys exhibit divergent patterns of structure clustering and low-density urbanism, distributions of agricultural infrastructure, and economic practices during the Classic period. Beyond meeting basic subsistence needs, agricultural production in multiple areas permitted surpluses likely for the purposes of tribute, taxation, and marketing. More broadly, this research highlights the strengths of HDBSCAN to the archaeological study of settlement distributions when compared to more commonly applied methods of density-based cluster analysis.

Assessing Classic Maya multi-scalar household inequality in southern Belize

Inequality is present to varying degrees in all human societies, pre-modern and contemporary. For archaeological contexts, variation in house size reflects differences in labor investments and serves as a robust means to assess wealth across populations small and large. The Gini coefficient, which measures the degree of concentration in the distribution of units within a population, has been employed as a standardized metric to evaluate the extent of inequality. Here, we employ Gini coefficients to assess wealth inequality at four nested socio-spatial scales-the micro-region, the polity, the district, and the neighborhood-at two medium size, peripheral Classic Maya polities located in southern Belize. We then compare our findings to Gini coefficients for other Classic Maya polities in the Maya heartland and to contemporaneous polities across Mesoamerica. We see the patterning of wealth inequality across the polities as a consequence of variable access to networks of exchange. Different forms of governance played a role in the degree of wealth inequality in Mesoamerica. More autocratic Classic Maya polities, where principals exercised degrees of control over exclusionary exchange networks, maintained high degrees of wealth inequality compared to most other Mesoamerican states, which generally are characterized by more collective forms of governance. We examine how household wealth inequality was reproduced at peripheral Classic Maya polities, and illustrate that economic inequity trickled down to local socio-spatial units in this prehispanic context.

Documentation of Archaeology-Specific Workflow for Airborne LiDAR Data Processing

Airborne LiDAR is a widely accepted tool for archaeological prospection. Over the last decade an archaeology-specific data processing workflow has been evolving, ranging from raw data acquisition and processing, point cloud processing and product derivation to archaeological interpretation, dissemination and archiving. Currently, though, there is no agreement on the specific steps or terminology. This workflow is an interpretative knowledge production process that must be documented as such to ensure the intellectual transparency and accountability required for evidence-based archaeological interpretation. However, this is rarely the case, and there are no accepted schemas, let alone standards, to do so. As a result, there is a risk that the data processing steps of the workflow will be accepted as a black box process and its results as “hard data”. The first step in documenting a scientific process is to define it. Therefore, this paper provides a critical review of existing archaeology-specific workflows for airborne LiDAR-derived topographic data processing, resulting in an 18-step workflow with consistent terminology. Its novelty and significance lies in the fact that the existing comprehensive studies are outdated and the newer ones focus on selected aspects of the workflow. Based on the updated workflow, a good practice example for its documentation is presented.