Groundwater Transit Times Research Articles

Steady-state flow under surface recharge through unconfined aquifer with depth-decaying hydraulic conductivity: Analytical results

This paper presents some analytical results related to groundwater flow through an unconfined horizontal aquifer induced by areal recharge on the top and a decaying hydraulic conductivity. Under steady-state conditions, and in the context of Boussinesq equation, exact analytical solutions describing drainage into one or two penetrating channels at the aquifer boundaries have been derived in closed forms. The analytical expressions have been used to investigate the existence of water divide inside the unconfined aquifer. Based on the derived analytical solutions, a closed-form analytical formula is developed for calculating groundwater transit time within Dupuit-type flow system. The analytical transit time expression is simply expressed in terms of common mathematical functions. Analytical results are verified against a numerical model and show excellent agreements. The analytical transit time is validated against a numerical one obtained from a two-dimensional model representing a rectangular unconfined aquifer. Based on the explicit closed-form solution, a fully analytical inverse problem to estimate the decay parameter is proposed using a least-square objective function. The solution of the inverse problem simply corresponds to the solution of a nonlinear algebraic equation obtained from the definition of the objective function and the form of the analytical solution.

Complex hydrology and variability of nitrogen sources in a karst watershed.

Streams draining karst areas with rapid groundwater transit times may respond relatively quickly to nitrogen reduction strategies, but the complex hydrologic network of interconnected sinkholes and springs is challenging for determining the placement and effectiveness of management practices. This study aims to inform nitrogen reduction strategies in a representative agricultural karst setting of the Chesapeake Bay watershed (Fishing Creek watershed, Pennsylvania) with known elevated nitrate contamination and a previous documented groundwater residence time of less than a decade. During baseflow conditions, streamflow did not increase with drainage area. Headwaters and the main stem lost substantial flow to sinkholes until eventually discharging along large springs downstream. Seasonal hydrologic conditions shift the flow and nitrogen load spatially among losing and gaining stream sections. A compilation of nitrogen source inputs with the geochemistry and the pattern of enrichment of δ15N and δ18O suggest that the nitrogen in streams and springs during baseflow represents a mixture of manure, fertilizer, and wastewater sources with low potential for denitrification. The pH and calcite saturation index increased along generalized flow paths from headwaters to springs and indicate shorter groundwater residence times in baseflow during the spring versus summer. Given the substantial investment in management practices, fixed monitoring sites could incorporate synoptic water sampling to properly monitor long-term progress and help inform management actions in karst watersheds. Although karst watersheds have the potential to respond to nitrogen reduction strategies due to shorter groundwater residence times, high nitrogen inputs, effectiveness of conservation practices, and release of legacy nutrients within the karst cavities could confound progress of water quality goals.

Spatial Variation in Transit Time Distributions of Groundwater Discharge to a Stream Overlying the Northern High Plains Aquifer, Nebraska, USA

AbstractGroundwater transit time distributions (TTDs) describe the spectrum of flow‐weighted apparent ages of groundwater from aquifer recharge to discharge. Regional‐scale TTDs in stream baseflow are often estimated from numerical models with limited calibration from groundwater sampling and suggest much younger groundwater discharge than has been observed by discrete age‐dating techniques. We investigate both local and regional‐scale groundwater TTDs in the Upper Middle Loup watershed (5,440 km2) overlying the High Plains Aquifer in the Nebraska Sand Hills, USA. We determined flow‐weighted apparent ages of groundwater discharging through the streambed at 88 discrete points along a 99 km groundwater‐dominated stream segment using 3H, noble gases, 14C, and groundwater flux measurements at the point‐scale (<7.6 cm diameter). Points were organized in transects across the stream width (3–10 points per transect) and transects were clustered in five sampling areas (10–610 m in stream length) located at increasing distances along the stream. Groundwater apparent ages ranged from 0 to 8,200 years and the mean groundwater transit time along the 99 km stream is >3,000 years. TTDs from upstream sampling areas were best fit by distributions with a narrow range of apparent ages, but when older groundwater from downstream sampling areas is included, the regional TTD is scale dependent and the distribution is better described by a gamma model (α ≈ 0.4) which accommodates large fractions of millennial‐aged groundwater. Observations indicate: (a) TTDs can exhibit spatial variability within a watershed and (b) watersheds can discharge larger fractions of old groundwater (>1,000 years) than commonly assumed.

EVALUATION OF AQUEOUS GASTROPOD SHELLS AS GROUNDWATER RADIOCARBON PROXIES ACROSS SPECIES AND SITES

ABSTRACTWe evaluate carbonate gastropod shells as 14C proxies for groundwater discharge at springs. Groundwater 14C is commonly used to estimate groundwater transit times, and a carbonate shell proxy would present a different way of collecting groundwater 14C data. Specifically, we test the hypothesis that in exclusively groundwater-fed spring systems, water 14C is preserved in carbonate shells at multiple sites, species, and water 14C. We first present isotopic and water temperature variability over several years at three spring sites in Utah. We then compare the 14C of contemporaneously collected water, sediment, and shells of benthic gastropods (Melanoides tuberculata, Pyrgulopsis pilsbryana, and Physella gyrina). We show that water and shell 14C activities at each site are correlated (slope = 1.00, R2 = 0.999, n = 22). These results support the hypothesis that 14C from groundwater is preserved in carbonate shells, and that aqueous gastropods a viable groundwater 14C proxy. Finally, we describe the utility and limitations of using gastropod shells as a groundwater 14C proxy.

Toward Understanding of Long‐Term Nitrogen Transport and Retention Dynamics Across German Catchments

AbstractLong‐term nitrogen (N) transport and retention dynamics across catchments are not well understood. Using a process‐based model for 89 German catchments, results across study catchments reveal that most N surplus (during 1950–2014) was removed by denitrification (mean ± standard deviation: 58 ± 15%) while the remaining fraction was mostly stored in the soil (14% ± 11%). The mean groundwater transit times in these catchments varied from 3.2 to 20.3 years. These results indicate that past N inputs could continue to affect surface and groundwater quality in the coming years. We identified four catchment groups with distinct archetypal N transport and retention dynamics, which are linked to the catchments' climate, topographic, and geological conditions. Overall, our results shed light on long‐term N dynamics in German catchments and how they are linked to catchment characteristics, emphasizing the role of long‐term N accumulation and transport for water quality management and evaluation programs.

Modeling groundwater transit time distributions and means across a Nebraska watershed: Effects of heterogeneity in the aquifer, riverbed, and recharge parameters

Groundwater transit time distributions (TTDs), the travel times through the aquifer from recharge at the water table to discharge at the surface water body, provide critical information on the timescales of hydrologic response of subsurface flow systems. We investigated the effects of spatial patterns of recharge, aquifer heterogeneity, and systematic variation in riverbed hydraulic conductivity on the mean transit times (MTTs) and TTDs of groundwater discharge to the Upper Middle Loup River (UMLR) and headwaters using a 3D-steady-state MODFLOW USG-MODPATHD3DU model. This 5436 km2 watershed overlies the High Plains Aquifer in the Sand Hills of Nebraska, USA. Modeled MTTs differed by up to three orders of magnitude from upstream to downstream in a 158 km section of the UMLR under varying recharge, aquifer, and riverbed heterogeneity scenarios. The simulated MTTs ranged from 1 to 397 years for upstream sites and 820 to 7968 years for the downstream sites. The TTDs at upstream sites were dominated by young groundwater from shallow flow paths and were sensitive to changes in riverbed hydraulic conductivity. Recharge parameterization had greater influence on the shape of the TTDs and magnitude of MTTs at the downstream sites, where much older groundwater discharged to the UMLR. Overall, spatial trends in transit times under varying model scenarios provided important information for refined conceptualization and calibration of future numerical models.

An analytical approach for urban groundwater transit time distributions accounting for the effect of stormwater infiltration system

The stormwater infiltration system (SIS) has been recognized as a feasible way to manage stormwater and restore streamflow. As a lumped representation of the flow and transport processes, the groundwater transit time distribution (TTD) is highly relevant to this restoration goal. However, the influence of the SIS on groundwater TTDs has rarely been studied. In this paper, we propose a novel analytical model for groundwater TTDs in regional urban aquifers accounting for the impact of the SIS. This analytical solution is developed for a homogeneous 2-D cross-section of an idealized unconfined aquifer following the Dupuit-Forchheimer assumption. The analytical solution is verified through the comparison with simulated TTDs using the Lagrangian particle tracking algorithm. The global sensitivity analysis reveals that the stormwater infiltration rate, the number of infiltration facilities, the size of the urban area, and the evapotranspiration are the most important factors for groundwater mean transit times (MTTs). Although groundwater discharge and MTT can be restored to the pre-urban one, groundwater TTDs will inevitably deviate from the pre-urban ones due to the alteration from diffuse recharge to infiltration-focused recharge, indicating that the urban groundwater flow regime cannot be fully restored to the pre-urban conditions. Overall, as a first-order approximation of the real-world urban groundwater system, the proposed analytical model provides decision-making support for urban water resources management and the design of the SIS.

A novel analytical model for the transit time distributions in urban groundwater systems

With the ongoing rapid urbanization across the globe, its interference with groundwater resources is critical to freshwater sustainability. The groundwater transit time distribution (TTD) lumps the flow and transport processes of a regional groundwater system and therefore characterizes the aquifer’s resilience to nonpoint-source contamination. However, the influence of large-scale urban areas on the regional groundwater TTD is not clear. This study proposed a novel analytical model for groundwater TTDs accounting for the effects of impervious urban structures. After the verification against results from particle tracking, we apply this analytical expression to investigate how the position and spatial extent of the urban area change the TTDs from the pre-urban ones. The sensitivity analysis suggests that urban areas tend to increase both the mean and the variance of groundwater transit times. Moreover, modeling results highlight the importance of the spatial relationship between the urban area and the aquifer in determining the urban groundwater TTDs. For aquifers intersected by a local urban area, mean transit time (MTT) is dominated by the horizontal extent of the urban area (i.e., the impervious area), whereas for aquifers intersected by a regional urban area, MTT is strongly controlled by the vertical extent of the urban area in addition to its horizontal size. Being computationally efficient, the proposed analytical model can aid decision-making in urban freshwater resources management and urban planning.

The 3H/3He Groundwater Age-Dating Method And Applications

Groundwater age-dating is an important tool for quantifying and managing water resources. Groundwater age is the elapsed time between recharge (at the land surface or water table) and the time when groundwater is sampled. If groundwater is sampled at the point of discharge from an aquifer, then the age represents the groundwater transit time. Groundwater that has recharged in recent decades is considered young groundwater. In many areas, the quality and quantity of young groundwater has been impacted by human activities and groundwater age-dating is useful for quantifying current and historical water and contaminant fluxes into and through aquifers. This review is focused on the tritium-helium (3H/3He) method, which is a robust and widely applied age-dating technique for young groundwater. We present the development of the 3H/3He method and practical considerations for sampling groundwater in shallow unconfined aquifers. Along the way, we highlight available resources: (1) educational software for building intuition around groundwater age-dating and selection of sampling sites and (2) software that can be used to calculate 3H/3He age from noble gas and 3H data. We also highlight strengths and potential uncertainties associated with the method. For example, while other age-dating techniques require a known historical record of tracer concentration in the atmosphere, the 3H/3He age-dating technique does not require such historical records. However, the 3H/3He method requires measurement of two tracers to produce a groundwater age estimate (“apparent age” or “tracer age”). Precise measurement of 3H and noble gases, plus careful analysis of noble gas data to calculate the tritiogenic 3He (i.e., the portion of 3He derived from decay of 3H in the aquifer) is required to calculate the groundwater apparent age. Sampling for noble gases is sometimes challenging and requires specialized sample containers and technique. We also introduce basic sampling methods in this review but highlight how practitioners should work closely with a noble gas laboratory to obtain the correct containers and assess field conditions and/or the overall feasibility of projects. Lastly, the review highlights recent applications of the 3H/3He method, including recharge rate estimation, characterization of contaminant input histories for aquifers, quantifying groundwater transit times by sampling at aquifer discharge points, and the use of isotope data to constrain and inform numerical and statistical models of groundwater and contaminant movement in the subsurface.

Geomorphological Controls on Groundwater Transit Times: A Synthetic Analysis at the Hillslope Scale

AbstractWe investigated how geomorphological structures shape Transit Time Distributions (TTDs) in shallow aquifers. Extensive three‐dimensional simulations were performed to determine the TTDs for synthetic convergent, straight, and divergent hillslopes with a constant slope. The uniform recharge applied on top of the aquifer is transferred to the receiving stream through steady‐state groundwater flows, return flows and saturation excess overland flows. Without seepage, TTDs evolve from uniform‐ to power law‐like‐ distributions depending on the average distance of the groundwater volume to the river (barycenter). Remarkably, the coefficient of variation (ratio of the standard deviation to the mean) of the TTDs scales linearly with the barycenter in agreement with a theoretical prediction based on three analytical approximations derived for specific cases. With seepage, the TTD has three separate modes corresponding to rapid saturation excess overland flows, to the intermediate flow paths ending in seepage area and to the slower flow paths going all the way to a discharge in the river. The coefficient of variation additionally depends on the extent of the seepage area. For a natural hillslope in the crystalline basement of Normandy (France), the same synthetic analysis demonstrates that the coefficient of variation is not only determined by the extent of the seepage zone but also by its structure in relation to the local and global geomorphological organization. The results suggest the possibility to assess the variability of transit times by combining geomorphological analysis, surface soil saturation observations and environmental tracers.

Hillslopes in Headwaters of Qinghai‐Tibetan Plateau as Hotspots for Subsurface Dissolved Organic Carbon Processing During Permafrost Thaw

AbstractClimate warming has accelerated thawing of northern permafrost, resulting in changes to the supply of dissolved organic carbon (DOC) to inland waters with uncertain fate. Extensive surface water and groundwater interactions occur in alpine permafrost watersheds and likely influence DOC processing differently than systems with limited interactions. Here, we quantify and characterize DOC in waters collected from eight types of water sampled across a small (25 km2) alpine (elevation 2,960–4,820 m a.s.l) watershed in the Qinghai‐Tibetan Plateau (QTP) containing variably degraded permafrost. Three types of water (thermokarst ponds, red mud gully, and seepage‐I) contained high DOC concentrations (5.2–22.6 mg L−1, n = 38), with carbon contributions predominantly from frozen soil meltwater. Spatial patterns of DOC in stream (0.3–4.8 mg L−1, n = 41), and subsurface waters (0.4–3.8 mg L−1, n = 34), all contained frozen soil meltwater carbon as constrained by δ18O and electrical conductivity, reflecting surface‐groundwater exchanges in the upper‐, middle‐, and lower stretches of the watershed. Further, patterns of increasing DOC loss in subsurface waters with decreased proportions of protein‐like organic matter and specific UV absorbance at 254 nm, suggest subsurface microbial processing. Using previously established biodegradation DOC kinetics (0.06 day−1) from the QTP, the groundwater transit time is estimated to be between 6 and 20 days based on DOC loss changes of 32% and 74% for July and September, respectively. Mass balance of DOC inputs and export fluxes demonstrate nearly half of all DOC was lost in this small watershed, indicating hillslopes are hotspots for DOC processing, with subsurface environments playing a key role.

A multi-tracer approach for understanding the functioning of heterogeneous phreatic coastal aquifers in humid tropical zones

ABSTRACT In many areas around the world, phreatic coastal aquifers are vulnerable to contamination and saltwater intrusion. This study aims to explore the hydrodynamic functioning of the coastal phreatic aquifer of the Urabá region, Colombia, based on the identification of hydrogeochemical processes and groundwater residence and transit times. Hydrodynamic data showed a fast response of the piezometric levels to precipitation with a maximum time lag of 48 h. Mean annual recharge was estimated at 19% with respect to precipitation. Another recharge mode is linked to infiltration from surface water, according to the hydrodynamic results. The main geochemical process identified is hydrolysis due to dissolved carbon dioxide, inducing silicate weathering and resulting in predominant HCO3-Ca-Mg facies. Monitoring of δ18O and δ2H compositions combined with 3H, 14C, chlorofluorocarbons (CFCs) and SF6 contents showed a fast recharge in the aquifer. Groundwater dating indicates recent waters mixed with older (>40 years) waters.

New perspectives on the use of 224Ra/228Ra and 222Rn/226Ra activity ratios in groundwater studies

The naturally occurring Ra isotopes (223Ra; T1/2 = 11.4 d, 224Ra; T1/2 = 3.66 d, 226Ra; T1/2 = 1,600 y, and 228Ra; T1/2 = 5.7 y) and Rn (222Rn; T1/2 = 3.82 d) have been widely applied as environmental tracers. The application of these radioactive tracers has mainly been restricted to the evaluation of oceanographic and land-ocean interaction processes, although in recent years their use has also been extended to the study of groundwater systems. In this context, the activity ratios of 224Ra/228Ra and 222Rn/226Ra can be instrumental in providing key information on groundwater transit times in aquifers and those processes governing groundwater discharge into the coastal sea (often referred to as Submarine Groundwater Discharge or SGD). This work evaluates the potential use of these activity ratios as proxies for investigating groundwater systems through an advective transport model that integrates the radionuclides involved in these activity ratios (224Ra, 228Ra, 226Ra, and 222Rn) and their immediate parents into a single formulation. The results provided by the transport model indicate that the main factors controlling the 224Ra/228Ra and 222Rn/226Ra activity ratios are the alpha recoil supply, the retardation factor of Ra, and the groundwater transit times. The advective transport model and the activity ratios are used to present novel applications that interrelate the disciplines of hydrogeology and coastal oceanography. The main applications include the determination of groundwater transit times and the assessment of pathways and end-members related to submarine groundwater discharge processes. These applications were tested in a Mediterranean coastal aquifer.

Recharge and residence times of groundwater in hyper arid areas: The confined aquifer of Calama, Loa River Basin, Atacama Desert, Chile

Groundwater recharge in hyper arid areas often depends on surface water infiltration and diffuse recharge of highly evaporated precipitation only contribute under favorable conditions. This happens in the Calama basin two-aquifer system, in the Central Andean area of northern Chile. A conceptual model of the groundwater system and its relationship with the Loa River is defined. We focus on the confined aquifer of the Calama basin, combining hydrodynamic, hydrogeochemical and isotopic methods. Radiocarbon (14C) activity data of dissolved inorganic carbon (DIC), in conjunction with chemical data, are applied to evaluate groundwater residence time within the confined aquifer. The Loa River recharges the Calama basin aquifers in its northeastern part, with water that has chemical and isotopic characteristics inherited from the arid environment and volcanic rocks in its upper basin. In the central and northeastern part of the confined aquifer, minor variations in chloride concentration suggest that the deep aquifer is well confined. The δ18O and δ2H values in groundwater of the confined aquifer show an increasing isotopic fractionation from the recharge area (around −10‰ δ18O) to those in the discharge area (between −8.5‰ and −8‰) in the southwestern part of the aquifer. The 14C activity continuously decreases down flow from the recharge by the Loa River. Adjusted DIC radiocarbon ages indicate a groundwater travel time between 1500 and 4000 years in the confined aquifer of Calama. Despite the limitations and uncertainties of radiocarbon in DIC to estimate groundwater transit times for the confined aquifer and considering complementary chemical and isotopic constraints, the DIC 14C provides acceptable values. The approach may be applicable in other confined aquifers in hyper-arid climates in which the formation of aquifer systems linked to river damming by geological action took place. This information is needed for sound management of the scarce groundwater resources.

Effects of Tidally Varying Salinity on Groundwater Flow and Solute Transport: Insights From Modelling an Idealized Creek Marsh Aquifer

AbstractMost existing numerical research on tide‐induced groundwater dynamics assumes a constant surface water salinity on the seaward boundary (constant salinity case). Few studies have investigated the influence of tidally varying salinity on shallow groundwater dynamics in coastal aquifers (tidal salinity case). We compiled field observations of tidally varying salinity in multiple estuaries across the eastern coast of China and a tidal creek in North Inlet‐Winyah Bay, United States. Numerical simulations were then conducted to explore the effect of tidally varying salinity on groundwater flow and salt transport in an idealized creek marsh aquifer. Results showed that the upper saline plume and classical saltwater wedge appeared in all cases, but the salinity in the saltwater wedge was diluted in the tidal salinity cases. Notably, groundwater transit times were shorter in the tidal salinity case than in the constant salinity case, especially under the creek bottom. Quantitative analyses indicated that tidally varying salinity significantly enhanced surface water‐groundwater exchange, increasing submarine groundwater discharge by 10% and the total inflow of surface water across the water‐sediment interface by 7%. As the density of groundwater differs from that of the overlying surface water, fingered saltwater flow formed in sediments under the creek bottom, leading to some small local water circulation cells. These small cells reduced groundwater transit times and almost doubled the water exchange rate. Coupling the density‐dependent flow to a simplified nitrogen reaction network revealed that the tidally varying salinity may have the potential to influence nitrogen biogeochemical transformations that modify nitrogen loads prior to discharge.

Effect of depth-dependent hydraulic conductivity and anisotropy on transit time distributions

The paper presents exact closed-form analytical solutions describing a two-dimensional profile confined groundwater flow induced by areal recharge in a laterally bounded aquifer with anisotropic permeability, which decreases with depth. The mathematical formulation of the problem for a rectangular flow domain in terms of hydraulic head and stream function are free of the limitation of the Dupuit–Forchheimer assumption. Two different types of outflow boundary conditions associated with the aquifer discharge area, Dirichlet and Neumann, are explored. Analytical solutions with respect to stream function allow examining the distribution of recharge over depth (between contouring streamlines) for a variety of flow parameter combinations. The solutions are extended to allow the groundwater transit time distribution (TTD) to be calculated. It was found that the dependence of the transit time function on hydraulic conductivity anisotropy and depth-decay coefficients may exhibit non-monotonic behavior. The mathematical models introduced in the article are accompanied by computational simulations.

Multi-tracing of recharge seasonality and contamination in groundwater: A tool for urban water resource management

In this study, sources of recharge and contamination in urban groundwater and in groundwater underneath a forest in the same aquifer were determined and compared. Data on hydro-chemical parameters and stable isotopes of water were collected in urban and forest springs in the Kharkiv region, Ukraine, over a period of 12 months. Groundwater transit time and precipitation contribution were calculated using hydrogeological data and stable isotopes of water to delineate groundwater recharge conditions. Hydro-chemical data, stable isotopes and emerging contaminants were used to trace anthropogenic groundwater recharge and approximate sewage and tap water contributions to the aquifer. The results indicated that each spring had unique isotopic signatures that could be explained by recharge conditions, groundwater residence time, and specific mixing patterns with sewage and water leaks. Elevated nitrate content, stable isotopes of nitrate, and the presence of emerging pollutants (mainly illicit drugs) in most of the urban springs confirmed mixing of urban groundwater with sewage leaks. These leaks amounted to up to 25% of total recharge and exhibited seasonal variations in some springs. Overall, the results show that urban groundwater receives variable seasonal contributions of anthropogenic components that increase the risk to the environment and human health, and reduce its usability for drinking water production. The multi-tracing approach presented can be useful for other cities worldwide that have similar problems of poor water management and inadequate sewage and water supply infrastructure.

How Aquifer Characteristics of a Watershed Affect Transit Time Distributions of Groundwater.

Article Impact Statement: The selection of lumped parameter models for groundwater transit time distributions should consider aquifer characteristics.

Using tracer-derived groundwater transit times to assess storage within a high-elevation watershed of the upper Colorado River Basin, USA

Previous watershed assessments have relied on annual baseflow to evaluate the groundwater contribution to streams. To quantify the volume of groundwater in storage, additional information such as groundwater mean transit time (MTT) is needed. This study determined the groundwater MTT in the West Fork Duchesne watershed in Utah (USA) with lumped-parameter modeling of environmental tracers (SF6, CFCs, and 3H/3He) from 21 springs. Approximately 30% of the springs exhibited an exponential transit time distribution (TTD); the remaining ~70% were best characterized by a piston-flow TTD. The flow-weighted groundwater MTT for the West Fork watershed is about 40 years with approximately 20 years in the unsaturated zone. A cumulative distribution of these ages revealed that most of the groundwater is between 30 and 50 years old, suggesting that declining recharge associated with 5–10-year droughts is less likely to have a profound effect on this watershed compared with systems with shorter MTTs. The estimated annual baseflow of West Fork stream flow based on chemical hydrograph separation is ~1.7 × 107 m3/year, a proxy for groundwater discharge. Using both MTT and groundwater discharge, the volume of mobile groundwater stored in the watershed was calculated to be ~6.5 × 108 m3, or ~20 m thickness of active groundwater storage and recharge of ~0.09 m/year (assuming porosity = 15%). Future watershed-scale assessments should evaluate groundwater MTT, in addition to annual baseflow, to quantify groundwater storage and more accurately assess watershed susceptibility to drought, groundwater extraction, and land-use change.

Application of tritium in precipitation and baseflow in Japan: a case study of groundwater transit times and storage in Hokkaido watersheds

Abstract. In this study, we demonstrate the application of tritium in precipitation and baseflow to estimate groundwater transit times and storage volumes in Hokkaido, Japan. To establish the long-term history of tritium concentration in Japanese precipitation, we used tritium data from the global network of isotopes in precipitation and from local studies in Japan. The record developed for Tokyo area precipitation was scaled for Hokkaido using tritium values for precipitation based on wine grown at Hokkaido. Then, tritium concentrations measured with high accuracy in river water from Hokkaido, Japan, were compared to this scaled precipitation record and used to estimate groundwater mean transit times (MTTs). A total of 16 river water samples in Hokkaido were collected in June, July, and October 2014 at 12 locations with altitudes between 22 and 831 m above sea level and catchment areas between 14 and 377 km2. Measured tritium concentrations were between 4.07 (± 0.07) TU and 5.29 (± 0.09) TU in June, 5.06 (± 0.09) TU in July, and between 3.75 (± 0.07) TU and 4.85 (± 0.07) TU in October. We utilised TracerLPM (Jurgens et al., 2012) for MTT estimation and introduced a Visual Basic module to automatically simulate tritium concentrations and relative errors for selected ranges of MTTs, exponential–piston ratios, and scaling factors of tritium input. Using the exponential (70 %) piston flow (30 %) model (E70 %PM), we simulated unique MTTs for seven river samples collected in six Hokkaido headwater catchments because their low tritium concentrations were no longer ambiguous. These river catchments are clustered in similar hydrogeological settings of Quaternary lava as well as Tertiary propylite formations near Sapporo city. However, nine river samples from six other catchments produced up to three possible MTT values with E70 % PM due to the interference by the tritium from the atmospheric hydrogen bomb testing 5–6 decades ago. For these catchments, we show that tritium in Japanese groundwater will reach natural levels in a decade, when one tritium measurement will be sufficient to estimate a unique MTT. Using a series of tritium measurements over the next few years with 3-year intervals will enable us to estimate the correct MTT without ambiguity in this period. These unique MTTs will allow estimation of groundwater storage volumes for water resources management during droughts and improvement of numerical model simulations. For example, the groundwater storage ranges between 0.013 and 5.07 km3 with saturated water thickness from 0.2 and 24 m. In summary, we emphasise three important points from our findings: (1) one tritium measurement is already sufficient to estimate MTTs for some Japanese catchments, (2) the hydrogeological settings control the tritium transit times of subsurface groundwater storage during baseflow, and (3) in the future, one tritium measurement will be sufficient to estimate MTTs in most Japanese watersheds.