A Salinity Module for SWAT to Simulate Salt Ion Fate and Transport 2 at the Watershed Scale 3 4 Ryan T
Total Page:16
File Type:pdf, Size:1020Kb
Hydrol. Earth Syst. Sci. Discuss., https://doi.org/10.5194/hess-2018-614 Manuscript under review for journal Hydrol. Earth Syst. Sci. Discussion started: 22 January 2019 c Author(s) 2019. CC BY 4.0 License. 1 A Salinity Module for SWAT to Simulate Salt Ion Fate and Transport 2 at the Watershed Scale 3 4 Ryan T. Bailey1*, Saman Tavakoli-Kivi1, Xiaolu Wei1 5 1 Department of Civil and Environmental Engineering, Colorado State University, 1372 Campus Delivery, Fort Collins, CO, 6 80523-1372, United States. 7 8 *Correspondence to: Ryan Bailey ([email protected]) 9 10 Abstract. Salinity is one of the most common water quality threats in river basins and irrigated regions worldwide. However, no 11 available numerical models simulate all major processes affecting salt ion fate and transport at the watershed scale. This study 12 presents a new salinity module for the SWAT model that simulates the fate and transport of 8 major salt ions (SO4, Ca, Mg, Na, 13 K, Cl, CO3, HCO3) in a watershed system. The module accounts for salt transport in surface runoff, soil percolation, lateral flow, 14 groundwater, and streams, and equilibrium chemistry reactions in soil layers and the aquifer. The module consists of several new 15 subroutines that are imbedded within the SWAT modelling code and one input file containing soil salinity and aquifer salinity 16 data for the watershed. The model is applied to a 732 km2 salinity-impaired irrigated region within the Arkansas River Valley in 17 southeastern Colorado, and tested against root zone soil salinity, groundwater salt ion concentration, groundwater salt loadings to 18 the river network, and in-stream salt ion concentration. The model can be a useful tool in simulating baseline salinity transport 19 and investigating salinity best management practices in watersheds of varying spatial scales worldwide. 20 21 1 Introduction 22 Salinity is one of the most common water quality threats in river basins and irrigated regions worldwide. Sustainability of 23 crop production in irrigated areas in semi-arid and arid areas is threatened by over-irrigation, poor quality of irrigation water 24 (high salinity), inadequate drainage, shallow saline groundwater, and salinization of soil and underlying groundwater, all of 25 which can lead to decreasing crop yield. Of the estimated 260 million ha of irrigated land worldwide, approximately 20-30 26 million ha (7-12%) is salinized (Tanji and Kielen, 2002), with a loss of 0.25 to 0.5 million ha each year globally. Approximately 27 8.8 million ha in western Australia alone may be lost to production by the year 2050 (NLWRA, 2001), and 25% of the Indus 28 River basin is affected by high salinity. Within the western United States, 27-28% of irrigated land has experienced sharp 29 declines in crop productivity due to high salinity (Umali, 1993; Tanji and Kielen, 2002), thereby rendering irrigated-induced 30 salinity as the principal water quality problem in the semi-arid regions of the western United States. 31 Salinization of soil and groundwater systems is caused by both natural processes and human-made activities. Salt naturally 32 can be dissolved from parent rock and soil material, with salt minerals (e.g. gypsum CaSO4, halite NaCl) dissolving to mobile 2+ - + - 33 ions such as Ca , SO4 , Na , and Cl . In addition, salt ions can accumulate in the shallow soil zone due to waterlogging, which is 34 a result of over-irrigating and irrigating in areas with inadequate drainage. Salts moving up into the soil zone can become evapo- 35 concentrated due to the removal of pure water by crop roots. Soil water salinization leads to a decrease in osmotic potential, i.e. 36 the potential for water to move from soil to the crop root cells via osmosis, leading to a decrease in crop production. 37 Numerical models have been used extensively to assess saline conditions, simulate salt movement across landscapes and 38 within soil profiles, predict salt build-up and movement in the root zone, and investigate the impact of best management 39 practices (Oosterbaan, 2005; Schoups et al., 2005; Burkhalter and Gates, 2006; Singh and Panda, 2012). Available models that 40 either have inherent salinity modules or can be applied to salinity transport problems include UNSATCHEM (Šimůnek and 1 Hydrol. Earth Syst. Sci. Discuss., https://doi.org/10.5194/hess-2018-614 Manuscript under review for journal Hydrol. Earth Syst. Sci. Discussion started: 22 January 2019 c Author(s) 2019. CC BY 4.0 License. 41 Suarez, 1994), HYDRUS linked with UNSATCHEM (Šimůnek et al., 2012); DRAINMOD, LEACHC (Wagenet and Hutson, 42 1987), SAHYSMOD (Oosterbaan, 2005; Singh and Panda, 2012), CATSALT, and MT3DMS (Burkhalter and Gates, 2006). 43 Whereas several of these models include major ion chemistry for salt ions (e.g. precipitation-dissolution, cation exchange, 44 complexation) (UNSATCHEM, HYDRUS), their application typically is limited to small field-scale or soil-profile domains (e.g. 45 Kaledhonkar and Keshari, 2006; Schoups et al., 2006; Kaledhonkar et al., 2012; Rasouli et al., 2013). Conversely, models such 46 as SAHYSMOD and MT3DMS have been applied to regional-scale problems, but lack the reaction chemistry and treat salinity 47 as a conservative solute. SAHYSMOD uses seasonal water and salt balance components for large-scale systems on a seasonal 48 time step (Singh and Panda, 2012). MT3DMS is a finite-difference contaminant transport groundwater model that uses 49 MODFLOW output for groundwater flow rates, but does not include salt ion solution chemistry (Burkhalter and Gates, 2006). 50 Schoups et al. (2005) used a hydro-salinity model that couples MODHMS with UNSATCHEM to simulate subsurface salt 51 transport and storage in a 1,400 km2 region of the San Joaquin Valley, California. The model, however, does not consider 52 salinity transport in surface runoff or salt transport in streams, limiting results to soil salinity and groundwater. Currently, there is 53 no model that simulates salt transport in all major hydrologic pathways (surface runoff, soil percolation and leaching, 54 groundwater flow, streamflow) at the watershed-scale that also considers important solution reaction chemistry. Such a model is 55 important for assessing watershed-scale and basin-scale salt movement and investigating the impact of large-scale salinity 56 remediation schemes. 57 The objective of this paper is to present a salinity transport modeling code that can be used to simulate the fate and transport 58 of the major ions (SO4, Ca, Mg, Na, K, Cl, CO3, HCO3) in a watershed hydrologic system. The salinity module is implemented 59 within the SWAT modeling code, and thereby salt transport pathways include surface runoff, percolation, soil later flow, 60 groundwater flow and streamflow. The soil water and groundwater concentration of each salt ion is also affected by equilibrium 61 chemistry reactions: precipitation-dissolution, complexation, and cation exchange. The use of the model is demonstrated through 62 application to a 732 km2 region of the Lower Arkansas River Valley (LARV) in southeastern Colorado, an irrigated alluvial 63 valley in which soil and groundwater salinization has occurred over the past few decades. The model is tested against salt ion and 64 total dissolved solids (TDS) concentration in surface water (Arkansas River and its tributaries), groundwater (from a network of 65 monitoring wells), and soil water (from a large dataset of soil salinity measurements). The salinity module for SWAT can be 66 applied to any watershed to simulate baseline conditions and to test the effect of best management practices on watershed 67 salinity. 68 69 2 Development of the SWAT Salinity Module 70 This section provides a brief overview of the SWAT model, followed by a description of the SWAT salinity module. Sect. 3 71 demonstrates the use of the salinity module to a regional-scale irrigated stream-aquifer system in the Lowe Arkansas River 72 Valley, Colorado. 73 2.1 The SWAT Model 74 The SWAT (Soil and Water Assessment Tool, Arnold et al., 1998) hydrologic model simulates water flow, nutrient 75 mass transport and sediment mass transport at the watershed scale. It is a continuous, daily time-step, basin-scale, distributed- 76 parameter watershed model that simulates water flow and nutrient (nitrogen, phosphorus) transport in surface runoff, soil 77 percolation, soil later flow, groundwater flow and discharge to streams, and streamflow. The watershed is divided into subbasins, 78 which are then further divided into multiple unique combinations (Hydrologic Response Units HRUs) of land use, soil type and 79 topographic slope for which detailed water and nutrient mass balance calculations are performed. Routing algorithms route water 80 and nutrient mass through the stream network to the watershed outlet. SWAT has been applied to hundreds of watersheds and 2 Hydrol. Earth Syst. Sci. Discuss., https://doi.org/10.5194/hess-2018-614 Manuscript under review for journal Hydrol. Earth Syst. Sci. Discussion started: 22 January 2019 c Author(s) 2019. CC BY 4.0 License. 81 river basins worldwide to assess water supply and nutrient contamination under baseline conditions (Abbaspour et al., 2015) and 82 scenarios of land use change (Zhao et al., 2016; Zuo et al., 2016; Napoli et al., 2017), best management practices (Arabi et al., 83 2006; Maringanti et al., 2009; Ullrich and Volk, 2009; Dechmi and Skhiri, 2013), and climate change (Jyrkama and Sykes, 2007; 84 Ficklin et al., 2009; Tweed et al., 2009; Haddeland et al., 2010; Brown et al., 2015). However, it has not yet been applied to 85 salinity issues. 86 2.2 Salinity Module for SWAT 87 The new SWAT salinity module simulates the fate and transport of 8 major salt ions (SO4, Ca, Mg, Na, K, Cl, CO3, HCO3) 88 via surface runoff, soil later flow, soil percolation and leaching, groundwater flow, and streamflow, subject to chemical reactions 89 such as precipitation-dissolution, complexation, and cation exchange within soil layers and the alluvial aquifer.