Abstract
The Atacama Desert is the driest and oldest desert on Earth. Despite the abundance evidence for long-term landscape stability, there are subtle signs of localised fluvial erosion and deposition since the onset of hyperaridity in the rock record. In the dry core of the Atacama Desert, pluvial episodes allowed antecedent drainage to incise into uplifting fault scarps, which in turn generated sinuous to meandering channels. Incision of ancient alluvial fan surfaces occurred during intermittent fluvial periods, albeit without signs of surface erosion. Fluvial incision during predominantly hyperarid climate periods is evident from these channels in unconsolidated alluvium. The absence of dense vegetation to provide bank stability and strength led us to investigate the potential role of regionally ubiquitous CaSO4-rich surface cover. This has enabled the preservation of Miocene surfaces and we hypothesize that it provided the required bank stability by adding strength to the upper decimetre to meter of incised alluvium to allow high sinuosity of stream channels to form during pluvial episodes in the Quaternary.
Similar content being viewed by others
Introduction
Hyperarid conditions have prevailed in the dry core of the Atacama Desert since the Early Miocene1,2,3,4. There are several inter-related reasons for the aridification; (1) its position within the area of the subtropical high pressure belt5, (2) upwelling of the cold Humboldt current in conjunction with the evolution and existence of a persistent atmospheric inversion layer trapping moist Pacific air below 1000 m a. s. l.6,7, (3) the orogenic rainshadow effect of the Andes to the east5,6,7,8, as well as (4) the continental rainshadow effect of moisture originating from the Atlantic. The timing and location of the onset of aridity, and its evolution are still a matter of debate. For instance, there are strong indications that parts of the hyperarid Coastal Cordillera and Central Depression (19–23°S) experienced the earliest onset of hyperaridity during the Early Miocene1,2,4. Long-term hyperarid conditions and the absence of significant precipitation enabled the evolution of thick surface covers of evaporites and dust, sourced from atmospheric deposition9,10,11. Relict surfaces are covered by decimetre to meter thick calcium-sulphate (CaSO4) rich soils or colluvium11 and are unable to sustain significant macroscopic vegetation. However, there is evidence that the predominantly hyperarid conditions were repeatedly interrupted by wetter—though still (hyper-) arid—periods2 which were capable of initiating fluvial incision and erosion2,3,12.
Fluvial systems in arid-hyperarid environments are characterized by large alluvial fans and small incised channels13, as well as braided channel networks, where ephemeral drainage dominates13,14. They form in response to rare, but intense precipitation events15. Ephemeral streams are often unstable and their channel form quickly alter in response to precipitation events, which cause a dynamic expression by channel migration, low temporal stability, high gradients and mostly a predominance of bedload15,16. Assuming that pluvial periods in the Atacama Desert are ephemeral, the observation of single thread sinuous to meandering channels in the hyperarid core of the Atacama Desert is surprising (Fig. 1). Sinuous to meandering channel forms are reported in modern poorly vegetated arid environments16,17, however, only rarely in ancient hyperarid environments16,17,18, such as the core of the Atacama Desert. This suggests that they are either fossil, formed during wetter periods prior, or the development of high sinuosity fluvial systems was enabled by other processes.
In this study we bracket the timing of alluvial fan abandonment and channel incision in the Atacama Desert using cosmogenic nuclide exposure dating of quartz clasts (10Be, 26Al and 21Ne) and tephrochronology. This frames the paleoenvironmental conditions during channel formation and provides causal mechanism for the development of sinuous to meandering stream channels in a hyperarid environment. We show that fluvial incision was coeval with long term (hyper-) aridity and that the ubiquitous CaSO4-rich soil and crust helped to preserve landscapes and promoted high sinuosity of fluvial channels.
Alluvial fans and meander in the hyperarid Atacama Desert
Seemingly inactive highly sinuous to meandering stream channels are etched into the landscape in the eastern part of the northern Costal Cordillera, 30 km NE of Iquique in the Pampa de la Perdiz (Fig. 1). The E-W channels cut into the Alto Hospicio gravels (Fig. 1). These unconsolidated and partially salt-cemented upper Oligocene to Pliocene gravels, sands, silts and clays were deposited as sheet-, debris and locally sand and mud flows19,20. The absence of large outcrops made it difficult to study the internal composition and structure of the Alto Hospicio gravels. However, they are likely comparable to active alluvial fan system in the Precordillera and Central Depression to the East21,22. Alluvial fan surfaces dip moderately (mean ~ 4°) to the west and consist of erosional products of the elevated (max. ~ 1700 m) Huara Intrusive Complex (HIC, 140–136 Ma19) to the east. All alluvial fan surfaces are covered by a several decimetres to meter thick CaSO4-rich soil and crust sourced from atmospheric deposition9,10,23. Compressional tectonic forces caused a reverse reactivation of former normal faults during the Miocene along N–S and E–W running fault scarps19. Tectonic activity resulted in the beheading of the large depositional system to the west of the N-S-trending scarp. The fluvial channels were the final spilling pathways across the scarp, before they were cut from their source area (Figs. 1, 2). Fluvial channels are incised up to ~ 14 m into the crest of the thrust fault, indicating that initially incision kept pace with uplift. The tectonically truncated drainages terminate in endorheic basins filled by pan deposits to the east of the scarp24. All truncated channels crossing the scarp display distinct meanders with a sinuosity between 1.02 and 1.55 (see Supplementary Dataset), with highest values near the crest of the fault. The rims of the channel slopes on the outside of bends are commonly steeper than adjacent channel slopes due to induration by the regionally pervasive CaSO4 cover. The fluvial drainage system in the alluvial deposits can be best explained as confined sinuous to meandering channel. Although the fossil channels indicate a rather immature stage of evolution (i.e. floodplain did not form and the number of cut-offs is limited), they preserve original depositional fluvial features, such as mid-channel bars in straight sections (Fig. 2D). This is best observed at sites that are distal to the crest of the fault. Close to the crest, tectonic deformation and slumping of the channel slopes mostly obscure the original channel morphology. Occasionally upper CaSO4-rich layers protrude at channel slopes and cantilever failures as a consequence of the undermining of more erodible sediments beneath, indicating a significant difference in the strength of both materials. Fluvial channel bottoms are in some places covered by pristine tephra (Fig. 2D) which are covered by a friable CaSO4-rich cover. Channel slopes are close to the angle of repose of the Alto Hospicio gravels and/or their colluvium. No tributaries originate from the adjacent flat alluvial fan surfaces, indicating that runoff was dominantly sourced from the catchments in HIC (Fig. 1). The latter indicates that precipitation events during fluvial episodes, to the same extent affecting the entire area, did not modify fluvial surfaces on directly adjacent fan surfaces.
Miocene surfaces and Plio-Pleistocene fluvial incision
Quartz clasts were sampled from alluvial fan surfaces and fluvial channel beds, as were two tephra layers from one channel bed, to determine the chronology of landscape evolution and constrain the environment under which sinuous to meandering channels evolved (Figs. 1, 2). The oldest ages (21Ne exposure ages 22–25 Ma) are from the alluvial fan surfaces to the east of the thrust fold (Fig. 3B), and point to long-term surface stability. This is in line with the notion that the region experienced predominantly (hyper-) arid conditions since the early Miocene1,2,4,25. Younger exposure age clusters are consistent with localized debris flow deposition onto the alluvial fans during the Mid-Miocene and a cessation of deposition after 5–6 Ma (Fig. 3), with the latter providing a maximum age-limit for the tectonic truncation of fan deposits from their source area. 10Be concentrations of fan-surface materials are saturated, confirming that exposure lasted longer than 4–5 Myr. Based on the ‘born at the surface model’ by Wells et al.26, an accretionary mantle of ‘soil-modified dust’ can build up beneath clasts previously exhumed to, or deposited onto, the surface. Therefore, the accumulation of atmospheric deposits took place after the initial exposure of the clasts. Any significant fluvial event will remove the clasts and erode the accretionary gypsic soil beneath. Long-term persistent aridity and the absence of any significant surface activity promoted the accumulation of CaSO4 and other clastic aeolian material since at least the deposition of the clasts and caused the preservation of a thick CaSO4 surface cover from atmospheric deposition9,10,23.
Tectonic ponding of the drainage to the east of the scarp presumably assisted in creating hydrologic conditions favourable for fluvial erosion (continuous flow, smoothing of peak discharges) further downstream before further tectonic uplift eventually beheaded these channels from their catchment source, creating fully closed endorheic basins. The tectonically created basins might acted as settling ponds for sediment transported during rare precipitation events, allowing bedload-free and sediment-poor water to overflow the tectonic barrier and subsequently cause incision. Rare and localised storm events in the Atacama Desert generated precipitation rates of 24–30 mm/h27. Based on disc infiltrometer measurements, Pfeiffer et al.27 estimated that precipitation has to exceed on average 78 mm/h for hillslopes to initiate infiltration/saturation-excess overland flow. Based on the recent occurrence of gypsic soils in the Atacama Desert mean annual precipitation probably did not exceed ~ 30 mm/year9,23. A sediment record from the southernmost endorheic clay pan confirms that fluvial deposition was restricted to these closed basins24. Besides removing the bulk of sediments from potential debris flows, the basins likely served as hydrologic buffers by flattening peak discharge and temporarily mimicking continuous flow during spilling events. Tectonic steepening of the alluvial fan slope near the crest may have promoted the initiation of, or further amplified28 high sinuosity29,30, ensuring that the flow remained confined.
Fluvial channels incised into the antecedent alluvial fan surface. This was at first capable of keeping pace with tectonic uplift but was later outpaced. 10Be and 26Al exposure ages from individual channel bed clasts and correlative dating of the channel tephra deposits within indicates that the drainages were terminally truncated at 640–800 ka in the southernmost catchment and ~ 185 ka in the northernmost catchment (Fig. 3). Based on the age-depth model of the southernmost clay pan climate record24 and the thickness of accumulated sediments it can be estimated that this pan received its current endorheic fill in the last 650–850 kyr. From the age of the alluvial fan surfaces (termination of fan-deposition) and channel deposits (termination of sediment transport, full cut-off), the sinuous to meandering channels incised sometime between the latest Miocene and late Quaternary. During this time, climatic conditions were arid to hyperarid on the Andean fore slope to the east of the study area3,4. Since the areas to the west of the Andean fore slope are generally drier (e.g.1,2,3,4), we also assume the study area was arid to hyperarid.
During this time, intermittent locally-sourced precipitation events must have been capable of incising the channels into the flat alluvial fans without modifying the surface. Vegetation cover is widely considered as one major prerequisite for adding strength to channel banks and slopes, so permitting high sinuosity of channels to develop in alluvial systems31. However, the absence of vegetation in the region requires another mechanism or process to enable the formation of these channel planforms by adding strength to channel banks and slopes.
CaSO4-rich soil as landscape forming agent
Preventing surface erosion
The old exposure ages of fan-surface clasts bear witness to hyperaridity and stability of the substrate since the Early Miocene. These conditions are a prerequisite for the accumulation and preservation of CaSO4-rich soils (Fig. 4). Long-term accumulation, preservation and modification of atmospheric CaSO4 deposition has led to cemented soils that cover much of the Atacama Desert9,10,11,23,32,33,34, that can be up to several meters thick11,23,34. Based on field observations, we ascertain a thickness of at least several decimetres to meters for our field site. Pedogenic processes which cause induration and cementation, include a combination of hydromorphic and illuvial processes34,35. These can be accentuated by diagenetic and surface exhumation36 and modified by the variable interaction of meteoric water or moisture; causing dissolution, redistribution and reprecipitation in the CaSO4-rich soil. The latter causes growth of CaSO4 and other soluble salts, leading to the cementation of the soil and underlying alluvial sediments (Fig. 4,27,34,35,36). These transformations are facilitated by moisture originating as rain, fog or dew. Dissolution and reprecipitation of CaSO4 phases by meteoric waters is verifiable through a poikilitic texture23,34,36. Surface dissolution and subsurface reprecipitation creates a thick, clast poor layer of almost purely CaSO4-rich material11,23,34,35,36. Anhydrite may be formed in the presence of highly saline solutions10,32,37 or due to dehydration of primary gypsum37. The uppermost soil layer, which is covered by a friable crust 11, usually remains powdery (‘chusca’ sensu Ericksen11), whereas concretions occur at depth (Fig. 4,9,11). Quantitative determination of the atmospheric deposition rate and soil inflation are difficult and rare from the Atacama Desert. Wang, et al.10 derived rates of ~ 35 cm/Ma soil accumulation based on meteoric 10Be data, which is in line with data from Placzek, et al.38. Soil accumulation rates using tephra deposits indicated also higher deposition rates of up 1 m/Ma32. The powdery layer promotes infiltration and prevents overland flow of any precipitation27,39,40. If the precipitation rate is below the infiltration capacity of the soil, overland flow will be prevented in favour of complete infiltration15. The infiltrated moisture drives dissolution and reprecipitation at depth, causing the induration of subsurface sediments27,34,39,40. Jordan et al.40 found indications that flat angle hillslopes and surfaces (< 7°, based on 90 m SRTM) experienced no measurable surface modifications, thus implying surficial covers of CaSO4 rich soils may hinder surface activity by absorbing rain. Precipitation falling the higher elevated catchments of the HIC were able to produce concentrated fluvial flows. Our study illustrates the role of this process in buffering the local erosional capacity of precipitation over millions of years, allowing the preservation of Miocene surfaces in direct proximity to Quaternary fluvial incisions.
Adding strength: modulating fluvial incision
The development of sinuosity of a fluvial system, depends on several factors. Among them arguably the strength of parent rock or alluvium to resist fluvial erosion (flow resistance) dominates41. Strength can be provided or enhanced by sediment cohesiveness, bedrock strength or dense vegetation41,42,43,44. Riparian vegetation or a dense vegetation cover and bedrock strength can be clearly excluded here. Alternative mechanisms that may increase flow resistance in the absence of vegetation are: (1) permafrost45, (2) the cohesion of fine-grained sediments18,44,46,47,48, and (3) chemical cementation synchronous with channel flow18,47,49. These mechanisms may have promoted high sinuosity and meandering in arid to hyperarid environments in pre-vegetational streams of early Phanerozoic and Precambrian, and as potential analogues for sinuous to meandering channels on Mars45,46,47,48,50,51,52,53. None of these possibilities applies directly to this study. The alluvium of the Altos de Hospicio gravels might be the source of fine-grained sediments, however, strong prevailing winds presumably did not permit a significant deposition within the channels, which could be responsible for increasing strength and bank stability by cohesion. Fluvial systems with high sinuosity typically range between two extremes, bedrock/ice restricted channels with low amounts of alluvium involved, called incised or entrenched or channels fully developed in an alluvial setting (54, and citations therein), previously deposited in or near them (i.e. alluvial meander). Confined sinuous to meandering channels, as recognized in our study area, are therefore restricted in their lateral growth (forming/creating larger floodplains, cut-offs, etc.) by increased flow resistance of the sediment41,54. Hillslope failures by toe slope erosion and undercutting indicate that this increased flow resistance is only restricted to the upper part of the alluvium that is indurated by CaSO4. The observed channel forms therefore might be inherited, created during initial incision. A lack of discharge has probably hindered the migration of loops since the initial incision.
The case described here points to an additional mechanism that may promote higher sinuosity in a non-vegetated environment by adding strength to the sediments prior to incision in arid to hyperarid environments and dominance of atmospheric deposition. We posit that the regionally ubiquitous CaSO4-rich soils (Figs. 3, 4) that evolved from atmospheric deposition9,10,11,32 and subsequent pedogenic modification (e.g.9,10,37) provide the necessary strength to form high sinuosity of fluvial channels, restricting lateral erosion. This cementation of the upper sediments might mimic bedrock strength, ground-ice, or vegetation, causing higher sinuosity, although being confined, restricting/preventing the evolution to a fully meandering system in the first place. Carbonate as cement in this kind of hyperarid soil can be excluded as dominant ingredient, as ‘its formation is not chemically favoured in this hyperarid setting’32. CaSO4-wedges9,11,55 (‘sand wedges’56) and pressure regimes inside the soil due to dissolution and recrystallization may cause fracturing, similar to ignimbrites/bedrock, which promotes incision and channel formation, reflecting similar bedrock characteristics. Due to limited outcrop situation, we cannot fully exclude that cohesion or partly cementation of alluvial deposits additionally caused increased shear strength. However, the occurrence of a protruding CaSO4-rich soil and crust at the channels’ rims points CaSO4 cementation as the dominating factor for increased shear strength. We assume that the sinuous to meandering channel planforms were created during the initiation of fluvial incision into the strengthened upper meter or so of alluvium, and that subsequent fluvial flow was largely confined by the initial planform. Rapid channel incision into the CaSO4 cemented alluvium may have initially shaped the channel pattern with a subsequent inherited sinuosity. Fluvial periods presumably were not frequent or long enough to enable the evolution of floodplains, or the antecedent drainage was disconnected from the catchments by tectonic activity. Observed cantilever failures caused by the strength difference of the upper cemented alluvium and the weakly cemented and unconsolidated alluvium of the Altos de Hospicio gravels19 underlines the effect that pedogenic CaSO4 modification of sediments can have in the Atacama Desert. The stable pedogenically indurated substrate apparently provides the required sediment strength to maintain high sinuosity of channel forms. Capped duricrust are known to have a protective effect, resistant enough to preserve ancient landscapes15,35. The observation of sulphate-capped ridges from the Atacama Desert, resembling inverted fluvial stream channel indicate that sulphate cementation is able to increase the erosion resistance against fluvial and aeolian processes57. Studies of gypcrete in Egypt36 supports the protective nature of gypsum crusts against erosion and weathering, consolidating surface layers by cementation. The features observed here, such as bank failure/slumping by toe-slope erosion/cantilever failures of undercut banks, are equivalent to those commonly observed in environments with vegetated channel-banks58 (Fig. 2G,H).
The notion that antecedent encrustation/cementation, i.e., independent from and prior to channel flow, is a potential mechanism for surface induration and cementation that preserves landscape and adds strength to sediments, setting the foundation to enable high sinuosity in previously unconsolidated sediments in long-term predominantly hyperarid environments. Enduring atmospheric deposition and subsequent pedogenic cementation, i.e. the evolution of duricrusts, might have relevance to other arid to hyperarid climates and for the geomorphology of Mars, where similar channel forms are evident45,46,47,48,50,51,52,53. The Atacama Desert serves as a natural laboratory for potential surface processes on Mars due to the broadly similar climatic conditions and geochemical/mineralogical compositions32,57,59. Sulphates (albeit MgSO460) are abundant on Martian surface60,61,62,63,64,65,66 and atmospheric deposition of sulphates from volcanic origin67,68 may give rise to surface encrustations that are similar to those found in the Atacama9,10,11,27. Mangold et al.69 have already indicated that Martian surfaces with polyhydrated sulphates in Candor Chasma are less eroded than the surfaces with monohydrated sulphates, which might be attributed to the formation of a duricrusts69,70.
The rates of and conditions for CaSO4 induration of sediments, and the build-up of CaSO4-rich soils, in the Atacama Desert will be the focus of future research. In particular, the conditions for cementation (state and volume of moisture involved, infiltration depth, and resulting cementation thickness) and their potential effects increasing sediment strength will be a focus. Experiments and models have to be conducted to further constrain this hypothesis and its effect on fluvial geomorphology in hyperarid settings, which is beyond the scope of this manuscript.
Methods
We used geological field observations and a high-resolution Digital Elevation Model (DEM) of the study area, based on Pléiades 1B data (see Supplementary Dataset), to create swath profiles. Optical satellite imagery (Pléiades 1B Multispectral Image) was additionally used to identify geomorphologic features, such as the extent of alluvial fans and channels. Tephra layers in the southern main channel were geochemically characterized along with already dated and published tephra deposits from the local vicinity. Details about the lab procedure are outlined in the Supplementary Dataset. The chronology of the depositional systems (alluvial fans) and fluvial incision is reconstructed using exposure dating utilizing cosmogenic nuclides (10Be, 26Al, 21Ne). Exposure ages were determined for quartz clasts sampled during expeditions in 2014, 2015, and 2017. Clasts were sampled from flat fan surfaces (HU15-18, HU17-01, HU15-15), incised fluvial channels (HU14-05, HU15-16, HU17-09, HU15-14), and a channel terrace (HU17-04). Surface quartz clasts exhibit a reddish-brown desert varnish and are mostly angular. In cases of ‘kernsprung’, indicated by a localized cluster (usually < 2 m diameter) of quartz fragments, only one fragment per cluster was sampled and no other samples were collected within a 5 m radius. Quartz clasts from channels are partly rounded and have only a light desert varnish. The size of sampled clasts ranges between 2 and 5 cm. The abundance of quartz on the surfaces and channel beds is very low. Attempts to retrieve quartz from the subsurface for depth-profiling were impractical due to the CaSO4 indurated crust. To identify and correct for possible pre-exposure and exhumation we utilize the multiple clast approach2,71,72, using 5–7 clasts per sampling site. Further details to sampling sites and sampling details are given in the Supplementary Dataset. Quartz samples were crushed, sieved to retain the 250–710 µm fraction and etched several times in a dilute HF-HNO3 mixture73. Splits of the etched samples were used for 10Be, 26Al, and 21Ne analysis. 10Be and 26Al sample preparation and measurement followed the single stacked column approach detailed in74. Samples were measured at the CologneAMS75. Ne isotope analyses were conducted at SUERC (UK)76 and in Cologne (Germany)77. Further details about sampling preparation and measurement procedure are outlined in the Supplementary Dataset. Exposure ages were calculated using the LSDn scaling scheme of Lifton et al.78, as implemented in version 3 of ‘the online calculators formerly known as the CRONUS-Earth online calculators’ https://hess.ess.washington.edu/math/v3/v3_age_in.html described in79; see also Supplement. We assume a time integrated linear uplift of 40 m/Ma since 23 Ma, and report uplift-corrected ages accordingly (see Supplement).
Data availability
All data generated or analysed during this study are included in this published article (and its Supplementary Information files).
References
Dunai, T. J., Lopez, G. A. G. & Juez-Larre, J. Oligocene-Miocene age of aridity in the Atacama Desert revealed by exposure dating of erosion-sensitive landforms. Geology 33, 321–324. https://doi.org/10.1130/g21184.1 (2005).
Ritter, B. et al. Neogene fluvial landscape evolution in the hyperarid core of the Atacama Desert. Sci. Rep. 8, 13952. https://doi.org/10.1038/s41598-018-32339-9 (2018).
Jordan, T. E., Kirk-Lawlor, N. E., Blanco, P. N., Rech, J. A. & Cosentino, N. J. Landscape modification in response to repeated onset of hyperarid paleoclimate states since 14 Ma, Atacama Desert Chile. Geol. Soc. Am. Bull. 14, 15. https://doi.org/10.1130/b30978.1 (2014).
Evenstar, L. et al. Geomorphology on geologic timescales: Evolution of the late Cenozoic Pacific paleosurface in Northern Chile and Southern Peru. Earth-Sci. Rev. 171, 1–27 (2017).
Houston, J. Variability of precipitation in the Atacama Desert: Its causes and hydrological impact. Int. J. Climatol. 26, 2181–2198. https://doi.org/10.1002/joc.1359 (2006).
Garreaud, R. D., Molina, A. & Farias, M. Andean uplift, ocean cooling and Atacama hyperaridity: A climate modeling perspective. Earth Planet. Sci. Lett. 292, 39–50. https://doi.org/10.1016/j.epsl.2010.01.017 (2010).
Garreaud, R. D., Vuille, M., Compagnucci, R. & Marengo, J. Present-day South American climate. Palaeogeogr. Palaeoclimatol. Palaeoecol. 281, 180–195. https://doi.org/10.1016/j.palaeo.2007.10.032 (2009).
Houston, J. & Hartley, A. J. The Central Andean west-slope rainshadow and its potential contribution to the origin of hyper-aridity in the Atacama desert. Int. J. Climatol. 23, 1453–1464 (2003).
Rech, J. A. et al. Massive middle Miocene gypsic paleosols in the Atacama Desert and the formation of the Central Andean rain-shadow. Earth Planet. Sci. Lett. 506, 184–194 (2019).
Wang, F. et al. Beryllium-10 concentrations in the hyper-arid soils in the Atacama Desert, Chile: Implications for arid soil formation rates and El Niño driven changes in Pliocene precipitation. Geochim. Cosmochim. Acta 160, 227–242. https://doi.org/10.1016/j.gca.2015.03.008 (2015).
Ericksen, G. E. Geology and Origin of the Chilean Nitrate Deposits. Report No. 1188, 37 (USGS, Washington, 1981).
Binnie, S. et al. The origins and implications of paleochannels in hyperarid, tectonically active regions: The northern Atacama Desert, Chile. Glob. Planet. Change 185, 103083. https://doi.org/10.1016/j.gloplacha.2019.103083 (2020).
Tooth, S. Process, form and change in dryland rivers: A review of recent research. Earth Sci. Rev. 51, 67–107 (2000).
Tooth, S. & Nanson, G. C. Equilibrium and nonequilibrium conditions in dryland rivers. Phys. Geogr. 21, 183–211 (2000).
Griffiths, J., Fookes, P., Goudie, A. & Stokes, M. Processes and landforms in deserts. Geol. Soc. Lond. Eng. Geol. Spec. Publ. 25, 33–95 (2012).
Billi, P., Demissie, B., Nyssen, J., Moges, G. & Fazzini, M. Meander hydromorphology of ephemeral streams: Similarities and differences with perennial rivers. Geomorphology 319, 35–46 (2018).
Santos, M. G. et al. Meandering rivers in modern desert basins: Implications for channel planform controls and prevegetation rivers. Sediment. Geol. 385, 1–14 (2019).
Ielpi, A. Morphodynamics of meandering streams devoid of plant life: Amargosa River, Death Valley, California. GSA Bull. 131, 782–802 (2019).
Vásquez, P. & Sepúlveda, F. Cartas Iquique y Pozo Almonte - Región de Tarapacá No. 161–163 Escala 1:100.000. Carta Geológica de Chile Serie Geología Básica (2013).
Marquardt, R., Marinovic, S. & Muñoz, T. Geología de las ciudades de Iquique y Alto Hospicio, región de Tarapacá, Escala 1: 25.000. (2008).
Morgan, A. et al. Sedimentology and climatic environment of alluvial fans in the martian Saheki crater and a comparison with terrestrial fans in the Atacama Desert. Icarus 229, 131–156 (2014).
Kiefer, E., Dorr, M. J., Ibbeken, H. & Gotze, H. J. Gravity-based mass balance of an alluvial fan giant: The Arcas Fan, Pampa del Tamarugal, Northern Chile. Rev. Geol. Chile 24, 165–185 (1997).
Rech, J. A., Quade, J. & Hart, W. S. Isotopic evidence for the source of Ca and S in soil gypsum, anhydrite and calcite in the Atacama Desert, Chile. Geochim. Cosmochim. Acta 67, 575–586 (2003).
Diederich, J. L. et al. A 68 ka precipitation record from the hyperarid core of the Atacama Desert in northern Chile. Glob. Planet. Change 184, 103054 (2020).
Carizzo, D., González, G. & Dunai, T. J. Constricción neógena en la Cordillera de la Costa, norte de Chile: Neotectónica y datación de superficies con 21Ne cosmogénico. Rev. Geol. Chile 35, 1–38 (2008).
Wells, S. G., McFadden, L. D., Poths, J. & Olinger, C. T. Cosmogenic 3He surface exposure dating of stone pavements. Geology 23, 613–616 (1995).
Pfeiffer, M. et al. Century scale rainfall in the absolute Atacama Desert: Landscape response and implications for past and future rainfall. Quatern. Sci. Rev. 254, 106797 (2021).
Dente, E., Lensky, N. G., Morin, E. & Enzel, Y. From straight to deeply incised meandering channels: Slope impact on sinuosity of confined streams. Earth Surf. Process. Landf. 46, 1041–1054 (2021).
Hooke, J. M. River Meandering. (2020).
Schumm, S. A. River Variability and Complexity (Cambridge University Press, 2007).
Ielpi, A., Lapôtre, M. G., Gibling, M. R. & Boyce, C. K. The impact of vegetation on meandering rivers. Nat. Rev. Earth Environ. 3, 165–178 (2022).
Ewing, S. A. et al. A threshold in soil formation at Earth’s arid-hyperarid transition. Geochim. Cosmochim. Acta 70, 5293–5322 (2006).
Rech, J. A., Currie, B. S., Michalski, G. & Cowan, A. M. Neogene climate change and uplift in the Atacama Desert, Chile. Geology 34, 761–764. https://doi.org/10.1130/g22444.1 (2006).
Hartley, A. J. & May, G. Miocene gypcretes from the Calama Basin Northern Chile. Sedimentology 45, 351–364 (1998).
Watson, A. Desert gypsum crusts as plaeoenvironmental indicators: A micropetrographic study of crusts from southern Tunisia and the central Namib Desert. J. Arid Environ. 15, 19–42 (1988).
Aref, M. A. Classification and depositional environments of Quaternary pedogenic gypsum crusts (gypcrete) from east of the Fayum Depression, Egypt. Sediment. Geol. 155, 87–108 (2003).
Voigt, C., Klipsch, S., Herwartz, D., Chong, G. & Staubwasser, M. The spatial distribution of soluble salts in the surface soil of the Atacama Desert and their relationship to hyperaridity. Glob. Planet. Change 184, 103077 (2020).
Placzek, C., Quade, J., Rech, J. A., Patchett, P. & de Arce, C. P. Geochemistry, chronology and stratigraphy of Neogene tuffs of the Central Andean region. Quat. Geochronol. 4, 22–36 (2009).
May, S. M. et al. Origin and timing of past hillslope activity in the hyper-arid core of the Atacama Desert-The formation of fine sediment lobes along the Chuculay Fault System, Northern Chile. Glob. Planet. Change 184, 103057 (2020).
Jordan, T. et al. XIV Congreso Geologico Chileno (La Serena).
Lazarus, E. D. & Constantine, J. A. Generic theory for channel sinuosity. Proc. Natl. Acad. Sci. USA 110, 8447–8452 (2013).
Tal, M. & Paola, C. Dynamic single-thread channels maintained by the interaction of flow and vegetation. Geology 35, 347–350 (2007).
Braudrick, C. A., Dietrich, W. E., Leverich, G. T. & Sklar, L. S. Experimental evidence for the conditions necessary to sustain meandering in coarse-bedded rivers. Proc. Natl. Acad. Sci. USA 106, 16936–16941 (2009).
Howard, A. D. How to make a meandering river. Proc. Natl. Acad. Sci. USA 106, 17245–17246 (2009).
Fairén, A., Davies, N. S. & Squyres, S. 44th Lunar and Planetary Science Conference, Abstract.
Lapôtre, M. G., Ielpi, A., Lamb, M. P., Williams, R. M. & Knoll, A. H. Model for the formation of single-thread rivers in barren landscapes and implications for pre-Silurian and martian fluvial deposits. J. Geophys. Res. 124, 2757–2777 (2019).
Matsubara, Y. et al. River meandering on Earth and Mars: A comparative study of Aeolis Dorsa meanders, Mars and possible terrestrial analogs of the Usuktuk River, AK, and the Quinn River, NV. Geomorphology 240, 102–120 (2015).
McMahon, W. J. & Davies, N. S. The shortage of geological evidence for pre-vegetation meandering rivers. Fluvial Meanders Sediment. Prod. Rock Rec. 48, 119–148 (2018).
Gibling, M. R. & Rust, B. R. Ribbon sandstones in the Pennsylvanian Waddens Cove Formation, Sydney Basin, Atlantic Canada: The influence of siliceous duricrusts on channel-body geometry. Sedimentology 37, 45–66 (1990).
Kereszturi, Á. Fluvial Geomorphology of Mars: Background to Separate Biogenic and Abiogenic Effects and to Identify Climate Change Related Features. (2015).
Lapotre, M. G. A. & Ielpi, A. AGU Fall Meeting Abstracts.
Lapôtre, M. G. & Ielpi, A. The pace of fluvial meanders on Mars and implications for the western delta deposits of Jezero crater. AGU Adv. 1, e2019AV000141 (2020).
Ielpi, A. & Lapôtre, M. G. Barren meandering streams in the modern Toiyabe Basin of Nevada, USA, and their relevance to the study of the pre-vegetation rock record. J. Sediment. Res. 89, 399–415 (2019).
Allen, J. Free meandering channels and lateral deposits. Sediment. Struct. 2, 53–100 (1982).
Zinelabedin, A., Riedesel, S., Reimann, T., Ritter, B. & Dunai, T. J. Testing the potential of using coarse-grain feldspars for post-IR IRSL dating of calcium sulphate-wedge growth in the Atacama Desert. Quat. Geochronol. 71, 101341 (2022).
Sager, C., Airo, A., Arens, F. L. & Schulze-Makuch, D. New type of sand wedge polygons in the salt cemented soils of the hyper-arid Atacama Desert. Geomorphology 373, 107481 (2021).
Williams, R. M. et al. Inverted channel variations identified on a distal portion of a bajada in the central Atacama Desert, Chile. Geomorphology 393, 107925 (2021).
Merritt, D. M. Reciprocal Relations between Riparian Vegetation Fluvial Landforms and Channel Processes (Academic Press, 2020).
Azua-Bustos, A., González-Silva, C. & Fairén, A. G. The Atacama Desert in Northern Chile as an analog model of Mars. Front. Astron. Space Sci. 8, 810426. https://doi.org/10.3389/fspas (2022).
Ehlmann, B. L. & Edwards, C. S. Mineralogy of the Martian surface. Annu. Rev. Earth Planet. Sci. 42, 291–315 (2014).
Bibring, J.-P. et al. Global mineralogical and aqueous Mars history derived from OMEGA/Mars Express data. Science 312, 400–404 (2006).
Christensen, M. O., Hamilton, V., Edwards, C., Wray, J. & Anderson, F. S. Aqueous Mineral Deposits in an Ancient, Channeled, Equatorial Terrain. PR (2008).
Osterloo, M. et al. Chloride-bearing materials in the southern highlands of Mars. Science 319, 1651–1654 (2008).
Mustard, J. F. et al. Hydrated silicate minerals on Mars observed by the Mars Reconnaissance Orbiter CRISM instrument. Nature 454, 305–309 (2008).
Poulet, F. et al. Phyllosilicates on Mars and implications for early Martian climate. Nature 438, 623–627 (2005).
Sharp, R. P. & Malin, M. C. Surface geology from Viking landers on Mars: A second look. Geol. Soc. Am. Bull. 95, 1398–1412 (1984).
Settle, M. Formation and deposition of volcanic sulfate aerosols on Mars. J. Geophys. Res. 84, 8343–8354 (1979).
Franz, H. B., King, P. L. & Gaillard, F. Volatiles in the Martian Crust 119–183 (Elsevier, 2019).
Mangold, N. et al. Spectral and geological study of the sulfate-rich region of West Candor Chasma, Mars. Icarus 194, 519–543 (2008).
Robertson, K. & Bish, D. Constraints on the distribution of CaSO4· nH2O phases on Mars and implications for their contribution to the hydrological cycle. Icarus 223, 407–417 (2013).
Binnie, A. et al. Accelerated late quaternary uplift revealed by 10 Be exposure dating of marine terraces, Mejillones Peninsula, northern Chile. Quat. Geochronol. 36, 12–27 (2016).
Farbod, Y. et al. Spatial variations in late Quaternary slip rates along the Doruneh Fault System (Central Iran). Tectonics 35, 386–406 (2016).
Kohl, C. & Nishiizumi, K. Chemical isolation of quartz for measurement of in-situ-produced cosmogenic nuclides. Geochim. Cosmochim. Acta 56, 3583–3587. https://doi.org/10.1016/0016-7037(92)90401-4 (1992).
Binnie, S. A. et al. Separation of Be and Al for AMS using single-step column chromatography. Nucl. Instrum. Methods Phys. Res. Sect. B 361, 397–415 (2015).
Dewald, A. et al. CologneAMS, a dedicated center for accelerator mass spectrometry in Germany. Nucle. Instrum. Methods Phys. Res. Sect. B 294, 18–23. https://doi.org/10.1016/j.nimb.2012.04.030 (2013).
Codilean, A. T. et al. Single-grain cosmogenic 21Ne concentrations in fluvial sediments reveal spatially variable erosion rates. Geology 36, 159–162 (2008).
Ritter, B., Vogt, A. & Dunai, T. J. Technical Note: Noble gas extraction procedure and performance of the Cologne Helix MC Plus multi-collector noble gas mass spectrometer for cosmogenic neon isotope analysis. Geochronology 3(2), 421–431 (2021).
Lifton, N., Sato, T. & Dunai, T. J. Scaling in situ cosmogenic nuclide production rates using analytical approximations to atmospheric cosmic-ray fluxes. Earth Planet. Sci. Lett. 386, 149–160. https://doi.org/10.1016/j.epsl.2013.10.052 (2014).
Balco, G., Stone, J. O., Lifton, N. A. & Dunai, T. J. A complete and easily accessible means of calculating surface exposure ages or erosion rates from (10)Be and (26)Al measurements. Quat. Geochronol. 3, 174–195. https://doi.org/10.1016/j.quageo.2007.12.001 (2008).
Acknowledgements
Thanks go to Elena Voronina, Damián López, and Tomasz Góral for help processing samples for 10Be and 26Al analysis at the University of Cologne, and Luigia Di Nicola for the Ne isotope analysis at SUERC. For supporting DEM analysis, we want to thank Joel Mohren. Finally, we would like to thank Eduardo Campos and colleagues at the Universidad Catolica del Norte at Antofagasta for their patient and essential support. Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Projektnummer 268236062 – SFB 1211.
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Contributions
B.R., field-work, sample preparation, TCN analysis, data evaluation, manuscript writing. J.L.D. field-work, sample preparation, TCN analysis, S.A.B. data evaluation, manuscript writing. F.M.S. TCN 21Ne analysis, data evaluation. V.W. tephra analysis. A.B. DEM generation. T.J.D. fieldwork, data evaluation, manuscript writing. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ritter, B., Diederich-Leicher, J.L., Binnie, S.A. et al. Impact of CaSO4-rich soil on Miocene surface preservation and Quaternary sinuous to meandering channel forms in the hyperarid Atacama Desert. Sci Rep 12, 17951 (2022). https://doi.org/10.1038/s41598-022-22787-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-22787-9
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.