Modeling Spatial Soil Water Dynamics in a Tropical Floodplain , East Africa

Analyzing the spatial and temporal distribution of soil moisture is critical for ecohydrological processes and for sustainable water management studies in wetlands. The characterization of soil moisture dynamics and its influencing factors in agriculturally used wetlands pose a challenge in data-scarce regions such as East Africa. High resolution and good-quality time series soil moisture data are rarely available and gaps are frequent due to measurement constraints and device malfunctioning. Soil water models that integrate meteorological conditions and soil water storage may significantly overcome limitations due to data gaps at a point scale. The purpose of this study was to evaluate if the Hydrus-1D model would adequately simulate soil water dynamics at different hydrological zones of a tropical floodplain in Tanzania, to determine controlling factors for wet and dry periods and to assess soil water availability. The zones of the Kilombero floodplain were segmented as riparian, middle, and fringe along a defined transect. The model was satisfactorily calibrated (coefficient of determination; R2 = 0.54–0.92, root mean square error; RMSE = 0.02–0.11) on a plot scale using measured soil moisture content at soil depths of 10, 20, 30, and 40 cm. Satisfying statistical measures (R2 = 0.36–0.89, RMSE = 0.03–0.13) were obtained when calibrations for one plot were validated with measured soil moisture for another plot within the same hydrological zone. Results show the transferability of the calibrated Hydrus-1D model to predict soil moisture for other plots with similar hydrological conditions. Soil water storage increased towards the riparian zone, at 262.8 mm/a while actual evapotranspiration was highest (1043.9 mm/a) at the fringe. Overbank flow, precipitation, and groundwater control soil moisture dynamics at the riparian and middle zone, while at the fringe zone, rainfall and lateral flow from mountains control soil moisture during the long rainy seasons. In the dry and short rainy seasons, rainfall, soil properties, and atmospheric demands control soil moisture dynamics at the riparian and middle zone. In addition to these factors, depths to groundwater level control soil moisture variability at the fringe zone. Our results support a better understanding of groundwater-soil water interaction, and provide references for wetland conservation and sustainable agricultural water management.


Introduction
Analyzing the spatial and temporal distribution of soil moisture content is essential for understanding ecohydrological wetland processes and for the development of sustainable agricultural water management wetland studies [1,2].Soil moisture content is an important variable in regulating and predicting a range of hydrological processes like flooding [3], plant water availability [4][5][6], and soil water-groundwater interaction [7].Soil moisture content is equally crucial in the hydrological cycle for the estimation of the plot-catchment scale water balance [8].Therefore, a better spatial and temporal assessment of soil moisture content yields valuable insights into hydrological processes and supports weather forecasting as well as climate projections [1,9,10].
Wetlands rely on water, whether derived from precipitation, surface, or sub-surface water, to fulfill their ecosystem functions such as nutrient retention, flood mitigation, and food provisioning through agricultural support [11,12].Indeed, knowledge of hydrology is indispensable for wetland conservation and agricultural production [13][14][15].In East Africa, floodplain wetlands are increasingly being utilized for agricultural production development projects, as a key intervention to synergistically achieve food security in the region [16].This is due to their relatively large size, high soil nutrient stocks, and prolonged periods of soil water availability [17][18][19].This is also evident for the floodplain wetland in the catchment of Kilombero river, the most important tributary of the Rufiji river basin, representing one of the largest basins of Tanzania.The catchment is, however, increasingly used for agricultural production, especially during periods of inundation in the long rainy season.According to References [20,21], rain-fed rice production has recently increased in the Kilombero floodplain.In fact, the wetland supplies about 9% of all rice produced in Tanzania.It is noteworthy that the Kilombero river catchment is also targeted for agricultural expansion through the project Southern Agricultural Growth Corridor Of Tanzania (SAGCOT) [22], a large-scale agricultural intensification project.In this context, assessment and an improved understanding of the amount and spatial dynamics of soil moisture is fundamental for crop growth, land use planning, and the evaluation of the potential impacts of agricultural management interventions on the hydrology and other ecosystem services of the floodplain.
Spatial and temporal soil moisture dynamics in floodplain wetlands are influenced by a range of factors; thus, assessing these factors through field observations is difficult and time-consuming.Soil moisture exhibits tremendous spatial heterogeneity even over small catchments, resulting from the interplay of hydrological, biological, and meteorological processes [23].In order to better understand the spatial pattern of soil moisture for wetland conservation and land use planning, good-quality time series data is required.However, there is a paucity of reliable field-scale measured data, and the remotely sensed soil moisture data is too coarse in the spatial resolution and refers to 5-cm depth soil moisture [24].Often, data gaps for the available soil moisture time series exist due to measurement constraints, device malfunctioning, or site inaccessibility during flooding periods.Moreover, during dry conditions, soil cracking affects the measuring moisture sensors and, therefore, data quality [25].Statistical and data-driven methods are an option for filling in missing values in distributed soil moisture datasets.Alternatively, models that integrate meteorological conditions, soil water storage, land use, and depth to groundwater table [26] may significantly overcome data gaps at a point scale, which are inevitable under field conditions.Once these models are well calibrated and validated, assessment of the impacts of agricultural management on the ecosystem functioning of wetlands is feasible.
In this study, a widely used soil water model, Hydrus-1D [27], which simulates one-dimensional water, heat, and solute transport in variably saturated and unsaturated media, was applied to simulate the spatial soil moisture dynamics in the Kilombero floodplain.Models such as MODular Hydrologic Modelling System (MODHMS) [28], European Hydrological System (MIKE SHE) [29], and Hydrus-2D and -3D [30] can be adopted to simulate soil water dynamics.However, they are data demanding, require large computer resources, and are difficult to adapt to the investigated study area because of complex boundary conditions.Hydrus-1D, in comparison, is less data demanding and can easily be adopted at the local scale as well as at the East African regional scale, and has also been proven to be an effective tool for evaluating various water and solute fluxes in agricultural fields with different crops and various irrigation schemes [31,32].Hydrus-1D has been successfully applied in numerous studies for predicting soil moisture content and movement under different conditions [33][34][35][36].In spite of these previous efforts, there is still limited research on the assessment of the spatial soil moisture dynamics at the different hydrological zones of a tropical wetland using Hydrus-1D, especially for data-scarce tropical regions such as East Africa.Most of these studies [31,[34][35][36] have been conducted in temperate regions, yet tropical regions continue to face data scarcity and quality challenges, especially in the recently degrading wetland ecosystems.The similarity of this study to previous works is the type of data (soil moisture content, depth to groundwater level, soil physical properties, leaf area index, and climatic data) used for model calibration and simulation of soil water dynamic and vertical fluxes.Many studies have described and explained the factors that influence soil water dynamics [6, 37,38] and soil moisture retention [16,39,40] at the wetland-catchment scale in East Africa.However, these studies are tailored to inland valley wetlands, yet floodplain wetlands exhibit different hydrological behavior compared to the inland valley wetlands.Therefore, there is paucity of sufficient scientific information for proper floodplain wetland agricultural management planning in the region.This calls for more research on these floodplains to improve the scientific knowledge base for decision-making.The purpose of the study is to conduct a thorough study of the soil water dynamic based on field measurements so as to evaluate the agricultural usability of the different hydrological zones of the tropical floodplain wetland.This contributes to soil and land management as it provides the required scientific basis for decision-making and agricultural planning.Accordingly, the study has the following objectives (i) to test the applicability of the Hydrus-1D model for simulating spatial soil moisture content dynamics along different hydrological zones of the floodplain; (ii) to determine factors controlling soil moisture content dynamics at each hydrological zone during drying and wetting periods; and (iii) to assess soil water availability for each hydrological zone.In addition, a sensitivity analysis on selected water flux parameters, with varying depth to groundwater levels, was conducted for each hydrological zone.This was performed to understand how changes in groundwater levels, which may be caused by water extraction and catchment scale processes, will influence local water availability and water fluxes.

Study Area
The Kilombero floodplain wetland is part of the Rufiji basin and is situated in Kilombero and Ulanga districts, Morogoro region, southern Tanzania.The floodplain lies between latitudes 7.654 • and 10.023 • S and longitudes 34.563 • and 37.797 • E (Figure 1a), with elevation ranging from 250 to 300 m above sea level (a.s.l).The Kilombero catchment covers an area of 40,240 km 2 , including the broad floodplain of the main Kilombero river with a spatial extent of about 7967 km 2 [41], in which the study area of about 96 km 2 is found (Figure 1).The floodplain is seasonally flooded and drained by the Kilombero river.The Kilombero river is an important tributary catchment of the Rufiji river basin, contributing 62% of the total annual average basin water flow [42].A number of perennial and seasonal streams drain into the Kilombero River from the Udzungwa Mountains in the northwest and the Mahenge Mbarika escarpment in the southeast with elevation ranging from 200 to 2500 m a.s.l.[41,43].
The Kilombero floodplain is part of the Southern Agricultural Growth Corridor Of Tanzania (SAGCOT) [22], an area earmarked for future investments in agricultural development.The area forms a large triangle of some 500 km 2 that is bordered by three national parks.Besides a range of agro-industrial uses (mainly sugar cane), the cropped portion is dominated by rain-fed rice (wet season) and several irrigation schemes allowing for year-round irrigated rice production.In fact, the wetland supplies about 9% of all rice produced in Tanzania compared to the Mkomazi floodplain in the sub-humid lowlands west of Usambara Mountains.The area is characterized by diverse land uses, land use intensity gradients, and interactions between large-and small-scale crop farmers, landless herders, and urban populations.
The climate is sub-humid tropical with annual rainfall ranging from 1200 to 1400 mm and a mean annual temperature of about 23-25 • C [44].The rainfall is distributed in a bi-modal structure.The long rainy season is from March to May, while the short rainy season spans from October to December [41].The flood peaks are experienced mainly during the long rainy season, though they can also happen as early as January.
The geological substrate in the floodplain includes mainly a sedimentary basin, forming a seasonal alluvial floodplain dominated by Fluvisols soils [45].However, according to Reference [46], sediments constitute alluvial fans, river deposits, and colluvial materials.The alluvial fans mark the transition between the floodplain and the fringes, with high infiltration rates and high soil moisture variability [46].The soils tend to crack upon drying and experience ponding once the soil is swollen and cracks are closed.Within the study area, groundwater flow is directed from north to south towards Kilombero river throughout the year [46].Groundwater recharge mainly occurs due to rainwater infiltration in the study area and the mountainous areas in the northern catchment.
The valley is predominantly characterized by Miombo woodlands, mainly Brachystegia species and grasses such as reed (Phragmites mauritianus), guinea grass (Panicum maximum), and elephant grass (Penisetum purpureum) [21].Agriculture is the main economic activity in the valley.Agricultural production is practiced at both subsistence and commercial scales.Lowland rain-fed rice is the dominating food and cash crop cultivated in the valley; however, other crops such as maize and peas are grown, especially during the short rainy season.
The climate is sub-humid tropical with annual rainfall ranging from 1200 to 1400 mm and a mean annual temperature of about 23-25 °C [44].The rainfall is distributed in a bi-modal structure.The long rainy season is from March to May, while the short rainy season spans from October to December [41].The flood peaks are experienced mainly during the long rainy season, though they can also happen as early as January.
The geological substrate in the floodplain includes mainly a sedimentary basin, forming a seasonal alluvial floodplain dominated by Fluvisols soils [45].However, according to Reference [46], sediments constitute alluvial fans, river deposits, and colluvial materials.The alluvial fans mark the transition between the floodplain and the fringes, with high infiltration rates and high soil moisture variability [46].The soils tend to crack upon drying and experience ponding once the soil is swollen and cracks are closed.Within the study area, groundwater flow is directed from north to south towards Kilombero river throughout the year [46].Groundwater recharge mainly occurs due to rainwater infiltration in the study area and the mountainous areas in the northern catchment.
The valley is predominantly characterized by Miombo woodlands, mainly Brachystegia species and grasses such as reed (Phragmites mauritianus), guinea grass (Panicum maximum), and elephant grass (Penisetum purpureum) [21].Agriculture is the main economic activity in the valley.Agricultural production is practiced at both subsistence and commercial scales.Lowland rain-fed rice is the dominating food and cash crop cultivated in the valley; however, other crops such as maize and peas are grown, especially during the short rainy season.

Study Design
The valley was subdivided into three hydrological zones (riparian, middle, and fringe, see Figure 1a) based on the origin of flooding water, its extent, and duration.In total, 462 Landsat 7 Enhanced Thematic Mapper Plus (ETM+) and eight Operational Land Imager (OLI) surface reflectance images [47,48] across six tiles from the years 2013 to 2015 with a cloud coverage <80% were processed.Their multi-temporal statistical values were classified in a Random Forest approach to delineate riverine flooded areas and generate inundation extents for seasonally flooded zones during the long rainy season (Figure 1a).The non-flooded zone was defined as fringe, characterized by precipitationand groundwater table-induced flooding [49].The middle zone is partially flooded, whereas the riparian zone is nearly completely flooded during the long rainy season.Further spatial differentiation of the hydrological zones included the application of the TerraSAR-X digital elevation model with a resolution of 12 m (TanDEM-X-12 m) to derive the absolute height and the distance from the river [50].The lower parts of the transect are flooded longer and with a higher water depth, thus the riparian zone was close to the river with a lower absolute height altitude (Figure 1a).

Hydrological Data Monitoring
Hydrological data variables were monitored for a period of 17 months (March 2015-June 2016).

Meteorological Data
Daily rainfall values were measured in a tipping bucket rain gauge with a 0.2-mm resolution [51], installed at the fringe zone of the floodplain (Figure 1).Temperature, relative humidity, wind speed, and solar radiation were obtained from an automatic weather station situated close to the study area (Figure 1a).

Soil Moisture and Groundwater Level Monitoring
Soil moisture measurements were conducted at the three defined hydrological zones.At each hydrological zone, at least four soil moisture access tubes (Figure 1b) were installed for a Frequency-Domain-Reflectometry (FDR) profile probe type PR2 [52] to measure soil moisture content every two days.Other moisture measurement devices such as soil Hydra probe sensors [51], decagon 5TE soil moisture sensors [53], Theta probes, and tensiometers can be adopted for field soil moisture measurements.In this study, the profile probe used was a portable device which takes simultaneous measurements at depths of 10, 20, 30, and 40 cm within the access tube.The FDR profile probe type PR2 was used to measure the profile soil moisture up to soil depths of 40 cm and 100 cm, consisting of a sealed polycarbonate rod with a diameter of about 25 mm, electronic sensors arranged at fixed intervals along its length, and presenting a measurement accuracy of ±0.06 m 3 /m 3 at 0 to 40 • C under generalized soil calibration in normal soils [52].FDR is based on a unique relation between the volumetric water content and the dielectric constant of mineral soils.According to References [25,54,55], a Delta-T profile probe shows measurement difficulties in high clay soils; however, results are improved with site-specific calibration.In this study, site-specific calibrations were not conducted since the calibration of FDR sensors to tropical wetland soil conditions have been established by Reference [37], who approved their applicability with an error of 0.07 m 3 /m 3 .Therefore, the measured moisture results from the probes are reliable for tropical wetland conditions.Soil moisture readings were manually recorded with a handheld (HH2) moisture meter connected to a profile probe sensor [56].
To monitor the depth to groundwater level in the study area, three 10-cm-diameter piezometer pipes were installed, one at each hydrological zone (Figure 1b).The depth to groundwater level was measured hourly using pressure data loggers [57] installed in the piezometers.

Soil Properties Sampling
Soil physical properties were characterized at each installed access tube.Composite soil samples were collected at depths of 10, 20, 30, and 40 cm using a soil auger to analyze the soil particle size distribution and soil organic carbon (SOC).Soil particle size distribution was measured using the laser method [58], while SOC was determined using the modified Blackley Walkley wet method [59].Three undisturbed soil samples (250 cm 3 ) for each sampled soil depth were taken at each site for saturated hydraulic conductivity (K sat ) measurements using the constant head method with a laboratory permeameter [60].Bulk density was determined with the core method using the same samples as those employed for the saturated hydraulic conductivity.These soil properties were used to drive the model setup and calibration.

Model Description
Simulations of soil water dynamics were performed using the Windows-based one-dimensional Hydrus-1D v4.09 software package [27] developed by PC-Progress s.r.o in Czech Republic.The model simulates variably saturated-unsaturated water flow by numerically solving Richard's equation [61]: where θ is the volumetric soil water content (L 3 /L 3 ), t is time (T), x is the vertical space coordinate (L), K is the hydraulic conductivity, h is the water pressure head (L), and S (h) is a water sink term accounting for root water uptake (L 3 /(L 3 •T)).Unsaturated soil hydraulic properties were represented by the parameterization given by [62]: where θ r and θ s are residual and saturated water contents (L 3 /L 3 ), respectively, K s is the saturated hydraulic conductivity (L/T 1 ), S e is the effective saturation, α (L −1 ), and l and n represent empirical shape parameters (the inverse air-entry point, pore connectivity, and pore-size distribution parameter, respectively), m = 1 − 1/n.The sink term S (h) in Equation ( 1) is defined as the volume of water derived from a unit volume of soil per unit time by plant roots.It accounts for actual root water uptake equivalent to actual transpiration, calculated by the model using Feddes equation [63]: where α (h) is the root water uptake stress response function, and S p is a potential root water uptake rate (T −1 ).Root water uptake is assumed to be zero close to saturation (i.e., wetter than the "anaerobiosis point", h 1 ), and dryer than the wilting point pressure head (h 4 ).Root water uptake is considered to be at optimum between pressure heads h 2 and h 3 , while root water uptake increases (or decreases) linearly with h when the pressure heads are between h 3 and the wilting point (or between h 2 ).Potential root water uptake is equal to root water uptake rate during periods of no water stress (when α (h) = 1).

Initial and Boundary Conditions
The initial conditions in the model were set to pressure heads in equilibrium with the measured depth to groundwater level, which was taken as a lower boundary condition.When the water table was above the soil surface, the pressure head was fixed to full saturation to facilitate modeling.The upper boundary condition was specified as an atmospheric boundary condition with a surface layer (to permit water build-up due to flooding) using daily rainfall, potential evapotranspiration, and leaf area index (LAI).Potential evapotranspiration (ET p ) was calculated by the Food and Agriculture Organization (FAO) Penman-Monteith equation [64].A number of methods, for example, the water balance method [65], the Priestley-Taylor method [66], the Bowen-ratio method [67], and the Hargreaves equation [68], can be used for the calculation of potential evapotranspiration.According to Reference [64], the FAO Penman-Monteith method is physically based and can be used to predict ET p in a wide range of locations and climates, even in data-scarce situations like Sub-Saharan Africa.The method uses meteorological data such as maximum and minimum temperature, relative humidity, wind speed, and solar radiation.These data were obtained from the automatic weather station located within the study site.Time series for meteorological conditions used for the simulation period are depicted in Figure 2. Actual evaporation and transpiration were directly computed by the model based on the given soil moisture conditions and root water uptake functions.
Water 2018, 10, x FOR PEER REVIEW 7 of 27

Initial and Boundary Conditions
The initial conditions in the model were set to pressure heads in equilibrium with the measured depth to groundwater level, which was taken as a lower boundary condition.When the water table was above the soil surface, the pressure head was fixed to full saturation to facilitate modeling.The upper boundary condition was specified as an atmospheric boundary condition with a surface layer (to permit water build-up due to flooding) using daily rainfall, potential evapotranspiration, and leaf area index (LAI).Potential evapotranspiration (ETp) was calculated by the Food and Agriculture Organization (FAO) Penman-Monteith equation [64].A number of methods, for example, the water balance method [65], the Priestley-Taylor method [66], the Bowen-ratio method [67], and the Hargreaves equation [68], can be used for the calculation of potential evapotranspiration.According to Reference [64], the FAO Penman-Monteith method is physically based and can be used to predict ETp in a wide range of locations and climates, even in data-scarce situations like Sub-Saharan Africa.The method uses meteorological data such as maximum and minimum temperature, relative humidity, wind speed, and solar radiation.These data were obtained from the automatic weather station located within the study site.Time series for meteorological conditions used for the simulation period are depicted in Figure 2. Actual evaporation and transpiration were directly computed by the model based on the given soil moisture conditions and root water uptake functions.Leaf area index (LAI) values were measured at different growth stages of rain-fed lowland rice using a LI-COR area meter, an indirect method.Direct methods like the planimetrical approach and the gravimetric procedure can be used, but are time-consuming and destructive [69].The LI-COR area meter is fast and uses a fisheye lens to project a nearly hemispheric image of the canopy and sky onto a ringed detector.This enables quick measurement of the leaf area index and foliage orientation over a large area [69].In this study, LAI values were linearly interpolated between the measurement dates to obtain time series values.

Model Calibration and Validation
The model was calibrated for one year (March 2015 to March 2016), by comparing the simulated and field-measured daily soil moisture content at the different soil depths (10,20,30, and 40 cm), and different plots for each hydrological zone.Plots PR_1, PR_3, PR_4, and PR_8 for the riparian zone, PM_3, PM_5, and PM_7 for the middle zone, and PF_4, PF_8, PF_10, and PF_12 for fringe zone were used for calibration (Figure 1b).The one-year period included the drying and wetting conditions, allowing us to fit the model at a wide range of soil moisture conditions in the floodplain.The soil material layer/horizon was defined in accordance with depth intervals used for monitoring soil moisture content in a profile of 200 cm.A profile of 200 cm was defined because the observed groundwater table is within this range.
Hydrus-1D was calibrated using site-specific boundary conditions and field-measured soil moisture content.Van Genuchten hydraulic parameters, which describe the soil retention and unsaturated hydraulic conductivity functions, were optimized through inverse modeling with no Leaf area index (LAI) values were measured at different growth stages of rain-fed lowland rice using a LI-COR area meter, an indirect method.Direct methods like the planimetrical approach and the gravimetric procedure can be used, but are time-consuming and destructive [69].The LI-COR area meter is fast and uses a fisheye lens to project a nearly hemispheric image of the canopy and sky onto a ringed detector.This enables quick measurement of the leaf area index and foliage orientation over a large area [69].In this study, LAI values were linearly interpolated between the measurement dates to obtain time series values.

Model Calibration and Validation
The model was calibrated for one year (March 2015 to March 2016), by comparing the simulated and field-measured daily soil moisture content at the different soil depths (10,20,30, and 40 cm), and different plots for each hydrological zone.Plots PR_1, PR_3, PR_4, and PR_8 for the riparian zone, PM_3, PM_5, and PM_7 for the middle zone, and PF_4, PF_8, PF_10, and PF_12 for fringe zone were used for calibration (Figure 1b).The one-year period included the drying and wetting conditions, allowing us to fit the model at a wide range of soil moisture conditions in the floodplain.The soil material layer/horizon was defined in accordance with depth intervals used for monitoring soil moisture content in a profile of 200 cm.A profile of 200 cm was defined because the observed groundwater table is within this range.
Hydrus-1D was calibrated using site-specific boundary conditions and field-measured soil moisture content.Van Genuchten hydraulic parameters, which describe the soil retention and unsaturated hydraulic conductivity functions, were optimized through inverse modeling with no hysteresis.In the Hydrus-1D model, inverse parameter estimation employs a relatively simple, gradient-based, local optimization algorithm based on the Marquardt-Levenberg method [27].The aim of this method is to determine the best estimate of the hydraulic model parameters [61].Hydraulic parameters estimated via the Rosetta pedotransfer functions [70] using measured soil particle size distribution (percentage clay, silt, and sand) and bulk density were used as initial estimates, in addition to the measured saturated hydraulic conductivity (Table 1).Alpha, n, and saturated water content parameters were fitted first, since the Hydrus-1D model could optimize only 15 parameters at a time.Residual water content and the L parameter were the last to be estimated.The values of K sat were reviewed and re-adjusted accordingly.During the rice growing period, Feddes' parameters for the tropical climate were taken from [71]: h 1 = 100 cm, h 2 = 55 cm, h 3 = −160 cm, h 4 = −250 cm, h 4 = −16,000 cm and the root distribution with a maximum of 40 cm were specified according to Reference [72].After rice harvest (i.e., during the dry and the short rainy seasons, in which the sites were under fallow), we used the Feddes' parameters for grass implemented in Hydrus-1D.
The calibrated model from each plot was validated for one year, with measured soil moisture data for the neighboring plots at each hydrological zone.The aim of this approach was to test the transferability of the calibrated model to predict spatial soil moisture for sites with similar biophysical properties and climate conditions.At each hydrological zone, two calibrated plots were validated with measured moisture data of the neighboring plots.At the riparian zone, plot PR_1 and PR_4 were validated with measured soil moisture data from plot PR_3 (88 m distant) and PR_1 (75 m distant), respectively.At the middle zone, plot PM_3 and PM_7 were validated with measured soil moisture values from plot PM_7 (400 m distant) and plot PM_5 (102 m distant), respectively.At the fringe, plot PF_10 and PF_8 were validated with measured soil moisture data from plot PF_12 (109 m distant), and PF_4 (225 m distant), respectively (Figure 1b).

Criterion for Model Calibration Performance
The agreement between predicted and observed moisture data was evaluated by (i) the coefficient of determination (R 2 ); (ii) the Nash-Sutcliffe model efficiency coefficient (NSE); and (iii) the root mean square error (RMSE).Kling-Gupta efficiency (KGE) [73] and percent bias (PBIAS) [74] are other evaluation approaches which can be used for evaluating model quality.In this study, we used R 2 , NSE, and RMSE.R 2 describes the proportion of the variance in measured data explained by the model.R 2 ranges from 0 to 1, with values greater than 0.5 considered to be acceptable [75].The NSE is a dimensionless quality measure used to determine the relative magnitude of the residual variance between the simulated and measured data variance [76].NSE ranges between −∞ and 1.0.An NSE of 1.0 indicates a perfect fit between the simulated and observed data [75].RMSE is an error index statistic which measures the difference between simulated and observed values.The lower the RMSE, the better the performance of the model.
where O i and P i are the observed and predicted moisture values, respectively, O and P are the mean of observed and predicted moisture values, and n is the number of observations.

Sensitivity Analysis
A sensitivity analysis was performed to analyze the impact of changes in the depth to groundwater level (dgwl) on root water uptake, evaporation, and change in soil water storage parameters.Evapotranspiration (transpiration and evaporation) is the major cause of water loss in the wetlands [36,77], while soil water storage is the most important state variable for crop growth.Therefore, understanding the impact of changes in the depth to groundwater level is essential to assess the water balance of the study area for sustainable agricultural and wetland management planning.Considering data for one year (July 2015-June 2016) as a reference, the mean dgwl was increased and decreased by 15%, 30%, 45%, and 60%.To evaluate the sensitivity of each parameter, the sensitivity index (SI) for each percentage change in dgwl was calculated according to de Roo (1993) cited in [78]: With P i = model output with a 15%, 30%, 45%, and 60% increase in the mean dgwl; P m = model output with a related 15%, 30%, 45%, and 60% decrease in the mean dgwl; P = model output with base dgwl.
A number of sensitivity analysis methods exist, for example, variance-based methods, correlation and regression analysis methods, and Monte Carlo methods [79].We used the sensitivity index method for this study because of its simplicity and because the sensitivity is clearly attributed to one model parameter.

Simulation of Spatial Soil Water Dynamics at Different Scales
Hyrus-1D was set up, calibrated, and validated at the riparian, middle, and fringe hydrological zones of the floodplain.

Calibration Results
A comparison of model performances and quality of statistical measures are presented in Table 2 for the riparian zone.Time series soil moisture content was well captured by the model for the four calibrated individual plots.A good agreement between modeled and observed soil moisture content as indicated by high R 2 (0.67-0.92),NSE (0.63-0.91), and low RMSE (0.03-0.06) was obtained (Table 2).
The best fit of van Genuchten hydraulic parameters for each plot is shown in Table 3.The calibrated hydraulic parameters showed slight differences among the plots at soil depths of 10, 20, 30, and 40 cm.The final calibrated saturated hydraulic conductivity was higher than the measured values (Tables 1 and 3).Soil moisture content dynamics for both the rainy and dry seasons were well reproduced by the model for all of the soil layers.Measured and modeled soil moisture content values increased with soil depth.Flooding from both river overbank flow and local rainfall occurred during the long rainy season (March to May) and in the months of January to February 2016, which experienced heavy rains, correlated to El Niño-Southern Oscillation (ENSO) [80,81] and the Indian Ocean zonal mode (IOZM) [82].
At the riparian zone, the depth to groundwater level increased on average up to 1.6 m above the soil surface during the long rainy seasons; thus, the soils were fully saturated.During these periods, the model was able to reproduce soil moisture content trends and showed constant values (indicating soil saturation) until the floods receded (Figure 3).During flood recession (when the depth to groundwater level increased below the surface i.e., after June 2015), both modeled and measured soil moisture content showed a similar decreasing trend.For the short rainy season, the depth to groundwater level continued to increase below the surface until it reached 150 cm.Interestingly, measured and modeled soil moisture content at the different soil depths responded positively to rainfall events.Therefore, rainfall controlled soil moisture dynamics during this period (Figure 3).In addition, soil moisture dynamics could be explained by the differences observed in calibrated hydraulic parameters along the soil depths.The calibrated soil hydraulic parameters showed slight differences among the plots by visual comparison (Table 3).Thus, within the hydrological zones, there were relatively similar soil moisture patterns among the plots.This was depicted in the measured and modeled soil moisture content for the individual plots in the riparian zone, which showed almost similar trends throughout the different seasons (Figure 3).At the riparian zone, the depth to groundwater level increased on average up to 1.6 m above the soil surface during the long rainy seasons; thus, the soils were fully saturated.During these periods, the model was able to reproduce soil moisture content trends and showed constant values (indicating soil saturation) until the floods receded (Figure 3).During flood recession (when the depth to groundwater level increased below the surface i.e., after June 2015), both modeled and measured soil moisture content showed a similar decreasing trend.For the short rainy season, the depth to groundwater level continued to increase below the surface until it reached 150 cm.Interestingly, measured and modeled soil moisture content at the different soil depths responded positively to rainfall events.Therefore, rainfall controlled soil moisture dynamics during this period (Figure 3).In addition, soil moisture dynamics could be explained by the differences observed in calibrated hydraulic parameters along the soil depths.The calibrated soil hydraulic parameters showed slight differences among the plots by visual comparison (Table 3).Thus, within the hydrological zones, there were relatively similar soil moisture patterns among the plots.This was depicted in the measured and modeled soil moisture content for the individual plots in the riparian zone, which showed almost similar trends throughout the different seasons (Figure 3).Overall, there was a good agreement between modeled and measured soil moisture content after calibration (R 2 = 0.59-0.89,NSE = 0.51-0.88,and RMSE = 0.02-0.10cm 3 cm 3 ) (Table 4).The best fit of the van Genuchten hydraulic parameters for individual plots are shown in Table 5.There was almost no difference observed between the calibrated hydraulic parameters among the plots, except for saturated hydraulic conductivity.Hence, small variations between the modeled and measured soil moisture content values were noted among the plots.Overall, there was a good agreement between modeled and measured soil moisture content after calibration (R 2 = 0.59-0.89,NSE = 0.51-0.88,and RMSE = 0.02-0.10cm 3 /cm 3 ) (Table 4).The best fit of the van Genuchten hydraulic parameters for individual plots are shown in Table 5.There was almost no difference observed between the calibrated hydraulic parameters among the plots, except for saturated hydraulic conductivity.Hence, small variations between the modeled and measured soil moisture content values were noted among the plots.The depth to groundwater level at the middle hydrological zone increased up to approximately 100 cm above the soil surface, causing full saturation of the soil layers during the long rainy season (March-June) (Figure 4).In the dry and short rainy seasons, the depth to groundwater level gradually declined up to 300 cm below the surface.Moreover, the model was capable of simulating the seasonal soil moisture content variations.During the short rainy season, both simulated and measured soil moisture content strongly responded to rainfall events in spite of the continued declining depth to groundwater level (Figure 4).Therefore, rainfall and soil properties played an important role in controlling soil moisture dynamics during dry and short rainy seasons at the middle zone.Modeled soil moisture dynamics during this period were well captured, except for some plots in which the simulated soil moisture showed high amplitude from the measured soil moisture content.Overall, the simulations were in agreement with the measurements.
Figure 5 and Table 6 show the calibrated results of soil moisture and model performance, respectively at the fringe zone of the floodplain.The modeled soil moisture content depicted a trend relatively similar to the measured soil moisture over time at the 10-cm depth intervals among the plots.The R 2 among the calibrated plots varied from 0.54 to 0.87, RMSE ranged from 0.03 to 0.11 cm 3 /cm 3 , and NSE varied from <0 to 0.85 (Table 6).Although modeled values were satisfactory, the model somewhat over-and underestimated soil moisture content at all of the plots (Figure 5).
Water 2018, 10, x FOR PEER REVIEW 13 of 27 The depth to groundwater level at the middle hydrological zone increased up to approximately 100 cm above the soil surface, causing full saturation of the soil layers during the long rainy season (March-June) (Figure 4).In the dry and short rainy seasons, the depth to groundwater level gradually declined up to 300 cm below the surface.Moreover, the model was capable of simulating the seasonal soil moisture content variations.During the short rainy season, both simulated and measured soil moisture content strongly responded to rainfall events in spite of the continued declining depth to groundwater level (Figure 4).Therefore, rainfall and soil properties played an important role in controlling soil moisture dynamics during dry and short rainy seasons at the middle zone.Modeled soil moisture dynamics during this period were well captured, except for some plots in which the simulated soil moisture showed high amplitude from the measured soil moisture content.Overall, the simulations were in agreement with the measurements.
Figure 5 and Table 6 show the calibrated results of soil moisture and model performance, respectively at the fringe zone of the floodplain.The modeled soil moisture content depicted a trend relatively similar to the measured soil moisture over time at the 10-cm depth intervals among the plots.The R 2 among the calibrated plots varied from 0.54 to 0.87, RMSE ranged from 0.03 to 0.11 cm 3 /cm 3 , and NSE varied from <0 to 0.85 (Table 6).Although modeled values were satisfactory, the model somewhat over-and underestimated soil moisture content at all of the plots (Figure 5).Likewise, the calibrated van Genuchten hydraulic parameters showed slight variations at 10-cm depth intervals within and across the individual plots.This could be due to the low heterogeneity in soil texture (clay, sand, and silt content) and bulk density within the fringe position.Saturated hydraulic conductivity greatly varied along the four soil depths within and across the plots (Table 7), which can be explained by the effect of soil macrospores, the cracking, and the occurrence of slickensides as a consequence of swell-and-shrink behavior, which were commonly encountered in clay-rich horizons during sampling.Modeled and measured soil moisture time series showed a positive response to rainfall events, although high dynamics were observed during the wetting and drying periods (October to December 2016) for all of the soil depths across the plots.The monitored plots were flooded for a short time during the rainy season (March to June during the study) and in the short heavy rains (January to February 2016) (Figure 5).The depth to groundwater level gradually increased up to approximately 50 cm above the soil surface, saturating the soil layers.Thereafter, the depth to groundwater level declined to approximately 120 cm below the surface during the dry season (July to September 2015).However, in comparison to the riparian and middle zones, the variations in the depth to groundwater level were lower at the fringe zone.Equally important, the model was able to capture the seasonal soil moisture variability at the different plots amidst the existing data gaps in the measured soil moisture content.Also, the discrepancies between the simulated and measured soil moisture time series at the fringe were not different from those observed at the riparian and middle zones across the calibrated plots.Unlike at the riparian and middle zones, the depth to groundwater level and soil moisture content exhibited a positive response to single rainfall events during the wetting and drying periods (October to December 2016) at the fringe zone.Thus, the depth to groundwater level, rainfall, and soil properties control soil moisture dynamics at this zone.
Additionally, we tested the hypothesis that the model can simulate soil moisture dynamics averaged from individual plots, after aggregating the plot-specific calibrated hydraulic parameters within the hydrological zone.Therefore, the calibrated van Genuchten hydraulic parameters for each plot from the respective hydrological zones were averaged at each soil depth (Table 8).The Hydrus-1D model was set up with these averaged hydraulic parameters to simulate soil moisture content for each hydrological zone.The simulations were evaluated by comparing the modeled soil moisture content with the averaged measured soil moisture content at 10-cm depth increments.Figure 6 shows the typical results.
Water 2018, 10, x FOR PEER REVIEW 15 of 27 periods (October to December 2016) at the fringe zone.Thus, the depth to groundwater level, rainfall, and soil properties control soil moisture dynamics at this zone.Additionally, we tested the hypothesis that the model can simulate soil moisture dynamics averaged from individual plots, after aggregating the plot-specific calibrated hydraulic parameters within the hydrological zone.Therefore, the calibrated van Genuchten hydraulic parameters for each plot from the respective hydrological zones were averaged at each soil depth (Table 8).The Hydrus-1D model was set up with these averaged hydraulic parameters to simulate soil moisture content for each hydrological zone.The simulations were evaluated by comparing the modeled soil moisture content with the averaged measured soil moisture content at 10-cm depth increments.Figure 6 shows the typical results.Overall, the modeled soil moisture time series showed good agreement with the averaged measured soil moisture content.There were some discrepancies between the modeled and measured soil moisture content values.This could be explained by the low variations in the hydraulic parameters observed at each individual plot of the respective hydrological zones.Generally, the model performance was good at the different hydrological zones.For example, R 2 varied from 0.78 to 0.90, NSE = 0.74-0.87,and RMSE ranged from 0.03 to 0.06 at the riparian zone; at the middle zone, Overall, the modeled soil moisture time series showed good agreement with the averaged measured soil moisture content.There were some discrepancies between the modeled and measured soil moisture content values.This could be explained by the low variations in the hydraulic parameters observed at each individual plot of the respective hydrological zones.Generally, the model performance was good at the different hydrological zones.For example, R 2 varied from 0.78 to 0.90, NSE = 0.74-0.87,and RMSE ranged from 0.03 to 0.06 at the riparian zone; at the middle zone, R 2 = 0.62-0.90,NSE = 0.54-0.86,and RMSE = 0.03-0.06;and for the fringe zone, R 2 = 0.75-0.82,RMSE = 0.03-0.07,and NSE ranged from 0.37 to 0.75 (Table 9).Table 9. Hydrus-1D performance for soil moisture simulation after using averaged hydraulic parameters calculated from all of the calibrated plots for each hydrological zone.We randomly validated the calibrated plots against the neighboring plots with measured soil moisture content data at soil depths of 10 and 30 cm for each hydrological zone.In doing so, we were able to test the transferability of model parameters calibrated at one plot to predict spatial soil moisture dynamics of another plot located in the same hydrological zone.As expected, model accuracy was higher for calibration than validation.Nevertheless, the modeled results are reasonable with low variability against the measured soil moisture content at each plot along the different soil depths by visual comparison (Figure 7).Validation resulted in R 2 = 0.64-0.89and RMSE = 0.06-0.07for the riparian zone; R 2 = 0.36-0.82and RMSE = 0.07-0.13for the middle zone; and R 2 = 0.52-0.91 and RMSE = 0.04-0.10for the fringe zone (Table 10).The goodness-of-fit values indicated satisfactory agreement between modeled values and field measurements for both rainy and dry seasons and demonstrated the acceptable accuracy of model simulations.Thus, the results supporting the hypothesis of the model's feasibility to simulate spatial soil water dynamics within sites of similar topography, as well as other biophysical properties once it is well calibrated.

Soil
Water 2018, 10, x FOR PEER REVIEW 16 of 27 R 2 = 0.62-0.90,NSE = 0.54-0.86,and RMSE = 0.03-0.06;and for the fringe zone, R 2 = 0.75-0.82,RMSE = 0.03-0.07,and NSE ranged from 0.37 to 0.75 (Table 9).We randomly validated the calibrated plots against the neighboring plots with measured soil moisture content data at soil depths of 10 and 30 cm for each hydrological zone.In doing so, we were able to test the transferability of model parameters calibrated at one plot to predict spatial soil moisture dynamics of another plot located in the same hydrological zone.As expected, model accuracy was higher for calibration than validation.Nevertheless, the modeled results are reasonable with low variability against the measured soil moisture content at each plot along the different soil depths by visual comparison (Figure 7).Validation resulted in R 2 = 0.64-0.89and RMSE = 0.06-0.07for the riparian zone; R 2 = 0.36-0.82and RMSE = 0.07-0.13for the middle zone; and R 2 = 0.52-0.91 and RMSE = 0.04-0.10for the fringe zone (Table 10).The goodness-of-fit values indicated satisfactory agreement between modeled values and field measurements for both rainy and dry seasons and demonstrated the acceptable accuracy of model simulations.Thus, the results supporting the hypothesis of the model's feasibility to simulate spatial soil water dynamics within sites of similar topography, as well as other biophysical properties once it is well calibrated.

Soil Water Availability at the Different Hydrological Zones
To assess soil water availability, annual soil water balance was established at each hydrological zone for a period of one year (July 2015 to June 2016).Modeled water fluxes showed high temporal variability (Figure 8) and distinct differences in annual values (Table 11) among the three hydrological zones.Results showed positive and negative cumulative bottom flux among the three hydrological zones for the simulation period.Negative bottom flux indicated drainage from the profile and groundwater recharge.Moreover, positive bottom fluxes were associated with capillary rise and groundwater discharge.The highest positive cumulative bottom flux was observed at the riparian zone (296.2 mm), and the lowest value was observed at the middle zone (116.2 mm) of the floodplain during the simulation period (Figure 8).
The largest annual change in soil water storage was observed at the riparian zone.An increase of about 262.8 mm was calculated at the riparian zone, followed by the middle zone (232.9 mm), and fringe zone (120.1 mm) (Table 11).During the short (October to December 2015) and long (March to June 2016) rainy seasons, as well as the months of January and February 2016, soil water storage increased due to the filling up of the soil profile by the infiltration process.While in the dry season (July to September 2015), soil water storage gradually decreased at all of the hydrological zones.It should be noted that, at the middle zone, the decrease reached values below zero (Figure 8).
Cumulative root water uptake was highest at the fringe zone, followed by the middle and riparian zones (Figure 8).Annual cumulative root water uptake reached 201.0 mm at the fringe, 183.3 mm at the middle, and 183.5 mm at the riparian zones.Annual cumulative evaporation was substantially higher at the fringe (842.9 mm) than the riparian (785.9 mm) and middle (712.8 mm) zones during the simulation period.Therefore, actual evapotranspiration was highest at the fringe zone (1043.9mm), contributing to 76.5% of the water loss from precipitation, followed by the riparian zone (969.4 mm, 71% water loss), and the middle zone (896.1 mm, 65.7% water loss) at each hydrological zone (Table 11).

Sensitivity Analysis
Sensitivity analysis was performed to analyze the impact of changes in the depth to groundwater level on root water uptake, evaporation, and change in soil water storage.Figure 9 and Table 12 show the results of the sensitivity analysis.The sensitivity analysis revealed that among the water flux components, evaporation was the most sensitive parameter to changes in the dgwl, followed by change in soil water storage at all of the hydrological zones (Table 12 and Figure 9).Root water uptake was sensitive at the middle zone, while it was insensitive to changes in the dgwl at the fringe and riparian zones (Figure 9 and Table 12).Change in soil water storage and evaporation were highly sensitive at the fringe followed by the middle zone.The most sensitive parameter at the fringe and middle zones was change in soil water storage, while evaporation was highly sensitive at the fringe zone for all of the changes in the dgwl.Sensitivity indices increased with increase in changes to the dgwl for all of the water flux components at the three hydrological zones (Table 12).
At the fringe zone, cumulative evaporation and change in soil water storage had the maximum annual values of 1057.7 mm and 128.2 mm, higher than their respective references when the dgwl decreased by 60% and 30%, respectively.The minimum values of 696.4 mm and 65.0 mm for cumulative evaporation and change in soil water storage, respectively, were observed when the dgwl increased by 60% from the reference (Figure 9).In summary, the sensitivity analysis results indicate that very shallow or very deep groundwater levels will have great impacts on evaporation and change in soil water storage.

Sensitivity Analysis
Sensitivity analysis was performed to analyze the impact of changes in the depth to groundwater level on root water uptake, evaporation, and change in soil water storage.Figure 9 and Table 12 show the results of the sensitivity analysis.The sensitivity analysis revealed that among the water flux components, evaporation was the most sensitive parameter to changes in the dgwl, followed by change in soil water storage at all of the hydrological zones (Table 12 and Figure 9).Root water uptake was sensitive at the middle zone, while it was insensitive to changes in the dgwl at the fringe and riparian zones (Figure 9 and Table 12).Change in soil water storage and evaporation were highly sensitive at the fringe followed by the middle zone.The most sensitive parameter at the fringe and middle zones was change in soil water storage, while evaporation was highly sensitive at the fringe zone for all of the changes in the dgwl.Sensitivity indices increased with increase in changes to the dgwl for all of the water flux components at the three hydrological zones (Table 12).
At the fringe zone, cumulative evaporation and change in soil water storage had the maximum annual values of 1057.7 mm and 128.2 mm, higher than their respective references when the dgwl decreased by 60% and 30%, respectively.The minimum values of 696.4 mm and 65.0 mm for cumulative evaporation and change in soil water storage, respectively, were observed when the dgwl increased by 60% from the reference (Figure 9).In summary, the sensitivity analysis results indicate that very shallow or very deep groundwater levels will have great impacts on evaporation and change in soil water storage.

Discussion
The variations between calibrated and measured saturated hydraulic conductivity can be partly explained by the existence of macro pores [83,84] due to soil cracking, as well as by the sampling and measurement errors [85].The discrepancies between measured and modeled soil moisture content are probably indicative of soil heterogeneity, model uncertainty, and measurement errors, which are inevitable in field conditions and the model setup [34,86].Overall, the calibrated soil water model Hydrus-1D simulated the spatial soil moisture dynamics at the different hydrological zones in a reasonable manner.Therefore, the model can be used to predict spatial soil moisture dynamics as well as to close soil moisture data gaps within sites of similar biophysical properties, and rainfall in the floodplain.The results also demonstrate that Hydrus-1D is an effective tool for evaluating soil water dynamics [31,32,35,87,88] and hence would be acceptable for performing scenario simulations.The findings of the current study corroborate with those of Reference [33], who successfully applied a calibrated Hydrus-1D model with soil moisture measurements from one site to predict soil moisture conditions at another site, given that the biophysical properties and rainfall were similar.The results of Reference [89] further confirm the ability of the Hydrus model to reproduce spatial soil moisture content dynamics in floodplains, although the study applied the Hydrus-2D version.
Our observations and modeling results showed a distinct spatial pattern of the soil moisture time series along the defined hydrological zones.Soil moisture content increased from the surface layer (10-cm depth) to the bottom layer (40-cm depth).Furthermore, soil moisture content and change in soil water storage increased towards the riparian zone.Since all of the hydrological zones have similar rainfall levels, land use, and are situated a few kilometers apart from each other, the dynamics in soil moisture content and change in soil water storage can be explained by the depth to groundwater level according to the distance from the river, overbank flow, atmospheric influence, and soil properties.
Factors such as flooding from river overbank flow, rainfall, and seasonally shallow depth to groundwater level replenish soil water storage and soil moisture content at the riparian and middle zones during the long rainy season.Conversely, the fringe zone does not experience flooding from river overbank flow due to the distance from the river.Rainfall and lateral flow from the mountains are the main sources of water causing flooding [49] and replenishment of the soil water storage at the fringe zone.During the dry and short rainy seasons, soil moisture dynamics at the riparian and middle zones is mainly controlled by atmospheric influence, rainfall, soil properties, and the influence of the river water level on the depth to groundwater.Furthermore, the riparian and middle zones are directly connected to the river [49].Therefore, a decrease or increase in the river water level directly impacts the depth to groundwater level in these zones.Similar results on the direct impact of the river water level on the depth to groundwater level in floodplains have been observed by References [90][91][92][93].Soil moisture dynamics in the dry and short rainy seasons is independent of the influence from the depth to groundwater level at the riparian and middle zones of the floodplain, as indicated by both the modeled and measured time series.Similar results from Reference [43] indicated that decreasing low flows were directly linked to the surface water-groundwater interaction in the floodplain during the dry season, altering the depth to groundwater level dynamics and hence the soil water availability in the floodplain of Kilombero valley.The soil water storage after flood recession at the riparian and middle zones provides conditions for the cultivation of short-cycled upland crops like vegetables following the lowland rice; however, at the fringe zone, the soil water storage and the low depth to groundwater level throughout the year provides an opportunity for the cultivation of deep-rooted upland crops.
Furthermore, the increasing clay and soil organic carbon content towards the riparian zone could partly explain the high soil water storage observed at this zone.Our results show high soil moisture content with increasing clay content in soils, which is in line with studies by References [33,94].From the theoretical point of view, the observed increasing soil organic carbon and clay content towards the riparian zone would imply high residual and saturated water content [95,96].Subsequently, the riparian zone of the floodplain has a high agricultural production potential.
Atmospheric influence, especially of the potential evapotranspiration on soil moisture dynamics, is more pronounced at the fringe zone.The highest actual evapotranspiration at the fringe zone could partly explain the lowest soil moisture content observed at this zone.Actual evapotranspiration is always lower than the potential evapotranspiration at all of the hydrological zones, thus indicating limited water availability [97,98].The high sensitivity of actual evaporation and soil water storage to changes in the depth to groundwater level was indicative of the influence of the groundwater table.Comparing the different hydrological zones, the highest sensitivity of actual evaporation at the fringe zone could be partly explained by the relatively low depth to groundwater level observed during the study throughout the year, as a result of direct recharge and lateral flow from the surrounding mountains [49].For this reason, to reduce the loss of soil water through evaporation, soil and water management strategies like mulching can be used for sustainable agricultural production at the local scale.
The temporal soil moisture variability at the fringe zone is mainly determined by the depth to the groundwater (more so during the dry and short rainy seasons, during which soil moisture was directly responsive to changes in the depth to groundwater), soil properties, precipitation, and atmospheric influence.Similar findings were observed by Reference [99], who conducted a study in the alluvial floodplain of the Baratz Lake watershed, in the northwest of Sardinia, Italy.Their study showed that groundwater dynamics in the floodplain are an important source of temporal variability in soil moisture and vertical recharge processes in the alluvial floodplain.On the contrary, at the riparian and middle zones, soil moisture content is less responsive to the depth to groundwater level, which continuously declines as soil moisture directly responds to rain events during the short rainy season.This finding is supported by the sensitivity analysis results, which showed that variations in the depth to groundwater level have a greater impact on soil water storage at the fringe than at middle and riparian zones of the floodplain.Therefore, agricultural practices that might cause the depth to groundwater level to decline further have to be implemented with caution, especially at the fringe zones, whose soil moisture dynamics, soil water storage, and evapotranspiration are strongly controlled by the groundwater level.Our sensitivity analysis results indicate that an increase in the depth to groundwater level below the surface will decrease soil water storage and actual evapotranspiration at all of the hydrological zones.
The observed high positive cumulative bottom flux at riparian and fringe zones compared to the middle zone could be due to the contribution of groundwater flow from overbank flow towards the riparian zones and lateral flow from the mountains to the fringes.For a detailed understanding of the interacting hydrological processes among the river, floodplain, and its headwaters, the current Hydrus-1D model has to be extended to a 2D model.This will explicitly account for the lateral flows and subsurface runoff, which are not taken into account by the current model used in this study.Nonetheless, the model is adequate as the first attempt in the Kilombero floodplain to achieve the study objectives.

Conclusions
The Hydrus-1D model was calibrated and validated using daily field measurements of soil moisture content, over a period of one year in the floodplain of the Kilombero catchment, southern Tanzania.Simulated soil moisture content from the model fit quite well with the measured values at all of the four soil depths for all of the hydrological zones.Calibrated soil water model at one site can be used to predict spatial soil moisture content, as well as to close data gaps at other sites within the hydrological zone, with similar biophysical properties and rainfall levels.
Soil water storage increases towards the riparian zone.The highest evapotranspiration is observed at the fringe, followed by the riparian and middle zones.The riparian zone has the highest capillary rise and groundwater discharge, followed by the fringe and middle zones.The fringe zone maintains a relatively low depth to groundwater level throughout the year due to the lateral flow contribution from mountains.
Flooding from river overbank flow, local rainfall, and groundwater contribution replenish soil water storage during the long rainy season at the riparian and middle zones.A combination of local rainfall and lateral flow from mountains are the main sources of water supply at the fringe zone, replenishing soil water storage during the long rainy season.During the dry season, soil properties and atmospheric influence control change in soil water storage and soil moisture dynamics at the riparian and middle zones, whereas, in addition to these factors, the depth to groundwater and lateral flow from mountains play an important role in controlling soil moisture dynamics at the fringe zone.In the short rainy seasons, local rainfall and soil properties control soil moisture dynamics at the riparian and middle zones.Local rainfall, soil properties, lateral flow from mountains, and depth to groundwater level control soil moisture dynamics at the fringe zone.
Actual evaporation and change in soil water storage are highly sensitive to changes in the depth to groundwater level at all of the hydrological zones.It is worth noting that the sensitivity of actual evaporation and change in soil water storage is highest at the fringe zone.
The results help gain a better understanding of the groundwater-soil water interactions, and may provide references for wetland conservation and sustainable agricultural water management.
The calibrated and validated model can be used for climate and land use scenario studies.This study has shown that local-scale water fluxes and soil moisture availability are driven by catchment-scale processes.Flooding of the riparian zone depends on catchment-scale processes, which are sensitive to land use and climate change [43].Lateral fluxes and groundwater availability at the fringe will be influenced by land use and climate change.However, the one of the limitations of the study was it inability to describe the lateral flow; henceforth, in analyzing future water availability, the 1D modeling approach used in this study has to be linked to a 3D catchment-scale model that is able to provide the lower boundary conditions required to run Hydrus-1D as well as account for the lateral flows within the floodplain.

Figure 1 .
Figure 1.(a) Location of study area andoverview of hydrological zones and instrumentation; and (b) detailed instrumentation at each hydrological zone used for modeling.

Figure 1 .
Figure 1.(a) Location of study area andoverview of hydrological zones and instrumentation; and (b) detailed instrumentation at each hydrological zone used for modeling.

Figure 2 .
Figure 2. Daily rainfall and ET0 values used during the simulation period.

Figure 2 .
Figure 2. Daily rainfall and ET 0 values used during the simulation period.

Water 2018 ,
10, x FOR PEER REVIEW 11 of 27

Figure 3 .
Figure 3. Measured and modeled soil moisture content for calibration at the four plots in the riparian zone.(a,b) plot PR_4; (c,d) plot PR_1; (e,f) plot PR_8; (g,h) plot PR_3.dgwl, measured depth to groundwater level.

Figure 4
Figure4presents the time series measured depth to groundwater level and calibration results of soil moisture content at the middle hydrological zone of the floodplain.Results revealed that the modeled soil moisture content exhibited variation tendencies quite similar to the measured data at the different soil depths in all of the plots.Modeled and measured soil moisture content increased with increase in soil depth.Overall, there was a good agreement between modeled and measured soil moisture content after calibration (R 2 = 0.59-0.89,NSE = 0.51-0.88,and RMSE = 0.02-0.10cm 3 cm 3 ) (Table4).The best fit of the van Genuchten hydraulic parameters for individual plots are shown in Table5.There was almost no difference observed between the calibrated hydraulic parameters among the plots, except for saturated hydraulic conductivity.Hence, small variations between the modeled and measured soil moisture content values were noted among the plots.

Figure 3 .
Figure 3. Measured and modeled soil moisture content for calibration at the four plots in the riparian zone.(a,b) plot PR_4; (c,d) plot PR_1; (e,f) plot PR_8; (g,h) plot PR_3.dgwl, measured depth to groundwater level.

Figure 4
Figure4presents the time series measured depth to groundwater level and calibration results of soil moisture content at the middle hydrological zone of the floodplain.Results revealed that the modeled soil moisture content exhibited variation tendencies quite similar to the measured data at the different soil depths in all of the plots.Modeled and measured soil moisture content increased with increase in soil depth.Overall, there was a good agreement between modeled and measured soil moisture content after calibration (R 2 = 0.59-0.89,NSE = 0.51-0.88,and RMSE = 0.02-0.10cm 3 /cm 3 ) (Table4).The best fit of the van Genuchten hydraulic parameters for individual plots are shown in Table5.There was almost no difference observed between the calibrated hydraulic parameters among the plots, except for saturated hydraulic conductivity.Hence, small variations between the modeled and measured soil moisture content values were noted among the plots.

Figure 4 .
Figure 4. Measured and modeled soil moisture content after calibration at the four plots for the middle zone.(a,b) plot PM_3; (c,d) plot PM_5; (e,f) plot PM_7.dgwl; measured depth to groundwater level.

Figure 4 .
Figure 4. Measured and modeled soil moisture content after calibration at the four plots for the middle zone.(a,b) plot PM_3; (c,d) plot PM_5; (e,f) plot PM_7.dgwl; measured depth to groundwater level.

Figure 5 .
Figure 5. Measured and modeled soil moisture content after calibration at the four plots for the fringe zone.(a,b) plot PF_12; (c,d) plot PF_10; (e,f) plot PF_4; (g,h) plot PF_8.dgwl; measured depth to groundwater level.

Figure 5 .
Figure 5. Measured and modeled soil moisture content after calibration at the four plots for the fringe zone.(a,b) plot PF_12; (c,d) plot PF_10; (e,f) plot PF_4; (g,h) plot PF_8.dgwl; measured depth to groundwater level.

Figure 6 .
Figure 6.Comparison of measured (averaged from all plots) and modeled soil moisture from the averaged calibrated hydraulic parameters at each hydrological zone.(a-c) Riparian zone; (d-f) middle zone; (g-i) fringe zone.dgwl: measured depth to groundwater level.

Figure 6 .
Figure 6.Comparison of measured (averaged from all plots) and modeled soil moisture from the averaged calibrated hydraulic parameters at each hydrological zone.(a-c) Riparian zone; (d-f) middle zone; (g-i) fringe zone.dgwl: measured depth to groundwater level.

Figure 9 .
Figure 9. Sensitivity analysis for different water fluxes on changes in the depth to groundwater level; dgwl, depth to groundwater level.

Figure 9 .
Figure 9. Sensitivity analysis different water fluxes on changes in the depth to groundwater level; dgwl, depth to groundwater level.

Table 1 .
Measured soil physical properties at the three hydrological zones of the floodplain.

Table 2 .
Statistical measures of the Hydrus-1D model performance for simulations of soil moisture content at the riparian zone after calibration.

Table 3 .
Calibrated van Genuchten hydraulic parameters for the different plots at the riparian zone of the floodplain.

Table 4 .
Statistical measures of the Hydrus-1D model performance for simulations of soil moisture content at the middle zone after calibration.

Table 5 .
Calibrated van Genuchten hydraulic parameters for the different plots at the middle zone of the floodplain.

Table 4 .
Statistical measures of the Hydrus-1D model performance for simulations of soil moisture content at the middle zone after calibration.

Table 5 .
Calibrated van Genuchten hydraulic parameters for the different plots at the middle zone of the floodplain.

Table 6 .
Statistical measures of the Hydrus-1D model performance for simulations of soil moisture content at the fringe zone after calibration.

Table 6 .
Statistical measures of the Hydrus-1D model performance for simulations of soil moisture content at the fringe zone after calibration.

Table 7 .
Calibrated van Genuchten hydraulic parameters from the different plots at the fringe zone of the floodplain.

Table 8 .
Averaged calibrated hydraulic parameters at each hydrological zone.

Table 8 .
Averaged calibrated hydraulic parameters at each hydrological zone.

Table 9 .
Hydrus-1D performance for soil moisture simulation after using averaged hydraulic parameters calculated from all of the calibrated plots for each hydrological zone.

Table 10 .
Hydrus-1D performance for soil moisture simulation after validation.

Table 11 .
Annual water balance from the Hydrus-1D method for each hydrological zone.Gi, capillary rise/groundwater discharge; ΔS, change in soil water storage.

Table 11 .
Annual water balance from the Hydrus-1D method for each hydrological zone.ET p , potential evapotranspiration; ET a , actual evapotranspiration; G o , groundwater recharge; G i , capillary rise/groundwater discharge; ∆S, change in soil water storage.

Table 12 .
Sensitivity indices (SI) for different water fluxes calculated after changes in the dgwl.
Numbers in bold are the most sensitive parameters at each hydrological zone.

Table 12 .
Sensitivity indices (SI) for different water fluxes calculated after changes in the dgwl.
Numbers in bold are the most sensitive parameters at each hydrological zone.