Since 1900, Semarang City has been meeting its industrial water needs by pumping groundwater through its underlying aquifers. The trend toward exploiting groundwater resources has driven the number of deep wells and their production capacity to increase, and therefore leads to the water table to drop from time to time, which has been marked as one of the primary causes of land subsidence there. The main aim of the current study was to numerically model the temporal and spatial evolution of groundwater table under excess abstraction so that a groundwater management strategy can be accordingly drawn up for ensuing the sustainability of groundwater resources in the future. A series of numerical simulations were carried out to take into account hydrogeological data, artificial and natural discharges of deep wells, and boundary effects in Semarang City. The groundwater modeling is calibrated under two flow conditions of the steady state from 1970 to 1990 and the transient state from 1990 to 2005 for six observation wells distributed in Semarang City. Four scenarios that reflect potential management strategies were developed, and then their effectiveness was systematically investigated. The results of our study indicate that the implementation of proper groundwater control management and measure is able to restore the groundwater level to rise back in Semarang City, and in turn achieve the sustainability of groundwater resources.

This study was carried out to assess land subsidence due to excessive groundwater abstraction in the northern region of Semarang City by integrating the application of both numerical models and geodetic measurements, particularly those based on the synthetic aperture radar interferometry (InSAR) technique. Since 1695, alluvial deposits caused by sedimentations have accumulated in the northern part of Semarang City, in turn resulting in changes in the coastline and land use up to the present. Commencing in 1900, excessive groundwater withdrawal from deep wells in the northern section of Semarang City has exacerbated natural compaction and aggravated the problem of land subsidence. In the current study, a groundwater model equivalent to the hydrogeological system in this area was developed using MODFLOW to simulate the hydromechanical coupling of groundwater flow and land subsidence. The numerical computation was performed starting with the steady-state flow model from the period of 1970 to 1990, followed by the model of transient flow and land subsidence from the period of 1990 to 2010. Our models were calibrated with deformation data from field measurements collected from various sources (e.g., leveling, GPS, and InSAR) for simulation of land subsidence, as well as with the hydraulic heads from observation wells for simulation of groundwater flow. Comparison of the results of our numerical calculations with recorded observations led to low RMSEs, yet high R² values, mathematically indicating that the simulation outcomes are in good agreement with monitoring data. The findings in the present study also revealed that land subsidence arising from groundwater pumping poses a serious threat to the northern part of Semarang City. Two groundwater management measures are proposed and the future development of land subsidence is accordingly projected until 2050. Our study shows quantitatively that the greatest land subsidence occurs in Genuk District, with a magnitude of 36.8 mm/year. However, if the suggested groundwater management can be implemented, the rate and affected area of land subsidence can be reduced by up to 59% and 76%, respectively.

Land subsidence is a global threat to coastal areas worldwide, including the North Java coastal area. Of many areas experiencing land subsidence in North Java, the rate of land subsidence in Pekalongan has matched the high subsidence rates usually found in big cities. The rate of land subsidence in Pekalongan far exceeds the sea-level rise, resulting in a looming threat of land loss. The devastating impacts of land subsidence are the manifestation of its subsurface movement. Therefore, it is essential to understand the subsurface to elucidate the mechanism of land subsidence. Previous studies on land subsidence in Pekalongan are mainly related to subsidence rate monitoring and have not elaborated on the subsurface condition. This paper reviews the Pekalongan subsurface geology based on available literature to provide insight into the land subsidence problem. The results revealed that the land subsidence occurs in the recent alluvial plain of Pekalongan, consisting of a 30-70 m thick compressible deposit. Possible mechanisms of land subsidence arise from natural compaction, over-exploitation of confined groundwater, and increased built areas. As the seismicity of the study area is low, tectonic influence on land subsidence is considered negligible. It is expected that the offshore, nearshore, and swamp deposits are still naturally compacting. As the surface water supply is minimal, over-exploitation of groundwater resources from the deltaic and Damar Formation aquifers occurs. In the end, future research direction is proposed to reduce the impacts of the subsidence hazard.

Semarang-Demak and other cities along the coast of North Java are vulnerable to land subsidence. The presence of saline groundwater in the coastal region is thought to affect the high subsidence rate, in this case the compressibility of the aquitard layer. We aimed to analyze the effects of groundwater salinity on the compression characteristics of the Semarang-Demak clay using physical analysis. Methods included the determination of groundwater salinity, clay mineralogy and fabrics, and consolidation tests under various salinity conditions. The Semarang-Demak clay is dominated by smectite of high activity, and saline clay exists at the depth of 10 to 35 m. Consolidation tests revealed that the increase insalinity increases the average consolidation rate and hydraulic conductivity up to 42% and 37.5%, respectively. Clay fabric analysis showed that the groundwater salinity modified the interconnectivity of pores by changing the fabric into parallel alignments, facilitating faster porewater dissipation, hence the clay is more readily compressed. These findings are useful for explaining the mechanism of the fast-subsiding coastal plains of North Java.

The Semarang-Demak plain has experienced intense human intervention over the last 40 years, thereby causing land subsidence. This study aims to assess long-term conditions in the study area using the drivers-pressuresstate- impacts-response (DPSIR) framework to mitigate land subsidence. Methods include analysis of land subsidence, socioeconomic, surface, and subsurface data, as well as spatial analysis. Results show that rapid population growth and economic activities are major driving forces, manifesting as pressures exerted from overexploitation of groundwater, increasing building and infrastructure loads, and decreasing non-built areas. Groundwater overexploitation reduced the artesian pressure in the 1980s, forming depression cones of the groundwater level from 5 to 30 m below mean sea level. From 1984 to the present, the constructed areas have increased more than tenfold, with Semarang City possessing the most densely built area. Based on our findings, we propose responses consisting of surface water utilization, spatial building regulation, and rigorous groundwater and land subsidence monitoring. Moreover, we encourage the strengthening of law enforcement and inter-sectoral management to ensure the successful land subsidence mitigation.