Calculating flux to predict future cave radon concentrations

Matt D. Rowberry, Xavi Marti, Carlos Frontera, Marco J. Van De Wiel, Miloš Briestenský

Research output: Contribution to journalArticle

4 Citations (Scopus)
32 Downloads (Pure)

Abstract

Cave radon concentration measurements reflect the outcome of a perpetual competition which pitches flux against ventilation and radioactive decay. The mass balance equations used to model changes in radon concentration through time routinely treat flux as a constant. This mathematical simplification is acceptable as a first order approximation despite the fact that it sidesteps an intrinsic geological problem: the majority of radon entering a cavity is exhaled as a result of advection along crustal discontinuities whose motions are inhomogeneous in both time and space. In this paper the dynamic nature of flux is investigated and the results are used to predict cave radon concentration for successive iterations. The first part of our numerical modelling procedure focuses on calculating cave air flow velocity while the second part isolates flux in a mass balance equation to simulate real time dependence among the variables. It is then possible to use this information to deliver an expression for computing cave radon concentration for successive iterations. The dynamic variables in the numerical model are represented by the outer temperature, the inner temperature, and the radon concentration while the static variables are represented by the radioactive decay constant and a range of parameters related to geometry of the cavity. Input data were recorded at Driny Cave in the Little Carpathians Mountains of western Slovakia. Here the cave passages have developed along splays of the NE–SW striking Smolenice Fault and a series of transverse faults striking NW–SE. Independent experimental observations of fault slip are provided by three permanently installed mechanical extensometers. Our numerical modelling has revealed four important flux anomalies between January 2010 and August 2011. Each of these flux anomalies was preceded by conspicuous fault slip anomalies. The mathematical procedure outlined in this paper will help to improve our understanding of radon migration along crustal discontinuities and its subsequent exhalation into the atmosphere. Furthermore, as it is possible to supply the model with continuous data, future research will focus on establishing a series of underground monitoring sites with the aim of generating the first real time global radon flux maps.

Publisher Statement: NOTICE: this is the author’s version of a work that was accepted for publication in Journal of Environmental Radioactivity. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Journal of Environmental Radioactivity, [157, (2016)] DOI: 10.1016/j.jenvrad.2016.02.023

© 2016, Elsevier. Licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International http://creativecommons.org/licenses/by-nc-nd/4.0/
Original languageEnglish
Pages (from-to)16-26
Number of pages10
JournalJournal of Environmental Radioactivity
Volume157
DOIs
Publication statusPublished - 4 Mar 2016

Fingerprint

Caves
Radon
radon
cave
Fluxes
Fault slips
radioactive decay
Radioactivity
fault slip
anomaly
radioactivity
Publications
mass balance
discontinuity
cavity
Exhalation
Slovakia
Peer Review
Dilatometers
Temperature

Keywords

  • Cave radon concentration
  • Cave radon flux
  • Cave ventilation
  • Radioactive decay
  • Fault slip
  • Numerical modelling

Cite this

Calculating flux to predict future cave radon concentrations. / Rowberry, Matt D.; Marti, Xavi; Frontera, Carlos; Van De Wiel, Marco J.; Briestenský, Miloš.

In: Journal of Environmental Radioactivity, Vol. 157, 04.03.2016, p. 16-26.

Research output: Contribution to journalArticle

Rowberry, Matt D. ; Marti, Xavi ; Frontera, Carlos ; Van De Wiel, Marco J. ; Briestenský, Miloš. / Calculating flux to predict future cave radon concentrations. In: Journal of Environmental Radioactivity. 2016 ; Vol. 157. pp. 16-26.
@article{094f603bed5b4731b1c29b4e55b838bd,
title = "Calculating flux to predict future cave radon concentrations",
abstract = "Cave radon concentration measurements reflect the outcome of a perpetual competition which pitches flux against ventilation and radioactive decay. The mass balance equations used to model changes in radon concentration through time routinely treat flux as a constant. This mathematical simplification is acceptable as a first order approximation despite the fact that it sidesteps an intrinsic geological problem: the majority of radon entering a cavity is exhaled as a result of advection along crustal discontinuities whose motions are inhomogeneous in both time and space. In this paper the dynamic nature of flux is investigated and the results are used to predict cave radon concentration for successive iterations. The first part of our numerical modelling procedure focuses on calculating cave air flow velocity while the second part isolates flux in a mass balance equation to simulate real time dependence among the variables. It is then possible to use this information to deliver an expression for computing cave radon concentration for successive iterations. The dynamic variables in the numerical model are represented by the outer temperature, the inner temperature, and the radon concentration while the static variables are represented by the radioactive decay constant and a range of parameters related to geometry of the cavity. Input data were recorded at Driny Cave in the Little Carpathians Mountains of western Slovakia. Here the cave passages have developed along splays of the NE–SW striking Smolenice Fault and a series of transverse faults striking NW–SE. Independent experimental observations of fault slip are provided by three permanently installed mechanical extensometers. Our numerical modelling has revealed four important flux anomalies between January 2010 and August 2011. Each of these flux anomalies was preceded by conspicuous fault slip anomalies. The mathematical procedure outlined in this paper will help to improve our understanding of radon migration along crustal discontinuities and its subsequent exhalation into the atmosphere. Furthermore, as it is possible to supply the model with continuous data, future research will focus on establishing a series of underground monitoring sites with the aim of generating the first real time global radon flux maps.Publisher Statement: NOTICE: this is the author’s version of a work that was accepted for publication in Journal of Environmental Radioactivity. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Journal of Environmental Radioactivity, [157, (2016)] DOI: 10.1016/j.jenvrad.2016.02.023{\circledC} 2016, Elsevier. Licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International http://creativecommons.org/licenses/by-nc-nd/4.0/",
keywords = "Cave radon concentration, Cave radon flux, Cave ventilation, Radioactive decay, Fault slip, Numerical modelling",
author = "Rowberry, {Matt D.} and Xavi Marti and Carlos Frontera and {Van De Wiel}, {Marco J.} and Miloš Briestensk{\'y}",
year = "2016",
month = "3",
day = "4",
doi = "10.1016/j.jenvrad.2016.02.023",
language = "English",
volume = "157",
pages = "16--26",
journal = "Journal of Environmental Radioactivity",
issn = "0265-931X",
publisher = "Elsevier",

}

TY - JOUR

T1 - Calculating flux to predict future cave radon concentrations

AU - Rowberry, Matt D.

AU - Marti, Xavi

AU - Frontera, Carlos

AU - Van De Wiel, Marco J.

AU - Briestenský, Miloš

PY - 2016/3/4

Y1 - 2016/3/4

N2 - Cave radon concentration measurements reflect the outcome of a perpetual competition which pitches flux against ventilation and radioactive decay. The mass balance equations used to model changes in radon concentration through time routinely treat flux as a constant. This mathematical simplification is acceptable as a first order approximation despite the fact that it sidesteps an intrinsic geological problem: the majority of radon entering a cavity is exhaled as a result of advection along crustal discontinuities whose motions are inhomogeneous in both time and space. In this paper the dynamic nature of flux is investigated and the results are used to predict cave radon concentration for successive iterations. The first part of our numerical modelling procedure focuses on calculating cave air flow velocity while the second part isolates flux in a mass balance equation to simulate real time dependence among the variables. It is then possible to use this information to deliver an expression for computing cave radon concentration for successive iterations. The dynamic variables in the numerical model are represented by the outer temperature, the inner temperature, and the radon concentration while the static variables are represented by the radioactive decay constant and a range of parameters related to geometry of the cavity. Input data were recorded at Driny Cave in the Little Carpathians Mountains of western Slovakia. Here the cave passages have developed along splays of the NE–SW striking Smolenice Fault and a series of transverse faults striking NW–SE. Independent experimental observations of fault slip are provided by three permanently installed mechanical extensometers. Our numerical modelling has revealed four important flux anomalies between January 2010 and August 2011. Each of these flux anomalies was preceded by conspicuous fault slip anomalies. The mathematical procedure outlined in this paper will help to improve our understanding of radon migration along crustal discontinuities and its subsequent exhalation into the atmosphere. Furthermore, as it is possible to supply the model with continuous data, future research will focus on establishing a series of underground monitoring sites with the aim of generating the first real time global radon flux maps.Publisher Statement: NOTICE: this is the author’s version of a work that was accepted for publication in Journal of Environmental Radioactivity. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Journal of Environmental Radioactivity, [157, (2016)] DOI: 10.1016/j.jenvrad.2016.02.023© 2016, Elsevier. Licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International http://creativecommons.org/licenses/by-nc-nd/4.0/

AB - Cave radon concentration measurements reflect the outcome of a perpetual competition which pitches flux against ventilation and radioactive decay. The mass balance equations used to model changes in radon concentration through time routinely treat flux as a constant. This mathematical simplification is acceptable as a first order approximation despite the fact that it sidesteps an intrinsic geological problem: the majority of radon entering a cavity is exhaled as a result of advection along crustal discontinuities whose motions are inhomogeneous in both time and space. In this paper the dynamic nature of flux is investigated and the results are used to predict cave radon concentration for successive iterations. The first part of our numerical modelling procedure focuses on calculating cave air flow velocity while the second part isolates flux in a mass balance equation to simulate real time dependence among the variables. It is then possible to use this information to deliver an expression for computing cave radon concentration for successive iterations. The dynamic variables in the numerical model are represented by the outer temperature, the inner temperature, and the radon concentration while the static variables are represented by the radioactive decay constant and a range of parameters related to geometry of the cavity. Input data were recorded at Driny Cave in the Little Carpathians Mountains of western Slovakia. Here the cave passages have developed along splays of the NE–SW striking Smolenice Fault and a series of transverse faults striking NW–SE. Independent experimental observations of fault slip are provided by three permanently installed mechanical extensometers. Our numerical modelling has revealed four important flux anomalies between January 2010 and August 2011. Each of these flux anomalies was preceded by conspicuous fault slip anomalies. The mathematical procedure outlined in this paper will help to improve our understanding of radon migration along crustal discontinuities and its subsequent exhalation into the atmosphere. Furthermore, as it is possible to supply the model with continuous data, future research will focus on establishing a series of underground monitoring sites with the aim of generating the first real time global radon flux maps.Publisher Statement: NOTICE: this is the author’s version of a work that was accepted for publication in Journal of Environmental Radioactivity. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Journal of Environmental Radioactivity, [157, (2016)] DOI: 10.1016/j.jenvrad.2016.02.023© 2016, Elsevier. Licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International http://creativecommons.org/licenses/by-nc-nd/4.0/

KW - Cave radon concentration

KW - Cave radon flux

KW - Cave ventilation

KW - Radioactive decay

KW - Fault slip

KW - Numerical modelling

U2 - 10.1016/j.jenvrad.2016.02.023

DO - 10.1016/j.jenvrad.2016.02.023

M3 - Article

VL - 157

SP - 16

EP - 26

JO - Journal of Environmental Radioactivity

JF - Journal of Environmental Radioactivity

SN - 0265-931X

ER -