Introduction
Lacustrine sediments are formed by accumulated organic and inorganic particles derived from different bedrock sources usually transported by rivers under continuous and episodic sediment loads [1]-[2]. Each lacustrine setting provides a unique record of the geographic, climatic, and geological variations occurring over time in a basin [3]-[7]. Proxies such as grain size, mineralogy, isotopes, organic matter, dating, and geochemistry provide valuable information to reconstruct the environmental record of such lacustrine environments [8]-[9]. This environmental record includes the analysis of sediment bedrock sources, mineral phases, sediment accumulation rates, weathering intensity, heavy metal thresholds, and chemical processes to identify past environmental disturbances [10]-[12].
The climate of the last millennium has been widely investigated worldwide. However, some tropical lacustrine settings remain unexplored. That is the case of swamps in the Atrato River delta in Colombia, in which paleoenvironmental reconstructions are limited to characterizations based on water, surface sediments, or shallow cores (< 80 cm) [13]-[15]. This study focuses on the flood plain area of the Atrato River, considered one of the most important rivers of western Colombia, with more than 20 lacustrine systems spread along its vast flood plain [13]. The Marriaga Swamp is the northernmost lacustrine landform of the Atrato River delta (8° 06' 70" N; 76° 58' 28" W; 15 m.a.s.l.) with a surface area of 12.47 km2 and a maximum water depth (measured during this research) of 4 m (Fig. 1). The swamp is connected to two main streamflows: the Atrato River to the west and the Caño Tarena to the north. The area has a tropical climate corresponding to the Intertropical Convergence Zone (iTCz) with average annual rainfall from 2,500 to 3,000 mm [16]-[17]. The average temperature is estimated at 33.8 °C. Humidity ranges from 80 to 86 % throughout the year [18].
The principal human activity in the area is fishing by the citizens of a small village named Marriaga. This village is located on the east side of the swamp, along the inlet feeds from Atrato River, just in the main entrance towards the swamp.
Bedrock in the area consists of Cretaceous rocks of igneous origin (basalts, andesites, and granitoids) cropping to the west and south of the Marriaga Swamp in the Baudó-Darien mountain ridge [19]-[20]. A broad group of Eocene sedimentary rocks outcrop towards the west of the Marriaga Swamp and comprise cherts, detrital limestones, mudstones, carbonaceous sandstones, and marls [21]. The Quaternary deposits of the vast alluvial plain are composed of clay and silt sediments [22] surrounding the Marriaga Swamp.
This research aims to determine the depositional processes and climatic impacts of the Marriaga Swamp, northern Colombia, in the last millennium, using quantitative data obtained from XRF geochemical analyses. It is the first approach to understanding the Marriaga Swamp's paleoenvironmental features with major and trace elements interpretations within the framework of the Atrato delta.
Materials and Methods
We drilled a sediment core at the center of the Marriaga Swamp with a Livingstone Bolivia corer. The core consists of two overlapping sections with a total length of 220 cm, retrieved at the central area of the swamp where the depth of water to surface sediments was 4 m. We sliced the core into 16 sections, between 5 and 10 cm long each, according to visually identified stratigraphic changes. We chose eight samples distributed on top, middle, and bottom sections of the sequence for detailed characterizations.
The physical properties of the sediment layers were initially determined with visual and tactile inspections of color, texture, and organic matter content. The geochemical analyses followed the patterns related to the three segments previously defined by their sediment characteristics, age, and organic carbon (OC) content. We measured variations of major and trace elements and determined correlations from which we analyze dry and humid periods, biological productivity, carbonate precipitation, fluvial-marine sediment input fluctuations, and human impacts.
We performed the geochemical analyses at the Universidad de Antioquia Spectrometry Laboratory, using an ARL OPTIM'X WDRRF consisted of grounding samples in an S-1000 RETCH Centrifugal Ball Mill with a zirconium oxide vessel at 3,000 rpm for 5 minutes. The grounded powder was calcinated at 550 °C for three hours to determine Loss of Ignition (LOI) and calculate OC. The geochemical data included concentrations of ten oxides (SiO2, Al2O3, Fe2O3, MgO, Na2O, CaO, K2O, TiO2, P2O5, and MnO) and 22 trace elements (Cl, S, Ba, V, Cu, Zr, Sm, Cr, Zn, Sr, W, Tb, Ni, Co, La, Rb, Br, Y, Sc, Ga, and Cd). Measurements were made with the standard and duplicate samples (Table 1).
Samples | S1 | S2 | S3 | S4 | S5 | S6 | S7 | S8 |
---|---|---|---|---|---|---|---|---|
LOI | 17,83 | 16,33 | 15,20 | 14,38 | 13,49 | 21,22 | 21,61 | 19,35 |
OC | 10,34 | 9,47 | 8,82 | 8,34 | 7,82 | 12,31 | 12,54 | 11,22 |
Depth | 5,00 | 33,00 | 61,00 | 82,00 | 113,00 | 162,00 | 201,00 | 219,00 |
SiO2 | 41,62 | 41,16 | 45,68 | 45,25 | 46,41 | 44,27 | 46,79 | 52,39 |
20,37 | 17,89 | 17,83 | 18,13 | 18,64 | 19,06 | 19,16 | 19,41 | |
Fe2O3 | 11,52 | 11,27 | 10,56 | 10,53 | 10,50 | 9,14 | 8,88 | 8,38 |
MgO | 3,38 | 2,86 | 3,31 | 3,44 | 3,37 | 3,58 | 3,94 | 4,04 |
Na2O | 0,94 | 1,16 | 1,39 | 1,50 | 1,31 | 2,08 | 2,43 | 2,50 |
CaO | 1,43 | 1,48 | 1,75 | 1,77 | 1,71 | 1,71 | 2,57 | 2,76 |
K2O | 0,90 | 0,98 | 1,25 | 1,33 | 1,32 | 1,40 | 1,37 | 1,51 |
TiO2 | 0,65 | 0,67 | 0,76 | 0,78 | 0,76 | 0,78 | 0,79 | 0,74 |
P2O5 | 0,32 | 0,19 | 0,09 | 0,17 | 0,10 | 0,27 | 0,26 | 0,15 |
MnO | 0,40 | 0,19 | 0,11 | 0,10 | 0,11 | 0,09 | 0,12 | 0,11 |
Cl | 606 | 3320 | 3950 | 4140 | 2490 | 6880 | 7260 | 4550 |
S | 2630 | 3840 | 13100 | 18900 | 17000 | 6370 | 2610 | 976 |
Ba | 589 | 429 | 497 | 514 | 502 | 810 | 937 | 977 |
V | 270 | 214 | 206 | 210 | 205 | 253 | 215 | 217 |
Cu | 273 | 261 | 269 | 246 | 254 | 253 | 220 | 207 |
Zr | 123 | 165 | 233 | 197 | 177 | 304 | 380 | 456 |
Sm | 146 | 175 | 150 | 152 | 199 | 170 | 121 | 0 |
Cr | 104 | 138 | 180 | 179 | 194 | 174 | 187 | 166 |
Zn | 148 | 134 | 144 | 138 | 142 | 157 | 126 | 122 |
Sr | 111 | 159 | 202 | 200 | 192 | 201 | 260 | 274 |
W | 108 | 103 | 112 | 105 | 90 | 61 | 0 | 74 |
Tb | 129 | 121 | 123 | 107 | 153 | 122 | 111 | 94 |
Ni | 111 | 89 | 100 | 111 | 115 | 120 | 113 | 112 |
Ar | 178 | 186 | 146 | 0 | 111 | 149 | 0 | 0 |
Co | 65 | 71 | 81 | 80 | 72 | 78 | 61 | 38 |
La | ND | ND | ND | ND | ND | ND | ND | 65 |
Rb | 27 | 36 | 38 | 44 | 35 | 51 | 37 | 40 |
Br | 22 | 48 | 42 | 29 | 26 | 39 | 23 | 14 |
Y | 22 | 24 | 27 | 25 | 25 | 26 | 28 | 16 |
Sc | 22 | 20 | 26 | 0 | 21 | 0 | 0 | 29 |
Ga | 11 | 20 | 13 | 16 | 17 | 26 | 20 | 18 |
Cd | ND | ND | ND | ND | ND | ND | ND | 65 |
Note. ND: Not detected.
Source: Own elaboration.
The ages at the top, middle, and bottom sections were determined with the Infrared Stimulated Luminescence (IRSL) technique in the Institute of Earth and Environmental Sciences laboratory at the Albert-Ludwigs-University of Freiburg, Germany [23]. The basic sample was prepared in the Universidad EAFIT photoluminescence laboratory in Medellín, Colombia. The core was split into two parts, selecting three samples on the top, middle, and bottom of the sedimentary sequence. A first dry process was carried out at 40 °C for 24 hours to obtain the water content. The robust material was prepared at the University of Freiburg red-light adapted laboratory. Before an HCl (30 %) and H2O2 (30 %) treatment, several sieve stages were completed to remove carbonates and organic matter.
Both coarse and fine fractions were used for age determination procedures.
Results and Discussion
Physical Features of Sediments
The 220 cm sedimentary sequence consists of laminated grayish to brownish fine sand, silt, and clay layers, intercalated with peat layers (Fig. 2). The ages estimated with the IRSL technique indicate 1,130 ± 90 a at the bottom of the core, around 870 ± 70 a in the middle section, and around 530 ± 40 a on the top. The core was divided into three segments, according to sediment features and estimated ages: bottom segment (Seg. i: 220-135 cm; 1,130 ± 90 to 870 ± 70 a), middle segment (Seg. ii: 135 cm-70 cm; 870 ± 70 to 530 ± 40 a), and upper segment (Seg. iii: 70-0 cm; 530 ± recent a).
Seg. i involves the coarser sediments of brownish fine sand to silty sediments, interbedded with black laminations enriched with organic matter and mica. Seg. ii comprises grey-brown alternations of silt and clay sediments with thin sandy laminations, scarce organic matter, and mica content. Seg. iii comprises greyish silt with intercalations of clay with poor content of organic matter.
Organic matter content is present in the core, with values ranging between 7.8 and 12.5 %. In the oldest segment (Seg. i), the organic matter content is considerably high, as evinced by four peaty layers. These layers indicate exposure conditions that allowed vegetation cropping between approx. 1,130 ± 90 and 870 ± 70 a. In Seg. ii (between 870 ± 70 and 530 ± 40 a), the organic matter is low and relates to a sedimentary sequence of clay that indicates subaerial and redox conditions. The coarser sediments in the most recent segment (Seg. iii) indicate the influence of high-energy flows feeding the swamp. Seg. iii is characterized by an increase in organic matter content towards the surface, approx. 530 ± 40 a to present. However, Seg. iii resembles Seg. ii in terms of grain size and the absence of clear peaty layers, except for milli-metric laminations.
Geochemistry
Major elements
Si, Al, Fe, Mg, Na, Ca, and K are associated with a siliciclastic provenance (Fig. 3). Seg. i (lower part of the sequence) has high Si, Mg, Na, Ca, and K concentrations, whereas Seg. ii (middle of the sequence) shows an increment of Fe and a reduction of Mg and Na. Seg. iii (upper sequence) exhibits high Al and Mg concentrations and reductions in Si and Ca. The patterns of Ca and Mg are more complex, probably due to the variety of bedrock sources with detrital and endogenic inputs.
High concentrations of major elements (Si, Al, K, Ti, and Fe) represent detrital contributions related to flooding events; high values of Ca demonstrate sediment reworking during dry periods or swamp level fluctuations [24], whereas Mg, Cl, and Na indicate salinity intrusions [25]-[26].
Trace elements
The extended analysis of trace elements exposes several paleoenvironmental features. Dry and humid periods, for instance, are related with (Zr+Rb)/ Sr and Ca/Ti ratios; biological productivity is represented by Mn/Fe, OC/Ti, and Mg/Ca ratios (Fig. 4). Carbonate precipitation employs Ba/Al, Sr/Al, and Ca/Al ratios (Fig. 5), and the possible anthropogenic impact is associated with Cu, Ni, Zn, Co, and Cr (Fig. 6).
The (Zr+Rb)/Sr ratio reveals terrigenous inputs, whereas the Ca/Ti ratio indicates carbonate production [26] and the relation with drought in humid conditions. High (Zr+Rb)/ Sr and Ca/Ti imply possible drought conditions observed in Seg. i (220 to 160 cm) from 1,130 ± 90 and 870 ± 70 a, with a significant fall towards more recent sediments. A relatively constant humid environment is consistent with these low ratios observed in the upper 120 cm of the core (Seg. ii and Seg. iii).
Biological productivity was important from 1,130 ± 90 to 870 ± 70 a. The middle to low Mn/Fe ratio in Seg. i indicates some intermittent, short periods of paleo-redox conditions (220-160 cm). Low values of Mn/Fe at 120 to 60 cm in Seg. ii indicate more oxygenated conditions from 870 ± 70 to 530 ± 40 a. The rise of biological productivity remains from 60 cm to the surface. Important OC/Ti ratios show high concentrations in Seg. i and Seg. iii with a significant decrease at Seg. ii. The Mg/Ca is relatively high at Seg. iii, close to the surface. The fine lamination of sediments preserved in upper and lower sediments indicates anoxic conditions in the Marriaga Swamp in Seg. i and Seg. iii
Carbonate precipitation (Fig. 5.), interpreted from Ba/Al, Sr/Al, and Ca/Al ratios, shows high Ba/Al profile values in Seg. i (220-120 cm), with a notable increase in the sequence's upper sediments, suggest authigenic processes. The presence of ion-rich waters and chemical precipitation is interpreted from Sr/Al and may indicate a saline, brackish stage in the Marriaga Swamp from 1,130 ± 90 to 530 ± 40 a (220-60 cm), with clear decreasing values from 60 cm to surface.
The brackish state is consistent with the fact that the Marriaga Swamp could be affected by sea-water inputs during high tidal conditions. High values of Ca/Al in the lower sequence (220-60 cm) correspond to half of Seg. i. Moderate and constant values from 160 to 60 cm are compatible with half of Seg. i and whole Seg. ii. A sharp decrease of Ca/ Al occurs in the upper sequence (Seg. iii), consistent with more fluvial inputs.
Heavy metal enrichments have been typically associated with a human impact [25], but they could also be associated with natural thresholds due to bedrock origins [28]-[29]. We considered that the high values found in the Marriaga Swamp profile (Fig. 6) are a combination between natural thresholds from basic bedrock weathering and specific anthropogenic effects associated with mining activities upstream the Atrato River that can incorporate mined sediments during floods into the Marriaga Swamp [30].
Cu, Ni, Zn, Co, and Cr concentrations were normalized with average Upper Continental Crust (UCC) concentrations [31]-[32] and compared to regional background data (Fig. 7) from stream sediments sampled by the Geological Survey of Colombia upstream the Marriaga Swamp in the Atrato River [33]. The purpose of this normalization was to determine if the heavy metal concentrations in the Marriaga Swamp showed extreme thresholds.
The normalization of heavy metals against the UCC shows that all concentrations are several times higher than the UCC values, which points to the possible anomalous enrichment of Cu, Ni, Zn, Co, and Cr in the swamp. The values for the regional background are also higher compared to the UCC and conservatively close to those found in the Marriaga Swamp sediments, except for Cr, which is higher for the regional records. Natural leaching could be one of the causes of the high values in the Marriaga Swamp sediments, which is also related to the high Chemical Index of Alteration (CIA) and the dilution effect of silica that promotes heavy metal dilutions [34].
If human activities were the most substantial reason for heavy metal enrichment in the Marriaga Swamp, the concentrations would be higher only in the upper sediments (Seg. ii) after 500 a, resulting from colonial mining [35]. However, the values in the profile (Fig. 7) demonstrate a constant presence along with the sequences. Since intense human presence or activities are not reported in the area before 500 a in the late Holocene (corresponding to Seg. i and part of Seg. ii: ~1130 ± 70 to 500 a), it is unlikely heavy metals are only a contribution of human sources.
Chemical Index of Alteration (CIA)
The CIA is used as a proxy of weathering grade and paleoclimatic signal to distinguish humid and arid conditions [36]-[38]. Te CIA has been employed to examine the relationship between some major elements and assume the disintegration of feldspars and further formation of clay minerals during the weathering processes [39]. We used the CIA to determine the production of detrital sediments and their contributions to sediment accumulation in the lacustrine setting of the Marriaga Swamp. Te CIA is the molecular percentage of principal (oxides) elements (Equation 1):
Where Al2O3, Na2O, and K2O are the mole percentage and CaO* is the calcium content of silicates.
An increase in CIA values is generally attributed to climatic fluctuations from cold to warm and humid [40]-[41]. For the Marriaga Swamp, the CIA is in a medium range for the lower sedimentary sequence (220-162 cm) and high over 162 cm to near-surface (Fig. 8). The CIA increases within the medium range from bottom to upper sequences, indicating climatic changes from 1,130 ± 90 a to present, from more arid and maybe colder conditions in late Holocene to humid and warmer conditions nowadays.
The CIA ranges hint at low weathering from 0 to 40, moderate from 40 to 70, and high from 70 to 100 [40]. The Marriaga Swamp CIA values range from 65 to 82, with an average of 72, which means the Marriaga Swamp sediments are highly weathered.
Statistical analysis
We applied a Pearson correlation to identify the statistical significance association (r > 0.70; p-value > 0.05) between OC, depth, major elements, and trace elements (Table 2). The results can be clustered in three main groups: 1) OC has a positive correlation with Na and Ba, concerning organic features and biological productivity, 2) Si, Fe, Mg, Na, Ca, K, Ba, Cl, and Cu are positively correlated with depth, concerning the interpretation of bedrock components, and 3) Cr, Zn, and Co, related with heavy metal inputs.
OC | Depth | Si | Al | Fe | Mg | Na | Ca | K | Ti | P | Mn | Cl | S | Ba | Cu | Cr | Zn | Ni | Co | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
OC | 1,00 | |||||||||||||||||||
Depth | 0,60 | 1,00 | ||||||||||||||||||
Si | 0,18 | 0,81 | 1,00 | |||||||||||||||||
Al | 0,54 | 0,18 | 0,10 | 1,00 | ||||||||||||||||
Fe | -0,67 | -0,98 | -0,80 | -0,20 | 1,00 | |||||||||||||||
Mg | 0,62 | 0,86 | 0,80 | 0,51 | -0,87 | 1,00 | ||||||||||||||
Na | 0,71 | 0,97 | 0,75 | 0,18 | -0,99 | 0,87 | 1,00 | |||||||||||||
Ca | 0,52 | 0,88 | 0,87 | 0,21 | -0,87 | 0,88 | 0,90 | 1,00 | ||||||||||||
K | 0,27 | 0,89 | 0,84 | -0,10 | -0,87 | 0,74 | 0,83 | 0,73 | 1,00 | |||||||||||
Ti | 0,14 | 0,67 | 0,53 | -0,31 | -0,64 | 0,52 | 0,62 | 0,47 | 0,86 | 1,00 | ||||||||||
P | 0,67 | -0,05 | -0,43 | 0,69 | 0,00 | 0,15 | 0,06 | -0,10 | -0,36 | -0,31 | 1,00 | |||||||||
Mn | 0,01 | -0,64 | -0,56 | 0,61 | 0,61 | -0,29 | -0,58 | -0,44 | -0,84 | -0,85 | 0,59 | 1,00 | ||||||||
Cl | 0,65 | 0,76 | 0,36 | -0,17 | -0,79 | 0,53 | 0,83 | 0,57 | 0,69 | 0,75 | 0,10 | -0,70 | 1,00 | |||||||
S | -0,78 | -0,25 | -0,05 | -0,58 | 0,32 | -0,31 | -0,38 | -0,37 | 0,17 | 0,44 | -0,60 | -0,40 | -0,16 | 1,00 | ||||||
Ba | 0,85 | 0,88 | 0,66 | 0,54 | -0,91 | 0,92 | 0,92 | 0,85 | 0,63 | 0,38 | 0,34 | -0,23 | 0,64 | -0,60 | 1,00 | |||||
Cu | -0,49 | -0,89 | -0,83 | -0,19 | 0,86 | -0,83 | -0,89 | -0,95 | -0,74 | -0,46 | 0,05 | 0,47 | -0,57 | 0,33 | -0,82 | 1,00 | ||||
Cr | -0,09 | 0,61 | 0,55 | -0,49 | -0,54 | 0,36 | 0,50 | 0,43 | 0,82 | 0,93 | -0,58 | -0,92 | 0,61 | 0,53 | 0,20 | -0,42 | 1,00 | |||
Zn | -0,12 | -0,46 | -0,59 | 0,05 | 0,41 | -0,42 | -0,48 | -0,76 | -0,29 | -0,06 | 0,25 | 0,21 | -0,18 | 0,31 | -0,40 | 0,77 | -0,17 | 1,00 | ||
Ni | 0,37 | 0,56 | 0,40 | 0,57 | -0,53 | 0,69 | 0,46 | 0,31 | 0,57 | 0,52 | 0,29 | -0,17 | 0,30 | 0,09 | 0,56 | -0,34 | 0,33 | 0,23 | 1,00 | |
Co | -0,44 | -0,55 | -0,66 | -0,55 | 0,55 | -0,64 | -0,56 | -0,75 | -0,27 | 0,18 | -0,10 | -0,11 | -0,04 | 0,67 | -0,69 | 0,75 | 0,16 | 0,71 | -0,15 | 1,00 |
Source: Own elaboration
A Principal Component Analysis (PCA) of the variables OC, depth, major elements, and trace elements (Fig. 9) reveal that PC-1 explains 54 % of the total variance and is positively linked to Al, OC, Ba, Na, Mg, Ca, Si, Cl, K, Ti, Cl, and depth and negatively linked to Fe, Cu, Zn, Co, and S. PC-2 accounts for 24.1 % of the total variance and is positively linked to Mg and P and negatively linked to Cu, Zn, Co, S, Cr, Ti, K, Cl, and Si.
The positive values of PC-1 indicate more detrital and allochthonous components, whereas the negative values denote more anoxic conditions and heavy metal components, which also could justify the humid and dry alternations inside the Marriaga Swamp. Positive values of PC-2 may imply an endogenic precipitation process, supported by the Mg and P, whereas negative values of PC-2 suggest heavy metal precipitation, combined with saline-brackish swamp regimes.
Geochemical signatures in northwestern Colombia during the Late Holocene
The Marriaga Swamp sedimentary sequence analyses contribute and complement the paleoenvironmental reconstructions for northwestern Colombia [42]-[43]. Tis research specifies the alternation of humid and dry conditions in the Atrato River delta setting in the last millennium. Tese climatic alternations are consistent with the migration of the ITCZ during the Holocene in the north of South America [44]-[46].
Based on geochemical proxies such as the oscillations of Si, Al, K, and Ti (representing detrital inputs) and Ba/Al and Sr/Al ratios, we identify different authigenic lacustrine and allochthonous processes in the Marriaga Swamp. The lower sequence (Seg. iii), for instance, is characterized by the predominance of a dry (arid) phase, where the high OC represents the potential for lamination of fine sand (coarsest sediments of the sequence) [47]. Sediment reworking and carbonate precipitations also occur at the lower end of the sequence, as illustrated by the Ca/Al ratio and the high concentrations of heavy metals (Cu, Ni, Co, Zn, and Cr).
On the other hand, Seg. ii represents a transition period between dry to moderately humid conditions. Tis part of the sedimentary sequence exhibits a light decrease in detrital inputs (Si, Al, K, Ti, and Fe) and a reduction of OC and carbonate production, represented by a low, Ca/Ti ratio and the deposition of finer sediments (wide clay layer interbedded with some coarser lamination towards 90 cm deep).
The upper sequence (Seg. iii) is the clearest sedimentary section representing more humid conditions and fluvial influence, with a significant drop of detrital Si, rise of Al, and moderate presence of Na and Ca. A rise in biological productivity is inferred by high Fe concentrations and Mg/Ca ratios. However, the tidal influence is represented by high Ba/Al ratios towards the lowest part of the sequence (< 25 cm). Concentrations of heavy metals also prove an increase of Cu, Ni, and Zn for Seg. iii, possibly related to human activities such as mining, whereas the reduction of Cr and Co could be interpreted as the influence of intense weathering and dilution of natural sources for these two metals.
Conclusions
The sedimentary sequence and the geochemistry of the Marriaga Swamp expose significant environmental changes that occurred during the last 1130 a in the Atrato River delta. Detrital contributions portray flooding events, generally associated with high concentrations of Si, Al, K, Ti, and Fe in sediments. In contrast, high concentrations of Ca reveal the possibility of sediments reworking during dry periods and level fluctuations in the Marriaga Swamp. Biological productivity manifests in the high contents of organic carbon and high carbonate precipitation (Mg/Ca). Higher Ca/Al, Sr/Al, and Ba/Al ratios are also indicative of endogenic calcite precipitation. The presence of high concentrations of heavy metals (Cu, Zn, Co, Ni, and Cr) along the entire sequence of the Marriaga Swamp sediments imply that their origin is likely a combination of natural bedrock sources and human impacts, especially towards the more recent deposits. High regional background concentrations for the same metals suggest natural enrichments of these heavy metals due to dilution processes in the sources that feed the swamp. The CIA is high along the sedimentary sequence, with a remarkable increase towards the surface, which points to a transition from arid to more humid conditions. The statistical correlations show three main clusters that distinguish among organic and biological productivity, bedrock source components, and heavy metal inputs.