<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing with OASIS Tables v3.0 20080202//EN" "journalpub-oasis3.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:oasis="http://docs.oasis-open.org/ns/oasis-exchange/table" dtd-version="3.0">
  <front>
    <journal-meta>
<journal-id journal-id-type="publisher">BG</journal-id>
<journal-title-group>
<journal-title>Biogeosciences</journal-title>
<abbrev-journal-title abbrev-type="publisher">BG</abbrev-journal-title>
<abbrev-journal-title abbrev-type="nlm-ta">Biogeosciences</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub">1726-4189</issn>
<publisher><publisher-name>Copernicus Publications</publisher-name>
<publisher-loc>Göttingen, Germany</publisher-loc>
</publisher>
</journal-meta>

    <article-meta>
      <article-id pub-id-type="doi">10.5194/bg-14-921-2017</article-id><title-group><article-title>Multi-frequency electrical impedance tomography as a non-invasive tool to characterize and monitor crop root systems</article-title>
      </title-group><?xmltex \runningtitle{EIT on root systems}?><?xmltex \runningauthor{M. Weigand and A. Kemna}?>
      <contrib-group>
        <contrib contrib-type="author" corresp="yes" rid="aff1">
          <name><surname>Weigand</surname><given-names>Maximilian</given-names></name>
          <email>mweigand@geo.uni-bonn.de</email>
        <ext-link>https://orcid.org/0000-0003-0510-1938</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Kemna</surname><given-names>Andreas</given-names></name>
          
        </contrib>
        <aff id="aff1"><institution>Department of Geophysics, University of Bonn, Meckenheimer
Allee 176, 53115 Bonn, Germany</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">Maximilian Weigand (mweigand@geo.uni-bonn.de)</corresp></author-notes><pub-date><day>28</day><month>February</month><year>2017</year></pub-date>
      
      <volume>14</volume>
      <issue>4</issue>
      <fpage>921</fpage><lpage>939</lpage>
      <history>
        <date date-type="received"><day>23</day><month>April</month><year>2016</year></date>
           <date date-type="rev-request"><day>23</day><month>August</month><year>2016</year></date>
           <date date-type="rev-recd"><day>27</day><month>January</month><year>2017</year></date>
           <date date-type="accepted"><day>28</day><month>January</month><year>2017</year></date>
      </history>
      <permissions>
<license license-type="open-access">
<license-p>This work is licensed under a Creative Commons Attribution 3.0 Unported License. To view a copy of this license, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/3.0/">http://creativecommons.org/licenses/by/3.0/</ext-link></license-p>
</license>
</permissions><self-uri xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017.html">This article is available from https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017.html</self-uri>
<self-uri xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017.pdf">The full text article is available as a PDF file from https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017.pdf</self-uri>


      <abstract>
    <p>A better understanding of root–soil interactions and associated processes is
essential in achieving progress in crop breeding and management, prompting
the need for high-resolution and non-destructive characterization methods. To
date, such methods are still lacking or restricted by technical constraints,
in particular the charactization and
monitoring of root growth and function in the field. A promising technique in
this respect is electrical impedance tomography (EIT), which utilizes
low-frequency (<inline-formula><mml:math id="M1" display="inline"><mml:mo>&lt;</mml:mo></mml:math></inline-formula> 1 kHz)-
electrical conduction- and
polarization properties in an imaging framework. It is well established that
cells and cell clusters exhibit an electrical polarization response in
alternating electric-current fields
due to electrical double layers which form at cell membranes. This double
layer is directly related to the electrical surface properties of the
membrane, which in turn are influenced by nutrient dynamics (fluxes and
concentrations on both sides of the membranes). Therefore, it can be assumed
that the electrical polarization properties of roots are inherently related
to ion uptake and translocation processes in the root systems. We hereby
propose broadband (mHz to hundreds of Hz) multi-frequency EIT as a
non-invasive methodological approach for the monitoring and physiological,
i.e., functional, characterization of crop root systems. The approach
combines the spatial-resolution
capability of an imaging method with the diagnostic potential of
electrical-impedance spectroscopy.
The capability of multi-frequency EIT to characterize and monitor crop root
systems was investigated in a rhizotron laboratory experiment, in which the
root system of oilseed plants was monitored in a water–filled rhizotron,
that is, in a nutrient-deprived environment. We found a low-frequency
polarization response of the root system, which enabled the successful
delineation of its spatial extension. The magnitude of the overall
polarization response decreased along with the physiological decay of the
root system due to the stress situation. Spectral polarization parameters, as
derived from a pixel-based Debye decomposition analysis of the
multi-frequency imaging results, reveal systematic changes in the spatial and
spectral electrical response of the root system. In particular, quantified
mean relaxation times (of the order of 10 ms) indicate changes in the length scales on which the polarization
processes took place in the root system, as a response to the prolonged
induced stress situation. Our results demonstrate that broadband EIT is a
capable, non-invasive method to image root system extension as well as to
monitor changes associated with the root physiological processes. Given its
applicability on both laboratory and field scales, our results suggest an
enormous potential of the method for the structural and functional imaging of root systems for various
applications. This particularly holds for the field scale, where
corresponding methods are highly desired but to date are lacking.</p>
  </abstract>
    </article-meta>
  </front>
<body>
      

<sec id="Ch1.S1" sec-type="intro">
  <title>Introduction</title>
      <p>Interest in and development of non-invasive methods for the structural and
functional characterization as well as the monitoring of root systems and the
surrounding rhizosphere have substantially increased in recent years
<xref ref-type="bibr" rid="bib1.bibx42" id="paren.1"><named-content content-type="pre">e.g.,</named-content><named-content content-type="post">and references therein</named-content></xref>. This trend is
driven mostly by the need to improve crop management and breeding techniques,
and to reduce fertilizer usage <xref ref-type="bibr" rid="bib1.bibx41" id="paren.2"><named-content content-type="pre">e.g.,</named-content></xref>. In this
context, various non-invasive methods for the investigation and
characterization of crop root systems have been proposed <xref ref-type="bibr" rid="bib1.bibx65" id="paren.3"><named-content content-type="pre">for a
comprehensive overview of current methods, both for laboratory and field
studies, see</named-content></xref>. These methods include light transmission
tomography <xref ref-type="bibr" rid="bib1.bibx72" id="paren.4"><named-content content-type="pre">e.g.,</named-content></xref>, X-ray computer tomography
<xref ref-type="bibr" rid="bib1.bibx39 bib1.bibx72" id="paren.5"><named-content content-type="pre">e.g.,</named-content></xref>, neutron radiography
<xref ref-type="bibr" rid="bib1.bibx99" id="paren.6"><named-content content-type="pre">e.g.,</named-content></xref>, magnetic resonance imaging
<xref ref-type="bibr" rid="bib1.bibx68" id="paren.7"><named-content content-type="pre">e.g.,</named-content><named-content content-type="post">and references therein</named-content></xref>, electrical
resistivity tomography (ERT) <xref ref-type="bibr" rid="bib1.bibx65" id="paren.8"><named-content content-type="pre">e.g.,</named-content></xref>, electrical
capacitance measurements, and electrical impedance spectroscopy (EIS)
<xref ref-type="bibr" rid="bib1.bibx65 bib1.bibx6" id="paren.9"><named-content content-type="pre">see</named-content><named-content content-type="post">and references
therein</named-content></xref>. However, most of these
methods cannot, or only under special circumstances, be used at the field scale, or they lack sensitivity to structural or physiological features of the rhizosphere <xref ref-type="bibr" rid="bib1.bibx65" id="paren.10"><named-content content-type="pre">e.g.,</named-content></xref>.</p>
      <p>Electrical methods, including both tomographic and spectroscopic approaches,
are gaining importance for root research due to their universal applicability at different scales and the
recognized potential to provide pertinent information on root systems via
their electrical properties. Advances in measurement accuracy
<xref ref-type="bibr" rid="bib1.bibx102" id="paren.11"/> and large-scale deployments
<xref ref-type="bibr" rid="bib1.bibx46 bib1.bibx62" id="paren.12"><named-content content-type="pre">e.g.,</named-content></xref> allow imaging studies
with high spatial and temporal resolution at both laboratory and field scales
<xref ref-type="bibr" rid="bib1.bibx49 bib1.bibx85" id="paren.13"><named-content content-type="pre">see, e.g.,</named-content></xref>.</p>
      <p>Electrical resistance measurements on root systems have been related to root
age <xref ref-type="bibr" rid="bib1.bibx7" id="paren.14"/>, to absorbing root surfaces of trees
<xref ref-type="bibr" rid="bib1.bibx8 bib1.bibx20" id="paren.15"/>, and to surface area in contact
with the ambient solution <xref ref-type="bibr" rid="bib1.bibx18" id="paren.16"/>. The measured resistances
are usually interpreted by means of equivalent electrical circuit models of
the root–soil continuum, and relations to biological properties are analyzed
in terms of the circuit model parameters. Electrical imaging applications on
crop root systems, however, are relatively rare. ERT has been used to map
root zones
<xref ref-type="bibr" rid="bib1.bibx1 bib1.bibx4 bib1.bibx5 bib1.bibx2 bib1.bibx81" id="paren.17"/>
and to monitor water content in maize fields
<xref ref-type="bibr" rid="bib1.bibx86 bib1.bibx12" id="paren.18"/> and under an apple orchard
<xref ref-type="bibr" rid="bib1.bibx16" id="paren.19"/>. <xref ref-type="bibr" rid="bib1.bibx98" id="text.20"/> showed that ERT can be
used in the field to indirectly phenotype root systems by monitoring water
content distributions over time.</p>
      <p>As pointed out by <xref ref-type="bibr" rid="bib1.bibx91" id="text.21"/>, resistance methods for root
characterization suffer from an inherent ambiguity of effective conductivity
(or resistivity), making interpretation difficult. Polarization properties,
on the other hand, provide valuable additional information, in particular if
their spectral variation is explored. In geophysics, corresponding
measurement approaches are referred to as induced polarization (IP) or
spectral induced polarization
(SIP) methods, since the polarization is provoked by an impressed electric
field. A wide range of studies have investigated electrical polarization
properties of plant root systems, mostly in terms of capacitances, using
alternating-current measurements at a particular frequency
<xref ref-type="bibr" rid="bib1.bibx92 bib1.bibx21 bib1.bibx32 bib1.bibx27 bib1.bibx9 bib1.bibx29" id="paren.22"><named-content content-type="pre">e.g.,</named-content></xref>.
Correlations of varying strength have been found between measured
capacitances and root (dry) mass, root surface, and various attributes
associated with physiological processes such as root development. For
example, <xref ref-type="bibr" rid="bib1.bibx33" id="text.23"/> used an improved measurement setup to
investigate the relation of electrical capacitances to root mass, root
surface area, and root length in soil experiments. For an overview of studies
using electrical capacitance measurements on root systems, we refer the
reader to <xref ref-type="bibr" rid="bib1.bibx54" id="text.24"/>. While in the above-mentioned studies
single-frequency capacitance measurements were used, more recent studies also
focused on the analysis of spectral measurements covering a broad frequency
range, in terms of both capacitances <xref ref-type="bibr" rid="bib1.bibx70" id="paren.25"/> and
impedances
<xref ref-type="bibr" rid="bib1.bibx70 bib1.bibx19 bib1.bibx100 bib1.bibx25 bib1.bibx77" id="paren.26"/>.</p>
      <p>Research has also been conducted on electrical properties at the cellular
scale, including electrical surface properties of cell membranes, also called
plasma membranes <xref ref-type="bibr" rid="bib1.bibx51 bib1.bibx94" id="paren.27"><named-content content-type="pre">e.g.,</named-content></xref>. An
electrical double layer (EDL) forms at an electrically charged surface in
contact with an electrolyte <xref ref-type="bibr" rid="bib1.bibx63" id="paren.28"><named-content content-type="pre">e.g.,</named-content></xref>. This EDL gives
rise to electrical polarizability <xref ref-type="bibr" rid="bib1.bibx64" id="paren.29"><named-content content-type="pre">e.g.,</named-content></xref>,
which can be measured with EIS or electrical impedance tomography (EIT).
Accordingly, variations in the EDL characteristics related to structural or
functional changes in the root system should manifest in electrical impedance
measurements. According to <xref ref-type="bibr" rid="bib1.bibx53" id="text.30"/> cell walls can be
assumed to be near ionic equilibrium with the surrounding electrolyte and
thus do not contribute to the formation of EDLs in biomaterial.</p>
      <p>Imaging of IP or SIP parameters has, so far, to our knowledge, not been
applied to the field of root research. However, various applications in
near-surface petro- and biogeophysics have been successful. For example,
spectral (i.e., multi-frequency) EIT was used to map subsurface hydrocarbon
contamination at an industrial site <xref ref-type="bibr" rid="bib1.bibx35" id="paren.31"/> and to
monitor uranium precipitation induced by bacterial injections within the
frame of contaminated site remediation <xref ref-type="bibr" rid="bib1.bibx37" id="paren.32"/> –
both studies demonstrating the field-scale applicability of the method for
subsurface (bio)geochemical characterization. <xref ref-type="bibr" rid="bib1.bibx66" id="text.33"/>
applied EIT to investigate fungus infestation of trees; however, in the
imaging they did not take the spectral variation into account.</p>
      <p>In the present work, we propose broadband (mHz – kHz) multi-frequency EIT as an
imaging tool for the physiological, i.e., functional, characterization of crop
root systems. This novel approach for functional root imaging combines the
spatial resolution benefits of EIT with the diagnostic capability of EIS, and
builds upon instrumentation and processing tools that have been developed in
recent years.  Analogous to the now widely accepted interpretation of SIP
signatures of soils and rocks in terms of textural and mineral surface
characteristics, we hypothesize that the SIP response of crop root systems,
which is imaged with the proposed methodology, is directly related to
physicochemical processes in the vicinity of electrical double layers forming
in association with root physiological activity (e.g., nutrient uptake) at
specific scales of the root system.</p>
      <p>Besides the spatial delineation and monitoring of active root zones in terms
of polarization magnitude, we aim at the analysis of the imaged SIP response
in terms of relaxation times, which provides information on the spatial
length scale at which the underlying processes occur. Relaxation times are
determined using the Debye decomposition scheme, a phenomenological model
that can describe a wide variety of SIP signatures
<xref ref-type="bibr" rid="bib1.bibx69 bib1.bibx96" id="paren.34"><named-content content-type="pre">e.g.,</named-content></xref>. A similar procedure
to analyse SIP signatures is also proposed by
<xref ref-type="bibr" rid="bib1.bibx70" id="text.35"/> for the analysis of SIP signatures
measured on root systems.</p>
      <p>To demonstrate the proposed methodology, we conducted a laboratory experiment
on oilseed plants grown in hydroponic conditions. The plants were placed in a rhizotron
container filled with tap water and monitored using multi-frequency EIT in
the course of prolonged nutrient deficiency. The recovered spectral
electrical signatures at various time steps were analyzed with regard to
total polarization strength and dominant relaxation timescales, and
qualitatively related to the macroscopic reaction of the root system to the
induced stress situation.</p>
      <p>The next section shortly reviews electrical measurements on, and the underlying polarization properties of root systems. Then, the geophysical
methods used in the presented study are described, followed by the
experimental setup and data acquisition/processing steps. The last two
sections present the results and discuss methodological and biological
aspects of the experiment.</p>
</sec>
<sec id="Ch1.S2">
  <?xmltex \opttitle{Electrical properties and measurements of\hack{\break} root systems}?><title>Electrical properties and measurements of<?xmltex \hack{\break}?> root systems</title>
      <p>This section develops our working hypotheses regarding the electrical
polarization of crop root systems. A more detailed description of the EDL is
given and linked to the measurement methodology. Moreover, we shortly review
previous works on the small-scale (cells and cell suspensions) polarization of
biomatter and the approaches used to analyze polarization measurements on
whole root systems.</p>
<sec id="Ch1.S2.SS1">
  <title>Electrical double layer polarization</title>
      <p>Electrical conduction properties of soils are primarily determined by
electrolytic soil water conductivity, i.e., ion concentration and mobility,
and the interface conduction processes at water–mineral interfaces. Electrical
polarization properties originate mainly in ion accumulation processes in
constrictions of the pore network and at water–mineral interfaces. If
surfaces are electrically charged and in contact with an electrolyte such as those found at mineral grain surfaces or cell membranes, electrical double
layers (EDLs) form, which comprise the so-called Stern layer of bound
counterions and the so-called diffusive layer. The latter forms in
the equilibrium of electromigrative and diffusive ion fluxes, and is characterized
by ion concentration gradients. The EDL is affected by external electric
fields, manifesting an induced polarization (IP), and takes a finite time
(relaxation time) to reach equilibrium again once an impressed external field
is turned off <xref ref-type="bibr" rid="bib1.bibx63" id="paren.36"><named-content content-type="pre">e.g.,</named-content></xref>. Models of both Stern layer
polarization (the build-up of counterion concentration gradients in the Stern
layer in the direction of the external electric field)
<xref ref-type="bibr" rid="bib1.bibx84 bib1.bibx59" id="paren.37"><named-content content-type="pre">e.g.,</named-content></xref> and diffuse layer
polarization (the build-up of counterion and coion concentration gradients in the
diffuse layer in the direction of the external electric field)
<xref ref-type="bibr" rid="bib1.bibx30 bib1.bibx34" id="paren.38"><named-content content-type="pre">e.g.,</named-content></xref> have been developed, as
well as models encompassing both the Stern layer and the diffuse layer
<xref ref-type="bibr" rid="bib1.bibx64 bib1.bibx75" id="paren.39"><named-content content-type="pre">e.g.,</named-content></xref>. In a porous
system such as soil, diffuse layer polarization is also referred to as
membrane polarization since the resultant ion concentration gradients, for
instance along a pore constriction, have an effect similar to an
ion-selective membrane <xref ref-type="bibr" rid="bib1.bibx17" id="paren.40"><named-content content-type="pre">e.g.,</named-content></xref>. Strength and
relaxation behavior of EDL polarization are influenced
by background ion concentration in the pore water and surface charge density,
among other factors
<xref ref-type="bibr" rid="bib1.bibx63" id="paren.41"><named-content content-type="pre">e.g.,</named-content></xref>. Importantly, the relaxation time relates
to the spatial length scale of the polarization process and the ionic
diffusion coefficient in the EDL, which may be different for Stern layer and
diffuse layer polarization <xref ref-type="bibr" rid="bib1.bibx64" id="paren.42"><named-content content-type="pre">e.g.,</named-content></xref>. The
relationship between relaxation time and characteristic length scale for
induced polarization in soils and sediments has been investigated in many
studies
<xref ref-type="bibr" rid="bib1.bibx89 bib1.bibx15 bib1.bibx55 bib1.bibx78 bib1.bibx80" id="paren.43"><named-content content-type="pre">e.g.,</named-content></xref>.</p>
</sec>
<sec id="Ch1.S2.SS2">
  <title>Electrical measurements</title>
      <p>Electrical methods measure the conduction and polarization properties of a
medium. In the frequency domain, the measured quantity is the complex-valued
impedance, with the real (ohmic) part accounting for conduction, and the
imaginary part accounting for polarization (capacitive) effects.</p>
      <p>The electrical impedance, <inline-formula><mml:math id="M2" display="inline"><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover></mml:math></inline-formula>, at a measurement (angular) frequency
<inline-formula><mml:math id="M3" display="inline"><mml:mi mathvariant="italic">ω</mml:mi></mml:math></inline-formula> is defined as the ratio of the complex voltage <inline-formula><mml:math id="M4" display="inline"><mml:mover accent="true"><mml:mi>U</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover></mml:math></inline-formula> to the
current <inline-formula><mml:math id="M5" display="inline"><mml:mover accent="true"><mml:mi>I</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover></mml:math></inline-formula>, and can be represented by a real part <inline-formula><mml:math id="M6" display="inline"><mml:mrow><mml:msup><mml:mi>Z</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> and an
imaginary part <inline-formula><mml:math id="M7" display="inline"><mml:mrow><mml:msup><mml:mi>Z</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>:
            <disp-formula id="Ch1.E1" content-type="numbered"><mml:math id="M8" display="block"><mml:mrow><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mover accent="true"><mml:mi>U</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:mover accent="true"><mml:mi>I</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:msup><mml:mi>Z</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo><mml:mo>+</mml:mo><mml:mi>j</mml:mi><mml:msup><mml:mi>Z</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msup><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
          with <inline-formula><mml:math id="M9" display="inline"><mml:mi>j</mml:mi></mml:math></inline-formula> denoting the imaginary unit. The inverse of the impedance is the
admittance <inline-formula><mml:math id="M10" display="inline"><mml:mrow><mml:mover accent="true"><mml:mi>Y</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>/</mml:mo><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:msup><mml:mi>Y</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo><mml:mo>+</mml:mo><mml:mi>j</mml:mi><mml:msup><mml:mi>Y</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msup><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula>.
Impedances, or admittances, can be translated to effective material properties
by means of a (real-valued) geometrical factor <inline-formula><mml:math id="M11" display="inline"><mml:mi>K</mml:mi></mml:math></inline-formula>, which takes into account
the geometric dimensions of the measurement (in particular electrode
positions):

                <disp-formula specific-use="align" content-type="numbered"><mml:math id="M12" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E2"><mml:mtd/><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mtext>a</mml:mtext></mml:msub><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mo>=</mml:mo><mml:mi>K</mml:mi><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi>K</mml:mi><mml:mrow><mml:mover accent="true"><mml:mi>Y</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E3"><mml:mtd/><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">σ</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mtext>a</mml:mtext></mml:msub><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mover accent="true"><mml:mi>Y</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo></mml:mrow><mml:mi>K</mml:mi></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:mi>K</mml:mi><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mn mathvariant="normal">1</mml:mn><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mtext>a</mml:mtext></mml:msub><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            with <inline-formula><mml:math id="M13" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">ρ</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mtext>a</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M14" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">σ</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mtext>a</mml:mtext></mml:msub></mml:mrow></mml:math></inline-formula> being the apparent
complex resistivity and apparent complex conductivity, respectively. These
quantities are referred to as “apparent” because they actually only represent the
true properties if the medium under investigation is homogeneous. Otherwise,
they represent an effective (average) value. Spatial discrimination of
electrical properties can be achieved by the use of multiple measurements
with different electrode locations, which also form the basis for tomographic
processing (inversion), i.e., imaging.</p>
      <p>Impedance measurements can be conducted using only two electrodes for a
combined current- and voltage
measurement, or by using four electrodes (quadrupole measurements, also
called four-point spreads) with separate current and voltage electrode pairs.
In the latter case, the contact impedance between electrode and medium, which
becomes significant towards lower measurement frequencies, has practically no
influence on the voltage measurement <xref ref-type="bibr" rid="bib1.bibx10" id="paren.44"><named-content content-type="pre">e.g.,</named-content></xref>.</p>
</sec>
<sec id="Ch1.S2.SS3">
  <title>Polarization of biomatter</title>
      <p>Polarization phenomena of biomatter are commonly classified into three
frequency regions with different polarization sources, namely the <inline-formula><mml:math id="M15" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula>,
<inline-formula><mml:math id="M16" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula>, and <inline-formula><mml:math id="M17" display="inline"><mml:mi mathvariant="italic">γ</mml:mi></mml:math></inline-formula> regions
<xref ref-type="bibr" rid="bib1.bibx83 bib1.bibx74" id="paren.45"><named-content content-type="pre">e.g.,</named-content></xref>. While overlapping, the
low-frequency <inline-formula><mml:math id="M18" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> polarization is thought to extend into the lower kHz
range, followed by the <inline-formula><mml:math id="M19" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula> polarization in the range up to about <inline-formula><mml:math id="M20" display="inline"><mml:mn mathvariant="normal">100</mml:mn></mml:math></inline-formula> MHz, and joined by the <inline-formula><mml:math id="M21" display="inline"><mml:mi mathvariant="italic">γ</mml:mi></mml:math></inline-formula> polarization at
higher frequencies <xref ref-type="bibr" rid="bib1.bibx76" id="paren.46"><named-content content-type="pre">e.g.,</named-content></xref>. Controlled by the
mobility of the charge carriers, the <inline-formula><mml:math id="M22" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> range is assumed to be
associated with electrochemical polarization (i.e., the build-up and
relaxation of ionic concentration gradients such as those found in EDLs, in
an electric, time-variable field); the <inline-formula><mml:math id="M23" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula> range by the Maxwell–Wagner
polarization of composite media
<xref ref-type="bibr" rid="bib1.bibx73 bib1.bibx74" id="paren.47"><named-content content-type="pre">e.g.,</named-content></xref>; and the
<inline-formula><mml:math id="M24" display="inline"><mml:mi mathvariant="italic">γ</mml:mi></mml:math></inline-formula> range by molecular, ionic, and atomic polarization. The different
processes lead to different current flow paths within biomatter for different
frequencies <xref ref-type="bibr" rid="bib1.bibx76" id="paren.48"/>. These observations have been primarily
made on cell suspensions and various kinds of tissue,
which exhibit structures
that are much more homogeneous than fully developed plant and root systems.
Polarization processes in plant roots are assumed to originate among others
in the cell membranes, the apoplast and the symplast <xref ref-type="bibr" rid="bib1.bibx77" id="paren.49"/>.
The frequency dependence of published multi-frequency measurements
<xref ref-type="bibr" rid="bib1.bibx70 bib1.bibx77" id="paren.50"><named-content content-type="pre">e.g.,</named-content></xref> indicates
multiple length scales and associated structures as the origin of electrical
polarization responses.</p>
      <p>On a cellular, or multi-cellular level, much work has been conducted to gain
information about the electrical surface characteristics of cells (cell
structures). A Gouy–Chapman-Stern model relating surface charges to external ion concentrations has been
formulated and subsequently improved
<xref ref-type="bibr" rid="bib1.bibx53 bib1.bibx51 bib1.bibx94" id="paren.51"/>. Using this
model, ion activity at membrane surfaces can be computed and analyzed for the
investigation of physiological effects. These and subsequent studies
regarding ion toxicity and related surface electric potential have provided
further evidence that certain surface potentials can be linked to
physiological states and processes, e.g., ion availability and uptake
<xref ref-type="bibr" rid="bib1.bibx93 bib1.bibx94 bib1.bibx95 bib1.bibx52" id="paren.52"/>.
<xref ref-type="bibr" rid="bib1.bibx61" id="text.53"/> estimated the electric potential at rice-root surfaces
of macroscopic root segments using measurements of the electrokinetic
zeta-potential. The zeta-potential is the experimentally accessible electric
potential at a distance from the surface where slipping in the electrolyte
occurs upon a flow-driving pressure gradient.</p>
      <p>The EDL is the source of polarization responses in the low-frequency range
usually measured with EIS/EIT. It is sensitive to physiological processes
that affect ion (nutrient) availability in the vicinity of, and ion fluxes
across, charged cell membranes. The key function of roots is the uptake of
water and nutrients, which is highly dependent on nutrient availability,
demand, and stress factors
<xref ref-type="bibr" rid="bib1.bibx23 bib1.bibx28 bib1.bibx43" id="paren.54"><named-content content-type="pre">e.g.,</named-content></xref>.
Nutrient availability can influence water (and nutrient) transport within
plant systems
<xref ref-type="bibr" rid="bib1.bibx24 bib1.bibx67" id="paren.55"><named-content content-type="pre">e.g.,</named-content></xref>, and
nutrient availability within roots can fluctuate in response to certain
depletion situations <xref ref-type="bibr" rid="bib1.bibx13" id="paren.56"><named-content content-type="pre">e.g.,</named-content></xref>.
Furthermore, the distribution of stress hormones such as
abscisic acid (ABA) increases in response to
stress situations, possibly inducing the aforementioned reactions
<xref ref-type="bibr" rid="bib1.bibx82" id="paren.57"><named-content content-type="pre">e.g.,</named-content></xref>. The formation and properties of
large-scale ion-selective structures such as endodermis and hypodermis are
also directly influenced by the growth environment, and can change in
response to external stress factors <xref ref-type="bibr" rid="bib1.bibx43" id="paren.58"/>. In addition,
<xref ref-type="bibr" rid="bib1.bibx27" id="text.59"/> noted that electrical polarization effects originate
in the “active” parts of a root system only, which change according to age,
nutrient availability, and other stress factors <xref ref-type="bibr" rid="bib1.bibx6" id="paren.60"><named-content content-type="pre">see
also</named-content></xref>.</p>
      <p>The majority of studies concerning full root systems work with equivalent
electrical circuit models to describe the measured signals of various
biostructures  <xref ref-type="bibr" rid="bib1.bibx76" id="paren.61"><named-content content-type="pre">see</named-content><named-content content-type="post">and references therein</named-content></xref>. The scale
and composition of these models vary considerably. For example,
<xref ref-type="bibr" rid="bib1.bibx27" id="text.62"/> equates root segments with cylindrical capacitors,
whose conducting plates are formed by the inner xylem and the fluid surrounding
the root segment, with root cortex acting as a dielectric.
<xref ref-type="bibr" rid="bib1.bibx56" id="text.63"/> proposed a simplified model of cell polarization in
root systems, in which the cell membrane acts as a dielectric between the
conducting inner and outer regions of the cells, thus representing a classical
capacitor. Equivalent circuit representations inherently depend on the assumed
flow paths of the electric current. For instance, impedance measurements using
the stem as one pole for current injection and the medium surrounding the roots
as the other pole  <xref ref-type="bibr" rid="bib1.bibx22 bib1.bibx29 bib1.bibx77" id="paren.64"><named-content content-type="pre">as frequently being done,
e.g.,</named-content></xref> force the
current to cross all radial layers of the roots.  However, even for stem
injection, equivalent circuit models considerably simplify the true electrical
processes in the root and root–rhizosphere system, and it is questionable
whether these models can be transferred between different experimental setups
<xref ref-type="bibr" rid="bib1.bibx27 bib1.bibx70 bib1.bibx29" id="paren.65"><named-content content-type="pre">as evident from the large number of slightly different models that were
proposed,
e.g.,</named-content></xref>.  A
purely phenomenological analysis is made by <xref ref-type="bibr" rid="bib1.bibx77" id="text.66"/>, who use a
classification approach to analyze spectral impedance data measured on pine
roots infested with mycorrhizal fungi.</p>
</sec>
<sec id="Ch1.S2.SS4">
  <title>Working hypotheses</title>
      <p>We propose to describe and interpret low-frequency (<inline-formula><mml:math id="M25" display="inline"><mml:mo>&lt;</mml:mo></mml:math></inline-formula> 1 kHz) polarization
processes in biomatter using concepts similar to those established for soils
and rocks in recent years, under the assumption that the observed responses
originate from the polarization of EDLs present in the biomatter.
Accordingly, it should be possible to link the polarization magnitude to the
average EDL thickness, which depends on the electric potential drop between
the charged surface/membrane and the background ambient electrolyte, and link
characteristic relaxation times to the length scales at which the
polarization processes take place.</p>
      <p>Given the current observations and understanding of electrical polarization
processes in biomatter, as reviewed in the previous section, our hypotheses
are as follows:</p>
      <p><list list-type="order">
            <list-item>

      <p>The magnitude of the low-frequency polarization response of roots is
related to the overall surface area comprised by EDLs in the
root–rhizosphere system, including the inner root structure. EDLs may
form at Casparian strips (e.g., hypodermis and endodermis)
and plasma membranes.</p>
            </list-item>
            <list-item>

      <p>The characteristic relaxation times of the low-frequency polarization
response of roots provide information on the length scales at which the
polarization processes take place. While it is not clear to what
extent a discrimination of specific polarization processes (e.g., plasma
membrane polarization and polarization of the hypodermis) is possible,
changes in the relaxation times should indicate changes in the length
scale of the polarizing structures.</p>
            </list-item>
            <list-item>

      <p>EDLs in the inner root system are influenced by ions (nutrients) in
the sap fluid, EDLs at the outer root surface are influenced by ion
concentrations in the external fluid. Thus, physiological processes
that influence the availability, usage, and translocation of ions
directly influence the low-frequency polarization response.</p>
            </list-item>
            <list-item>

      <p>Spectral EIT is a suitable non-invasive method to image and monitor
magnitude and characteristic relaxation times of the low-frequency
polarization response of root systems.</p>
            </list-item>
          </list></p>
      <p>In the present study, we address the second part of hypothesis three, as well as
hypothesis four. Hypotheses one and two are based on the synthesis of existing
works, but can neither be validated nor invalidated by the present study.</p>
</sec>
</sec>
<sec id="Ch1.S3">
  <title>Material and methods</title>
<sec id="Ch1.S3.SS1">
  <title>Electrical impedance tomography</title>
      <p>The EIS (or SIP) method involves the measurement of impedances at multiple
frequencies (usually in the mHz to kHz range). It can be extended by
utilizing electrode arrays consisting of tens to hundreds of electrodes to
collect numerous, spatially distributed four-point impedance measurements.
From these data sets images of the complex conductivity (or its inverse,
complex resistivity) can be computed using tomographic inversion algorithms
<xref ref-type="bibr" rid="bib1.bibx48 bib1.bibx26" id="paren.67"><named-content content-type="pre">e.g.,</named-content></xref>. This method is called
complex conductivity (or complex resistivity) imaging or electrical
impedance tomography (EIT), and refers to both single- and multi-frequency
(spectral) approaches. EIT images are characterized by a spatially variable
resolution, which decreases with increasing distance from the electrodes
<xref ref-type="bibr" rid="bib1.bibx3 bib1.bibx38 bib1.bibx14 bib1.bibx26" id="paren.68"><named-content content-type="pre">e.g.,</named-content></xref>.
The method's primary fields of application are in near-surface geophysics
<xref ref-type="bibr" rid="bib1.bibx14 bib1.bibx26 bib1.bibx79" id="paren.69"><named-content content-type="pre">e.g.,</named-content></xref> and
medical imaging <xref ref-type="bibr" rid="bib1.bibx11" id="paren.70"><named-content content-type="pre">e.g.,</named-content></xref>.</p>
      <p>Spectral EIT measurements presented in this study were conducted using the
40-channel EIT-40 impedance tomograph <xref ref-type="bibr" rid="bib1.bibx102" id="paren.71"/>, which
was configured in a monitoring setup to acquire up to seven EIT data sets
(frames) from a mini rhizotron per day.</p>
<sec id="Ch1.S3.SS1.SSS1">
  <title>2-D forward modeling</title>
      <p>Synthetic impedance data, required in the tomographic inversion process, were
modeled using the finite-element (FE) forward modeling code of
<xref ref-type="bibr" rid="bib1.bibx48" id="text.72"/>. The code solves the Poisson equation for a 2-D complex
conductivity distribution and 2-D source currents in a domain of given thickness
<xref ref-type="bibr" rid="bib1.bibx36" id="paren.73"/>. At the boundaries of the 2-D modeling
domain, no-flow Neumann conditions
are imposed, which do not allow any current flow out of the modeling
domain. Details of the
implementation can be found in <xref ref-type="bibr" rid="bib1.bibx48" id="text.74"/>.</p>
      <p>A sketch of the FE grid (also used for the inversion and presentation of
imaging results) resembling the rhizotron is shown in Fig. <xref ref-type="fig" rid="Ch1.F1"/>b along with the position of
38 electrodes. The grid consists of 60 elements in the <inline-formula><mml:math id="M26" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> direction, and 157
elements in the <inline-formula><mml:math id="M27" display="inline"><mml:mi>z</mml:mi></mml:math></inline-formula> direction (9420 elements in total).</p>
</sec>
<sec id="Ch1.S3.SS1.SSS2">
  <title>2-D tomographic inversion</title>
      <p>Complex conductivity images at multiple measurement frequencies were computed
using the smoothness-constraint inversion code of <xref ref-type="bibr" rid="bib1.bibx48" id="text.75"/>. The
code computes the distribution of complex conductivity <inline-formula><mml:math id="M28" display="inline"><mml:mover accent="true"><mml:mi mathvariant="italic">σ</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover></mml:math></inline-formula>
(expressed in either magnitude (<inline-formula><mml:math id="M29" display="inline"><mml:mrow><mml:mo>|</mml:mo><mml:mover accent="true"><mml:mi mathvariant="italic">σ</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mo>|</mml:mo></mml:mrow></mml:math></inline-formula>) and phase (<inline-formula><mml:math id="M30" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula>) or real
component (<inline-formula><mml:math id="M31" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula>) and imaginary component
(<inline-formula><mml:math id="M32" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>))-in the 2-D (<inline-formula><mml:math id="M33" display="inline"><mml:mrow><mml:mi>x</mml:mi><mml:mo>,</mml:mo><mml:mi>z</mml:mi></mml:mrow></mml:math></inline-formula>)
image plane from the given set of complex transfer impedances
<inline-formula><mml:math id="M34" display="inline"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>)–(expressed in magnitude
(<inline-formula><mml:math id="M35" display="inline"><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub><mml:mo>|</mml:mo></mml:mrow></mml:math></inline-formula>) and phase (<inline-formula><mml:math id="M36" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">φ</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>)) under the constraint of maximum
model smoothness. Log-transformed impedances and log-transformed complex
conductivities (of the individual elements of the grid) are used as data and
model parameters, respectively, in the inversion.</p>
      <p>The iterative, Gauss–Newton type of inversion scheme minimizes an objective
function composed of measures of data misfit and model roughness. The data
misfit is weighted by the (real-valued) magnitudes of individual,
complex-valued data errors <xref ref-type="bibr" rid="bib1.bibx48" id="paren.76"><named-content content-type="pre">see</named-content></xref>, which, however, are
dominated by the resistance errors since the phase values are relatively
small for the measurements considered here. Therefore the resistance error
model <xref ref-type="bibr" rid="bib1.bibx58" id="paren.77"/>

                  <disp-formula id="Ch1.E4" content-type="numbered"><mml:math id="M37" display="block"><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mi mathvariant="normal">Δ</mml:mi><mml:mo>|</mml:mo><mml:msub><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub><mml:mo>|</mml:mo><mml:mo>=</mml:mo><mml:mi>a</mml:mi><mml:mo>|</mml:mo><mml:msub><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub><mml:mo>|</mml:mo><mml:mo>+</mml:mo><mml:mi>b</mml:mi></mml:mrow></mml:math></disp-formula>

            can be used in the complex inversion for the weighting of complex data
(including the phase), with <inline-formula><mml:math id="M38" display="inline"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi><mml:mo>|</mml:mo><mml:msub><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub><mml:mo>|</mml:mo></mml:mrow></mml:math></inline-formula> being the error of impedance
magnitude (resistance) <inline-formula><mml:math id="M39" display="inline"><mml:mrow><mml:mo>|</mml:mo><mml:msub><mml:mover accent="true"><mml:mi>Z</mml:mi><mml:mo mathvariant="normal" stretchy="false">^</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub><mml:mo>|</mml:mo></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M40" display="inline"><mml:mi>a</mml:mi></mml:math></inline-formula> a relative error contribution, and <inline-formula><mml:math id="M41" display="inline"><mml:mi>b</mml:mi></mml:math></inline-formula>
an absolute error contribution. For more details on the inversion scheme we
refer the reader to <xref ref-type="bibr" rid="bib1.bibx48" id="text.78"/>. The inversion is separately performed for each frequency
of the given data set.</p>
</sec>
</sec>
<sec id="Ch1.S3.SS2">
  <title>Debye decomposition</title>
      <p>The Debye decomposition (DD) approach
<xref ref-type="bibr" rid="bib1.bibx90 bib1.bibx60" id="paren.79"><named-content content-type="pre">e.g.,</named-content></xref> was used to
analyze the complex conductivity spectra recovered from the multi-frequency
EIT inversion results. The approach yields integral parameters describing the
spectral characteristics of the SIP signature. The complex conductivity
spectrum is represented as a superposition of a large number of Debye
relaxation terms at relaxation times <inline-formula><mml:math id="M42" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> (suitably distributed over the
range implicitly defined by the data frequency limits; see
<xref ref-type="bibr" rid="bib1.bibx96" id="altparen.80"/>):
            <disp-formula id="Ch1.E5" content-type="numbered"><mml:math id="M43" display="block"><mml:mrow><mml:mover accent="true"><mml:mi mathvariant="italic">σ</mml:mi><mml:mo stretchy="false" mathvariant="normal">^</mml:mo></mml:mover><mml:mo>(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>)</mml:mo><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="italic">σ</mml:mi><mml:mi mathvariant="normal">∞</mml:mi></mml:msub><mml:mfenced close="]" open="["><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mi>k</mml:mi></mml:munder><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>+</mml:mo><mml:mi>j</mml:mi><mml:mi mathvariant="italic">ω</mml:mi><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
          with <inline-formula><mml:math id="M44" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">σ</mml:mi><mml:mi mathvariant="normal">∞</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> being the (real-valued) conductivity in the
high-frequency limit, and <inline-formula><mml:math id="M45" display="inline"><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> the <inline-formula><mml:math id="M46" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula>th chargeability, describing the
relative weight of the <inline-formula><mml:math id="M47" display="inline"><mml:mi>k</mml:mi></mml:math></inline-formula>th Debye relaxation term in the decomposition. The
chargeabilities <inline-formula><mml:math id="M48" display="inline"><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> at the different relaxation times <inline-formula><mml:math id="M49" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> form a
relaxation time distribution (RTD), from which the following descriptive
parameters are computed <xref ref-type="bibr" rid="bib1.bibx69" id="paren.81"><named-content content-type="pre">e.g.,</named-content></xref>:</p>
      <p><list list-type="bullet">
            <list-item>

      <p>The normalized total chargeability <inline-formula><mml:math id="M50" display="inline"><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula> is
a measure of the overall polarization reflected in the spectrum
<xref ref-type="bibr" rid="bib1.bibx87 bib1.bibx96" id="paren.82"><named-content content-type="pre">e.g.,</named-content></xref>:

                      <disp-formula id="Ch1.E6" content-type="numbered"><mml:math id="M51" display="block"><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:msub><mml:mi mathvariant="italic">σ</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mi>k</mml:mi></mml:munder><mml:msub><mml:mi>m</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>

                with <inline-formula><mml:math id="M52" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">σ</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> being the (real-valued) conductivity in the low-frequency
limit.</p>
            </list-item>
            <list-item>

      <p>The mean logarithmic relaxation time <inline-formula><mml:math id="M53" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi mathvariant="normal">mean</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>
represents a weighted mean of the RTD:
                  <disp-formula id="Ch1.E7" content-type="numbered"><mml:math id="M54" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi mathvariant="normal">mean</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi>exp⁡</mml:mi><mml:mfenced open="(" close=")"><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mo>∑</mml:mo><mml:mi>k</mml:mi></mml:msub><mml:msub><mml:mi>m</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mi>log⁡</mml:mi><mml:mo>(</mml:mo><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mo>)</mml:mo></mml:mrow><mml:mrow><mml:msub><mml:mo>∑</mml:mo><mml:mi>k</mml:mi></mml:msub><mml:msub><mml:mi>m</mml:mi><mml:mi>k</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle></mml:mfenced><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula></p>
            </list-item>
            <list-item>

      <p>The uniformity parameter <inline-formula><mml:math id="M55" display="inline"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mn mathvariant="normal">60</mml:mn><mml:mo>,</mml:mo><mml:mn mathvariant="normal">10</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> describes the
frequency dispersion of the spectrum:
                  <disp-formula id="Ch1.E8" content-type="numbered"><mml:math id="M56" display="block"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mn mathvariant="normal">60</mml:mn><mml:mo>,</mml:mo><mml:mn mathvariant="normal">10</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mn mathvariant="normal">60</mml:mn></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mn mathvariant="normal">10</mml:mn></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
                with <inline-formula><mml:math id="M57" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mn mathvariant="normal">10</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M58" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mn mathvariant="normal">60</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> being the relaxation times at which the
cumulative chargeability reaches 10 and 60 %, respectively, of the
total chargeability sum.</p>
            </list-item>
          </list></p>
      <p>The implementation of <xref ref-type="bibr" rid="bib1.bibx96" id="text.83"/> was used for the DD analysis.
The iterative inversion scheme balances between (error-weighted) data fitting
and smoothing requirements.</p>
</sec>
<sec id="Ch1.S3.SS3">
  <title>Experimental setup</title>
<sec id="Ch1.S3.SS3.SSS1">
  <title>Rhizotron</title>

      <?xmltex \floatpos{t}?><fig id="Ch1.F1"><caption><p><bold>(a)</bold> Experimental setup of plant root systems in water
filled rhizotron. <bold>(b)</bold> Corresponding finite-element grid used for
electrical modeling and imaging. Red dots indicate position of electrodes.
Red arrows indicate the ascending order of electrode numbering, some of which
are marked using the notation E1–E38.</p></caption>
            <?xmltex \igopts{width=170.716535pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f01.png"/>

          </fig>

      <p>The experiment was conducted using a mini-rhizotron with the dimensions of
30 cm width, 78 cm height, and 2 cm thickness, and a transparent front
plate (Fig. <xref ref-type="fig" rid="Ch1.F1"/>). The front of the rhizotron is equipped
with 38 brass pins of 5 mm diameter as electrodes, which do not extend into
the rhizotron's inner volume. A growth lamp was installed above the rhizotron
and turned on during daylight hours.</p>
</sec>
<sec id="Ch1.S3.SS3.SSS2">
  <title>Plant treatment</title>
      <p>Oilseed plants had been grown in nutrient solution prior to the experiment.
To increase the root mass, three plants were tied together and centrally
placed at the top of the rhizotron (Fig. <xref ref-type="fig" rid="Ch1.F1"/>), which had been filled with tap water before. No water was added to the rhizotron during
the experiment, nor was the water in the rhizotron disturbed in any way.
Thus, it is possible that at some point anaerobic conditions manifested.</p>
</sec>
<sec id="Ch1.S3.SS3.SSS3">
  <title>Data acquisition</title>
      <p>Over the course of 3 ḋays, 21 EIT data sets at 35 frequencies between
0.46 Hz and 45 kHz were collected in regular intervals, starting right
after the placement of the plant system in the rhizotron. A total of 1158
quadrupoles were measured for each data set, involving 74 individual current
injections (plus 767 reciprocal configurations, where current and voltage
electrode pairs are interchanged, for quality assessment), requiring less
than 4 h acquisition time. These quadrupoles consisted mostly of skip-0 and
skip-2 (numbers of electrodes between the two electrodes used for current
injection and voltage measurement, respectively) dipole–dipole
configurations, as well as quadrupoles with current electrodes on opposite
sides (left and right) of the rhizotron and skip-0 voltage readings.</p>
</sec>
</sec>
<sec id="Ch1.S3.SS4">
  <title>Data processing</title>
<sec id="Ch1.S3.SS4.SSS1">
  <title>Selection of impedance data</title>
      <p>The inversion scheme assumes normally distributed and uncorrelated data
errors and is very sensitive to outliers
<xref ref-type="bibr" rid="bib1.bibx57" id="paren.84"><named-content content-type="pre">e.g.,</named-content></xref>. Outliers are usually associated with
low signal-to-noise ratios or systematic errors due to missing or bad
electrode contacts. Outliers can either be removed from the data set prior to
inversion or accounted for by sophisticated, “robust” inversion schemes
<xref ref-type="bibr" rid="bib1.bibx57" id="paren.85"/>. In these robust schemes, the weighting of
individual data points is iteratively adapted, which can lead to a reduction
of spatial resolution as well as recovered contrast in the imaging results.
However, usually this does not change the qualitative results of the
inversion. In the present study, we sought to analyze data across the
frequency and time domains, which requires a careful and consistent analysis
of the inversion data. Thus, to prevent introducing unnecessary variations
between time steps and frequencies, we opted to remove outliers using the
criteria described below and use individual, but consistent, data weighting
schemes for all measurements.</p>
      <p>The measured impedance data (also referred to as “raw data”, in contrast to
complex conductivity data recovered from the imaging results, referred to as
“intrinsic data”) were screened (filtered) for outliers and faulty data
according to multiple criteria: First, outliers were identified for each
frequency and time step and removed from the data set. Due to the underlying
physical principles, EIT measurements usually do not show strong variations
when electrode positions are only slightly shifted. The exception here is
measurements with electrodes located close to the plant stem system, where a
very localized anomalous response is expected in the data. Accordingly, care
was taken not to remove these data as outliers.
Second, following this selection process,
only impedance spectra were kept that retained more than 90 % of the
original data points below 300 Hz and showed consistency over several time
steps. Third, to avoid errors due to
electromagnetic coupling effects
<xref ref-type="bibr" rid="bib1.bibx71 bib1.bibx101" id="paren.86"><named-content content-type="pre">e.g.,</named-content></xref>, only data below 220 Hz
were considered for the imaging. Fourth, measurements at 50 Hz were discarded due to powerline noise.</p>
      <p>The applied data selection criteria resulted in small variations in the
number of measurements actually used for the inversions for the different
time steps, ranging between 530 and 555 measurements per data set and
frequency. The average injected current strength of the measurements at each
time step increased slightly over time from approximately 1.0 to 1.2 mA.</p>
</sec>
<sec id="Ch1.S3.SS4.SSS2">
  <?xmltex \opttitle{Correction of impedance data for imperfect\hack{\break} 2-D situations}?><title>Correction of impedance data for imperfect<?xmltex \hack{\break}?> 2-D situations</title>
      <p>Since the electrodes do not extend across the entire depth (i.e., horizontal
direction perpendicular to the image plane) of the rhizotron, the electric
current and potential field distributions in the rhizotron are not perfectly
2-D, as is assumed in the forward modeling. Therefore measurements were
conducted on a rhizotron solely filled with tap water of known conductivity.
By comparing the latter with the apparent conductivity
(Eq. <xref ref-type="disp-formula" rid="Ch1.E3"/>) derived from the measured impedance and the
numerically determined geometric factor (obtained from running the forward
model for a homogeneous case) for each measurement configuration, correction
factors were computed and applied to all measured impedances.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F2"><caption><p>EIT inversion results (real component of complex conductivity) for
measurements on a rhizotron filled with water of known conductivity
(375 <inline-formula><mml:math id="M59" display="inline"><mml:mi mathvariant="normal">µ</mml:mi></mml:math></inline-formula>S cm<inline-formula><mml:math id="M60" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> <inline-formula><mml:math id="M61" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> <inline-formula><mml:math id="M62" display="inline"><mml:mrow><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1.43</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> S cm<inline-formula><mml:math id="M63" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>)
without <bold>(a)</bold> and with <bold>(b)</bold> correction of the impedance data
for the imperfect 2-D situation.</p></caption>
            <?xmltex \igopts{width=156.490157pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f02.pdf"/>

          </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F3"><caption><p>EIT inversion results (real component of complex conductivity)
obtained from synthetic data using a modeling grid in the inversion with
a <bold>(a)</bold> 1.5 cm lower, <bold>(b)</bold> identical, and <bold>(c)</bold> 1.5 cm higher position of the top
boundary compared to the grid (forward model) used to compute synthetic measurments.
The forward model was homogeneously parameterized with a conductivity
distribution of 0.1 S m<inline-formula><mml:math id="M64" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>. Only the upper part of the modeling domain
(rhizotron) is shown. Electrode positions and measurement configurations are
the same for all three cases.</p></caption>
            <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f03.png"/>

          </fig>

      <p>From a theoretical point of view, these correction factors are independent of
the conductivity value of a homogeneous distribution since measured
resistance and resistivity (inverse of conductivity) are linearly related for
a homogeneous distribution. Significant changes in the correction factors
can only occur for strong spatial conductivity variations, in particular
across the thickness of the rhizotron (2-D/3-D effects). However, even if
present, such effects in the correction factors would primarily result in
inaccuracies in the inverted conductivity magnitude image, while the
conductivity-phase image and also the DD-derived spectral parameters (total
chargeability, relaxation time) are relatively robust against magnitude,
i.e., correction factor, errors. We therefore, for the small to moderate
conductivity variations observed in the experiment in the upper region of the
rhizotron, assume that the conducted calibration survey, i.e., one universal
set of correction factors, was actually sufficient.</p>
      <p>In Fig. <xref ref-type="fig" rid="Ch1.F2"/> the effect of this correction procedure
on the EIT inversion result is shown. Without correction the obtained image
exhibits an artificial pattern (Fig. <xref ref-type="fig" rid="Ch1.F2"/>a), while
with correction a practically homogeneous distribution is recovered, in
agreement with the conductivity of the tap water
(Fig. <xref ref-type="fig" rid="Ch1.F2"/>b).</p>
      <p>The inversion was conducted using the error parameter values <inline-formula><mml:math id="M65" display="inline"><mml:mrow><mml:mi>a</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.5</mml:mn></mml:mrow></mml:math></inline-formula> %
and <inline-formula><mml:math id="M66" display="inline"><mml:mrow><mml:mi>b</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.012</mml:mn><mml:mi mathvariant="normal">Ω</mml:mi></mml:mrow></mml:math></inline-formula> (see Eq. <xref ref-type="disp-formula" rid="Ch1.E4"/>). These values were found
to be appropriate and were also used in the inversions, of which the results
are shown in the following.</p>
</sec>
<sec id="Ch1.S3.SS4.SSS3">
  <title>Adaptation of modeling domain to changing water table</title>

      <?xmltex \floatpos{t}?><fig id="Ch1.F4"><caption><p>Photographs of the oilseed plants during the experiment:
<bold>(a)</bold> close-up of the root systems in the rhizotron container,
<bold>(b)</bold> day 1 and <bold>(c)</bold> day 3. The colored dots in
<bold>(a)</bold> indicate the approximate position at which intrinsic
signatures, recovered from tomographic inversion results, are investigated;
red dot: stem area; blue dot: fine root area.</p></caption>
            <?xmltex \igopts{width=236.157874pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f04.jpg"/>

          </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F5" specific-use="star"><caption><p>Temporal evolution of raw data spectra, plotted as real <bold>(a, b)</bold> and imaginary <bold>(c, d)</bold> components of apparent complex conductivity
(Eq. <xref ref-type="disp-formula" rid="Ch1.E3"/>) for two example measurement configurations:
<bold>(a, c)</bold> current injection between electrodes 3 and 4, and voltage
measurement between electrodes 5 and 6 (quadrupole located directly above the
root system, i.e., response “with roots”); <bold>(b, d)</bold> current
injection between electrodes 34 and 35 and voltage measurement between
electrodes 36 and 37 (quadrupole located in an area relatively far away from
the root system, i.e., “water-only” response). For electrode numbering, see
Fig. <xref ref-type="fig" rid="Ch1.F1"/>b. Blue indicates early measurements, while
later ones are shown in red. Values of <inline-formula><mml:math id="M67" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> that lie below the
measurement accuracy of the system (1 mrad phase shift at 1 kHz for water;
see <xref ref-type="bibr" rid="bib1.bibx102" id="altparen.87"/>) are indicated by gray areas. </p></caption>
            <?xmltex \igopts{width=298.753937pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f05.png"/>

          </fig>

      <p>Due to evaporation and root water uptake the water table fell by ca. 2 cm
over the course of the monitoring experiment. This was not problematic in
terms of electrode contact as electrodes always remained in the water.
However, the changing water table has to be accounted for in the EIT
inversion by means of an adapted modeling domain, i.e., by adapting the
position of the top boundary of the FE modeling grid, where no-flow
conditions are assumed. Otherwise, as we checked in numerical experiments,
significant artifacts appear in the inversion results
(Fig. <xref ref-type="fig" rid="Ch1.F3"/>).</p>
      <p>From the known water tables at the beginning and end of the experiment, and the
average time when each EIT data set was collected, the positions of the top
boundary of the individual grids used for the inversion of each data set were
determined by linear interpolation.</p><?xmltex \hack{\newpage}?>
</sec>
<sec id="Ch1.S3.SS4.SSS4">
  <title>Analysis of spectral imaging results</title>
      <p>The spectral imaging results were analyzed by means of pixel-wise application
of the Debye decomposition scheme. As water exhibits no significant
polarization response in the examined frequency range, the area free of roots
from 20 to 78 cm depth of the rhizotron was used to quantify a
<inline-formula><mml:math id="M68" display="inline"><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula> threshold value below which polarization is
considered insignificant. Based on this threshold value all images of the spectral
parameters obtained from the Debye decomposition, including the top 20 cm of
the rhizotron, were partitioned into pixels with and without significant
polarization. The observed polarization can be fully attributed to the root
system (no polarization is expected from the surrounding water in the
examined frequency range) and thus the corresponding pixels delineate
polarizable areas of the root system, which we refer to as the root pixel
group.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F6" specific-use="star"><caption><p>Single-frequency inversion results in terms of real <bold>(b, d)</bold>
and imaginary <bold>(c, e)</bold> components of complex conductivity for
1 Hz <bold>(b, c)</bold> and 70 Hz <bold>(c, e)</bold> at the first time step. The
photograph of the root system at this time <bold>(a)</bold> shows the same area
of the rhizotron as the inversion results.</p></caption>
            <?xmltex \igopts{width=298.753937pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f06.png"/>

          </fig>

      <p>To analyze the temporal evolution of the overall root system polarization (in
terms of normalized total chargeability) the root pixel group was determined
for the first time step, and then kept fixed for the following time steps.
Relaxation times, however, can only be reliably extracted from SIP signatures
if they show significant polarization. Therefore, for the relaxation time
analysis (in terms of mean relaxation time and uniformity parameter) the root
pixel groups were determined for each time step individually.</p><?xmltex \hack{\newpage}?>
</sec>
</sec>
</sec>
<sec id="Ch1.S4">
  <title>Results</title>
<sec id="Ch1.S4.SS1">
  <title>Physiological response</title>
      <p>Photographs of the plant systems at the beginning and the end of the
experiment are shown in Fig. <xref ref-type="fig" rid="Ch1.F4"/>. As is evident from
the photographs, the plants significantly reacted to the nutrient stress
situation (and possibly anaerobic conditions) and degraded over time. The
root systems extended down to a depth of approximately 13 cm, with an
approximate maximum lateral extension of 13 cm (see Fig. <xref ref-type="fig" rid="Ch1.F1"/>).</p>
</sec>
<sec id="Ch1.S4.SS2">
  <title>Impedance spectra</title>
      <p>In Fig. <xref ref-type="fig" rid="Ch1.F5"/> the temporal evolution of the raw data spectra
in terms of apparent complex conductivity (Eq. <xref ref-type="disp-formula" rid="Ch1.E3"/>) is
shown for two example measurement configurations: a quadrupole with electrode
pairs on both sides of the rhizotron located directly above the root system,
i.e., with sensitivity to the root system and a quadrupole with electrodes
from the horizontal electrode line at 37 cm depth, i.e., located relatively
far away from the root system and thus sensitive only to the water. The real
component of apparent complex conductivity shows a smooth, consistent
behavior across the time and frequency domains for both responses, “with
roots” and “water-only” (Fig. <xref ref-type="fig" rid="Ch1.F5"/>a, b). However, the
conductivity decreases for the quadrupole around the root system, while it
increases in the “water-only” quadrupole. The imaginary components, i.e.,
the polarization responses, with roots (Fig. <xref ref-type="fig" rid="Ch1.F5"/>c) are
also consistent and show changes, especially in the lower-frequency range,
over time. The water-only measurements, on the other hand, only exhibit
negligible polarization responses, likely dominated by measurement errors and
noise (Fig. <xref ref-type="fig" rid="Ch1.F5"/>d). The polarization magnitudes, on the one
hand, lie well below the signal threshold that can be reliably measured with
the EIT-40 system <xref ref-type="bibr" rid="bib1.bibx102" id="paren.88"/>; the jittery shape of these
signatures (Fig. <xref ref-type="fig" rid="Ch1.F5"/>d) is attributed to the logarithmic
scale of the plot. On the other hand, measured root signatures lie clearly
above the measurement threshold of the system (Fig. <xref ref-type="fig" rid="Ch1.F5"/>c).</p>
</sec>
<sec id="Ch1.S4.SS3">
  <title>Single-frequency imaging results</title>
      <p>The spatial variability of the
electrical response was assessed using the complex conductivity imaging
results, i.e., <inline-formula><mml:math id="M69" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M70" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>, at the first time step for the two
frequencies 1 and 70 Hz (Fig. <xref ref-type="fig" rid="Ch1.F6"/>). Only weak
variations in the real component (in-phase conductivity) can be observed at
the location of the root system (Fig. <xref ref-type="fig" rid="Ch1.F6"/>b, d).
However, a significant polarization response in the imaginary component
(Fig. <xref ref-type="fig" rid="Ch1.F6"/>c, e) coincides with the extension of the root
system. The frequency dependence previously found in the apparent complex
conductivity spectra (see Fig. <xref ref-type="fig" rid="Ch1.F5"/>) is also revealed in
the imaging results, with a stronger response at 70 than at 1 Hz. It
manifests both in signal strength and in the spatial extension of the
polarizable anomaly associated with the root system.</p>
</sec>
<sec id="Ch1.S4.SS4">
  <title>Complex conductivity spectra recovered from imaging results</title>

      <?xmltex \floatpos{t}?><fig id="Ch1.F7" specific-use="star"><caption><p>Intrinsic complex conductivity spectra (in terms of real component
<inline-formula><mml:math id="M71" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula> and imaginary component <inline-formula><mml:math id="M72" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) for all time steps (indicated
by color of the curves) recovered from the multi-frequency inversion results
at different locations in the rhizotron: <bold>(a, d)</bold> stem area (<inline-formula><mml:math id="M73" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula>
position: 20.25 cm; <inline-formula><mml:math id="M74" display="inline"><mml:mi>z</mml:mi></mml:math></inline-formula> position: <inline-formula><mml:math id="M75" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>0.75 cm); <bold>(b, e)</bold> bottom of
the root system (<inline-formula><mml:math id="M76" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> position: 20.25 cm; <inline-formula><mml:math id="M77" display="inline"><mml:mi>z</mml:mi></mml:math></inline-formula> position <inline-formula><mml:math id="M78" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>6.25 cm);
<bold>(c, f)</bold> water-only location (<inline-formula><mml:math id="M79" display="inline"><mml:mi>x</mml:mi></mml:math></inline-formula> position: 20.25 cm; <inline-formula><mml:math id="M80" display="inline"><mml:mi>z</mml:mi></mml:math></inline-formula> position:
<inline-formula><mml:math id="M81" display="inline"><mml:mo>-</mml:mo></mml:math></inline-formula>50.25 cm). Positions for <bold>(a and d)</bold> and <bold>(b and e)</bold> are also
indicated in Fig. <xref ref-type="fig" rid="Ch1.F4"/>a. Values of <inline-formula><mml:math id="M82" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> that lie
below the measurement accuracy of the system (1 mrad phase shift at 1 kHz
for water; see <xref ref-type="bibr" rid="bib1.bibx102" id="altparen.89"/>) are indicated by gray
areas.</p></caption>
          <?xmltex \igopts{width=298.753937pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f07.pdf"/>

        </fig>

      <p>Complex conductivity spectra were extracted from the multi-frequency imaging
results at three locations: near the stem area of the root system, from the
lower area of the root system, and from the lower half of the rhizotron,
where no root segments were present. The two locations near the root system
are indicated in Fig. <xref ref-type="fig" rid="Ch1.F4"/>a. Thus, these locations
represent areas with relatively thick root segments, thin root
segments, and no root segments (water only) at all.
Figure <xref ref-type="fig" rid="Ch1.F7"/> shows the temporal evolution of the spectral
response for the three locations. The real component of complex conductivity
(<inline-formula><mml:math id="M83" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mo>′</mml:mo></mml:msup></mml:mrow></mml:math></inline-formula>) increased over the course of the experiment for all three
locations (Fig. <xref ref-type="fig" rid="Ch1.F7"/>a–c). The imaginary component of
complex conductivity (<inline-formula><mml:math id="M84" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula>) reveals a frequency-dependent polarization
response at the root segment locations for all time steps
(Fig. <xref ref-type="fig" rid="Ch1.F7"/>d, e). The polarization magnitude decreases
over time, and changes in the shape of the spectra can be observed for later
time steps. These changes are most pronounced for the location near the stem
(Fig. <xref ref-type="fig" rid="Ch1.F7"/>a). The polarization signatures recovered at
the bottom of the rhizotron (water only) show an amount that is almost 2 orders of magnitude
smaller (Fig. <xref ref-type="fig" rid="Ch1.F7"/>f) compared to those in the
root system areas; they contain more noise and do not exhibit a clear frequency
trend.</p>
</sec>
<sec id="Ch1.S4.SS5">
  <title>Debye decomposition of recovered complex conductivity spectra</title>
      <p>The Debye decomposition scheme was applied to the complex conductivity
spectra recovered from multi-frequency EIT to quantify the overall
polarization (normalized total chargeability <inline-formula><mml:math id="M85" display="inline"><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula>)
and the characteristic relaxation time (mean relaxation time
<inline-formula><mml:math id="M86" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi mathvariant="normal">mean</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>), as well as the uniformity parameter <inline-formula><mml:math id="M87" display="inline"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mn mathvariant="normal">60</mml:mn><mml:mo>,</mml:mo><mml:mn mathvariant="normal">10</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>. By
means of this analysis, the intrinsic spectra can be assessed with respect to
the magnitude and shape of the polarization response for all pixels at each
time step. Figure <xref ref-type="fig" rid="Ch1.F8"/> shows a decomposition result of
a pixel signature from the stem area for the
first time step, corresponding to the spectrum plotted in
Fig. <xref ref-type="fig" rid="Ch1.F7"/>a, d at “0 h” (dark blue curve). The complex
conductivity spectrum was fitted by means of 96 Debye relaxation terms
(Fig. <xref ref-type="fig" rid="Ch1.F8"/>a), yielding a relaxation time distribution
(RTD) (Fig. <xref ref-type="fig" rid="Ch1.F8"/>b), from which <inline-formula><mml:math id="M88" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi mathvariant="normal">mean</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>,
<inline-formula><mml:math id="M89" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mn mathvariant="normal">10</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, and <inline-formula><mml:math id="M90" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mn mathvariant="normal">60</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> can be determined
(Fig. <xref ref-type="fig" rid="Ch1.F8"/>a, b). We note that <inline-formula><mml:math id="M91" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi mathvariant="normal">mean</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>
does not coincide with the RTD peak, which only happens if the RTD shows a
perfect symmetry (in log scale), which is not the case here.</p>
</sec>
<sec id="Ch1.S4.SS6">
  <title>Images of spectral parameters obtained from Debye decomposition</title>
      <p>Images depicting the DD-derived total polarization
(<inline-formula><mml:math id="M92" display="inline"><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula>) results of the complex conductivity spectra
(obtained from multi-frequency EIT) for selected time steps are presented in
Fig. <xref ref-type="fig" rid="Ch1.F9"/>. The extension of the root system against the
surrounding water (characterized by low polarization) is clearly delineated
in the images, and a continuous decrease in polarization strength is observed
over time.</p>
      <p>For further analysis, the complex conductivity spectra (also referred to as
pixel spectra) were classified into two categories, with and without root
segments, as described in Sect. <xref ref-type="sec" rid="Ch1.S3.SS4.SSS4"/>. The resulting
“root system spectra” were then processed separately, and care was taken
that the selected spectra exhibit a sufficiently strong and consistent
polarization response to allow a reliable relaxation time analysis.
Figure <xref ref-type="fig" rid="Ch1.F10"/> shows the comparison of the
<inline-formula><mml:math id="M93" display="inline"><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula> results for the first time step with the
extension of the root system according to the photograph. The root system
area reconstructed from the spectral EIT results shows a good agreement with
the known outer boundaries of the root system. Systematic changes in the
overall root system response were analyzed by averaging the
<inline-formula><mml:math id="M94" display="inline"><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula> pixel values in the root system zone
(Fig. <xref ref-type="fig" rid="Ch1.F11"/>). This average polarization response shows a
steady decrease over time.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F8"><caption><p>Debye decomposition (DD) of the recovered complex conductivity spectrum
for a pixel from the stem area with maximum polarization response (see.
Figs. <xref ref-type="fig" rid="Ch1.F5"/> and <xref ref-type="fig" rid="Ch1.F7"/>):
<bold>(a)</bold> complex conductivity (gray: real component; black: imaginary
component) from spectral EIT (dots) and fitted DD response (solid curves);
<bold>(b)</bold> corresponding relaxation time distribution. Vertical gray solid
lines indicate <inline-formula><mml:math id="M95" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi mathvariant="normal">mean</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, and the dashed vertical lines indicate
<inline-formula><mml:math id="M96" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mn mathvariant="normal">10</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M97" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mn mathvariant="normal">60</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>, respectively.</p></caption>
          <?xmltex \igopts{width=184.942913pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f08.pdf"/>

        </fig>

      <p>Images of the DD-derived mean relaxation time <inline-formula><mml:math id="M98" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi mathvariant="normal">mean</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are
presented in Fig. <xref ref-type="fig" rid="Ch1.F12"/> for selected time steps. Spatial
variations within the root system zone can be observed for each time step, as
well as changes between time steps. A general trend from larger
relaxation times (up to 18 ms) to smaller relaxation times (down to 9 ms)
over the course of the experiment is noticeable. Corresponding images of the uniformity
parameter <inline-formula><mml:math id="M99" display="inline"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mn mathvariant="normal">60</mml:mn><mml:mo>,</mml:mo><mml:mn mathvariant="normal">10</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> are shown in Fig. <xref ref-type="fig" rid="Ch1.F13"/>. Observed
variations within images and between time steps indicate changes in the shape
of the pixel spectra. Values approaching 1 indicate a stronger spectral
dispersion, i.e., a focusing of the spectral polarization response in a
narrower frequency band.</p>
</sec>
</sec>
<sec id="Ch1.S5">
  <title>Discussion</title>
      <p>The following discussion is divided into two parts: the biological discussion
of the experiment and the assessment of the geophysical methodology for crop
root investigations.</p>
<sec id="Ch1.S5.SS1">
  <title>Biological interpretation</title>
      <p>In the course of the experiment, a conductivity increase is observed
(Fig. <xref ref-type="fig" rid="Ch1.F7"/>a–c), which originated at the bottom of the
rhizotron and continuously migrated upwards (see S1 in the supplement for
conductivity images of the whole rhizotron). This spatial pattern rules out
the root system as the cause of the conductivity increase, and we attribute
it to the dissolution and subsequent upwards diffusion of impurities at the
bottom of the rhizotron frame. While a certain influence on the spectral
parameter results (Figs. <xref ref-type="fig" rid="Ch1.F9"/> and
<xref ref-type="fig" rid="Ch1.F12"/>) is possible, we believe the impact to be
relatively small for two reasons: First, the observed time evolution of
<inline-formula><mml:math id="M100" display="inline"><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula> shows changes more or less centered around
the stem region, not following the distribution pattern of the conductivity
increase; second, we observe changes in the spectral behavior of the
signatures (in terms of <inline-formula><mml:math id="M101" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi mathvariant="normal">mean</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>), which can not be explained
by an increase in conductivity. However, in future experiments the background
conditions should be monitored and the plant only inserted once equilibrium
of the system has been reached.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F9"><caption><p>Spatial distribution of DD-derived parameter
<inline-formula><mml:math id="M102" display="inline"><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula> for selected time steps. The top boundary is
adjusted according to the estimated water table for each measurement time.</p></caption>
          <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f09.png"/>

        </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F10"><caption><p>Comparison of root extension inferred from the
photograph <bold>(a)</bold>, indicated by overlaid solid lines, and <bold>(b)</bold>
<inline-formula><mml:math id="M103" display="inline"><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula> results for the first time step.
In <bold>(b)</bold> only pixels from the root zone, i.e., pixels with
a polarization response above the identified <inline-formula><mml:math id="M104" display="inline"><mml:mrow><mml:msup><mml:mi mathvariant="italic">σ</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:math></inline-formula> threshold value, are plotted.</p></caption>
          <?xmltex \igopts{width=184.942913pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f10.pdf"/>

        </fig>

      <p>The physiological response of the root system to the imposed nutrient
deprivation and possible anaerobic conditions is reflected by a decreasing
overall polarization response (Figs. <xref ref-type="fig" rid="Ch1.F9"/> and
<xref ref-type="fig" rid="Ch1.F11"/>). Note that it is highly unlikely that the dropping
water table caused this decrease in polarization, as more current is forced
through the main bulk of the root system with the dropping water table. Thus,
it should actually increase the polarization response, assuming that this
response does not change due to physiological reactions in the root system.
Note that the water table always remains above the actual root system and
only drops in the stem area. We attribute the observed decrease in
polarization to a general weakening of the EDLs present in the root system.
The cause of this EDL weakening can be manifold and can not be isolated in
this study. In the following, we shortly discuss two (possibly superimposing)
approaches to interpretation.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F11"><caption><p>Mean value of DD-derived total chargeability
<inline-formula><mml:math id="M105" display="inline"><mml:mrow><mml:msubsup><mml:mi>m</mml:mi><mml:mi mathvariant="normal">tot</mml:mi><mml:mi mathvariant="normal">n</mml:mi></mml:msubsup></mml:mrow></mml:math></inline-formula> plotted vs. time after start of the
experiment. Average values were computed based on all pixels belonging to the
root system zone.</p></caption>
          <?xmltex \igopts{width=184.942913pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f11.png"/>

        </fig>

      <p>The first approach is to consider ion (nutrient) concentration in the fluid
phase of the EDL. Shoot–root systems
represent a hydraulically connected system whose water potential is primarily
controlled by water transpiration at the leaves. In the case of intact
hydraulic connectivity in the plant, a decrease in water potential due to
transpiration causes water and solute uptake by the roots, and water and
solute flow from the roots to the leaves <xref ref-type="bibr" rid="bib1.bibx88" id="paren.90"/>. While
it is possible that solutes were taken up by the plant in our experiment, no
nutrients were available to the plant. Accordingly, for our experiment we
expected that the root regions farther away from the stem, that is, the
regions with the highest water potentials, were depleted of nutrients first,
as the available nutrients were translocated to the stem and the leaf areas,
following the (negative) water potential gradient. Plants can sense, and
react to, stress situations and can initiate changes in solute transport and
hydraulic membrane conductivity properties
<xref ref-type="bibr" rid="bib1.bibx24 bib1.bibx82" id="paren.91"><named-content content-type="pre">e.g.,</named-content></xref>, especially if
dynamic responses are considered. As a result, without being able to pinpoint
the exact cause, the ion concentration could have decreased in the vicinity
of cell membranes in these depleted areas, leading to a weakening of
associated EDLs, in turn implying a decreased polarization response. The
second approach is to consider changes in the electrical surface
characteristics of cell membranes in reaction to the imposed physiological
stress situation. These surface characteristics could have changed in
response to certain active triggers, such as stress hormones, or as a result
of balancing processes across the membranes. In light of these
considerations, the stem should retain a more stable electrical polarization
since it is likely not affected as much by physiological responses as the
other, smaller, root segments (in terms of larger nutrient storage
capability, better air availability, and general metabolic activity). This is
consistent with the time-lapse imaging results, which show more stable
polarization responses in the stem area (Fig. <xref ref-type="fig" rid="Ch1.F9"/>).</p>
      <p>Another indication of the physiological stress response is the changes in
the shape of the spectra (Figs. <xref ref-type="fig" rid="Ch1.F12"/> and
<xref ref-type="fig" rid="Ch1.F13"/>). Relative changes in the relaxation time
contributions suggest changes in the underlying structures that control the
polarization response at certain time steps. These changes might be related
to new or ceased ion fluxes and their varying pathways within the root
systems, as well as to varying surface charges at various structures such as
the endodermis. If these structures change, or break down, in response to
stress situations, corresponding changes in the electrical properties can be
expected. However, given its spatial resolution limits, EIT does not allow one to
distinguish these different structures directly.</p>
      <p>Assuming that relaxation times can be linked to length scales of the
underlying polarization processes, the observed signatures indicate multiple
polarizable structures. However, the methodology applied here prevents
further investigations in this direction. In contrast to most of the existing
studies, we did not inject current directly in the stem, and correspondingly
the explicit current pathways are much less defined in our approach. This
prevents (at this stage) a simple formulation of an equivalent lumped
electrical circuit model. Comparison of measurements using the procedure
presented here with a stem-injection approach could, however, help to
elucidate the origin of polarization and its length-scale characteristics.
Current injection into the stem forces the current to flow through the root
system and through all radial layers of the root segments, and thus a
stronger polarization response from inside the root segment can be expected,
as well as the polarization of additional membrane structures. Additional
experiments could focus on establishing relationships between recovered
spectral polarization parameters and root-specific parameters, such as
surface area and root length density. The use of sophisticated electrical
models, coupled to existing macroscopic root development and nutrient uptake
models <xref ref-type="bibr" rid="bib1.bibx31 bib1.bibx45" id="paren.92"><named-content content-type="pre">e.g.,</named-content></xref>, could
provide further insight to identify the key processes that control the
electrical polarization signatures, of root systems. While this study focused
on establishing the EIT methodology for crop root research, in future
studies, physiological plant parameters such as leaf transpiration rates,
which could help to identify the key processes that control electrical
polarization signatures, should also
be monitored.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F12"><caption><p>Spatial distribution of the DD-derived parameter,
<inline-formula><mml:math id="M106" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">τ</mml:mi><mml:mi mathvariant="normal">mean</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, for selected time steps. Only pixels belonging to the
root system zone are plotted. Masked (white) pixels were classified as water.
The top boundary is adjusted according to the estimated water table for each
measurement time.</p></caption>
          <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f12.png"/>

        </fig>

      <?xmltex \floatpos{t}?><fig id="Ch1.F13"><caption><p>Spatial distribution of the DD-derived parameter, <inline-formula><mml:math id="M107" display="inline"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mrow><mml:mn mathvariant="normal">60</mml:mn><mml:mo>,</mml:mo><mml:mn mathvariant="normal">10</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>, for
selected time steps. Only pixels belonging to the root system zone are
plotted. Masked (white) pixels were classified as water. The top boundary is
adjusted according to the estimated water table for each measurement time.</p></caption>
          <?xmltex \igopts{width=241.848425pt}?><graphic xlink:href="https://bg.copernicus.org/articles/14/921/2017/bg-14-921-2017-f13.png"/>

        </fig>

      <p>As already pointed out by <xref ref-type="bibr" rid="bib1.bibx76" id="text.93"/>, single-frequency
measurements are of limited value to determine electrical polarization
properties of root systems, both in terms of spatial distribution and polarization strength, and this
finding is supported by our spectral EIS results
(Figs. <xref ref-type="fig" rid="Ch1.F5"/>, <xref ref-type="fig" rid="Ch1.F7"/>,
<xref ref-type="fig" rid="Ch1.F12"/>, and <xref ref-type="fig" rid="Ch1.F13"/>). This becomes even more
obvious when interpreting the polarization signature as an EDL response,
which typically exhibits a strong frequency dependence. It should be noted
that the frequency range analyzed here in an imaging framework (0.46 to
220 Hz) does not cover the full bandwidth that could be, in principle,
measured with the presented setup, and corresponding advances are within easy
reach <xref ref-type="bibr" rid="bib1.bibx44" id="paren.94"><named-content content-type="pre">e.g.,</named-content></xref>. However, reliable measurements
at lower and higher frequencies will require careful adaptations in
measurement and data processing procedures.</p>
      <p>The classification of image pixels into two classes cannot, and should not,
be treated as a universal analysis procedure. For the simple conditions in
this experiment a clear distinction between the root system area and the surrounding
medium could be made, which facilitated the assessment of the method (e.g.,
Fig. <xref ref-type="fig" rid="Ch1.F11"/>), and can potentially be used for further
experiments with root systems in aqueous solutions. However, the primary
results of this study do not rely on this specific classification, and
likewise soil-based experiments could be conducted with the measurement
setup.</p>
      <p>This study does not involve any kind of granular substrate and thus excludes
possible influences from such a background material. In fact, significant
additional electrical polarization can be expected when soil surrounds the
root system, which will superimpose on the root system response. Organic
matter and micorrhiza may also contribute to the overall electrical
signature. Finally, a variable water content can significantly influence the
electrical response of the soil and the root system, either directly by
influencing present EDLs or indirectly by inducing physiological processes
such as nutrient uptake, which in turn can affect the electrical signatures
of the EDLs.</p>
</sec>
<sec id="Ch1.S5.SS2">
  <title>Geophysical methodology</title>
      <p>The observed polarization response of the root system is relatively weak and
its measurement requires a corresponding measurement instrument accuracy. This accuracy is provided by the EIT-40 tomograph that was used
in this study <xref ref-type="bibr" rid="bib1.bibx102" id="paren.95"/>. The high accuracy of the
instrument was also recently demonstrated in an imaging study on soil columns
<xref ref-type="bibr" rid="bib1.bibx47" id="paren.96"/>.</p>
      <p>If fixed data weighting is used, which we believe would produce more reliable
and consistent results for multi-frequency time-lapse data, data selection,
i.e., filtering, becomes a relevant step in the processing pipeline before
the inversion and subsequent spectral analysis. While it is common to remove
outliers from geophysical data prior to inversion, filtering becomes
challenging if multiple time steps are to be analyzed consistently.
The number of retained data points varied slightly between time steps,
although the same filtering criteria were applied. This can be explained by
data noise and varying contact impedances at the electrodes. However, data
quality was sufficient enough to produce consistent imaging results for all
time steps and frequencies, as is evident from the impedance spectra
(Fig. <xref ref-type="fig" rid="Ch1.F5"/>).</p>
      <p>Another important issue is the data processing flow in the imaging framework,
coupled with the spectral analysis based on the Debye decomposition. The
inversion algorithm produces spatially smooth images; however, the images
were computed for each frequency separately, and thus no smooth variation
between adjacent frequencies is enforced in the inversion, although
physically expected. Corresponding inversion algorithms have been developed
recently <xref ref-type="bibr" rid="bib1.bibx50 bib1.bibx40" id="paren.97"/> and could lead to a
further improvement of the multi-frequency imaging results. However, a
similar constraint is introduced by the Debye decomposition, where smoothness
is imposed along the relaxation time axis, which directly corresponds to the
frequency axis. Minor noise components can thus be expected to be smoothed
out both spatially and spectrally.</p>
      <p>EIT applications in pseudo-2-D rhizotron containers require specific
processing steps. The determination and testing of correction factors
accounting for modeling errors due to an imperfect 2-D situation
(Fig. <xref ref-type="fig" rid="Ch1.F2"/>) are as important as the correct
representation of the rhizotron in terms of the FE grid underlying the
inversion process (Fig. <xref ref-type="fig" rid="Ch1.F3"/>). Not taking these aspects
into account can produce artifacts in the imaging results that can easily be
misinterpreted in biological terms. Similarly, data errors should not be
underestimated, as this can also produce misleading imaging results when data
are overfitted <xref ref-type="bibr" rid="bib1.bibx49" id="paren.98"><named-content content-type="pre">e.g.,</named-content></xref>. Contrarily, an
overestimation of data errors can mask information present in the data. Among
other analysis steps, raw (impedance) data and imaging (complex conductivity) data
should be checked for consistency and plausibility by taking into account the
much lower spatial resolution of the raw data (see
Figs. <xref ref-type="fig" rid="Ch1.F5"/> and <xref ref-type="fig" rid="Ch1.F7"/>).</p>
      <p>Electrical imaging results exhibit a spatially variable resolution, which
usually decreases as the distance from the electrodes increases. One could
question such a method's usefulness if the resolution cannot be clearly
determined. Nonetheless, even limited spatial information allows for a
distinction of polarizing and non-polarizing regions in the investigated
object. This is not possible with spectroscopic measurements, and
so it is even more difficult to analyze spatially distributed root systems with measurements such as these.  We suspect that some of the reported
inconsistencies in electrical capacitance relationships with biological
parameters <xref ref-type="bibr" rid="bib1.bibx54" id="paren.99"><named-content content-type="pre">e.g.,</named-content><named-content content-type="post">and references therein</named-content></xref> can be
ascribed to missing spatial information in the measurement data. The resolution
of EIT is not sufficient to image microscopic current flow paths in the root
system, but the imaged macroscopic electrical properties can be compared for
different regions of the root system, for instance the (older) top part of the
root system compared to the younger lower part. Future improvements in
experimental setups (electrode distribution and spacing) and measurement
configurations will most probably lead to increased spatial resolution.</p>
      <p>Another advantage of the EIT approach presented in this study is the
possibility of arbitrarily placed electrodes (as long as the resulting
geometrical arrangement allows for a sufficient measurement coverage of the
root system), in contrast to using stem electrodes as commonly done in previous
studies.  If the stem of a plant is used to inject current into the root
system, measurements, and resulting correlations to biological parameters, are
highly sensitive to the electrode position above the stem base
<xref ref-type="bibr" rid="bib1.bibx27 bib1.bibx70" id="paren.100"/>. Another problem is that
electrodes can not be inserted into the stem if damage to the plant is to be
avoided. However, injections can also be realized by use of non-invasive
clamps.</p>
      <p>The timescale of the physiological response to be monitored is also important
for the experimental design. It took approximately 3.5 h to complete a
single time frame of the spectral EIT measurements presented here. Thus,
physiological processes taking place on a shorter time span cannot be
resolved. Reducing the data acquisition time can be achieved by reducing
either the number of low-frequency measurements or the number of current
injections. This can result in a loss of spectral and spatial resolution if
measurement configurations are not suitably optimized to compensate for the
lost number of measurements.</p>
</sec>
</sec>
<sec id="Ch1.S6" sec-type="conclusions">
  <title>Conclusions</title>
      <p>The goal of this study was to investigate and establish spectral (i.e.,
multifrequency) EIT as a non-invasive tool for the characterization and
monitoring of crop root systems. Based on working hypotheses derived from the
state of science in the involved fields, including geophysics and plant
science, we designed and conducted a controlled experiment in which the root
systems of oilseed plants were monitored in a 2-D, water-filled rhizotron
container. Since water does not exhibit a significant polarization response
in the considered frequency range, the observed electrical polarization
response could be attributed to the root systems.</p>
      <p>The spectral EIT imaging results revealed a low-frequency polarization
response of the root system, which enabled the successful delineation of the
spatial extension of the root system. Based on a pixel-based Debye
decomposition analysis of the spectral imaging results, we found a mean
relaxation time of the root system's polarization signature in the covered
frequency range of the order of 10 ms, corresponding to a frequency of the
order of 15 Hz. Importantly, upon ongoing nutrient deprivation (with
possibly anaerobic conditions), the magnitude of the overall polarization
response steadily decreased and the spectral characteristics systematically
changed, indicating changes in the length scales on which the polarization
processes took place in the root system. The spectral EIT imaging results
could be explained by the macroscopically observed and expected physiological
response of the plant to the imposed stress situation. The identification of
the root structures and processes controlling the root electrical signatures,
however, was beyond the scope of this study given the inherent spatial
resolution limits of EIT. Nonetheless the recovered electrical signatures
could be used in the future to develop and calibrate improved macroscopic
root electrical models which incorporate microscopic processes.</p>
      <p>We showed, for the first time (to the best of our knowledge), that spectral
EIT is a capable, non-invasive method to image root system extension as well
as to monitor changes associated with root physiological processes. Given the
applicability of the method at both the laboratory and field scale, our results
suggest an enormous potential of spectral EIT for the structural and
functional imaging of root systems for various applications. In particular, at
the field scale, non-invasive methods for root system characterization and
imaging are lacking and EIT seems to be a very promising method to fill this
gap. In future studies we will aim to further prove the suitability of
spectral EIT to monitor physiological responses in different situations and
to different stimuli, at both laboratory and field scales.</p>
</sec>
<sec id="Ch1.S7">
  <title>Data availability</title>
      <p>Measured raw data, electrical imaging data, spectral analysis results, and
Python scripts used to generate the figures can be accessed at
<ext-link xlink:href="http://dx.doi.org/10.5281/zenodo.260087" ext-link-type="DOI">10.5281/zenodo.260087</ext-link> (<xref ref-type="bibr" rid="bib1.bibx97" id="altparen.101"/>).</p>
</sec>

      
      </body>
    <back><app-group>
        <supplementary-material position="anchor"><p><bold>The Supplement related to this article is available online at <inline-supplementary-material xlink:href="http://dx.doi.org/10.5194/bg-14-921-2017-supplement" xlink:title="pdf">doi:10.5194/bg-14-921-2017-supplement</inline-supplementary-material>.</bold></p></supplementary-material>
        </app-group><notes notes-type="competinginterests">

      <p>The authors declare that they have no conflict of
interest.</p>
  </notes><ack><title>Acknowledgements</title><p>Parts of this work were funded by the Deutsche Forschungsgemeinschaft (DFG)
in the framework of the project “Non-destructive characterization and
monitoring of root structure and function at the rhizotron and field scale
using spectral electrical impedance tomography” (KE 1138/1-1) and the
collaborative research center “Patterns in soil-vegetation-atmosphere
systems: monitoring, modeling and data assimilation” (SFB/TR 32). We are
especially grateful to Egon Zimmermann and Matthias Kelter for valuable
discussions regarding the measurement setup. We also thank Achim Walter,
Johannes Pfeifer and Kerstin Nagel for technical support and discussions on
root physiology in the initial phase of the work. The authors would also like
to thank the three anonymous referees for their very constructive reviews,
which lead to substantial improvements of the
paper.
<?xmltex \hack{\newline}?><?xmltex \hack{\newline}?> Edited by: P. Stoy<?xmltex \hack{\newline}?> Reviewed by: three
anonymous referees</p></ack><ref-list>
    <title>References</title>

      <ref id="bib1.bibx1"><label>al Hagrey(2007)</label><mixed-citation>al Hagrey, S.: Geophysical imaging of root-zone, trunk, and moisture
heterogeneity, J. Exp. Bot., 58, 839–854,
<ext-link xlink:href="http://dx.doi.org/10.1093/jxb/erl237" ext-link-type="DOI">10.1093/jxb/erl237</ext-link>, 2007.</mixed-citation></ref>
      <ref id="bib1.bibx2"><label>al Hagrey and Petersen(2011)</label><mixed-citation>al Hagrey, S. and Petersen, T.: Numerical and experimental mapping of small
root zones using optimized surface and borehole resistivity tomography,
Geophysics, 76, G25–G35, <ext-link xlink:href="http://dx.doi.org/10.1190/1.3545067" ext-link-type="DOI">10.1190/1.3545067</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bibx3"><label>Alumbaugh and Newman(2000)</label><mixed-citation>Alumbaugh, D. and Newman, G.: Image appraisal for 2-D and 3-D
electromagnetic inversion, Geophysics, 65, 1455–1467,
<ext-link xlink:href="http://dx.doi.org/10.1190/1.1444834" ext-link-type="DOI">10.1190/1.1444834</ext-link>, 2000.</mixed-citation></ref>
      <ref id="bib1.bibx4"><label>Amato et al.(2008)Amato, Basso, Celano, Bitella, Morelli, and
Rossi</label><mixed-citation>Amato, M., Basso, B., Celano, G., Bitella, G., Morelli, G., and Rossi, R.: In
situ detection of tree root distribution and biomass by multi-electrode
resistivity imaging, Tree Physiol., 28, 1441–1448,
<ext-link xlink:href="http://dx.doi.org/10.1093/treephys/28.10.1441" ext-link-type="DOI">10.1093/treephys/28.10.1441</ext-link>, 2008.</mixed-citation></ref>
      <ref id="bib1.bibx5"><label>Amato et al.(2009)Amato, Bitella, Rossi, Gómez, Lovelli, and
Gomes</label><mixed-citation>Amato, M., Bitella, G., Rossi, R., Gómez, J. A., Lovelli, S., and Gomes,
J. J. F.: Multi-electrode 3D resistivity imaging of alfalfa root zone, Eur.
J. Agron., 31, 213–222, <ext-link xlink:href="http://dx.doi.org/10.1016/j.eja.2009.08.005" ext-link-type="DOI">10.1016/j.eja.2009.08.005</ext-link>,
2009.</mixed-citation></ref>
      <ref id="bib1.bibx6"><label>Anderson and Hopmans(2013)</label><mixed-citation>Anderson, S. and Hopmans, J. (Eds.): Soil-water-root Processes: Advances in
Tomography and Imaging, Soil Sci. Soc. Am.,
<ext-link xlink:href="http://dx.doi.org/10.2136/sssaspecpub61" ext-link-type="DOI">10.2136/sssaspecpub61</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx7"><label>Anderson and Higinbotham(1976)</label><mixed-citation>Anderson, W. and Higinbotham, N.: Electrical resistances of corn root
segments,
Plant Physiol., 57, 137–141, <ext-link xlink:href="http://dx.doi.org/10.1104/pp.57.2.137" ext-link-type="DOI">10.1104/pp.57.2.137</ext-link>, 1976.</mixed-citation></ref>
      <ref id="bib1.bibx8"><label>Aubrecht et al.(2006)Aubrecht, Staněk, and
Koller</label><mixed-citation>Aubrecht, L., Staněk, Z., and Koller, J.: Electrical measurement of the
absorption surfaces of tree roots by the earth impedance method: 1. Theory,
Tree Physiol., 26, 1105–1112, <ext-link xlink:href="http://dx.doi.org/10.1093/treephys/26.9.1105" ext-link-type="DOI">10.1093/treephys/26.9.1105</ext-link>, 2006.</mixed-citation></ref>
      <ref id="bib1.bibx9"><label>Aulen and Shipley(2012)</label><mixed-citation>Aulen, M. and Shipley, B.: Non-destructive estimation of root mass using
electrical capacitance on ten herbaceous species, Plant Soil, 355,
41–49, <ext-link xlink:href="http://dx.doi.org/10.1007/s11104-011-1077-3" ext-link-type="DOI">10.1007/s11104-011-1077-3</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bibx10"><label>Barsoukov and Macdonald(2005)</label><mixed-citation>Barsoukov, E. and Macdonald, J. (Eds.): Impedance spectroscopy: theory,
experiment, and applications, Wiley-Interscience, <ext-link xlink:href="http://dx.doi.org/10.1002/0471716243" ext-link-type="DOI">10.1002/0471716243</ext-link>,
2005.</mixed-citation></ref>
      <ref id="bib1.bibx11"><label>Bayford(2006)</label><mixed-citation>Bayford, R.: Bioimpedance tomography (electrical impedance tomography),
Annual
Reviews Biomedical Engineering, 8, 63–91,
<ext-link xlink:href="http://dx.doi.org/10.1146/annurev.bioeng.8.061505.095716" ext-link-type="DOI">10.1146/annurev.bioeng.8.061505.095716</ext-link>, 2006.</mixed-citation></ref>
      <ref id="bib1.bibx12"><label>Beff et al.(2013)Beff, Günther, Vandoorne, Couvreur, and
Javaux</label><mixed-citation>Beff, L., Günther, T., Vandoorne, B., Couvreur, V., and Javaux, M.:
Three-dimensional monitoring of soil water content in a maize field using
Electrical Resistivity Tomography, Hydrol. Earth Syst. Sci., 17, 595–609,
<ext-link xlink:href="http://dx.doi.org/10.5194/hess-17-595-2013" ext-link-type="DOI">10.5194/hess-17-595-2013</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx13"><label>Benlloch-González et al.(2010)Benlloch-González, Fournier,
and Benlloch</label><mixed-citation>Benlloch-González, M., Fournier, J., and Benlloch, M.: K<inline-formula><mml:math id="M108" display="inline"><mml:msup><mml:mi/><mml:mo>+</mml:mo></mml:msup></mml:math></inline-formula> deprivation
induces xylem water and K<inline-formula><mml:math id="M109" display="inline"><mml:msup><mml:mi/><mml:mo>+</mml:mo></mml:msup></mml:math></inline-formula> transport in sunflower: evidence for a
co-ordinated control, J. Exp. Bot., 61, 157–164,
<ext-link xlink:href="http://dx.doi.org/10.1093/jxb/erp288" ext-link-type="DOI">10.1093/jxb/erp288</ext-link>, 2010.</mixed-citation></ref>
      <ref id="bib1.bibx14"><label>Binley and Kemna(2005)</label><mixed-citation>Binley, A. and Kemna, A.: DC resistivity and induced polarization methods,
in:
Hydrogeophysics, edited by: Rubin, Y. and Hubbard, S. S.,
Springer, the
Netherlands, 129–156, <ext-link xlink:href="http://dx.doi.org/10.1007/1-4020-3102-5_5" ext-link-type="DOI">10.1007/1-4020-3102-5_5</ext-link>, 2005.</mixed-citation></ref>
      <ref id="bib1.bibx15"><label>Binley et al.(2005)Binley, Slater, Fukes, and
Cassiani</label><mixed-citation>Binley, A., Slater, L., Fukes, M., and Cassiani, G.: Relationship between
spectral induced polarization and hydraulic properties of saturated and
unsaturated sandstone, Water Resour. Res., 41, 1–13,
<ext-link xlink:href="http://dx.doi.org/10.1029/2005WR004202" ext-link-type="DOI">10.1029/2005WR004202</ext-link>, 2005.</mixed-citation></ref>
      <ref id="bib1.bibx16"><label>Boaga et al.(2013)Boaga, Rossi, and Cassiani</label><mixed-citation>Boaga, J., Rossi, M., and Cassiani, G.: Monitoring soil-plant interactions in
an apple orchard using 3D electrical resistivity tomography, Procedia
Environmental Sciences, 19, 394–402, <ext-link xlink:href="http://dx.doi.org/10.1016/j.proenv.2013.06.045" ext-link-type="DOI">10.1016/j.proenv.2013.06.045</ext-link>,
2013.</mixed-citation></ref>
      <ref id="bib1.bibx17"><label>Bücker and Hördt(2013)</label><mixed-citation>Bücker, M. and Hördt, A.: Analytical modelling of membrane
polarization
with explicit parametrization of pore radii and the electrical double layer,
Geophys. J. Int., 194, 804–813, <ext-link xlink:href="http://dx.doi.org/10.1093/gji/ggt136" ext-link-type="DOI">10.1093/gji/ggt136</ext-link>,
2013.</mixed-citation></ref>
      <ref id="bib1.bibx18"><label>Cao et al.(2010)Cao, Repo, Silvennoinen, Lehto, and
Pelkonen</label><mixed-citation>Cao, Y., Repo, T., Silvennoinen, R., Lehto, T., and Pelkonen, P.: An
appraisal
of the electrical resistance method for assessing root surface area, J. Exp. Bot., 61, 2491–2497, <ext-link xlink:href="http://dx.doi.org/10.1093/jxb/erq078" ext-link-type="DOI">10.1093/jxb/erq078</ext-link>, 2010.</mixed-citation></ref>
      <ref id="bib1.bibx19"><label>Cao et al.(2011)Cao, Repo, Silvennoinen, Lehto, and
Pelkonen</label><mixed-citation>Cao, Y., Repo, T., Silvennoinen, R., Lehto, T., and Pelkonen, P.: Analysis of
the willow root system by electrical impedance spectroscopy, J.
Exp. Bot., 62, 351–358, <ext-link xlink:href="http://dx.doi.org/10.1093/jxb/erq276" ext-link-type="DOI">10.1093/jxb/erq276</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bibx20"><label>Čermák et al.(2006)</label><mixed-citation>Čermák, J., Ulrich, R., Staněk, Z., Koller, J., and Aubrecht,
L.: Electrical measurement of tree root absorbing surfaces by the earth
impedance method: 2. Verification based on allometric relationships and
root severing experiments, Tree Physiol., 26, 1113–1121,
<ext-link xlink:href="http://dx.doi.org/10.1093/treephys/26.9.1113" ext-link-type="DOI">10.1093/treephys/26.9.1113</ext-link>, 2006.</mixed-citation></ref>
      <ref id="bib1.bibx21"><label>Chloupek(1972)</label><mixed-citation>Chloupek, O.: The relationship between electric capacitance and some other
parameters of plant roots, Biol. Plantarum, 14, 227–230,
<ext-link xlink:href="http://dx.doi.org/10.1007/BF02921255" ext-link-type="DOI">10.1007/BF02921255</ext-link>, 1972.</mixed-citation></ref>
      <ref id="bib1.bibx22"><label>Chloupek(1976)</label><mixed-citation>Chloupek, O.: Die Bewertung des Wurzelsystems von Senfpflanzen auf Grund der
dielektrischen Eigenschaften und mit Rücksicht auf den Endertrag,
Biol. Plantarum, 18, 44–49, <ext-link xlink:href="http://dx.doi.org/10.1007/BF02922333" ext-link-type="DOI">10.1007/BF02922333</ext-link>, 1976.</mixed-citation></ref>
      <ref id="bib1.bibx23"><label>Claassen and Barber(1974)</label><mixed-citation>Claassen, N. and Barber, S.: A method for characterizing the relation between
nutrient concentration and flux into roots of intact plants, Plant
Physiol., 54, 564–568, <ext-link xlink:href="http://dx.doi.org/10.1104/pp.54.4.564" ext-link-type="DOI">10.1104/pp.54.4.564</ext-link>, 1974.</mixed-citation></ref>
      <ref id="bib1.bibx24"><label>Clarkson et al.(2000)Clarkson, Carvajal, Henzler, Waterhouse, Smyth,
Cooke, and Steudle</label><mixed-citation>Clarkson, D., Carvajal, M., Henzler, T., Waterhouse, R., Smyth, A., Cooke,
D.,
and Steudle, E.: Root hydraulic conductance: diurnal aquaporin expression and
the effects of nutrient stress, J. Exp. Bot., 51, 61–70,
<ext-link xlink:href="http://dx.doi.org/10.1093/jexbot/51.342.61" ext-link-type="DOI">10.1093/jexbot/51.342.61</ext-link>, 2000.</mixed-citation></ref>
      <ref id="bib1.bibx25"><label>Cseresnyés et al.(2013)Cseresnyés, Tünde, Végh,
Anton, and Rajkai</label><mixed-citation>Cseresnyés, I., Tünde, T., Végh, K., Anton, A., and Rajkai, K.:
Electrical impedance and capacitance method: A new approach for detection of
functional aspects of arbuscular mycorrhizal colonization in maize, European
J. Soil Biol., 54, 25–31, <ext-link xlink:href="http://dx.doi.org/10.1016/j.ejsobi.2012.11.001" ext-link-type="DOI">10.1016/j.ejsobi.2012.11.001</ext-link>,
2013.</mixed-citation></ref>
      <ref id="bib1.bibx26"><label>Daily et al.(2005)Daily, Ramirez, Binley, and
LaBrecque</label><mixed-citation>Daily, W., Ramirez, A., Binley, A., and LaBrecque, D.: Electrical resistance
tomography – theory and practice, in: Near-Surface Geophysics, edited by:
Butler, D., Vol. 154, chap. 17, Society of Exploration
Geophysicists, Tulsa, OK, 525–550, <ext-link xlink:href="http://dx.doi.org/10.1190/1.9781560801719.ch17" ext-link-type="DOI">10.1190/1.9781560801719.ch17</ext-link>, 2005.</mixed-citation></ref>
      <ref id="bib1.bibx27"><label>Dalton(1995)</label><mixed-citation>Dalton, F.: In-situ root extent measurements by electrical capacitance
methods,
Plant Soil, 173, 157–165, <ext-link xlink:href="http://dx.doi.org/10.1007/BF00155527" ext-link-type="DOI">10.1007/BF00155527</ext-link>, 1995.</mixed-citation></ref>
      <ref id="bib1.bibx28"><label>Delhon et al.(1995)Delhon, Gojon, Tillard, and
Passama</label><mixed-citation>Delhon, P., Gojon, A., Tillard, P., and Passama, L.: Diurnal regulation of
NO<inline-formula><mml:math id="M110" display="inline"><mml:mrow><mml:msubsup><mml:mi/><mml:mn mathvariant="normal">3</mml:mn><mml:mo>-</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula> uptake in soybean plants I. Changes in NO<inline-formula><mml:math id="M111" display="inline"><mml:mrow><mml:msubsup><mml:mi/><mml:mn mathvariant="normal">3</mml:mn><mml:mo>-</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula> influx, efflux, and N
utilization in the plant during the day/night cycle, J. Exp.
Bot., 46, 1585–1594, <ext-link xlink:href="http://dx.doi.org/10.1093/jxb/46.10.1585" ext-link-type="DOI">10.1093/jxb/46.10.1585</ext-link>, 1995.</mixed-citation></ref>
      <ref id="bib1.bibx29"><label>Dietrich et al.(2013)Dietrich, Bengough, Jones, and
White</label><mixed-citation>Dietrich, R., Bengough, A., Jones, H., and White, P.: Can root electrical
capacitance be used to predict root mass in soil?, Ann. Bot., 112,
457–464, <ext-link xlink:href="http://dx.doi.org/10.1093/aob/mct044" ext-link-type="DOI">10.1093/aob/mct044</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx30"><label>Dukhin and Shilov(1974)</label><mixed-citation>
Dukhin, S. S. and Shilov, V. N.: Dielectric Phenomena and the Double Layer in
Disperse Systems and Polyelectrolytes, Wiley, New York, 1974.</mixed-citation></ref>
      <ref id="bib1.bibx31"><label>Dunbabin et al.(2013)Dunbabin, Postma, Schnepf, Pagès, Javaux,
Wu, Leitner, Chen, Rengel, and Diggle</label><mixed-citation>Dunbabin, V., Postma, J., Schnepf, A., Pagès, L., Javaux, M., Wu, L.,
Leitner, D., Chen, Y., Rengel, Z., and Diggle, A.: Modelling root-soil
interactions using three-dimensional models of root growth, architecture and
function, Plant Soil, 372, 93–124, <ext-link xlink:href="http://dx.doi.org/10.1007/s11104-013-1769-y" ext-link-type="DOI">10.1007/s11104-013-1769-y</ext-link>,
2013.</mixed-citation></ref>
      <ref id="bib1.bibx32"><label>Dvořák et al.(1981)</label><mixed-citation>Dvořák, M., Černohorská, J., and Janáček, K.:
Characteristics of current passage through plant tissue, Biol. Plantarum,
23, 306–310, <ext-link xlink:href="http://dx.doi.org/10.1007/BF02895374" ext-link-type="DOI">10.1007/BF02895374</ext-link>, 1981.</mixed-citation></ref>
      <ref id="bib1.bibx33"><label>Ellis et al.(2013)Ellis, Murray, and Kavalieris</label><mixed-citation>Ellis, T., Murray, W., and Kavalieris, L.: Electrical capacitance of bean
(Vicia faba) root systems was related to tissue density – a test for the
Dalton Model, Plant Soil, 366, 575–584,
<ext-link xlink:href="http://dx.doi.org/10.1007/s11104-012-1424-z" ext-link-type="DOI">10.1007/s11104-012-1424-z</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx34"><label>Fixman(1980)</label><mixed-citation>Fixman, M.: Charged Macromolecules in External Fields. 2. Preliminary
Remarks
on the Cylinder, Macromolecules, 13, 711–716, <ext-link xlink:href="http://dx.doi.org/10.1021/ma60075a043" ext-link-type="DOI">10.1021/ma60075a043</ext-link>,
1980.</mixed-citation></ref>
      <ref id="bib1.bibx35"><label>Flores Orozco et al.(2012a)Flores Orozco, Kemna,
Oberdörster, Zschornack, Leven, Dietrich, and
Weiss</label><mixed-citation>Flores Orozco, A., Kemna, A., Oberdörster, C., Zschornack, L., Leven, C.,
Dietrich, P., and Weiss, H.: Delineation of subsurface hydrocarbon
contamination at a former hydrogenation plant using spectral induced
polarization imaging, J. Contam. Hydrol., 136, 131–144,
<ext-link xlink:href="http://dx.doi.org/10.1016/j.jconhyd.2012.06.001" ext-link-type="DOI">10.1016/j.jconhyd.2012.06.001</ext-link>, 2012a.</mixed-citation></ref>
      <ref id="bib1.bibx36"><label>Flores Orozco et al.(2012b)Flores Orozco, Kemna, and
Zimmermann</label><mixed-citation>Flores Orozco, A., Kemna, A., and Zimmermann, E.: Data error quantification
in
spectral induced polarization imaging, Geophysics, 77, E227–E237,
<ext-link xlink:href="http://dx.doi.org/10.1190/geo2010-0194.1" ext-link-type="DOI">10.1190/geo2010-0194.1</ext-link>, 2012b.</mixed-citation></ref>
      <ref id="bib1.bibx37"><label>Flores Orozco et al.(2013)Flores Orozco, Williams, and
Kemna</label><mixed-citation>Flores Orozco, A., Williams, K., and Kemna, A.: Time-lapse spectral induced
polarization imaging of stimulated uranium bioremediation, Near Surf.
Geophys., 11, 531–544, <ext-link xlink:href="http://dx.doi.org/10.3997/1873-0604.2013020" ext-link-type="DOI">10.3997/1873-0604.2013020</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx38"><label>Friedel(2003)</label><mixed-citation>Friedel, S.: Resolution, stability and efficiency of resistivity tomography
estimated from a generalized inverse approach, Geophys. J.
Int., 153, 305–316, <ext-link xlink:href="http://dx.doi.org/10.1046/j.1365-246X.2003.01890.x" ext-link-type="DOI">10.1046/j.1365-246X.2003.01890.x</ext-link>, 2003.</mixed-citation></ref>
      <ref id="bib1.bibx39"><label>Gregory et al.(2003)Gregory, Hutchison, Read, Jenneson, Gilboy, and
Morton</label><mixed-citation>Gregory, P., Hutchison, D., Read, D., Jenneson, P., Gilboy, W., and Morton,
E.:
Non-invasive imaging of roots with high resolution X-ray micro-tomography,
in: Roots: The Dynamic Interface between Plants and the Earth,
Springer, 351–359, <ext-link xlink:href="http://dx.doi.org/10.1023/A:1026179919689" ext-link-type="DOI">10.1023/A:1026179919689</ext-link>, 2003.</mixed-citation></ref>
      <ref id="bib1.bibx40"><label>Günther and Martin(2016)</label><mixed-citation>Günther, T. and Martin, T.: Spectral two-dimensional inversion of
frequency-domain induced polarization data from a mining slag heap, J. Appl. Geophys., 135, 436–448,  <ext-link xlink:href="http://dx.doi.org/10.1016/j.jappgeo.2016.01.008" ext-link-type="DOI">10.1016/j.jappgeo.2016.01.008</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bibx41"><label>Heege(2013)</label><mixed-citation>
Heege, H.: Precision in Crop Farming: Site Specific Concepts and Sensing
Methods: Applications and Results, Springer Science &amp; Business Media, Springer, the
Netherlands,
2013.</mixed-citation></ref>
      <ref id="bib1.bibx42"><label>Heřmanská et al.(2015)</label><mixed-citation>Heřmanská, A., Středa, T., and Chloupek, O.: Improved wheat
grain yield by a new method of root selection, Agron. Sustain.
Dev., 35, 195–202, <ext-link xlink:href="http://dx.doi.org/10.1007/s13593-014-0227-4" ext-link-type="DOI">10.1007/s13593-014-0227-4</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bibx43"><label>Hose et al.(2001)Hose, Clarkson, Steudle, Schreiber, and
Hartung</label><mixed-citation>Hose, E., Clarkson, D., Steudle, E., Schreiber, L., and Hartung, W.: The
exodermis: a variable apoplastic barrier, J. Exp. Bot., 52,
2245–2264, <ext-link xlink:href="http://dx.doi.org/10.1093/jexbot/52.365.2245" ext-link-type="DOI">10.1093/jexbot/52.365.2245</ext-link>, 2001.</mixed-citation></ref>
      <ref id="bib1.bibx44"><label>Huisman et al.(2015)Huisman, Zimmermann, Esser, Haegel, Treichel, and
Vereecken</label><mixed-citation>Huisman, J., Zimmermann, E., Esser, O., Haegel, F., Treichel, A., and
Vereecken, H.: Evaluation of a novel correction procedure to remove
electrode impedance effects from broadband SIP measurements, J.
Appl. Geophys., 135, 466–473,  <ext-link xlink:href="http://dx.doi.org/10.1016/j.jappgeo.2015.11.008" ext-link-type="DOI">10.1016/j.jappgeo.2015.11.008</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bibx45"><label>Javaux et al.(2013)Javaux, Couvreur, Vanderborght, and
Vereecken</label><mixed-citation>Javaux, M., Couvreur, V., Vanderborght, J., and Vereecken, H.: Root water
uptake: From three-dimensional biophysical processes to macroscopic modeling
approaches, Vadose Zone J., 12, 1–16, <ext-link xlink:href="http://dx.doi.org/10.2136/vzj2013.02.0042" ext-link-type="DOI">10.2136/vzj2013.02.0042</ext-link>,
2013.</mixed-citation></ref>
      <ref id="bib1.bibx46"><label>Johnson et al.(2012)Johnson, Slater, Ntarlagiannis, Day-Lewis, and
Elwaseif</label><mixed-citation>Johnson, T., Slater, L., Ntarlagiannis, D., Day-Lewis, F., and Elwaseif, M.:
Monitoring groundwater-surface water interaction using time-series and
time-frequency analysis of transient three-dimensional electrical resistivity
changes, Water Resour. Res., 48, W07506, <ext-link xlink:href="http://dx.doi.org/10.1029/2012WR011893" ext-link-type="DOI">10.1029/2012WR011893</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bibx47"><label>Kelter et al.(2015)Kelter, Huisman, Zimmermann, Kemna, and
Vereecken</label><mixed-citation>Kelter, M., Huisman, J., Zimmermann, E., Kemna, A., and Vereecken, H.:
Quantitative imaging of spectral electrical properties of variably saturated
soil columns, J. Appl. Geophys., 123, 333–344,
<ext-link xlink:href="http://dx.doi.org/10.1016/j.jappgeo.2015.09.001" ext-link-type="DOI">10.1016/j.jappgeo.2015.09.001</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bibx48"><label>Kemna(2000)</label><mixed-citation>Kemna, A.: Tomographic inversion of complex resistivity – theory and
application, PhD thesis, Ruhr-Universität Bochum,
<ext-link xlink:href="http://dx.doi.org/10.1111/1365-2478.12013" ext-link-type="DOI">10.1111/1365-2478.12013</ext-link>, 2000.</mixed-citation></ref>
      <ref id="bib1.bibx49"><label>Kemna et al.(2012)Kemna, Binley, Cassiani, Niederleithinger, Revil,
Slater, Williams, Flores Orozco, Haegel, Hoerdt, Kruschwitz, Leroux, Titov,
and Zimermann</label><mixed-citation>Kemna, A., Binley, A., Cassiani, G., Niederleithinger, E., Revil, A., Slater,
L., Williams, K., Flores Orozco, A., Haegel, F., Hoerdt, A., Kruschwitz, S.,
Leroux, V., Titov, K., and Zimermann, E.: An overview of the spectral induced
polarization method for near-surface applications, Near Surf. Geophys.,
10, 453–468, <ext-link xlink:href="http://dx.doi.org/10.3997/1873-0604.2012027" ext-link-type="DOI">10.3997/1873-0604.2012027</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bibx50"><label>Kemna et al.(2014)Kemna, Huisman, Zimmermann, Martin, Zhao, Treichel,
Flores Orozco, and Fechner</label><mixed-citation>Kemna, A., Huisman, J., Zimmermann, E., Martin, R., Zhao, Y., Treichel, A.,
Flores Orozco, A., and Fechner, T.: Broadband electrical impedance tomography
for subsurface characterization using improved corrections of electromagnetic
coupling and spectral regularization, in: Tomography of the Earth's Crust:
From Geophysical Sounding to Real-Time Monitoring, Springer, 1–20, <ext-link xlink:href="http://dx.doi.org/10.1007/978-3-319-04205-3_1" ext-link-type="DOI">10.1007/978-3-319-04205-3_1</ext-link>,
2014.</mixed-citation></ref>
      <ref id="bib1.bibx51"><label>Kinraide(1994)</label><mixed-citation>
Kinraide, T.: Use of a Gouy-Chapman-Stern model for membrane-surface
electrical potential to interpret some features of mineral rhizotoxicity,
Plant Physiol., 106, 1583–1592, 1994.</mixed-citation></ref>
      <ref id="bib1.bibx52"><label>Kinraide and Wang(2010)</label><mixed-citation>Kinraide, T. and Wang, P.: The surface charge density of plant cell membranes
(<inline-formula><mml:math id="M112" display="inline"><mml:mi mathvariant="italic">σ</mml:mi></mml:math></inline-formula>): an attempt to resolve conflicting values for intrinsic <inline-formula><mml:math id="M113" display="inline"><mml:mi mathvariant="italic">σ</mml:mi></mml:math></inline-formula>,
J. Exp. Bot., 61, 2507–2518, <ext-link xlink:href="http://dx.doi.org/10.1093/jxb/erq082" ext-link-type="DOI">10.1093/jxb/erq082</ext-link>,
2010.</mixed-citation></ref>
      <ref id="bib1.bibx53"><label>Kinraide et al.(1998)Kinraide, Yermiyahu, and
Rytwo</label><mixed-citation>Kinraide, T., Yermiyahu, U., and Rytwo, G.: Computation of surface electrical
potentials of plant cell membranes correspondence to published zeta
potentials from diverse plant sources, Plant Physiol., 118, 505–512,
<ext-link xlink:href="http://dx.doi.org/10.1104/pp.118.2.505" ext-link-type="DOI">10.1104/pp.118.2.505</ext-link>, 1998.</mixed-citation></ref>
      <ref id="bib1.bibx54"><label>Kormanek et al.(2015)</label><mixed-citation>Kormanek, M., Gła̧b, T., and Klimek-Kopyra, A.: Modification of the
tree root
electrical capacitance method under laboratory conditions, Tree Physiol.,
36, 121–127, <ext-link xlink:href="http://dx.doi.org/10.1093/treephys/tpv088" ext-link-type="DOI">10.1093/treephys/tpv088</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bibx55"><label>Kruschwitz et al.(2010)Kruschwitz, Binley, Lesmes, and
Elshenawy</label><mixed-citation>Kruschwitz, S., Binley, A., Lesmes, D., and Elshenawy, A.: Textural controls
on
low-frequency electrical spectra of porous media, Geophysics, 75,
WA113, <ext-link xlink:href="http://dx.doi.org/10.1190/1.3479835" ext-link-type="DOI">10.1190/1.3479835</ext-link>, 2010.</mixed-citation></ref>
      <ref id="bib1.bibx56"><label>Kyle et al.(1999)Kyle, Chan, and Minchinton</label><mixed-citation>Kyle, A., Chan, C., and Minchinton, A.: Characterization of three-dimensional
tissue cultures using electrical impedance spectroscopy, Biophys. J.,
76, 2640–2648, <ext-link xlink:href="http://dx.doi.org/10.1016/S0006-3495(99)77416-3" ext-link-type="DOI">10.1016/S0006-3495(99)77416-3</ext-link>, 1999.</mixed-citation></ref>
      <ref id="bib1.bibx57"><label>LaBrecque and Ward(1990)</label><mixed-citation>
LaBrecque, D. and Ward, S.: Two-dimensional cross-borehole resistivity model
fitting, Geotechnical and Environmental Geophysics, 1, 51–57, 1990.</mixed-citation></ref>
      <ref id="bib1.bibx58"><label>LaBrecque et al.(1996)LaBrecque, Ramirez, Daily, Binley, and
Schima</label><mixed-citation>LaBrecque, D., Ramirez, A., Daily, W., Binley, A., and Schima, S.: ERT
monitoring of environmental remediation processes, Meas. Sci. Technol., 7,
375–383, <ext-link xlink:href="http://dx.doi.org/10.1088/0957-0233/7/3/019" ext-link-type="DOI">10.1088/0957-0233/7/3/019</ext-link>, 1996.</mixed-citation></ref>
      <ref id="bib1.bibx59"><label>Leroy et al.(2008)Leroy, Revil, Kemna, Cosenza, and
Ghorbani</label><mixed-citation>Leroy, P., Revil, A., Kemna, A., Cosenza, P., and Ghorbani, A.: Complex
conductivity of water-saturated packs of glass beads, J. Colloid
Interf. Sci., 321, 103–117, <ext-link xlink:href="http://dx.doi.org/10.1016/j.jcis.2007.12.031" ext-link-type="DOI">10.1016/j.jcis.2007.12.031</ext-link>, 2008.</mixed-citation></ref>
      <ref id="bib1.bibx60"><label>Lesmes and Frye(2001)</label><mixed-citation>Lesmes, D. and Frye, K.: Influence of pore fluid chemistry on the complex
conductivity and induced polarization responses of Berea sandstone, J. Geophys. Res., 106, 4079–4090, <ext-link xlink:href="http://dx.doi.org/10.1029/2000JB900392" ext-link-type="DOI">10.1029/2000JB900392</ext-link>, 2001.</mixed-citation></ref>
      <ref id="bib1.bibx61"><label>Li et al.(2015)Li, Liu, Zheng, and Xu</label><mixed-citation>Li, Z., Liu, Y., Zheng, Y., and Xu, R.: Zeta potential at the root surfaces
of
rice characterized by streaming potential measurements, Plant Soil, 386,
237–250, <ext-link xlink:href="http://dx.doi.org/10.1007/s11104-014-2259-6" ext-link-type="DOI">10.1007/s11104-014-2259-6</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bibx62"><label>Loke et al.(2013)Loke, Chambers, Rucker, Kuras, and
Wilkinson</label><mixed-citation>Loke, M., Chambers, J., Rucker, D., Kuras, O., and Wilkinson, P.: Recent
developments in the direct-current geoelectrical imaging method, J.
Appl. Geophys., 95, 135–156, <ext-link xlink:href="http://dx.doi.org/10.1016/j.jappgeo.2013.02.017" ext-link-type="DOI">10.1016/j.jappgeo.2013.02.017</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx63"><label>Lyklema(2005)</label><mixed-citation>
Lyklema, J.: Fundamentals of interface and colloid science: solid-liquid
interfaces, Vol. 2, Academic press, London, UK, 2005.</mixed-citation></ref>
      <ref id="bib1.bibx64"><label>Lyklema et al.(1983)Lyklema, Dukhin, and
Shilov</label><mixed-citation>Lyklema, J., Dukhin, S., and Shilov, V.: The relaxation of the double layer
around colloidal particles and the low-frequency dielectric dispersion: Part
I. Theoretical considerations, Journal of Electroanalytical Chemistry and
Interfacial Electrochemistry, 143, 1–21,
<ext-link xlink:href="http://dx.doi.org/10.1016/S0022-0728(83)80251-4" ext-link-type="DOI">10.1016/S0022-0728(83)80251-4</ext-link>, 1983.</mixed-citation></ref>
      <ref id="bib1.bibx65"><label>Mancuso(2012)</label><mixed-citation>Mancuso, S.: Measuring roots: an updated approach, Springer Science &amp;
Business
Media,
<ext-link xlink:href="http://dx.doi.org/10.1007/978-3-642-22067-8" ext-link-type="DOI">10.1007/978-3-642-22067-8</ext-link>, Springer, Heidelberg, Dordrecht, London, New
York, 2012.</mixed-citation></ref>
      <ref id="bib1.bibx66"><label>Martin and Günther(2013)</label><mixed-citation>Martin, T. and Günther, T.: Complex resistivity tomography (CRT) for
fungus
detection on standing oak trees, Eur. J. Forest Res., 132,
765–776, <ext-link xlink:href="http://dx.doi.org/10.1007/s10342-013-0711-4" ext-link-type="DOI">10.1007/s10342-013-0711-4</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx67"><label>Martínez-Ballesta et al.(2011)</label><mixed-citation>Martínez-Ballesta, M., Rodríguez-Hernández, M.,
Alcaraz-López, C., Mota-Cadenas, C., Muries, B., and Carvajal, M.: Plant
hydraulic conductivity: The aquaporins contribution, Hydraulic
conductivity – issues, determination and applications, Rijeka: In Tech,
103–21, <ext-link xlink:href="http://dx.doi.org/10.5772/18580" ext-link-type="DOI">10.5772/18580</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bibx68"><label>Metzner et al.(2015)Metzner, Eggert, van Dusschoten, Pflugfelder,
Gerth, Schurr, Uhlmann, and Jahnke</label><mixed-citation>Metzner, R., Eggert, A., van Dusschoten, D., Pflugfelder, D., Gerth, S.,
Schurr, U., Uhlmann, N., and Jahnke, S.: Direct comparison of MRI and X-ray
CT technologies for 3D imaging of root systems in soil: potential and
challenges for root trait quantification, Plant Meth., 11, 17,
<ext-link xlink:href="http://dx.doi.org/10.1186/s13007-015-0060-z" ext-link-type="DOI">10.1186/s13007-015-0060-z</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bibx69"><label>Nordsiek and Weller(2008)</label><mixed-citation>Nordsiek, S. and Weller, A.: A new approach to fitting induced-polarization
spectra, Geophysics, 73, F235–F245, <ext-link xlink:href="http://dx.doi.org/10.1190/1.2987412" ext-link-type="DOI">10.1190/1.2987412</ext-link>, 2008.</mixed-citation></ref>
      <ref id="bib1.bibx70"><label>Ozier-Lafontaine and Bajazet(2005)</label><mixed-citation>Ozier-Lafontaine, H. and Bajazet, T.: Analysis of root growth by impedance
spectroscopy (EIS), Plant Soil, 277, 299–313,
<ext-link xlink:href="http://dx.doi.org/10.1007/s11104-005-7531-3" ext-link-type="DOI">10.1007/s11104-005-7531-3</ext-link>, 2005.</mixed-citation></ref>
      <ref id="bib1.bibx71"><label>Pelton et al.(1978)Pelton, Ward, Hallof, Sill, and
Nelson</label><mixed-citation>Pelton, W., Ward, S., Hallof, P., Sill, W., and Nelson, P.: Mineral
discrimination and removal of inductive coupling with multifrequency IP,
Geophysics, 43, 588–609, <ext-link xlink:href="http://dx.doi.org/10.1190/1.1440839" ext-link-type="DOI">10.1190/1.1440839</ext-link>, 1978.</mixed-citation></ref>
      <ref id="bib1.bibx72"><label>Pierret et al.(2003)</label><mixed-citation>Pierret, A., Doussan, C., Garrigues, E., and Mc Kirby, J.: Observing plant
roots in their environment: current imaging options and specific contribution
of two-dimensional approaches, Agronomie, 23, 471–479,
<ext-link xlink:href="http://dx.doi.org/10.1051/agro:2003019" ext-link-type="DOI">10.1051/agro:2003019</ext-link>, 2003.</mixed-citation></ref>
      <ref id="bib1.bibx73"><label>Prodan and Prodan(1999)</label><mixed-citation>
Prodan, C. and Prodan, E.: The dielectric behaviour of living cell
suspensions,
J. Phys. D, 32, 335–343, 1999.</mixed-citation></ref>
      <ref id="bib1.bibx74"><label>Prodan et al.(2008)Prodan, Prodan, and Miller Jr</label><mixed-citation>Prodan, E., Prodan, C., and Miller Jr., J.: The dielectric response of
spherical
live cells in suspension: an analytic solution, Biophys. J., 95,
4174–4182, <ext-link xlink:href="http://dx.doi.org/10.1529/biophysj.108.137042" ext-link-type="DOI">10.1529/biophysj.108.137042</ext-link>, 2008.</mixed-citation></ref>
      <ref id="bib1.bibx75"><label>Razilov and Dukhin(1995)</label><mixed-citation>
Razilov, I. and Dukhin, S.: Simultaneous influence of concentration
polarization of the diffuse layer and polarization of the Stern layer
according to the mechanism of bound counterions at arbitrary magnitudes of
the relaxation parameter, Colloid J. Russ. Acad., 57, 364–371, 1995.</mixed-citation></ref>
      <ref id="bib1.bibx76"><label>Repo et al.(2012)Repo, Cao, Silvennoinen, and
Ozier-Lafontaine</label><mixed-citation>Repo, T., Cao, Y., Silvennoinen, R., and Ozier-Lafontaine, H.: Electrical
impedance spectroscopy and roots, in: Measuring Roots, edited by: Mancuso, S., Springer,
25–49, <ext-link xlink:href="http://dx.doi.org/10.1007/978-3-642-22067-8_2" ext-link-type="DOI">10.1007/978-3-642-22067-8_2</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bibx77"><label>Repo et al.(2014)Repo, Korhonen, Laukkanen, Lehto, and
Silvennoinen</label><mixed-citation>Repo, T., Korhonen, A., Laukkanen, M., Lehto, T., and Silvennoinen, R.:
Detecting mycorrhizal colonisation in Scots pine roots using electrical
impedance spectra, Biosyst. Eng., 121, 139–149,
<ext-link xlink:href="http://dx.doi.org/10.1016/j.biosystemseng.2014.02.014" ext-link-type="DOI">10.1016/j.biosystemseng.2014.02.014</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bibx78"><label>Revil and Florsch(2010)</label><mixed-citation>Revil, A. and Florsch, N.: Determination of permeability from spectral
induced
polarization in granular media, Geophys. J. Int., 181,
1480–1498, <ext-link xlink:href="http://dx.doi.org/10.1111/j.1365-246X.2010.04573.x" ext-link-type="DOI">10.1111/j.1365-246X.2010.04573.x</ext-link>, 2010.</mixed-citation></ref>
      <ref id="bib1.bibx79"><label>Revil et al.(2012)Revil, Karaoulis, Johnson, and
Kemna</label><mixed-citation>Revil, A., Karaoulis, M., Johnson, T., and Kemna, A.: Review: Some
low-frequency electrical methods for subsurface characterization and
monitoring in hydrogeology, Hydrogeol. J., 20, 617–658,
<ext-link xlink:href="http://dx.doi.org/10.1007/s10040-011-0819-x" ext-link-type="DOI">10.1007/s10040-011-0819-x</ext-link>, 2012.</mixed-citation></ref>
      <ref id="bib1.bibx80"><label>Revil et al.(2014)Revil, Florsch, and Camerlynck</label><mixed-citation>Revil, A., Florsch, N., and Camerlynck, C.: Spectral induced polarization
porosimetry, Geophys. J. Int., 198, 1016–1033,
<ext-link xlink:href="http://dx.doi.org/10.1093/gji/ggu180" ext-link-type="DOI">10.1093/gji/ggu180</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bibx81"><label>Rossi et al.(2011)Rossi, Amato, Bitella, Bochicchio, Ferreira Gomes,
Lovelli, Martorella, and Favale</label><mixed-citation>Rossi, R., Amato, M., Bitella, G., Bochicchio, R., Ferreira Gomes, J.,
Lovelli,
S., Martorella, E., and Favale, P.: Electrical resistivity tomography as a
non-destructive method for mapping root biomass in an orchard, Eur.
J. Soil Sci., 62, 206–215,
<ext-link xlink:href="http://dx.doi.org/10.1111/j.1365-2389.2010.01329.x" ext-link-type="DOI">10.1111/j.1365-2389.2010.01329.x</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bibx82"><label>Schraut et al.(2005)Schraut, Heilmeier, and
Hartung</label><mixed-citation>Schraut, D., Heilmeier, H., and Hartung, W.: Radial transport of water and
abscisic acid (ABA) in roots of Zea mays under conditions of nutrient
deficiency, J. Exp. Bot., 56, 879–886,
<ext-link xlink:href="http://dx.doi.org/10.1093/jxb/eri080" ext-link-type="DOI">10.1093/jxb/eri080</ext-link>, 2005.</mixed-citation></ref>
      <ref id="bib1.bibx83"><label>Schwan(1957)</label><mixed-citation>Schwan, H.: Electrical Properties of Tissue and Cell Suspensions, Vol. 5 of
Advances in Biological and Medical Physics, Elsevier, 147–209,
<ext-link xlink:href="http://dx.doi.org/10.1016/B978-1-4832-3111-2.50008-0" ext-link-type="DOI">10.1016/B978-1-4832-3111-2.50008-0</ext-link>,
1957.</mixed-citation></ref>
      <ref id="bib1.bibx84"><label>Schwarz(1962)</label><mixed-citation>Schwarz, G.: A theory of the low-frequency dielectric dispersion of colloidal
particles in electrolyte solution, J. Phys. Chem., 66,
2636–2642, <ext-link xlink:href="http://dx.doi.org/10.1021/j100818a067" ext-link-type="DOI">10.1021/j100818a067</ext-link>, 1962.</mixed-citation></ref>
      <ref id="bib1.bibx85"><label>Singha et al.(2014)Singha, Day-Lewis, Johnson, and
Slater</label><mixed-citation>Singha, K., Day-Lewis, F., Johnson, T., and Slater, L.: Advances in
interpretation of subsurface processes with time-lapse electrical imaging,
Hydrol. Process., 29, 1549–1576,  <ext-link xlink:href="http://dx.doi.org/10.1002/hyp.10280" ext-link-type="DOI">10.1002/hyp.10280</ext-link>, 2014.</mixed-citation></ref>
      <ref id="bib1.bibx86"><label>Srayeddin and Doussan(2009)</label><mixed-citation>Srayeddin, I. and Doussan, C.: Estimation of the spatial variability of root
water uptake of maize and sorghum at the field scale by electrical
resistivity tomography, Plant Soil, 319, 185–207,
<ext-link xlink:href="http://dx.doi.org/10.1007/s11104-008-9860-5" ext-link-type="DOI">10.1007/s11104-008-9860-5</ext-link>, 2009.</mixed-citation></ref>
      <ref id="bib1.bibx87"><label>Tarasov and Titov(2013)</label><mixed-citation>Tarasov, A. and Titov, K.: On the use of the Cole-Cole equations in spectral
induced polarization, Geophys. J. Int., 195, 352–356,
<ext-link xlink:href="http://dx.doi.org/10.1093/gji/ggt251" ext-link-type="DOI">10.1093/gji/ggt251</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx88"><label>Tinker and Nye(2000)</label><mixed-citation>Tinker, P. and Nye, P.: Solute movement in the rhizosphere, Oxford University
Press, Inc., New York, <ext-link xlink:href="http://dx.doi.org/10.1046/j.1365-2389.2001.00418-2.x" ext-link-type="DOI">10.1046/j.1365-2389.2001.00418-2.x</ext-link>,
2000.</mixed-citation></ref>
      <ref id="bib1.bibx89"><label>Titov et al.(2002)Titov, Komarov, Tarasov, and
Levitski</label><mixed-citation>Titov, K., Komarov, V., Tarasov, V., and Levitski, A.: Theoretical and
experimental study of time domain-induced polarization in water-saturated
sands, J. Appl. Geophys., 50, 417–433,
<ext-link xlink:href="http://dx.doi.org/10.1016/S0926-9851(02)00168-4" ext-link-type="DOI">10.1016/S0926-9851(02)00168-4</ext-link>, 2002.</mixed-citation></ref>
      <ref id="bib1.bibx90"><label>Uhlmann and Hakim(1971)</label><mixed-citation>Uhlmann, D. and Hakim, R.: Derivation of distribution functions from
relaxation
data, J. Phys. Chem. Solids, 32, 2652–2655,
<ext-link xlink:href="http://dx.doi.org/10.1016/S0022-3697(71)80114-2" ext-link-type="DOI">10.1016/S0022-3697(71)80114-2</ext-link>, 1971.</mixed-citation></ref>
      <ref id="bib1.bibx91"><label>Urban et al.(2011)Urban, Bequet, and Mainiero</label><mixed-citation>Urban, J., Bequet, R., and Mainiero, R.: Assessing the applicability of the
earth impedance method for in situ studies of tree root systems, J.
Exp. Bot., 62, 1857–1869, <ext-link xlink:href="http://dx.doi.org/10.1093/jxb/erq370" ext-link-type="DOI">10.1093/jxb/erq370</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bibx92"><label>Walker(1965)</label><mixed-citation>Walker, J.: Electrical AC resistance and capacitance of Zea mays L, Plant
Soil, 23, 270–274, <ext-link xlink:href="http://dx.doi.org/10.1007/BF01358354" ext-link-type="DOI">10.1007/BF01358354</ext-link>, 1965.</mixed-citation></ref>
      <ref id="bib1.bibx93"><label>Wang et al.(2009)Wang, Zhou, Li, and Li</label><mixed-citation>Wang, P., Zhou, D.-M., Li, L.-Z., and Li, D.-D.: What role does cell membrane
surface potential play in ion-plant interactions, Plant Signaling &amp;
Behavior, 4, 42–43, <ext-link xlink:href="http://dx.doi.org/10.4161/psb.4.1.7270" ext-link-type="DOI">10.4161/psb.4.1.7270</ext-link>, 2009.</mixed-citation></ref>
      <ref id="bib1.bibx94"><label>Wang et al.(2011)Wang, Kinraide, Zhou, K., and
Peijnenburg</label><mixed-citation>Wang, P., Kinraide, T., Zhou, D., K., P. M., and Peijnenburg, W.: Plasma
membrane surface potential: dual effects upon ion uptake and toxicity, Plant
Physiol., 155, 808–820, <ext-link xlink:href="http://dx.doi.org/10.1104/pp.110.165985" ext-link-type="DOI">10.1104/pp.110.165985</ext-link>, 2011.</mixed-citation></ref>
      <ref id="bib1.bibx95"><label>Wang et al.(2013)Wang, Kinraide, Wang, Zhou, and
Hao</label><mixed-citation>Wang, Y.-M., Kinraide, T., Wang, P., Zhou, D., and Hao, X.: Modeling
rhizotoxicity and uptake of Zn and Co singly and in binary mixture in wheat
in terms of the cell membrane surface electrical potential, Environ.
Sci. Technol., 47, 2831–2838, <ext-link xlink:href="http://dx.doi.org/10.1021/es3022107" ext-link-type="DOI">10.1021/es3022107</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx96"><label>Weigand and Kemna(2016)</label><mixed-citation>Weigand, M. and Kemna, A.: Debye decomposition of time-lapse spectral induced
polarisation data, Comput. Geosci., 86, 34–45,
<ext-link xlink:href="http://dx.doi.org/10.1016/j.cageo.2015.09.021" ext-link-type="DOI">10.1016/j.cageo.2015.09.021</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bibx97"><label>Weigand and Kemna(2017)</label><mixed-citation>Weigand, M. and Kemna, A.: Data and results for manuscript “Multi-frequency
electrical impedance tomography as a non-invasive tool to characterize and
monitor crop root systems”, Data set, Zenodo, <ext-link xlink:href="http://dx.doi.org/10.5281/zenodo.260087" ext-link-type="DOI">10.5281/zenodo.260087</ext-link>,
2017.</mixed-citation></ref>
      <ref id="bib1.bibx98"><label>Whalley et al.(2017)Whalley, Binley, Watts, Shanahan, Dodd, Ober,
Ashton, Webster, White, and Hawkesford</label><mixed-citation>Whalley, W., Binley, A., Watts, C., Shanahan, P., Dodd, I., Ober, E., Ashton,
R., Webster, C., White, R., and Hawkesford, M. J.: Methods to estimate
changes in soil water for phenotyping root activity in the field, Plant
Soil,  1–16, <ext-link xlink:href="http://dx.doi.org/10.1007/s11104-016-3161-1" ext-link-type="DOI">10.1007/s11104-016-3161-1</ext-link>, 2017.
</mixed-citation></ref><?xmltex \hack{\newpage}?>
      <ref id="bib1.bibx99"><label>Willatt et al.(1978)Willatt, Struss, and Taylor</label><mixed-citation>Willatt, S., Struss, R., and Taylor, H.: In situ root studies using neutron
radiography, Agron. J., 70, 581–586,
<ext-link xlink:href="http://dx.doi.org/10.2134/agronj1978.00021962007000040016x" ext-link-type="DOI">10.2134/agronj1978.00021962007000040016x</ext-link>, 1978.</mixed-citation></ref>
      <ref id="bib1.bibx100"><label>Zanetti et al.(2011)Zanetti, Weller, Vennetier, and
Mériaux</label><mixed-citation>Zanetti, C., Weller, A., Vennetier, M., and Mériaux, P.: Detection of
buried tree root samples by using geoelectrical measurements: a laboratory
experiment, Plant Soil, 339, 273–283, <ext-link xlink:href="http://dx.doi.org/10.1007/s11104-010-0574-0" ext-link-type="DOI">10.1007/s11104-010-0574-0</ext-link>,
2011.</mixed-citation></ref>
      <ref id="bib1.bibx101"><label>Zhao et al.(2013)Zhao, Zimmermann, Huisman, Treichel, Wolters, van
Waasen, and Kemna</label><mixed-citation>Zhao, Y., Zimmermann, E., Huisman, J., Treichel, A., Wolters, B., van Waasen,
S., and Kemna, A.: Broadband EIT borehole measurements with high phase
accuracy using numerical corrections of electromagnetic coupling effects,
Meas. Sci. Technol., 24, 085005,
<ext-link xlink:href="http://dx.doi.org/10.1088/0957-0233/24/8/085005" ext-link-type="DOI">10.1088/0957-0233/24/8/085005</ext-link>, 2013.</mixed-citation></ref>
      <ref id="bib1.bibx102"><label>Zimmermann et al.(2008)Zimmermann, Kemna, Berwix, Glaas, and
Vereecken</label><mixed-citation>Zimmermann, E., Kemna, A., Berwix, J., Glaas, W., and Vereecken, H.: EIT
measurement system with high phase accuracy for the imaging of spectral
induced polarization properties of soils and sediments, Meas. Sci.
Technol., 19, 094010, <ext-link xlink:href="http://dx.doi.org/10.1088/0957-0233/19/9/094010" ext-link-type="DOI">10.1088/0957-0233/19/9/094010</ext-link>, 2008.</mixed-citation></ref>

  </ref-list><app-group content-type="float"><app><title/>

    </app></app-group></back>
    <!--<article-title-html>Multi-frequency electrical impedance tomography as a non-invasive tool to characterize and monitor crop root systems</article-title-html>
<abstract-html><p class="p">A better understanding of root–soil interactions and associated processes is
essential in achieving progress in crop breeding and management, prompting
the need for high-resolution and non-destructive characterization methods. To
date, such methods are still lacking or restricted by technical constraints,
in particular the charactization and
monitoring of root growth and function in the field. A promising technique in
this respect is electrical impedance tomography (EIT), which utilizes
low-frequency ( &lt;  1 kHz)-
electrical conduction- and
polarization properties in an imaging framework. It is well established that
cells and cell clusters exhibit an electrical polarization response in
alternating electric-current fields
due to electrical double layers which form at cell membranes. This double
layer is directly related to the electrical surface properties of the
membrane, which in turn are influenced by nutrient dynamics (fluxes and
concentrations on both sides of the membranes). Therefore, it can be assumed
that the electrical polarization properties of roots are inherently related
to ion uptake and translocation processes in the root systems. We hereby
propose broadband (mHz to hundreds of Hz) multi-frequency EIT as a
non-invasive methodological approach for the monitoring and physiological,
i.e., functional, characterization of crop root systems. The approach
combines the spatial-resolution
capability of an imaging method with the diagnostic potential of
electrical-impedance spectroscopy.
The capability of multi-frequency EIT to characterize and monitor crop root
systems was investigated in a rhizotron laboratory experiment, in which the
root system of oilseed plants was monitored in a water–filled rhizotron,
that is, in a nutrient-deprived environment. We found a low-frequency
polarization response of the root system, which enabled the successful
delineation of its spatial extension. The magnitude of the overall
polarization response decreased along with the physiological decay of the
root system due to the stress situation. Spectral polarization parameters, as
derived from a pixel-based Debye decomposition analysis of the
multi-frequency imaging results, reveal systematic changes in the spatial and
spectral electrical response of the root system. In particular, quantified
mean relaxation times (of the order of 10 ms) indicate changes in the length scales on which the polarization
processes took place in the root system, as a response to the prolonged
induced stress situation. Our results demonstrate that broadband EIT is a
capable, non-invasive method to image root system extension as well as to
monitor changes associated with the root physiological processes. Given its
applicability on both laboratory and field scales, our results suggest an
enormous potential of the method for the structural and functional imaging of root systems for various
applications. This particularly holds for the field scale, where
corresponding methods are highly desired but to date are lacking.</p></abstract-html>
<ref-html id="bib1.bib1"><label>al Hagrey(2007)</label><mixed-citation>
al Hagrey, S.: Geophysical imaging of root-zone, trunk, and moisture
heterogeneity, J. Exp. Bot., 58, 839–854,
<a href="http://dx.doi.org/10.1093/jxb/erl237" target="_blank">doi:10.1093/jxb/erl237</a>, 2007.
</mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>al Hagrey and Petersen(2011)</label><mixed-citation>
al Hagrey, S. and Petersen, T.: Numerical and experimental mapping of small
root zones using optimized surface and borehole resistivity tomography,
Geophysics, 76, G25–G35, <a href="http://dx.doi.org/10.1190/1.3545067" target="_blank">doi:10.1190/1.3545067</a>, 2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>Alumbaugh and Newman(2000)</label><mixed-citation>
Alumbaugh, D. and Newman, G.: Image appraisal for 2-D and 3-D
electromagnetic inversion, Geophysics, 65, 1455–1467,
<a href="http://dx.doi.org/10.1190/1.1444834" target="_blank">doi:10.1190/1.1444834</a>, 2000.
</mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>Amato et al.(2008)Amato, Basso, Celano, Bitella, Morelli, and
Rossi</label><mixed-citation>
Amato, M., Basso, B., Celano, G., Bitella, G., Morelli, G., and Rossi, R.: In
situ detection of tree root distribution and biomass by multi-electrode
resistivity imaging, Tree Physiol., 28, 1441–1448,
<a href="http://dx.doi.org/10.1093/treephys/28.10.1441" target="_blank">doi:10.1093/treephys/28.10.1441</a>, 2008.
</mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>Amato et al.(2009)Amato, Bitella, Rossi, Gómez, Lovelli, and
Gomes</label><mixed-citation>
Amato, M., Bitella, G., Rossi, R., Gómez, J. A., Lovelli, S., and Gomes,
J. J. F.: Multi-electrode 3D resistivity imaging of alfalfa root zone, Eur.
J. Agron., 31, 213–222, <a href="http://dx.doi.org/10.1016/j.eja.2009.08.005" target="_blank">doi:10.1016/j.eja.2009.08.005</a>,
2009.
</mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>Anderson and Hopmans(2013)</label><mixed-citation>
Anderson, S. and Hopmans, J. (Eds.): Soil-water-root Processes: Advances in
Tomography and Imaging, Soil Sci. Soc. Am.,
<a href="http://dx.doi.org/10.2136/sssaspecpub61" target="_blank">doi:10.2136/sssaspecpub61</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>Anderson and Higinbotham(1976)</label><mixed-citation>
Anderson, W. and Higinbotham, N.: Electrical resistances of corn root
segments,
Plant Physiol., 57, 137–141, <a href="http://dx.doi.org/10.1104/pp.57.2.137" target="_blank">doi:10.1104/pp.57.2.137</a>, 1976.
</mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>Aubrecht et al.(2006)Aubrecht, Staněk, and
Koller</label><mixed-citation>
Aubrecht, L., Staněk, Z., and Koller, J.: Electrical measurement of the
absorption surfaces of tree roots by the earth impedance method: 1. Theory,
Tree Physiol., 26, 1105–1112, <a href="http://dx.doi.org/10.1093/treephys/26.9.1105" target="_blank">doi:10.1093/treephys/26.9.1105</a>, 2006.
</mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>Aulen and Shipley(2012)</label><mixed-citation>
Aulen, M. and Shipley, B.: Non-destructive estimation of root mass using
electrical capacitance on ten herbaceous species, Plant Soil, 355,
41–49, <a href="http://dx.doi.org/10.1007/s11104-011-1077-3" target="_blank">doi:10.1007/s11104-011-1077-3</a>, 2012.
</mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>Barsoukov and Macdonald(2005)</label><mixed-citation>
Barsoukov, E. and Macdonald, J. (Eds.): Impedance spectroscopy: theory,
experiment, and applications, Wiley-Interscience, <a href="http://dx.doi.org/10.1002/0471716243" target="_blank">doi:10.1002/0471716243</a>,
2005.
</mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>Bayford(2006)</label><mixed-citation>
Bayford, R.: Bioimpedance tomography (electrical impedance tomography),
Annual
Reviews Biomedical Engineering, 8, 63–91,
<a href="http://dx.doi.org/10.1146/annurev.bioeng.8.061505.095716" target="_blank">doi:10.1146/annurev.bioeng.8.061505.095716</a>, 2006.
</mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>Beff et al.(2013)Beff, Günther, Vandoorne, Couvreur, and
Javaux</label><mixed-citation>
Beff, L., Günther, T., Vandoorne, B., Couvreur, V., and Javaux, M.:
Three-dimensional monitoring of soil water content in a maize field using
Electrical Resistivity Tomography, Hydrol. Earth Syst. Sci., 17, 595–609,
<a href="http://dx.doi.org/10.5194/hess-17-595-2013" target="_blank">doi:10.5194/hess-17-595-2013</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>Benlloch-González et al.(2010)Benlloch-González, Fournier,
and Benlloch</label><mixed-citation>
Benlloch-González, M., Fournier, J., and Benlloch, M.: K<sup>+</sup> deprivation
induces xylem water and K<sup>+</sup> transport in sunflower: evidence for a
co-ordinated control, J. Exp. Bot., 61, 157–164,
<a href="http://dx.doi.org/10.1093/jxb/erp288" target="_blank">doi:10.1093/jxb/erp288</a>, 2010.
</mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>Binley and Kemna(2005)</label><mixed-citation>
Binley, A. and Kemna, A.: DC resistivity and induced polarization methods,
in:
Hydrogeophysics, edited by: Rubin, Y. and Hubbard, S. S.,
Springer, the
Netherlands, 129–156, <a href="http://dx.doi.org/10.1007/1-4020-3102-5_5" target="_blank">doi:10.1007/1-4020-3102-5_5</a>, 2005.
</mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>Binley et al.(2005)Binley, Slater, Fukes, and
Cassiani</label><mixed-citation>
Binley, A., Slater, L., Fukes, M., and Cassiani, G.: Relationship between
spectral induced polarization and hydraulic properties of saturated and
unsaturated sandstone, Water Resour. Res., 41, 1–13,
<a href="http://dx.doi.org/10.1029/2005WR004202" target="_blank">doi:10.1029/2005WR004202</a>, 2005.
</mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>Boaga et al.(2013)Boaga, Rossi, and Cassiani</label><mixed-citation>
Boaga, J., Rossi, M., and Cassiani, G.: Monitoring soil-plant interactions in
an apple orchard using 3D electrical resistivity tomography, Procedia
Environmental Sciences, 19, 394–402, <a href="http://dx.doi.org/10.1016/j.proenv.2013.06.045" target="_blank">doi:10.1016/j.proenv.2013.06.045</a>,
2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>Bücker and Hördt(2013)</label><mixed-citation>
Bücker, M. and Hördt, A.: Analytical modelling of membrane
polarization
with explicit parametrization of pore radii and the electrical double layer,
Geophys. J. Int., 194, 804–813, <a href="http://dx.doi.org/10.1093/gji/ggt136" target="_blank">doi:10.1093/gji/ggt136</a>,
2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>Cao et al.(2010)Cao, Repo, Silvennoinen, Lehto, and
Pelkonen</label><mixed-citation>
Cao, Y., Repo, T., Silvennoinen, R., Lehto, T., and Pelkonen, P.: An
appraisal
of the electrical resistance method for assessing root surface area, J. Exp. Bot., 61, 2491–2497, <a href="http://dx.doi.org/10.1093/jxb/erq078" target="_blank">doi:10.1093/jxb/erq078</a>, 2010.
</mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>Cao et al.(2011)Cao, Repo, Silvennoinen, Lehto, and
Pelkonen</label><mixed-citation>
Cao, Y., Repo, T., Silvennoinen, R., Lehto, T., and Pelkonen, P.: Analysis of
the willow root system by electrical impedance spectroscopy, J.
Exp. Bot., 62, 351–358, <a href="http://dx.doi.org/10.1093/jxb/erq276" target="_blank">doi:10.1093/jxb/erq276</a>, 2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>Čermák et al.(2006)</label><mixed-citation>
Čermák, J., Ulrich, R., Staněk, Z., Koller, J., and Aubrecht,
L.: Electrical measurement of tree root absorbing surfaces by the earth
impedance method: 2. Verification based on allometric relationships and
root severing experiments, Tree Physiol., 26, 1113–1121,
<a href="http://dx.doi.org/10.1093/treephys/26.9.1113" target="_blank">doi:10.1093/treephys/26.9.1113</a>, 2006.
</mixed-citation></ref-html>
<ref-html id="bib1.bib21"><label>Chloupek(1972)</label><mixed-citation>
Chloupek, O.: The relationship between electric capacitance and some other
parameters of plant roots, Biol. Plantarum, 14, 227–230,
<a href="http://dx.doi.org/10.1007/BF02921255" target="_blank">doi:10.1007/BF02921255</a>, 1972.
</mixed-citation></ref-html>
<ref-html id="bib1.bib22"><label>Chloupek(1976)</label><mixed-citation>
Chloupek, O.: Die Bewertung des Wurzelsystems von Senfpflanzen auf Grund der
dielektrischen Eigenschaften und mit Rücksicht auf den Endertrag,
Biol. Plantarum, 18, 44–49, <a href="http://dx.doi.org/10.1007/BF02922333" target="_blank">doi:10.1007/BF02922333</a>, 1976.
</mixed-citation></ref-html>
<ref-html id="bib1.bib23"><label>Claassen and Barber(1974)</label><mixed-citation>
Claassen, N. and Barber, S.: A method for characterizing the relation between
nutrient concentration and flux into roots of intact plants, Plant
Physiol., 54, 564–568, <a href="http://dx.doi.org/10.1104/pp.54.4.564" target="_blank">doi:10.1104/pp.54.4.564</a>, 1974.
</mixed-citation></ref-html>
<ref-html id="bib1.bib24"><label>Clarkson et al.(2000)Clarkson, Carvajal, Henzler, Waterhouse, Smyth,
Cooke, and Steudle</label><mixed-citation>
Clarkson, D., Carvajal, M., Henzler, T., Waterhouse, R., Smyth, A., Cooke,
D.,
and Steudle, E.: Root hydraulic conductance: diurnal aquaporin expression and
the effects of nutrient stress, J. Exp. Bot., 51, 61–70,
<a href="http://dx.doi.org/10.1093/jexbot/51.342.61" target="_blank">doi:10.1093/jexbot/51.342.61</a>, 2000.
</mixed-citation></ref-html>
<ref-html id="bib1.bib25"><label>Cseresnyés et al.(2013)Cseresnyés, Tünde, Végh,
Anton, and Rajkai</label><mixed-citation>
Cseresnyés, I., Tünde, T., Végh, K., Anton, A., and Rajkai, K.:
Electrical impedance and capacitance method: A new approach for detection of
functional aspects of arbuscular mycorrhizal colonization in maize, European
J. Soil Biol., 54, 25–31, <a href="http://dx.doi.org/10.1016/j.ejsobi.2012.11.001" target="_blank">doi:10.1016/j.ejsobi.2012.11.001</a>,
2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib26"><label>Daily et al.(2005)Daily, Ramirez, Binley, and
LaBrecque</label><mixed-citation>
Daily, W., Ramirez, A., Binley, A., and LaBrecque, D.: Electrical resistance
tomography – theory and practice, in: Near-Surface Geophysics, edited by:
Butler, D., Vol. 154, chap. 17, Society of Exploration
Geophysicists, Tulsa, OK, 525–550, <a href="http://dx.doi.org/10.1190/1.9781560801719.ch17" target="_blank">doi:10.1190/1.9781560801719.ch17</a>, 2005.
</mixed-citation></ref-html>
<ref-html id="bib1.bib27"><label>Dalton(1995)</label><mixed-citation>
Dalton, F.: In-situ root extent measurements by electrical capacitance
methods,
Plant Soil, 173, 157–165, <a href="http://dx.doi.org/10.1007/BF00155527" target="_blank">doi:10.1007/BF00155527</a>, 1995.
</mixed-citation></ref-html>
<ref-html id="bib1.bib28"><label>Delhon et al.(1995)Delhon, Gojon, Tillard, and
Passama</label><mixed-citation>
Delhon, P., Gojon, A., Tillard, P., and Passama, L.: Diurnal regulation of
NO<sub>3</sub><sup>−</sup> uptake in soybean plants I. Changes in NO<sub>3</sub><sup>−</sup> influx, efflux, and N
utilization in the plant during the day/night cycle, J. Exp.
Bot., 46, 1585–1594, <a href="http://dx.doi.org/10.1093/jxb/46.10.1585" target="_blank">doi:10.1093/jxb/46.10.1585</a>, 1995.
</mixed-citation></ref-html>
<ref-html id="bib1.bib29"><label>Dietrich et al.(2013)Dietrich, Bengough, Jones, and
White</label><mixed-citation>
Dietrich, R., Bengough, A., Jones, H., and White, P.: Can root electrical
capacitance be used to predict root mass in soil?, Ann. Bot., 112,
457–464, <a href="http://dx.doi.org/10.1093/aob/mct044" target="_blank">doi:10.1093/aob/mct044</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib30"><label>Dukhin and Shilov(1974)</label><mixed-citation>
Dukhin, S. S. and Shilov, V. N.: Dielectric Phenomena and the Double Layer in
Disperse Systems and Polyelectrolytes, Wiley, New York, 1974.
</mixed-citation></ref-html>
<ref-html id="bib1.bib31"><label>Dunbabin et al.(2013)Dunbabin, Postma, Schnepf, Pagès, Javaux,
Wu, Leitner, Chen, Rengel, and Diggle</label><mixed-citation>
Dunbabin, V., Postma, J., Schnepf, A., Pagès, L., Javaux, M., Wu, L.,
Leitner, D., Chen, Y., Rengel, Z., and Diggle, A.: Modelling root-soil
interactions using three-dimensional models of root growth, architecture and
function, Plant Soil, 372, 93–124, <a href="http://dx.doi.org/10.1007/s11104-013-1769-y" target="_blank">doi:10.1007/s11104-013-1769-y</a>,
2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib32"><label>Dvořák et al.(1981)</label><mixed-citation>
Dvořák, M., Černohorská, J., and Janáček, K.:
Characteristics of current passage through plant tissue, Biol. Plantarum,
23, 306–310, <a href="http://dx.doi.org/10.1007/BF02895374" target="_blank">doi:10.1007/BF02895374</a>, 1981.
</mixed-citation></ref-html>
<ref-html id="bib1.bib33"><label>Ellis et al.(2013)Ellis, Murray, and Kavalieris</label><mixed-citation>
Ellis, T., Murray, W., and Kavalieris, L.: Electrical capacitance of bean
(Vicia faba) root systems was related to tissue density – a test for the
Dalton Model, Plant Soil, 366, 575–584,
<a href="http://dx.doi.org/10.1007/s11104-012-1424-z" target="_blank">doi:10.1007/s11104-012-1424-z</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib34"><label>Fixman(1980)</label><mixed-citation>
Fixman, M.: Charged Macromolecules in External Fields. 2. Preliminary
Remarks
on the Cylinder, Macromolecules, 13, 711–716, <a href="http://dx.doi.org/10.1021/ma60075a043" target="_blank">doi:10.1021/ma60075a043</a>,
1980.
</mixed-citation></ref-html>
<ref-html id="bib1.bib35"><label>Flores Orozco et al.(2012a)Flores Orozco, Kemna,
Oberdörster, Zschornack, Leven, Dietrich, and
Weiss</label><mixed-citation>
Flores Orozco, A., Kemna, A., Oberdörster, C., Zschornack, L., Leven, C.,
Dietrich, P., and Weiss, H.: Delineation of subsurface hydrocarbon
contamination at a former hydrogenation plant using spectral induced
polarization imaging, J. Contam. Hydrol., 136, 131–144,
<a href="http://dx.doi.org/10.1016/j.jconhyd.2012.06.001" target="_blank">doi:10.1016/j.jconhyd.2012.06.001</a>, 2012a.
</mixed-citation></ref-html>
<ref-html id="bib1.bib36"><label>Flores Orozco et al.(2012b)Flores Orozco, Kemna, and
Zimmermann</label><mixed-citation>
Flores Orozco, A., Kemna, A., and Zimmermann, E.: Data error quantification
in
spectral induced polarization imaging, Geophysics, 77, E227–E237,
<a href="http://dx.doi.org/10.1190/geo2010-0194.1" target="_blank">doi:10.1190/geo2010-0194.1</a>, 2012b.
</mixed-citation></ref-html>
<ref-html id="bib1.bib37"><label>Flores Orozco et al.(2013)Flores Orozco, Williams, and
Kemna</label><mixed-citation>
Flores Orozco, A., Williams, K., and Kemna, A.: Time-lapse spectral induced
polarization imaging of stimulated uranium bioremediation, Near Surf.
Geophys., 11, 531–544, <a href="http://dx.doi.org/10.3997/1873-0604.2013020" target="_blank">doi:10.3997/1873-0604.2013020</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib38"><label>Friedel(2003)</label><mixed-citation>
Friedel, S.: Resolution, stability and efficiency of resistivity tomography
estimated from a generalized inverse approach, Geophys. J.
Int., 153, 305–316, <a href="http://dx.doi.org/10.1046/j.1365-246X.2003.01890.x" target="_blank">doi:10.1046/j.1365-246X.2003.01890.x</a>, 2003.
</mixed-citation></ref-html>
<ref-html id="bib1.bib39"><label>Gregory et al.(2003)Gregory, Hutchison, Read, Jenneson, Gilboy, and
Morton</label><mixed-citation>
Gregory, P., Hutchison, D., Read, D., Jenneson, P., Gilboy, W., and Morton,
E.:
Non-invasive imaging of roots with high resolution X-ray micro-tomography,
in: Roots: The Dynamic Interface between Plants and the Earth,
Springer, 351–359, <a href="http://dx.doi.org/10.1023/A:1026179919689" target="_blank">doi:10.1023/A:1026179919689</a>, 2003.
</mixed-citation></ref-html>
<ref-html id="bib1.bib40"><label>Günther and Martin(2016)</label><mixed-citation>
Günther, T. and Martin, T.: Spectral two-dimensional inversion of
frequency-domain induced polarization data from a mining slag heap, J. Appl. Geophys., 135, 436–448,  <a href="http://dx.doi.org/10.1016/j.jappgeo.2016.01.008" target="_blank">doi:10.1016/j.jappgeo.2016.01.008</a>, 2016.
</mixed-citation></ref-html>
<ref-html id="bib1.bib41"><label>Heege(2013)</label><mixed-citation>
Heege, H.: Precision in Crop Farming: Site Specific Concepts and Sensing
Methods: Applications and Results, Springer Science &amp; Business Media, Springer, the
Netherlands,
2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib42"><label>Heřmanská et al.(2015)</label><mixed-citation>
Heřmanská, A., Středa, T., and Chloupek, O.: Improved wheat
grain yield by a new method of root selection, Agron. Sustain.
Dev., 35, 195–202, <a href="http://dx.doi.org/10.1007/s13593-014-0227-4" target="_blank">doi:10.1007/s13593-014-0227-4</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib43"><label>Hose et al.(2001)Hose, Clarkson, Steudle, Schreiber, and
Hartung</label><mixed-citation>
Hose, E., Clarkson, D., Steudle, E., Schreiber, L., and Hartung, W.: The
exodermis: a variable apoplastic barrier, J. Exp. Bot., 52,
2245–2264, <a href="http://dx.doi.org/10.1093/jexbot/52.365.2245" target="_blank">doi:10.1093/jexbot/52.365.2245</a>, 2001.
</mixed-citation></ref-html>
<ref-html id="bib1.bib44"><label>Huisman et al.(2015)Huisman, Zimmermann, Esser, Haegel, Treichel, and
Vereecken</label><mixed-citation>
Huisman, J., Zimmermann, E., Esser, O., Haegel, F., Treichel, A., and
Vereecken, H.: Evaluation of a novel correction procedure to remove
electrode impedance effects from broadband SIP measurements, J.
Appl. Geophys., 135, 466–473,  <a href="http://dx.doi.org/10.1016/j.jappgeo.2015.11.008" target="_blank">doi:10.1016/j.jappgeo.2015.11.008</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib45"><label>Javaux et al.(2013)Javaux, Couvreur, Vanderborght, and
Vereecken</label><mixed-citation>
Javaux, M., Couvreur, V., Vanderborght, J., and Vereecken, H.: Root water
uptake: From three-dimensional biophysical processes to macroscopic modeling
approaches, Vadose Zone J., 12, 1–16, <a href="http://dx.doi.org/10.2136/vzj2013.02.0042" target="_blank">doi:10.2136/vzj2013.02.0042</a>,
2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib46"><label>Johnson et al.(2012)Johnson, Slater, Ntarlagiannis, Day-Lewis, and
Elwaseif</label><mixed-citation>
Johnson, T., Slater, L., Ntarlagiannis, D., Day-Lewis, F., and Elwaseif, M.:
Monitoring groundwater-surface water interaction using time-series and
time-frequency analysis of transient three-dimensional electrical resistivity
changes, Water Resour. Res., 48, W07506, <a href="http://dx.doi.org/10.1029/2012WR011893" target="_blank">doi:10.1029/2012WR011893</a>, 2012.
</mixed-citation></ref-html>
<ref-html id="bib1.bib47"><label>Kelter et al.(2015)Kelter, Huisman, Zimmermann, Kemna, and
Vereecken</label><mixed-citation>
Kelter, M., Huisman, J., Zimmermann, E., Kemna, A., and Vereecken, H.:
Quantitative imaging of spectral electrical properties of variably saturated
soil columns, J. Appl. Geophys., 123, 333–344,
<a href="http://dx.doi.org/10.1016/j.jappgeo.2015.09.001" target="_blank">doi:10.1016/j.jappgeo.2015.09.001</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib48"><label>Kemna(2000)</label><mixed-citation>
Kemna, A.: Tomographic inversion of complex resistivity – theory and
application, PhD thesis, Ruhr-Universität Bochum,
<a href="http://dx.doi.org/10.1111/1365-2478.12013" target="_blank">doi:10.1111/1365-2478.12013</a>, 2000.
</mixed-citation></ref-html>
<ref-html id="bib1.bib49"><label>Kemna et al.(2012)Kemna, Binley, Cassiani, Niederleithinger, Revil,
Slater, Williams, Flores Orozco, Haegel, Hoerdt, Kruschwitz, Leroux, Titov,
and Zimermann</label><mixed-citation>
Kemna, A., Binley, A., Cassiani, G., Niederleithinger, E., Revil, A., Slater,
L., Williams, K., Flores Orozco, A., Haegel, F., Hoerdt, A., Kruschwitz, S.,
Leroux, V., Titov, K., and Zimermann, E.: An overview of the spectral induced
polarization method for near-surface applications, Near Surf. Geophys.,
10, 453–468, <a href="http://dx.doi.org/10.3997/1873-0604.2012027" target="_blank">doi:10.3997/1873-0604.2012027</a>, 2012.
</mixed-citation></ref-html>
<ref-html id="bib1.bib50"><label>Kemna et al.(2014)Kemna, Huisman, Zimmermann, Martin, Zhao, Treichel,
Flores Orozco, and Fechner</label><mixed-citation>
Kemna, A., Huisman, J., Zimmermann, E., Martin, R., Zhao, Y., Treichel, A.,
Flores Orozco, A., and Fechner, T.: Broadband electrical impedance tomography
for subsurface characterization using improved corrections of electromagnetic
coupling and spectral regularization, in: Tomography of the Earth's Crust:
From Geophysical Sounding to Real-Time Monitoring, Springer, 1–20, <a href="http://dx.doi.org/10.1007/978-3-319-04205-3_1" target="_blank">doi:10.1007/978-3-319-04205-3_1</a>,
2014.
</mixed-citation></ref-html>
<ref-html id="bib1.bib51"><label>Kinraide(1994)</label><mixed-citation>
Kinraide, T.: Use of a Gouy-Chapman-Stern model for membrane-surface
electrical potential to interpret some features of mineral rhizotoxicity,
Plant Physiol., 106, 1583–1592, 1994.
</mixed-citation></ref-html>
<ref-html id="bib1.bib52"><label>Kinraide and Wang(2010)</label><mixed-citation>
Kinraide, T. and Wang, P.: The surface charge density of plant cell membranes
(<i>σ</i>): an attempt to resolve conflicting values for intrinsic <i>σ</i>,
J. Exp. Bot., 61, 2507–2518, <a href="http://dx.doi.org/10.1093/jxb/erq082" target="_blank">doi:10.1093/jxb/erq082</a>,
2010.
</mixed-citation></ref-html>
<ref-html id="bib1.bib53"><label>Kinraide et al.(1998)Kinraide, Yermiyahu, and
Rytwo</label><mixed-citation>
Kinraide, T., Yermiyahu, U., and Rytwo, G.: Computation of surface electrical
potentials of plant cell membranes correspondence to published zeta
potentials from diverse plant sources, Plant Physiol., 118, 505–512,
<a href="http://dx.doi.org/10.1104/pp.118.2.505" target="_blank">doi:10.1104/pp.118.2.505</a>, 1998.
</mixed-citation></ref-html>
<ref-html id="bib1.bib54"><label>Kormanek et al.(2015)</label><mixed-citation>
Kormanek, M., Gła̧b, T., and Klimek-Kopyra, A.: Modification of the
tree root
electrical capacitance method under laboratory conditions, Tree Physiol.,
36, 121–127, <a href="http://dx.doi.org/10.1093/treephys/tpv088" target="_blank">doi:10.1093/treephys/tpv088</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib55"><label>Kruschwitz et al.(2010)Kruschwitz, Binley, Lesmes, and
Elshenawy</label><mixed-citation>
Kruschwitz, S., Binley, A., Lesmes, D., and Elshenawy, A.: Textural controls
on
low-frequency electrical spectra of porous media, Geophysics, 75,
WA113, <a href="http://dx.doi.org/10.1190/1.3479835" target="_blank">doi:10.1190/1.3479835</a>, 2010.
</mixed-citation></ref-html>
<ref-html id="bib1.bib56"><label>Kyle et al.(1999)Kyle, Chan, and Minchinton</label><mixed-citation>
Kyle, A., Chan, C., and Minchinton, A.: Characterization of three-dimensional
tissue cultures using electrical impedance spectroscopy, Biophys. J.,
76, 2640–2648, <a href="http://dx.doi.org/10.1016/S0006-3495(99)77416-3" target="_blank">doi:10.1016/S0006-3495(99)77416-3</a>, 1999.
</mixed-citation></ref-html>
<ref-html id="bib1.bib57"><label>LaBrecque and Ward(1990)</label><mixed-citation>
LaBrecque, D. and Ward, S.: Two-dimensional cross-borehole resistivity model
fitting, Geotechnical and Environmental Geophysics, 1, 51–57, 1990.
</mixed-citation></ref-html>
<ref-html id="bib1.bib58"><label>LaBrecque et al.(1996)LaBrecque, Ramirez, Daily, Binley, and
Schima</label><mixed-citation>
LaBrecque, D., Ramirez, A., Daily, W., Binley, A., and Schima, S.: ERT
monitoring of environmental remediation processes, Meas. Sci. Technol., 7,
375–383, <a href="http://dx.doi.org/10.1088/0957-0233/7/3/019" target="_blank">doi:10.1088/0957-0233/7/3/019</a>, 1996.
</mixed-citation></ref-html>
<ref-html id="bib1.bib59"><label>Leroy et al.(2008)Leroy, Revil, Kemna, Cosenza, and
Ghorbani</label><mixed-citation>
Leroy, P., Revil, A., Kemna, A., Cosenza, P., and Ghorbani, A.: Complex
conductivity of water-saturated packs of glass beads, J. Colloid
Interf. Sci., 321, 103–117, <a href="http://dx.doi.org/10.1016/j.jcis.2007.12.031" target="_blank">doi:10.1016/j.jcis.2007.12.031</a>, 2008.
</mixed-citation></ref-html>
<ref-html id="bib1.bib60"><label>Lesmes and Frye(2001)</label><mixed-citation>
Lesmes, D. and Frye, K.: Influence of pore fluid chemistry on the complex
conductivity and induced polarization responses of Berea sandstone, J. Geophys. Res., 106, 4079–4090, <a href="http://dx.doi.org/10.1029/2000JB900392" target="_blank">doi:10.1029/2000JB900392</a>, 2001.
</mixed-citation></ref-html>
<ref-html id="bib1.bib61"><label>Li et al.(2015)Li, Liu, Zheng, and Xu</label><mixed-citation>
Li, Z., Liu, Y., Zheng, Y., and Xu, R.: Zeta potential at the root surfaces
of
rice characterized by streaming potential measurements, Plant Soil, 386,
237–250, <a href="http://dx.doi.org/10.1007/s11104-014-2259-6" target="_blank">doi:10.1007/s11104-014-2259-6</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib62"><label>Loke et al.(2013)Loke, Chambers, Rucker, Kuras, and
Wilkinson</label><mixed-citation>
Loke, M., Chambers, J., Rucker, D., Kuras, O., and Wilkinson, P.: Recent
developments in the direct-current geoelectrical imaging method, J.
Appl. Geophys., 95, 135–156, <a href="http://dx.doi.org/10.1016/j.jappgeo.2013.02.017" target="_blank">doi:10.1016/j.jappgeo.2013.02.017</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib63"><label>Lyklema(2005)</label><mixed-citation>
Lyklema, J.: Fundamentals of interface and colloid science: solid-liquid
interfaces, Vol. 2, Academic press, London, UK, 2005.
</mixed-citation></ref-html>
<ref-html id="bib1.bib64"><label>Lyklema et al.(1983)Lyklema, Dukhin, and
Shilov</label><mixed-citation>
Lyklema, J., Dukhin, S., and Shilov, V.: The relaxation of the double layer
around colloidal particles and the low-frequency dielectric dispersion: Part
I. Theoretical considerations, Journal of Electroanalytical Chemistry and
Interfacial Electrochemistry, 143, 1–21,
<a href="http://dx.doi.org/10.1016/S0022-0728(83)80251-4" target="_blank">doi:10.1016/S0022-0728(83)80251-4</a>, 1983.
</mixed-citation></ref-html>
<ref-html id="bib1.bib65"><label>Mancuso(2012)</label><mixed-citation>
Mancuso, S.: Measuring roots: an updated approach, Springer Science &amp;
Business
Media,
<a href="http://dx.doi.org/10.1007/978-3-642-22067-8" target="_blank">doi:10.1007/978-3-642-22067-8</a>, Springer, Heidelberg, Dordrecht, London, New
York, 2012.
</mixed-citation></ref-html>
<ref-html id="bib1.bib66"><label>Martin and Günther(2013)</label><mixed-citation>
Martin, T. and Günther, T.: Complex resistivity tomography (CRT) for
fungus
detection on standing oak trees, Eur. J. Forest Res., 132,
765–776, <a href="http://dx.doi.org/10.1007/s10342-013-0711-4" target="_blank">doi:10.1007/s10342-013-0711-4</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib67"><label>Martínez-Ballesta et al.(2011)</label><mixed-citation>
Martínez-Ballesta, M., Rodríguez-Hernández, M.,
Alcaraz-López, C., Mota-Cadenas, C., Muries, B., and Carvajal, M.: Plant
hydraulic conductivity: The aquaporins contribution, Hydraulic
conductivity – issues, determination and applications, Rijeka: In Tech,
103–21, <a href="http://dx.doi.org/10.5772/18580" target="_blank">doi:10.5772/18580</a>, 2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib68"><label>Metzner et al.(2015)Metzner, Eggert, van Dusschoten, Pflugfelder,
Gerth, Schurr, Uhlmann, and Jahnke</label><mixed-citation>
Metzner, R., Eggert, A., van Dusschoten, D., Pflugfelder, D., Gerth, S.,
Schurr, U., Uhlmann, N., and Jahnke, S.: Direct comparison of MRI and X-ray
CT technologies for 3D imaging of root systems in soil: potential and
challenges for root trait quantification, Plant Meth., 11, 17,
<a href="http://dx.doi.org/10.1186/s13007-015-0060-z" target="_blank">doi:10.1186/s13007-015-0060-z</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib69"><label>Nordsiek and Weller(2008)</label><mixed-citation>
Nordsiek, S. and Weller, A.: A new approach to fitting induced-polarization
spectra, Geophysics, 73, F235–F245, <a href="http://dx.doi.org/10.1190/1.2987412" target="_blank">doi:10.1190/1.2987412</a>, 2008.
</mixed-citation></ref-html>
<ref-html id="bib1.bib70"><label>Ozier-Lafontaine and Bajazet(2005)</label><mixed-citation>
Ozier-Lafontaine, H. and Bajazet, T.: Analysis of root growth by impedance
spectroscopy (EIS), Plant Soil, 277, 299–313,
<a href="http://dx.doi.org/10.1007/s11104-005-7531-3" target="_blank">doi:10.1007/s11104-005-7531-3</a>, 2005.
</mixed-citation></ref-html>
<ref-html id="bib1.bib71"><label>Pelton et al.(1978)Pelton, Ward, Hallof, Sill, and
Nelson</label><mixed-citation>
Pelton, W., Ward, S., Hallof, P., Sill, W., and Nelson, P.: Mineral
discrimination and removal of inductive coupling with multifrequency IP,
Geophysics, 43, 588–609, <a href="http://dx.doi.org/10.1190/1.1440839" target="_blank">doi:10.1190/1.1440839</a>, 1978.
</mixed-citation></ref-html>
<ref-html id="bib1.bib72"><label>Pierret et al.(2003)</label><mixed-citation>
Pierret, A., Doussan, C., Garrigues, E., and Mc Kirby, J.: Observing plant
roots in their environment: current imaging options and specific contribution
of two-dimensional approaches, Agronomie, 23, 471–479,
<a href="http://dx.doi.org/10.1051/agro:2003019" target="_blank">doi:10.1051/agro:2003019</a>, 2003.
</mixed-citation></ref-html>
<ref-html id="bib1.bib73"><label>Prodan and Prodan(1999)</label><mixed-citation>
Prodan, C. and Prodan, E.: The dielectric behaviour of living cell
suspensions,
J. Phys. D, 32, 335–343, 1999.
</mixed-citation></ref-html>
<ref-html id="bib1.bib74"><label>Prodan et al.(2008)Prodan, Prodan, and Miller Jr</label><mixed-citation>
Prodan, E., Prodan, C., and Miller Jr., J.: The dielectric response of
spherical
live cells in suspension: an analytic solution, Biophys. J., 95,
4174–4182, <a href="http://dx.doi.org/10.1529/biophysj.108.137042" target="_blank">doi:10.1529/biophysj.108.137042</a>, 2008.
</mixed-citation></ref-html>
<ref-html id="bib1.bib75"><label>Razilov and Dukhin(1995)</label><mixed-citation>
Razilov, I. and Dukhin, S.: Simultaneous influence of concentration
polarization of the diffuse layer and polarization of the Stern layer
according to the mechanism of bound counterions at arbitrary magnitudes of
the relaxation parameter, Colloid J. Russ. Acad., 57, 364–371, 1995.
</mixed-citation></ref-html>
<ref-html id="bib1.bib76"><label>Repo et al.(2012)Repo, Cao, Silvennoinen, and
Ozier-Lafontaine</label><mixed-citation>
Repo, T., Cao, Y., Silvennoinen, R., and Ozier-Lafontaine, H.: Electrical
impedance spectroscopy and roots, in: Measuring Roots, edited by: Mancuso, S., Springer,
25–49, <a href="http://dx.doi.org/10.1007/978-3-642-22067-8_2" target="_blank">doi:10.1007/978-3-642-22067-8_2</a>, 2012.
</mixed-citation></ref-html>
<ref-html id="bib1.bib77"><label>Repo et al.(2014)Repo, Korhonen, Laukkanen, Lehto, and
Silvennoinen</label><mixed-citation>
Repo, T., Korhonen, A., Laukkanen, M., Lehto, T., and Silvennoinen, R.:
Detecting mycorrhizal colonisation in Scots pine roots using electrical
impedance spectra, Biosyst. Eng., 121, 139–149,
<a href="http://dx.doi.org/10.1016/j.biosystemseng.2014.02.014" target="_blank">doi:10.1016/j.biosystemseng.2014.02.014</a>, 2014.
</mixed-citation></ref-html>
<ref-html id="bib1.bib78"><label>Revil and Florsch(2010)</label><mixed-citation>
Revil, A. and Florsch, N.: Determination of permeability from spectral
induced
polarization in granular media, Geophys. J. Int., 181,
1480–1498, <a href="http://dx.doi.org/10.1111/j.1365-246X.2010.04573.x" target="_blank">doi:10.1111/j.1365-246X.2010.04573.x</a>, 2010.
</mixed-citation></ref-html>
<ref-html id="bib1.bib79"><label>Revil et al.(2012)Revil, Karaoulis, Johnson, and
Kemna</label><mixed-citation>
Revil, A., Karaoulis, M., Johnson, T., and Kemna, A.: Review: Some
low-frequency electrical methods for subsurface characterization and
monitoring in hydrogeology, Hydrogeol. J., 20, 617–658,
<a href="http://dx.doi.org/10.1007/s10040-011-0819-x" target="_blank">doi:10.1007/s10040-011-0819-x</a>, 2012.
</mixed-citation></ref-html>
<ref-html id="bib1.bib80"><label>Revil et al.(2014)Revil, Florsch, and Camerlynck</label><mixed-citation>
Revil, A., Florsch, N., and Camerlynck, C.: Spectral induced polarization
porosimetry, Geophys. J. Int., 198, 1016–1033,
<a href="http://dx.doi.org/10.1093/gji/ggu180" target="_blank">doi:10.1093/gji/ggu180</a>, 2014.
</mixed-citation></ref-html>
<ref-html id="bib1.bib81"><label>Rossi et al.(2011)Rossi, Amato, Bitella, Bochicchio, Ferreira Gomes,
Lovelli, Martorella, and Favale</label><mixed-citation>
Rossi, R., Amato, M., Bitella, G., Bochicchio, R., Ferreira Gomes, J.,
Lovelli,
S., Martorella, E., and Favale, P.: Electrical resistivity tomography as a
non-destructive method for mapping root biomass in an orchard, Eur.
J. Soil Sci., 62, 206–215,
<a href="http://dx.doi.org/10.1111/j.1365-2389.2010.01329.x" target="_blank">doi:10.1111/j.1365-2389.2010.01329.x</a>, 2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib82"><label>Schraut et al.(2005)Schraut, Heilmeier, and
Hartung</label><mixed-citation>
Schraut, D., Heilmeier, H., and Hartung, W.: Radial transport of water and
abscisic acid (ABA) in roots of Zea mays under conditions of nutrient
deficiency, J. Exp. Bot., 56, 879–886,
<a href="http://dx.doi.org/10.1093/jxb/eri080" target="_blank">doi:10.1093/jxb/eri080</a>, 2005.
</mixed-citation></ref-html>
<ref-html id="bib1.bib83"><label>Schwan(1957)</label><mixed-citation>
Schwan, H.: Electrical Properties of Tissue and Cell Suspensions, Vol. 5 of
Advances in Biological and Medical Physics, Elsevier, 147–209,
<a href="http://dx.doi.org/10.1016/B978-1-4832-3111-2.50008-0" target="_blank">doi:10.1016/B978-1-4832-3111-2.50008-0</a>,
1957.
</mixed-citation></ref-html>
<ref-html id="bib1.bib84"><label>Schwarz(1962)</label><mixed-citation>
Schwarz, G.: A theory of the low-frequency dielectric dispersion of colloidal
particles in electrolyte solution, J. Phys. Chem., 66,
2636–2642, <a href="http://dx.doi.org/10.1021/j100818a067" target="_blank">doi:10.1021/j100818a067</a>, 1962.
</mixed-citation></ref-html>
<ref-html id="bib1.bib85"><label>Singha et al.(2014)Singha, Day-Lewis, Johnson, and
Slater</label><mixed-citation>
Singha, K., Day-Lewis, F., Johnson, T., and Slater, L.: Advances in
interpretation of subsurface processes with time-lapse electrical imaging,
Hydrol. Process., 29, 1549–1576,  <a href="http://dx.doi.org/10.1002/hyp.10280" target="_blank">doi:10.1002/hyp.10280</a>, 2014.
</mixed-citation></ref-html>
<ref-html id="bib1.bib86"><label>Srayeddin and Doussan(2009)</label><mixed-citation>
Srayeddin, I. and Doussan, C.: Estimation of the spatial variability of root
water uptake of maize and sorghum at the field scale by electrical
resistivity tomography, Plant Soil, 319, 185–207,
<a href="http://dx.doi.org/10.1007/s11104-008-9860-5" target="_blank">doi:10.1007/s11104-008-9860-5</a>, 2009.
</mixed-citation></ref-html>
<ref-html id="bib1.bib87"><label>Tarasov and Titov(2013)</label><mixed-citation>
Tarasov, A. and Titov, K.: On the use of the Cole-Cole equations in spectral
induced polarization, Geophys. J. Int., 195, 352–356,
<a href="http://dx.doi.org/10.1093/gji/ggt251" target="_blank">doi:10.1093/gji/ggt251</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib88"><label>Tinker and Nye(2000)</label><mixed-citation>
Tinker, P. and Nye, P.: Solute movement in the rhizosphere, Oxford University
Press, Inc., New York, <a href="http://dx.doi.org/10.1046/j.1365-2389.2001.00418-2.x" target="_blank">doi:10.1046/j.1365-2389.2001.00418-2.x</a>,
2000.
</mixed-citation></ref-html>
<ref-html id="bib1.bib89"><label>Titov et al.(2002)Titov, Komarov, Tarasov, and
Levitski</label><mixed-citation>
Titov, K., Komarov, V., Tarasov, V., and Levitski, A.: Theoretical and
experimental study of time domain-induced polarization in water-saturated
sands, J. Appl. Geophys., 50, 417–433,
<a href="http://dx.doi.org/10.1016/S0926-9851(02)00168-4" target="_blank">doi:10.1016/S0926-9851(02)00168-4</a>, 2002.
</mixed-citation></ref-html>
<ref-html id="bib1.bib90"><label>Uhlmann and Hakim(1971)</label><mixed-citation>
Uhlmann, D. and Hakim, R.: Derivation of distribution functions from
relaxation
data, J. Phys. Chem. Solids, 32, 2652–2655,
<a href="http://dx.doi.org/10.1016/S0022-3697(71)80114-2" target="_blank">doi:10.1016/S0022-3697(71)80114-2</a>, 1971.
</mixed-citation></ref-html>
<ref-html id="bib1.bib91"><label>Urban et al.(2011)Urban, Bequet, and Mainiero</label><mixed-citation>
Urban, J., Bequet, R., and Mainiero, R.: Assessing the applicability of the
earth impedance method for in situ studies of tree root systems, J.
Exp. Bot., 62, 1857–1869, <a href="http://dx.doi.org/10.1093/jxb/erq370" target="_blank">doi:10.1093/jxb/erq370</a>, 2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib92"><label>Walker(1965)</label><mixed-citation>
Walker, J.: Electrical AC resistance and capacitance of Zea mays L, Plant
Soil, 23, 270–274, <a href="http://dx.doi.org/10.1007/BF01358354" target="_blank">doi:10.1007/BF01358354</a>, 1965.
</mixed-citation></ref-html>
<ref-html id="bib1.bib93"><label>Wang et al.(2009)Wang, Zhou, Li, and Li</label><mixed-citation>
Wang, P., Zhou, D.-M., Li, L.-Z., and Li, D.-D.: What role does cell membrane
surface potential play in ion-plant interactions, Plant Signaling &amp;
Behavior, 4, 42–43, <a href="http://dx.doi.org/10.4161/psb.4.1.7270" target="_blank">doi:10.4161/psb.4.1.7270</a>, 2009.
</mixed-citation></ref-html>
<ref-html id="bib1.bib94"><label>Wang et al.(2011)Wang, Kinraide, Zhou, K., and
Peijnenburg</label><mixed-citation>
Wang, P., Kinraide, T., Zhou, D., K., P. M., and Peijnenburg, W.: Plasma
membrane surface potential: dual effects upon ion uptake and toxicity, Plant
Physiol., 155, 808–820, <a href="http://dx.doi.org/10.1104/pp.110.165985" target="_blank">doi:10.1104/pp.110.165985</a>, 2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib95"><label>Wang et al.(2013)Wang, Kinraide, Wang, Zhou, and
Hao</label><mixed-citation>
Wang, Y.-M., Kinraide, T., Wang, P., Zhou, D., and Hao, X.: Modeling
rhizotoxicity and uptake of Zn and Co singly and in binary mixture in wheat
in terms of the cell membrane surface electrical potential, Environ.
Sci. Technol., 47, 2831–2838, <a href="http://dx.doi.org/10.1021/es3022107" target="_blank">doi:10.1021/es3022107</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib96"><label>Weigand and Kemna(2016)</label><mixed-citation>
Weigand, M. and Kemna, A.: Debye decomposition of time-lapse spectral induced
polarisation data, Comput. Geosci., 86, 34–45,
<a href="http://dx.doi.org/10.1016/j.cageo.2015.09.021" target="_blank">doi:10.1016/j.cageo.2015.09.021</a>, 2016.
</mixed-citation></ref-html>
<ref-html id="bib1.bib97"><label>Weigand and Kemna(2017)</label><mixed-citation>
Weigand, M. and Kemna, A.: Data and results for manuscript “Multi-frequency
electrical impedance tomography as a non-invasive tool to characterize and
monitor crop root systems”, Data set, Zenodo, <a href="http://dx.doi.org/10.5281/zenodo.260087" target="_blank">doi:10.5281/zenodo.260087</a>,
2017.
</mixed-citation></ref-html>
<ref-html id="bib1.bib98"><label>Whalley et al.(2017)Whalley, Binley, Watts, Shanahan, Dodd, Ober,
Ashton, Webster, White, and Hawkesford</label><mixed-citation>
Whalley, W., Binley, A., Watts, C., Shanahan, P., Dodd, I., Ober, E., Ashton,
R., Webster, C., White, R., and Hawkesford, M. J.: Methods to estimate
changes in soil water for phenotyping root activity in the field, Plant
Soil,  1–16, <a href="http://dx.doi.org/10.1007/s11104-016-3161-1" target="_blank">doi:10.1007/s11104-016-3161-1</a>, 2017.

</mixed-citation></ref-html>
<ref-html id="bib1.bib99"><label>Willatt et al.(1978)Willatt, Struss, and Taylor</label><mixed-citation>
Willatt, S., Struss, R., and Taylor, H.: In situ root studies using neutron
radiography, Agron. J., 70, 581–586,
<a href="http://dx.doi.org/10.2134/agronj1978.00021962007000040016x" target="_blank">doi:10.2134/agronj1978.00021962007000040016x</a>, 1978.
</mixed-citation></ref-html>
<ref-html id="bib1.bib100"><label>Zanetti et al.(2011)Zanetti, Weller, Vennetier, and
Mériaux</label><mixed-citation>
Zanetti, C., Weller, A., Vennetier, M., and Mériaux, P.: Detection of
buried tree root samples by using geoelectrical measurements: a laboratory
experiment, Plant Soil, 339, 273–283, <a href="http://dx.doi.org/10.1007/s11104-010-0574-0" target="_blank">doi:10.1007/s11104-010-0574-0</a>,
2011.
</mixed-citation></ref-html>
<ref-html id="bib1.bib101"><label>Zhao et al.(2013)Zhao, Zimmermann, Huisman, Treichel, Wolters, van
Waasen, and Kemna</label><mixed-citation>
Zhao, Y., Zimmermann, E., Huisman, J., Treichel, A., Wolters, B., van Waasen,
S., and Kemna, A.: Broadband EIT borehole measurements with high phase
accuracy using numerical corrections of electromagnetic coupling effects,
Meas. Sci. Technol., 24, 085005,
<a href="http://dx.doi.org/10.1088/0957-0233/24/8/085005" target="_blank">doi:10.1088/0957-0233/24/8/085005</a>, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib102"><label>Zimmermann et al.(2008)Zimmermann, Kemna, Berwix, Glaas, and
Vereecken</label><mixed-citation>
Zimmermann, E., Kemna, A., Berwix, J., Glaas, W., and Vereecken, H.: EIT
measurement system with high phase accuracy for the imaging of spectral
induced polarization properties of soils and sediments, Meas. Sci.
Technol., 19, 094010, <a href="http://dx.doi.org/10.1088/0957-0233/19/9/094010" target="_blank">doi:10.1088/0957-0233/19/9/094010</a>, 2008.
</mixed-citation></ref-html>--></article>
