Lithuanian Journal of Physics, Vol. 61, No. 1, pp. 53-65 (2021) © Lietuvos mokslų akademija, 2021
Received 28 August 2020; revised 6 January 2021; accepted 12 January 2021
In this study, modelled groundwater actual flow times in intermediate and deep aquifers, covered by regional scale impermeable aquitards, were compared with 4He and 81Kr age dating results. To improve the reliability of the steady state 3D groundwater flow model, the isotopic ages of deep groundwater were compared to the MODPATH modelled travel times. The highest helium values in groundwater reservoirs coincide with fault zones in the crystalline basement and sedimentary cover near Rapakivi granite massifs. Insights into isotope-geochemical anomalies of the Baltic Artesian Basin intermediate and deep groundwater support their main distribution peculiarities in the flow path towards the Baltic Sea coast lowland and seabed depression as the regional groundwater discharge area.
Keywords: deep groundwater dating, isotope-geochemistry, flow modelling
PACS: 91.67.Qr, 91.67.Rx, 92.40.Kf
A new approach of deep aquifer investigation in the Baltic Artesian Basin (BAB) was focused on using a wide complex of environmental isotopes, noble gases, groundwater chemistry, and hydraulic data via coupled analysis of groundwater flow numerical models. The main aim was to analyse crustal fracturing and hydraulic migration properties of sedimentary cover related to the fault systems and aquifer reservoir role in noble gas and fluid distribution. In this study the modelled particle travel time in intermediate and deep aquifers, which are covered by thick impermeable aquitards, was compared with 4He and 81Kr dating results. To improve the reliability of 3D groundwater flow model verification, the isotopic ages of deep groundwater were compared to the MODPATH particle tracking modelled times.
In 2016-2019 Vilnius University participated in the IAEA Research Program F33022. Following this project work plan, additional 6 Lithuani-an intermediate-deep boreholes with their depth varying from 300 up to 1011 m were studied. Stable isotopes were analysed at the Institute of Geology, Tallinn University of Technology, Estonia. In 2016-2019, at the Institute of Geosciences of Vilnius University, numerical 3D steady-state and transient groundwater models of the BAB were created and calibrated with groundwater head and hydraulic parameter data of wells. The modelled area is approximately 500 000 km2, grid size 5 × 5 km, total cell number 614823, and has 31 separated layers from the crystalline basement aquifer up to the Quaternary aquifer system. In 2013-2017, radiokrypton age was estimated for the BAB by cooperation of the international team for samples of seven boreholes of deep groundwater using the atom trap trace analysis (ATTA) method .
The BAB is a multistoried geological structure of aquifers and aquitards (Fig. 1). In western and northern marginal parts, it is inundated by the Baltic Sea that is one of the main groundwater discharge areas. The vertical sequence of groundwater bodies forms three hydrogeological units or zones - active, delayed and stagnant, separated by regional scale aquitards and formed boundary conditions for flow direction during the geotectonic development of the sedimentary basin. This separation strongly determines the isotope-geochemical composition of groundwater. Geodynamic loading processes, hydrogeochemical interactions between groundwater and rocks, groundwater flow features from the meteoric recharge area to the discharge area, and partitioning processes of isotopic composition have affected the evolution of stable oxygen-18 isotope ratio values and its locali-zation sites inside of these zones. The zone of an active water exchange extends from the inland meteoric water recharge area to the coastal submarine discharge area which is separated in the nearshore with a transitional zone, where the inversion of groundwater heads takes place. In the inland part, meteoric water infiltration predominates downwards up to depths of 400 to 450 m. Here the fresh carbonate type groundwater reaches a total dissolved solids (TDS) value up to 2-4 g/L. The zone of delayed groundwater is located at depths from 0.5 up to 1.5 km. The carbonate-chloride and sulfate type groundwater of this zone varies in TDS from 5 to 60-100 g/L. Below the 1.8 km depth, stagnant chloride type groundwater with more than 100 g/L in TDS is formed. Thus, the three main zones according to δ18O values distribution in the abovementioned groundwater zones were distinguished in the BAB.
In the active zone near surface, dominant δ18O values are from -8.2 up to -11.6%thou that are related to the modern meteoric water recharge and by features of lateral partitioning of the δ18O isotope ratio of precipitation across the region . Zones with more positive δ18O values (from -8.2 to -10.5%thou) are located mostly in the southwestern part of the Baltic Region. In Estonia, extremely depleted zones with δ18O values from -14 to -22.5%thou are observed reflecting palaeore-charge during the cold climate Pleistocene time. These differences in δ18O of groundwater are due to changes in global paleoclimatic conditions and their impact on the groundwater formation during the Pleistocene and Holocene . The δ18O values of groundwater in most aquifers in North Poland, Kaliningrad District of Russia, Lithuania and Latvia territory range from -7.7 to -13.9%thou [6, 7]. However, in the Island of Gotland (Sweden Homocline) at the same depth groundwater has significantly higher δ18O values ranging from -5.7 to -6.1%thou .
The delayed zone is located in the western part of BAB where the TDS of groundwater gradually rises from 5-7 up to 90-100 g/L. Under the impact of lithostatic load and tectonic activities, the compaction flow of this zone is moving toward the periphery part of the artesian basin. This zone includes three hydrogeochemical facies of groundwater: bicarbonate chloride calcium sodium, chloride magnesium sodium (marine) and chloride sodium. The average isotopic ratio δ18O values are -6 to -9%thou. In this zone, the δ18O values are formed like transitional media between the motion of shallow active and deep stagnant groundwater and mixing adaption.
Stagnant zone brines are under a closed system thermodynamically in complete equilibrium with sedimentary rocks of the artesian basin. Their salinity increases up to 140-300 g/L and the value of δ18O is enriched (-3 to -4.5%thou). The temperature at the top of Cambrian-Ordovician rocks varies from 98°C in the southwestern Lithuania coast area to 21°C at the Jelgava site located close to Riga city. The geothermal anomaly promotes the processes of ultrafiltration, dehydration and ion-exchange because of enrichment of deep groundwater with calcium ions, i.e. the water of chloride calcium fa-cies is formed.
Thus, three main zones according to δ18O distribution in groundwater can be distinguished: 1) a near subsurface zone with δ18O average values from -10.0 to -11.6%thou (close to local meteoric water), that are related to the modern meteoric water percolation on watersheds of highlands (eastern part of the Baltic Region); 2) zones up to 150-250 m depth with δ18O values from -10.5 to -12.5%thou (more negative values than modern mete-ogenous water); 3) a zone with more enriched values δ18O >-7%thou (more positive values than modern and Late Pleistocene time meteogenous water). The differences in δ18O of groundwater reflect changes in paleoclimatic conditions and their impact on the groundwater formation from the Early Pleistocene to Holocene.
During the sampling for radiokrypton dating, groundwater is being transferred through the degassing device, extracted gas is collected in a steel container [1, 9]. The krypton from bulk gas was separated and analysed with the noble gas mass spectrometric system at the University of Bern and the 81Kr/Kr isotope ratio was determined by using the ATTA-3 instrument in the Laboratory for Radiokrypton Dating, Argonne National Laboratory .
It is noteworthy that using the 81Kr/Kr isotope ratio for the groundwater age estimation for the BAB is possible only up to the lowest part of the intermedium zone at 1 km depth because the downward radiokrypton isotope does not have fixed tracks. So, for the deepest BAB segments, we could only use 4He based age evaluation methods.
Helium in samples was analysed in the Isotope Hydrology Laboratory of the International Atomic Energy Agency (IAEA) using the methods described by Refs.  and . A mass spectrometer MM5400 and two quadrupole mass spectrometers (QMS) along with a sample extraction system were used to determine helium quantities in groundwater sampled during the 2017 campaign. For the IAEA equipment helium was collected in copper tubes through which a constant flow of groundwater was ensured. Through clamping both ends of the copper tube, the groundwater sample is being taken without any contact with air .
The helium gas concentration of boreholes was analysed also using an INGEM-1 device at Vilnius University. The INGEM-1 analyzer measurement is based on the helium gas diffusion and absorption process through the quartz membrane induced by titanium cathode in the high voltage electric discharge in the magnetic field . Groundwater samples from the boreholes were collected into 0.33 L glass bottles. The sample bottles were sealed with rubber corks. To remove air bubbles during the sample sealing a thin metal thread was used by removing it out of the bottle while the cork was pushed finally in.
For the 4He age calculation (Eq. (1)), the average uranium (440 ppm) and thorium (1500 ppm) content in the crystalline basement at the northern part of the BAB was used . The correction (Eq. (2)) for the helium production to the calculation of accumulation rate  was made using the average rock density (2.5 g/cm3), the void ratio calculated for the porosity of reservoir aquifers (n = 0.15) and release factor ΛHe = 1, corrected in the fluid phase helium accumulation rate J and helium production J' from the rock rate. Helium age (a) is the ratio of helium content in aquifer groundwater (cm3STP/g) and corrected helium production rate J (cm3STP/ (g × a)) :
Deep and intermediate aquifer investigation of the BAB was carried out using a wide complex of environmental isotopes, noble gases, groundwater chemistry and hydraulic data via coupled analysis of groundwater flow numerical models. Helium distribution peculiarities in the BAB were studied for 40 years [16, 17]. Regional hydraulic ages were assessed by using particle tracking from the well sites up to the recharge area endpoints (Table 1, Fig. 2). In the case of a few flow paths crossing in the site, two most significant recharge areas were picked: Riga Well No. 50194 and Likėnai, Aukstaitija Well No. 21965. Thus, two particle track ages were obtained (Table 1). The average value of both hydraulic ages is presented in the diagram of dating method comparison (Fig. 3). The deep groundwater seepage velocity estimated by particle tracking for the BAB for intermediate depth up to 500 metres is 5E-3 m/a and for 1 km depth it is ~8E-4 m/a. These results are significantly lower than modelled by the BAB steady state condition with Darcy hydraulic conductivity of ~1 m/a . The modelled particle travel time in the aquifers confined with regional scale aquitards was compared with the 4He and 81Kr dating results. The groundwater reservoirs where elevated 4He values were observed coincide with the steep fault blocks that transect the basin sedimentary cover.
The radiokrypton age was estimated in deep and intermediate BAB aquifer groundwater samples collected from seven boreholes in 2013-2017 . The Lithuanian groundwater sampling took place during the campaigns in 2013 collecting noble gas for radiokrypton analysis in 3 boreholes: Ignalina, Klaipėda and Genčiai with cooperation from the Tallinn University of Technology and the Bern University. Helium was not measured in the Genčiai site sample. During the 2017 campaign helium and other noble gas were collected from 7 wells in Lithuania (in cooperation with IAEA). The results of Birstonas and Telsiai sites are not included in this study. The helium content data of Riga, Häädemeeste and Värska were previously published .
Additionally, the dissolved 4He concentration data were obtained during previous fieldwork sampling and from publications [1, 16, 17].
The cross-section of the BAB multilayered aquifers system is a setting of the regional aquitard-aquifer framework. The regional scale impermeable aquitard acts as the main separation boundary for helium migration in deep and intermediate flows. Helium produced and emitted by basement rock and aquifer matrix migrates in vertical and lateral directions toward the surface and periphery of the basin. According to these trajectories, helium accumulation in the traps results in a cumulative scheme like apparent aggregation. The traditional binary mixing model could be used for the two aquifer units separated by regional scale impermeable aquitard calculation. A leakage from beneath is not significant because permeability of the regional scale aquitard is low. The helium anomalies are influenced by its reservoir's external/internal sources and require correction.
|Well No.||Site||System||Depth, m||Particle travel time (modelled), ka recharge from south/east||4He age, ka||81Kr age, ka|
Note. D is Devonian, O-Cm is Ordovician-Cambrian.
Diﬃculties of such estimation of helium content are related to the optimization of leakage parameters and the estimation of release rates of rock radiogenic alpha particles in situ. According to our preliminary assessment based on the BAB numerical flow model, for the interface condition between the deep/intermediate zone end-members f1, f2 values may vary in intervals of 0.15 and 0.85, and for the intermediate/shallow interface zone in intervals of 0.33 and 0.66, respectively. There f1 is the helium upwelling inflow part through aquitard and f2 = (1-f1) is the in situ accumulation/depletion rate part. This analysis requires more detailed modelling with data in the vicinities of faults.
For the 4He age estimation by Eqs. (1) and (2) many diﬃculties are often related to lambda factor use because a multilayered aquifer system in many cases dramatically changes the effective porosity in regard to double porosities with fracturing and diffusion coeﬃcients of rocks.
δ2H and δ18O stable isotopes were measured in the Laboratory of Mass Spectrometry at the Department of Geology, Tallinn University of Technology. Major and trace element analyses were performed in Lithuanian and Estonian accredited laboratories.
A selective database was compiled from the data collected during the research of the BAB and previous studies. Criteria for data selection were to pick samples that are exclusive in the context of other BAB groundwater. In the case of intermediate and deep aquifers, additional criteria for groundwater dating results (4He, 81Kr, or particle travel time) were necessary. All data were divided into groups considering groundwater depth (modern-shallow, intermediate, deep), location, and structural features of the crystalline basement. Saturation indi-ces for calcite, gypsum and halite were modelled using the software Phreeqc, wateq4 database. In case the data was not available, the initial conditions were applied: temperature 10°C, electron activity (pe) 0, density 1 g/mL. All stable oxygen and deuterium isotope data are expressed in the Vienna Standard Mean Ocean Water (VSMOW) system.
Radiokrypton is meteoric in origin and has the same downward-lateral flow trajectory as a modelled particle track. This statement is applicable for the mul-tilayered aquifer systems and confined aquifer isolated from above by local aquitards where a shallow and intermediate groundwater flow direction up to 0.5-1 km depth is predominantly lateral from the recharge area toward local discharge places on land and regional discharge sites on the Baltic Sea lowland and offshore. The BAB hydraulic conductivity values of such aquitards, which confine aquifers from meteoric water leakage on a regional scale, are below 5E-7 m/d . Here the main source of emanation of 4He occurs due to the decay of uranium and thorium in minerals of the crystalline basement rock. Helium diffuses through fractures and tectonic faults to sedimentary reservoirs where it accumulates. Cases, where vertical local uplift blocks confine groundwater aquifers (offset from 100 to 600 m), could prevent lateral migration of fluid. At that boundary as such on the contact with the regional impermeable aquitard, the accumulation rate of crustal origin noble gases has increased significantly. The helium amount in the deep and at the lower part of intermediate groundwater of the BAB vary, respectively: 4.6E-8 up to 9.0E-4 cm3 STP/g . High helium accumulation rate in groundwater is observed in the periphery margin of the Estonian Homocline and Belarus-Mazurian Massif, Polish-Lithuanian Trough, Liepaja-Riga- Pskov Ridge, and with Baltic Sea Depression related Rapakivi granite massifs in the basement . The mentioned tectonic structures are bounded by Paldisk-Pskov, Middle Estonian, Liepaja-Sal-dus-Riga and West Lithuanian fault zones. These fault zones separate regions where groundwater movement rates are relatively fast from stagnant in the deepest parts of the BAB. The average velocity in the peripheral part up to the intermediate depth of the BAB varies from 0.2 to 1 m/a.
The helium age calculation shows that older groundwater may originate in an intermediate depth aquifer compared to that located deeper (Tables 2, 3). It contradicts groundwater chemical composition and dating results obtained using other techniques, therefore should be considered while interpreting helium ages. Apparent aggregated helium age in the multilayered aquifer matrix unit is the sum of all groundwater helium rates released from sources in a particular site (Tables 2, 3). During the Pleistocene period, very specific migration conditions through the sedimentary cover existed. Permafrost screening conditions, ice loading, etc. drastically changed the flow path. One of the problems is that helium bulk loss from reservoirs is practically unproved. These circumstances plausibly formed many replacements for helium anomalies and in general established a highly elevated content of helium gas in the Baltic Sea shore vicinity.
|System||He, cm3STP/g||Age, ka|
|P2||4.56E-8||Modern by 14C data|
|D3||Low helium content||Modern by 14C data|
|D2||Regional scale intermediate depth aquitard||-|
|S||Regional scale deep aquitard||-|
|Σ||Total aggregated 2.0E-3|
Note. T is Triassic, P is Permian, D is Devonian, S is Silu-rian, O-Cm is Ordovician-Cambrian.
|System||He, cm3STP/g||Age, ka|
|Cmln||Intermediate depth regional scale aquitard||-|
|Σ||Total aggregated 2.78E-4|
Note. D is Devonian, S-O is Silurian-Ordovician, O-Cm is Ordovician-Cambrian, Cmln is Cambrian Lontova, V2vr is Vendian Voronka, V2kt is Vendian Kotlin, V2gd is Vendian Gdov.
Interesting results of the U-He relationship to identify the 4He source within the Canadian Que-bec Region reveal that glaciation time fracturing of periglacial aquifers induced high helium release rates . The coupled model between the activity of 234U/238U isotope ratio, fractionation and radio-genic helium excesses was developed. It suggests a process within the aquifer to explain the 4He excesses, providing a complementary approach to the hypothesis of external sources of helium. According to the authors, the U-He relationship analysis suggests a common enhanced radiogenic 4He release process in groundwater by 234U α-recoil and helium diffusion: these release rates are between 1000 and 30,000 times higher than the local U and Th steady state production rate. Simulated 234U/238U activity ratio evolution as a function of groundwater residence time, based on the measured for the ratio maximum value, was calibrated by the 14C adjusted age 6.7 ka . In the northwestern part of BAB, the Cambrian-Vendian aquifer, which is close to Rapakivi granite massifs, the groundwater at the interface boundary with the basement groundwater 234U/238U activity ratio varies from 3 to 26 . The groundwater of Cl-HCO3-Ca-Na type has increased the TDS content, mostly 0.8-1.2 g/L, and the permafrost model for that formation was build up [6, 20]. Applying  the simulation build up trend-line for residence time, the Cambrian-Vendian groundwater age by the 234U/238U activity ratio may be evaluated as 30-80 ka. It is noteworthy that in these groundwater facies the stable oxygen-18 isotope ratio value is strongly depleted (δ18O from -20 to -22.5%thou) with a low 14C content .
Two cases of helium distribution and apparent aggregate helium age are presented in this study. The distribution of helium content (cm3STP/g) in the Klaipėda-Palanga site aquifers is the following: O-Cm 9.0E-4, D1-2 8.1E-4, D2 3.1E-4, P2 4.6E-8. The sum of helium content in the Klaipėda site is 20.2E-4 cm3STP/g (Table 2). The value of aggregate helium age could suggest more reliable dating of the deepest aquifer which in this case is O-Cm. Perhaps the real groundwater age is more likely about 1.5 Ma for this aquifer system than 0.66 Ma (by conventional helium dating for lower aquifer) considering that radiokrypton measured in the overlaying D1-2 aquifer dates groundwater up to 1.2 Ma.
|Well No.||pH||Eh, mV||δ18O, %thou||δ2H, %thou||Cl-, mg/L||SO42-, mg/L||HCO3-, mg/L||Na+, mg/L||K+, mg/L||Mg2+, mg/L||Ca2+, mg/L||TDS, mg/L||SI calcite||SI halite|
Note. SI is saturation index.
The case of Värska, in southern Estonia, presents an even more extreme case, where the deepest aquifer (of V2gd) groundwater contains more than 40 times less helium than overlaying aquifers (Table 3). The helium amount in the aquifer of Värska (cm3STP/g) is the following: V2gd 3.1E-6, V2vr 1.2E-4, O-Cm 1.1E-4, D1-2 4.8E-5. The sum of helium content in the Värska site is 2.8E-4 cm3STP/g. There the O-Cm aquifer has an in situ source for the helium rate released from the Alum Shale Formation radiogenic rocks (Table 3). A similar estimation of the age of the deepest V2gd aquifer groundwater could be deduced applying the aggregate helium age of 0.2 Ma. The conventional helium method dates groundwater to 0.002 Ma, almost 100 times younger. It is necessary to mention that this approach is not necessarily accurate. The radiokrypton age in the Värska O-Cm aquifer dates groundwater up to 0.6 Ma, that is 3 times older than the aggregate helium age. The most likely reason for that is a loss of helium which escapes to the atmosphere due to insuﬃcient aqui-tard confining capacity or ice period influences by depletion. In such cases the sum of accumulated helium is not equal to its total bulk emanation.
Groundwater residence times in intermediate and deep aquifers in general increase with depth (Fig. 3). The oldest groundwater age is obtained using a radiokrypton tracer. The particle tracking and 4He dating results show a good correlation and are significantly younger as compared with the radiokrypton method in most sites. Linear equations of the age-depth correlation are presented in Fig. 3.
The greatest differences of groundwater dating results are at 1 km depth where the 81Kr obtained age is more than 2 times greater compared to helium and particle tracking methods. These aquifers are well confined with the regional scale aquitard of Silurian rocks. The differences between ages could be explained by the reduced helium amount in the aquifer. It is most likely that helium escapes from deep aquifers through tectonic faults and diffuses upwards. The results of modelled particle tracking resident times are slightly higher than 4He and much lower than radiokrypton. The ra-diokrypton and particle track paths of shallow and intermedium aquifer are moved from the Baltic Highlands on lateral flow towards discharge sites in the Baltic Sea seabed area. The modelled aquifers are digitized without any geological structures able to retard a flow (for example, vertical faults) which can immobilize or slow groundwater motion, therefore the radiokrypton dating shows older groundwater compared to particle track.
The particle travel time and 4He age are very similar in the Lithuanian intermediate aquifers (Table 1). The intermediate depth Estonian and Latvian groundwaters show a very young 4He age compared to the radiokrypton and particle travel time. Structural features of the intermediate aquifers located in the Lithuanian-Poland Trough may result in a higher accumulation of helium than on the Saldus-Riga-Pskov Ridge (Latvia) and southern Estonian High. An uneven distribution of uranium and thorium in the crystalline basement could partly explain a low helium content in eastern Latvia, Estonia and Lithuania as well, yet this hypothesis is inconclusive because an elevated content of radioactive elements is distributed sporadically in many places [14, 16, 17]. High helium values in groundwater reservoirs coincide with the crystalline basement deformation zone near Rapakivi granite massifs [6, 14, 16, 17]. Two sites in the Lower-Middle Devonian aquifer system of East Lithuania show 3 times younger 4He age than the particle track: Likėnai (Well No. 11978) and Rokiskis (Well No. 12349).
The groundwater of deep aquifers in Lithuania and Latvia (Fig. 4(a)) are usually enriched with δ2H and δ18O isotopes in respect to current modern water in shallow aquifers (δ18O from -11.2 to -10.5%thou). Utmost enrichment is in deep brines, where δ18O varies between -5.8 to -4.4 per mil (Table 4). A deuterium excess could be explained by the evaporation factor which took place during the brine evolution. Another hypothesis could suggest that groundwater isotope fractionation occurs through ultrafiltration through thick aquitards. The deep groundwater age is from a few hundred thousand to more than a million years old that is suﬃcient for a significant stable isotope fractionation due to ultrafiltration . The evaporation trend and enrichment of stable isotopes are clear in Fig. 4(a). Stable isotope values of groundwater collected near the tectonic faults fall exactly between modern freshwater and brine, as they do at one intermediate aquifer site (Well No. 47543, Palanga). This supports the opinion that groundwater in these hydrogeological conditions is formed utilizing binary mixing, brine groundwater discharge into freshwater aquifers [21-26].
A few sites of intermediate aquifer groundwater located in northern Lithuania are slightly depleted with stable isotopes (δ18O from -11.6 to -12.2%thou). A minor deuterium excess is observable. A few hypotheses could be forwarded to explain this phenomenon. Depletion of stable isotopes could be a result of glacial time precipitation recharge and meltwater injection or caused by permafrost induced Rayleigh distillation [27-29].
Radiocarbon dating of shallow groundwater in the northern Estonia Cm-V aquifer system suggests formation time during the Late Pleistocene (17 to 33 ka) [30-32]. At sites where modern groundwater recharge takes place, groundwater is significantly younger (7 to 10 ka). Radiokryp-ton measured in the intermediate depth aquifer sites Värska and Häädemeeste dates groundwater 550 and 408 ka, respectively (Table 1). These ages correspond to the Middle Pleistocene. A relation between groundwater age and its stable isotope content is presented in Fig. 4(b). Three major groups are established: groundwater affected in the Late Pleistocene (δ18O -23 to -17%thou), the Middle (and possibly Early) Pleistocene groundwater (δ18O from -17 to -12%thou) and modern recharge (δ18O > -12%thou). Due to the fact that most of the Cm-V and O-Cm groundwater date back to the Pleistocene it is safe to assume that global glaciation and climate change processes could have affected its isotope hydrochemistry. Multiple formation scenarios are suggested to explain the origin of the Cm-V and O-Cm groundwater isotope geochemistry: glacial time precipitation recharge and meltwater injection, seawater intrusion, modern groundwater recharge and cryogenic alteration [6, 8, 19, 20, 33-37].
A new approach of deep aquifer investigation should be focused on coupled analysis using a wide complex of environmental isotopes, noble gases, groundwater chemistry, and hydraulic data of numerical flow models.
The biggest accumulation rates of helium in the BAB deep aquifers are located near the Bal-tic Sea coast and islands. In this area, intermediate and deep groundwater flow is discharged to the Baltic Sea depression and confined from meteoric water percolation. Deep and intermediate flows are separated by regional scale aquitards and helium leakages are possible only through the faults system network in the sedimentary cover. Vertical fault block dislocations in many cases may reduce or interrupt regional groundwater lateral flow paths. Depending on sedimentary cover parameters and aquifer features, the radio-genic helium from the crystalline basement has good conditions to accumulate below regional scale aquitards at reservoir traps. Due to its atmospheric origin, the radiokrypton pathway should be similar to the modelled meteoric recharged particle travel trajectory in the multilayered confined aquifer system. Yet, the difference between radiokrypton ages and particle travel times is significant. The absence of fault zones in the numerical model may be substituted for the correction of hydraulic parameters in layers to reflect the natural hydrogeological internal boundary setting more accurately. Including fault zone geometry and other parameters into the numerical model is a necessary step for further analysis of the coupled radiokrypton, 4He, and groundwater age obtained by particle tracking.
In this study, modelled groundwater actual flow times in intermediate and deep aquifers, covered by impermeable aquitards, were compared with the 4He and 81Kr dating results of the Baltic Artesian Basin groundwater. The isotope-geochemistry data helped to calibrate and preliminarily verify earlier developed steady state 3D groundwater flow models. The elevated helium content in the aquifers coincides with deformation zones in the crystalline basement and sedimentary cover close to Rapakivi granite massifs. Atmospheric in origin, the 81Kr flow path should conform to the particle travel trajectory and show similar groundwater age results. Yet, the 81Kr age and particle travel times are significantly different. The numerical model lacks structural elements such as faults, therefore correction of hydraulic parameters must be made to reflect natural conditions. The formation, age and origin of the Baltic Region groundwater anomalies are supported by the environmental isotope-geochemical data.
The work has been performed in a joint effort by Research Program F33022 between Vilnius University and IAEA. The authors are grateful for the noble gas analysis to Takuya Matsumoto, Isotope Hy-drologic Section, IAEA, Vienna, Austria.
 C. Gerber, R. Vaikmäe, W. Aeschbach, A. Babre, W. Jiang, M. Leuenberger, Z.T. Lu, R. Mokrik,P. Müller, V. Raidla, et al., Using 81Kr and noble gases to characterize and date groundwater and brines in the Baltic Artesian Basin on the one-million-year timescale, Geochim. Cosmochim. Acta 205, 187-210 (2017), https://doi.org/10.1016/j.gca.2017.01.033
 V. Juodkazis, Hydrogeological Map of the Pre-quaternary Deposits of the Soviet Baltic Republics (Ministry of Geology of the USSR, 1980).
 J. Virbulis, U. Bethers, T. Saks, J. Sennikovs, and A. Timuhins, Hydrogeological model of the Baltic Artesian Basin, Hydrogeol. J. 21, 845-862 (2013), https://doi.org/10.1007/s10040-013-0970-7
 J. Pärn, V. Raidla, R. Vaikmäe, T. Martma, J. Ivask, R. Mokrik, and K. Erg, The recharge of glacial meltwater and its influence on the geochemical evolution of groundwater in the Ordovician-Cambrian aquifer system, northern part of the Baltic Artesian Basin, Appl. Geochem. 72, 125-135 (2016), https://doi.org/10.1016/j.apgeochem.2016.07.007
 V. Raidla, Z. Kern, J. Pärn, A. Babre, K. Erg, J. Ivask, A. Kalvāns, B. Kohán, M. Lelgus, T. Martma, et al., A δ18O isoscape for the shallow groundwater in the Baltic Artesian Basin, J. Hydrol. 542, 254-267 (2016), https://doi.org/10.1016/j.jhydrol.2016.09.004
 R. Mokrik, The Paleohydrogeology of the Baltic Basin (Vilnius University Publishing House, Vilnius, 2003).
 K. Rozanski and A. Zuber, Glacial infiltration in Europe - myth or reality, Przegląd Geol. 48, 796- 803 (2000).
 R. Mokrik, The Palaeohydrogeology of the Baltic Basin. Vendian and Cambrian (Tartu University Press, 1997).
 Z.T. Lu, P. Schlosser, W.M. Smethie, N.C. Sturchio, T.P. Fischer, B.M. Kennedy, R. Purtschert, J.P. Severinghaus, D.K. Solomon, T. Tanhua, and R. Yokochi, Tracer applications of noble gas ra-dionuclides in the geosciences, Earth-Sci. Rev. 138, 196-214 (2014), https://doi.org/10.1016/j.earscirev.2013.09.002
 A. Suckow, M. Gröning, M. Jaklitsch, L.-H. Han, and P. Aggarwal, in: Proceedings of 4th Mini Conference on Noble Gases in the Hydrosphere and Natural Gas Reservoirs (Potsdam, 2007) pp. 115, https://doi.org/10.2312/GFZ.mga.050
 T. Matsumoto, D.K. Solomon, L. Araguás-Araguás, and P. Aggarwal, The IAEA's coordinated research project on estimation of groundwater recharge and discharge by using the tritium, helium-3 dating technique, Geochem. J. 51, 385-390 (2017), https://doi.org/10.2343/geochemj.2.0500
 W. Aeschbach-Hertig, Groundwater Sampling for Helium/Noble Gases (2020), https://www.iup.uni-heidelberg.de/
 L.E. Levina, V.V. Demanov, V.E. Stadnik, Yu.A. Galaskokov, V.V. Ermakov, N.N. Gar taye-nov, and I.N. Yanitskii, A measuring device for helium analysis in geological investigations, EI. VEAMS Series 10, 1-22 (1975) [in Russian].
 R. Raudsep, in: Geology and Mineral Resources of Estonia, eds. A. Raukas, A. Teedumäe (Estonian Academy Publishers, Tallinn, 1997) pp. 436, https://geoloogia.info/
 T. Torgersen and M. Stute, in: Isotope Methods for Dating Old Groundwater (International Atomic Energy Agency, Vienna, 2013) pp. 376, https://www-pub.iaea.org/MTCD/Publications/PDF/Pub1587_web.pdf
 V. Juodkazis and K. Tibar, Helium in groundwater on the northern flank of the Baltic Artesian Basin, Int. Geol. Rev. 31, 736-743(1989).
 R. Mokrik, V. Puura, T. Floden, and R. Petkevičius, Peculiarities of helium distribution in the Baltic Basin, Litosfera 6, 121-123 (2002).
 P. Mejean, D.L. Pinti, B. Ghaleb, and M. Larocque, Fracturing-induced release of radiogenic 4He and 234U into groundwater during the last de-glaciation: An alternative source to crustal helium fluxes in periglacial aquifers, J. Am. Water Resour. Assoc. 53, 5677-5689 (2017), https://doi.org/10.1002/2016WR020014
 M. Yezhova, V. Polyakov, A. Tkachenko, L. Sa vits-ki, and V. Belkina, Paleowaters of North Estonia and their influence on changes of resources and quality of fresh groundwaters of large coastal water supplies, Geology 19, 37-40 (1996).
 R. Mokrik, Pecularities of the formation of the isotopic composition of underground waters on the southern slope of the Baltic shield, Geologija 19, 16-25 (1996).
 V. Juodkazis, Regional Hydrogeology of the Baltic Region (Mokslas, Vilnius, 1989) [in Russian].
 V. Juodkazis, Regional Hydrogeology Foundations (Vilnius University Publishing House, Vilnius, 2003) [in Lithuanian].
 J. Mazeika, Regularities of Radionuclide Migration and Transformation in Lithuanian Geological Environment, Habilitation Thesis (1999).
 A. Zuzevičius, J. Mazeika, and V. Baltrūnas, A model of brackish groundwater formation in the Nemunas River Valley, Geologija 60, 63-75 (2007).
 A. Zuzevičius, The groundwater dynamics in the southern part of the Baltic Artesian Basin during the Late Pleistocene, Baltica 23, 1-12 (2010).
 M. Gregorauskas, K. Kausinis, M. Bujanauskas, V. Samalavičius, and R. Mokrik, Cenomanio-apatinės kreidos sluoksnio pozeminio vandens isteklių ir hidrocheminių anomalijų mode-linis įvertinimas, Geologija. Geografija 3, 73-79 (2017), https://doi.org/10.6001/geol-geogr.v3i2.3538 [in Lithuanian].
 I. Clark and P. Fritz, Environmental Isotopes in Hydrogeology (Lewis Publishers, Boca Raton, 1997), https://doi.org/10.1201/9781482242911
 J.R. O'Neil, Hydrogen and oxygen isotope fractionation between ice and water, J. Phys. Chem. 72(10), 3683-3684 (1968), https://doi.org/10.1021/j100856a060
 M. Lehmann and U. Siegenthaler, Equilibrium oxygen- and hydrogen-isotope fractionation between ice and water, J. Glaciol. 37, 23-26 (1991), https://doi.org/10.3189/S0022143000042751
 J. Banys, V. Juodkazis, and R. Mokrik, Regional regularities of radiocarbon distribution in groundwaters of the Baltic artesian basin, Water Resour. 6, 243-248 (1979) [in Russian].
 J. Pärn, S. Affolter, J. Ivask, S. Johnson, K. Kir-si mäe, M. Leuenberger, T. Martma, V. Raidla, S. Schloemer, H. Sepp, R. Vaikmäe, and K. Wal-rae vens, Redox zonation and organic matter oxidation in palaeogroundwater of glacial origin from the Baltic Artesian Basin, Chem. Geol. 488, 149-161 (2018), https://doi.org/10.1016/j.chemgeo.2018.04.027
 V. Raidla, K. Kirsimäe, R. Vaikmäe, E. Kaup, and T. Martma, Carbon isotope systematics of the Cambrian-Vendian aquifer system in the northern Baltic Basin: Implications to the age and evolution of groundwater, Appl. Geochem. 27, 2042-2052 (2012), https://doi.org/10.1016/j.apgeochem.2012.06.005
 A. Sterckx, J.M. Lemieux, and R. Vaikmäe, Representing glaciations and subglacial processes in hydrogeological models: A numerical investigation, Geofluids 2017, 1-12 (2017), https://doi.org/10.1155/2017/4598902
 A. Sterckx, J.-M. Lemieux, and R. Vaikmäe, Assessment of paleo-recharge under the Fenno-scandian Ice Sheet and its impact on regional groundwater flow in the northern Baltic Artesian Basin using a numerical model, Hydrogeol. J. 26, 2793-2810 (2018), https://doi.org/10.1007/s10040-018-1838-7
 A. Babre, A. Kalvāns, K. Popovs, I. Retiķe, A. Dē-liņa, R. Vaikmäe, and T. Martma, Pleistocene age paleo-groundwater inferred from water-stable isotope values in the central part of the Baltic Artesian Basin, Isot. Environ. Health Stud. 52, 706-725 (2016), https://doi.org/10.1080/10256016.2016.1168411
 A. Delina, A. Kalvans, T. Saks, U. Bethers, and V. Valdis, Highlights of Groundwater Research in the Baltic Artesian Basin (University of Latvia, Riga, 2012).
 R. Vaikmäe, L. Vallner, H.H. Loosli, P.C. Blaser, and M. Juillard-Tardent, in: Palaeowaters in Coastal Europe: Evolution of Groundwater Since the Late Pleistocene, Special Publications, eds. W.M. Edmunds, C.J. Milne, Vol. 189 (Geological Society, London, 2001) pp. 17-27, https://doi.org/10.1144/GSL.SP.2001.189.01.03
Tyrime buvo palygintas modeliuotas pozeminio vandens srauto judėjimo laikas vidutiniuose ir giliuo-se vandeninguose sluoksniuose su 4He ir 81Kr izotopais nustatytu amziumi Baltijos artezinio baseino (BAB) pozeminiame vandenyje. Izotopiniai ir hidrogeochemi-niai duomenys padėjo kalibruoti ir patikrinti anksčiau sukurtus nuostovios filtracijos 3D pozeminio vandens srauto modelius. Padidėjęs helio kiekis vandeninguose sluoksniuose sutampa su kristalinio pamato ir nuosė-dinės dangos lūzių zonomis, esančiomis salia rapakivio granito masyvų Baltijos jūros pakrantėje. Atmosferinės kilmės 81Kr izotopo trajektorija turėtų atitikti mode-liuotą dalelių judėjimo kryptį ir parodyti panasius po-zeminio vandens amziaus rezultatus. Gauti duomenys atskleidė, kad 81Kr izotopo ir modeliuotų dalelių amziai reiksmingai skiriasi. Skaitmeniniame modelyje trūksta vidinių geologinių struktūrinių elementų, tokių kaip lūzių zonos, todėl norint atspindėti natūralias sąlygas būtina koreguoti hidraulinius modelio parametrus. Remiantis izotopiniais ir hidrogeocheminiais duome-nimis, patvirtinta BAB vandens helio anomalijų kilmė, formavimosi procesai ir amzius.