Management of nitrogen nutrient flow under legume practice from La Teresa basin to the Batabanó Gulf

Nitrogen management in a karst basin is examined through the establishment of legumes Leucaena leucocephala and Phaseolus vulgaris in pasture and crop areas respectively, with the objective of enriching the nutrient flow interrupted by the construction of dams and canals built during the 1990s, which have altered the natural patterns of hydrological variation and transport to the sea. The different soil treatments show low values of overland fluxes to the sea between 0.9623 and 0.9624 kg ha−1 yr−1 for crops, discretely increasing by 0.02% in pasture and by 0.28% for crop + pasture according to the geomorphic potential of the basin. For land use conditions, a positive economic nitrogen balance was estimated in the range of 55–1560 for crops and 22–2400 USD ha−1 yr−1 for pastures, the latter showing positive balances at the lowest fixation rates, confirming that cover management with legumes maximizes gains and minimizes losses for both agricultural and marine ecosystems, providing a source of indirect fertilization downstream where agricultural areas are located. In the flat karst region, geomorphological potential is the main factor in predicting the dynamics of sediment movement. Leguminous land use has a weak effect on alleviating pressures on the island's oligotrophic marine ecosystems. In both crop and pasture, a positive annual economic balance of N was found in the range of 55–2400 USD ha−1 yr−1. In the flat karst region, geomorphological potential is the main factor in predicting the dynamics of sediment movement. Leguminous land use has a weak effect on alleviating pressures on the island's oligotrophic marine ecosystems. In both crop and pasture, a positive annual economic balance of N was found in the range of 55–2400 USD ha−1 yr−1.


Introduction
The waters of the Gulf of Batabanó are oligotrophic waters [1,2] under strong historical pressure, due to the execution of various hydro-technical works during the 1970s (reservoirs) for the development of agricultural programs [3,4], in which the Dique Sur (South Dike) stands out. These actions have notably influenced nitrogen and phosphorus inputs from the land to the sea, as these, act as an obstacle, along with the drastic reduction of fertilizer imports during the so-called Special Period in Cuba (1991)(1992)(1993)(1994)(1995)(1996) as shown in Table 1 and Fig. 1 [5,6].
Natural factors such as the low geomorphic potential for runoff, typical of a karstic plain, combined with the effects of engineering works built along the plain, have led to a reduction in the amounts of nutrients reaching the sea from the land. Hence, this brings about a decrease in the natural rate of biogeochemical cycles.) These results agree with those of Baisre [4,7,8] and Claro et al. [9] who find a decrease in total biomass (Fig. 2) and fish length in the gulf since the early 1990s. This declining trend coincides with the completion of the Dique Sur (South Dike) between 1985 and 1991 [10] and the completion of the last dam in 1991 [11].
Studies carried out in Cuba on soil management to date do not contemplate its influence on the marine ecosystem. These studies focus on soil management for its fertility, either via compost, or legumes, among other practices such as those carried out by [12][13][14][15] and others. In addition, research methods for comparing agricultural systems (conventional and organic) are established as well as highlighting which ones have better socio-economic benefits and less soil degradation [16][17][18]. From the marine science side, there are no studies that delve into this topic, focusing only on the coastal zone at the mangrove-marine or beach-marine water interface [19][20][21]. This apparent disconnection is perhaps due to the fact that the waters of the Cuban gulfs do not present great degrees of eutrophication and the volumes of nutrients that reach the waters due to agriculture are low, but they leave aside that the problem of oligotrophication of the gulfs can be alleviated or maintained in a stable situation avoiding its depletion owed to agricultural changes.  1960 1963 1966 1969 1972 1975 1978 1981 1984 1987 1990 1993 1996 1999 Inorganic N ferƟlizers (Gt) In the last decade, there have been no studies that manage soil and oligotrophic marine water resources and their interrelationships in the territory. This type of approach is the basis for knowing the main variables that determine the system and the physical characteristics of the study area or/and type of soil management. In this research, soil management focused on the potential rates of nitrogen fixation, the surface extent of this practice in the basin, and the drift of nitrogen to the sea in each of the different combinations of the two elements mentioned above. Knowledge of these elements allows estimating the cost-benefit balance of the soil resource under legume management. All of the above will contribute to better management, planning, and creation of better policies to manage Cuba's natural resources.

Legume choice, Phaseolus vulgaris, and Leucaena leucocephala
The choice of crops was made based on the existing specialized literature, social and cultural factors, and experience with the plant. Bean (Phaseolus vulgaris) was chosen in the areas of temporary crop use, since the insertion would not present alteration in eating habits, in addition to complementing the culinary cultural heritage of this crop. Leucaena leucocephala for pasture areas because of the large number of research studies carried out on the country. The beans present in a countless number of species, with great nutritional values of 20.2-23.4 g per 100 g of protein [22]. According to the study by González Robaina et al. [23], beans in Cuba is a crop that requires little irrigation demand (280-380 mm for 1.6-2.98 t ha −1 ), considering this a rainfed crop. Yields greater than 100% can be expected from this crop when the necessary water is applied to satisfy its water demand. Several times this demand can be solved with the rains, as can be seen in the results of Flores et al. [24] in Zacatecas, Mexico. Another study was conducted by Torres and Mendoza [25] in Nicaragua, where six different treatments were used for the bean crop, fertilization (non-application, chemical, and organic), weeds (control and non-control) and was subjected to the comparison of two tillage systems: conventional and conservation tillage. It was concluded that the best yields were found in the fertilizer application treatment with weed control and traditional tillage, but not the one with the best cost-benefit, being the control treatment with weed control.
The use of leguminous plants is a common practice in livestock farming, either to feed them directly or as fertilizer due to their capacity for nitrogen fixation or as the stubble that is incorporated into the soil. Based on research conducted in the last decade on the island, it is proposed to use L. leucocephala cv. Cunningham as a legume for fertilizing pasture fields in a silvopastoral system [29,30].
According to the literature, their fixation rate ranges from 90 to 170 kg N ha −1 yr −1 [41], 100 to 500 kg N ha −1 yr −1 [42], and 300 kg N ha −1 yr −1 [43]. The average rate of transfer of this nitrogen to the associated pasture in each of these studies was 30%, reaching 50% in some cases, as in the first study referenced. In this context, we used the methodology and principles used in the research carried out by Febles et al. [44]. This reveals and highlights the importance of applying soil erosion models of the areal domain in karst environments of intermittent surface streams; which, although they do not offer an exact simulation of real situations, allow obtaining approximate results, according to the assumed simplifications and applying them at plot or basin scale. This article aims to manage the dynamics of nutrient fluxes under leguminous management in order to obtain the maximum benefit for soil and marine resources.
The present study has been divided as follows. Section 2.1 presents the physical-geographical characteristics of the study area. The following subsection 2.2 shows the models used for the simulation and provides the reasons for the choice of legumes. Section 3.0 presents the results achieved, which are subdivided into Sect.

Geographical location
The research was carried out in the La Teresa Karst basin located on the southern coast of Cuba (Fig. 3); belonging to the Southern Havana-Matanzas Karst plain (Table 2). Limited to the north by the coordinates 82° 7′ 0″ W and 22° 55′ 0″ N, to the east 82° 15′ 0″ W and 22° 47′ 0″ N; to the west 82° 0′ 0″ W and 22° 55′ 0″ N and to the south by the coastline.
The basin was mapped following the watershed lines, according to Strahler's methodology [45]. It also has a surface hydrography with intermittent currents.

Climatic regime
According to the Köppen-Geiger Climate Classification, 1991-2020 [46], this can be considered tropical with alternating humidity, with a long rainy period and a less humid one. In general, the rainy season extends from May to October and the less humid period extends from November to April (Fig. 4).
It is estimated that the amount of precipitation that falls annually in the province of Mayabeque ranges from 1200 to 1600 mm [47] (Fig. 4), and the average temperature ranges between 24 and 26 °C.

Geological-geomorphological characteristics
The geomorphology corresponds to a cumulative plain with a flat terrace system (Fig. 3) consisting of carbonate rocks and marsh deposits [48].

Soils and land use
The predominant soils correspond to the Red Ferralitic Grouping [49,50], basically used for various crops and pastures (Fig. 5).
Red Ferralitic soils, despite having high clay content, have low plasticity, high capillary elevation, and low hygroscopicity, which favors their use for a wide range of plant species [52]. In terms of soil use, sugarcane is the only crop that continues to be cultivated under traditional agricultural management.

Nitrogen fixation rate
As discussed above (Sect. 1.1), nitrogen fixation values of 30, 40, 50, and 125 kg N ha −1 yr −1 for crops and 90, 100, 300, and 500 kg N ha −1 yr −1 for pasture were taken, covering the low, intermediate, and high scenarios. These practices are used in 10, 20, 30, and 40% of the area employed in these uses, with 40% being the highest range.

Soil nutrient content
The concentration of nutrients ( Fig. 6), concerning N in the soil of the basin, is medium (0.14-0.25%) [53,54] with higher values towards the center of the territory.
The P concentration in the soils is moderate to very high (10-17 moderate, 17-25 high, greater than 25 ppm very high) [54,58], with higher values in the southern part. Finally, K concentration in the soils is low (less than 44 ppm) [54,58,59].

Soil erosion
Mean annual erosion, according to Febles et al. [44], is 0.6 ± 0.1 (t ha −1 ) (Fig. 7), being classified as very low according to Geler [60]. The transport capacity of surface flow (G) was the limiting factor of erosion.
The low rates are the result of a gentle slope and the mechanical composition of the predominantly clayey soils (85%), which have a high degree of resistance to detachment by raindrops. The highest erosion rates were estimated in the cultivated sectors, both in short and prolonged periods.

Nutrient flux and export of nutrients to the sea
According to Febles et al., 2020, the highest export rates in the basin were recorded in those sectors occupied by various crops on clay-textured soils (Fig. 5A, B). These sectors are also the most eroded ( Fig. 7) and with the greatest loss of nutrients. Sub-basins No. 6 and 10 are those with the highest values (Fig. 8).
The zonal-spatial pattern in nutrient export is higher in the middle sector of the basin and lower in the other two sectors, associated with the effectiveness of sediment traps in the upper sector and the low geomorphic potential near the coastline (Fig. 8).
Nutrient export of nitrogen and phosphorus was 35 ± 1 t (1.0 ± 0.02 kg ha −1 yr −1 ); 2.3 ± 0.4 t (0.064 ± 0.01 kg ha −1 yr −1 ) respectively. These are classified as low when compared to other studies, none of these developed in karst regions.   These low values are a consequence of the low erosion rates of the Red Ferralitic and Alfisol soils since they present moderate to high values ( Fig. 6) of nutrients [44].

Methodology used
The management principles proposed by Febles et al. [44]  were applied to determine the optimal management of soils in the La Teresa Basin under legume management, which exert the least possible pressure on marine ecosystems. These indicate maintaining the speed of biogeochemical cycles and increasing the concentration of nutrients in soils. This study focused only on one chemical element of fundamental importance for agriculture, nitrogen, through the organic management of legumes. Ten thousand possible scenarios were simulated of nitrogen fixation by legumes in different surface extensions and the export of nutrients to the Gulf of Batabanó in different spatial configurations taken at random. The selected areas of use were those currently occupied by temporary crops and pastures (Fig. 4B) since they allow flexible management as they are not strategic crops such as sugar cane or permanent crops. The simulation was performed in a Python console, under the methodology of Febles et al. [44] and Febles et al. [61] the first for estimating nutrient fluxes and the second for estimating soil particle fluxes (Fig. 9).

Calculation of nutrient fluxes
To calculate the amount of soil eroded (Fig. 9), the Morgan, Morgan and Finney (MMF) parametric erosion model [62] was applied, based on the dynamics of erosion processes, soil properties, climatic characteristics, and conditions of use, which made it possible to evaluate the magnitude of losses.
The equation of Verstraeten and Poesen [63] was used to calculate nutrient export. This equation was designed for use in small basins to quantify sediment-bound nutrient transfer based on measurements of sediment volumes and nutrient content. It is most effective when erosion rates are medium to low, such as those characteristics of a karst plain.
MMF model: (1) E = R 11.9 + 8.7 ⋅ log 10 I The final prediction of the soil loss is made by comparing the values obtained of G and F. The lowest of them is taken as annual soil loss rate.
For the estimation of nutrient losses, the assumption was that the basin discharges by a single sector into the sea. Equation (7) with slight modifications (8) was used to estimate the effects of the amount of nutrients and sediments coming from the upper part of the basin through 13 small sub-basins that were mapped (Fig. 3).

Pastures
Land used for pasture occupies 21% of the basin area, mainly located in the upper portion of the basin. The average soil nitrogen fixation concentration in the different treatments was within the range of 1.07-1070 ppm with a standard deviation of ± 2.23-141 ppm, with a maximum fixation of 1619 ppm (Fig. 10).
Significant differences were found in all treatments in the Conover-Iman test (Fig. 10). The export to the Gulf of Batabanó ranged between 0.9623 and 0.9624 kg ha −1 yr −1 , with the maximum value coinciding with the treatment of 500 kg ha −1 at 40% of the area under legume management and less than 90 kg ha −1 at 10% (Figs. 11, 12).

Crops
The land used for cultivation occupies 27% of the surface of the basin, mainly located in the middle portion of the basin. The average soil nitrogen fixation concentration in the different treatments was within the range of 0.3-1017 ppm with a standard deviation of ± 0.85-69 ppm, with a maximum fixation of 1301 ppm (Fig. 13).
Significant differences were found in all treatments in the Conover-Iman test (Fig. 13). Export to the gulf ranged between 0.962379 and 0.962474 kg ha −1 yr −1 ,  (Figs. 14, 15). The dynamics of nutrient export when combining the spatial configurations of pasture and temporary crops in which higher export is achieved is shown in Fig. 16. The average difference between using 40% and 100% of the area is 0.0003 kg ha −1 yr −1 of N under the maximum capture rate reported in the literature.

Costs and revenues
Land, in terms of agriculture, turns out to be the most significant single asset. Soil degradation is a relevant phenomenon in agricultural activity; therefore, it is necessary to take it into account. The negative effects of degraded soil on the economy are more severe in countries that strongly depend on agriculture [65,66] and there is no doubt that the agricultural sector is one of the most significant in the Cuban economic activity [67].
For a correct estimate of its devaluation, it is necessary to consider the number of ways in which the phenomenon of degradation can be found. Bearing in mind that, it is common for more than one of these degradation variants to subsist in the same geographical area, these can be: excessive extraction of nutrients, saturation of salts, contamination with chemical products, soil erosion, vegetation degradation, logging of forests for farmland, water use, etc.. This process causes a decrease in the productive land capacity, decreasing potential yields. This decrease must be reflected by accounting in the financial statements [68].
Due to the complexity of estimating the depreciation of the soil resource, this research only offers an approximation of the costs produced by the export of nutrients caused by water erosion of the soils of the basin and inputs produced by the use of legume management. In estimating these costs, the value of fertilizers in the market was used (Table 3) and the spatial configurations that produced the highest export value (Figs. 12, 15).
The highest nutrient loss costs were found in sub-basin number 6 (Supplemental material III) between the range of 44-60 USD ha −1 yr −1 in pastures and 44-62 USD ha −1 yr −1 in crops. The income rates in each of the treatments, both in crops and pastures, vary according to the spatial configuration. These incomes have zero values in the sub-basins where these land uses are not present, mainly those found in the lower zone of La Teresa (Figs. 3, 5B).  The general annual capital flow balance shows that the pasture treatments are more profitable than the crop treatments (Table 4), only taking into account the input and output of N, leaving out of the analysis the value of beans in the national and international market. The pasture treatments had a positive balance of 22 USD ha −1 yr −1 , which was less than the crop treatment of 55 USD ha −1 yr −1 ; but the latter, used more area at higher rates. To see these costs in situ would be to go against the systemic model of nature whereby everything that is lost in one, place ends up redeposited in another. In this case, what is lost by erosion in the middle and upper basin fertilizes the lower basin; and what escapes from the basin, in general, reaches the sea where it is used by marine biota as an external nutrient input. Therefore, proper accounting becomes complicated when nature only understands cycles and not losses. But these approaches serve to address the problem from another reference system, much more appropriate when estimating the economic benefits of implementing ecological measures [70].

Nitrogen fixation in soils of the La Teresa basin
A distinctive feature of legumes is their ability to develop root nodules and fix N 2 in symbiosis with compatible Rhizobium, which brings about increased nutrient properties of soils [71]. The data show that all legume treatments are effective in increasing soil fertility in the basin, being higher as the rate of legume fixation increases.
However, these rates do not progress linearly over time. Specialized research shows that as more nitrogen becomes available in soils, these rates begin to decrease, since nutrient uptake directly from the soil is thermodynamically more feasible than from the symbiont [72−74]. There is evidence that legumes use more photosynthates for N assimilation if the N comes from N 2 fixation compared to soil N uptake [75−79]. Increasing any form of N in the soil or nutrient solution reduces root nodule numbers and N 2 fixation rates [80,81].
Nitrogen-fixing bacteria are equipped with an enzyme complex called nitrogenase. Nitrogenase is a protein enzyme, composed of two metalloproteins: ferroprotein and molybdoferroprotein. This enzyme catalyzes the primary nitrogen fixation reaction, in other words, the transformation of nitrogen present in the atmosphere into ammonia. The nitrogenase enzyme complex is irreversibly inactivated in the presence of oxygen [82].
The supply of oxygen (O 2 ) to the interior of the nodule is a key factor in the regulation of nodule activity [83]. This regulation can be either by excess, which could jeopardize the destruction of nitrogenase, or by defect, preventing the production of ATP necessary for fixation [84]. This problem presented by legumes is not a major limitation if   These optimum states may vary from sector to sector, due to the characteristics of management and crops that have been used, with density being the main indicator of loss of pore space.
The soils of the basin currently do not present any P limitation, but a continued use of this management under high N fixation rates can cause a decrease in this element, causing imbalances in long-term biogeochemical cycles. A study by Graham and Vance [85], showed that P deficiency constitutes a limitation for N 2 fixation in legumes, particularly in soils with low P availability, which has been associated with slow plant growth and a decrease in N 2 fixation [86,87]. In another study by Robson et al. [88] in soils with low P availability, legumes dependent on N 2 fixation, have responded positively to P fertilization and have shown a higher N content in shoots and roots.

Nitrogen export
The amount of nutrients that reach the sea in each of the different treatments in the ten thousand simulations, do not present high values. The difference between the use of 40% relative to 100% of the surface destined to pasture is 0.004% and 0.005% for crops. In the union of both treatments, the values rise quite conservatively by 0.28% for crop + pasture at 40% with respect to pasture or crops separately, and 0.29% for 100% of the area under such management (Fig. 16).
The non-increase in nutrient export values concerning the increase in nutrient concentration in soils due to N 2 fixation produced by legumes, is mainly because of the physical characteristics of the basin. The convex geomorphological layout influences the overall behavior of sediment transport, with sedimentation prevailing in the more depressed sections. The most notorious inflection point of this convexity is the flattening of the slope in the lower sector, which can be appreciated from 19 km away in Fig. 17, corresponding to point C on the map in Fig. 3.
The dynamics governing this particular spatial configuration is the fact that the horizontal component of gravity decreases along the less steep slopes, which decreases the runoff velocity and thus lengthens the runoff generation time [89,90]. Therefore, the sub-basins of the lower part of La Teresa act as a geomorphological barrier for sediments on their way to the sea. Figure 17 also shows the configuration exerted by the karstic environment of the basin in a series of "U" and "V" shaped karstic depressions that act as geomorphological barriers arranged transversally to the slopes, intercepting runoff, as well as soil erosion that transports nutrients [44,91]. The mechanical composition of the soils is predominantly clayey (85%), which has a high degree of resistance to detachment by raindrops [92,93]. Therefore, texture is another factor to take into account when analyzing the arrival of nutrients to the sea. Febles et al. [44] 500 kg ha-1 125 kg ha-1 125 kg ha-1+500 kg ha-1 kg N ha -1 año -1 40% surface 100% surface Another variable to consider that influences the dynamics of nutrient export from La Teresa basin is the spatial arrangement of the areas of use, with exports being greater as the distance to the coastline shortens. The results show that the difference across spatial configurations and between the two land uses, which achieve the highest export values, only differs in the sixth significant digit after zero (Fig. 16). Both uses cover similar areas (21% pasture and 27% temporary crops (Fig. 5B) and markedly different nutrient fixation rates 90-500 kg ha −1 yr −1 for pasture and 30-125 kg ha −1 yr −1 for crops. However, since the cropland area is mainly located in the middle sector of the basin, as opposed to the pasture area in the upper sector, it manages to maintain a similar level of exports since it requires less total geomorphic potential to reach the coastline (Fig. 5B).

Soil use management in the basin
Given the above, the basin management proposed in accordance with the topic of soil management with leguminous plants so that it does not exert more pressure on the oligotrophic ecosystems of the Gulf of Batabanó is as follows: This research does not require the use of this agroecological practice to increase yields in cultivated areas; however, it strongly recommends its use in these regions. In the case of applying this management, it should be done under the economic risks of having less land available for the cultivation of other crops during one stage of the year. Being a recommended practice, it gives autonomy to farmers to design/redesign the strategy that best suits their harvest and crop calendars, being legume management one of the many strategies that can be followed by this sector.
This research does encourage legume management in pasture areas, if feasible to the greatest extent possible. First, this management is suggested because the use of Leucaena leucocephala does not require such demanding care as would be required by a legume whose destination is the agricultural market or self-consumption, in addition to the fact that its life cycle is longer. This management also showed a positive economic balance, under low fixation rates depending on its extension or moderate to high N 2 fixation rates regardless of its extension (Table 4). Besides, the use of this management in these regions and its location in the upper zone of La Teresa, provides an indirect fertilization of the lower zones of these regions where the area of temporary or permanent crops is mostly concentrated (Fig. 5B). This would not only benefit the livestock sector, but also the agricultural sector.
The forested areas near the coastline (Fig. 5B) should not be disturbed in order to increase the export of nutrients to the shelf. This is a natural barrier against extreme weather events by attenuating storm surges and windstorms [94]. This arboreal strip also protects from marine aerosol, which, according to Morcillo et al., [95] can have concentrations of 100 mg Cl m −2 at 1 km and 9 mg Cl m −2 at 5 km. Although, these data are subject to wind speed, obstacles, topography configuration and the concentration of sodium chloride in the seawater, an adverse effect on soils and infrastructures due to corrosion. Another advantage of this forest buffer is that they slow down or increase the shoreline as they retain sediment and stabilize the soil in the areas between high and low tide (intertidal zones) reducing erosion during storms and floods [96].
This study does not recommend a rate of N 2 fixation by the legume, because even the highest rates reported in agriculture do not achieve a notable export of nutrients. These fixation rates are not continuous over time for physiological reasons of the symbionts, and may vary from year to year and within the same year. This legume management is recommended to use as widely as possible in the basin.

Representativeness of this result to the rest of the southern zone of Cuba
Cuba is an island with a central watershed and a particular tilting where the northern zone rises, and the southern  zone submerges [97]. The geomorphology of this southern zone is mainly terraced (Fig. 18), with a slight convexity (Fig. 19), flattening several kilometers before reaching the coastline. These topographic characteristics, being exposed to the same type of climate and a mainly karstic lithology, produce similar soil groups in these regions (Fig. 20). Therefore, it can be observed that around the Ana María and Batabanó Gulfs, there is a predominance of Ferralitic soils. In the upper part, near the watershed, there are Brown soils or Fersialitic soils that have similar textural characteristics (silty).
The only area of these two gulfs that stands out for being quite dissimilar in this aspect, is the quartzite plain of Pinar del Río. Its texture is primarily sandy and possibly low in nutrients due to the little electric interaction that exists between its aggregates, in addition to the fact that an erosion much greater may occur than in the other sectors.
Therefore, it may not be surprising to extrapolate several of the results obtained in the La Teresa basin to almost Mainly those of physical-geographical origin such as geomorphic potential, slope shape, dissection caused by karstic relief, and similarity in soil cover at the group level. Another aspect in which both areas are similar is that their shorelines are generally covered by a transition of mangrove forest with vegetation more in line with automorphic surface terrain. The exception to this rule is Dique Sur (South Dike) where the mangrove forest has deteriorated.

Conclusions
The soil management under legume management that maximizes gains and minimizes losses for both agricultural and marine ecosystems in the basin is in the areas used for pasture since the legume used does not require as demanding care as a legume for human consumption. This management also showed a positive economic balance at low fixation rates (from 55 to 1560 USD ha −1 yr −1 for crops and from 22 to 2400 USD ha −1 yr −1 for pasture). Its location in the highlands of La Teresa provides indirect fertilization of the lowlands. In agricultural regions, the use of this management is left to the personal discretion of the farmers because of the benefits in productivity and field yields. The maximum nutrient concentration in La Teresa by legume management was 1.07-1070 ± 2-141 ppm in pastures and 0.3-1017 ± 0.85-69 ppm in crops. The amount of nutrients that reached the Gulf of Batabanó in each of the different treatments did not show high values (0.9623-0.9624 kg ha −1 yr −1 for pastures; in the case of crops 0.962379-0.962474 kg ha −1 yr −1 ). Regarding the spatial arrangement of the areas of use, it follows from the results that the shorter the distance to the coast, the higher the export. The main variable regulating this export is the geomorphic potential of a convex basin, with sedimentation predominating in the lower locations. These results demonstrate that through legume management it is not possible to increase the export of the nutrient nitrogen to the oligotrophic waters of the Gulf of Batabanó without causing major modifications to the ecosystem. Therefore, the future scope of research should be the management of this fishing zone, adequate to these new conditions of low nutrient export values from land and damming.
Funding The authors did not receive support from any organization for the submitted work.

Conflict of interest
On behalf of all authors, the corresponding author states that there is no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.