Flow and Transport Numerical Model of a Coastal Aquifer Based on the Hydraulic Importance of a Dyke and Its Impact on Water Quality. Manglaralto—Ecuador

: Coastal aquifers are part of the natural resources contributing to local development and promote resilience in the most vulnerable communities near the sea. Manglaralto, an Ecuadorian coastal parish, is affected by water resource scarcity. The increase in salinity and deterioration of the water quality is generated by the local and ﬂoating population’s demand, causing an increase in the Total Dissolved Solids (TDS) concentrations and decreasing the aquifer’s piezometric levels. The aim is to establish a numerical model of ﬂow and transport of the Manglaralto coastal aquifer by using hydrogeological data and Visual Transin software, relating the hydraulic importance of a dyke’s design (“tape”) and its impact on the quality of the water. The methodology is (i) hydrogeological database analysis, (ii) the system’s recharge concerning the soil water balance, (iii) the boundary conditions of the ﬂow and transport model and, (iv) the results and validation of the numerical simulation. The results conﬁgure the importance of the coastal aquifer’s artiﬁcial recharge in the area where the tape is located, as reﬂected in the increase in piezometric levels and the decrease in salinity in wells near the sea. In conclusion, the numerical model of ﬂow and transport allows expanding the knowledge of the variation of the piezometric levels and TDS concentrations over time, the importance of recharge in the hydrogeological system’s operation, and correct community management resilience and projection to sustainable development.


Introduction
An aquifer is a natural reserve of freshwater corresponding to 1.7% of the Earth's total and about 30.1% of Earth's freshwater except for the Poles [1], a natural space that must be protected and monitored by several organizations. The wetlands convention (RAMSAR), one of the first global and environmental conventions adopted since 1971, its objective is the conservation of surface and groundwater hydrological systems for the protection of coastal aquifers due to seawater's advance is a complex process that can be analyzed using numerical models, such as the case of Querença-Silves in Portugal and Andhra Pradesh in India, where the hydraulic system's behavior and the correct management of the aquifer are highlighted [31,32].
The numerical models of flow and transport are used within scientific research to establish the hydraulic behavior in advance of seawater intrusion, validate the hypotheses raised in hydrogeological conceptual models, as well as the evolution in time of piezometric levels and concentrations of groundwater, for the proper management and planning of the water resources (in this case, coastal aquifers) necessary for the development of communities [33][34][35].
Santa Elena is a coastal province of Ecuador, recognized for its geological characteristics, mineral resources, historical geoscientific knowledge, and mainly for its tourist contribution to local development [36]. It has problems in distributing water to meet its communes' demand [37]. It has two primary sources of water; in the southern area, which corresponds to the urban part, the water resources are obtained through the Daule-Santa Elena diversion from the Guayas river basin, while, in the northern region, where the rural sectors are, it is supplied mainly from coastal aquifers [38].
In the Manglaralto parish, the primary source of water is underground, extracted from the coastal aquifer associated with the river of the same name [39], and in turn, it is considered a geosite of great relevance in the Geopark Project of the Santa Elena Peninsula [36]. The daily recharge occurs through the infiltration of surface water (rivers and streams that cross the limits of the aquifer), in addition to high-intensity rainfall in winter (December to May frequently) and small drizzles in June, July, and August. One of the difficulties present is the material that composes it, which is very porous and permeable, which causes the fast flow of groundwater to the sea. The Regional Administrative Board of Drinking Water of Manglaralto (JAAPMAN) is in charge of supplying the vital liquid to six communes (Figure 1), through 12 wells drilled. Although it is possible to supply most of the inhabitants, the rest of the population uses clandestine wells, "albarradas" (artisanal dams), tankers, and other sources [18]. For this reason, the scope of water in Manglaralto is limited and does not cover the necessary and economic activities of the area [39], identifying the sustainable management of the aquifer as a problem [40].
Water 2021, 13, x FOR PEER REVIEW 3 of 23 processes that alter the physical properties of an aquifer [29,30]. Pollution in coastal aquifers due to seawater's advance is a complex process that can be analyzed using numerical models, such as the case of Querença-Silves in Portugal and Andhra Pradesh in India, where the hydraulic system's behavior and the correct management of the aquifer are highlighted [31,32]. The numerical models of flow and transport are used within scientific research to establish the hydraulic behavior in advance of seawater intrusion, validate the hypotheses raised in hydrogeological conceptual models, as well as the evolution in time of piezometric levels and concentrations of groundwater, for the proper management and planning of the water resources (in this case, coastal aquifers) necessary for the development of communities [33][34][35].
Santa Elena is a coastal province of Ecuador, recognized for its geological characteristics, mineral resources, historical geoscientific knowledge, and mainly for its tourist contribution to local development [36]. It has problems in distributing water to meet its communes' demand [37]. It has two primary sources of water; in the southern area, which corresponds to the urban part, the water resources are obtained through the Daule-Santa Elena diversion from the Guayas river basin, while, in the northern region, where the rural sectors are, it is supplied mainly from coastal aquifers [38].
In the Manglaralto parish, the primary source of water is underground, extracted from the coastal aquifer associated with the river of the same name [39], and in turn, it is considered a geosite of great relevance in the Geopark Project of the Santa Elena Peninsula [36]. The daily recharge occurs through the infiltration of surface water (rivers and streams that cross the limits of the aquifer), in addition to high-intensity rainfall in winter (December to May frequently) and small drizzles in June, July, and August. One of the difficulties present is the material that composes it, which is very porous and permeable, which causes the fast flow of groundwater to the sea. The Regional Administrative Board of Drinking Water of Manglaralto (JAAPMAN) is in charge of supplying the vital liquid to six communes (Figure 1), through 12 wells drilled. Although it is possible to supply most of the inhabitants, the rest of the population uses clandestine wells, "albarradas" (artisanal dams), tankers, and other sources [18]. For this reason, the scope of water in Manglaralto is limited and does not cover the necessary and economic activities of the area [39], identifying the sustainable management of the aquifer as a problem [40].  Manglaralto has come to be considered a natural laboratory, given the community management for 40 years in the development, use, and protection of the aquifer and academic and scientific participation, through various research projects and multiple undertakings about this resource. JAAPMAN is responsible for using groundwater and considers the management of wastewater with stabilization ponds and green filters for reforestation, using the soil as a natural filter in the commune of Montañita [41]. On the other hand, the issue of water and its management has recently been treated as an object of greater interest by many government researchers thanks to community participation, accompanied by ancestral knowledge, which have been key factors in solving problems regarding obtaining and distributing water, ensuring that the residents develop techniques to recharge the aquifer [18]. Within this participation, there are other significant parties involved, such as the International Atomic Energy Organization (IAEA) and the ESPOL Polytechnic University, which together with the community have managed to meet the growing water demand, rescuing ancestral knowledge by implementing tapes (artisanal dykes), addressing the overexploitation of the aquifer and empowering communities in water management [42]. However, in the geographical environment, problems remain due to tourism growth and climate change.
The main factors that affect and increase the water problem in the study area are marine currents and weather phenomena such as the warm El Niño current and the cold Humboldt current, which affect the intensity and regularity of rainfall [43,44]. In addition to having a semiarid climate [45], increased seawater intrusion is notoriously damaging water quality [38]. One of the solutions implemented is the rescue of ancestral knowledge with a technical-artisan tape [8] completed at the end of 2018. "Tape" is an ancestral-cultural term that the community members assign to the accumulation of rocks and sediments in certain strategic parts of the riverbed, damming the water or accumulate with the consequent artificial recharge aquifer [18]. By December 2019, it was evident how the coastal aquifer affected the communities' socioeconomic development. The JAAPMAN provided the water service 24 h a day without interruption, positively influencing tourism in the sector, unlike in the past, where the service was only available at certain times.
As of March 2020, due to the global pandemic, another inevitable challenge is presented for Manglaralto, focused on overcoming the difficulties of proper water management. Studies carried out by [46] indicate that the groundwater system in Manglaralto is resilient where the community, despite global circumstances, has managed to adequately manage the water to avoid overexploitation of the aquifer without ruling out the household water. The aquifer is essential for these communities' resilience in the face of the global crisis, even though it is an area affected by water scarcity related to climate change [18]. It has been possible to provide adequate hygiene, which is fundamental for each of the inhabitant's very life and the resurgence of the coastal area, due to control and dedication by JAAPMAN towards its clients, despite the decrease in its leading economic resource, tourism [47].
The high demand problem due to socioeconomic relationships and climate change incidence by El Niño and La Niña leads us to ask ourselves the following question: is it possible to determine the advance of seawater intrusion by the influence of a tape on hydraulics behavior and groundwater density through a numerical flow and transport model?
This research aims to establish a numerical model of flow and transport of the Manglaralto coastal aquifer, based on the conceptual model of the hydrogeological system, through the use of hydrogeological data and the application of Visual Transin software, relating the hydraulic importance of the design of a tape and its impact on water quality. This analysis would help contribute to the resilience and sustainable development of the Manglaralto coastal aquifer. It will allow us to predict how the tape works and its importance in artificial recharge, positively affecting the biodiversity and geotourism of this Ecuadorian coastal zone.

Materials and Methods
The hydrogeological conceptual model encompasses all the specific data and scientific knowledge of the study site ( Figure 2), managing to discretize the data domain and behavior of the processes that influence water flow and the transport of solutes [48]. The water table in the middle and upper zones follows the geological survey guidelines carried out in the Manglaralto riverbed, indicating that in the middle zone, the average power of the aquifer is 10 m. In comparison, in the upper zone, the average power is 5 m. It is assumed that the water table is from 2 m below the level of the riverbed.

Materials and Methods
The hydrogeological conceptual model encompasses all the specific data and scientific knowledge of the study site ( Figure 2), managing to discretize the data domain and behavior of the processes that influence water flow and the transport of solutes [48]. The water table in the middle and upper zones follows the geological survey guidelines carried out in the Manglaralto riverbed, indicating that in the middle zone, the average power of the aquifer is 10 m. In comparison, in the upper zone, the average power is 5 m. It is assumed that the water table is from 2 m below the level of the riverbed. Before the methodological process that this investigative work entails, the vital work of the scientific information carried out in the study area must be considered. These data are stored using various works carried out by the community and academic institutions, such as geoelectric prospecting, workshops with the community, and recording the physical-chemical parameters and the static-dynamic levels of groundwater wells. These studies have a link with various projects (Table 1), such as:  Before the methodological process that this investigative work entails, the vital work of the scientific information carried out in the study area must be considered. These data are stored using various works carried out by the community and academic institutions, such as geoelectric prospecting, workshops with the community, and recording the physicalchemical parameters and the static-dynamic levels of groundwater wells. These studies have a link with various projects (Table 1), such as: Darcy's law governs the movement of underground flow in porous media [49]. The expressed conditions of variable density are in the supplementary material section.
The Visual Transin 4.0 and Transdens software solve the flow equation for a saturated and unsaturated zone, steady and transient, for transport, including a scattering-diffusion and advective process. The flow and transport equations can be developed decoupled (in constant density) or coupled (in variable density). These codes solve the flow and transport equations using finite elements in space and finite differences weighted in time.
This research work focuses on analyzing the scientific information and data collected from the study area, bearing in mind the zone's geological and hydrogeological conditions, showing the reality of the coastal aquifer's physical environment in a numerical model of flow and transport. A logical sequence is stipulated that starts from the conceptual model of the aquifer be studied. Subsequently, with the entire hydrogeological database and the zone's hydrological cycle data, a soil water balance is carried out. It can be observed in the following methodological scheme, in the Figure 3.

Phase I
The compilation of information from the 11 study wells helps to ensure good order and technical formatting of the hydrogeological information (hydraulic, physical-chemical, and geological) managing to visualize in a general way the behavior of each well according to the month and year, and thus its direct relationship with piezometric levels,

Phase I
The compilation of information from the 11 study wells helps to ensure good order and technical formatting of the hydrogeological information (hydraulic, physical-chemical, and geological) managing to visualize in a general way the behavior of each well according to the month and year, and thus its direct relationship with piezometric levels, daily extraction, precipitations, and temperatures in the sector. These wells provide integrated hydrogeological information of the coastal aquifer, correlated with geology and geophysics, as shown in the conceptual model ( Figure 2).
During this phase, the database of each well is analyzed. This information has been registered and reported monthly by JAAPMAN and the Center for Research and Projects Applied to Earth Sciences (CIPAT-ESPOL, acronym in Spanish), such as Total Dissolved Solids (TDS), salinity, temperature, resistivity, piezometric levels, and daily pumping extractions. The information on precipitation and temperature in the basin is from the Manglaralto meteorological station (Code: M0619).

Phase II
The Manglaralto coastal aquifer covers approximately 508 hectares and corresponds to 3.84% of the Manglaralto river basin in the surface territory. Permeable and porous materials predominate in this sector, identified as an alluvial terrace [50]. Figure 4 shows the geological formations that directly influence the alluvial aquifer, consisting mainly of gravel, sand, silt, and clay materials [38]. canic rocks a clay matrix, and gray or green silicified argillites, unlike the Tablazo Formation, where sequences of conglomerates and sands are evident. The Socorro Formation (Ancon Group) is stratigraphically like flysch with turbidites and thin layers of shales and greenish sandstones [51]. The Zapotal Member that belongs to the Tosagua Formation contains significant amounts of chloride due to continuous contact with seawater and contains sandstone, conglomerates, and shales. The transmissivity of the hydrogeological system is determined considering a theoretical permeability of 40 m/day (m/d), according to the gravel, sand, silt and clay materials present in the alluvial terrace, with an average saturated thickness of 20 m, 10 m for the middle zone and 5 m for the high zone.
In the soil-water balance, characteristic parameters such as field capacity, wilting point, precipitation, and temperature entered into the Easy Balance software were considered, estimating the aquifer recharge and Potential Evapotranspiration (PET). Recharge is a fundamental parameter in the simulation process, given its importance in the behavior of piezometric levels and water quality in these free and coastal aquifers.
It was essential for the field capacity and the wilting point to consider the clay soil according to its loamy texture, represented in volumetric content from percentage data and soil density [52,53]. Due to the presence of weeds in Manglaralto, it was possible to The formations are sequences of sedimentary and volcano-sedimentary rocks, and the ages range from the Upper Cretaceous to the Lower Miocene, followed by alluvial deposits from the Holocene [38].
This basement belongs to the Cayo Formation, which has powerful sequences of volcanic-sedimentary rocks. It is recognized for presenting graywacke sandstones and volcanic rocks a clay matrix, and gray or green silicified argillites, unlike the Tablazo Formation, where sequences of conglomerates and sands are evident. The Socorro Formation (Ancon Group) is stratigraphically like flysch with turbidites and thin layers of shales and greenish sandstones [51]. The Zapotal Member that belongs to the Tosagua Formation contains significant amounts of chloride due to continuous contact with seawater and contains sandstone, conglomerates, and shales.
The transmissivity of the hydrogeological system is determined considering a theoretical permeability of 40 m/day (m/d), according to the gravel, sand, silt and clay materials present in the alluvial terrace, with an average saturated thickness of 20 m, 10 m for the middle zone and 5 m for the high zone.
In the soil-water balance, characteristic parameters such as field capacity, wilting point, precipitation, and temperature entered into the Easy Balance software were considered, estimating the aquifer recharge and Potential Evapotranspiration (PET). Recharge is a fundamental parameter in the simulation process, given its importance in the behavior of piezometric levels and water quality in these free and coastal aquifers.
It was essential for the field capacity and the wilting point to consider the clay soil according to its loamy texture, represented in volumetric content from percentage data and soil density [52,53]. Due to the presence of weeds in Manglaralto, it was possible to estimate a root thickness of 0.25 m and an initial humidity of 0. The design hyetograph of the Manglaralto basin indicates that the laminator value is 10 mm [54,55] Finally, there is a useful reserve of 0 mm, and the initial reserve minor of 25 mm, used to calculate the PET and the recharge. Considering that the soil's humidity and the vegetal layer are optimal, the information recorded on rainfall and temperatures from May 2018 to September 2020 is used through the Thornthwaite method to calculate the PET [56].

Phase III
Considering the basin's hydrogeological system's conceptualization, it could represent it through two layers; the upper layer corresponds to the aquifer, while the lower layer represents the basement ( Figure 5). The first layer corresponds to the alluvial layer, classified into upper, middle, and lower, the latter strategically dividing into eight sub-zones according to the lateral recharge imposed on the model. classified into upper, middle, and lower, the latter strategically dividing into eight subzones according to the lateral recharge imposed on the model. The transmissivities, the storage coefficient, porosity, and dispersivity of the aquifer correspond to the first layer's three main zones. In the second layer (Figure 6), the aquifer's basement, in all its domain, presents a single zone of these parameters. In the basement, it is considered that there is no flow (Q = 0) because it is mainly made up of shale, silt and clay, for which a prescribed concentration condition is established.  The transmissivities, the storage coefficient, porosity, and dispersivity of the aquifer correspond to the first layer's three main zones. In the second layer (Figure 6), the aquifer's basement, in all its domain, presents a single zone of these parameters. In the basement, it is considered that there is no flow (Q = 0) because it is mainly made up of shale, silt and clay, for which a prescribed concentration condition is established.     The initial conditions used in the constant density model were the resolution of the first time step's flow and transport problem, since the model's numerical stability allowed it to go directly to the transitory model. This model was constructed prior to undertaking the variable density modelling.
The variable density model's initial conditions were the solution to the constant density flow and transport problem and by using the equivalent levels of fresh water in the constant and variable density models.
This work was carried out from 1 May 2018, to 30 September 2020; the simulation considered 29 months'-worth of data. It was possible to perform spatial discretization with a mesh of 952 nodes and 2068 elements. The largest mesh refining was carried out in the aquifer's lower area because it is an area of high hydraulic activity and groundwater discharge and is the sector where seawater enters.

Phase IV
In this section, the numerical model of flow and transport was obtained, that was subsequently analyzed to expand the knowledge of the aquifer piezometric levels and TDS concentrations, since this is related to the effect of the natural and artificial recharge produced by the water reservoir in the tape, its concentration in wells close to the sea due to marine intrusion and wells within the basin due to the saline basement's influence.

Results
The methodology used in this work provides the following results:

Soil Water Balance
Using Easy Balance 4.0 software, it was possible to obtain the soil-water balance; thus determining the PET and the recharge that enters the aquifer, using the simulation period's average monthly temperature.
March was the month with the highest PET (mm) in all years studied, with 200. 16, 197.55, and 193.70 mm, respectively, in 2018, 2019 and 2020, ( Table 4). The year with the highest value of this parameter is 2019 with 1596.95, while in 2020 until September, it is 1303.27 (Table 5); the PET in 2020 is higher compared to the years 2018 and 2019 in the same period. Before estimating the recharge, the soil's characteristic parameters in the study basin needed for its determination were determined ( Table 6). The recharge calculated for this period was 37 mm, which corresponds to 11.4% of the rainfall; in 2019, the recharge's highest value was 20.3 mm (Table 5). Table 4 shows the months that influenced the value of this recharge. It is essential to mention that in addition to the characteristics of the soil (Table 6), the PET and rainfall influence recharge, especially in the months of February-March, which is when the highest rainfall occurs, while in July-August, low-intensity precipitations known as "garuas" (in Spanish) occur, which infiltrate the alluvial terrace and recharge the aquifer. The data and results of the PET and recharge estimation are in Table S1 of the complementary material.

Mass Balance
The mass balance in the flow model in the transient regime and variable density conditions indicates that the inputs to the hydrogeological system are generated by recharge due to the infiltration of precipitation in the alluvial layer (2.35 hm 3 /y); lateral recharge at the edge of the aquifer (32.21 hm 3 /y); in the reservoir and the water network is represented by a mixed condition or leaked (3.68 hm 3 /y), and by a constant level that represents seawater entering the aquifer (1 hm 3 /y).
The outputs of the hydrogeological system are by pumping (34.01 hm 3 /y), in drainage in the highest aquifer area (2.63 hm 3 /y), and by the constant level that represents the aquifer discharge to the sea (1.22 hm 3 /y).

Calibration
Regarding the calibration of the flow model, the following statistical parameters intervened: (i) on the mean absolute error (ε) of the piezometric levels, estimation of the sum of the absolute differences between the calculated levels (H cal ), and the number of total observations (H obs ), (ii) the Root Mean Square Error (RMSE), and (iii) the standard deviation. The mean absolute error was 0.0082, RMSE 0.12, and the standard deviation of 0.5.
The calibration strategy consisted of a manual stage based on trial and error criteria and an automatic one provided by the Visual Transin 4.0 software, which minimizes an objective function that it considers and the piezometric adjustment between measured and calculated values.
During this process, the zoning of the lower zone of the aquifer was modified due to inconsistencies in the reproduction of the measured levels, zoning it in eight subzones of transmissivity and storage coefficients in the first layer.
The piezometric adjustment between calculated and measured levels indicates that the calibration of the parameters presents a good fit, with a linear correlation coefficient (R 2 ) of 0.9818 (Figure 7), which indicates the robustness of the numerical representation due to the conditions contour of the natural zone.
due to the conditions contour of the natural zone.
The calibrated transmissivity and storage coefficient parameters in the flow model, both in constant and variable density, are shown in Table 7. The model's calibration is based on statistical analysis with an adjustment of the piezometric levels measured in the field and those calculated in the simulation, obtaining values that are very close to each other, achieving an adjustment with reality. In the transport model, the calibration was manual based on trial and error criteria. The calibration parameters were hydrodynamic dispersion and porosity, which remained constant at 0.2. These parameters were taken into account in the first and second layers since a saline concentration source was in the basement.  Table 7. The model's calibration is based on statistical analysis with an adjustment of the piezometric levels measured in the field and those calculated in the simulation, obtaining values that are very close to each other, achieving an adjustment with reality. In the transport model, the calibration was manual based on trial and error criteria. The calibration parameters were hydrodynamic dispersion and porosity, which remained constant at 0.2. These parameters were taken into account in the first and second layers since a saline concentration source was in the basement.

Lateral Recharge
The surface runoff from the rest of the basin is represented by recharging the aquifer's edge, called lateral recharge. It is a significant water contribution to the system, being estimated by Visual Transin 4.0, which considers the transmissivity, storage, and pumping parameters, the that fit the measured levels of the 11 wells used for this simulation.
The conceptualization of lateral recharge provided the simulation with calculated levels very close to those measured, approaching the reality of this hydrogeological system's physical environment. The aquifer's edge was divided into segments, highlighting that much of the data used is from the lower area due to the study wells' locations.

Constant Density Flow Model
The aquifer's upper zone is 5 m thick, while it is 10 m thick in the middle zone. The lower zone, where the porous and permeable aquifer materials' greatest potency is, reaches thicknesses of approximately 30 m. The calibrated transmissivity and storage coefficient parameters indicate that the average permeability throughout the alluvial layer is approximately 40 m/d, with a storage coefficient of 0.2.
The calculated piezometric follows the measured piezometric. The model reproduces the hydraulic characteristics observed in the field (Figure 8). The wells closest to the sea generally have an excellent piezometric adjustment. However, during specific periods corresponding to the beginning of the 2019 rainfall season (wet), it is impossible to reproduce the observed piezometric in the wells P2 and P3, while the P4 presents the best fit. The wells inside the basins P5, P6, P7, P8, P9, and P10 present the present work's best piezometric adjustment. In wells P11 and P12, the calculated versus measured piezometric levels are acceptable (Figure 8).

Transport Model
In well P2, the field concentrations show abrupt variations, going from 1000 mg/L to 6600 mg/L in one month, making the reproduction of these concentrations by Transden

Transport Model
In well P2, the field concentrations show abrupt variations, going from 1000 mg/L to 6600 mg/L in one month, making the reproduction of these concentrations by Transden 2.0 complex. Wells P3, P4, P5, P6, and P7, present a reasonable adjustment in concentrations (TDS). Then the transport reproduces acceptably. In well P8 and P9, the calculated concentrations are below those measured, with a moderately acceptable fit. In wells P11 and P12, the calculated concentrations are above those measured, indicating that these wells have the highest saline concentration in the basin's interior because they are located very close to the middle zone, which has concentrations above 1500 mg/L in TDS (Figure 9).

Variable Density
Flow Model In all wells, except for wells P11 and P12, the calculated piezometric is close to the measurement, leading to an excellent piezometric fit ( Figure 10). In wells P11 and P12, the calculated levels are below those measured. However, the fit can be considered relatively acceptable. This model, under conditions of variable density, in a general way, reproduces the conceptualization of the underground flow in which the eleven wells present hydrodynamic activity due to pumping, causing the aquifer levels to be depressed without re-

Variable Density Flow Model
In all wells, except for wells P11 and P12, the calculated piezometric is close to the measurement, leading to an excellent piezometric fit ( Figure 10). In wells P11 and P12, the calculated levels are below those measured. However, the fit can be considered relatively acceptable. This model, under conditions of variable density, in a general way, reproduces the conceptualization of the underground flow in which the eleven wells present hydrodynamic activity due to pumping, causing the aquifer levels to be depressed without recharge due to the absence of rainfall, which means that these levels do not recover in a short time. Transport Model In this model, the variability of the density of the flows present in the basin, which affect the aquifer, is considered, generally presents a good fit (Figure 11), similar to the concentrations calculated in the transport model with constant density. The sources of concentration of the underground flow are the sea and the basement or substrate. Wells P11 and P12 have a higher concentration within the basin; this can occur due to a mixture of recharge water with concentrated water, a product of the substrate's materials. The high variability in well P2 is influenced by its proximity to the shoreline, indicating seawater intrusion. The hydrograms of the flow and transport model of the 11 wells studied in both variable and constant density are in the supplementary material section.

Transport Model
In this model, the variability of the density of the flows present in the basin, which affect the aquifer, is considered, generally presents a good fit (Figure 11), similar to the concentrations calculated in the transport model with constant density. The sources of concentration of the underground flow are the sea and the basement or substrate. Wells P11 and P12 have a higher concentration within the basin; this can occur due to a mixture of recharge water with concentrated water, a product of the substrate's materials. The high variability in well P2 is influenced by its proximity to the shoreline, indicating seawater intrusion. The hydrograms of the flow and transport model of the 11 wells studied in both variable and constant density are in the supplementary material section. Water 2021, 13, x FOR PEER REVIEW 17 of 23 Figure 11. Variable density transport model and hydrograms of wells P2, P6, P8, and P12.
In Figure 12, the system analyzed in this study is synthesized, considering its geographical and hydrogeological parts, visualizing the hydraulic effect of the infiltration produced by the tape near well P2, as well as the advance of marine intrusion into the basin and the importance of the tape in the piezometric concentrations (TDS) water and piezometric levels of the aquifer. The water concentration sources in the basement and the coastal edge by seawater affect groundwater quality for human consumption and supply. Figure 11. Variable density transport model and hydrograms of wells P2, P6, P8, and P12.
In Figure 12, the system analyzed in this study is synthesized, considering its geographical and hydrogeological parts, visualizing the hydraulic effect of the infiltration produced by the tape near well P2, as well as the advance of marine intrusion into the basin and the importance of the tape in the piezometric concentrations (TDS) water and piezometric levels of the aquifer. The water concentration sources in the basement and the coastal edge by seawater affect groundwater quality for human consumption and supply.

Discussion
The hydrograms, both of the constant density flow model and the variable density flow model, obtained from piezometric levels measured in the field and those calculated by Visual Transin 4.0, in general, reproduce the hydrodynamic behavior of the hydrogeological system. The piezometric levels of wells P2, P3, and P4 are influenced by the tape because, in this area, the piezometric levels manage to stay between zero m.a.s.l. and six m.a.s.l. Without this structure, the interval above the levels could not be maintained and would decrease considerably. The tape generates an artificial recharge, which hydraulically causes more significant infiltration to occur, achieving the piezometric levels to increase due to this structure. The closeness or distance of the tape influences the piezometric levels of the water, analogous to what occurs in the numerical model of [57], simulating the effects of artificial recharge in the riverbed, which then infiltrates into the aquifer, increasing the levels in the proximity and forming a barrier against seawater intrusion.
At the points where it is not possible to reproduce the piezometric levels, it is necessary to redefine the parameters of natural recharge and pumping in the wells within the basin, while in the wells near the sea, in addition to these parameters, the hydraulics generated by tape must be reviewed; this is because as the infiltration increases in this area, the piezometric levels' behavior will depend on the volume in which it infiltrates and periods during which this occurs. The situations are similar to the piezometric levels produced by the effects of different artificial recharge methods around the world presented by [19].
The transport models for constant and variable density show an acceptable adjustment of the concentrations. The wells inside the basin present their concentrations due to the materials that make up the aquifer substrate, which has been subject to regressions and transgressions of the sea towards the continent [38], causing evaporitic materials in the basement, which form a solution in the water. In the aquifer's lower zone

Discussion
The hydrograms, both of the constant density flow model and the variable density flow model, obtained from piezometric levels measured in the field and those calculated by Visual Transin 4.0, in general, reproduce the hydrodynamic behavior of the hydrogeological system. The piezometric levels of wells P2, P3, and P4 are influenced by the tape because, in this area, the piezometric levels manage to stay between zero m.a.s.l. and six m.a.s.l. Without this structure, the interval above the levels could not be maintained and would decrease considerably. The tape generates an artificial recharge, which hydraulically causes more significant infiltration to occur, achieving the piezometric levels to increase due to this structure. The closeness or distance of the tape influences the piezometric levels of the water, analogous to what occurs in the numerical model of [57], simulating the effects of artificial recharge in the riverbed, which then infiltrates into the aquifer, increasing the levels in the proximity and forming a barrier against seawater intrusion.
At the points where it is not possible to reproduce the piezometric levels, it is necessary to redefine the parameters of natural recharge and pumping in the wells within the basin, while in the wells near the sea, in addition to these parameters, the hydraulics generated by tape must be reviewed; this is because as the infiltration increases in this area, the piezometric levels' behavior will depend on the volume in which it infiltrates and periods during which this occurs. The situations are similar to the piezometric levels produced by the effects of different artificial recharge methods around the world presented by [19].
The transport models for constant and variable density show an acceptable adjustment of the concentrations. The wells inside the basin present their concentrations due to the materials that make up the aquifer substrate, which has been subject to regressions and transgressions of the sea towards the continent [38], causing evaporitic materials in the basement, which form a solution in the water. In the aquifer's lower zone where the study wells are located, the aquifer presents some sustainability, with an approximate power of 30 m and a saturated thickness of approximately 20 m. However, when the recharge produced by the infiltration of precipitation decreases, the saturated thickness also decreases, leading to the pumping approaching this substrate with higher concentration, producing a mixture with the aquifer water.
In the wells P2, P3, and P4, where the concentration sources are from the sea's substrate, there is evidence that these wells reach high concentrations: mainly P2, which has the highest concentrations in this system (up to 9200 mg/L). The well P12 is at the limit of the lower and middle zones of the Manglaralto river basin (see Figure 2), where the water in the mix (aquifer-basement) has the highest concentration of TDS. Furthermore, the well's depth (approximately 28 m) is very close to the basement with the highest salt concentration. The interior wells present TDS concentrations of up to 1337 mg/L, indicating the mixing of the aquifer water with saline material in the basement. On the other hand, these wells' concentrations decrease considerably in higher rainfall periods (January-May) and the garuas (July-August), which occur throughout the basin. These rain events and the tape facilitate water reservoirs that reinforce the artificial recharge that influences these wells.
The concentrations obtained by the numerical simulation indicate that marine intrusion is present in the basin's interior up to approximately 1 km from the coastline. The transition zone, seawater-groundwater aquifer, is inferred around wells P2, P3, and P4. The advance of seawater intrusion to wells' locations is due to the high pumping periods, a necessary anthropic activity that highlights the communities' development, as indicated [57] and [58], which reinforces the results obtained in the present study. The dry periods related to the decrease in recharge, which is where the importance of the tape becomes evident since, without it, the advance front of the marine intrusion would be found further into the basin, causing these wells to increase their concentration (TDS) and with this decreasing there would be a corresponding decrease in the quality of the water. That this happens is due to the prolonged contribution of infiltration by the reservoir's presence, causing a decrease in concentrations.

Conclusions
The numerical model of flow and transport under constant density and variable density shows the tape's hydraulic and environmental importance. The piezometric levels increase considerably in the closest sites in reservoir zones since artificial recharge occurs due to infiltration from the reservoir, causing the levels to predominate for a longer time. With this, the groundwater reduces the concentration (TDS), improving its quality.
The groundwater concentration inside the basin is due to the dissolution of materials present in the substrate, which mix with recharge water from the infiltration of rainfall by pumping effects. The wells closest to the sea have their concentrations due to seawater intrusion and the aquifer substratum. Therefore, in the hydrogeological system, there are two sources of concentration of groundwater.
The numerical simulation shows that marine intrusion is present, and that the advancing front into the basin is approximately in the surroundings of wells P2, P3, and P4. The tape positively impacts the forward front since the recharge is more significant in this area, preventing seawater from entering the basin's interior. Therefore, without this artificial recharge, the concentration would increase towards the other wells. The hydraulic concrete structure has been designed based on the ancestral knowledge, used by the commune members, using artisanal dams, and with the technical-scientific contribution of the academy, it has been possible to design a dam that helps in the storage of surface water and in underground recharge. An added value of the tape located in the area closest to the coastline is seawater intrusion control.  Figure S41: Well 9 hydrogram (transport model-variable density), Figure S42: Well 10 hydrogram (transport model-variable density), Figure S43: Well 11 hydrogram (transport model-variable density), Figure S44: Well 12 hydrogram (transport model-variable density), Table S1: Monthly data on precipitation and values obtained in recharge and PET, Table S2:Annual data of precipitation and values obtained in recharge and PET. Funding: The ESPOL Social Link Project funded this research: "Resilience in water management, in the face of COVID-19, Manglaralto". CODE: PG03-PY20-03 and the IAEA Project: "Improvement of coastal aquifers' management by studying the recharge rate of the alluvial aquifer of the Manglaralto river basin". CODE: ECU/7/005.