The integration of thermal infrared imaging,discharge measurements and numerical simulation to quantify the relative contributions of freshwater inflows to small estuaries in Atlantic Canada

Nutrient fluxes from developed catchments are often a significant factor in the declining water quality and ecological functioning in estuaries. Determining the relative contributions of surface water and groundwater discharge to nutrient‐sensitive estuaries is required because these two pathways may be characterized by different nutrient concentrations and temporal variability, and may thus require different remedial actions. Quantifying the volumetric discharge of groundwater, which may occur via diffuse seepage or springs, remains a significant challenge. In this contribution, the total discharge of freshwater, including groundwater, to two small nutrient‐sensitive estuaries in Prince Edward Island (Canada) is assessed using a unique combination of airborne thermal infrared imaging, direct discharge measurements in streams and shoreline springs, and numerical simulation of groundwater flow. The results of the thermal infrared surveys indicate that groundwater discharge occurs at discrete locations (springs) along the shoreline of both estuaries, which can be attributed to the fractured sandstone bedrock aquifer. The discharge measured at a sub‐set of the springs correlates well with the area of the thermal signal attributed to each discharge location and this information was used to determine the total spring discharge to each estuary. Stream discharge is shown to be the largest volumetric contribution of freshwater to both estuaries (83% for Trout River estuary and 78% for McIntyre Creek estuary); however, groundwater discharge is significant at between 13% and 18% of the total discharge. Comparison of the results from catchment‐scale groundwater flow models and the analysis of spring discharge suggest that diffuse seepage to both estuaries comprises only about 25% of the total groundwater discharge. The methods employed in this research provide a useful framework for determining the relative volumetric contributions of surface water and groundwater to small estuaries and the findings are expected to be relevant to other fractured sandstone coastal catchments in Atlantic Canada. Copyright © 2009 John Wiley & Sons, Ltd.     

The Interface Configuration of the Fresh‐/Dead Sea Water – Theory and Measurements

The high‐density Dead Sea water (1.235 g/cm³) forms a special interface configuration with the fresh groundwater resources of its surrounding aquifers. The fresh groundwater column beneath its surroundings is around one tenth of its length compared to oceanic water. This fact alone indicates the vulnerability of the fresh groundwater resources to the impacts of changes in the Dead Sea level and to saltwater migration. Ghyben‐Herzberg and Glover equations were used to calculate the volumes of water in coastal aquifers which were replaced by freshwater due to the interface seaward migration as a result of the drop in the level of the Dead Sea. For that purpose, the dynamic equation of Glover approach has been integrated to accommodate that type of interface readjustment. The calculated amounts of freshwater which substituted salt Dead Sea water due to the migration of interface are 3.21 · 10¹¹ m³, from a Dead Sea level of –392 m to τ411 m below sea level. The average porosity of coastal aquifers was calculated to range from 2.8 to 2.94%. Geoelectric sounding measurements showed that areas underlying the coastal aquifers formerly occupied by the Dead Sea water are gradually becoming flushed and occupied by freshwater. The latter is becoming salinized due to the residuals of Dead Sea water in the aquifer matrix, the present salinity of which is lower than that of the Dead Sea water. At the same time salt dissolution from the Lisan Marl formation is causing collapses along the shorelines in the form of sinkholes, tens of meters in diameter and depth.     

Groundwater connections under a barrier beach: A case study in the Venice Lagoon

The flow of groundwater beneath barrier islands has been cited as a possible pathway for salt water and chemical exchange between a protected embayment and the open sea. Evidence is presented that identifies an exchange of groundwater through a highly permeable paleoinlet along the barrier beach of Cavallino, which separates the northern Venice Lagoon from the Adriatic Sea. We utilized both point measurements of submarine groundwater discharge (SGD) and a geophysical investigation of the subsurface resistivity to analyze the movement of saline groundwater. Discharge of groundwater and associated nutrients, was higher at the site of a former inlet than at a similar site along the barrier and modulated by the difference in tidal water level between the lagoon and Adriatic Sea. If the measured conditions are typical, storm surge barriers could potentially result in a saline groundwater flow of up to 1.5×10⁶ m³ d⁻¹ into the lagoon.     

Evaluating the source and seasonality of submarine groundwater discharge using a radon-222 pore water transport model

Pore water radon (²²²Rn) distributions from Indian River Lagoon, Florida, are characterized by three zones: a lower zone where pore water ²²²Rn and sediment-bound radium (²²⁶Ra) are in equilibrium and concentration gradients are vertical; a middle zone where ²²²Rn is in excess of sediment-bound ²²⁶Ra and concentration gradients are concave-downward; and an upper zone where ²²²Rn concentration gradients are nearly vertical. These ²²²Rn data are simulated in a one-dimensional numerical model including advection, diffusion, and non-local exchange to estimate magnitudes of submarine groundwater discharge components (fresh or marine). The numerical model estimates three parameters, fresh groundwater seepage velocity, irrigation intensity, and irrigation attenuation, using two Monte Carlo (MC) simulations that (1) ensure the minimization algorithm converges on a global minimum of the merit function and the parameter estimates are consistent within this global minimum, and (2) provide 90% confidence intervals on the parameter estimates using the measured ²²²Rn activity variance. Model estimates of seepage velocities and discharge agree with previous estimates obtained from numerical groundwater flow models and seepage meter measurements and show the fresh water component decreases offshore and varies seasonally by a factor of nine or less. Comparison between the discharge estimates and precipitation patterns suggests a mean residence time in unsaturated and saturated zones on the order of 5 to 7 months. Irrigation rates generally decrease offshore for all sampling periods. The mean irrigation rate is approximately three times greater than the mean seepage velocity although the ranges of irrigation rates and seepage velocities are the same. Possible mechanisms for irrigation include density-driven convection, wave pumping, and bio-irrigation. Simulation of both advection and irrigation allows the separation of submarine groundwater discharge into fresh groundwater and (re)circulated lagoon water.     

The water balance of the Dead Sea: an integrated approach

The Dead Sea is the lowest spot on Earth. It is a closed saline lake located in the middle of the Jordan Rift Valley between Lake Tiberias and the Red Sea. Its major tributaries are the Jordan River itself and the Dead Sea side wadis. The Dead Sea has a unique ecosystem and its water has curative, industrial and recreational significance. The level of the Dead Sea has been continuously falling since the early 1930s at an average rate of 0·7 m per year. The water level, as of February 1998, is about 410·9 m below mean sea level. In this paper, a water balance model is developed for the Dead Sea by considering different hydrological components of this water balance, including precipitation, runoff, evaporation and groundwater flow. This model is calibrated based on historical levels of the Dead Sea. Different scenarios are investigated, including the proposed Dead Sea–Red Sea Canal. This project is supposed to halt the shrinking of the Dead Sea and restore it to pre‐1950 levels in the next century. Copyright © 2000 John Wiley & Sons, Ltd.     

Water balance estimation under the challenge of data scarcity in a hyperarid to Mediterranean region

Water budget analyses are important for the evaluation of the water resources in semiarid and arid regions. The lack of observed data is the major obstacle for hydrological modelling in arid regions. The aim of this study is the analysis and calculation of the natural water resources of the Western Dead Sea subsurface catchment, one which is highly sensitive to rainfall resulting in highly variable temporal and spatial groundwater recharge. We focus on the subsurface catchment and subsequently apply the findings to a large‐scale groundwater flow model to estimate the groundwater discharge to the Dead Sea. We apply a semidistributed hydrological model (J2000g), originally developed for the Mediterranean, to the hyperarid region of the Western Dead Sea catchment, where runoff data and meteorological records are sparsely available. The challenge is to simulate the water budget, where the localized nature of extreme rainstorms together with sparse runoff data results in few observed runoff and recharge events. To overcome the scarcity of climate input data, we enhance the database with mean monthly rainfall data. The rainfall data of 2 satellites are shown to be unsuitable to fill the missing rainfall data due to underrepresentation of the steep hydrological gradient and temporal resolution. Hydrological models need to be calibrated against measured values; hence, the absence of adequate data can be problematic. Therefore, our calibration approach is based on a nested strategy of diverse observations. We calculate a direct surface runoff of the Western Dead Sea surface area (1,801 km²) of 3.4 mm/a and an average recharge (36.7 mm/a) for the 3,816 km² subsurface drainage basin of the Cretaceous aquifer system.     

Earthquake-related Changes in Groundwater Levels at the Dogo Hot Spring,Japan

The Dogo hot spring, situated in Matsuyama City, Ehime Prefecture, Japan, is one of the oldest and most famous hot springs in Japan. The groundwater level or discharge at the spring decreased four times during the past eight or nine Nankai earthquakes. These are large interplate earthquakes that have occurred repeatedly in the western part of the Nankai Trough at intervals of 100–200 years since A.D. 684. To clarify the mechanism of these earthquake-related changes in the water level at the spring, we analyzed groundwater-level data recorded at the spring immediately after the 1946 Nankai earthquake and over the period from 1985 to 2006. We detected the other nine postseismic increases in groundwater level and no decreases, except for a large decrease of 11.4 m related to the 1946 Nankai earthquake. The increases were probably caused by ground-shaking, while the decrease was caused by a change in coseismic volumetric strain. These results lead to the following explanation of the recorded earthquake-related changes in the groundwater level at the Dogo hot spring. Both coseismic changes in volumetric strain and ground-shaking can lead to postseismic changes in groundwater pressure. The increase in groundwater pressure arising from ground-shaking is generally greater than the change in pressure associated with changes in coseismic volumetric strain; however, at the time of the Nankai earthquakes, the spring experiences a large increase in coseismic volumetric strain, leading to a considerably larger decrease in the groundwater level than the increase associated with ground-shaking. Therefore, the groundwater level at the Dogo hot spring usually increases at times of relatively large earthquakes, although the groundwater level or discharge decreases in the case of the Nankai earthquakes.     

The enhancing effect on tidal signals of a submarine spring connected to a semi-infinite confined aquifer

Abstract

Submarine springs play an important role in submarine groundwater discharge (SGD). To investigate the effects of these springs on the propagation of tidal signals in coastal confined aquifers, this paper considers a general coastal aquifer system with a submarine spring on the seabed where the length of the aquifer's offshore extent is finite and its submarine outlet is covered by an impermeable outlet-capping. An approximate analytical solution is obtained for describing the tidal head fluctuations in the aquifer. Solution analyses indicate that the error of the approximate analytical solution is negligible when both distances from the spring hole to the coastline and to the submarine outlet-capping are much greater than the radius of the spring hole. Sensitivity tests are conducted to investigate the effects of hydraulic properties, tidal and spring geometric configuration parameters on the tidal signal propagation in the inland aquifer. For aquifers with infinite offshore length, or without submarine springs, existing solutions in the literature are obtained. The comparison of groundwater head fluctuations for the cases with and without a submarine spring demonstrate the enhancing effect of the submarine spring on tidal signal propagation in the inland aquifer. Three situations that fit our model assumptions are given for future potential applications. A hypothetical example is used to show the possibility of identifying a spring's location using the present analytical solution together with tidal signals observed from inland wells.

Editor D. Koutsoyiannis; Associate editor Y. Guttmann

Citation Xia, Y.Q., Li, H.L., Yang, Y., and Huang, W., 2012. Enhancing effect on tidal signals of a submarine spring related to a semi-infinite confined aquifer. Hydrological Sciences Journal, 57 (6), 1231–1248.     

Changes in the Dead Sea Level and their Impacts on the Surrounding Groundwater Bodies

In this paper the reaction of the salt‐/freshwater interface due to the changes in the Dead Sea level are elaborated at in details by using the inflows into the Dead Sea, the outflows due to evaporation losses and artificial discharges, and the hydrographic registrations of the Dead Sea level. The analyses show that the interface seaward migration resulted in a groundwater discharge of around 423 Mio m³ per meter drop in the level of the Dead Sea in the period 1994–1998 and of around 525 Mio m³/m in the period 1930–1937. The additional amount of groundwater joining the Dead Sea due to the interface seaward migration was 51 Mio m³ per one square kilometer of shrinkage in the area of the Dead Sea in the period 1930–1937 and 91 Mio m³/km² in the period 1994–1998. The riparian states of the Dead Sea are nowadays loosing 370 Mio m³/a of freshwater to the Dead Sea through the interface readjustment mechanisms as a result of their over exploitation of waters which formerly fed the Dead Sea.     

Remotely sensed estimation of water discharge into the rapidly dwindling Dead Sea

Abstract

The Dead Sea is shrinking as its water level drops at the alarming rate of about 1 m year^-1. The Dead Sea is important to the economies of Israel and Jordan due to the extracted minerals (primarily potassium, also magnesium and bromide). It is also central to regional tourism. It is the lowest place on Earth and its endorheic, saline basin attracts international research in various disciplines. Additional to the Lower Jordan River, the discharge of which has been decreased to a small fraction of its original value, fresh to brackish springs are the main source of water to the rapidly dwindling Dead Sea. Although the existence of these springs has been known for decades, until recently estimates of spring water discharge into the Dead Sea were scarce. In this study, we developed a methodology of water discharge estimation for channels incising into a lacustrine bed using remotely sensed data and a single hydraulic geometry variable, water-surface width. Based on calibration of over 400 in situ measurements and simultaneous hydrometric data from aerial images, width–discharge rating curves were established, tested, and found suitable (r² = 0.92, p = 0.001) for the estimation of water discharge, with ±5% uncertainty. Furthermore, we used these relationships to estimate retrospectively the temporal changes in water discharge of seven main channels traceable in historical aerial images. The reconstructed trend reveals a major (63%) reduction in average freshwater inflow between 1990 and 2006. Our results maximize the use of water-surface width information from aerial imagery, and suggest applicability to areas experiencing rapid exposure of sea/lake bed and consequent access difficulties in in situ discharge monitoring.

Editor Z.W. Kundzewicz     

Effects of beach slope breaks on nearshore groundwater dynamics

Interactions between fresh groundwater and seawater affect significantly the nearshore pore water flow, which in turn influences the fate of nutrients and contaminants in coastal aquifers prior to discharge to the marine environment. Field investigations and numerical simulations were carried out to examine the groundwater dynamics in the intertidal zone of a carbonate sandy aquifer on the tropical island of Rarotonga, Cook Islands. The study site was featured by distinct cross‐shore slope breaks on the beach surface. Measured pore water salinities revealed different distributions under the influences of different beach profiles, inland heads, and tidal oscillations. Fresh groundwater was found to discharge around a beach slope break located in the middle area of the intertidal zone. The results indicate a strong interplay between the slope break beach morphology and tidal force in controlling the nearshore groundwater flow and solute transport. The fresh groundwater discharge location was largely determined by the beach morphology in combination with the tidal force. The nearshore groundwater flow can be very sensitive to beach slope breaks, which induce local circulation and flow instabilities. As slope breaks are a common feature of beaches around the world, these results have important, general implications for future studies of nutrients transport and transformations in nearshore aquifers and associated fluxes via submarine groundwater discharge.     

Does the Actual Drop in Dead Sea Level Reflect the Development of Water Sources Within its Drainage Basin?

As a part of Jordan’s efforts to quantify the effect of the Dead Sea level decline on the precious groundwater resources of the surrounding aquifers, the authors analyzed the historic or predevelopment inflows and outflows of the Dead Sea basin and the resulting water balance which included precipitation, evaporation, surface‐ and groundwaters. The predevelopment situation was taken as the point of departure for the sake of this study. Furthermore, the present situation was analyzed in an attempt to quantify the groundwater inflows into the Dead Sea as a result of drop in the Dead Sea level. The groundwater component and the corresponding saltwater/freshwater interface were taken as the variables to balance the levels of the sea that would have been reached without the contribution of the uncontrolled groundwater inflows as a result of the salt/freshwater interface seaward migration. The present day water balance that includes all the water diversion projects from all riparians indicates serious declines in the Dead Sea level. The effects of the present day level declines on the fresh groundwater/saltwater interface indicate that considerable amounts of groundwater are driven into the Sea as a result of the seaward migration of the freshwater/saline water interface.     

Variations in turbulent kinetic energy at a point source submarine groundwater discharge in a reef lagoon

The influence of sea level variations due to tides and wave setup on turbulent kinetic energy (TKE) was observed at a point source submarine groundwater discharge in a fringing coral reef lagoon. Tidal and wave setup variations modulated speed, TKE, TKE dissipation, and water temperature and salinity at the buoyant jet. The primary driver of jet TKE and speed variations was tides, while wave setup was a minor contributor. An inverse relationship between surface elevation and TKE was explained with an exponential equation based on sea level variations. During low tides, peak jet speeds (up to 0.3 m s^?1) and TKE per unit mass (up to 0.4 m² s^?2) were observed. As high tide approached, the jet produced minimum TKE of ~0.003 m² s^?2 and TKE dissipation ranged from 2 to 8×10^?4 m² s^?3. This demonstrated the sensitivity of the jet discharge to tides despite the small tidal range (<20 cm). Jet temperatures and salinities displayed semidiurnal oscillations with minimum salinity and temperature values during maximum discharge. Jet salinities increased throughout low tides while temperatures decreased. This pattern suggested the jet conduit was connected to a stratified cavity within the aquifer containing cool fresh water over cool salty water. As low tides progressed, jet outflow increased in salinity because of the mixing within the conduit, while lower jet temperatures suggested water coming from further or deeper in the aquifer. The presence of such a cavity has been recently confirmed by divers.     

Numerical simulation of 2D buoyant jets in ice-covered and temperature-stratified water

A two-dimensional (2D) unsteady simulation model is applied to the problem of a submerged warm water discharge into a stratified lake or reservoir with an ice cover. Numerical simulations and analyses are conducted to gain insight into large-scale convective recirculation and flow processes in a cold waterbody induced by a buoyant jet. Jet behaviors under various discharge temperatures are captured by directly modeling flow and thermal fields. Flow structures and processes are described by the simulated spatial and temporal distributions of velocity and temperature in various regions: deflection, recirculation, attachment, and impingement. Some peculiar hydrothermal and dynamic features, e.g. reversal of buoyancy due to the dilution of a warm jet by entraining cold ambient water, are identified and examined. Simulation results show that buoyancy is the most important factor controlling jet behavior and mixing processes. The inflow boundary is treated as a liquid wall from which the jet is offset. Similarity and difference in effects of boundaries perpendicular and parallel to flow, and of buoyancy on jet attachment and impingement, are discussed. Symmetric flow configuration is used to de-emphasize the Coanda effect caused by offset.     

Quantitative estimation of submarine groundwater discharge using airborne thermal infrared data acquired at two different tidal heights

Submarine groundwater discharge (SGD) plays an important role in coastal biogeochemical processes and hydrological cycles, particularly off volcanic islands in oligotrophic oceans. However, the spatial and temporal variations of SGD are still poorly understood owing to difficulty in taking rapid SGD measurements over a large scale. In this study, we used four airborne thermal infrared surveys (twice each during high and low tides) to quantify the spatiotemporal variations of SGD over the entire coast of Jeju Island, Korea. On the basis of an analytical model, we found a linear positive correlation between the thermal anomaly and squares of the groundwater discharge velocity and a negative exponential correlation between the anomaly and water depth (including tide height and bathymetry). We then derived a new equation for quantitatively estimating the SGD flow rates from thermal anomalies acquired at two different tide heights. The proposed method was validated with the measured SGD flow rates using a current meter at Gongcheonpo Beach. We believe that the method can be effectively applied for rapid estimation of SGD over coastal areas, where fresh groundwater discharge is significant, using airborne thermal infrared surveys.     

Estimation of the inland extending length of the freshwater–saltwater interface in coastal unconfined aquifers

ABSTRACT

The inland extending length of the freshwater–saltwater interface toe is useful in studies of seawater intrusion in coastal areas. The submarine fresh groundwater discharge in coastal zones is affected not only by hydraulic conductivity and hydraulic gradient of the aquifer, but also by the position of the interface. Two observation wells at different distances from the coast are required to calculate the fresh groundwater flow rate in coastal unconfined aquifers. By considering that the submarine groundwater discharge is equal to the groundwater flow rate, the length of the interface toe extending inland can be estimated when the groundwater flow is at a steady-flow state. Aquifers with horizontal and sloping confined beds and without/with unique surface vertical infiltration are considered. Examples used to illustrate the application of these methods indicate that the inland extending lengths of the interface toe in aquifers with vertical surface infiltration are much shorter than those in aquifers without vertical surface infiltration, and the length of the interface in aquifers with a horizontal confining lower bed are smaller than those in aquifers with a confining lower bed sloping towards the sea. The extent of the interface on the northwestern coast near the city of Beihai in southern Guangxi, China, on 18 January 2013 was estimated as 471–478 m.

Editor M.C. Acreman Associate editor not assigned     

Response analysis of a submarine cable under fault movement

Based on the performance of submarine cables in past earthquakes, an analytical method to determine cable performance under seabed fault movement is proposed in this paper. First, common types of earthquake damage to submarine cables are summarized, which include seabed displacement induced by fault movement, submarine landslides and seabed soil liquefaction, etc. The damage is similar to damage observed to buried pipelines following land earthquakes. The Hengchun earthquake of Dec. 26, 2006 is used as a case study. The M7.2 earthquake occurred in the South China Sea at 20:26 Beijing Time, and caused 14 international submarine cables to sever and break. The results show that the proposed method predicts damage similar to that observed in the Hengchun earthquake. Based on parametric studies of the influence of the water depth and the magnitude of the submarine earthquake, countermeasures to prevent damage to submarine cables are proposed.     

A natural tracer investigation of the hydrological regime of Spring Creek Springs, the largest submarine spring system in Florida

This work presents results from a nearly two-year monitoring of the hydrologic dynamics of the largest submarine spring system in Florida, Spring Creek Springs. During the summer of 2007 this spring system was observed to have significantly reduced flow due to persistent drought conditions. Our examination of the springs revealed that the salinity of the springs' waters had increased significantly, from 4 in 2004 to 33 in July 2007 with anomalous high radon (²²²Rn, t_1/2=3.8 days) in surface water concentrations indicating substantial saltwater intrusion into the local aquifer. During our investigation from August 2007 to May 2009 we deployed on an almost monthly basis a continuous radon-in-water measurement system and monitored the salinity fluctuations in the discharge area. To evaluate the springs' freshwater flux we developed three different models: two of them are based on water velocity measurements and either salinity or ²²²Rn in the associated surface waters as groundwater tracers. The third approach used only salinity changes within the spring area. The three models showed good agreement and the results confirmed that the hydrologic regime of the system is strongly correlated to local precipitation and water table fluctuations with higher discharges after major rain events and very low, even reverse flow during prolong droughts. High flow spring conditions were observed twice during our study, in the early spring and mid-late summer of 2008. However the freshwater spring flux during our observation period never reached that reported from a 1970s value of 4.9×10⁶ m³/day. The maximum spring flow was estimated at about 3.0×10⁶ m³/day after heavy precipitation in February-March 2008. As a result of this storm (total of 173 mm) the salinity in the spring area dropped from about 27 to 2 in only two days. The radon-in-water concentrations dramatically increased in parallel, from about 330 Bq/m³ to about 6600 Bq/m³. Such a rapid response suggests a direct connection between the deep and the surficial aquifers.     

Simulation of submarine ground water discharge to a marine estuary: Biscayne Bay, Florida

Variable density ground water flow models are rarely used to estimate submarine ground water discharge because of limitations in computer speed, data availability, and availability of a simulation tool that can minimize numerical dispersion. This paper presents an application of the SEAWAT code, which is a combined version of MODFLOW and MT3D, to estimate rates of submarine ground water discharge to a coastal marine estuary. Discharge rates were estimated for Biscayne Bay, Florida, for the period from January 1989 to September 1998 using a three-dimensional, variable density ground water flow and transport model. Hydrologic stresses in the 10-layer model include recharge, evapotranspiration, ground water withdrawals from municipal wellfields, interactions with surface water (canals in urban areas and wetlands in the Everglades), boundary fluxes, and submarine ground water discharge to Biscayne Bay. The model was calibrated by matching ground water levels in monitoring wells, baseflow to canals, and the position of the 1995 salt water intrusion line. Results suggest that fresh submarine ground water discharge to Biscayne Bay may have exceeded surface water discharge during the 1989, 1990, and 1991 dry seasons, but the average discharge for the entire simulation period was only approximately 10% of the surface water discharge to the bay. Results from the model also suggest that tidal canals intercept fresh ground water that might otherwise have discharged directly to Biscayne Bay. This application demonstrates that regional scale variable density models are potentially useful tools for estimating rates of submarine ground water discharge.     

Modelling tidal signals enhanced by a submarine spring in a coastal confined aquifer extending under the sea

Submarine springs discharge offshore groundwater from confined aquifers extending under the sea. The effects of these springs on the propagation of tidal oscillations in coastal confined aquifers are not known. This paper presents an approximate analytical solution of tidal head fluctuations in a confined aquifer with one submarine spring. The aquifer is assumed to extend in all directions infinitely. The spring is represented by a permeable round column on the seabed, which penetrates completely the impermeable layer overlying the confined aquifer. The error of the approximate solution is negligible if the distance from the spring to the coastline is much greater than the radius of the permeable column representing the spring. Through a hypothetical example, we demonstrate that it is possible to identify the spring’s location using tidal signals observed from inland wells. Tidal groundwater head fluctuations from three inland observation wells at least are needed to determine the 5 model parameters, including the location (2 parameters), the radius of the permeable column representing the spring, the diffusivity of the aquifer, and the tidal loading efficiency of the system.