M.: Automated calculation of vertical pore-water flux from field temperature time series using the VFLUX method and computer program

Keywords: Heat tracing Groundwater-surface water interaction Hyporheic MATLAB Dynamic Harmonic Regression Distributed temperature sensing s u m m a r y Heat is a useful tracer for quantifying groundwater-surface water interaction, but analyzing large amounts of raw thermal data has many challenges. We present a computer program named VFLUX, written in the MATLAB computing language, for processing raw temperature time series and calculating vertical water flux in shallow sub-surface-water systems. The step-by-step workflow synthesizes several recent advancements in signal processing, and adds new techniques for calculating flux rates with large numbers of temperature records from high-resolution sensor profiles. The program includes functions for quantitatively evaluating the ideal spacing between sensor pairs, and for performing error and sensitivity analyses for the heat transport model due to thermal parameter uncertainty. The program synchronizes and resamples temperature data from multiple sensors in a vertical profile, isolates the diurnal signal from each time series and extracts its amplitude and phase angle information using Dynamic Harmonic Regression (DHR), and calculates vertical water flux rates between multiple sensor pairs using heat transport models. Flux rates are calculated every 1-to-2 h using four similar analytical methods. One or more ''sliding analysis windows'' can be used to automatically identify any number of variably spaced sensor pairs for flux calculations, which is necessary when a single vertical profile contains many sensors, such as in a high-resolution fiber-optic distributed temperature sensing (DTS) profile. We demonstrate the new method by processing two field temperature time series datasets collected using discrete temperature sensors and a high-resolution DTS profile. The analyses of field data show vertical flux rates significantly decreasing with depth at high-spatial resolution as the sensor profiles penetrate shallow, curved hyporheic flow paths, patterns which may have been obscured without the unique analytical abilities of VFLUX