The application of Monte Carlo modelling to quantify in situ hydrogen and associated element production in the deep subsurface

The subsurface production, accumulation, and cycling of hydrogen (H2), and cogenetic elements such as sulfate (SO4 2-) and the noble gases (e.g., 4He, 40Ar) remains a critical area of research in the 21st century. Understanding how these elements generate, migrate, and accumulate is essential in terms of developing hydrogen as an alternative low-carbon energy source and as a basis for helium exploration which is urgently needed to meet global demand of this gas used in medical, industrial, and research fields. Beyond this, understanding the subsurface cycles of these compounds is key for investigating chemosynthetically-driven habitability models with relevance to the subsurface biosphere and the search for life beyond Earth. The challenge is that to evaluate each of these critical element cycles requires quantification and accurate estimates of production rates. The natural variability and intersectional nature of the critical parameters controlling production for different settings (local estimates), and for the planet as a whole (global estimates) are complex. To address this, we propose for the first time a Monte Carlo based approach which is capable of simultaneously incorporating both random and normally distributed ranges for all input parameters. This approach is capable of combining these through deterministic calculations to determine both the most probable production rates for these elements for any given system as well as defining upper and lowermost production rates as a function of probability and the most critical variables. This approach, which is applied to the Kidd Creek Observatory to demonstrate its efficacy, represents the next-generation of models which are needed to effectively incorporate the variability inherent to natural systems and to accurately model H2, 4He, 40Ar, SO4 2- production on Earth and beyond.