Age Tracers Research Articles

Historical Trends in Ocean Heat, Carbon, Salinity, and Oxygen Simulations: Impact of a Changing Ocean Transport

AbstractExamination of historical simulations from CMIP6 models shows substantial pre‐industrial to present‐day changes in ocean heat (ΔH), salinity (ΔS), oxygen (ΔO2), dissolved inorganic carbon (ΔDIC), chlorofluorocarbon‐12 (ΔCFC12), and sulfur hexafluoride (ΔSF6). The spatial structure of the changes and the consistency among models differ among tracers: ΔDIC, ΔCFC12, and ΔSF6 all are largest near the surface, are positive throughout the thermocline with weak changes below, and there is good agreement among the models. In contrast, the largest ΔH, ΔS, and ΔO2 are not necessarily at the surface, their sign varies within the thermocline, and there are large differences among models. These differences between the two groups of tracers are linked to climate‐driven changes in the ocean transport, with this tracer “redistribution” playing a significant role in changes in ΔH, ΔS, and ΔO2 but not the other tracers. The spatial structure, and differences between models, of changes in age tracers are consistent with ΔH, ΔS, and ΔO2, supporting the hypothesis that redistribution plays a major role for these tracers. Further, the impact of the vertical displacement of isopycnals (heave) plays a major role in the differing impact of redistribution between the two groups, with this process causing insignificant changes to ΔDIC, ΔCFC12, and ΔSF6 due to their weak spatial gradients. A similar multi‐tracer analysis of observations could provide insights into the relative role of the addition and redistribution of tracers in the ocean.

Estimated in-situ carbon sequestration rates in a weathered silicate basin, southwestern Idaho, U.S.A.

Silicate weathering can induce calcite precipitation from groundwater, enabling carbon dioxide (CO2) sequestration in the critical zone (CZ), which acts as a net carbon sink with significant implications for the global carbon budget. In weathered silicates, secondary calcite dissolution accompanies precipitation-dissolution reactions, and it is unclear how calcite dissolution affects CO2 consumption in natural settings. At the Reynolds Creek Experimental Watershed - Critical Zone Observatory (RCEW-CZO), southwestern Idaho, USA, we estimate in-situ carbon sequestration rates in a semi-arid weathered silicate aquifer using hydrochemical compositions and age tracers from 6 springs and 10 wells. We delineate water-rock interactions by using observed groundwater chemistry to model open system carbon evolution, evapoconcentration in wells, silicate weathering, and formation of clays along groundwater flowpaths. We suggest carbonate precipitation under closed system conditions in deep groundwater, as calcite saturation is reached and CZ CO2 drops to just 41 % of initial concentrations in older waters. In a closed system, we estimate approximately 9 % of CZ CO2 would precipitate, indicating that on-going water-rock interactions in our weathered silicate system appear to drive continued carbon sequestration. Carbon sequestration rates via silicate weathering may help to explain a missing C sink observed at the RCEW-CZO and in other weathered silicate basins.

Design and performance of ELSA v2.0: an isochronal model for ice-sheet layer tracing

Abstract. We provide a detailed description of the ice-sheet layer age tracer Englacial Layer Simulation Architecture (ELSA) – a model that uses a straightforward method to simulate the englacial stratification of large ice sheets – as an alternative to Eulerian or Lagrangian tracer schemes. ELSA's vertical axis is time, and individual layers of accumulation are modeled explicitly and are isochronal. ELSA is not a standalone ice-sheet model but requires unidirectional coupling to another model providing ice physics and dynamics (the “host model”). Via ELSA's layer tracing, the host model’s output can be evaluated throughout the interior using ice core or radiostratigraphy data. We describe the stability and resolution dependence of this coupled modeling system using simulations of the last glacial cycle of the Greenland ice sheet using one specific host model. Key questions concern ELSA's design to maximize usability, which includes making it computationally efficient enough for ensemble runs, as well as exploring the requirements for offline forcing of ELSA with output from a range of existing ice-sheet models. ELSA is an open source and collaborative project, and this work provides the foundation for a well-documented, flexible, and easily adaptable model code to effectively force ELSA with (any) existing full ice-sheet model via a clear interface.

Strontium isotopes in the atmosphere, geosphere and hydrosphere: Developing a systematic “fingerprinting” framework of rocks and water in sedimentary basins in eastern Australia

Understanding the connection between aquifers, aquitards, and groundwater-dependant ecosystems remains a key challenge when developing a conceptual hydrogeological model. The aim of this study was to develop a systematic strontium isotope (87Sr/86Sr) fingerprinting framework of rocks and water within the sedimentary Surat and Clarence-Moreton basins (SCM basins) in eastern Australia – an area of extensive coal seam gas development and high potential for aquifer and groundwater-surface water connectivity. To do this, new groundwater samples (n = 298) were collected, analyzed and integrated with published data (n = 154) from the basins' major sedimentary, volcanic and alluvial aquifers, including the major coal seam gas target, the Walloon Coal Measures. Samples were also analyzed from rainfall (n = 2) and surface water (n = 40). In addition, rock core samples (n = 39) from exploration and stratigraphic wells were analyzed to determine the range of Sr isotope composition from host rocks. The analyses of cores demonstrate a distinct and systematic contrast in 87Sr/86Sr between different hydrogeological units. This confirms that all major hydrogeological units have a narrow range with unique 87Sr/86Sr population characteristics that are useful for guiding conceptual model development. Comparison with selected hydrochemical and groundwater age tracers (14C and 36Cl) suggests only limited changes of 87Sr/86Sr from recharge beds to the deeper parts of the basins or with a decrease in natural 14C and 36Cl tracer content along flow paths. Stream sampling during baseflow conditions confirms that 87Sr/86Sr in surface waters are similar to those of the underlying bedrock formations. We demonstrated that 87Sr/86Sr analyses of rocks and water provide a powerful hydrostratigraphic and chemostratigraphic fingerprinting framework in the SCM basins, enabling reliable assessments of plausible aquifer and groundwater-surface water interconnectivity pathways. Applied in other complex multi-aquifer sedimentary basins in Australia, and globally, a similar approach can help to constrain conceptual hydrogeological models and facilitate improved water resource management.

A Framework for Constraining Ocean Mixing Rates and Overturning Circulation from Age Tracers

Abstract The age of seawater refers to the amount of time that has elapsed since that water encountered the surface. This age measures the ventilation rate of the ocean, and the spatial distribution of age can be influenced by multiple processes, such as overturning circulation, ocean mixing, and air–sea exchange. In this work, we aim to gain new quantitative insights about how the ocean’s age tracer distribution reflects the strength of the meridional overturning circulation and diapycnal diffusivity. We propose an integral constraint that relates the age tracer flow across an isopycnal surface to the geometry of the surface. With the integral constraint, a relationship between the globally averaged effective diapycnal diffusivity and the meridional overturning strength at an arbitrary density level can be inferred from the age tracer concentration near that level. The theory is tested in a set of idealized single-basin simulations. A key insight from this study is that the age difference between regions of upwelling and downwelling, rather than any single absolute age value, is the best indicator of overturning strength. The framework has also been adapted to estimate the strength of abyssal overturning circulation in the modern North Pacific, and we demonstrate that the age field provides an estimate of the circulation strength consistent with previous studies. This framework could potentially constrain ocean circulation and mixing rates from age-like realistic tracers (e.g., radiocarbon) in both past and present climates. Significance Statement The age of seawater—the local mean time since local water from different pathways was last at the surface—is a valuable indicator of ocean circulation and the transport time scale of heat and carbon. We introduce a novel constraint that relates total age flow across a density surface to its geometry, which provides new insights into constraining ocean circulation and mixing rates from age-like realistic tracers (e.g., radiocarbon).

Carbon evolution and mixing effects on groundwater age calculations in fractured basalt, southwestern Idaho, U.S.A.

Using hydrochemical and isotopic compositions of springs and wells, we trace carbon from critical zone carbon dioxide (CO2) into groundwater of the semi-arid Reynolds Creek Experimental Watershed - Critical Zone Observatory, southwestern Idaho, USA. Dissolved inorganic carbon (DIC) concentrations, pH and stable isotope tracers of carbon for DIC (δ13CDIC), are used to show that most groundwater evolves under open system conditions, moving carbon into the groundwater and acting as a carbon sink. However, one sample (−10.94‰ δ13CDIC, 6,350 14C years before present (yrs. BP)) may have evolved under closed system conditions with a higher partial pressure of critical zone CO2 than present-day soils. By characterizing the carbon cycle, we show that (1) carbon evolution is primarily under open-system conditions, (2) shallow groundwater samples are generally less mixed and more recent (10 to 70 3H yrs. BP) than deeper groundwater samples (1,469 to 6,350 14C yrs. BP), and (3) the older portion of the groundwater may be even older than the calculated 14C ages, as indicated by the mixing of age tracers in intermediate wells. Our global conception of the deep critical zone should include carbon cycling of critical zone CO2 in old groundwater. Characterizing the deep critical zone in a semi-arid weathered silicate watershed improves our global understanding of carbon, nutrient and water cycling.

Can near-to-mid Infrared Spectral Energy Distribution Quantitatively Trace Protoplanetary Disk Evolution?

Infrared (IR) spectral energy distribution (SED) is the major tracer of protoplanetary disks. It was recently proposed to use the near-to-mid IR (or K-24) SED slope α defined between 2 and 24μm as a potential quantitative tracer of disk age. We critically examine the viability of this idea and confront it with additional statistics of IR luminosities and SED shapes. We point out that, because the statistical properties of most of the complicated physical factors involved in disk evolution are still poorly understood in a quantitative sense, the only viable way is to assume them to be random so that an idealized “average disk” can be defined, which allows the α histogram to trace its age. We confirm that the statistics of the zeroth order (luminosity), first order (slope α), and second order characteristics (concavity) of the observed K-24 SEDs indeed carry useful information upon the evolutionary processes of the “average disk”. We also stress that intrinsic diversities in K-24 SED shapes and luminosities are always large at the level of individual stars so that the application of the evolutionary path of the “average disk” to individual stars must be done with care. The data of most curves in plots are provided on GitHub (Disk-age package https://github.com/starage/disk-age/).

Deciphering the Ages of Saline Water in the Baltic Sea by Anthropogenic Radiotracers

AbstractThe slow water renewal endows the Baltic Sea a strong retention of pollutants/nutrients. Constraining water age is a practical way to depict the transport timescales for water masses and accompanying soluble substances. Although the water ages in the Baltic Sea have been resolved by 3D ocean models 20 year ago, the simulated results have not been verified. In this work, we exploited two anthropogenic radionuclides (129I and 236U) to constrain the ages of inflowing North Sea saline waters into the Baltic Sea by using the Tracer Age and Transit Time Distribution approaches. Our results indicate that the Baltic Sea has a highly stratified structure with a more diffusive bottom circulation and a more advective surface circulation. The circulation timescale for the saline water in the Baltic Sea was estimated to be >30 years from the bottom of Arkona Basin to the surface of central Baltic Sea. This work demonstrates the power of anthropogenic radiotracers in investigating circulation timescale and mixing processes in the Baltic Sea and provides the first observation‐based proof for the multi‐decadal retention of (radioactive) pollutants within the Baltic Sea.

The dominance and growth of shallow groundwater resources in continuous permafrost environments

Water is a limited resource in Arctic watersheds with continuous permafrost because freezing conditions in winter and the impermeability of permafrost limit storage and connectivity between surface water and deep groundwater. However, groundwater can still be an important source of surface water in such settings, feeding springs and large aufeis fields that are abundant in cold regions and generating runoff when precipitation is rare. Whether groundwater is sourced from suprapermafrost taliks or deeper regional aquifers will impact water availability as the Arctic continues to warm and thaw. Previous research is ambiguous about the role of deep groundwater, leading to uncertainty regarding Arctic water availability and changing water resources. We analyzed chemistry and residence times of spring, stream, and river waters in the continuous permafrost zone of Alaska, spanning the mountains to the coastal plain. Water chemistry and age tracers show that surface waters are predominately sourced from recent precipitation and have short (<50 y) subsurface residence times. Remote sensing indicates trends in the areal extent of aufeis over the last 37 y, and correlations between aufeis extent and previous year summer temperature. Together, these data indicate that surface waters in continuous permafrost regions may be impacted by short flow paths and shallow suprapermafrost aquifers that are highly sensitive to climatic and hydrologic change over annual timescales. Despite the lack of connection to regional aquifers, continued warming and permafrost thaw may promote deepening of the shallow subsurface aquifers and creation of shallow taliks, providing some resilience to Arctic freshwater ecosystems.

Mechanisms of global ocean ventilation age change during the last deglaciation

Abstract. Marine radiocarbon (14C) is widely used to trace deep-ocean circulation, providing insight into the atmosphere–ocean exchange of CO2 during the last deglaciation. Evidence shows a significantly depleted Δ14C in the glacial deep ocean, suggesting an increased ventilation age at the Last Glacial Maximum (LGM). In this study, using two transient simulations with tracers of 14C and ideal age (IAGE), we found that the oldest ventilation age is not observed at the LGM. In contrast, the models show a modestly younger ventilation age during the LGM compared to the present day. The global mean ventilation ages averaged below 1 km are approximately 800 (630) years and 930 (2000) years at the LGM and in the present day, respectively, in two simulations. This younger glacial ventilation age is mainly caused by the stronger glacial Antarctic Bottom Water (AABW) transport associated with sea ice expansion. Notably, the ocean ventilation age is significantly older predominantly in the deep Pacific during deglaciation compared to the age at the LGM, with global mean ventilation ages peaking at 1900 and 2200 years around 14–12 ka in two simulations, primarily due to the weakening of AABW transport.

TESS photometry and CAOS spectroscopy of six eclipsing binaries with Am components

In this paper, we present the results of a comprehensive study of six eclipsing binaries whose components are confirmed or suspected Am stars. By combining long-term, high-resolution Catania Astrophysical Observatory Spectropolarimeter (CAOS) spectroscopy and TESS photometry, we have been able to accurately obtain the orbital parameters of each system as well as the atmospheric parameters of its components. We performed an in-depth chemical analysis and provided chemical abundances of C, O, Na, Mg, Si, Ca, Sc, Ti, Cr, Mn, Fe, Ni, Zn, Sr, Y, Zr, and Ba. From the solution of the light and radial curves, we have determined the masses, radii, and temperatures with a good degree of accuracy. We observe apsidal motion in the eccentric system HD 216429, in which the Rossiter-McLaughlin effect is also noted. We inferred the age of our targets by fitting isochrones on the HR diagram and find that both components in each system are properly described with the same isochrone, which reinforces our results. Furthermore, dynamical and evolutionary masses, independently obtained, show an excellent agreement. According to the out-of-eclipse variability shown in their TESS light curves and their position on the Hertzsprung-Russell diagram, we claim that the stars HD 42954 (as δ Sct type) and HD 151604 (γ Dor) have a pulsating nature. Based on the chemical analysis, we corroborate that four of the systems studied here are formed by Am stars, while in the remaining ones (HD 126031 and HD 216429) only the primary component exhibits a peculiar composition. Additionally, the age distribution found in Am stars supports their suitability as age tracers in stellar populations.

Timing and Source of Recharge to the Columbia River Basalt Groundwater System in Northeastern Oregon.

Recharge to and flow within the Columbia River Basalt Group (CRBG) groundwater flow system of northeastern Oregon were characterized using isotopic, gas, and age-tracer samples from wells completed in basalt, springs, and stream base flow. Most groundwater samples were late-Pleistocene to early-Holocene; median age of well samples was 11,100 years. The relation between mean groundwater age and completed well depth across the eastern portion of the study area was similar despite differences in precipitation, topographic position, incision, thickness of the sedimentary overburden, and CRBG geologic unit. However, the lateral continuity in groundwater age was disrupted across large regional fault zones indicating these structures are substantial impediments to groundwater flow from the high-precipitation uplands to adjacent lower-precipitation and lower-elevation portions of the study area. Recharge rates calculated from the age-depth relations were <3 mm/yr and independent of the modern precipitation gradient across the study area. The age-constrained recharge rates to the CRBG groundwater system are considerably smaller than previously published estimates and highlight the uncertainty of prevailing models used to estimate recharge to the CRBG groundwater system across the Columbia Plateau in Oregon and Washington. Age tracer and isotopic evidence indicate recharge to the CRBG groundwater system is an exceedingly slow and localized process.

Temporal changes in iodine-129 and radiocesium in the Canada Basin in the Arctic Ocean between 1993 and 2020

The Arctic Ocean plays an important role in global climate and global warming through freshwater and heat exchange with subarctic waters. A better understanding of circulation time scales in the Arctic Ocean is essential to predict changes in climate and biogeochemical cycling in the Arctic Ocean. 129I and 137Cs, which have been discharged from the nuclear fuel-reprocessing facilities, have been employed to determine the time scale of the circulation in the Arctic Ocean. However, its temporal change has not been understood well. In 2017, 2019, and 2020, we measured 129I and 137Cs in the Canada Basin in the Arctic Ocean. Using our new and historical data, we discuss temporal changes in the circulation in the basin between 1993 and 2020. The tracer ages derived from the 129I/137Cs ratio indicate that the transport of the Atlantic water into the Canada Basin was accelerated in 2020. This is consistent with results of recent studies that indicated the intensified inflow of the Atlantic water into the eastern Arctic Ocean in the late 2010s, which is termed “atlantification”. Our results confirmed the “atlantification” in the Canada Basin by temporal changes in the transient tracers for the first time.

The epoch of the Milky Way’s bar formation: dynamical modelling of Mira variables in the nuclear stellar disc

ABSTRACT A key event in the history of the Milky Way is the formation of the bar. This event affects the subsequent structural and dynamical evolution of the entire Galaxy. When the bar formed, gas was likely rapidly funnelled to the centre of the Galaxy settling in a star-forming nuclear disc. The Milky Way bar formation can then be dated by considering the age distribution of the oldest stars in the formed nuclear stellar disc. In this highly obscured and crowded region, reliable age tracers are limited, but bright, high-amplitude Mira variables make useful age indicators as they follow a period–age relation. We fit dynamical models to the proper motions of a sample of Mira variables in the Milky Way’s nuclear stellar disc region. Weak evidence for inside-out growth and both radial and vertical dynamical heating with time of the nuclear stellar disc is presented, suggesting that the nuclear stellar disc is dynamically well-mixed. Furthermore, for Mira variables around a ∼350-d period, there is a clear transition from nuclear stellar disc-dominated kinematics to background bar-bulge-dominated kinematics. Using a Mira variable period–age relation calibrated in the solar neighbourhood, this suggests the nuclear stellar disc formed in a significant burst in star formation $(8\pm 1)\, \mathrm{Gyr}$ ago, although the data are also weakly consistent with a more gradual formation of the nuclear stellar disc at even earlier epochs. This implies a relatively early formation time for the Milky Way bar ($\gtrsim 8\, \mathrm{Gyr}$), which has implications for the growth and state of the young Milky Way and its subsequent history.

New insights into the flow dynamics of a deep freshwater aquifer in the semi-arid and saline Cuvelai-Etosha Basin, Northern Namibia: Results of a multi-environmental tracer study

Study regionA paleo-megafan system of the Cubango River in the northern parts of the semi-arid Cuvelai-Etosha Basin, shared by Angola and Namibia. It hosts a deep freshwater aquifer, the so-called Kalahari-Ohangwena 2 (KOH-2), with the potential to resolve the imminent regional water supply shortages. Study focusHydrogeochemical and multi-environmental tracer studies incorporating the use of age tracers 14C, 36Cl, 81Kr and 4He to determine the age of groundwater and provide insights into the flow dynamics of the KOH-2. New hydrological insights for the regionStable water isotopes and noble gas thermometry show that in a period with higher rainfall and recharge, temperatures were at least 3 – 4 °C lower than today. Several arguments led to the conclusion that younger groundwater, possibly of an age of 35,000 years, is mixed with ancient saline pore water. These include: 1) the correlation of measured 36Cl and 81Kr ratios, as well as 4He concentrations, using a binary mixing model, and 2) the substantial variation in 81Kr ages, ranging from 40,000 to 170,000 years, over relatively short distances—a phenomenon challenging to explain by advective groundwater flow equations. Consequently, the ages derived from 81Kr measurements serve as indicators of the extent of freshening and therefore describe mixing ages rather than absolute travel times.

Chronological Study of Coal‐seam Water and its Implication on Gas Production in the South Qinshui Basin

AbstractThe knowledge of the residence time of formation water is fundamental to understanding the subsurface flow and hydrological setting. To better identify the origin and evolution of coal seam water and its impact on gas storage and production, this study collected coalbed methane co‐produced water in the southeast Qinshui Basin and detected chemical and isotopic compositions, especially 36Cl and 129I concentrations. The calculated tracer ages of 129I (5.2–50.6 Ma) and 36Cl (0.13–0.76 Ma) are significantly younger than the age of coal‐bearing formation (Pennsylvanian ‐ Cisuralian), indicating freshwater recharge after coal deposition. The model that utilises 129I/I and 36Cl/Cl ratios to constrain the timing of recharge and the proportion of recharge water reveals that over 60% of pre‐anthropogenic meteoric water entered coal seams since 10 Ma and mixed with residue initial deposition water, corresponding to the basin inversion in Cenozoic. The spatial distribution of major ion concentrations reveals the primary recharge pathway for meteoric water from coal outcrops at the eastern margin to the basin center. This study demonstrates the occurrence of higher gas production rates from wells that accept water recharge in recent times and suggests the possible potential of the non‐stagnant zones for high gas production.

Lagrangian and Eulerian modelling of river plumes in the Great Barrier Reef system, Australia

In this paper, we describe methods to verify the adequacy and accuracy of Lagrangian particles from a Lagrangian model to reproduce the concentrations of a passive tracer from an Eulerian-model in river plumes. The modelling simulates plumes from two major rivers discharging in the Great Barrier Reef (GBR), Australia, under real-world scenarios. The study has been a part of a major project to aid in the protection of the GBR system from the impacts of extreme events and climate change. We employ the Regional Ocean Modelling System (ROMS) activated with its built-in Lagrangian model, and forced with wind fields from global models and recorded river volume discharges. The ROMS-Lagrangian model tracks the Lagrangian particles using the spatially interpolated velocities computed on the Eulerian ROMS three-dimensional (3D) grid. The Lagrangian particles are released in the river in proportion to the measured river volume flux. We apply a novel technique that exploits Voronoi polygon areas to convert Lagrangian particle separation into a concentration field. This facilitates comparison with the passive tracer concentrations driven by the Eulerian velocities computed on the ROMS 3D grid. We evaluate the fate of Lagrangian particles activated on a 4-km grid resolution with those of equivalent Eulerian tracer concentrations. For validation, we compare the Lagrangian particles with the Eulerian passive tracer simulated using a higher model grid resolution of 500-m, and found that the computed Lagrangian particle concentrations showed similar overall patterns compared to the passive tracer concentrations from both the coarse 4-km, and high-resolution 500-m Eulerian models. The spatial extent of the particles was in better agreement with the coarse 4-km model than with the higher resolution 500-m model. The clustering of particles resulted in structures at finer scales than the Eulerian model, over the high particle density areas that compared well against observations. The Lagrangian particles are able to capture the general river plume patterns seen in the satellite MODIS images and the satellite derived Chlorophyll-a. We further compare the online ROMS model with an offline particle tracker, OceanPARCELS and determined the similarity of results from both online and offline methods. Overall, our study demonstrates the viability to use Lagrangian particles to reproduce the Eulerian tracer properties, which can be combined with particle properties include tracer age providing to enhance the guidance on the distribution and concentration of effluent from localised flooding river events.

In vivo distribution of cerebrospinal fluid tracer in human upper spinal cord and brain stem.

BACKGROUNDIntrathecal injection is an attractive route through which drugs can be administered and directed to the spinal cord, restricted by the blood-spinal cord barrier. However, in vivo data on the distribution of cerebrospinal fluid (CSF) substances in the human spinal cord are lacking. We conducted this study to assess the enrichment of a CSF tracer in the upper cervical spinal cord and the brain stem.METHODSAfter lumbar intrathecal injection of a magnetic resonance imaging (MRI) contrast agent, gadobutrol, repeated blood samples and MRI of the upper cervical spinal cord, brain stem, and adjacent subarachnoid spaces (SAS) were obtained through 48 hours. The MRI scans were then analyzed for tracer distribution in the different regions and correlated to age, disease, and amounts of tracer in the blood to determine CSF-to-blood clearance.RESULTSThe study included 26 reference individuals and 35 patients with the dementia subtype idiopathic normal pressure hydrocephalus (iNPH). The tracer enriched all analyzed regions. Moreover, tracer enrichment in parenchyma was associated with tracer enrichment in the adjacent SAS and with CSF-to-blood clearance. Clearance from the CSF was delayed in patients with iNPH compared with younger reference patients.CONCLUSIONA CSF tracer substance administered to the lumbar thecal sac can access the parenchyma of the upper cervical spinal cord and brain stem. Since CSF-to-blood clearance is highly individual and is associated with tracer level in CSF, clearance assessment may be used to tailor intrathecal treatment regimes.FUNDINGSouth-Eastern Norway Regional Health and Østfold Hospital Trust supported the research and publication of this work.

A novel Eulerian model based on central moments to simulate age and reactivity continua interacting with mixing processes

Abstract. In geoscientific models, simulating the properties associated with particles in a continuum can serve many scientific purposes, and this has commonly been addressed using Lagrangian models. As an alternative approach, we present an Eulerian method here: diffusion–advection–reaction type partial differential equations are derived for centralized moments, which can describe the distribution of properties associated with chemicals in reaction–transport models. When the property is age, the equations for centralized moments (unlike non-central moments) do not require terms to account for aging, making this method suitable for modeling age tracers. The properties described by the distributions may also represent kinetic variables affecting reaction rates. In practical applications, continuous distributions of ages and reactivities are resolved to simulate organic matter mineralization in surficial sediments, where macrofaunal and physical mixing processes typically dominate transport. In test simulations, mixing emerged as the predominant factor shaping reactivity and age distributions. Furthermore, the applications showcase the method's aptitude for modeling continua in mixed environments while also highlighting practical considerations and challenges.

Estimation of groundwater age distributions from hydrochemistry: comparison of two metamodelling algorithms in the Heretaunga Plains aquifer system, New Zealand

Abstract. Groundwater age or residence time is important for identifying flow and contaminant pathways through groundwater systems. Typically, groundwater age and age distributions are inferred via lumped parameter models based on measured age tracer concentrations. However, due to cost and time constraints, age tracers are usually only sampled at a small percentage of the wells in a catchment. This paper describes and compares two methods to increase the number of groundwater age data points and assist with validating age distributions inferred from lumped parameter models. Two machine learning techniques with different strengths were applied to develop two independent metamodels that each aim to establish relationships between the hydrochemical parameters and the modelled groundwater age distributions in one test catchment. Ensemble medians from the best model realisations per age distribution percentile were used for comparison with the results from traditional lumped parameter models based on age tracers. Results show that both metamodelling techniques predict age distributions from hydrochemistry with good correspondence to traditional lumped parameter model (LPM)-derived age distributions. Therefore, these techniques can be used to assist with the interpretation of lumped parameter models where age tracers have been sampled, and they can also be applied to predict groundwater age distributions for wells in a similar hydrogeological regime that have hydrochemistry data available but no age tracer data.