An areal numerical simulation has been carried out to investigate effects of below-sea-level (BSL) excavation in the raw material site of a cement factory in Turkey. A finite element model (585 nodes with 534 elements) is formed to solve for the head distribution in the quarry site upon quarry operation planned to be implemented in the near future. The model is calibrated to the field conditions and appropriate boundary conditions and the physical parameters are obtained to be used in future prediction studies. After a successful calibration the model is run to estimate the water levels and the discharge rates required during below sea level quarry operations. Above sea level (ASL) and below sea level (BSL) operations are simulated and water level contour maps are obtained both for above sea level (ASL) production for the 2000-2030 period, and for each BSL (-10m, -20m, -30m) production periods, which would totally take 13 years. Estimation show that the proposed model runs properly and it calculates the water levels and discharge rates accurately for probable future quarry operations. It is clear that quarry operations would not create a serious problem in terms of water discharge from the quarry site.