Aquifer parameter and zone structure estimation using kernel-based fuzzy c-means clustering and genetic algorithm M. Tamer Ayvaz a , Halil Karahan a , Mustafa M. Aral b, * a Department of Civil Engineering, Pamukkale University, TR-20070 Denizli, Turkey b Multimedia Environmental Simulations Laboratory, School of Civil and Environmental Engineering, Georgia Institute of Technology, Atlanta, GA 30332-0355, USA Received 29 January 2007; received in revised form 23 May 2007; accepted 27 June 2007 KEYWORDS Inverse problems; Parameter estimation; Groundwater modeling; Kernel based fuzzy c- means clustering; Zone structure; Genetic algorithm Summary In this study, we propose an inverse solution algorithm through which both the aquifer parameters and the zone structure of these parameters can be determined based on a given set of observations on piezometric heads. In the zone structure identification problem, kernel-based fuzzy c-means (KFCM) clustering method is used. The association of the zone structure with the transmissivity distribution is accomplished through a coupled simulation–optimization model. In the optimization model, genetic algorithm (GA) is used due to its efficiency in finding global or near global optimum solutions. Since the solution is based on the GA procedures, the optimization process starts with a randomly generated initial solution. Thus, there is no need to define an initial estimate of the solution. This is an advantage when compared to other studies reported in the literature. Further, the objective function used in the optimization model does not include a reference to field transmissivity data, which is another advantage of the proposed methodology. Numerical examples are provided to demonstrate the performance of the proposed algorithm. In the first example, transmissivity values and zone structures are determined for a known num- ber of zones in the solution domain. In the second example, optimum number of zones as well as the transmissivity values and the zone structures are determined iteratively. A sensitivity analysis is also performed to test the performance of the proposed solution algorithm based on the number of observation data necessary to solve the problem accu- rately. Numerical results indicate that the proposed algorithm is effective and efficient and may be used in the inverse parameter estimation problems when both parameter val- ues and zone structure are unknown. ª 2007 Elsevier B.V. All rights reserved. 0022-1694/$ - see front matter ª 2007 Elsevier B.V. All rights reserved. doi:10.1016/j.jhydrol.2007.06.018 * Corresponding author. Tel.: +1 404 894 2243; fax: +1 404 894 5111. E-mail address: maral@ce.gatech.edu (M.M. Aral). Journal of Hydrology (2007) 343, 240– 253 available at www.sciencedirect.com journal homepage: www.elsevier.com/locate/jhydrol