Hydrogeological Characterization of Coastal Aquifer on the Basis of Observed Sea Level and Groundwater Level Fluctuations: Neretva Valley Aquifer, Croatia

Abstract

Hydrogeological data availability is often limited to local areas where usual in situ tests or methods are applied (slug/bail or pumping tests, Electrical Resistivity Tomography (ERT)). Because most problems (e.g., saltwater intrusion mitigation) require problem analysis on larger scales (catchment or sub catchment), hydrogeological identification of global character is preferable. This work leads to the determination of aquifer hydrogeological parameters on the basis of observed sea level, groundwater piezometric head found inland, and barometric pressure. When applied to observed signals, the approach led efficiently to final hydrogeological characterization. After identification of dominant tidal constituents from observed signals, barometric efficiency was successfully determined. Following available information on geological settings, an appropriate conceptual model was applied and updated to count for polychromatic signals. Final determination of hydrogeological parameters relied on root mean square error (RMSE) minimization and led to determination of (i) presence of three stratigraphic units: unconfined sandy aquifer on the top, a confining layer made of clay, and a confined gravel layer; (ii) existence of the clay layer under the sea with a total length of 1400 m; (iii) a clay layer has been identified as confining one by both spectral analysis and determined leakance value; and (iv) estimated confined aquifer specific storage ranging from 2.87 × 10−6 to 4.98 × 10−6 (m−1), whereas hydraulic conductivity ranged from 7.0 × 10−4 to 7.5 × 10−3 (m s−1). Both range intervals corresponded to previous in situ findings conducted within the area of interest.