Excessive exploitation of groundwater has led to a significant subsidence in Darab plain, located in Fars province, Iran. The compaction of aquifer layers especially in fine-grain sediments leads to changes in aquifer parameters including hydraulic conductivity, specific yield, and compressibility which consequently changes the permeability and storage of the aquifer. Accordingly, a precise estimation of aquifer parameters are essential for both future water resources planning and management. The main objective of this paper is to determine the most optimum values of aquifer parameters based on the groundwater information and subsidence measurements spanning from 2010 to 2016. The proposed approach for inverse solution of groundwater flow is based on Lattice Boltzmann Method (LBM) integrated with Genetic Algorithm for optimization (GA). Subsidence measurements made by previous studies as well as the initial values of unknown parameters obtained from piezometric information are incorporated into the inverse modeling. The whole process of inverse modeling is repeated from 2010 to 2016 which leads to annually estimation of the aquifer parameters. Due to the compaction occurring in the aquifer system, a decreasing temporal trend is observed in the aquifer parameters in most parts of the plain. By fitting a function to time-dependent aquifer parameters, their corresponding values and consequently the amount of subsidence in 2017 are predicted. The small average relative error (~ 3.5%) between the predicted land subsidence and the measurements demonstrates the high performance of the proposed inverse modeling approach.