skip to main content


Title: Four millennia of geomorphic change and human settlement in the lower Usumacinta–Grijalva River Basin, Mexico

The lower Usumacinta–Grijalva River Basin contains one of the richest biodiversity landscapes of the Maya region. Our research is based on (1) an integrative literature review of the geomorphological and archaeological papers published about the lower Usumacinta–Grijalva River Basin and (2) topographic analysis of digital elevation models using a geographical information system to explore the relationship between past human settlement and landscape accessibility along the coastal plain of Tabasco. This work provides a new synthesis of previous research and proposes new models for the geomorphic evolution of the lower Usumacinta–Grijalva River Basin in the context of four millennia of human land use and settlement. For the evolution of the strand-plain of the Usumacinta and Grijalva rivers, there are two published geochronological models that provide different chronologies. We discuss here how both geochronological models encompass Pre-Columbian human settlement in the delta. Interestingly, we notice that one of them overlaps a possible high-magnitude flood event (or events) that drove large geomorphic change around 750 CE (1200 BP), with implications for settlement patterns and chronology. Based on topographical analysis of the eastern-distal sector of the Usumacinta–Grijalva delta, we propose a new model for the evolution of this area with implications for the human occupation during the Mesoamerican Terminal Classic and Early Postclassic on the delta. As one of the main conclusions, we propose that the Pom–Atasta water bodies predate much of the Usumacinta–Grijalva delta and the most recent phase of delta building overlays the original lagoon barriers, resulting in a geomorphic setting more attractive to local human occupation after the Terminal Classic period. According to one of the geochronological models of the delta, this dates to ca. 900 CE, preceding the establishment of nearby settlements such as Atasta.

 
more » « less
NSF-PAR ID:
10397941
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
SAGE Publications
Date Published:
Journal Name:
Progress in Physical Geography: Earth and Environment
Volume:
47
Issue:
2
ISSN:
0309-1333
Page Range / eLocation ID:
p. 227-248
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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. 
    more » « less
  2. Abstract River deltas all over the world are sinking beneath sea-level rise, causing significant threats to natural and social systems. This is due to the combined effects of anthropogenic changes to sediment supply and river flow, subsidence, and sea-level rise, posing an immediate threat to the 500–1,000 million residents, many in megacities that live on deltaic coasts. The Mississippi River Deltaic Plain (MRDP) provides examples for many of the functions and feedbacks, regarding how human river management has impacted source-sink processes in coastal deltaic basins, resulting in human settlements more at risk to coastal storms. The survival of human settlement on the MRDP is arguably coupled to a shifting mass balance between a deltaic landscape occupied by either land built by the Mississippi River or water occupied by the Gulf of Mexico. We developed an approach to compare 50 % L:W isopleths (L:W is ratio of land to water) across the Atchafalaya and Terrebonne Basins to test landscape behavior over the last six decades to measure delta instability in coastal deltaic basins as a function of reduced sediment supply from river flooding. The Atchafalaya Basin, with continued sediment delivery, compared to Terrebonne Basin, with reduced river inputs, allow us to test assumptions of how coastal deltaic basins respond to river management over the last 75 years by analyzing landward migration rate of 50 % L:W isopleths between 1932 and 2010. The average landward migration for Terrebonne Basin was nearly 17,000 m (17 km) compared to only 22 m in Atchafalaya Basin over the last 78 years (p\0.001), resulting in migration rates of 218 m/year (0.22 km/year) and\0.5 m/year, respectively. In addition, freshwater vegetation expanded in Atchafalaya Basin since 1949 compared to migration of intermediate and brackish marshes landward in the Terrebonne Basin. Changes in salt marsh vegetation patterns were very distinct in these two basins with gain of 25 % in the Terrebonne Basin compared to 90 % decrease in the Atchafalaya Basin since 1949. These shifts in vegetation types as L:W ratio decreases with reduced sediment input and increase in salinity also coincide with an increase in wind fetch in Terrebonne Bay. In the upper Terrebonne Bay, where the largest landward migration of the 50 % L:W ratio isopleth occurred, we estimate that the wave power has increased by 50–100 % from 1932 to 2010, as the bathymetric and topographic conditions changed, and increase in maximum storm-surge height also increased owing to the landward migration of the L:W ratio isopleth. We argue that this balance of land relative to water in this delta provides a much clearer understanding of increased flood risk from tropical cyclones rather than just estimates of areal land loss. We describe how coastal deltaic basins of the MRDP can be used as experimental landscapes to provide insights into how varying degrees of sediment delivery to coastal deltaic floodplains change flooding risks of a sinking delta using landward migrations of 50 % L:W isopleths. The nonlinear response of migrating L:W isopleths as wind fetch increases is a critical feedback effect that should influence human river-management decisions in deltaic coast. Changes in land area alone do not capture how corresponding landscape degradation and increased water area can lead to exponential increase in flood risk to human populations in low-lying coastal regions. Reduced land formation in coastal deltaic basins (measured by changes in the land:water ratio) can contribute significantly to increasing flood risks by removing the negative feedback of wetlands on wave and storm-surge that occur during extreme weather events. Increased flood risks will promote population migration as human risks associated with living in a deltaic landscape increase, as land is submerged and coastal inundation threats rise. These system linkages in dynamic deltaic coasts define a balance of river management and human settlement dependent on a certain level of land area within coastal deltaic basins (L). 
    more » « less
  3. null (Ed.)
    Abstract. The morphology of deltas is determined by the spatial extent and variability of the geomorphic processes that shape them. While in some cases resilient, deltas are increasingly threatened by natural and anthropogenic forces, such as sea level rise and land use change, which can drastically alter the rates and patterns of sediment transport. Quantifying process patterns can improve our predictive understanding of how different zones within delta systems will respond to future change. Available remotely sensed imagery can help, but appropriate tools are needed for pattern extraction and analysis. We present a method for extracting information about the nature and spatial extent of active geomorphic processes across deltas with 10 parameters quantifying the geometry of each of 1239 islands and the channels around them using machine learning. The method consists of a two-step unsupervised machine learning algorithm that clusters islands into spatially continuous zones based on the 10 morphological metrics extracted from remotely sensed imagery. By applying this method to the Ganges–Brahmaputra–Meghna Delta, we find that the system can be divided into six major zones. Classification results show that active fluvial island construction and bar migration processes are limited to relatively narrow zones along the main Ganges River and Brahmaputra and Meghna corridors, whereas zones in the mature upper delta plain with smaller fluvial distributary channels stand out as their own morphometric class. The classification also shows good correspondence with known gradients in the influence of tidal energy with distinct classes for islands in the backwater zone and in the purely tidally controlled region of the delta. Islands at the delta front under the mixed influence of tides, fluvial–estuarine construction, and local wave reworking have their own characteristic shape and channel configuration. The method is not able to distinguish between islands with embankments (polders) and natural islands in the nearby mangrove forest (Sundarbans), suggesting that human modifications have not yet altered the gross geometry of the islands beyond their previous “natural” morphology or that the input data (time, resolution) used in this study are preventing the identification of a human signature. These results demonstrate that machine learning and remotely sensed imagery are useful tools for identifying the spatial patterns of geomorphic processes across delta systems. 
    more » « less
  4. The morphology of deltas is determined by the spatial extent and variability of the geomorphic processes that shape them. While in some cases resilient, deltas are increasingly threatened by natural and anthropogenic forces, such as sea level rise and land use change, which can drastically alter the rates and patterns of sediment transport. Quantifying process patterns can improve our predictive understanding of how different zones within delta systems will respond to future change. Available remotely sensed imagery can help but appropriate tools are needed for pattern extraction and analysis. We present a method for extracting information about the nature and spatial extent of active geomorphic processes across deltas from ten parameters quantifying the geometry of each of 1239 islands and the channels around them using machine learning. The method consists of a two-step unsupervised machine learning algorithm that clusters islands into spatially continuous zones based on the ten morphological metrics extracted from remotely sensed imagery. By applying this method to the Ganges–Brahmaputra–Meghna Delta, we find that the system can be divided into six major zones. Classification results show that active fluvial island construction and bar migration processes are limited to relatively narrow zones along the main Ganges River and Brahmaputra and Meghna corridors, whereas zones in the mature upper delta plain, with smaller fluvial distributary channels stand out as their own morphometric class. The classification also shows good correspondence with known gradients in the influence of tidal energy with distinct classes for islands in the backwater zone and in the purely tidally-controlled region of the delta. Islands at the delta front, under the mixed influence of tides, fluvial-estuarine construction, and local wave reworking have their own characteristic shape and channel configuration. The method does not distinguish between islands with embankments (polders) and natural islands in the nearby mangrove forest (Sundarbans), suggesting that human modifications have not yet altered the gross geometry of the islands beyond their previous natural morphology. These results demonstrate that machine learning and remotely sensed imagery are useful tools for identifying the spatial patterns of geomorphic processes across delta systems. 
    more » « less
  5. Abstract

    Rivers that traverse the terrestrial‐marine interface may have lower reaches that are influenced by both terrestrial and marine processes. However, only a handful of studies have focused on how the interactions of fluvial and tidal processes translate to channel geomorphology, and those are largely from delta/distributary systems. Here we quantify channel properties along the fluvial‐tidal transition reach of a coastal plain river and provide insight into their origins. The study site is a 47 km long tidal, single‐thread freshwater section of a river at 29 to 76 river kilometers inland of the estuary mouth, upstream of the delta/distributary system, and with average riverbed slope of 10−4. Results show that a tidal wave approaching the study reach loses 15%–17% of its incident energy (per horizontal area) per kilometer of channel, and at 51 km upstream of the mouth the incident energy is reduced to <1%. Also, at or near 51 km we observed breaks in along‐channel trends of channel cross‐section geometry, bed grain size, sinuosity, channel bed and water surface slopes. We propose that fluvial‐tidal flow processes and corresponding geomorphic feedbacks are apparent as abrupt changes in channel properties that highlight the influence of tides, and these discontinuities may be endemic to fluvial‐tidal transition zones in general. How these transition reaches self‐adjust in response to climate change remains largely unexplored but these reaches are likely to become important geomorphic hotspots.

     
    more » « less