Modeling the impact of solute recycling on groundwater salinization under irrigated lands: A study of the Alto Piura aquifer, Peru

Summary Studies of groundwater quality in arid and semi-arid lands show that irrigation return flow is one of major factors contributing to aquifer salinization. Existing mathematical models do not account explicitly for solute recycling during irrigation on a daily scale, which is considered as an important salinization input. The main objective of this research was to develop a mathematical numerical model that can simulate impact of irrigation return flow by coupling water and solute fluxes at the soil surface with quality of water pumped from the aquifer. This was obtained with a Quasi-3D model representing flow in the vadose zone – aquifer system by series of 1D Richards’ equations in a variably-saturated zone and by a 3D flow equation in groundwater. The 3D advection-dispersion equation is solved in the entire domain. Concentration of irrigation water is calculated at each time step as a function of concentration of both surface water and groundwater extracted at specific locations. The model was applied to simulate the impact of irrigation on groundwater salinization of Alto Piura aquifer (Northern Peru) over thirty years. Three scenarios were considered: (i) use of flood irrigation and groundwater extraction (the present situation); (ii) increase of groundwater pumping by 50% compared to the first scenario; and (iii) transition from flood irrigation to drip irrigation, thus decreasing irrigation volume by around 60% compared to the first scenario. Results indicate that in different irrigation areas, the simulated increase rates of total dissolved solids in groundwater vary from 3–5 to 15–17 mg/L/year, depending on hydrogeological and hydrochemical conditions, volumes of water extracted, and proportion between surface water and groundwater applied. The transition from flood irrigation to drip irrigation decreases the negative impact of return flow on groundwater quality; however drip irrigation causes faster soil salinization compared to flood irrigation. Irrigation return coefficients were calculated in the order of 21–23% and 22–24% for the first and second scenarios, respectively.

[1]  Paul P. Wang,et al.  MT3DMS: A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User's Guide , 1999 .

[2]  Andres Alcolea,et al.  Inverse problem in hydrogeology , 2005 .

[3]  J. Cavero,et al.  Nitrate exported in drainage waters of two sprinkler-irrigated watersheds. , 2003, Journal of environmental quality.

[4]  L. Prunty,et al.  Lysimeter Study of Nitrogen Fertilizer and Irrigation Rates on Quality of Recharge Water and Corn Yield , 1991 .

[5]  E. Custodio,et al.  State of knowledge of coastal aquifer management in South America , 2010 .

[6]  S. P. Neuman.  Universal scaling of hydraulic conductivities and dispersivities in geologic media , 1990 .

[7]  J. Ayars,et al.  Occurrence and measurement of salinity stratification in shallow groundwater in the Murrumbidgee Irrigation Area, south-eastern Australia , 2006 .

[8]  Jirka Šimůnek,et al.  Reply to “Comment on ‘Evaluating Interactions between Groundwater and Vadose Zone Using the HYDRUS-based Flow Package for MODFLOW’” by Navin Kumar C. Twarakavi, Jirka Šimůnek, and Sophia Seo , 2009 .

[9]  T. Jang,et al.  Estimation of irrigation return flow from paddy fields considering the soil moisture , 2009 .

[10]  Y. Rubin,et al.  Estimating the hydraulic conductivity at the south oyster site from geophysical tomographic data using Bayesian Techniques based on the normal linear regression model , 2001 .

[11]  S. Siegel,et al.  Nonparametric Statistics for the Behavioral Sciences , 2022, The SAGE Encyclopedia of Research Design.

[12]  H. Riedwyl.  Goodness of Fit , 1967 .

[13]  Raghavan Srinivasan,et al.  Return‐flow assessment for irrigation command in the Palleru river basin using SWAT model , 2005 .

[14]  M. Schaap,et al.  ROSETTA: a computer program for estimating soil hydraulic parameters with hierarchical pedotransfer functions , 2001 .

[15]  P. S. Huyakorn,et al.  A composite modeling approach for subsurface transport of degrading contaminants from land-disposal sites , 1994 .

[16]  D. Suarez.  Impact of agricultural practices on groundwater salinity , 1989 .

[17]  Leonard F. Konikow,et al.  Modeling flow and chemical quality changes in an irrigated stream‐aquifer system , 1974 .

[18]  Jirka Šimůnek,et al.  Evaluating Interactions between Groundwater and Vadose Zone Using the HYDRUS‐Based Flow Package for MODFLOW , 2008 .

[19]  Mark Person,et al.  Assessment of Long-Term Salinity Changes in an Irrigated Stream-Aquifer System , 1985 .

[20]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[21]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[22]  P. Perrochet,et al.  Direct simulation of solute recycling in irrigated areas , 2006 .

[23]  Viacheslav Borisov,et al.  A quasi three-dimensional model for flow and transport in unsaturated and saturated zones: 1. Implementation of the quasi two-dimensional case , 1998 .

[24]  D. Legates,et al.  Evaluating the use of “goodness‐of‐fit” Measures in hydrologic and hydroclimatic model validation , 1999 .

[25]  Davide Maggi,et al.  Coupled SVAT-groundwater model for water resources simulation in irrigated alluvial plains , 2004, Environ. Model. Softw..

[26]  Arlen W. Harbaugh,et al.  MODFLOW-2000, The U.S. Geological Survey Modular Ground-Water Model - User Guide to Modularization Concepts and the Ground-Water Flow Process , 2000 .

[27]  S. Shapiro,et al.  An Analysis of Variance Test for Normality (Complete Samples) , 1965 .

[28]  R. Aragüés,et al.  Irrigation management and hydrosalinity balance in a semi-arid area of the middle Ebro river basin (Spain) , 2001 .

[29]  H. Pauwels,et al.  Solute recycling: An emerging threat to groundwater quality in southern India? , 2011 .

[30]  B. Richter,et al.  Identification of sources of ground-water salinization using geochemical techniques , 1991 .

[31]  V. Zlotnik,et al.  Using pedotransfer functions in vadose zone models for estimating groundwater recharge in semiarid regions , 2009 .

[32]  L. S. Pereira,et al.  Crop evapotranspiration : guidelines for computing crop water requirements , 1998 .

[33]  G. Vachaud,et al.  Field measurements of water and nitrogen losses under irrigated maize , 1994 .

[34]  Robert L. Street,et al.  A numerical model based on coupled one‐dimensional Richards and Boussinesq equations , 1974 .

[35]  K. Pruess,et al.  TOUGH2 User's Guide Version 2 , 1999 .

[36]  B. Scanlon,et al.  Assessing controls on diffuse groundwater recharge using unsaturated flow modeling , 2005 .

[37]  James C. McWilliams,et al.  1997-1998 El Niño off Peru: A numerical study , 2008 .

[38]  E. Small.  Climatic controls on diffuse groundwater recharge in semiarid environments of the southwestern United States , 2004 .

[39]  S. Mohan,et al.  Prediction of irrigation return flows through a hierarchical modeling approach , 2009 .

[40]  Jirka Šimůnek,et al.  Inverse Dual‐Permeability Modeling of Preferential Water Flow in a Soil Column and Implications for Field‐Scale Solute Transport , 2006 .

[41]  E. Salameh,et al.  Sustainable management of groundwater resources , 2009 .

[42]  Shaul Sorek,et al.  Quasi 3D modeling of water flow in vadose zone and groundwater , 2012 .

[43]  Mark Ross,et al.  Extinction Depth and Evapotranspiration from Ground Water under Selected Land Covers , 2007, Ground water.

[44]  B. Dewandel,et al.  An efficient methodology for estimating irrigation return flow coefficients of irrigated crops at watershed and seasonal scale , 2008 .

[45]  Richard W. Healy,et al.  Factors influencing ground-water recharge in the eastern United States , 2007 .

[46]  D. Quílez,et al.  Assessment of irrigation and environmental quality at the hydrological basin level , 2004 .

[47]  Seong Taek Yun,et al.  Regional hydrochemical study on salinization of coastal aquifers, western coastal area of South Korea , 2005 .

[48]  Georges Vachaud,et al.  Comments on A numerical model based on coupled one-dimensional Richards and Boussinesq equations b , 1975 .

[49]  Philippe Renard,et al.  The problem of salt recycling and seawater intrusion in coastal irrigated plains: an example from the Kiti aquifer (Southern Cyprus) , 2004 .

[50]  Richard G. Allen,et al.  PENMAN–MONTEITH EQUATION , 2005 .