Izeh Plain Subsidence Modeling Using MODFLOW Mathematical Code

The subsidence of the plains structure can be directly due to the drop in groundwater level and the destruction of aquifer alluvial tissue. In this study, in the Izeh plain, the MODFLOW mathematical model was performed by extracting a raster map of saturated aquifer subsidence in the ten-year period from 2008 to 2019 and the forecast period for three climatic scenarios between 2019 and 2029. The subsidence model was developed by developing the SUB package in the limited numerical code structure structure of MODFLOW. The quantitative groundwater flow model was calibrated using PEST automation code and sensitivity analysis. The results show that the mathematical model used to simulate Izeh aquifer has a relative relative error NRMSE, which confirms the ideal modeling after reviewing the validation process. An examination of the vertical changes in the structure of the land showed that in the long term of 20 years, with the climatic assumptions raised, the aquifer level will occur a total of up to 1.5 meters and a subsidence of at least 0.9 meters. If the aquifer alluvial structure is destroyed, this amount may be less than these figures, but in return, it is equal to the conditions without returning the groundwater aquifer feed.