DOI: 10.21952/wtr/1506943
title: Time domain-induced polarization geophysical data collected in the Rifle Floodplain, Colorado, USA

description: The IP data set was collected using a time domain-induced polarization (TDIP) method along 65 profiles of various lengths over the floodplain. Surface time domain?induced polarization (TDIP) data was collected to create the 3?D images of the complex electrical resistivity, in terms of magnitude and phase, which are associated with mineral precipitation and other lithological properties. The TDIP data was then used to estimate the spatial distribution of naturally reduced zones, which are known to be biogeochemical hotspots in floodplains having increased uranium and other redox-sensitive metal concentrations as well as high carbon contents. In the TDIP method, the transient decay of voltage is measured after current shut-off, typically in the form of an integral of decay curves over a predefined time window (so-called integral chargeability). TDIP measurements at the site were collected using the Syscal Iris Pro Switch equipment with a square-wave current injection, 50% duty cycle, and a pulse length of 2 s. The integral chargeability measurements were carried out using 20 windows during voltage decay between 240 and 1840 ms after current shut-off. Tomographic measurements were collected by deploying stainless steel electrodes with an electrode separation of 1.8 m and using a dipole-dipole ‘‘skip-2’’ and ‘‘skip-3’’ measuring protocol (i.e., for a dipoles length of 5.4 and 7.2 m, respectively). The sequence of dipole-dipole measurements was carefully arranged to (1) minimize unwanted electromagnetic coupling effects in the data, avoiding potential measurements with electrodes located inside the current dipole, (2) prevent voltage measurements using electrodes, which might be polarized due to previous current injection, and (3) increase the signal-to-noise ratio for an intended exploration depth of 8 m, i.e., the bottom of the aquifer. All measurements were collected as normal and reciprocal pairs for estimation of the data error. The IP measurements were collected with symmetric arrays (i.e., the measuring equipment placed at the center of the electrode array) with a maximum of 36 electrodes, considering that longer profiles revealed a significant increase in the normal-reciprocal misfit for the measurements of the decay curve, probably due to greater impact of electromagnetic coupling on the data. The TDIP data sets were inverted in a twodimensional domain along each transect using CRTomo, which is a smoothness-constraint inversion code based on a finite element algorithm. The resistivity and phase shift values at each pixel were then assigned at the corresponding point within the 3-D floodplain domain (the black rectangle Figure 1b) and used in the 3-D estimation. The TDIP inversion results provided the distribution of the complex resistivity, expressed in terms of its magnitude and phase-shift. The analysis of the normal-reciprocal misfit was used to estimate the relative error. We removed the measurement with smallest voltage difference (&lt;2 mV), representing about 2% variations in the data. For the inversion of TDIP measurements, chargeability values were linearly converted to frequency domain phase values (at the fundamental frequency of 0.125 Hz), by assuming a constant-phase response. This approach has been demonstrated to provide consistent results in previous studies