Next Article in Journal
The Possibility of Managed Aquifer Recharge (MAR) for Normal Functioning of the Public Water-Supply of Zagreb, Croatia
Next Article in Special Issue
The River Influence Controls Water Quality and Spatio-Temporal Microalgal Distribution in Pacific Estuaries (Padre Ramos and Salinas Grandes) of Nicaragua
Previous Article in Journal
Climate Change Impacts on Water Resources of the Largest Hydropower Plant Reservoir in Southeast Brazil
Previous Article in Special Issue
Regional Hydrogeochemical Evolution of Groundwater in the Ring of Cenotes, Yucatán (Mexico): An Inverse Modelling Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Particulate Organic Carbon in the Tropical Usumacinta River, Southeast Mexico: Concentration, Flux, and Sources

1
Programa de Posgrado en Ciencias del Mar y Limnología, Universidad Nacional Autónoma de México, Av. Universidad 3000, Alcaldía Coyoacán, C. P., Mexico City 04510, Mexico
2
Grupo de Investigación en Limnología Tropical, FES Iztacala, Universidad Nacional Autónoma de México, Av. de los Barrios No.1, Los Reyes Iztacala, Tlalnepantla 54090, Mexico
3
Instituto de Investigaciones en Ecosistemas y Sustentabilidad, Universidad Nacional Autónoma de México, AP 27-3, Santa María de Guido, Morelia, Michoacán 58090, Mexico
4
Departamento de Biogeoquímica y Ecología Microbiana, Museo Nacional de Ciencias Naturales (CSIC), Serrano 115 dpdo, 28006 Madrid, Spain
*
Author to whom correspondence should be addressed.
Water 2021, 13(11), 1561; https://doi.org/10.3390/w13111561
Submission received: 30 April 2021 / Revised: 25 May 2021 / Accepted: 29 May 2021 / Published: 31 May 2021

Abstract

:
Particulate organic carbon (POC) derived from inland water plays an important role in the global carbon (C) cycle; however, the POC dynamic in tropical rivers is poorly known. We assessed the POC concentration, flux, and sources in the Usumacinta, the largest tropical river in North America, to determine the controls on POC export to the Gulf of Mexico. We examined the Mexican middle and lower Usumacinta Basin during the 2017 dry (DS) and rainy (RS) seasons. The POC concentration ranged from 0.48 to 4.7 mg L−1 and was higher in the RS, though only in the middle basin, while remaining similar in both seasons in the lower basin. The POC was predominantly allochthonous (54.7 to 99.6%). However, autochthonous POC (phytoplankton) increased in the DS (from 5.1 to 17.7%) in both basins. The POC mass inflow–outflow balance suggested that floodplains supply (C source) autochthonous POC during the DS while retaining (C sink) allochthonous POC in the RS. Ranging between 109.1 (DS) and 926.1 t POC d−1 (RS), the Usumacinta River POC export to the Gulf of Mexico was similar to that of other tropical rivers with a comparable water discharge. The extensive floodplains and the “Pantanos de Centla” wetlands in the lowlands largely influenced the POC dynamics and export to the southern Gulf of Mexico.

1. Introduction

During the last two decades, several studies mentioned that inland aquatic ecosystems are important agents in the coupling of biogeochemical cycles between continents, the atmosphere, and oceans (e.g., [1,2]). Rivers play a mediatory function in carbon (C) storage, processing, and release among these compartments. Rivers also have effects on the regional C cycle since they couple landscapes through the lateral export from terrestrial ecosystems [3]. However, integral approaches to the carbon biogeochemistry associations between habitats remain poorly explored in inland water research [2].
Rivers annually deliver approximately 0.5 Pg of organic carbon (OC) to the world’s oceans, with nearly half arriving in particulate form [4,5,6]. The difference between rock-derived and biosphere-derived particulate organic carbon (POC) loads arriving into oceans from terrestrial systems through rivers contributes to the regulation of atmospheric and terrestrial C reservoirs over a geological timescale by sequestration in marine sediments [5,7,8,9]. The activity of primary producers in rivers may regulate the C inventory on a shorter scale. For example, the fluvial C inventory may be affected by algal in situ production, increasing the POC (autochthonous) and reducing the CO2 concentrations [10]. Therefore, a considerable fraction of the C dynamics in river ecosystems depends on the variability of POC within the C pool [11], including the temporal and spatial changes occurring in their sources (allochthonous vs. autochthonous).
Seasonal hydrologic variation often guides POC supply and flux in tropical rivers. While low-water discharge conditions enhance autochthonous production and decrease POC flux [12], high-water discharge conditions increase POC flux with the allochthonous terrestrial matter via erosion [13,14]. However, the seasonal effect on the flux and source of POC is not always clear in tropical rivers, as in some mountainous watersheds [15].
Studies on POC seasonal variations usually circumscribe to limited zones of the rivers. This approach prevents a comprehensive understanding of POC dynamics along the fluvial continuum, from headwaters downstream to the river’s mouth. For example, the presence of extensive floodplains and wetlands may diminish the POC through transformation to dissolved forms (mineralized as CO2) and may increase it via primary production activity. This affects the incoming and outgoing POC flux [16,17,18].
Studies regarding riverine POC flux in the tropics—including Latin America—are still limited [19]. Although available data suggest that tropical rivers are important for the transport of C to oceans, POC flux estimates are still uncertain due to the lack of data, the poor representation of spatial and temporal C dynamics, and indirect estimates, such as for suspended matter flux [8]. Therefore, providing information on the spatial variation in the amount and sources of fluvial POC in the tropics in a seasonal framework will help to unravel the functioning of tropical rivers in the regional and global C cycle.
The Usumacinta River is a transboundary river (Mexico–Guatemala) that is among the largest (i.e., in terms of water discharge, length, drainage basin) tropical rivers in North America and the longest in Mesoamerica [20]. Through its river basin, marked biophysical and hydrological shifts occur [21,22,23]: the Mexican portion of the middle basin is dominated by steep slopes and covered by well-preserved tropical rainforests, such as the Selva Lacandona [24], while a large expanse of floodplain wetlands and shallow lakes (Pantanos de Centla Biosphere Reserve) extend across the lower basin [25]. The strong extreme climate seasonality and complex mosaic of landscapes should exert strong control over both the river hydrology and POC dynamics in the tropical Central America region, but no information on this topic is currently available.
In this study, our primary goal was to determine the effect of hydrologic tropical seasons (rainy versus dry season) on fluvial POC throughout the basin, assessing its transport, processing, and export. To understand this effect, herein we addressed two main questions: (1) Do the POC sources (autochthonous and allochthonous), concentration, and flux change seasonally along the Usumacinta River and its main tributaries? (2) What effect do the floodplain and wetlands located in the lower basin of the Usumacinta River have on the seasonal POC export? We expected (i) a higher proportion of allochthonous POC in the high-relief, forested landscapes of the middle basin compared to the flat, lower basin where better conditions for autochthonous production may exist (lower water current and low turbidity); (ii) an increase in the POC concentration and flux along the river in the rainy season this is associated with a higher erosion and dragging of terrestrial matter; (iii) that floodplains and wetlands add POC to the fluvial export toward the Gulf of Mexico due to their potential algal development in the dry seasons while retaining POC during the rainy season.
Changes in POC concentration, water quality variables, and C/N ratios of particulate matter along with hydrological variables were analyzed throughout two consecutive seasons to unravel the mechanisms that were most likely to be controlling the POC flux and sources in the Usumacinta River basin. This study contributes to improving the understanding of tropical rivers as biogeochemical conduits of continental carbon. Furthermore, this contribution provides the first estimate of POC export from a Mexican river to the southern Gulf of Mexico.

2. Materials and Methods

2.1. Study Site

The Usumacinta River is located within the Grijalva–Usumacinta Hydrological Basin to the southeast of Mexico and north of Guatemala. The basin lies between 14°55′ and 18°35′ north and 91°20′ and 94°15′ west [26]. The Usumacinta basin covers ≈77,000 km2 of the Guatemalan and Mexican territories [27], while the Grijalva basin covers ≈57,000 km2 [28]. The maximum altitude in the basin is 3800 m [20].
The Usumacinta River originates in the Cuchumatanes mountain range (upper basin), Guatemala, and runs through the Mexican States of Chiapas and Tabasco (middle and lower basins). The adjacent Grijalva River, which is artificially connected to the Usumacinta River before the river mouth at the southern Gulf of Mexico [29], increases the total water discharge through the Usumacinta River’s mouth. The maximum order (Strahler method) of both rivers at their confluence is 7 [30]. The mean annual water discharge of both rivers is 100 × 109 m3, where nearly 64% of this corresponds to the Usumacinta River discharge [31].
The climate is tropical and wet, with a mean annual temperature of between 8 and 12 °C in the highlands (>1000 m altitude) and 26 and 36 °C in the lowlands (<1000 m) [32]. The rainy season (RS) takes place between May and October and the dry season (DS) between November and April [33,34]; the average precipitation for the seasons is 4500 and 200 mm, respectively [35]. During the RS, characteristic tropical storms increase the Usumacinta River’s water discharge, causing flooding in the lowlands. The geomorphological units that characterize the Usumacinta basin are the mountainous area in the highlands and the lowland coastal plain [20]. The mountainous area is composed of the Sierra de Los Cuchumatanes (limestone and dolomites from the Cretaceous), the Sierra Madre de Chiapas (diorites and granites from the Paleozoic), and the Altos de Chiapas (marine and continental carbonates from the Mesozoic with volcanic deposits from the Cenozoic). The coastal plain is a low-relief area formed by sedimentary Cenozoic rocks of alluvial origin [36].
Evergreens, coniferous, and oak forests covered 85% of the middle basin in the 2000s, but they have seen a near 14% reduction in area in the last 15 years [37,38]. The replacement of the natural vegetal cover with agricultural land, space for raising cattle, and urban development in the middle basin has caused landslides and floods during extreme rain events [39]. Mangrove forests, along with other hydrophytes, extend across the lower portion of the basin, covering 8% of the Usumacinta–Grijalva basin [40]. Submerged and floating macrophytes are also present in the lakes and floodplains connected to the Usumacinta’s main channel [41].

2.2. Sampling

Two sampling campaigns were carried out in the Mexican portion of the Usumacinta River (middle and lower basins); the first was conducted at the end of the DS in May 2017 and the second was conducted in the RS in November 2017, just after the rainiest month (October). We selected seven sites along the main channel and major tributaries of the Usumacinta River: 4 in the middle (M1 to M4) and 3 in the lower basin (LU, LG, and LPP; Figure 1).
M1 corresponds to the Lacantún River, which drains 15,772 km2 of the highlands of Chiapas. From this site, the course stretches 656 km to the river mouth. M2, the Tzendales River, was chosen as a reference site due to its undisturbed and pristine condition, as its water emerges and crosses through the protected area of the Lacandona rainforest [37]. M2’s watershed covers 1487 km2, and from this site, the river stretches 653 km to the river mouth. Both M1 and M2 are in the Montes Azules Biosphere Reserve tropical rainforest. M3 (Frontera Corozal) and M4 (Boca del Cerro) correspond to sites along the main channel of the Usumacinta River in the middle basin (Figure 1). They have watersheds of 46,564 and 51,064 km2 and lengths of 525 and 385 km to the river mouth, respectively. LU (Rivera Alta) is the main channel of the Usumacinta River in the lower basin. It stretches 25 km to the river mouth, while its watershed spans 70,350 km2. LG refers to the Grijalva River before its confluence with the Usumacinta River. We used it as a reference to the adjacent Grijalva watershed, which drains 57,272 km2 of land. This site is also 25 km away from the river mouth. LPP corresponds to the distributary San Pedro–San Pablo River; it is 7 km away from the river mouth and drains 72,569 km2. LPP represents the second and smaller river mouth of the Usumacinta River. The three lower basin sites are in the Ramsar wetland and Biosphere Reserve “Pantanos de Centla,” which is the largest (3027 km2) wetland area of North America and the third-largest in Latin America ([42]; Figure 1).
At each sampling site, a cross-section profile was recorded with a Garmin GPSMap 526S Echo Sounder. Vertical profiles (every meter) of temperature (T), dissolved oxygen concentration (DO), pH, electrical conductivity (K25), turbidity (Turb), and oxidation/reduction potential (ORP) were determined with a Hydrolab DS5 multiparameter sonde (OTT Hydromet, Kempten, Germany). The mean current velocity was measured with a Swoffer 3000 current meter (Swoffer, Sumner, Washington, USA). Three sub-superficial (1/3 maximum depth) water samples were collected along the cross-section of the river (center and 1/3 toward both sides across the river section) in all sites but M2; based on the M2 reduced width (40 m), only one water sample was taken at the center. The water discharge (Q, m3 s−1) was calculated using cross-section profiles and the measured mean river flow velocities and then normalized by the watershed area upstream of the sampling point (L km2 s−1). Water samples for the POC, total suspended solids (TSS), and chlorophyll a (Chl-a) concentrations were taken with a horizontal Van-Dorn-type sampling bottle.

2.3. Analytical Methods

The Chl-a concentration measurement followed the EPA 445.0 method [43]. Triplicate water samples were filtered through Whatman GF/F glass microfiber filters (nominal 0.7 µm pore size). A total of 10 mL of 90% acetone was added to each filter. Chl-a was extracted at 4 °C in dark conditions for 12–24 h. After centrifugation (2700 rpm/15 min), the Chl-a extracted from the supernatant was measured with a digital fluorometer TD 10-AU (Turner Designs, San Jose, CA, USA, country). To measure the TSS concentration, duplicate water samples (300 to 1000 mL) were pre-filtered through 100 µm pore mesh to remove “large swimmers” and were furthered retained in pre-combusted (550 °C/4 h) and pre-weighed Whatman GF/F filters. The filters were then dried for at least 48 h (50 °C) and weighed again.
Duplicate water samples (20 to 80 mL) were filtered through pre-combusted (550 °C/4 h) Whatman GF/F filters. Filters were further acidified (HCl 10%) to eliminate inorganic C [44], folded in clean pre-baked aluminum foil, and frozen until analysis. The C and N content in the filters was measured via elemental analysis (Elemental Analyzer Carlo Erba NC2100, CE Elantech, Lakewood, NJ, USA) to calculate the POC concentration and the particulate carbon-to-nitrogen ratio (POC/PN). The POC concentrations and Q data were used to calculate instant POC fluxes (FPOC, t C d−1): Q × POC concentration. We also calculated the POC yields using (YPOC, kg C km2 d−1) Q × POC concentration/watershed area.

2.4. Data Analysis: POC Sources, Mass Balance, and Statistical and Geographical Analyses

Chl-a, TSS, and POC/PN were used to provide insight into the POC sources (autochthonous vs. allochthonous). The contribution of POC to the total particulate matter load (POC/TSS) was used to determine the allochthonous origin of POC (i.e., rock or soil; [9,45]); we calculated the POC contribution as a percentage. We used the Chl-a concentration to obtain the algal carbon from phytoplankton (CALG). The CALG was estimated by considering that 1 mg of Chl-a = 40 mg of CALG [14,46]. To obtain the autochthonous proportion, the CALG to POC (CALG /POC) and TSS (CALG /TSS) proportions were calculated as percentages. The latter was used to determine the seasonal change in POC sources, along with the POC/TSS. The POC/PN ratio was used to provide insight into the relative signals of decomposition or type of particulate matter [47].
The mass inflow–outflow balance of POC in the lower basin was calculated according to [48] using the difference between the M4 site (input), which is the last site of the middle basin, and the two Usumacinta River mouths (LU and LPP, output) at the lower basin. Thus, the difference in POC fluxes between the outputs and inputs determines whether the zone stores (sink, positive values) or supplies (source, negative values) POC, but without any implication regarding the process or fate of the C [19].
The seasonal differences in POC concentrations, POC flux, and POC/PN were evaluated using t-tests and non-parametric Mann–Whitney U-tests; U-tests were performed using groups with the same frequency distribution shape. The spatial variability of water quality was examined using principal component analysis (PCA). The PCA was performed using a correlation matrix. The variables were standardized to find the unit sample variance, which is the sample mean subtracted from values and divided by the standard deviation. We computed the distance to the river mouth with a map projection using the WGS84 datum. The limits of the catchments and lengths of the sampling sites were determined with digital elevation models from the USGS dataset [49] using the software QGIS Geographic Information System [50] with GRASS (v. 7.8.5). Non-parametric Spearman rank-order correlation (ρ) and linear regressions were performed to detect associations between the POC and water quality variables. All the statistical procedures were carried out in Sigmaplot version 14.0. An alpha level of 0.05 was used for all statistical tests.

3. Results

3.1. River Water Quality Dynamics

The Usumacinta River had significant seasonal changes in Q and maximum depth (Zmax), with the highest values recorded in the RS. The Q ranged between 1 and 20 L km−2 s−1 (10 ± 6 L km−2 s−1, average ± SD) in the DS and varied from 5 to 123 L km−2 s−1 (72 ± 46 L km−2 s−1) in the RS (Table 1). M3 and M4 had the most significant seasonal changes in Q and Zmax (Table 1). The middle basin sites had a higher Q than the lower basin sites due to the watershed. In the DS, the Lacandona rainforest sites had the highest Q (M1 and M2; >10 L km−2 s−1), while M2, M3, and M4 had the highest Q in the RS (>100 L km−2 s−1). The LPP site had the lowest Q in both seasons.
The total Usumacinta River Q into the Gulf of Mexico considering both Usumacinta River mouths (LU and LPP) and the Grijalva River (LG) was 1037 m3 s−1 in the DS and 5695 m3 s−1 in the RS. The Usumacinta River (LU plus LPP) contributed 49% and 54% to the total Q in the DS and RS, while the Grijalva contributed 51% and 46%, respectively. However, the Grijalva had a slightly higher watershed-normalized Q than the Usumacinta in both seasons, with 9 L km−2 s−1 in the DS and 46 L km−2 s−1 in the RS against 7 and 44 L km−2 s−1, respectively.
The water temperature at all the sites was significantly higher in the DS than in the RS (t-test; p < 0.001). The DO was similar in both seasons; it ranged from 3.3 to 7.2 mg L−1 in the DS (5.5 ± 1.7 mg L−1) and from 1.4 to 8.5 mg L−1 in the RS (5.2 ± 2.3 mg L−1). The DO decreased from the uppermost site toward the river mouth in both seasons; the middle basin sites averaged above 6 mg L−1, while the lower basin averaged concentrations below 5 mg L−1. The pH in the Usumacinta River ranged from 6.9 to 8.1; it rose from 6.9 (middle) to 7.3 (low basin) in the DS and decreased from 8.1 (middle) to 7.5 (low basin) in the RS (Table 1).
Seawater intrusion influenced the K25 values in the DS due to a conspicuous saltwater wedge detected in the lower basin sites’ vertical profiles. The K25 in the RS was lower than in the DS. In the middle basin sites, K25 did not exceed 1000 µS cm−1 in both seasons (Table 1). Overall, the Usumacinta River behaved as an oxidizing environment. Seasonally, the river had lower ORP values in the DS (174 ± 26 mV (142 to 208 mV)) than in the RS (390 ± 31 mV (355 to 428 mV)). Turb ranged from 6 to 132 NTU (21 ± 29 NTU) in the DS. In this season, LPP noticeably displayed high values and variation (Table 1). In the RS, Turb ranged from 8 to 100 NTU (55 ± 21 NTU). The average turbidity was higher in the RS in all river sites, except in M2 and LPP (Table 1).
The Chl-a concentration along the Usumacinta River varied from 0.1 to 17.6 µg L−1 (4.2 ± 4.6 µg L−1). In general, the Chl-a was lower in the RS than in the DS; in the DS, the Chl-a concentration ranged from 0.8 to 20.2 µg L−1 (6.3 ± 5.5 µg L−1), while it ranged from 0.1 to 4.4 µg L−1 in the RS (2.0 ± 1.5 µg L−1). The TSS concentration varied seasonally in the Usumacinta River, and higher values were recorded in the RS. The mean TSS concentration in the RS (77 ± 27 µg L−1) nearly doubled in the DS (28 ± 21 µg L−1; Table 1).
As shown by the PCA (Figure 2), three principal components explained 84.6% of the river water quality variance. PC 1 correlated inversely with temperature and directly with ORP and TSS. PC 2 directly correlated with the Turb, Chl-a, and K25. PC 3 was inversely correlated with DO.
The PC 1 displayed the seasonality of the Usumacinta River. In the DS, the river was less oxidant, warmer, and had lower concentrations of TSS and a lower flow rate than in the RS. The PC 1 arrangement also displayed warmer and more oxidant water, with lower TSS concentrations in the lower basin in the RS, except in M2 (Figure 2a). The PC 2 and PC 3 arrangements separated the observations into middle and lower basins in both seasons. In them, it is shown that the lower basin had a higher Chl-a concentration and salinity and less oxygenated waters than the middle basin in both seasons (Figure 2b). The LPP in the DS had very high Turb, K25, and Chl-a concentrations (Table 1).

3.2. POC Relationship with Water Quality Variables

The POC concentration was positively correlated with the TSS concentration (linear regression, p < 0.001; Figure 3) and positively correlated with Turb (ρ = 0.39, p = 0.016). In contrast, water discharge was not correlated with POC (ρ = −0.16, p = 0.333), although TSS was positively correlated with Q (ρ = 0.47, p < 0.05). In the DS, the POC concentration was positively correlated with the Chl-a concentration (ρ = 0.83, p < 0.001).
The percentage of POC in the total particulate constituents (POC/TSS) varied from 1.3 to 11.1%. The DS recorded higher percentages of POC compared with the RS (t-test, p < 0.001). There was a non-linear relationship between the POC/TSS and the TSS concentration (Figure 3). M2 showed the highest POC/TSS (7.8%) in the RS due to having the lowest TSS concentration (16 mg L−1) compared to the rest of the sites (> 60 mg L−1; Table 1). The percentage of algal biomass of TSS (CALG/TSS) varied from 0.01 to 3.7%. It also exhibited a lower percentage in the RS, displaying a non-linear relationship with the TSS concentration (Figure 3).
The CALG/POC averaged 11.6 ± 10.2%. It averaged 17.7 ± 10.3% in the DS and 5.1 ± 4.6% in the RS. Although the CALG/POC was higher in the DS compared to in the RS (U-test, p < 0.001), it remained below 50% (Figure 4) in both seasons. The algal contribution to the POC was higher in the lower basin in both seasons (DS: t-test, p < 0.05; RS: U-test, p < 0.001; Figure 4).
The POC/PN ratio averaged 10 ± 3 for the DS and 13 ± 17 for the RS, but there were no significant differences between the seasons (U-test, p > 0.05). The POC/PN ratios in the DS were different in the middle and lower basins (Figure 5).

3.3. Variation of POC Concentration and Yield

The Usumacinta River’s POC concentration ranged from 0.48 to 4.7 mg L−1 (1.69 ± 1.03 mg L−1). Overall, it was significantly higher in the RS (2.06 ± 1.14 mg L−1) than in the DS (1.42 ± 0.86 mg L−1; U-test, p > 0.01). While the POC concentrations in M3, LU, and LPP remained similar in both seasons (t-test, p > 0.05), the POC concentrations in M1, M2, M4, and LG were between 50 and 80% higher in the RS (Figure 6).
In both seasons, the POC followed the same dynamic throughout the river. In the middle basin, the POC diminished downstream from M1 (>2.2 mg L−1) to the main channel (0.48 to 1.24 mg L−1 in M3 and M4). M2 showed similar lower POC values compared to the main channel (<1.63 mg L−1). In the lower basin, the POC increased, especially in the LPP site, which exhibited concentrations above 1.72 mg L−1.
The POC yield throughout the Usumacinta River ranged between 0.2 and 19.6 kg km2 d−1 (5.5 ± 6.1 kg km2 d−1), with lower values in the DS than in the RS (t-test, p < 0.05). The POC yield increased from approximately four times to thirty times in the RS (Table 2). In the DS, only the sampling sites located in the Lacandona rainforest (M1 and M2) surpassed 1 kg km2 d−1. The yield in the RS always exceeded 1 kg km2 d−1 across all sampling sites, but clearly decreased downriver (Figure 7). While in the DS, the Usumacinta River POC yields were similar to the Grijalva River values, in the RS, the high POC concentration of the Grijalva River resulted in a higher POC yield (4.3 ± 1.2 kg km2 d−1 in the Usumacinta versus 9.5 ± 1.5 kg km2 d−1 in the Grijalva).

3.4. Mass Balance of POC in the Usumacinta Lower Basin and POC Export to the Southern Gulf of Mexico

In the DS, the POC flux increased (3.7 times) from the middle basin (M4) to the lower basin (LU plus LPP). This represented an increase in POC of approximately 70% in the lower Usumacinta River basin. In contrast, the POC mass inflow–outflow balance in the RS was −114.9 t C d−1. Hence, the lower basin retained near one-fifth of the incoming POC flux from the middle basin (494 t C d−1; Table 3).
In the DS, the POC export to the southern Gulf of Mexico was 109.1 t C d−1 and it ranged from 95.7 to 124.5 t C d−1. The RS was 926.1 t C d−1, which was almost nine times higher than that in the DS (Table 3), and it varied from 676.8 to 1080.4 t C d−1. The relative contribution to the total export of each of the rivers located in the lower basin was similar in both seasons. In the RS and DS, the Grijalva River (LG) contributed 44% and 50%, the Usumacinta River contributed 40% and 33%, and the San Pedro–San Pablo River (LPP) contributed 16% and 8%, respectively.

4. Discussion

4.1. Seasonal Hydrologic Differences in the Usumacinta River Basin and its Implication in POC Sources

Water discharge in the Usumacinta River revealed a strong seasonal hydrological contrast during rainy and dry periods. The Usumacinta River discharge to the Gulf of Mexico increased by five times in the RS compared to the DS. This increase is noteworthy compared to that of other large rivers, such as the Amazon, Congo, Tana, Yellow, and Mississippi, which exhibited lower seasonal differences from 1.2 to 3.2 times [44,51,52,53,54].
The relationship observed between TSS and water discharge clearly indicates seasonal erosion by runoff, as cited in other basins worldwide [6,55]. The cyclonic storms occurring in the Usumacinta–Grijalva basin during the RS promote the episodic movement of large volumes of organic and mineral soil through the watersheds, resulting in hyper-concentrated and highly turbid flows in the main channel of the Usumacinta River (>5% sediment volume [36]).
TSS and turbidity were positively correlated with the POC concentrations, but surprisingly not with Q. Our results showed the positive correlations between two but not all three variables (Q, TSS, and POC) commonly found in fluvial systems, which usually indicate seasonal POC inputs from terrestrial allochthonous sources (e.g., soil erosion; [10,53,56,57]). Therefore, the correlation between POC and Q in the Usumacinta basin was possibly affected by the dilution effect between the river sites. Likewise, the moderate difference between seasons regarding the POC concentration in the lower basin, which was probably associated with the algal activity, also affected the POC and Q relationship.
The non-linear relationship found in the Usumacinta basin between POC/TSS and the TSS’s concentration has been observed in other studies worldwide (e.g., [5,58]). This relationship typically results from two different hydrodynamic processes for POC supply: (a) the RS is a period of high soil erosion—allochthonous POC—during flooding events (e.g., topsoil with higher POC/TSS to deep layers with low POC/TSS), (b) the low water level and low turbidity recorded during the DS promote high aquatic primary productivity—autochthonous POC (high POC/TSS) [45]. Most studies in tropical rivers using both basin-scale approaches (e.g., [53,59,60]) and river-site-specific sampling (e.g., [10]) mention the positive effect of algal productivity in the POC supply process.
On the one hand, an intensification of algal development occurred in the DS since the POC was directly related to Chl-a; the CALG/POC increased in all sampling sites and we also found that CALG/TSS explained a significant part of the variation in POC/TSS (≈ 40%, ρ = 0.625, p = 0.004). Additionally, the sampling sites with high POC/TSS but low CALG values suggest high amounts of terrestrial carbon (as litter or soil organic matter) arriving from runoff to the rivers. On the other hand, the POC dilution could be explained by (a) the addition of soil-mineral-derived suspended solids during the RS for runoff increases, and (b) the low POC/TSS values that are registered in the Usumacinta River. Similar findings have been cited in other tropical rivers (e.g., [15,18]). For example, some Asian rivers draining through highly eroded regions exhibited a POC/TSS ranging between 0.2 and 3.6%, reflecting, in the same way, the terrestrial sedimentary origin [13,57,61].
Our results show that the contrasting seasonal hydrology in the Usumacinta River strongly impacted the TSS and POC, as well as in its sources (allochthonous versus autochthonous). The inter-tropical convergence zone oscillation largely influences the climatic seasonality and rainfall pattern in the tropics [62]. However, despite the climate, the hydrologic seasonality in tropical rivers is not a rule [63]. The extreme seasonality in the Usumacinta basin in comparison with other tropical rivers may reflect additional forcing agents, such as tropical cyclones, which drive outstanding seasonal changes in the hydrologic conditions and sediment transport [36]. This intra-seasonal variability may also affect the POC supply to the river, but it needs to be explored in future research by means of more detailed sampling campaigns over the rainy season. Since tropical storm events and interannual climatic oscillations, such as ENSO and NAO, disproportionately affect precipitation and sediment transport [31,34], their effects on POC fluxes should be addressed in the next few years.

4.2. Spatial Variation of POC Sources

As expected, the middle basin of the Usumacinta River received more allochthonous POC than the lower basin because it is dominated by high-relief forested landscapes. On average and according to global data, POC/PN ratios greater than 15 usually indicate terrestrial vegetation sources [64] or poorly degraded plant constituents in tropical soils (e.g., tree litterfall and foliage [65]). Hence, the highest POC/PN and lowest CALG/POC ratios recorded in the uplands of the basin, which correspond to the most forested section (M1 and M2), during the RS may imply that terrestrial vegetation debris constituted the principal POC source in the Usumacinta River. These findings differ to some extent from those found in other tropical rivers flowing in highly productive forests combined with monsoon climate and steep slopes. While other studies indicate that the bulk of POC comes mainly from the topsoil (matured organic matter, POC/PN ratio between 6 and 15, e.g., [15,66]), here it seems that the fresh vegetation litter would be the principal POC source in forested areas. This aspect is probably associated with the high carbon storage capacity of rainforest soils [67].
It is possible that with the increase in primary producers in the DS that the allochthonous POC in the middle Usumacinta basin mixed with autochthonous sources. In the middle basin, the observed POC/PN ratio decrease in the DS may reflect the mixing of vascular (terrestrial) and non-vascular (phytoplankton) plants [68]. Moreover, the increase in the algal proportion in POC could also support this finding.
The pristine reference site M2 did not show seasonal variation in the TSS and turbidity. The trapping effect of the well-preserved Lacandona rainforest ecosystem could be responsible for reducing rainfall erosivity during storm events. The low concentrations of suspended solids in rivers that were least affected by farming suggest that the TSS had a significant plant litter fraction arriving from the soil surface [69]. Alternatively, since the POC/TSS ratio was high in both seasons, the falling litter of overhanging riparian vegetation may continuously in M2.
The increase in the algal proportion within the POC in the lower basin sites denotes a significant autochthonous contribution in both seasons. The POC/PN ratio also indicates an autochthonous contribution in the lower basin for both seasons since the ratios (<10) match the values reported for algae in diverse aquatic environments [65,69,70,71]. The low turbidity and water discharge, high residence times, and higher water temperatures favor algal development in the lower Usumacinta basin sampling sites. Other tropical rivers have higher algal growth under these environmental conditions during the drier stages [72]. Primary production in floodplains or lentic-like environments also enhances autochthonous organic matter inputs to the rivers [73]; for instance, the lower-basin floodplains of the Orinoco River display nutrient-rich waters fostering higher algal contributions [74]. There was an increase in Chl-a concentration in the dry season in floodplain wetlands and lakes connected to the Usumacinta River in the lower basin because of the high nutrient availability and large residence time [41]. Commonly, the autochthonous POC produced in floodplains during dry seasons reaches the rivers during wet seasons [12,75]. Thus, the high algal contribution observed in both seasons in this study could alternate between the river production and floodplain contribution.
In summary, the fluvial POC in the Usumacinta River was mostly allochthonous. The particulate organic matter arriving at the river differs seasonally in the forested section of the basin from a mixture composed of terrestrial vegetation and algal particles in the DS to one composed of fresh terrestrial vegetation with soil organic matter in the RS. The seasonality of the organic matter sources likely has implications for the Usumacinta River metabolism. Organic matter derived from soils is usually enriched in refractory C [76,77], while algae have a large proportion of labile substances, such as proteins and carbohydrates [78]. Thus, in the Usumacinta River, carbon metabolism is fueled in the DS by inputs of fresh labile organic matter, which is more common in high water temperatures; by contrast, the more recalcitrant nature of organic matter reaching the river during the RS hinders the detrital pathway.

4.3. POC Concentration and Yield Variabilities in the Usumacinta River Basin

The POC concentration and flux increase observed in the middle basin during the RS represent the allochthonous sources delivered through runoff. The C litterfall from tropical forests transported toward the rivers constitutes a significant portion of the total fluvial POC transport, particularly in the RS in forested headwaters compared with plateau streams [79]. Tropical forested areas are major contributors to organic carbon flux with regard to other landscape units [67]. Litterfall from the riparian vegetation could also be exceptionally high [10].
Downstream of the tributary M1, the POC concentration decreased in both seasons because of a dilution effect caused by the tributaries. Some tributaries at the middle basin exhibited higher TSS and Q loads but a low C content. Although the flow of the Chixoy River is restricted by a hydroelectric dam built in the Guatemalan highlands, which should retain a considerable amount of fluvial sediments from the upper basin, the sediment flux toward the Usumacinta River is still relatively high [80]. Other factors, such as land-use changes (e.g., deforestation and agriculture), strongly decrease the soil organic matter and increase soil erosion and sediment transport to the La Pasion and Chixoy rivers. The loss of forest cover in the Usumacinta middle basin is more intense than in the upper and lower basins [38].

4.4. POC Export to the Gulf of Mexico

According to the mass balance and different to what was expected, the POC flux in the lower basin depended on primary producers and the hydrology dynamics. Therefore, the lower basin during the DS acted as a POC source due to the high autochthonous contribution, while in the rainy season it acted as a POC sink. The flooding occurring during the RS overflowed the main channel of the Usumacinta River, and a large fraction of the POC settled in the adjacent floodplains. During the RS, the flow connectivity with adjacent lakes and wetlands increases strongly in the lower Usumacinta River [22], as is commonly observed in other seasonal systems [81,82]. Indeed, the POC retention in the lower basin accounts for 23% of the middle basin delivery.
POC inputs enhance the mineralization in floodplains and the algae increase primary production through the use of recycled organic matter [41]; subsequently, the deposited allochthonous POC becomes new autochthonous matter. This fresh organic matter is possibly incorporated into the river during the DS when the water level decreases and the flow returns to the main channel. For instance, the lower basin retained 142.8 t of the allochthonous POC derived from the middle basin. Subsequently, 15.5 t of the POC was incorporated into the river again as autochtho-nous matter in the DS. Other tributaries (e.g., San Pedro River) draining through ex-tensive wetlands in the lower basin, supply autochthonous POC that adds to the final export of the Usumacinta River to the Gulf of Mexico.
Our results reveal that phytoplankton growth represents between 9 and 18% of the POC produced in the lower basin. This fraction concurs with the addition of organic matter by aquatic producers in the lower Amazon River (17% [83]), reaching 12% in the floodplain wetlands connected to the main river [18].
The POC export from the tropical Usumacinta River to the Gulf of Mexico is highly seasonal. The POC export increased by almost nine times in the RS with regard to the DS, which is larger than that reported for other temperate, larger rivers, such as the Mississippi River in the United States (1.6 to 5 times [54]; Table 4), but closer to that of other tropical monsoonal rivers, such as the Western Ghats in India (12 times [60]). However, the POC export from the Usumacinta River was 35 times lower than that of the Amazon river export [18]. These values are consistent due to the fact that the Amazon River has the highest water discharge (6300 km3 yr−1) and drainage area (6300 × 10−3 km2) [84]. The POC export in the Usumacinta River (Q = 147 km3 yr−1) is similar to that of the Maroni River at the border between French Guiana and Suriname (Table 4), which has an annual Q of 63.7 km3 yr−1, a smaller drainage area (0.07 × 1011 vs. 0.11 × 1011 km2), and high forest cover [85]. Additionally, being smaller in terms of water discharge (~5 times) and drainage area (~30 times) than the Mississippi River, the Usumacinta River showed a noteworthy proportional POC delivery to the Gulf of Mexico. For instance, although the POC export calculated in this study was lower than the Mississippi River POC transport (~1013 t C d−1 in the DS and ~5287 t C d−1 in the RS; [54]), it represented up to ~10% of the Mississippi River daily POC export in the DS and ~20% in the RS.
The POC concentration observed in the Usumacinta River is similar to or higher than that of other tropical rivers, such as the Maroni River [85] and Oyapock River, Senegal, Africa [12]. The Usumacinta River’s POC concentration nearly doubled the range of the Orinoco River, Venezuela, and Colombia, despite the water discharge differences (Table 4). Differently, the Usumacinta River POC concentration range is lower than that of the whitewater Amazon River tributaries [18]; the Mississippi River [54]; and some highly turbid tropical rivers with high TSS loads (>100 mg L−1), such as the Tana [86] and Yellow rivers [87].
Table 4. Water discharge (Q), POC concentrations, and flux of some rivers in the world.
Table 4. Water discharge (Q), POC concentrations, and flux of some rivers in the world.
RiverCountryZonePOCQRef.
(mg L−1)(t C d−1)(m3 s−1)
UsumacintaMexicoTrop0.48–4.7109–92623–5924This study
MaroniGuyanaTrop0.45–4.5617–1266178–4634a
OyapockGuyanaTrop0.64–3.169–279141–2564a
OrinocoVenezuelaTrop0.17–2.29200–7500500–70,000b
AmazonBrazilTrop0.27–26.84060–33,52328,000–209,000c
SenegalAfricaTrop0.2–4 ≤1700d
TanaKenyaTrop0.23–119.8 123–208e
MississippiUSATemp1.9–9.71013–528713,196–29,733f
YellowChinaTemp4.6–92.4 g
Trop: tropical; Temp: temperate. a [85]; b [88]; c [18]; d [12]; e [53]; f [54]; g [87]. POC and Q ranges were obtained from major tributaries and the main channel; POC flux at location closest to the river mouth.
The northern Gulf of Mexico receives an important amount of particulate matter and dissolved organic carbon, nutrients, and freshwater from the Mississippi River [89], leading to the development of hypoxic marine zones [90]. Hypoxia is an environmental phenomenon where the concentration of dissolved oxygen in the water column decreases to a level that can no longer support living aquatic organisms. Hypoxic areas, or “dead zones,” are likely to develop in the southern Gulf of Mexico due to the POC, nutrient, and freshwater delivery of the Usumacinta and Grijalva Rivers. The Usumacinta River mouths are highly eutrophicated (mesotrophic to eutrophic) in the upper water column, which consumes most oxygen derived from phytoplankton, causing lower oxygen saturations in coastal waters [91]. The information on the hypoxia and potential dead zone formation in the Usumacinta River influence zone should not be discarded and more investigatory studies need to be carried out on this topic [92].

5. Conclusions

The analysis of the seasonal dynamics of fluvial POC sources, concentration, and flux in the Usumacinta River basin relies on spatial differences and mechanisms that were determined using hydrogeomorphic features at the basin scale. The middle basin (steep slopes, forested) displayed a higher seasonal effect on POC concentration and source due to the potential erosion of terrestrial inputs compared to the lower basin (plateau, wetlands, and floodplains). On the other hand, the algal-derived autochthonous POC seemed to lead to a steady concentration throughout seasons in the lower basin. The POC flux difference between seasons, which was high compared with that of other large river basins, was in agreement with the hydrological contrast in the basin. The “Pantanos de Centla” wetlands located in the lowlands displayed a crucial alternate role in the Usumacinta River’s POC discharge into the Gulf of Mexico. The wetlands shifted between a sink zone of allochthonous POC in the RS (almost one-fifth of the POC transported downstream) and a source of algal fresh POC in the DS (nearly 70% of the final POC discharge). This seasonal effect caused around 20% of the total POC export to the ocean. The conservation of the large floodplains and wetlands covering most of the lower basin of the Usumacinta River is crucial for regulating the flow of POC to the southern Gulf of Mexico and maintaining a supply of labile organic matter to sustain marine productivity.

Author Contributions

Conceptualization, D.C.-L., J.A., D.C.-G., I.F.S.-R. and L.A.O.; Data curation, D.C.-L., D.C.-G. and I.F.S.-R.; Formal analysis, D.C.-L., J.A., D.C.-G., I.F.S.-R., F.G.-O. and S.S.-C.; Funding acquisition, J.A.; Investigation, D.C.-L., J.A., D.C.-G., I.F.S.-R. and L.A.O.; Methodology, D.C.-L., J.A., D.C.-G., I.F.S.-R. and L.A.O.; Project administration, J.A.; Resources, J.A.; Software, D.C.-L., D.C.-G. and I.F.S.-R.; Supervision, J.A.; Validation, D.C.-L., J.A., D.C.-G., I.F.S.-R., F.G.-O., S.S.-C. and L.A.O.; Writing—original draft, D.C.-L., J.A., D.C.-G., I.F.S.-R. and L.A.O.; Writing—review and editing, D.C.-L., J.A., D.C.-G., I.F.S.-R., F.G.-O., S.S.-C. and L.A.O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the FORDECYT-CONACYT project 273646 “Fortalecimiento de las capacidades científicas y tecnológicas para la gestión territorial sustentable de la Cuenca del Río Usumacinta y su Zona Marina de Influencia (CRUZMI), así como su adaptación ante el cambio climático” and UNAM-PAPIIT project IN216818 “Flujos de carbono, nutrientes y sedimentos en un sistema lótico tropical.” This research was also financially supported by the Programa de Investigación en Cambio Climático (PINCC 2020-2021) through project “Cuerpos acuáticos epicontinentales: papel en la dinámica del carbono y emisiones de gases de efecto invernadero en México.”

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author (J.A.) upon reasonable request.

Acknowledgments

We would like to thank the “Posgrado en Ciencias del Mar y Limnología,” UNAM, and CONACYT grants awarded to I.F.S.-R., D.C.-G. and J.D.C.-L. Jorge Ramírez and Julio Diaz helped with fieldwork and data collecting and processing. Natura y Ecosistemas Mexicanos AC foundation supported with logistic services at the Chajul Biological Station.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cole, J.J.; Prairie, Y.T.; Caraco, N.F.; McDowell, W.H.; Tranvik, L.J.; Striegl, R.G.; Duarte, C.M.; Kortelainen, P.; Downing, J.A.; Middelburg, J.J.; et al. Plumbing the Global Carbon Cycle: Integrating Inland Waters into the Terrestrial Carbon Budget. Ecosystems 2007, 10, 172–185. [Google Scholar] [CrossRef] [Green Version]
  2. Tranvik, L.J.; Cole, J.J.; Prairie, Y.T. The study of carbon in inland waters-from isolated ecosystems to players in the global carbon cycle. Limnol. Oceanogr. Lett. 2018, 3, 41–48. [Google Scholar] [CrossRef]
  3. Aufdenkampe, A.K.; Mayorga, E.; Raymond, P.A.; Melack, J.M.; Doney, S.C.; Alin, S.R.; Aalto, R.E.; Yoo, K. Riverine coupling of biogeochemical cycles between land, oceans, and atmosphere. Front. Ecol. Environ. 2011, 9, 53–60. [Google Scholar] [CrossRef] [Green Version]
  4. Schlesinger, W.H.; Melack, J.M. Transport of organic carbon in the world’s rivers. Tellus 1981, 33, 172–187. [Google Scholar] [CrossRef]
  5. Ludwig, W.; Probst, J.-L.L.; Kempe, S. Predicting the oceanic input of organic carbon by continental erosion. Global Biogeochem. Cycles 1996, 10, 23–41. [Google Scholar] [CrossRef] [Green Version]
  6. Li, M.; Peng, C.; Wang, M.; Xue, W.; Zhang, K.; Wang, K.; Shi, G.; Zhu, Q. The carbon flux of global rivers: A re-evaluation of amount and spatial patterns. Ecol. Indic. 2017, 80, 40–51. [Google Scholar] [CrossRef]
  7. Galy, V.; Peucker-Ehrenbrink, B.; Eglinton, T. Global carbon export from the terrestrial biosphere controlled by erosion. Nature 2015, 521, 204–207. [Google Scholar] [CrossRef] [PubMed]
  8. Huang, T.-H.; Fu, Y.-H.; Pan, P.-Y.; Chen, C.-T.A. Fluvial carbon fluxes in tropical rivers. Curr. Opin. Environ. Sustain. 2012, 4, 162–169. [Google Scholar] [CrossRef] [Green Version]
  9. Meybeck, M. Riverine transport of atmospheric carbon: Sources, global typology and budget. Water Air Soil Pollut. 1993, 70, 443–463. [Google Scholar] [CrossRef]
  10. Bouillon, S.; Yambélé, A.; Spencer, R.G.M.; Gillikin, D.P.; Hernes, P.J.; Six, J.; Merckx, R.; Borges, A.V. Organic matter sources, fluxes and greenhouse gas exchange in the Oubangui River (Congo River basin). Biogeosciences 2012, 9, 2045–2062. [Google Scholar] [CrossRef] [Green Version]
  11. Bouchez, J.; Galy, V.; Hilton, R.G.; Gaillardet, J.; Moreira-Turcq, P.; Pérez, M.A.; France-Lanord, C.; Maurice, L. Source, transport and fluxes of Amazon River particulate organic carbon: Insights from river sediment depth-profiles. Geochim. Cosmochim. Acta 2014, 133, 280–298. [Google Scholar] [CrossRef]
  12. Mbaye, M.L.; Gaye, A.T.; Spitzy, A.; Dähnke, K.; Afouda, A.; Gaye, B. Seasonal and spatial variation in suspended matter, organic carbon, nitrogen, and nutrient concentrations of the Senegal River in West Africa. Limnologica 2016, 57, 1–13. [Google Scholar] [CrossRef]
  13. Ziegler, A.D.; Benner, S.G.; Kunkel, M.L.; Phang, V.X.H.; Lupascu, M.; Tantasirin, C. Particulate carbon and nitrogen dynamics in a headwater catchment in Northern Thailand: Hysteresis, high yields, and hot spots. Hydrol. Process. 2016, 30, 3339–3360. [Google Scholar] [CrossRef]
  14. Shi, J.; Wang, B.; Wang, F.; Peng, X. Sources and fluxes of particulate organic carbon in the Wujiang cascade reservoirs, southwest China. Inland Waters 2018, 8, 141–147. [Google Scholar] [CrossRef]
  15. Lloret, E.; Dessert, C.; Lajeunesse, E.; Crispi, O.; Pastor, L.; Gaillardet, J.; Benedetti, M.F. Tropical oceanic islands carbon export Are small mountainous tropical watersheds of oceanic islands important for carbon export? Tropical oceanic islands carbon export. Biogeosci. Discuss. 2012, 9, 7117–7163. [Google Scholar] [CrossRef]
  16. Hedges, J.I.; Clark, W.A.; Quay, P.D.; Richey, J.E.; Devol, A.H.; Santos, M. Compositions and fluxes of particulate organic material in the Amazon River. Limnol. Oceanogr. 1986, 31, 717–738. [Google Scholar] [CrossRef] [Green Version]
  17. Richey, J.E.; Hedges, J.I.; Devol, A.H.; Quay, P.D.; Victoria, R.; Martinelli, L.; Forsberg, B.R. Biogeochemistry of carbon in the Amazon River. Limnol. Oceanogr. 1990, 35, 352–371. [Google Scholar] [CrossRef]
  18. Moreira-Turcq, P.; Seyler, P.; Guyot, J.L.; Etcheber, H. Exportation of organic carbon from the Amazon River and its main tributaries. Hydrol. Process. 2003, 17, 1329–1344. [Google Scholar] [CrossRef]
  19. Geeraert, N.; Omengo, F.O.; Tamooh, F.; Marwick, T.R.; Borges, A.V.; Govers, G.; Bouillon, S. Seasonal and inter-annual variations in carbon fluxes in a tropical river system (Tana River, Kenya). Biogeosci. Discuss. 2017, 14, 1–27. [Google Scholar] [CrossRef] [Green Version]
  20. Benke, A. Streams and Rivers of North America: Western, Northern and Mexican Basins. In Encyclopedia of Inland Waters; Likens, G.E., Ed.; Elsevier: Boston, MA, USA, 2009; pp. 414–424. ISBN 0120882531. [Google Scholar]
  21. Hudson, P.F.; Hendrickson, D.A.; Benke, A.C.; Varela-Rmero, A.; Rodiles-Hernández, R.; Minckley, W.L. Rivers of Mexico. In Rivers of North America; Benke, A., Cushing, C., Eds.; Elsevier: Burlington, MA, USA, 2005; pp. 1030–1084. [Google Scholar]
  22. Tapia-Silva, F.-O.; Contreras-Silva, A.-I.; Rosales-Arriaga, E.-R. Hydrological characterization of the Usumacinta River Basin towards the preservation of environmental services. ISPRS Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2015, XL-7/W3, 1505–1509. [Google Scholar] [CrossRef]
  23. Waibel, L. (Ed.) Sierra Madre de Chiapas, 1st ed.; Miguel Ángel Porrúa: Mexico City, Mexico, 1998; ISBN 968-842-716-7. [Google Scholar]
  24. Saavedra Guerrero, A.; López López, D.M.; Castellanos Fajardo, L.A. Descripción del medio físico de la cuenca media del río Usumacinta en México. In Conservación y Desarrollo Sustentable en la Selva Lacandona. 25 años de Actividades y Experiencias; Carabias, J., de La Maza, J., Cadena, R., Eds.; (coordinadores); Natura y Ecosistemas Mexicanos: Mexico City, Mexico, 2015; pp. 19–34. ISBN 9786079710200. [Google Scholar]
  25. Díaz Gallegos, J.R.; Mas, J.-F.; Velázquez, A.; Díaz-Gallegos, J.R.; Mas, J.-F.; Velázquez, A. Trends of tropical deforestation in Southeast Mexico. Singap. J. Trop. Geogr. 2010, 31, 180–196. [Google Scholar] [CrossRef]
  26. CONAGUA Programas de Medidas Preventivas y de Mitigación a la Sequía (PMPMS) por Consejo de Cuenca. Available online: https://www.gob.mx/cms/uploads/attachment/file/99961/PMPMS_CC_R_os_Grijalva_y_Usumacinta.pdf (accessed on 18 February 2020).
  27. CCGSS Servicios de Información Geoespacial Plataforma GeoWeb IDEGeo. Available online: https://idegeo.centrogeo.org.mx/layers/geonode:cuenca_17/layer_info_metadata (accessed on 18 March 2021).
  28. Plascencia Vargas, J.H.; González Espinosa, M.; Ramírez Marcial, N.; Musálem Castillejos, K. Características físico-bióticas de la cuenca del río Grijalva. In Montañas, Pueblos y Agua. Dimenciones y Realidades de la Cuenca Grijalva; González Espinosa, M., Brunel Manse, M.C., Eds.; El Colegio de la Frontera Sur: Mexico City, Mexico, 2014; Volume 1, pp. 29–79. ISBN 978-607-8429-01-1. [Google Scholar]
  29. Mendoza, A.; Soto-Cortes, G.; Priego-Hernandez, G.; Rivera-Trejo, F. Historical description of the morphology and hydraulic behavior of a bifurcation in the lowlands of the Grijalva River Basin, Mexico. CATENA 2019, 176, 343–351. [Google Scholar] [CrossRef]
  30. CONAGUA. Estadísticas del Agua en México, 2018th ed.; Secretaría de Medio Ambiente y Recursos Naturales, Comisión Nacional del Agua: Mexico City, Mexico, 2018. [Google Scholar]
  31. Arreguín-Cortés, F.I.; Rubio-Gutiérrez, H.; Domínguez-Mora, R.; de Luna-Cruz, F. Analysis of Floods in the Tabasco Plains from 1995–2010. Ecnología y Ciencias del Agua 2014, 5, 5–32. [Google Scholar]
  32. Soares, D.; García, A. La cuenca del río Usumacinta desde la Perspectiva del cambio Climático; Instituto Mexicano de Tecnología del Agua: Jiutepec, Mexico, 2017; ISBN 9786079368807. [Google Scholar]
  33. Grodsky, S.A.; Carton, J.A. The Intertropical Convergence Zone in the South Atlantic and the Equatorial Cold Tongue. J. Clim. 2003, 16, 723–733. [Google Scholar] [CrossRef] [Green Version]
  34. Muñoz-Salinas, E.; Castillo, M. Streamflow and sediment load assessment from 1950 to 2006 in the Usumacinta and Grijalva Rivers (Southern Mexico) and the influence of ENSO. Catena 2015, 127, 270–278. [Google Scholar] [CrossRef]
  35. Alcérreca-Huerta, J.C.; Callejas-Jiménez, M.E.; Carrillo, L.; Castillo, M.M. Dam implications on salt-water intrusion and land use within a tropical estuarine environment of the Gulf of Mexico. Sci. Total Environ. 2018, 652, 1102–1112. [Google Scholar] [CrossRef] [PubMed]
  36. Muñoz-Salinas, E.; Castillo, M.; Sanderson, D.; Kinnaird, T.; Cruz-Zaragoza, E. Using three different approaches of OSL for the study of young fluvial sediments at the coastal plain of the Usumacinta-Grijalva River Basin, southern Mexico. Earth Surf. Process. Landforms 2016, 41, 823–834. [Google Scholar] [CrossRef]
  37. de la Maza, J. Caracterización de la subcuenca del Lacantún. In Conservación y Desarrollo Sustentable de la Selva Lacandona; Natura y Ecosistemas Mexicanos, A.C.: Mexico City, Mexico, 2015; pp. 79–83. [Google Scholar]
  38. Peralta-Carreta, C.; Gallardo-Cruz, J.A.; Solórzano, J.V.; Hernandez-Gómez, M. Clasificación del uso de suelo y vegetación en áreas de pérdida de cobertura arbórea (2000–2016) en la cuenca del río Usumacinta. Madera y Bosques 2019, 25, 1–19. [Google Scholar] [CrossRef] [Green Version]
  39. García, A.G.; Kauffer, M.E. Las cuencas compartidas entre México, Guatemala y Belice: Un acercamiento a su delimitación y problemática general. Front. Norte 2011, 23, 131–162. [Google Scholar]
  40. Sánchez-Hernández, R.; De Dios Mendoza-Palacios, J.; De, J.C.; Reyes, C.; Mendoza Martínez, J.E.; Ramos-Reyes, R. Erosión potencial en la cuenca baja Grijalva-Usumacinta. Univ. y Cienc. 2013, 29, 153–161. [Google Scholar]
  41. Castillo, M.M. Suspended sediment, nutrients, and chlorophyll in tropical floodplain lakes with different patterns of hydrological connectivity. Limnologica 2020, 82, 125767. [Google Scholar] [CrossRef]
  42. Olivera, M. de los A.G.; Gil, G.O. Análisis de la flora de la Reserva de la Biosfera de los Pantanos de Centla, Tabasco, México. Univ. y Cienc. 2000, 15, 67–104. [Google Scholar]
  43. Arar, E.J.; Collins, G.B. Method 445.0 In Vitro Determination of Chlorophyll a and Pheophytin ain Marine and Freshwater Algae by Fluorescence; Environmental Protection Agency: Washington, DC, USA, 1997. [Google Scholar]
  44. Ran, L.; Lu, X.X.; Sun, H.; Han, J.; Li, R.; Zhang, J. Spatial and seasonal variability of organic carbon transport in the Yellow River, China. J. Hydrol. 2013, 498, 76–88. [Google Scholar] [CrossRef]
  45. Meybeck, M. Origins and Behaviors of Carbon Species in World Rivers. In Soil Erosion and Carbon Dynamics; Roose, E.J., Lal, R., Feller, C., Barthes, B., Stewart, B.A., Eds.; CRC Press: Boca Raton, FL, USA, 2005; pp. 209–238. [Google Scholar]
  46. Del Giorgio, P.A.; Gasol, J.M. Biomass Distribution in Freshwater Plankton Communities. Am. Nat. 1995, 146, 135–152. [Google Scholar] [CrossRef]
  47. Brust, G.E. Management strategies for organic vegetable fertility. In Safety and Practice for Organic Food; Elsevier: San Diego, CA, USA, 2019; pp. 193–212. ISBN 9780128120606. [Google Scholar]
  48. DelDuco, E.M.; Xu, Y.J. Dissolved carbon transport and processing in North America’s largest swamp river entering the Northern Gulf of Mexico. Water 2019, 11, 1395. [Google Scholar] [CrossRef] [Green Version]
  49. U.S. Geological Survey USGS EROS Archive—Digital Elevation—Shuttle Radar Topography Mission (SRTM). Available online: https://doi.org/10.5066/F7PR7TFT (accessed on 31 March 2021).
  50. QGIS.org QGIS Geographic Information System 2021. Available online: http://www.qgis.org (accessed on 21 March 2021).
  51. Benedetti, M.F.; Mounier, S.; Filizola, N.; Benaim, J.; Seyler, P. Carbon and metal concentrations, size distributions and fluxes in major rivers of the Amazon basin. Hydrol. Process. 2003, 17, 1363–1377. [Google Scholar] [CrossRef]
  52. Seyler, P.; Coynel, A.; Moreira-Turcq, P.; Etcheber, H.; Colas, C.; Orange, D.; Bricquet, J.-P.; Laraque, A.; Luc Guyot, J.; Olivry, J.-C.; et al. Organic Carbon Transported by the Equatorial Rivers: Example of Congo-Zaire and Amazon Basins. In Soil Erosion and Carbon Dynamics; CRC Press: Boca Raton, FL, USA, 2005; pp. 255–274. [Google Scholar]
  53. Tamooh, F.; Van Den Meersche, K.; Meysman, F.; Marwick, T.R.; Borges, A.V.; Merckx, R.; Dehairs, F.; Schmidt, S.; Nyunja, J.; Bouillon, S. Distribution and origin of suspended matter and organic carbon pools in the Tana River Basin, Kenya. Biogeosciences 2012, 9, 2905–2920. [Google Scholar] [CrossRef] [Green Version]
  54. Cai, Y.; Guo, L.; Wang, X.; Aiken, G. Abundance, stable isotopic composition, and export fluxes of DOC, POC, and DIC from the Lower Mississippi River during 2006-2008. J. Geophys. Res. Biogeosciences 2015, 120, 2273–2288. [Google Scholar] [CrossRef] [Green Version]
  55. Webster, J.R.; Benfield, E.F.; Golladay, S.W.; Hill, B.H.; Hornick, L.E.; Kazmierczak, R.F.; Perry, W.B. Experimental studies of physical factors affecting seston transport in streams1. Limnol. Oceanogr. 1987, 32, 848–863. [Google Scholar] [CrossRef]
  56. Zou, J. Geochemical characteristics and organic carbon sources within the upper reaches of the Xi River, southwest China during high flow. J. Earth Syst. Sci. 2017, 126, 6. [Google Scholar] [CrossRef] [Green Version]
  57. Tao, S.; Eglinton, T.I.; Zhang, L.; Yi, Z.; Montluçon, D.B.; McIntyre, C.; Yu, M.; Zhao, M. Temporal variability in composition and fluxes of Yellow River particulate organic matter. Limnol. Oceanogr. 2018, 63, S119–S141. [Google Scholar] [CrossRef]
  58. Ittekkot, V. Global trends in the nature of organic matter in river suspensions. Nature 1988, 332, 436–438. [Google Scholar] [CrossRef]
  59. Coynel, A.; Seyler, P.; Etcheber, H.; Meybeck, M.; Orange, D. Spatial and seasonal dynamics of total suspended sediment and organic carbon species in the Congo River. Global Biogeochem. Cycles 2005, 19, 1–17. [Google Scholar] [CrossRef]
  60. Kiran, S.; Reddy, K.; Gupta, H.; Badimela, U.; Reddy, D.V.; Kurakalva, R.M.; Kumar, D. Export of particulate organic carbon by the mountainous tropical rivers of Western Ghats, India: Variations and controls. Sci. Total Environ. 2021, 751, 142115. [Google Scholar] [CrossRef]
  61. Tao, F.-X.; Liu, C.-Q.; Li, S.-L. Source and flux of POC in two subtropical karstic tributaries with contrasting land use practice in the Yangtze River Basin. Appl. Geochemistry 2009, 24, 2102–2112. [Google Scholar] [CrossRef]
  62. Beck, H. Tropical ecology. Encycl. Ecol. 2018, 2, 671–678. [Google Scholar] [CrossRef]
  63. Syvitski, J.P.M.; Cohen, S.; Kettner, A.J.; Brakenridge, G.R. How important and different are tropical rivers?-An overview. Geomorphology 2014, 227, 5–17. [Google Scholar] [CrossRef]
  64. Kendall, C.; Silva, S.R.; Kelly, V.J. Carbon and nitrogen isotopic compositions of particulate organic matter in four large river systems across the United States. Hydrol. Process. 2001, 15, 1301–1346. [Google Scholar] [CrossRef]
  65. Cleveland, C.C.; Liptzin, D. C:N:P stoichiometry in soil: Is there a “Redfield ratio” for the microbial biomass? Biogeochemistry 2007, 85, 235–252. [Google Scholar] [CrossRef]
  66. Bird, M.I.; Robinson, R.A.J.; Oo, N.W.; Aye, M.M.; Lu, X.X.; Higgitt, D.L.; Swe, A.; Tun, T.; Lhaing Win, S.; Sandar Aye, K.; et al. A preliminary estimate of organic carbon transport by the Ayeyarwady (Irrawaddy) and Thanlwin (Salween) Rivers of Myanmar. Quat. Int. 2008, 186, 113–122. [Google Scholar] [CrossRef]
  67. Navarrete-Segueda, A.; Martínez-Ramos, M.; Ibarra-Manríquez, G.; Vázquez-Selem, L.; Siebe, C. Variation of main terrestrial carbon stocks at the landscape-scale are shaped by soil in a tropical rainforest. Geoderma 2018, 313, 57–68. [Google Scholar] [CrossRef]
  68. Meyers, P.A.; Ishiwatari, R. Lacustrine organic geochemistryman overview of indicators of organic matter sources and diagenesis in lake sediments. Org. Geochem. 1993, 20, 867–900. [Google Scholar] [CrossRef] [Green Version]
  69. Weiguo, L.; Zisheng, A.; Weijian, Z.; Head, M.; Delin, C. Carbon isotope and C/N ratios of suspended matter in rivers: An indicator of seasonal change in C4/C3 vegetation. Appl. Geochem. 2003, 18, 1241–1249. [Google Scholar] [CrossRef]
  70. Meyers, P.A. Preservation of elemental and isotopic source identification of sedimentary organic matter. Chem. Geol. 1994, 114, 289–302. [Google Scholar] [CrossRef] [Green Version]
  71. Onstad, G.D.; Canfield, D.E.; Quay, P.D.; Hedges, J.I. Sources of particulate organic matter in rivers from the continental usa: Lignin phenol and stable carbon isotope compositions. Geochim. Cosmochim. Acta 2000, 64, 3539–3546. [Google Scholar] [CrossRef]
  72. Kohler, T.J.; Heatherly, T.N.; El-Sabaawi, R.W.; Zandonà, E.; Marshall, M.C.; Flecker, A.S.; Pringle, C.M.; Reznick, D.N.; Thomas, S.A. Flow, nutrients, and light availability influence Neotropical epilithon biomass and stoichiometry. Freshw. Sci. 2012, 31, 1019–1034. [Google Scholar] [CrossRef]
  73. Ward, N.D.; Bianchi, T.S.; Medeiros, P.M.; Seidel, M.; Richey, J.E.; Keil, R.G.; Sawakuchi, H.O. Where Carbon Goes When Water Flows: Carbon Cycling across the Aquatic Continuum. Front. Mar. Sci. 2017, 4. [Google Scholar] [CrossRef] [Green Version]
  74. Lewis, W.M. Physical and Chemical Features of Tropical Flowing Waters. In Tropical Stream Ecology; Dudgeon, D., Ed.; Academic Press: London, UK, 2008; pp. 1–21. ISBN 9780120884490. [Google Scholar]
  75. Moreira-Turcq, P.; Bonnet, M.P.; Amorim, M.; Bernardes, M.; Lagane, C.; Maurice, L.; Perez, M.; Seyler, P. Seasonal variability in concentration, composition, age, and fluxes of particulate organic carbon exchanged between the floodplain and Amazon River. Global Biogeochem. Cycles 2013, 27, 119–130. [Google Scholar] [CrossRef]
  76. Jobbágy, E.G.; Jackson, R.B. The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol. Appl. 2000, 10, 423–436. [Google Scholar] [CrossRef]
  77. Keller, C.K. Carbon Exports from Terrestrial Ecosystems: A Critical-Zone Framework. Ecosystems 2019, 22, 1–15. [Google Scholar] [CrossRef]
  78. Ludwig, W. Continental Erosion and River Transport of Organic Carbon to the World’s Oceans; Mémoire, 9.; Sciences Géologiques, University Louis-Pasteur: Strasbourg, France, 1997. [Google Scholar]
  79. Selva, E.C.; Couto, E.G.; Johnson, M.S.; Lehmann, J. Litterfall production and fluvial export in headwater catchments of the southern Amazon. J. Trop. Ecol. 2007, 23, 329–335. [Google Scholar] [CrossRef] [Green Version]
  80. Nooren, K.; Hoek, W.Z.; Winkels, T.; Huizinga, A.; Van der Plicht, H.; Van Dam, R.L.; Van Heteren, S.; Van Bergen, M.J.; Prins, M.A.; Reimann, T.; et al. The Usumacinta–Grijalva beach-ridge plain in southern Mexico: A high-resolution archive of river discharge and precipitation. Earth Surf. Dyn. 2017, 5, 529–556. [Google Scholar] [CrossRef] [Green Version]
  81. Noe, G.B.; Hupp, C.R. Carbon, nitrogen, and phosphorus accumulation in floodplains of Atlantic Coastal Plain rivers, USA. Ecol. Appl. 2005, 15, 1178–1190. [Google Scholar] [CrossRef]
  82. Sutfin, N.A.; Wohl, E.E.; Dwire, K.A. Banking carbon: A review of organic carbon storage and physical factors influencing retention in floodplains and riparian ecosystems. Earth Surf. Process. Landforms 2016, 41, 38–60. [Google Scholar] [CrossRef]
  83. Hedges, J.I.; Mayorga, E.; Tsamakis, E.; McClain, M.E.; Aufdenkampe, A.; Quay, P.; Richey, J.E.; Benner, R.; Opsahl, S.; Black, B.; et al. Organic matter in Bolivian tributaries of the Amazon River: A comparison to the lower mainstream. Limnol. Oceanogr. 2000, 45, 1449–1466. [Google Scholar] [CrossRef] [Green Version]
  84. Milliman, J.D.; Farnsworth, K.L. River Discharge to the Coastal Ocean: A Global Synthesis; Cambridge University Press: Cambridge, UK, 2011; ISBN 9780511781247. [Google Scholar]
  85. Gallay, M.; Mora, A.; Martinez, J.-M.; Gardel, A.; Laraque, A.; Sarrazin, M.; Beaucher, E.; Doudou, J.-C.; Lagane, C. Dynamics and fluxes of organic carbon and nitrogen in two Guiana Shield river basins impacted by deforestation and mining activities. Hydrol. Process. 2018, 32, 17–29. [Google Scholar] [CrossRef] [Green Version]
  86. Tamooh, F.; Meysman, F.J.R.; Borges, A.V.; Marwick, T.R.; Van Den Meersche, K.; Dehairs, F.; Merckx, R.; Bouillon, S. Sediment and carbon fluxes along a longitudinal gradient in the lower Tana River (Kenya). J. Geophys. Res. G Biogeosci. 2014, 119, 1340–1353. [Google Scholar] [CrossRef] [Green Version]
  87. Liu, J.; Song, X.; Wang, Z.; Yang, L.; Sun, Z.; Wang, W. Variations of carbon transport in the Yellow River, China. Hydrol. Res. 2015, 46, 746–762. [Google Scholar] [CrossRef]
  88. Mora, A.; Laraque, A.; Moreira-Turcq, P.; Alfonso, J.A. Temporal variation and fluxes of dissolved and particulate organic carbon in the Apure, Caura and Orinoco rivers, Venezuela. J. S. Am. Earth Sci. 2014, 54, 47–56. [Google Scholar] [CrossRef]
  89. Wang, X.-C.; Chen, R.F.; Gardner, G.B. Sources and transport of dissolved and particulate organic carbon in the Mississippi River estuary and adjacent coastal waters of the northern Gulf of Mexico. Mar. Chem. 2004, 89, 241–256. [Google Scholar] [CrossRef]
  90. Liu, Z.; Xue, J. The Lability and Source of Particulate Organic Matter in the Northern Gulf of Mexico Hypoxic Zone. J. Geophys. Res. Biogeosciences 2020, 125, 1–16. [Google Scholar] [CrossRef]
  91. Signoret, M.; Monreal-Gómez, M.A.; Aldeco, J.; Salas-de-León, D.A. Hydrography, oxygen saturation, suspended particulate matter, and chlorophyll-a fluorescence in an oceanic region under freshwater influence. Estuar. Coast. Shelf Sci. 2006, 69, 153–164. [Google Scholar] [CrossRef]
  92. Cardoso-Mohedano, J.G.; Canales-Delgadillo, J.C.; Machain-Castillo, M.L.; Hernández-Hernández, J.G.; Sanchez-Cabeza, J.A.; Ruiz-Fernández, A.C.; Alonso-Rodríguez, R.; Gómez-Ponce, M.A.; Esqueda-Lara, K.; Merino-Ibarra, M.; et al. Absence of hypoxia events in the adjacent coastal waters of Grijalva-Usumacinta river, Southern Gulf of Mexico. Mar. Pollut. Bull. 2020, 156, 111174. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The Usumacinta River basin and sampling sites (elaborated using the INEGI database; datum: WGS84).
Figure 1. The Usumacinta River basin and sampling sites (elaborated using the INEGI database; datum: WGS84).
Water 13 01561 g001
Figure 2. Spatial and seasonal water quality variation in the Usumacinta River main channel and major tributaries, as displayed using PCA: (a) Representation of PC 1 against PC 2 and (b) representation of PC 1 against PC 3. The percentages indicate the explained variance of the PCs.
Figure 2. Spatial and seasonal water quality variation in the Usumacinta River main channel and major tributaries, as displayed using PCA: (a) Representation of PC 1 against PC 2 and (b) representation of PC 1 against PC 3. The percentages indicate the explained variance of the PCs.
Water 13 01561 g002
Figure 3. Variation in the organic carbon contribution to the total suspended solids in the Usumacinta River during the dry (DS) and rainy (RS) seasons: (a) relationship between the POC and TSS concentrations, (b) relative content of POC in TSS as a function of the TSS concentration, and (c) relative contribution of CALG to TSS as a function of the TSS concentration.
Figure 3. Variation in the organic carbon contribution to the total suspended solids in the Usumacinta River during the dry (DS) and rainy (RS) seasons: (a) relationship between the POC and TSS concentrations, (b) relative content of POC in TSS as a function of the TSS concentration, and (c) relative contribution of CALG to TSS as a function of the TSS concentration.
Water 13 01561 g003
Figure 4. Contribution of algal biomass to POC throughout the Usumacinta River basin during the dry (DS) and rainy (RS) seasons. Bars represent averages and the whiskers represent the standard deviations.
Figure 4. Contribution of algal biomass to POC throughout the Usumacinta River basin during the dry (DS) and rainy (RS) seasons. Bars represent averages and the whiskers represent the standard deviations.
Water 13 01561 g004
Figure 5. Mean algal biomass contribution to POC in the Usumacinta River as a function of the mean POC/PN ratio. The whiskers represent the standard deviations.
Figure 5. Mean algal biomass contribution to POC in the Usumacinta River as a function of the mean POC/PN ratio. The whiskers represent the standard deviations.
Water 13 01561 g005
Figure 6. Box and whisker plot showing the seasonal and spatial distribution of the POC concentration in the Usumacinta River basin. Boxplots represent the 25th and 75th percentiles of the concentration, and the dotted line within each box represents the mean. The whiskers represent the standard deviation.
Figure 6. Box and whisker plot showing the seasonal and spatial distribution of the POC concentration in the Usumacinta River basin. Boxplots represent the 25th and 75th percentiles of the concentration, and the dotted line within each box represents the mean. The whiskers represent the standard deviation.
Water 13 01561 g006
Figure 7. Seasonal and spatial distributions of POC yield in the Usumacinta River basin. Bars represent average yields and whiskers represent the standard deviations.
Figure 7. Seasonal and spatial distributions of POC yield in the Usumacinta River basin. Bars represent average yields and whiskers represent the standard deviations.
Water 13 01561 g007
Table 1. Average (first line) and standard deviation (second line) values of water quality variables in the Usumacinta River main channel and major tributaries in the dry and rainy seasons.
Table 1. Average (first line) and standard deviation (second line) values of water quality variables in the Usumacinta River main channel and major tributaries in the dry and rainy seasons.
SeasonRiver/Site QZmaxTDOpHK25ORPTurbChl-aTSS
Dry SeasonLacantún/M1 +X−3123.527.36.36.9 *395208464.749
σ 0<0.11.12280.57
Tzendales/M2 +X−231.527.87.27.1 *931202100.914
σ 001.00<110.11
Usumacinta/M3X−3964.230.76.57.1 *706197113.012
σ 0.1<0.11.1<1810.41
Usumacinta/M4X−36920.130.76.47.2 *83916971.97
σ 0.1<0.11.23310.41
Usumacinta/LUX−43212.529.93.37.4 *33,36815388.312
σ 0.92.11.118,8061430.82
Grijalva/LG +X−5271130.34.07.4 *27,61014264.040
σ 0.82.21.215,126511.015
San Pedro–San Pablo/LPP ─X−784.829.85.17.6 *57,09814213217.654
σ 0.30.81.111485852.211
TotalX−3058.229.65.57.2 *13,019174236.328
σ1876.62.21.91.619,86626365.520
Rainy SeasonLacantún/M1 +X−9486.322.28.58.1 *3743461000.4119
σ <0.1<0.11.1<11640.114
Tzendales/M2 +X−1733.823.48.08.1 *53935580.116
σ <0.1<0.11.0<1110.0
Usumacinta/M3X−571515.624.35.87.7 *316357550.876
σ <0.10.11.15910.16
Usumacinta/M4X−593435.224.96.67.8 *359428520.771
σ <0.1<0.11.0<1420.12
Usumacinta/L1X−270917.126.74.67.7 *361409443.463
σ <0.10.11.0<1280.37
Grijalva/LG +X−261812.327.91.47.3 *337381493.382
σ <0.10.11.111090.322
San Pedro–San Pablo/LPP ─X−3684.827.51.97.4 *394390433.769
σ 0.10.21.1144150.520
TotalX−26381425.45.27.7 *360390552.078
σ2395111.92.31.73731211.525
Q: water discharge (m3 s−1); Zmax: maximum depth (m); T: temperature (°C); DO: dissolved oxygen concentration (mg L−1); K25: electrical conductivity standardized at 25 °C (µS cm−1); ORP: oxidoreduction potential (mV); Turb: turbidity (NTU); Chl-a: chlorophyll-a concentration (µg L−1); TSS: total suspended solids (mg L−1); * geometric mean. + tributary; ─ distributary.
Table 2. Average (first line) and standard deviation (second line) values of POC concentration, POC flux, POC/TSS, POC/PN, and CALG/POC ratios in the sampling sites of the Usumacinta River in the dry and rainy seasons.
Table 2. Average (first line) and standard deviation (second line) values of POC concentration, POC flux, POC/TSS, POC/PN, and CALG/POC ratios in the sampling sites of the Usumacinta River in the dry and rainy seasons.
POCPOC/TSS
(%)
POC/PNCALG/POC
(%)
SeasonRiver/Site Concentration
(mg L−1)
Flux
(t C d−1)
Yield (kg C km2 d−1)
Dry SeasonLacantún/M1 +X−2.4766.74.25.1147.6
σ0.164.40.30.620.9
Tzendales/M2 +X−0.831.71.16.1114.2
σ0.030.10.0 1
Usumacinta/M3X−0.8629.50.67.51014.0
σ0.144.70.10.311.5
Usumacinta/M4X−0.5116.40.37.41214.9
σ0.051.70.00.622.1
Usumacinta/LUX−1.1643.40.69.5929.4
σ0.186.80.11.816.2
Grijalva/LG +X−1.0648.40.82.9615.1
σ0.052.10.00.912.9
San Pedro–San Pablo/LPP ─X−2.5717.30.24.5929.9
σ0.734.90.10.6013.4
TotalX−1.4233.31.13.71016.4
σ0.8621.11.43.224.8
Rainy SeasonLacantún/M1 +X−3.78309.319.63.2450.4
σ0.6654.33.40.3110.1
Tzendales/M2 +X−1.271912.87.8270.4
σ0.040.60.4 1
Usumacinta/M3X−1.07530.711.41.422.9
σ0.1363.61.40.000.1
Usumacinta/M4X−0.96494.79.71.423.1
σ0.0523.90.50.100.5
Usumacinta/LUX−1.303044.31.9512.7
σ0.3787.11.20.515.6
Grijalva/LG +X−2.42546.39.52.925.9
σ0.3885.61.50.611.6
San Pedro–San Pablo/LPP ─X−2.3875.91.03.856.3
σ0.247.70.11.321.6
TotalX−2.06323.09.81.9104.5
σ1.14198.96.02.1142.1
+ tributary; ─ distributary.
Table 3. Input–output balances of water discharge (Q) and POC in the Usumacinta lower basin and the total POC export to the Gulf of Mexico in the rainy (RS) and dry (DS) seasons of 2017. M4 is the output of the middle basin; LU: Usumacinta River; LG: Grijalva River, LPP: San Pedro–San Pablo River. See Figure 1 for the detailed locations of the sampling sites.
Table 3. Input–output balances of water discharge (Q) and POC in the Usumacinta lower basin and the total POC export to the Gulf of Mexico in the rainy (RS) and dry (DS) seasons of 2017. M4 is the output of the middle basin; LU: Usumacinta River; LG: Grijalva River, LPP: San Pedro–San Pablo River. See Figure 1 for the detailed locations of the sampling sites.
SiteQ
(m3 s−1)
POC Flux
(t C d−1)
DSRSDSRS
InputM4369593416.4494.7
OutputLU432270943.4304.0
LPP7836817.375.9
Lower Usumacinta basin balance(input-output)M4 − (LU + LPP)−1412857−44.3114.9
Usumacinta basin exportLU + LPP510307760.7379.9
Usumacinta–Grijalva basin exportLU + LG + LPP10375695109.1926.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cuevas-Lara, D.; Alcocer, J.; Cortés-Guzmán, D.; Soria-Reinoso, I.F.; García-Oliva, F.; Sánchez-Carrillo, S.; Oseguera, L.A. Particulate Organic Carbon in the Tropical Usumacinta River, Southeast Mexico: Concentration, Flux, and Sources. Water 2021, 13, 1561. https://doi.org/10.3390/w13111561

AMA Style

Cuevas-Lara D, Alcocer J, Cortés-Guzmán D, Soria-Reinoso IF, García-Oliva F, Sánchez-Carrillo S, Oseguera LA. Particulate Organic Carbon in the Tropical Usumacinta River, Southeast Mexico: Concentration, Flux, and Sources. Water. 2021; 13(11):1561. https://doi.org/10.3390/w13111561

Chicago/Turabian Style

Cuevas-Lara, Daniel, Javier Alcocer, Daniela Cortés-Guzmán, Ismael F. Soria-Reinoso, Felipe García-Oliva, Salvador Sánchez-Carrillo, and Luis A. Oseguera. 2021. "Particulate Organic Carbon in the Tropical Usumacinta River, Southeast Mexico: Concentration, Flux, and Sources" Water 13, no. 11: 1561. https://doi.org/10.3390/w13111561

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop