In this study, we develop a methodology to quantitatively integrate well resistivity logs and audio-magnetotelluric (AMT) sounding in the context of groundwater investigations. The experiments were conducted in a complex Quaternary depositional environment within a region of intensive groundwater use. A synthetic resistivity model was constructed from well resistivity log data and used for forward modelling to generate synthetic AMT responses, which were then inverted using Occam’s algorithm. The AMT-derived resistivity model was subsequently compared with the model obtained from the inversion of the field AMT sounding using the Pearson correlation coefficient (r) and the root mean square error (RMSE). A scalar shift factor (k) was introduced to optimize the match between both models. The comparison between the AMT-derived resistivity model and the well resistivity logs yielded a Pearson correlation coefficient of 0.669 and an RMSE of 0.1911. The optimal scalar shift factor was k = 1.1854, with a 95% confidence interval of [1.1470, 1.2246], indicating only a minor discrepancy. These results demonstrate that AMT can successfully recover a resistivity structure consistent with well resistivity logs. The proposed quantitative integration and validation workflow provides a robust framework to reduce the inherent ambiguity in AMT interpretation and highlights its potential as a direct hydrogeophysical tool for groundwater studies.

Groundwater is one of the main sources of water supply in tropical developing countries; however, its integrated management is often constrained by limited hydrogeological information and increasing anthropogenic pressures on aquifer systems. This study presents the numerical modeling of groundwater flow in the Neogene–Quaternary aquifer system of the Middle Magdalena Valley (Colombia), focusing on the rural area of Puerto Wilches, which is characterized by strong surface–groundwater interactions, particularly with the Yarirí wetland and the Magdalena River. A three-dimensional model was implemented and calibrated in FEFLOW v.8.1 under steady-state and transient conditions, integrating both primary and secondary data. The dataset included piezometric levels measured with water level meters and automatic loggers, hydrometeorological records, 21 physicochemical and microbiological parameters analyzed in 45 samples collected during three field campaigns under contrasting hydrological conditions, 79 pumping tests, detailed lithological columns from drilled wells, and complementary geological and geophysical models. The results indicate a predominant east–west groundwater flow from the Eastern Cordillera toward the Magdalena River, with seasonal recharge and discharge patterns controlled by the bimodal rainfall regime. Microbiological contamination (total coliforms in 69% of groundwater samples) and nitrate concentrations above 10 mg/L in 21% of wells were detected, mainly due to agricultural fertilizers and domestic wastewater infiltration. Particle tracking revealed predominantly horizontal flow paths, with transit times of up to 800 years in intermediate units of the Real Group and around 60 years in shallow Quaternary deposits, highlighting the differential vulnerability of the system to contamination. These findings provide scientific foundations for strengthening integrated groundwater management in tropical regions under agroindustrial and hydrocarbon pressures and emphasize the need to consolidate monitoring networks, promote sustainable agricultural practices, and establish preventive measures to protect groundwater quality.

Sustainable management of hydrogeological systems that supply water and exhibit high hydrologic complexity can be studied through pragmatic numerical modeling supported by field-constrained conceptualization. This study develops a local-scale three-dimensional groundwater flow numerical model using FEFLOW for the Barranca Lebrija settlement in Aguachica town, where the Lebrija River, the Musanda floodplain lake, and groundwater system converge. The numerical model incorporates: (i) the three-dimensional distribution of geological units and lithology; (ii) water level observations from the Musanda floodplain lake; (iii) stage records from the Lebrija River; (iv) boundary conditions and flux estimates inherited from a previous regional groundwater model; and (v) hydraulic heads from two monitoring wells and five community wells. Steady-state and transient conditions were calibrated, and a sensitivity analysis was performed to identify the parameters that most strongly control surface water–groundwater exchange. The simulations reproduce seasonal groundwater level trends and demonstrate the exchange pathways among the river, floodplain lake, and groundwater system. Results indicate dual behavior: during wet periods, flooding of the Musanda floodplain lake driven by high river levels seeps into the underlying aquifer, whereas in dry periods the floodplain lake reverses its role and becomes a principal discharge boundary. This local-scale, boundary-driven approach provides a computationally tractable framework to quantify SW–GW exchange in data-scarce tropical floodplains and supports monitoring design and water-supply management.