The Impact of Para Rubber Expansion on Streamflow and Other Water Balance Components of the Nam Loei River Basin, Thailand

At present, Para rubber is an economical crop which provides a high priced product and is in demand by global markets. Consequently, the government of Thailand is promoting the expansion of Para rubber plantations throughout the country. Traditionally, Para rubber was planted and grown only in the southern areas of the country. However, due to the Government's support and promotion as well as economic reasons, the expansion of Para rubber plantations in the northeast has increased rapidly. This support has occurred without accounting for suitable cultivation of Para rubber conditions, particularly in areas with steep slopes and other factors which have significant impacts on hydrology and water quality. This study presents the impacts of Para rubber expansion by applying the Soil and Water Assessment Tool (SWAT) hydrological model on the hydrology and water balance of the Nam Loei River Basin, Loei Province. The results showed that the displacement of original local field crops and disturbed forest land by Para rubber production resulted in an overall increase of evapotranspiration (ET) of roughly 3%. The major factors are the rubber canopy and precipitation. Moreover, the water balance results showed an annual reduction of about 3% in the basin average water yield, especially during the dry season.


Introduction
Zeigler et al. [1] estimated that over 500,000 ha of upland areas in southeast Asia had been converted to Para rubber (Heveabrasiliensis) production in southeast Asia by 2009 and that the land area devoted in the region to Para rubber production could double or triple by the year 2050 [1].Updated estimates for the same timeframe indicate that the expansion of total rubber production area in non-traditional Southeast Asia growing regions at >1,000,000 ha and that the production area could increase by a factor of four by 2050 [2].Much of the expansion is occurring in "marginal areas" that are vulnerable to increased soil erosion and other environmental problems [3].
Para rubber has become one of the most important economic crops in Thailand, which is now the largest exporter of Para rubber by volume worldwide [4].Para rubber production started in southern Thailand over a century ago [5] and has expanded greatly in that region since then due to favorable climatic conditions and land types.However, the government of Thailand has implemented policies to promote the expansion of Para rubber plantations throughout other areas of the country.Due to the Government's support and promotion as well as for economic reasons, the expansion of Para rubber plantations in the northeast has increased dramatically during the past decade.Continuing attractive prices have resulted in particularly rapid expansion of Para rubber plantations during the past few years, resulting an increase of nearly 500,000 ha (246,340 ha to 739,190 ha) between 2006 and 2015 [6][7][8][9], confirming earlier projections of greatly expanded production [1].The government support of Para rubber production in the northeast has occurred without adequate investigation of suitable cultivation conditions.This has resulted in Para rubber production occurring in areas with steep slopes, non-ideal climatic conditions, and other factors which have resulted in significant negative impacts on regional hydrology and water quality.
The northeastern region of Thailand consists of 20 provinces which cover a total area of about 170,226 km 2 or one-third of the country (Figure 1).Forest areas in the region are rapidly becoming degraded due to destruction of existing forest stands.This is occurring because of increased agricultural and Para rubber production to support the rapidly growing population, and burning during the summer to support wild game hunting.At present, the most extreme burning of forests in the country is occurring in north and northeast Thailand [10].
implemented policies to promote the expansion of Para rubber plantations throughout other areas of the country.Due to the Government's support and promotion as well as for economic reasons, the expansion of Para rubber plantations in the northeast has increased dramatically during the past decade.Continuing attractive prices have resulted in particularly rapid expansion of Para rubber plantations during the past few years, resulting an increase of nearly 500,000 ha (246,340 ha to 739,190 ha) between 2006 and 2015 [6][7][8][9], confirming earlier projections of greatly expanded production [1].The government support of Para rubber production in the northeast has occurred without adequate investigation of suitable cultivation conditions.This has resulted in Para rubber production occurring in areas with steep slopes, non-ideal climatic conditions, and other factors which have resulted in significant negative impacts on regional hydrology and water quality.
The northeastern region of Thailand consists of 20 provinces which cover a total area of about 170,226 km 2 or one-third of the country (Figure 1).Forest areas in the region are rapidly becoming degraded due to destruction of existing forest stands.This is occurring because of increased agricultural and Para rubber production to support the rapidly growing population, and burning during the summer to support wild game hunting.At present, the most extreme burning of forests in the country is occurring in north and northeast Thailand [10].The northeastern region of Thailand (Figure 1) has a total agricultural area of 15.90 million ha, of which 6.65 million ha are suitable for rubber plantations [11].To date, only a small portion of this potential area has been developed for Para rubber production although projections indicate greatly expanded production in the future.Investigations are urgently needed to determine how expanded rubber production in the northeast will impact environmental conditions in the region, especially rainwater, humidity, soil characteristics, hydrologic balance, flow regime and rock formation.Changes in soil quality can have a strong effect on the amount of drainable water, as well as physical, chemical and biological properties [12].Decision-makers and planners face difficult challenges in meeting water conservation objectives, and managing the engineering, socioeconomic and The northeastern region of Thailand (Figure 1) has a total agricultural area of 15.90 million ha, of which 6.65 million ha are suitable for rubber plantations [11].To date, only a small portion of this potential area has been developed for Para rubber production although projections indicate greatly expanded production in the future.Investigations are urgently needed to determine how expanded rubber production in the northeast will impact environmental conditions in the region, especially rainwater, humidity, soil characteristics, hydrologic balance, flow regime and rock formation.Changes in soil quality can have a strong effect on the amount of drainable water, as well as physical, chemical and biological properties [12].Decision-makers and planners face difficult challenges in meeting water conservation objectives, and managing the engineering, socioeconomic and environmental aspects of development and planning, related to Para rubber production in northeast Thailand.This is especially true in certain sub regions such as the Nam Loei River Basin (NLRB) in Loei province (Figure 1), where Para rubber production increased from 0.4% to 21.5% of the total land use between 2002 and 2015, resulting in an extremely volatile situation that is impacting the entire watershed.Hence, technical tools including the Soil and Water Assessment Tool (SWAT) watershed-scale water quality model [13][14][15][16] are needed to support in-depth hydrologic and environmental assessments of Para rubber production in the region.SWAT has been extensively tested for a wide range of environmental conditions and watershed scales [17][18][19][20] and has been used effectively in a number of land use change studies [21][22][23][24][25][26][27][28][29][30].Thus, the specific objectives of this research are to: (1) report the hydrologic impacts of the increased Para rubber production in the NRLB that occurred during 2002 to 2009, and 2009 to 2015; and (2) identification of inappropriate areas for rubber plantation and risks of landslide.

Description of Study Area
Loei Province covers 11,424 km 2 in the upper northeastern region of Thailand (Figure 1) and is the fifth largest province in the region.The NLRB drains 3915 km 2 from its combined upper basin and lower basin within Loei Province (Figure 1) and extends 231 km from the upper Phu Luang Range to its outlet.The Nam Puan is the major tributary of the upper basin, which is initially comprised of steep slopes but ultimately flows into a plain area where it joins the Nam Loei River within the Wang Sa Phung District.The main river of the lower basin is the Nam Loei River, which flows through the Muang District to the river plain within the Chiang Khan District to meet the Mekong River.
The average annual long-term rainfall and temperature is 1241 mm and 26.5 degree Celsius, respectively, for the NLRB [31].The range of monthly average minimum temperatures, maximum temperatures, and precipitation over the 30-year period of 1981 to 2010 are shown in Figure 2 for climate station 48353, which is located in the study region [32].However extended drought problems have resulted in streamflows of just 5% to 10% and 90% to 95% during the dry and rainy seasons, respectively, relative to annual average streamflow.Eight major groups comprise the spatial extent of soils in the NLRB, with the most dominant being the following three soil types: (1) the Slope complex (Sc) soil group which covers 44.3% of the basin, and represents a soil mixture in steep areas with >30% slopes that are generally characterized by forest, low permeability, and high risk of soil erosion; (2) the Wang hi (Wh) soil group, which covers 16.5% of the basin, represents soils derived from decay of various materials, and are characterized by fine-grain textures and high permeability; and (3) the Chiang Khan (Ch) soil group, which cover 14.4% of the basin, is derived from river sediments, and reflect sedimentary rock weathering and high permeability.Both the Wh and Ch soil types are prone to collapse or landslides in steep areas.The basin is further characterized by minimum, maximum and average elevations of 212 m, 1956 m, and 419 m, respectively.The land use of the basin, based on 2002 land use data [33], can be classified into 14 categories: corn (23.4%), disturbed forest (19.1%), forest-deciduous (12.5%), paddy field (12.3%), orchard (8.5%), sugarcane (5.9%), agricultural land-row crops (5.9%), field crop (4.8%), urban area (4.0%), miscellaneous land (1.7%), plantation (1.0%), water resources (0.4%), rubber tree (0.4%) and planted forest (0.2%).In 2002, the basin had a total Para rubber plantation area of just 762 ha.The plantation area then increased to 68,800 ha by 2009 and later to 100,000 ha by 2013, which met the expected goal established in a Para rubber production strategy for the NLRB [10].However, Para rubber production has continued to increase since 2013, reaching 129,280 ha by 2015.Currently, the 2016 provincial policy points to an expansion of an additional 100,000 ha of Para rubber plantations in the next five years [34], which will result in almost double the land area currently dedicated to Para rubber production in the NLRB by 2020.

Evaluation of Evapotranspiation (ET) at the Basin Scale
Evapotranspiration (ET) is a collective term that includes all processes by which water at the earth's surface is converted to water vapor.It includes evaporation from plant canopies, and sublimation and evaporation from soil.ET is usually the primary mechanism by which water is removed from a watershed.On average, ET equals about 62% of the precipitation that falls on landscapes across the globe except in Antarctica [35], where runoff exceeds ET.An accurate estimation of ET is critical in assessing the impact of climate and land use changes on water resources [36].The water losses through ET are significant in the hydrologic cycle; such losses are usually determined by estimating the availability of water through soil moisture or groundwater, the energy and drying power of the air, and/or via land cover and vegetation characteristics [37].A new method of determining ET for Para rubber trees has been developed in which the rubber tree ET is estimated by accounting for observed patterns of rubber root water uptake as affected by the plant's phrenology [38].Specifically, the method considers vegetation dynamics and corresponding water needs or evaporative demands.This contrasts with the traditional approach of estimating Para rubber ET, which neglects the increased water use during the dry season when both soil water content and canopy cover are minimal [38].It is expected here that the SWAT basin-scale hydrologic model will more accurately capture seasonal water balance and ET dynamics that are more consistent with recent research, especially for larger scale rubber expansion situations that exceed smaller stand levels.

Description of SWAT Model
SWAT is a public domain model jointly developed by the U.S. Department of Agriculture Agricultural Research Service (USDA-ARS) and Texas AgriLife Research, a unit within the Texas A&M University System [15][16][17].Watersheds are typically simulated in SWAT by delineating the respective watershed into subbasins and then further subdividing each subbasin into hydrologic response units (HRUs), which are non-spatial land units consisting of homogeneous topographic, soil type, land use, and management characteristics.Hydrologic cycling including precipitation inputs, surface runoff, infiltration into the soil profile, ET, lateral subsurface flow, and flow via other pathways is initially simulated at the HRU level.Nutrient cycling and transport, as well as sediment losses, are also simulated first at the HRU scale.The HRU-level hydrologic and pollutant outputs are then aggregated to the subbasin level and ultimately routed through the stream network to the watershed outlet.
SWAT first simulates atmospheric water demands to calculate the maximum, unstressed ET, commonly referred to as potential ET, before calculating the final actual ET.Three potential ET methods are included in SWAT that vary considerably in the amount of required inputs: (1) Penman-Monteith [39], which requires solar radiation, air temperature, relative humidity and wind speed; (2) Priestley-Taylor [40], which requires solar radiation, air temperature and relative humidity; and (3)

Evaluation of Evapotranspiation (ET) at the Basin Scale
Evapotranspiration (ET) is a collective term that includes all processes by which water at the earth's surface is converted to water vapor.It includes evaporation from plant canopies, and sublimation and evaporation from soil.ET is usually the primary mechanism by which water is removed from a watershed.On average, ET equals about 62% of the precipitation that falls on landscapes across the globe except in Antarctica [35], where runoff exceeds ET.An accurate estimation of ET is critical in assessing the impact of climate and land use changes on water resources [36].The water losses through ET are significant in the hydrologic cycle; such losses are usually determined by estimating the availability of water through soil moisture or groundwater, the energy and drying power of the air, and/or via land cover and vegetation characteristics [37].A new method of determining ET for Para rubber trees has been developed in which the rubber tree ET is estimated by accounting for observed patterns of rubber root water uptake as affected by the plant's phrenology [38].Specifically, the method considers vegetation dynamics and corresponding water needs or evaporative demands.This contrasts with the traditional approach of estimating Para rubber ET, which neglects the increased water use during the dry season when both soil water content and canopy cover are minimal [38].It is expected here that the SWAT basin-scale hydrologic model will more accurately capture seasonal water balance and ET dynamics that are more consistent with recent research, especially for larger scale rubber expansion situations that exceed smaller stand levels.

Description of SWAT Model
SWAT is a public domain model jointly developed by the U.S. Department of Agriculture Agricultural Research Service (USDA-ARS) and Texas AgriLife Research, a unit within the Texas A&M University System [15][16][17].Watersheds are typically simulated in SWAT by delineating the respective watershed into subbasins and then further subdividing each subbasin into hydrologic response units (HRUs), which are non-spatial land units consisting of homogeneous topographic, soil type, land use, and management characteristics.Hydrologic cycling including precipitation inputs, surface runoff, infiltration into the soil profile, ET, lateral subsurface flow, and flow via other pathways is initially simulated at the HRU level.Nutrient cycling and transport, as well as sediment losses, are also simulated first at the HRU scale.The HRU-level hydrologic and pollutant outputs are then aggregated to the subbasin level and ultimately routed through the stream network to the watershed outlet.
SWAT first simulates atmospheric water demands to calculate the maximum, unstressed ET, commonly referred to as potential ET, before calculating the final actual ET.Three potential ET methods are included in SWAT that vary considerably in the amount of required inputs: (1) Penman-Monteith [39], which requires solar radiation, air temperature, relative humidity and wind speed; (2) Priestley-Taylor [40], which requires solar radiation, air temperature and relative humidity; and (3) Hargreaves [41], which requires air temperature only [11].Multiple options are also provided in SWAT for simulating the partitioning of precipitation inputs between surface runoff and infiltration, as well as for some other processes simulated in the model.Complete theoretical and user input options are provided in the SWAT model documentation [36].A revision of SWAT version 2009 (SWAT2009) was used in conjunction with the ArcGIS SWAT (ArcSWAT) interface for this study [42].

Data Input Needs and Sources
Topographic, soil, land use, climate, and management data are key inputs required for simulating a watershed in SWAT.In-stream monitoring data are also important in regards to testing SWAT output.Topographic, soil, and land use are usually input into SWAT in the form of digital spatial layers that are overlaid within a Geographic Information System (GIS).All of these major data were prepared in input format files for the SWAT simulations including spatial topographic, soil and land use data (Table 1, Figures 3 and 4) obtained from various sources as described below.
Water 2016, 8, 1 5 of 19 Hargreaves [41], which requires air temperature only [11].Multiple options are also provided in SWAT for simulating the partitioning of precipitation inputs between surface runoff and infiltration, as well as for some other processes simulated in the model.Complete theoretical and user input options are provided in the SWAT model documentation [36].A revision of SWAT version 2009 (SWAT2009) was used in conjunction with the ArcGIS SWAT (ArcSWAT) interface for this study [42].

Data Input Needs and Sources
Topographic, soil, land use, climate, and management data are key inputs required for simulating a watershed in SWAT.In-stream monitoring data are also important in regards to testing SWAT output.Topographic, soil, and land use are usually input into SWAT in the form of digital spatial layers that are overlaid within a Geographic Information System (GIS).All of these major data were prepared in input format files for the SWAT simulations including spatial topographic, soil and land use data (Table 1, Figures 3 and 4) obtained from various sources as described below.1).
Spatial topographic data required for the SWAT application were obtained in the form of a digital elevation map (DEM) [43].These DEM data are characterized by a 30 m × 30 m (1:50,000 scale) resolution.The baseline land use data for year 2002 [44] and the spatial soil map were provided by the Land Development Department [45].The LDD soil data consist of the previously mentioned 8 major groups.All soil properties required for SWAT were surveyed from by the LDD in 2012 [46].Daily climate data were used, including precipitation, temperature, solar radiation, wind speed, and humidity data.These data were collected from rainfall gauges within the NLRB during the period from 1985 to 2015 [32].The daily precipitation data were obtained at 14 gauging stations (Figure 1), while the daily temperature, solar radiation, wind speed and humidity data were collected from a single major weather station (Figure 1) of the Thai Meteorological Department (TMD).Two  1).
Spatial topographic data required for the SWAT application were obtained in the form of a digital elevation map (DEM) [43].These DEM data are characterized by a 30 m × 30 m (1:50,000 scale) resolution.The baseline land use data for year 2002 [44] and the spatial soil map were provided by the Land Development Department [45].The LDD soil data consist of the previously mentioned 8 major groups.All soil properties required for SWAT were surveyed from by the LDD in 2012 [46].Daily climate data were used, including precipitation, temperature, solar radiation, wind speed, and humidity data.These data were collected from rainfall gauges within the NLRB during the period from 1985 to 2015 [32].The daily precipitation data were obtained at 14 gauging stations (Figure 1), while the daily temperature, solar radiation, wind speed and humidity data were collected from a single major weather station (Figure 1) of the Thai Meteorological Department (TMD).Two hydrologic runoff stations located in the basin (gauges Kh.28A and Kh.58A in Figure 1) that are maintained by the Royal Irrigation Department (RID) and have complete monthly runoff data were selected for model calibration and validation [47].hydrologic runoff stations located in the basin (gauges Kh.28A and Kh.58A in Figure 1) that are maintained by the Royal Irrigation Department (RID) and have complete monthly runoff data were selected for model calibration and validation [47].1).

Model Set Up
The NLRB was delineated into 19 subbasins and 389 hydrologic response units (HRUs) for the SWAT simulations using the ArcSWAT Interface [42].The delineation was performed as a function of the DEM-based surface topography, which resulted in the configuration of the 19 sub basins used in the SWAT simulations.The land use data were processed and reclassified to match the land use codes used in SWAT, resulting in the previously described 14 land use categories.The eight major soil types were also converted and reclassified to match the SWAT model soil formatting requirements.These land use and soil data were then overlaid with the DEM data within ArcSWAT to create the HRUs.The required weather data were incorporated or the simulations.The initial curve number values were assigned based on the land use type and soil hydrologic group for the average antecedent moisture condition of the runoff curve number method.The PET was computed by using the Penman-Monteith method.The overall model simulation scenario period covered a 25-year duration from 1985 to 2009 using a daily time step, with a shorter time period used for model calibration and validation as described below.

Sensitivity Analysis and SWAT Calibration and Validation
A sensitivity analysis is a useful procedure to determine which flow-related parameters are the most influential in impacting total streamflow for SWAT applications, following guidance reported for previous studies [16,[48][49][50][51]. Performing a sensitivity analysis further supports the calibration of SWAT and provides insight for the application of the model to other similar watersheds.In this study, the SWAT CUP software package [52] was used to perform an automatic sensitivity analysis of the impact of 19 different SWAT parameters on daily streamflow flow for the NLRB as described in the Results and Discussion section.
Following the sensitivity analysis, calibration and validation of SWAT was performed which is required to reduce uncertainty and increase confidence in its predictive abilities for the NRLB [16,49,50].The calibration process included both multisite and multivariable aspects as discussed in previous SWAT studies [16,50].The calibration procedure consisted of 3 stages: (1) replicating the long-term water balance over the calibration period; (2) accurately tracking the observed hydrograph shapes; and (3) obtaining an accurate comparison between observed and simulated flow duration curves.Calibration and validation of SWAT was performed by comparing the simulated monthly aggregated stream flows versus corresponding measured monthly stream flows at two hydrological gauge stations on the main stem of the Nam Loei River (Figure 1): Wang SaPhung (Kh.28A) in subbasin 14 and Ban FakLoei (Kh.58A) in subbasin 6. Calibration was performed for 1994 to 2004 while validation was conducted from 2005 to 2009.The parameters derived for the gauged catchments were then transferred to the ungauged catchments, based on proximity and similarities in land use and soil types which result in similar hydrological responses.
The accuracy of the model output variance was assessed using the Root Mean Squared Error (RMSE) and Nash-Sutcliffe Efficiency (NSE) statistics [53,54], which are expressed as follows in Equations ( 1) and ( 2): where Q obs are the observed values and Q sim are the simulated values at time/place i.
Values of RMSE can range from 0 to ∞ where a value of 0 indicates a perfect fit between the simulated data and counterpart measured data [53].The RMSE provides a measure of the difference between the measured and simulated values or residual variance [55].However, statistical results generated with the RMSE can result in significant model bias during model calibration, even when the error variances are small [56].Thus it is desirable to use additional statistical evaluation when using the RMSE such as the NSE.
The NSE estimates the magnitude of the simulated variance relative to the measured variance and how accurately a plot of the modeled versus measured data fit a 1:1 line.The NSE can vary from −∞ to 1, where 1 is a perfect fit and a negative value indicates that the average value of the measured data would provide a better prediction than the simulated data.Statistical criteria for judging the success of hydrological modeling results has been suggested [55,57] including an NSE value of at least 0.5 to achieve a satisfactory level for comparisons of aggregate monthly simulation output versus corresponding measured stream flow data.Following successful calibration and validation, the calibrated SWAT model was used to evaluate the Para rubber land use scenario (describe below) including evaluations of ET and water yield.

Development of Para Rubber Land Use Scenarios
Several factors need to be considered in developing the Para rubber land use scenarios for the NLRB.First, it is important to consider the optimal climatic, soil, slope, and other conditions that Para rubber should be grown under for the study region, especially in the context of typical current practices in which rubber plantations are being established on very high, vulnerable slopes.Second, crop parameters needed to be developed for Para rubber for this analysis.Third, 2002 baseline, 2009 scenario and 2015 scenario landuse layers had to be constructed in order to simulate the impact of expanded Para rubber production in the NLRB during the period of rapid production expansion.These aspects of the Para rubber scenario development are described below.

Optimal Environmental Conditions for Para Rubber Production
The rubber tree is native to the evergreen tropical rainforests which usually occur within 5 • latitude of the equator.The climate of this region is characterized by heavy rainfall and no distinct dry season [5].The optimal climatic conditions for Para rubber include rainfall of 1250 mm or more, evenly distributed throughout the year with no severe dry season and with 120-150 annual rainy days, and a temperature range of about 26-30 • C [34].In addition, the ideal elevation range for Para rubber growth is from sea level up to 600 m above mean sea level; the growth rate will decline at higher altitudes.Para rubber can grow on many soils, with the best options being well drained clayey and deep clay soils, but it can withstand physical conditions ranging from stiff clays with poor drainage to well drained sandy loams [5].The most suitable soil conditions for Para rubber production are: (1) planted at a depth which does not exceed 1 m in depth to allow for adequate future root penetration and growth; (2) soil textures that range between loamy sand to clay loam with good drainage; (3) no gravel or stone in the subsoil layer; and (4) the soil pH should range between 4.5 and 5.5 [34].
The ideal slope range for growing Para rubber trees is between 5% and 15%.Three categories can be classified for this ideal range for the Loei River basin region: (1) low land areas with slopes from 0% to 5%; (2) slopes ranging from 5% to 10% located in the plain areas; and (3) mild slopes that typically range from 10% to 15%.About 30% (117,433 ha) of the total area within the basin region meets these ideal slope characteristics [58].However, the majority of the current Para rubber production in the study region is currently occurring on much steeper slopes exceeding 15%, and is concentrated especially in an extreme slope range of 30% to 35%.Therefore, Para rubber production on slopes greater than 15% should be managed with terraces as shown in the photos and schematic in Figure 5.It was outside the scope of the present research to account for such terrace systems as part of the SWAT simulations reported in this study.

Para Rubber Crop Parameters
Crop growth in the SWAT model was set up for this study primarily using a heat unit schedule approach.However, rice and other "field crops" were simulated for two growing seasons during each year and thus both were simulated using specific planting and harvesting dates rather than heat unit scheduling.Double cropping of wet season rice and dry season rice were accounted for in the simulation; the second dry season crop is mostly grown in selected areas of Thailand that have sufficient irrigation water available.The planting schedules of wet season rice, dry season rice and field crops were simulated based on typical planting dates for each crop type in the region.The remaining operations were controlled by the fraction of heat units for each crop using the heat unit scheduling approach [36].
Para rubber is a perennial tree crop which grows year round; thus, annual planting and harvest and operations were not simulated for Para rubber.Para rubber crop parameters were not available in the ArcSWAT database for SWAT2009; thus, parameters had to be determined from measured data, inferred from existing parameters for other tree species in the ArcSWAT database, or determined from other sources.Values for 15 key Para rubber Tree (RUBR) crop parameters, and two other vegetation-influenced input parameters, which were used for the NLRB SWAT analysis, are listed in Table 2.The value selected for the maximum canopy height (CHTMX) were based on previously reported measurements [5] while the maximum rooting depth (RDMX) value is based on other measurements conducted by Thai scientists [58].The other crop parameters in Table 2 were derived mainly from existing tree crop parameters in the ArcSWAT database.The choice of curve number (CN2) value for Para rubber trees reflects a woodland condition consisting of a thin stand, poor cover, no mulch, and a soil type consistent with hydrologic soil group B drainage conditions [59].The Manning's n value for overland flow (OV_N) represents high runoff for timberland conditions [60].Minimum, maximum and average values are listed for each crop parameter in Table 2.These parameter ranges allow the user the flexibility to decrease or increase the crop parameter default values, typically by approximately 10%.
Finally, the daily rubber tree water requirement rate was calculated by using the Bowen Ratio method based on data collected for 10-year old rubber trees at an experimental site in Chachoengsao Province.Temperature and humidity sensors were installed within and above the rubber tree canopy at the site.Solar radiation and wind velocity sensors were also installed above the rubber tree canopy.Soil temperature and soil moisture sensors were installed at 0 m and 0.10 m depths.A S-shape regression had been proposed to describe the relation between crop coefficients (Kc) and the Julian date.Therefore, the S-shape regression can be applied to evaluate the water requirement of rubber trees at different locations beyond the experimental site [61].

Para Rubber Crop Parameters
Crop growth in the SWAT model was set up for this study primarily using a heat unit schedule approach.However, rice and other "field crops" were simulated for two growing seasons during each year and thus both were simulated using specific planting and harvesting dates rather than heat unit scheduling.Double cropping of wet season rice and dry season rice were accounted for in the simulation; the second dry season crop is mostly grown in selected areas of Thailand that have sufficient irrigation water available.The planting schedules of wet season rice, dry season rice and field crops were simulated based on typical planting dates for each crop type in the region.The remaining operations were controlled by the fraction of heat units for each crop using the heat unit scheduling approach [36].
Para rubber is a perennial tree crop which grows year round; thus, annual planting and harvest and operations were not simulated for Para rubber.Para rubber crop parameters were not available in the ArcSWAT database for SWAT2009; thus, parameters had to be determined from measured data, inferred from existing parameters for other tree species in the ArcSWAT database, or determined from other sources.Values for 15 key Para rubber Tree (RUBR) crop parameters, and two other vegetation-influenced input parameters, which were used for the NLRB SWAT analysis, are listed in Table 2.The value selected for the maximum canopy height (CHTMX) were based on previously reported measurements [5] while the maximum rooting depth (RDMX) value is based on other measurements conducted by Thai scientists [58].The other crop parameters in Table 2 were derived mainly from existing tree crop parameters in the ArcSWAT database.The choice of curve number (CN2) value for Para rubber trees reflects a woodland condition consisting of a thin stand, poor cover, no mulch, and a soil type consistent with hydrologic soil group B drainage conditions [59].The Manning's n value for overland flow (OV_N) represents high runoff for timberland conditions [60].Minimum, maximum and average values are listed for each crop parameter in Table 2.These parameter ranges allow the user the flexibility to decrease or increase the crop parameter default values, typically by approximately 10%.
Finally, the daily rubber tree water requirement rate was calculated by using the Bowen Ratio method based on data collected for 10-year old rubber trees at an experimental site in Chachoengsao Province.Temperature and humidity sensors were installed within and above the rubber tree canopy at the site.Solar radiation and wind velocity sensors were also installed above the rubber tree canopy.Soil temperature and soil moisture sensors were installed at 0 m and 0.10 m depths.A S-shape regression had been proposed to describe the relation between crop coefficients (Kc) and the Julian date.Therefore, the S-shape regression can be applied to evaluate the water requirement of rubber trees at different locations beyond the experimental site [61].

Land Use Change Scenarios
Table 3 shows the specific land use distributions for the 2002, 2009 and 2015 land use scenarios, and the percentage difference for each land use category for two time periods: (1) between 2002 and 2009; and (2) between 2009 and 2015.The major changes in land use between 2002 and 2009 included a nearly 11.5% increase in Para rubber production, a decline in corn production of 10%, a decrease in disturbed forest land of over 9.7%, and increases in evergreen and deciduous forest of 7.35% and 2.1%, respectively.The largest shift in land use between 2009 and 2015 was an increase in Para rubber production of 9.69%, which was primarily responsible for respective decreases of −3.57%, −2.99%, −1.93%, −1.84% and −1.77% of corn, disturbed forest land, sugarcane, paddy fields and orchards.

Sensitivity Analysis
The top five most sensitive parameters as ranked in Table 4 were: (1) ALPHA_BF, base flow alpha factor (days); (2) ESCO, soil evaporation compensation factor; (3) GQWMN, threshold depth of water in the shallow aquifer required for return flow to occur (mm); (4) CN2, initial SCS runoff curve number for moisture condition II; and (5) CH_K2, effective hydraulic conductivity in main channel alluvium (mm•h −1 ).The most influential parameters found in the sensitivity analysis are consistent with previously published summaries of the most widely used parameters in SWAT calibration [16,50].The results also underscore the importance of accurate spatial and temporal precipitation inputs [16,50].The choice of ALPHA_BF, CN2, ESCO and other parameters also varied (Table 4) between the two subbasins that drain to gauges Kh.28A and Kh.58A, respectively.

Sensitivity Analysis
The top five most sensitive parameters as ranked in Table 4 were: (1) ALPHA_BF, base flow alpha factor (days); (2) ESCO, soil evaporation compensation factor; (3) GQWMN, threshold depth of water in the shallow aquifer required for return flow to occur (mm); (4) CN2, initial SCS runoff curve number for moisture condition II; and (5) CH_K2, effective hydraulic conductivity in main channel alluvium (mm•h −1 ).The most influential parameters found in the sensitivity analysis are consistent with previously published summaries of the most widely used parameters in SWAT calibration [16,50].The results also underscore the importance of accurate spatial and temporal precipitation inputs [16,50].The choice of ALPHA_BF, CN2, ESCO and other parameters also varied (Table 4) between the two subbasins that drain to gauges Kh.28A and Kh.58A, respectively.

Model Calibration and Validation
The statistical results of comparing the simulated SWAT calibration and validation aggregated monthly streamflows versus corresponding measured streamflows are listed in Table 5 for both gauge sites (Figure 1).The graphical comparisons between the simulated and measured monthly streamflows for the calibration and validation period are shown in Figure 7 for gauge site Kh.58A.The four NSE values computed for the two gauges during the calibration and validation period all exceeded 0.5, indicating that SWAT produced satisfactory streamflow estimates per previously suggested criteria [55,57].The RMSE statistics were also all below 1.0 which further confirm that the SWAT streamflow estimates satisfactorily replicated the measured streamflows.The graphical comparisons between the simulated and measured aggregated monthly streamflows for the calibration and validation period at gauge site Kh.58A (Figure 7) show that SWAT accurately replicated most of the measured streamflow trends during the 11-year calibration period.However, several peak monthly streamflows were under predicted, especially in the last three years of the calibration period, which mirrors a tendency towards under prediction reported in a number of existing SWAT studies (e.g., [15]) and points to the need for further improvement of the SWAT hydrological algorithms, especially for Southeast Asia conditions.Similar graphical results occurred for the other three gauge site/time period combinations and thus are not reported here.
Water 2016, 8, 1 13 of 19 of existing SWAT studies (e.g., [15]) and points to the need for further improvement of the SWAT hydrological algorithms, especially for Southeast Asia conditions.Similar graphical results occurred for the other three gauge site/time period combinations and thus are not reported here.

Overall Water Balance Results for the Land Use Scenarios
Table 6 shows the overall long-term average annual water balance results for the entire NLRB predicted by SWAT for the 25-year 2002, 2009 and 2015 land use scenario simulations.The results of the scenarios reflect the increases in Para rubber production and other shifts in land use that occurred during the 2002 to 2009 and 2002 to 2015 time periods.Transmission losses were essentially negligible for all three land use scenarios, and the estimated combined lateral subsurface flow and groundwater flow were very similar between the three scenario simulations.However, the predicted ET increased by nearly 17 mm, and the predicted surface runoff and water yield decreased by similar amounts for the 2009 land use scenario as compared to the 2002 land use scenario.Similar, greater decreases in reduced surface runoff and water yield, relative to the 2009 land use scenario, were estimated for the 2015 land use scenario.However, the simulated ET decreased by almost 12 mm between the 2015 and 2009 land use scenarios.These results underscore the impacts of both the increased Para rubber production and shifting overall land use mixes in the NLRB between 2002 and 2015 (Table 3).3) resulted in predicted increases in ET in almost every month of the year,

Overall Water Balance Results for the Land Use Scenarios
Table 6 shows the overall long-term average annual water balance results for the entire NLRB predicted by SWAT for the 25 3).3) resulted in predicted increases in ET in almost every month of the year, relative to the baseline year of 2002 (Figure 8), except for March and April.The annual average ET increased about 3% (Table 6), with the highest percentage increases occurring during the dry season months of November to January and the lowest percentage increases occurring during the wet season months of May, September and October.The estimated water yield responses between the baseline and 2009 Para rubber expansion scenario resulted in the opposite trend, with water yields decreasing in most months, although slight increases occurred during the months of November, December and January.The predicted percentage water yield changes ranged from +0.6% in January to over −10% in March and April.The percentage declines in water yield predicted for the wet season were more constant as compared to the dry season water yield impacts, ranging from −1% in October to −7% in October.
The continued expansion of Para rubber production and other land use changes between 2009 and 2015 (Table 3) resulted in stronger shifts in the predicted monthly ET levels and in the annual hydrograph (Figure 8), with an earlier onset of the flood season and a decreased overall peak discharge in September.The estimated ET for the 2015 land use was higher during the dry season as compared to 2002 and 2009 (Figure 8), except for November, but the opposite trend occurred during the wet season, resulting in an overall decline of about 2% from 2009 to 2015 (Table 6).The overall average dry period water yield was predicted to be about 13% higher for the 2015 land use scenario as compared to the baseline during November to April.However, lower water yields were predicted during most of the wet period except for the months of July and August (Figure 8).In total, the average annual simulated water yield decreased almost 7% from the baseline to 2015 (Table 6).3) resulted in stronger shifts in the predicted monthly ET levels and in the annual hydrograph (Figure 8), with an earlier onset of the flood season and a decreased overall peak discharge in September.The estimated ET for the 2015 land use scenario was higher during the dry season as compared to 2002 and 2009 (Figure 8), except for November, but the opposite trend occurred during the wet season, resulting in an overall decline of about 2% from 2009 to 2015 (Table 6).The overall average dry period water yield was predicted to be about 13% higher for the 2015 land use scenario as compared to the baseline during November to April.However, lower water yields were predicted during most of the wet period except for the months of July and August (Figure 8).In total, the average annual simulated water yield decreased almost 7% from the baseline to 2015 (Table 6).7), as compared to the dry season.However, a more even distribution of ET occurred between the wet and the dry seasons, with slightly Table 7 provides comparisons of: (1) the long-term (25-year simulation period) average cumulative amounts of ET (mm) and water yield (mm) that occurred during the wet season, dry season, and annually; (2) the percentages of the cumulative amounts that occurred during the wet season and dry season; and (3) the percentage changes in annual average ET and WYLD that occurred due to the land use changes between 2002 and 2009, and 2009 and 2015.These results show that nearly 80% of the water yield occurred during the wet season during both the 2002 baseline and the 2009 and 2015 Para rubber expansion land use scenarios (Table 7), as compared to the dry season.However, a more even distribution of ET occurred between the wet and the dry seasons, with slightly higher levels occurring during the dry season.Slight shifts in the overall amounts that occurred between the two seasons were generally predicted the 2009 and 2015 conditions relative to the 2002 baseline; the 2015 land use scenario ET estimates resulted in the largest relative shift that was predicted between the two seasons (Table 7).
Various seasonal shifts were predicted to occur between 2002 baseline and the 2009 landuse scenarios, and between the 2009 and 2015 land use scenarios (Table 7).Relatively minor shifts were estimated between 2002 and 2009, with ET increasing by roughly 3% in both seasons versus 1% and 3% declines in water yield for the dry season and wet season, respectively.ET was predicted to decrease almost 10% during the wet season but increase nearly 6% during the dry season, between 2009 and 2015.In contrast, water yield was predicted to decrease during the wet season by about 3% and increase by 1% during the dry season, between the two landuse scenarios.The estimated average annual percentage change in ET from 2002 to 2009 was a decline of roughly 3% as compared to an increase in water yield of close to 3% for the same time period (Table 7).In contrast, the predicted average annual percentage change in both ET and water yield from 2009 to 2015 was a decline of roughly 2% (Table 7).

Discussion
As noted previously, substantial increases in rubber tree production are expected to occur in Southeast Asia by 2050.Preliminary research suggests that this massive land use change could exacerbate environmental problems in the region including increased soil erosion and sediment transport to surface water, degraded soil quality, decreased stream flow, risk of landslides and probable decreased soil carbon levels [1,2,62].The majority of the area where rapid and widespread land conversion to monoculture rubber plantation has occurred in continental SE Asia is also vulnerable to extreme climatic events including typhoons, frost or drought, which can greatly reduce or even destroy rubber production and further exacerbate environmental problems [3].Furthermore, these environmental problems may be magnified even more per projected future climate change in the region [3].Intensified Para rubber production has also negatively impacted biodiversity in production areas located in Thailand and other subregions in Southeast Asia [63].
The results of this study show that the expansion of rubber cultivation is resulting in a decreased volume of water in the rainy season in the NLRB in northeast Thailand.In addition, during the dry season the water content decreases, resulting in water shortages.The expansion of rubber production in inappropriate areas with slopes over 35% [51] is more likely to result in increased flash floods and landslides due to the intensive monoculture practices.This is especially true for the Nam Manh, Nam Phu and Nam Paow subbasins, where 3000 ha of Para rubber plantations have been introduced on landscapes with extremely high slopes.In addition, an additional 10,300 ha of Para rubber production exists in more moderate slope areas (15% to 35% slopes) that are distributed across the upper, middle and lower subregions of the NLRB.The Rubber Research Institute of Thailand has recommended that terrace systems should be used for rubber production on slopes >15% [34] but this recommendation is not being consistently followed.These specific problems occurring in Thailand, coupled with the previously described environmental problems that are increasing Southeast Asia, underscore the urgent need to develop measurement approaches, databases and modeling tools that can be used to investigate Para rubber production problems throughout the region.
The amount of water that is required to initiate and sustain Para rubber plant growth depends on many factors such as the type of plant species, the age of the plant and weather conditions such as wind speed, temperature and humidity [61].There are many ways to determine crop water use both directly and indirectly.The measure by instrumentation is more accurate, but cannot be conducted across large plantations.The application of mathematical models is a popular method and are widely used in the study and evaluation of water use and water requirements of crops, which saves time and cost.However, the use of mathematical models can require extensive input data depending on the type and format of the model.There are parameters that must be calibrated and verified to ensure the calculated results are close to the measured values.Furthermore, databases with appropriate parameter values may not be available for some models.Thus, those databases must be brought up-to-date to better meet the needs of end users.
The case study reported here that describes the expansion of Para rubber plantations in the NLRB using SWAT is an example of an application that requires the development of important plant input parameters.However, there is a need to further develop SWAT Para rubber input parameters that better account for rubber species and age, rubber stems and leaves of rubber trees.In addition, there is a need to improve the SWAT growth functions to be able to better account for the effects of rubber tree growth phenomena related to the rooting structure that depletes deeper soil layers and results in higher ET impacts, relative to traditional vegetation [36,64].Finally, an overall expanded Para rubber production and knowledge database is needed for Thailand and Southeast Asia in general.

Conclusions
The Soil and Water Assessment Tool (SWAT) model was applied to assess the impact of Para rubber expansion in the Nam Loei River Basin in this study The application has designated the land use in the year 2002 as the baseline versus historical Para rubber tree expansions in 2009 and 2015 as the land use change scenarios.The stream flow estimated by SWAT showed annual average stream flow of about 1580 MCM.The average stream flow occurring during the rainy season (June-November) was about 1264 MCM (80% of the average annual stream flow) and in the dry season about 316 MCM (20% of the average annual stream flow).The simulation was done using scenarios to assess the water balance in the hydrological process.The results of the simulations showed that the increased production of Para rubber, which replaced the original local field crop and disturbed forest land, resulted in an increase of ET of about 3% from 2002 to 2009.However, additional increased Para rubber production in combination with other land use shifts during 2009 to 2015 resulted in a predicted ET decrease of 2%.The major factors that influenced this result were the rubber canopy and precipitation.Moreover, runoff results reduced water balance in the basin by an annual average of about 3%, especially during the dry season.However, the effect of ET on water resources has increased complexity and uncertainty; the consideration of many parameters of Para rubber and a reflection on the past will help our understanding of the dynamic changes.The results of this study will help provide guidance for decision-making about land use allocation or zoning for suitable Para rubber area.In addition, this study will aid in the management and planning of water resources for the NLRB and other river basins located in northeastern Thailand.

Figure 1 .
Figure 1.Location of the Nam Loei River Basin (NLRB) within Loei Province and Loei Province within northeast Thailand.

Figure 1 .
Figure 1.Location of the Nam Loei River Basin (NLRB) within Loei Province and Loei Province within northeast Thailand.

Figure 2 .
Figure 2. Range of average monthly minimum temperatures, maximum temperatures and precipitation at climate station 48353 during the 30-year period 1981-2010 [30].

Figure 3 .
Figure 3. Elevation ranges and distribution of soil types for the Nam Loei River Basin (NLRB), which are based on 1:50,000 Digital Elevation and soil maps, respectively (Table1).

Figure 3 .
Figure 3. Elevation ranges and distribution of soil types for the Nam Loei River Basin (NLRB), which are based on 1:50,000 Digital Elevation and soil maps, respectively (Table1).
Notes: a DWR = Department of Water resources; LDD = Land Development Department; RID = Royal Irrigation Department; RTSD = Royal Thai Survey Department; TMD = Thai Meteorological Department.
Notes: a DWR = Department of Water resources; LDD = Land Development Department; RID = Royal Irrigation Department; RTSD = Royal Thai Survey Department; TMD = Thai Meteorological Department.

Figure 4 .
Figure 4. Distribution of land use in 2002 for the Nam Loei River Basin (NLRB) based on 1:50,000 land use/land cover map (Table1).

Figure 4 .
Figure 4. Distribution of land use in 2002 for the Nam Loei River Basin (NLRB) based on 1:50,000 land use/land cover map (Table1).

Figure 5 .
Figure 5. Photos and schematic showing terrace system needed to help mitigate excessive runoff and soil erosion that can occur when Para rubber trees are grown on high slopes >15%.

Figure 5 .
Figure 5. Photos and schematic showing terrace system needed to help mitigate excessive runoff and soil erosion that can occur when Para rubber trees are grown on high slopes >15%.
The distribution of Para rubber production is also shown in the 2002 baseline land use map, versus the 2009 and 2015 land use scenario maps, in Figure6.These distributions of Para rubber production areas further underscore the dramatic expansion of Para rubber tree plantations that occurred during the 14-year period of 2002 to 2015 in the NLRB.The effects of the increase in Para rubber production between the 2002 baseline and the two scenario years of 2009 and 2015 were accounted for in three separate scenario simulations performed in SWAT.The baseline scenario was first executed using the 2002 land use distribution (Table3) for a 25-year period (1985 to 2009).The 2009 and 2015 land use scenarios were then performed for the same 25-year period to provide a consistent basis of comparison versus the baseline scenario.Water 2016, 8, 1 11 of 19The distribution of Para rubber production is also shown in the 2002 baseline land use map, versus the 2009 and 2015 land use scenario maps, in Figure6.These distributions of Para rubber production areas further underscore the dramatic expansion of Para rubber tree plantations that occurred during the 14-year period of 2002 to 2015 in the NLRB.The effects of the increase in Para rubber production between the 2002 baseline and the two scenario years of 2009 and 2015 were accounted for in three separate scenario simulations performed in SWAT.The baseline scenario was first executed using the 2002 land use distribution (Table3) for a 25-year period (1985 to 2009).The 2009 and 2015 land use scenarios were then performed for the same 25-year period to provide a consistent basis of comparison versus the baseline scenario.

Figure 6 .
Figure 6.Spatial distribution of Para rubber production for the 2002 baseline versus the 2009 and 2015 land use scenarios in the Nam Loei River Basin (NLRB).

Figure 6 .
Figure 6.Spatial distribution of Para rubber production for the 2002 baseline versus the 2009 and 2015 land use scenarios in the Nam Loei River Basin (NLRB).

3. 4 .
Seasonal ET and Water Yield Responses The 25-year average monthly precipitation, ET levels and water yields simulated for the 2002, 2009 and 2015 NLRB land use scenarios during the 25-year simulation period (1985 to 2009) are shown in Figure 8.The initial increase in Para rubber plantations and other land use changes that had occurred by 2009 (Table
-year 2002, 2009 and 2015 land use scenario simulations.The results of the scenarios reflect the increases in Para rubber production and other shifts in land use that occurred during the 2002 to 2009 and 2002 to 2015 time periods.Transmission losses were essentially negligible for all three land use scenarios, and the estimated combined lateral subsurface flow and groundwater flow were very similar between the three scenario simulations.However, the predicted ET increased by nearly 17 mm, and the predicted surface runoff and water yield decreased by similar amounts for the 2009 land use scenario as compared to the 2002 land use scenario.Similar, greater decreases in reduced surface runoff and water yield, relative to the 2009 land use scenario, were estimated for the 2015 land use scenario.However, the simulated ET decreased by almost 12 mm between the 2015 and 2009 land use scenarios.These results underscore the impacts of both the increased Para rubber production and shifting overall land use mixes in the NLRB between 2002 and 2015 (Table

3. 4 .
Seasonal ET and Water Yield Responses The 25-year average monthly precipitation, ET levels and water yields simulated for the 2002, 2009 and 2015 NLRB land use scenarios during the 25-year simulation period (1985 to 2009) are shown in Figure 8.The initial increase in Para rubber plantations and other land use changes that had occurred by 2009 (Table (Table6), with the highest percentage increases occurring during the dry season months of November to January and the lowest percentage increases occurring during the wet season months of May, September and October.The estimated water yield responses between the baseline and 2009 Para rubber expansion scenario resulted in the opposite trend, with water yields decreasing in most months, although slight increases occurred during the months of November, December and January.The predicted percentage water yield changes ranged from +0.6% in January to over −10% in March and April.The percentage declines in water yield predicted for the wet season were more constant as compared to the dry season water yield impacts, ranging from −1% in October to −7% in October.The continued expansion of Para rubber production and other land use changes between 2009 and 2015 (Table

Figure 8 .
Figure 8.Average monthly precipitation, and average monthly water yield and evapotranspiration (ET) for the 2002 baseline, 2009 land use scenario and 2015 land use scenario, for the 25-year (1985 to 2009) simulation period for the Nam Loei River Basin (NLRB).
Table 7 provides comparisons of:(1) the long-term (25-year simulation period) average cumulative amounts of ET (mm) and water yield (mm) that occurred during the wet season, dry season, and annually; (2) the percentages of the cumulative amounts that occurred during the wet season and dry season; and (3) the percentage changes in annual average ET and WYLD that occurred due to the land use changes between 2002 and 2009, and 2009 and 2015.These results show that nearly 80% of the water yield occurred during the wet season during both the 2002 baseline and the 2009 and 2015 Para rubber expansion land use scenarios (Table

Figure 8 .
Figure 8.Average monthly precipitation, and average monthly water yield and evapotranspiration (ET) for the 2002 baseline, 2009 land use scenario and 2015 land use scenario, for the 25-year (1985 to 2009) simulation period for the Nam Loei River Basin (NLRB).

Table 1 .
Model input data sources for Nam Loei River Basin (NLRB).

Table 1 .
Model input data sources for Nam Loei River Basin (NLRB).

Table 2 .
Key Para rubber Tree (RUBR) crop parameters, and two other vegetation-influenced input parameters, that were used for the Nam Loei River Basin (NLRB) SWAT analysis.

Table 3 .
Distribution of land use for the 2002 and 2009 Nam Loei River Basin (NLRB) land use scenarios.

Table 4 .
Parameters ranges and results of the sensitivity analysis at the gauge stations located within the Nam Loei River Basin (NLRB).

Table 5 .
Calibration and validation results at streamflow gauge stations Kh.28A and Kh.58A a in the Nam Loei River Basin (NLRB).
Notes: a Locations of streamflow gauge stations shown in Figure1.

Table 6 .
Long-term (1985 to 2009) average annual water balance components for the 2002 and 2009 land use scenarios as estimated by SWAT for the entire Nam Loei River Basin (NLRB).
Notes: a Total water yield = surface runoff + groundwater flow + lateral flow − transmission loss.

Table 6 .
Long-term (1985 to 2009) average annual water balance components for the 2002 and 2009 land use scenarios as estimated by SWAT for the entire Nam Loei River Basin (NLRB).
Notes: a Total water yield = surface runoff + groundwater flow + lateral flow − transmission loss.

Table 7 .
Comparisons of 25-year cumulative average seasonal and annual ET and water yields, and total annual ET and water yield percentage changes that occurred due to the land use changes between 2002 and 2009, and 2009 and 2015, in the Nam Loei River Basin (NLRB).