In this study, we attempt to establish an alternative method for estimating the groundwater levels and the specific yields of an unconfined aquifer for the evaluation of potential groundwater reservoirs. We first converted the inverted resistivity into the normalized water content. Then, we inverted the parameters of the Brooks-Corey model from the vertical profiles of the water content by assuming that the suction head was in proportion to the elevation regarding a predefined base level. Lastly, we estimated the groundwater level, the theoretical specific yield, and the specific yield capacity from the Brooks-Corey parameters at every survey site in the study area. The contour maps of the time-lapse groundwater levels show that the groundwater flows downstream, with a higher hydraulic gradient near the river channel than in the area away from the main channel. We conclude that the estimated maximum specific yield capacities are consistent with that derived from the pumping tests in the nearby observation well. Additionally, the specific yield capacities are only three quarters to two thirds of the theoretical specific yields derived from the difference between the residual and saturated water contents in the Brooks-Corey model. We conclude that the distribution pattern of the specific yields had been subjected to the distribution of natural river sediments in the Minzu Basin, since the modern channel was artificially modified. Although we had to make some simple assumptions for the estimations, the results show that the surface resistivity surveys provide reasonable estimations of the hydraulic parameters for a preliminary assessment in an area with few available wells.

This study aims to apply geophysical methods to determine the Specific Yield (Sy) and Groundwater Level (GWL) in an unconfined aquifer of the Pingtung Plain in South Taiwan. Sy is an important hydraulic parameter for assessing groundwater potential. Obtaining specific yield for a large area is impractical due to the limited coverage and the high cost of the pumping test, which limits the potential evaluation of regional groundwater. Therefore, we used time-lapse Electrical Resistivity Imaging (ERI) to determine the Sy and GWL. Seasonal variations were considered when measuring time-lapse resistivity for five different months in 2019. We calculated the Sy and GWL from inverted resistivity data using empirical formulas and the soil–water characteristic curve (SWCC). We first used Archie’s law to calculate the relative saturation change with depth for each ERI profile, and then we used the Van Genuchten (VG) and Brooks–Corey (BC) empirical equations to estimate Sy and GWL. Finally, we compared the obtained GWL to the existing observation well to verify the findings of our study. The results showed that the VG and BC are able to predict Sy and GWL; however, the BC result is less consistent with the observation well result. In the study area, the dry season GWL ranged from 24.5 m to 35.2 m for the VG results and from 25.7 m to 35.5 m for the BC results. The wet season GWL ranged from 26.5 m to 38.9 m for the VG and from 26.4 m to 38.2 m for the BC results. The spatial distribution of the GWL shows a high gradient of GWL in the northeastern region, induced by significant proximal fan recharge. The determined spatial distribution of Sy varies from 0.15 to 0.21 for the VG and 0.14 to 0.20 for the BC results, indicating the study area has significant potential for groundwater resources. Therefore, nondestructive resistivity imaging can be used to aid in the determination of hydraulic parameters.

This work develops a rule curve-based conjunctive use management model for optimizing the operating rules for a lake–groundwater system with off-stream storage lakes. The proposed procedure is a simulation-optimization approach that embeds an Artificial Neural Network (ANN) instead of a groundwater numerical model into a genetic algorithm (GA). The direct physical exchange between lake water with groundwater is simulated using the ANN model, which is a reduced version of a full numerical model, MODFLOW with an LAK3 module. By applying the ANN model, the proposed procedure can reduce the computational burden that is induced by the nonlinear exchange. An operating rule-based optimal conjunctive use management model for the Gaopin Artificial lakes system in Taiwan was thus developed using the proposed framework. A set of optimal solutions involves rule curves and a discount ratio. Simulation results demonstrate that the embedded ANN model can accurately simulate the nonlinear exchange of a lake with groundwater. The embedded ANN model is less computationally complex than the numerical model. This work demonstrates a methodology for reducing the computational burden of the optimal conjunctive use management model that is associated with an internal nonlinear system by using the ANN reduced model. Specifically, the concept of, and results obtained using the developed operating rule-based model incorporating five artificial lakes and considering the nonlinear exchange of those lakes with the groundwater system provides a valuable practical reference for solving related conjunctive use problems.

Hydraulic conductivity (K) is crucial for groundwater studies and is conventionally obtained through pumping tests, which are costly due to well installation, resulting in a limited amount of data. Recent studies have addressed this limitation by estimating K through the integration of a pumping-test K and electrical resistivity measurements. While this approach is suitable for local areas, it cannot readily be applied to determine K fields for the composite fan delta. This study proposes and demonstrates an advanced method to estimate the K’s spatial distribution of a composite fan delta. The proposed method included data classification, linear regression, and kriging interpolation. Data classification was conducted using a physical-based zonation method. The K and formation factor (F) data pairs were classified into several groups. Linear regression was used to develop K-F mapping for each group. The regression results showed a good correlation between K and F in each group. These K-F mappings were verified by additional pumping tests. These results indicate that estimation errors were between 7 m/day and 58 m/day, and the correlation coefficient of each data group was greater than 0.8. Based on these regression equations and ordinary kriging method, the detailed K spatial distribution of the study area was derived. According to these results, the proposed method is efficient and economical to delineate K’s spatial features for regional groundwater systems and can benefit groundwater-related studies.