<?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" xml:lang="en" dtd-version="3.0"><?xmltex \makeatother\@nolinetrue\makeatletter?>
  <front>
    <journal-meta><journal-id journal-id-type="publisher">HESS</journal-id><journal-title-group>
    <journal-title>Hydrology and Earth System Sciences</journal-title>
    <abbrev-journal-title abbrev-type="publisher">HESS</abbrev-journal-title><abbrev-journal-title abbrev-type="nlm-ta">Hydrol. Earth Syst. Sci.</abbrev-journal-title>
  </journal-title-group><issn pub-type="epub">1607-7938</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/hess-23-3593-2019</article-id><title-group><article-title>A soil non-aqueous phase liquid (NAPL) flushing laboratory experiment based on measuring the dielectric properties of soil–organic mixtures via time domain reflectometry (TDR)</article-title><alt-title>A soil non-aqueous phase liquid (NAPL) flushing laboratory experiment</alt-title>
      </title-group><?xmltex \runningtitle{A soil non-aqueous phase liquid~(NAPL) flushing laboratory experiment}?><?xmltex \runningauthor{A.~Comegna et al.}?>
      <contrib-group>
        <contrib contrib-type="author" corresp="yes" rid="aff1">
          <name><surname>Comegna</surname><given-names>Alessandro</given-names></name>
          <email>alessandro.comegna@unibas.it</email>
        <ext-link>https://orcid.org/0000-0002-8501-5046</ext-link></contrib>
        <contrib contrib-type="author" corresp="no" rid="aff1">
          <name><surname>Coppola</surname><given-names>Antonio</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff2">
          <name><surname>Dragonetti</surname><given-names>Giovanna</given-names></name>
          
        </contrib>
        <contrib contrib-type="author" corresp="no" rid="aff3">
          <name><surname>Sommella</surname><given-names>Angelo</given-names></name>
          
        </contrib>
        <aff id="aff1"><label>1</label><institution>School of Agricultural Forestry, Food, and Environmental Sciences (SAFE), University of Basilicata, Potenza, Italy</institution>
        </aff>
        <aff id="aff2"><label>2</label><institution>Mediterranean Agronomic Institute, Land and Water Division, IAMB,
Bari, Italy</institution>
        </aff>
        <aff id="aff3"><label>3</label><institution>Department of Agriculture, University of Naples Federico II,
Naples, Italy</institution>
        </aff>
      </contrib-group>
      <author-notes><corresp id="corr1">Alessandro Comegna (alessandro.comegna@unibas.it)</corresp></author-notes><pub-date><day>4</day><month>September</month><year>2019</year></pub-date>
      
      <volume>23</volume>
      <issue>9</issue>
      <fpage>3593</fpage><lpage>3602</lpage>
      <history>
        <date date-type="received"><day>4</day><month>April</month><year>2019</year></date>
           <date date-type="rev-request"><day>29</day><month>April</month><year>2019</year></date>
           <date date-type="rev-recd"><day>25</day><month>July</month><year>2019</year></date>
           <date date-type="accepted"><day>6</day><month>August</month><year>2019</year></date>
      </history>
      <permissions>
        <copyright-statement>Copyright: © 2019 Alessandro Comegna et al.</copyright-statement>
        <copyright-year>2019</copyright-year>
      <license license-type="open-access"><license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><self-uri xlink:href="https://hess.copernicus.org/articles/23/3593/2019/hess-23-3593-2019.html">This article is available from https://hess.copernicus.org/articles/23/3593/2019/hess-23-3593-2019.html</self-uri><self-uri xlink:href="https://hess.copernicus.org/articles/23/3593/2019/hess-23-3593-2019.pdf">The full text article is available as a PDF file from https://hess.copernicus.org/articles/23/3593/2019/hess-23-3593-2019.pdf</self-uri>
      <abstract><title>Abstract</title>
    <p id="d1e121">The term non-aqueous phase liquid (NAPL) refers to a group of organic compounds with scarce solubility in water. They are the products of various human activities and may be accidentally introduced into the soil system. Given their toxicity level and high mobility, NAPLs constitute a serious geo-environmental problem. Contaminant distribution in the soil and groundwater contains fundamental information for the remediation of polluted
soil sites. The present research explored the possible employment of time
domain reflectometry (TDR) to estimate pollutant removal in a silt-loam soil that was primarily contaminated with a corn oil as a light NAPL and then flushed with different washing solutions. Known mixtures of soil and NAPL were prepared in the laboratory to achieve soil specimens with varying
pollution levels. The prepared soil samples were repacked into plastic
cylinders and then placed in testing cells. Washing solutions were then
injected upward into the contaminated sample, and both the quantity of
remediated NAPL and the bulk dielectric permittivity of the soil sample were
determined. The above data were also used to calibrate and validate a
dielectric model (the <inline-formula><mml:math id="M1" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> mixing model) which permits the volumetric
NAPL content (<inline-formula><mml:math id="M2" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>; m<inline-formula><mml:math id="M3" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">3</mml:mn></mml:msup></mml:math></inline-formula> m<inline-formula><mml:math id="M4" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>) within the contaminated sample to be determined and quantified during the different decontamination stages. Our results demonstrate that during a decontamination process, the TDR device is NAPL-sensitive: the dielectric permittivity of the medium increases as the NAPL volume decreases. Moreover, decontamination progression can be monitored using a simple (one-parameter) mixing model.</p>
  </abstract>
    </article-meta>
  </front>
<body>
      

      <?xmltex \hack{\newpage}?>
<sec id="Ch1.S1" sec-type="intro">
  <label>1</label><title>Introduction</title>
      <p id="d1e174">Soil and groundwater contamination with non-aqueous phase liquid (NAPL) from point or nonpoint sources
is a severe problem of considerable complexity (Fitts, 2002; Fetter, 1993).
The repercussions concern not only the deterioration of the soil's physical,
mechanical and chemical properties but also account for a potentially
severe hazard to the well-being of humans and other living species (Freeze,
2000).</p>
      <p id="d1e177">Soil flushing is the technical procedure used for treating polluted soils
with water, surfactants and co-solvents (such as methanol, ethanol and
propanols). Surfactant-enhanced flushing was developed from the conventional
pump-and-treat method. The success of this approach is related to the
capacity of such chemical compounds to greatly enhance the aqueous
solubility of oils (Pennell et al., 1994; Parnian and Ayatollahi, 2008).</p>
      <p id="d1e180">There is high interfacial tension between NAPL and water molecules that
makes water a non-efficient cleaning material in removing NAPL from the
soil. Instead, surfactants and co-solvent agents can promote the enhanced
removal of NAPL from the subsurface through mobilization and solubilization
(Martel et al., 1998; Rinaldi and Francisca, 2006; Parnian and Ayatollahi,
2008).</p>
      <p id="d1e183">Primary remediation entails the removal of the NAPL free phase by pumping.
This extraction mechanism returns appreciable effects if there is a region
of high NAPL saturation. After primary pumping, a considerable portion of
NAPL remains constrained within the soil as capillary forces overcome
viscous and buoyancy forces. This discontinuous<?pagebreak page3594?> NAPL phase is referred to as
trapped residual NAPL (or NAPL residual saturation), and its remediation is
referred to as secondary remediation (Parnian and Ayatollahi, 2008).
Residual NAPL is a long-term source of soil and groundwater pollution
(Mercer and Cohen, 1990; Troung Hong and Bettahar, 2000).</p>
      <p id="d1e187">To develop powerful decontamination procedures, the characterization of
polluted soils is required. Practices usually employed to characterize
polluted soil sites are coring, soil sampling and the installation of
monitoring wells for the collection of water samples from aquifers (Mercer
and Cohen, 1990). Since the above procedures are costly, different
dielectric techniques can be used to detect organic contaminants in soils.
The most widely accepted geophysical technique, based on the principle of
electromagnetic wave (EMW) propagation, is the ground-penetrating radar (GPR;
Knight, 2001). Redman et al. (1991) described some field experiments in the
application of GPR to detect NAPL plumes.</p>
      <p id="d1e190">Rinaldi and Francisca (2006) used a coaxial impedance dielectric
reflectometry (CIDR) technique to measure the complex dielectric
permittivity in sands contaminated by a paraffin oil. Their research into
the dielectric behavior of NAPL-contaminated soils during a decontamination
process mainly focused on the removal efficiency of different washing
solutions and on the spectral response of the contaminated medium during
the various tests conducted.</p>
      <p id="d1e193">Time domain reflectometry (TDR) is a further geophysical device based on electromagnetic wave (EMW)
principles that can also be used for this purpose (Endres and Redman,
1993; Redman and DeRyck, 1994; Mohamed and Said, 2005; Moroizumi and
Sasaki, 2006; Francisca and Montoro, 2012). Few experiments have been
conducted coupling the TDR technique and NAPL. In these studies estimation
of NAPLs using TDR measurements of dielectric properties relies greatly on
various mixing models relating the measured dielectric permittivity to the
volume fractions of the pore fluids and various soil phases such as solid,
water, air and NAPLs (van Dam et al., 2005).</p>
      <p id="d1e196">Some interesting results were achieved by Persson and Berndtsson (2002)
whilst investigating the influence of different LNAPLs (i.e., light NAPLs) on TDR measurements
in a homogeneous silica sand under saturated and unsaturated soil
conditions. Measurements of both dielectric permittivity and electrical
conductivity allowed a method to be developed (the two-step method) which
measured the dielectric properties of the system against the amount of NAPL
in soils. Comegna et al. (2016) developed a general TDR-based methodology
for evaluating the correlations between the dielectric response and the NAPL
content in variable saturated soils with different textures and pedological
characteristics.</p>
      <p id="d1e199">The purpose of this study was as follows: (i) to investigate a possible
extension of TDR technology to assess the effects of NAPL removal in soils
and (ii) revisit, on the basis of the acquired data and the experimental
results, a dielectric model to predict, “in real time”, the volumetric amounts of NAPL (<inline-formula><mml:math id="M5" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) within the contaminated soil during the decontamination process.</p>
</sec>
<sec id="Ch1.S2">
  <label>2</label><title>Theoretical concepts of TDR</title>
      <p id="d1e221">TDR is a geophysical technique employed to determine the dielectric
permittivity of liquids and solids (Ferrè and Topp, 2002, described this
method in detail). In general, the bulk dielectric permittivity is a complex
term (<inline-formula><mml:math id="M6" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">r</mml:mi><mml:mo>*</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula>), which may be expressed as follows (Robinson et al., 2003):
          <disp-formula id="Ch1.E1" content-type="numbered"><label>1</label><mml:math id="M7" display="block"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">r</mml:mi><mml:mo>*</mml:mo></mml:msubsup><mml:mo>=</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">r</mml:mi><mml:mo>′</mml:mo></mml:msubsup><mml:mo>-</mml:mo><mml:mi>j</mml:mi><mml:mfenced open="[" close="]"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">r</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msubsup><mml:mo>+</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi mathvariant="italic">σ</mml:mi><mml:mrow><mml:mi mathvariant="italic">ω</mml:mi><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M8" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">r</mml:mi><mml:mo>′</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula> is the real part of dielectric permittivity, which gives the energy stored in the dielectrics at a certain frequency and
temperature, and <inline-formula><mml:math id="M9" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">r</mml:mi><mml:mrow><mml:mo>′</mml:mo><mml:mo>′</mml:mo></mml:mrow></mml:msubsup></mml:mrow></mml:math></inline-formula> is the imaginary part due to
relaxations. The zero-frequency conductivity <inline-formula><mml:math id="M10" display="inline"><mml:mi mathvariant="italic">σ</mml:mi></mml:math></inline-formula>, the angle frequency <inline-formula><mml:math id="M11" display="inline"><mml:mi mathvariant="italic">ω</mml:mi></mml:math></inline-formula>, the imaginary number <inline-formula><mml:math id="M12" display="inline"><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msqrt></mml:mrow></mml:math></inline-formula> and the permittivity <inline-formula><mml:math id="M13" 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> in free space contribute to defining <inline-formula><mml:math id="M14" display="inline"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">r</mml:mi><mml:mo>*</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula>.</p>
      <p id="d1e375">When the frequency of a TDR cable tester ranges between 200 MHz and 1.5 GHz,
dielectric losses can be considered minimal (Heimovaara, 1994), and the bulk
dielectric permittivity <inline-formula><mml:math id="M15" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> (<inline-formula><mml:math id="M16" display="inline"><mml:mo lspace="0mm">≅</mml:mo></mml:math></inline-formula> the real part of
permittivity) of a probe of length <inline-formula><mml:math id="M17" display="inline"><mml:mi>L</mml:mi></mml:math></inline-formula> is determined from the propagation
velocity <inline-formula><mml:math id="M18" display="inline"><mml:mrow><mml:mi>v</mml:mi><mml:mo>(</mml:mo><mml:mo>=</mml:mo><mml:mn mathvariant="normal">2</mml:mn><mml:mi>L</mml:mi><mml:mo>/</mml:mo><mml:mi>t</mml:mi><mml:mo>)</mml:mo></mml:mrow></mml:math></inline-formula> of an electromagnetic wave along the wave guide across the investigated medium by the following expression:
          <disp-formula id="Ch1.E2" content-type="numbered"><label>2</label><mml:math id="M19" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">b</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mfenced close=")" open="("><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mi>c</mml:mi><mml:mi>v</mml:mi></mml:mfrac></mml:mstyle></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M20" display="inline"><mml:mrow><mml:mi>c</mml:mi><mml:mo>(</mml:mo><mml:mo>=</mml:mo><mml:mn mathvariant="normal">3</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mn mathvariant="normal">10</mml:mn><mml:mn mathvariant="normal">8</mml:mn></mml:msup></mml:mrow></mml:math></inline-formula> m s<inline-formula><mml:math id="M21" 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>) is the velocity of an electromagnetic
wave in vacuum (Topp et al., 1980) and <inline-formula><mml:math id="M22" display="inline"><mml:mi>t</mml:mi></mml:math></inline-formula> is travel time, i.e., the time
required by the generated signal to go back and forth through the TDR probe
of length <inline-formula><mml:math id="M23" display="inline"><mml:mi>L</mml:mi></mml:math></inline-formula> (m). This can be calculated as follows:
          <disp-formula id="Ch1.E3" content-type="numbered"><label>3</label><mml:math id="M24" display="block"><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mn mathvariant="normal">2</mml:mn><mml:mi>L</mml:mi></mml:mrow><mml:mi>c</mml:mi></mml:mfrac></mml:mstyle><mml:msqrt><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">b</mml:mi></mml:msub></mml:mrow></mml:msqrt><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
        The direct dependence of the signal's travel time <inline-formula><mml:math id="M25" display="inline"><mml:mi>t</mml:mi></mml:math></inline-formula> on soil dielectric
permittivity is expressed by Eq. (3).</p>
</sec>
<sec id="Ch1.S3">
  <label>3</label><title>Estimating volumetric NAPL content during a decontamination process in soils</title>
      <p id="d1e541">Dielectric mixing models, in their classical application, have been proposed
to estimate the bulk dielectric permittivity of a multi-phase medium, that
is, a combination of three or four dielectric phases, and to couple the
dielectric permittivity of the medium to the dielectric permittivity of each
single phase (Hilhorst, 1998). Recently, after analyzing the effects of
organic contaminants on soil dielectric properties, the above models were
further developed to estimate the dielectric properties of NAPL-polluted
soils (Redman et al.,<?pagebreak page3595?> 1991; Persson and Berndtsson, 2002; Francisca and
Montoro, 2012; Comegna et al., 2013a, 2016, 2017).</p>
      <p id="d1e544">Based on such models, in the present study, we analyze the possibility of
predicting the correlations between the volumetric contents of
NAPL (<inline-formula><mml:math id="M26" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) and the dielectric response (<inline-formula><mml:math id="M27" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">b</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) of contaminated soil during the progression of a steady-state remediation process.</p>
      <p id="d1e569">In the present research, we chose the so-called <inline-formula><mml:math id="M28" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> model (Birchack et
al., 1974; Knight and Endres, 1990; Roth et al., 1990):
          <disp-formula id="Ch1.E4" content-type="numbered"><label>4</label><mml:math id="M29" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">b</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msup><mml:mfenced close="]" open="["><mml:mrow><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>n</mml:mi></mml:munderover><mml:msub><mml:mi>V</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi>i</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow></mml:mfenced><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>/</mml:mo><mml:mi mathvariant="italic">α</mml:mi></mml:mrow></mml:msup><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M30" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the volume and <inline-formula><mml:math id="M31" 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> is the permittivity of
each component of the complex medium; the exponent <inline-formula><mml:math id="M32" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> is a fitting
parameter (<inline-formula><mml:math id="M33" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> varies between <inline-formula><mml:math id="M34" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> and 1), which may be related to the
internal structure of the investigated medium (Hilhorst, 1998; Coppola et
al., 2013, 2015). Under the following hypothesis, (i) the soil is homogeneous from a textural point of view, and (ii) when the soil porosity (<inline-formula><mml:math id="M35" display="inline"><mml:mi mathvariant="italic">ϕ</mml:mi></mml:math></inline-formula>) is constant, Eq. (4) was reformulated for our purposes.</p>
      <p id="d1e682">For mixtures of soil (<inline-formula><mml:math id="M36" display="inline"><mml:mi>s</mml:mi></mml:math></inline-formula>) saturated with a certain amount of washing solution (ws), in rearranging the model formulation of Rinaldi and Francisca (2006), the <inline-formula><mml:math id="M37" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> model yields the following:
          <disp-formula id="Ch1.E5" content-type="numbered"><label>5</label><mml:math id="M38" display="block"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mfenced close="]" open="["><mml:mrow><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>)</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">ws</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M39" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the soil-washing solution permittivity, and <inline-formula><mml:math id="M40" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M41" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">ws</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are the permittivities of soil particles and washing solutions, respectively. By the same token, for soil organic (<inline-formula><mml:math id="M42" display="inline"><mml:mi>s</mml:mi></mml:math></inline-formula>-NAPL) compounds at saturation, the <inline-formula><mml:math id="M43" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> model can be expressed as follows:
          <disp-formula id="Ch1.E6" content-type="numbered"><label>6</label><mml:math id="M44" display="block"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mfenced close="]" open="["><mml:mrow><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>)</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M45" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula> is the permittivity of the soil–NAPL mixture, and <inline-formula><mml:math id="M46" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the oil permittivity.</p>
      <p id="d1e884">A medium consisting of soil particles, washing solution and NAPL (<inline-formula><mml:math id="M47" display="inline"><mml:mi>s</mml:mi></mml:math></inline-formula>-ws-NAPL) can be viewed as a mix of soil-washing solution (Eq. 5) and soil–NAPL (Eq. 6):
          <disp-formula id="Ch1.E7" content-type="numbered"><label>7</label><mml:math id="M48" display="block"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:mfenced open="[" close="]"><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mi mathvariant="italic">β</mml:mi><mml:mo>)</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow></mml:mfenced><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M49" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula> is the relative volume of NAPL contained in the whole fluid
phase,
          <disp-formula id="Ch1.E8" content-type="numbered"><label>8</label><mml:math id="M50" display="block"><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">ws</mml:mi></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">f</mml:mi></mml:msub></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:math></disp-formula>
        where <inline-formula><mml:math id="M51" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">f</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is the volumetric fluid content (m<inline-formula><mml:math id="M52" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">3</mml:mn></mml:msup></mml:math></inline-formula> m<inline-formula><mml:math id="M53" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>), the sum of the volumetric washing solution content (<inline-formula><mml:math id="M54" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">ws</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) and volumetric NAPL
content (<inline-formula><mml:math id="M55" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>); <inline-formula><mml:math id="M56" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula> varies between 0 (i.e., a soil-washing solution mixture) and 1 (i.e., a soil–NAPL mixture).</p>
      <p id="d1e1080"><?xmltex \hack{\newpage}?>To estimate <inline-formula><mml:math id="M57" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>, Eq. (7) is first reformulated in terms of <inline-formula><mml:math id="M58" display="inline"><mml:mi mathvariant="italic">β</mml:mi></mml:math></inline-formula>:
          <disp-formula id="Ch1.E9" content-type="numbered"><label>9</label><mml:math id="M59" display="block"><mml:mtable rowspacing="0.2ex" columnspacing="1em" class="split" displaystyle="true" columnalign="right left"><mml:mtr><mml:mtd><mml:mrow><mml:mi mathvariant="italic">β</mml:mi></mml:mrow></mml:mtd><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow></mml:mfrac></mml:mstyle></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd/><mml:mtd><mml:mrow><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>)</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">ws</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:mfenced open="(" close=")"><mml:mrow><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>)</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">ws</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow></mml:mfenced><mml:mo>-</mml:mo><mml:mfenced close=")" open="("><mml:mrow><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>)</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow></mml:mfenced></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
        Substituting Eq. (8) into Eq. (9), and considering that for a saturated medium, the volumetric fluid content is equal to soil porosity
(i.e., <inline-formula><mml:math id="M60" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">f</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi></mml:mrow></mml:math></inline-formula>), <inline-formula><mml:math id="M61" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> can be calculated as follows:
          <disp-formula id="Ch1.E10" content-type="numbered"><label>10</label><mml:math id="M62" display="block"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:mo>(</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:mo>)</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">s</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:mi mathvariant="italic">ϕ</mml:mi><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">ws</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow><mml:mrow><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">ws</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>-</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>.</mml:mo></mml:mrow></mml:math></disp-formula>
        Equation (10) correlates the dependence of volumetric NAPL content with soil
porosity; <inline-formula><mml:math id="M63" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> can be estimated (within the contaminated soil) during the progression of a remediation process once the dielectric
permittivity of the soil-contaminated mixture (<inline-formula><mml:math id="M64" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>)
is known.</p>
</sec>
<sec id="Ch1.S4">
  <label>4</label><title>Materials and methods</title>
<sec id="Ch1.S4.SS1">
  <label>4.1</label><title>Soil and fluid properties</title>
      <p id="d1e1452">A silt-loam anthrosol (IUSS Working Group WRB, 2006) from the region of
Puglia (Italy) was used for this study. The soil texture was measured by
means of the hydrometer method (Day, 1965), while the Walkley–Black
procedure (Allison, 1965) was used to determine soil organic C content. The
method developed by Miller and Curtis (2006) was used to measure soil
electrical conductivity (EC<inline-formula><mml:math id="M65" display="inline"><mml:msub><mml:mi/><mml:mi mathvariant="normal">w</mml:mi></mml:msub></mml:math></inline-formula>), while soil pH was determined on the basis of a <inline-formula><mml:math id="M66" display="inline"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>:</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> soil-to-water ratio (Eckert, 1988). In textural terms, the soil comprised 15.7 % sand, 11.6 % clay and 72.4 % silt. Soil porosity was 0.57 %, organic content was 1.84 %, EC<inline-formula><mml:math id="M67" display="inline"><mml:msub><mml:mi/><mml:mi mathvariant="normal">w</mml:mi></mml:msub></mml:math></inline-formula> was 0.17 dS m<inline-formula><mml:math id="M68" 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> and soil pH was 8.40.</p>
      <p id="d1e1497">The NAPL employed for the laboratory tests was corn oil (<inline-formula><mml:math id="M69" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">3.2</mml:mn></mml:mrow></mml:math></inline-formula>; EC<inline-formula><mml:math id="M70" display="inline"><mml:mrow><mml:msub><mml:mi/><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.055</mml:mn></mml:mrow></mml:math></inline-formula> dS m<inline-formula><mml:math id="M71" 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> at 25 <inline-formula><mml:math id="M72" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>C) with a density of 0.905 g cm<inline-formula><mml:math id="M73" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula> (at 25 <inline-formula><mml:math id="M74" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>C). Three different removal solutions were employed for soil cleaning: (a) a first solution (referred to below as wd) composed of 99 % distilled water and 1 % commercial detergent (<inline-formula><mml:math id="M75" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">d</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">9.22</mml:mn></mml:mrow></mml:math></inline-formula>, at 25 <inline-formula><mml:math id="M76" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>C); (b) a second solution (wda no. 1) composed of 90 % distilled water, 1 % commercial detergent and 9 % methanol as co-solvent (<inline-formula><mml:math id="M77" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">alcohol</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">26.13</mml:mn></mml:mrow></mml:math></inline-formula>, at 25 <inline-formula><mml:math id="M78" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>C); and (c) a third solution (wda no. 2) composed of distilled water (85 %) with commercial detergent (1 %) and methanol (14 %). The dielectric permittivity of the washing solutions, measured at 25 <inline-formula><mml:math id="M79" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>C, was <inline-formula><mml:math id="M80" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">wd</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">75.04</mml:mn></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M81" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">wda</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mi mathvariant="normal">no</mml:mi><mml:mo>.</mml:mo><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">68.98</mml:mn></mml:mrow></mml:math></inline-formula> and
<inline-formula><mml:math id="M82" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">wda</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mi mathvariant="normal">no</mml:mi><mml:mo>.</mml:mo><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">65.92</mml:mn></mml:mrow></mml:math></inline-formula>, whereas the dielectric permittivity of
the tested soil saturated with each of the three cleaning solutions was
<inline-formula><mml:math id="M83" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">soil</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">wd</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">34.59</mml:mn></mml:mrow></mml:math></inline-formula>, <inline-formula><mml:math id="M84" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">soil</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">wda</mml:mi><mml:mspace linebreak="nobreak" width="0.25em"/><mml:mi mathvariant="normal">no</mml:mi><mml:mo>.</mml:mo><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">31.04</mml:mn></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M85" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">soil</mml:mi><mml:mo>+</mml:mo><mml:mi mathvariant="normal">wda</mml:mi><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mi mathvariant="normal">no</mml:mi><mml:mo>.</mml:mo><mml:mspace width="0.25em" linebreak="nobreak"/><mml:mn mathvariant="normal">2</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">30.10</mml:mn></mml:mrow></mml:math></inline-formula>.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F1" specific-use="star"><?xmltex \currentcnt{1}?><label>Figure 1</label><caption><p id="d1e1771">Experimental setup used in the NAPL removal experiments (from
Comegna et al., 2013c).</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://hess.copernicus.org/articles/23/3593/2019/hess-23-3593-2019-f01.png"/>

        </fig>

<?xmltex \hack{\newpage}?>
</sec>
<?pagebreak page3596?><sec id="Ch1.S4.SS2">
  <label>4.2</label><title>Measurement of dielectric permittivity of soil–NAPL-contaminated samples during soil remediation</title>
<sec id="Ch1.S4.SS2.SSS1">
  <label>4.2.1</label><title>Experimental setup</title>
      <p id="d1e1797">As illustrated in Fig. 1, the experimental layout consisted of the
following: (i) a Tektronix (model 1502C) cable tester, (ii) a three-rod probe, 14.5 cm long, with a wire diameter of 0.003 m and a wire spacing of 0.02 m, introduced vertically into the soil samples, (iii) a testing cell, 0.15 m high and 0.08 m in diameter, and (iv) a peristaltic pump used for upward movement of the washing solution.</p>
</sec>
<sec id="Ch1.S4.SS2.SSS2">
  <label>4.2.2</label><title>Sample preparation and testing procedures</title>
      <p id="d1e1808">Soil was oven dried at 105 <inline-formula><mml:math id="M86" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>C and passed through a 2 mm sieve.
Known amounts of soil and oil were mixed together, shaken, and then kept for
24 h in sealed plastic bags to avoid any evaporation and ensure a uniform distribution of oil within the sample and good oil adsorption by the soil matrix. The samples were then allocated to cylindrical boxes. With a view to achieving different degrees of (initial) soil contamination, volumetric NAPL content (<inline-formula><mml:math id="M87" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) was varied from 0.05 to 0.40 (in steps of 0.05). In all, each washing solution comprised eight oil-contaminated soil samples.</p>
      <p id="d1e1831">For all experiments, the soil samples were placed in the vessels in various
steps at a bulk density of 1.13 g cm<inline-formula><mml:math id="M88" display="inline"><mml:msup><mml:mi/><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">3</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>. During TDR measurements, the
soil samples were conserved at a temperature of 25 <inline-formula><mml:math id="M89" display="inline"><mml:msup><mml:mi/><mml:mo>∘</mml:mo></mml:msup></mml:math></inline-formula>C by using a
thermostat box. Remediation was performed using an upward flux of diverse
pore volumes <inline-formula><mml:math id="M90" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> of three washing solutions (wd, wda no. 1 and wda no. 2) supplied at the rate of 90 cm<inline-formula><mml:math id="M91" display="inline"><mml:msup><mml:mi/><mml:mn mathvariant="normal">3</mml:mn></mml:msup></mml:math></inline-formula> h<inline-formula><mml:math id="M92" 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>, corresponding to a Darcian velocity of
1.8 cm h<inline-formula><mml:math id="M93" 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>. After collection of the outflow from the soil columns, the
surnatant NAPL was separated from the washing solution and the quantity of
NAPL remediated from the soil was determined.</p>
      <p id="d1e1896">The obtained data series were employed to calibrate the proposed dielectric
model of Eq. (10). An independent dataset, obtained in the same manner as the calibration dataset, was used for model validation.</p>
</sec>
</sec>
<sec id="Ch1.S4.SS3">
  <label>4.3</label><title>Numerical indices for model performance evaluation</title>
      <p id="d1e1908">The goodness of Eq. (10) was evaluated using two different criteria: (i) the mean bias error (MBE) and (ii) the model efficiency (EF), computed according
to the following relations (Legates and McCabe Jr., 1999):

                <disp-formula specific-use="align" content-type="numbered"><mml:math id="M94" display="block"><mml:mtable displaystyle="true"><mml:mlabeledtr id="Ch1.E11"><mml:mtd><mml:mtext>11</mml:mtext></mml:mtd><mml:mtd><mml:mstyle displaystyle="true" class="stylechange"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle displaystyle="true" class="stylechange"/><mml:mi mathvariant="normal">MBE</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>N</mml:mi></mml:munderover><mml:mfenced close=")" open="("><mml:mrow><mml:msub><mml:mi>E</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>O</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mfenced></mml:mrow><mml:mi>N</mml:mi></mml:mfrac></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr><mml:mlabeledtr id="Ch1.E12"><mml:mtd><mml:mtext>12</mml:mtext></mml:mtd><mml:mtd><mml:mstyle class="stylechange" displaystyle="true"/></mml:mtd><mml:mtd><mml:mrow><mml:mstyle class="stylechange" displaystyle="true"/><mml:mi mathvariant="normal">EF</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn><mml:mo>-</mml:mo><mml:mstyle displaystyle="true"><mml:mfrac style="display"><mml:mrow><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>N</mml:mi></mml:munderover><mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>E</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi>O</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow><mml:mrow><mml:munderover><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow><mml:mi>N</mml:mi></mml:munderover><mml:msup><mml:mfenced open="(" close=")"><mml:mrow><mml:msub><mml:mi>O</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>-</mml:mo><mml:mover accent="true"><mml:mi>O</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover></mml:mrow></mml:mfenced><mml:mn mathvariant="normal">2</mml:mn></mml:msup></mml:mrow></mml:mfrac></mml:mstyle><mml:mo>,</mml:mo></mml:mrow></mml:mtd></mml:mlabeledtr></mml:mtable></mml:math></disp-formula>

            where <inline-formula><mml:math id="M95" display="inline"><mml:mrow><mml:msub><mml:mi>E</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> and <inline-formula><mml:math id="M96" display="inline"><mml:mrow><mml:msub><mml:mi>O</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are respectively the expected and the observed
value, <inline-formula><mml:math id="M97" display="inline"><mml:mover accent="true"><mml:mi>O</mml:mi><mml:mo mathvariant="normal">‾</mml:mo></mml:mover></mml:math></inline-formula> is the mean of the observed data, and <inline-formula><mml:math id="M98" display="inline"><mml:mi>N</mml:mi></mml:math></inline-formula> is the number
of observations.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F2" specific-use="star"><?xmltex \currentcnt{2}?><label>Figure 2</label><caption><p id="d1e2085">Selection of experimental relationships between the measured dielectric permittivity (<inline-formula><mml:math id="M99" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>) and number of pore volumes <inline-formula><mml:math id="M100" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> under the effect of different washing solutions: (i) water–detergent (wd) and (ii) water–detergent–alcohol (wda no. 1 and wda no. 2).</p></caption>
          <?xmltex \igopts{width=497.923228pt}?><graphic xlink:href="https://hess.copernicus.org/articles/23/3593/2019/hess-23-3593-2019-f02.png"/>

        </fig>

      <p id="d1e2121">MBE measures the differences between model-simulated data and measured values
(positive MBE values are used<?pagebreak page3597?> to indicate average overprediction, while negative values indicate underprediction). The model's ability to forecast <inline-formula><mml:math id="M101" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is described by parameter EF, according to which EF <inline-formula><mml:math id="M102" display="inline"><mml:mo>=</mml:mo></mml:math></inline-formula> 1 indicates perfect accord between predicted and measured data.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F3" specific-use="star"><?xmltex \currentcnt{3}?><label>Figure 3</label><caption><p id="d1e2145">Volume of NAPL recovered (<inline-formula><mml:math id="M103" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mrow><mml:mi mathvariant="normal">NAPL</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">Rem</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>) with respect to the initial volume of NAPL present in the soil sample (<inline-formula><mml:math id="M104" display="inline"><mml:mrow><mml:msub><mml:mi>V</mml:mi><mml:mn mathvariant="normal">0</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula>) of different washing solutions (wd, wda no. 1 and wda no. 2) for different experiments (<inline-formula><mml:math id="M105" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.15</mml:mn></mml:mrow></mml:math></inline-formula>, 0.20, 0.25, 0.30, 0.35 and 0.40).</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://hess.copernicus.org/articles/23/3593/2019/hess-23-3593-2019-f03.png"/>

        </fig>

</sec>
</sec>
<sec id="Ch1.S5">
  <label>5</label><title>Results and discussion</title>
<sec id="Ch1.S5.SS1">
  <label>5.1</label><title>Influence of washing solution on NAPL removal</title>
      <p id="d1e2212">Figure 2a–f, with reference to the most representative
experimental results, reveal the influence of pore volumes <inline-formula><mml:math id="M106" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula> on evaluated
bulk dielectric permittivity (<inline-formula><mml:math id="M107" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>) for the soil specimens initially polluted with oil. As the washing solution started to remove oil, the dielectric permittivity rose due to the larger dielectric
permittivity of the flushing mixture. As the remediation solution continued
to move upward, the rising rate of the dielectric permittivity decreased and
asymptotically approached a constant value. This steady value was smaller
than that observed when the soil specimens were completely saturated by only
the flushing solution (i.e., wd, wda no. 1 or wda no. 2), which in our tests
corresponds to the condition of a completely decontaminated soil. This
difference in values is undoubtedly due to oil confined in soil pores (i.e., NAPL residual saturation). For the same reason, residual saturation may
explain why insignificant oil remediation was observed for
<inline-formula><mml:math id="M108" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> values less than 0.15. This aspect may be explained by the fact that for low volumetric NAPL contents, the non-wetting fluid (oil) is disconnectedly distributed (i.e., immobile) in the soil samples, which means that <inline-formula><mml:math id="M109" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> is close to the limiting “residual value”,<?pagebreak page3598?> and thus NAPL loses its ability to move in the soil in response to a hydraulic gradient (i.e., capillary retention forces are greater than gravitational forces, which tend to immobilize the NAPL; Brost and DeVaull, 2000).</p>
      <p id="d1e2264">The NAPL volumes removed for different washing solutions and the initial
volumetric content of NAPL are compared in Fig. 3. For all the three
cleaning solutions adopted, the experiments ultimately demonstrate (for a
fixed <inline-formula><mml:math id="M110" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) the same results in terms of soil decontamination, and they show that NAPL removal increases with increasing <inline-formula><mml:math id="M111" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>. In some cases (i.e., <inline-formula><mml:math id="M112" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.15</mml:mn></mml:mrow></mml:math></inline-formula>, 0.20 and 0.30), contaminated samples flushed with the wda no. 1 solution yield slightly higher removal efficiency values compared to the samples flushed with wd and wda no. 2. Martel et al. (1998) suggest the need to investigate the best water–surfactant–alcohol combination in order to enhance NAPL solubilization in soil.</p>
</sec>
<sec id="Ch1.S5.SS2">
  <label>5.2</label><title>Model calibration and validation</title>
      <p id="d1e2312">For the model (Eq. 10) calibration methodology, with reference to the three washing solutions (wd, wda no. 1 and wda no. 2), we analyze the effect of the measured dielectric permittivity on volumetric NAPL content (<inline-formula><mml:math id="M113" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) in order to estimate the <inline-formula><mml:math id="M114" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> parameter of the model. The complete calibration dataset of estimated <inline-formula><mml:math id="M115" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> parameters is reported in Table 1. The <inline-formula><mml:math id="M116" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> parameter of the mixing model was determined, for a fixed <inline-formula><mml:math id="M117" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> value and washing solution, by an optimization procedure based on the least-squares technique and was kept constant for each of the remediation tests developed.</p>

<?xmltex \floatpos{t}?><table-wrap id="Ch1.T1"><?xmltex \currentcnt{1}?><label>Table 1</label><caption><p id="d1e2361">Estimated <inline-formula><mml:math id="M118" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> parameter of Eq. (10) for all three washing
solutions (wd, wda no. 1 and wda no. 2) and the volumetric NAPL content (<inline-formula><mml:math id="M119" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) tested.</p></caption><oasis:table frame="topbot"><?xmltex \begin{scaleboxenv}{.92}[.92]?><oasis:tgroup cols="8">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="left" colsep="1"/>
     <oasis:colspec colnum="3" colname="col3" align="center"/>
     <oasis:colspec colnum="4" colname="col4" align="center"/>
     <oasis:colspec colnum="5" colname="col5" align="center"/>
     <oasis:colspec colnum="6" colname="col6" align="center"/>
     <oasis:colspec colnum="7" colname="col7" align="center"/>
     <oasis:colspec colnum="8" colname="col8" align="center"/>
     <oasis:thead>
       <oasis:row>

         <oasis:entry colname="col1">Parameter</oasis:entry>

         <oasis:entry colname="col2">Washing</oasis:entry>

         <oasis:entry namest="col3" nameend="col8"><inline-formula><mml:math id="M120" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula></oasis:entry>

       </oasis:row>
       <oasis:row rowsep="1">

         <oasis:entry colname="col1"/>

         <oasis:entry colname="col2">solution</oasis:entry>

         <oasis:entry colname="col3"/>

         <oasis:entry colname="col4"/>

         <oasis:entry colname="col5"/>

         <oasis:entry colname="col6"/>

         <oasis:entry colname="col7"/>

         <oasis:entry colname="col8"/>

       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>

         <oasis:entry colname="col1" morerows="3"><inline-formula><mml:math id="M121" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula></oasis:entry>

         <oasis:entry colname="col2"/>

         <oasis:entry colname="col3">0.15</oasis:entry>

         <oasis:entry colname="col4">0.20</oasis:entry>

         <oasis:entry colname="col5">0.25</oasis:entry>

         <oasis:entry colname="col6">0.30</oasis:entry>

         <oasis:entry colname="col7">0.35</oasis:entry>

         <oasis:entry colname="col8">0.40</oasis:entry>

       </oasis:row>
       <oasis:row>

         <oasis:entry colname="col2">wd</oasis:entry>

         <oasis:entry colname="col3">0.45</oasis:entry>

         <oasis:entry colname="col4">0.30</oasis:entry>

         <oasis:entry colname="col5">0.49</oasis:entry>

         <oasis:entry colname="col6">0.65</oasis:entry>

         <oasis:entry colname="col7">0.67</oasis:entry>

         <oasis:entry colname="col8">0.55</oasis:entry>

       </oasis:row>
       <oasis:row>

         <oasis:entry colname="col2">wda no. 1</oasis:entry>

         <oasis:entry colname="col3">0.25</oasis:entry>

         <oasis:entry colname="col4">0.45</oasis:entry>

         <oasis:entry colname="col5">0.45</oasis:entry>

         <oasis:entry colname="col6">0.42</oasis:entry>

         <oasis:entry colname="col7">0.50</oasis:entry>

         <oasis:entry colname="col8">0.55</oasis:entry>

       </oasis:row>
       <oasis:row>

         <oasis:entry colname="col2">wda no. 2</oasis:entry>

         <oasis:entry colname="col3">0.20</oasis:entry>

         <oasis:entry colname="col4">0.45</oasis:entry>

         <oasis:entry colname="col5">0.30</oasis:entry>

         <oasis:entry colname="col6">0.45</oasis:entry>

         <oasis:entry colname="col7">0.55</oasis:entry>

         <oasis:entry colname="col8">0.52</oasis:entry>

       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup><?xmltex \end{scaleboxenv}?></oasis:table></table-wrap>

      <p id="d1e2559">A permittivity value of 3.70 was adopted for the solid phase. This value was
determined using the “immersion method”, which is commonly employed for estimating the <inline-formula><mml:math id="M122" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mi mathvariant="normal">s</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> of soils (Robinson and Friedman, 2003; Kameyama and Miyamoto, 2008; Comegna et al., 2013b).</p>
      <p id="d1e2574">For the sake of brevity, a selection of the experimental <inline-formula><mml:math id="M123" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> relationships (validation dataset) is reported in Fig. 4a–f. The data in Fig. 4 (except for Fig. 4e and f) show that some of the model-simulated values tend to overestimate the measured data. This behavior is mostly restricted to the beginning of the remediation process, when a rapid change in dielectric permittivity may be observed. This behavior was also verified in other tests (not shown here) and may be explained by invoking both NAPL properties such as liquid density; surface tension and viscosity; and soil properties including moisture content, relative permeability and soil porosity (Brost and DeVaull, 2000; Wang et al., 2013).</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F4" specific-use="star"><?xmltex \currentcnt{4}?><label>Figure 4</label><caption><p id="d1e2606"><bold>(a–f)</bold> Selection of observed (symbols) and modeled (dashed lines) volumetric NAPL content (<inline-formula><mml:math id="M124" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) versus dielectric permittivity (<inline-formula><mml:math id="M125" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi mathvariant="normal">s</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">ws</mml:mi><mml:mtext>-</mml:mtext><mml:mi mathvariant="normal">NAPL</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:math></inline-formula>), with reference to the three washing solutions (wd, wda no. 1 and wda no. 2) used during the remediation tests.</p></caption>
          <?xmltex \igopts{width=497.923228pt}?><graphic xlink:href="https://hess.copernicus.org/articles/23/3593/2019/hess-23-3593-2019-f04.png"/>

        </fig>

      <p id="d1e2648">Mercer and Cohen (1990) referred to the existence, in NAPL-contaminated
soils, of a “double fluid domain”, defined as the composition of the following: (i) mobile pools, which are NAPL-connected phases that move in the soil, and (ii) immobile residuals (i.e., low-permeability regions), which depend on small disconnected blobs or ganglia within<?pagebreak page3599?> the contaminated soil (see also Sect. 5.1 above). As long as the flushing continues, mobile pools are reduced and the oil tends increasingly to be trapped in the immobile areas. This means that, during soil cleaning, the capacity of non-wetting fluids to respond to gravitational forces gradually diminishes (Luckner et al., 1989). From a dielectric point of view, this mechanism may appear as a rapid dielectric permittivity increase (identified in Fig. 4 as a “fast oil mobility region”) within a few pore volumes. When this fast mobility mechanism is dominant, the predictions of Eq. (10) fail.</p>
      <p id="d1e2651">Another possible explanation for this discrepancy between the observed and
the predicted permittivity values may be linked to the propensity of
NAPL-water mixtures to form macroinclusions in the soil (Persson and
Berndtsson, 2002), which affected the initial pore-scale distribution of
NAPL, and thus the global dielectric response of the medium (Ferré et
al., 1996), during the first remediation stages.</p>
      <p id="d1e2654">However, since the phenomenon is mostly limited to the initial part of the
washing process, overall model effectiveness is not compromised, as also
shown in Table 2, which summarizes the goodness-of-fit statistical indices,
and in Fig. 5a–f, where the estimated <inline-formula><mml:math id="M126" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> from Eq. (10) and the known <inline-formula><mml:math id="M127" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> are illustrated in a series of <inline-formula><mml:math id="M128" display="inline"><mml:mrow><mml:mn mathvariant="normal">1</mml:mn><mml:mo>:</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula> scatter plots.</p>

      <?xmltex \floatpos{t}?><fig id="Ch1.F5" specific-use="star"><?xmltex \currentcnt{5}?><label>Figure 5</label><caption><p id="d1e2694"><bold>(a–f)</bold> Measured (Eq. 10) vs. known volumetric NAPL content (<inline-formula><mml:math id="M129" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) of contaminated soils, with reference to the
different remediation tests in Fig. 4.</p></caption>
          <?xmltex \igopts{width=341.433071pt}?><graphic xlink:href="https://hess.copernicus.org/articles/23/3593/2019/hess-23-3593-2019-f05.png"/>

        </fig>

<?xmltex \floatpos{t}?><table-wrap id="Ch1.T2"><?xmltex \currentcnt{2}?><label>Table 2</label><caption><p id="d1e2719">Model efficiency (EF) and mean bias error (MBE) statistical indices,
referring to measured and predicted (Eq. 10) volumetric NAPL content (<inline-formula><mml:math id="M130" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>), calculated in the range of model applicability (<inline-formula><mml:math id="M131" display="inline"><mml:mrow><mml:mn mathvariant="normal">0.15</mml:mn><mml:mo>≤</mml:mo><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>≤</mml:mo><mml:mn mathvariant="normal">0.40</mml:mn></mml:mrow></mml:math></inline-formula>).</p></caption><oasis:table frame="topbot"><?xmltex \begin{scaleboxenv}{.89}[.89]?><oasis:tgroup cols="9">
     <oasis:colspec colnum="1" colname="col1" align="left"/>
     <oasis:colspec colnum="2" colname="col2" align="center"/>
     <oasis:colspec colnum="3" colname="col3" align="right"/>
     <oasis:colspec colnum="4" colname="col4" align="left"/>
     <oasis:colspec colnum="5" colname="col5" align="center"/>
     <oasis:colspec colnum="6" colname="col6" align="right"/>
     <oasis:colspec colnum="7" colname="col7" align="left"/>
     <oasis:colspec colnum="8" colname="col8" align="center"/>
     <oasis:colspec colnum="9" colname="col9" align="right"/>
     <oasis:thead>
       <oasis:row>
         <oasis:entry colname="col1">Washing</oasis:entry>
         <oasis:entry rowsep="1" namest="col2" nameend="col3"><inline-formula><mml:math id="M132" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.15</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry rowsep="1" namest="col5" nameend="col6"><inline-formula><mml:math id="M133" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.20</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col7"/>
         <oasis:entry rowsep="1" namest="col8" nameend="col9"><inline-formula><mml:math id="M134" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.25</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">solution</oasis:entry>
         <oasis:entry colname="col2">EF</oasis:entry>
         <oasis:entry colname="col3">MBE</oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">EF</oasis:entry>
         <oasis:entry colname="col6">MBE</oasis:entry>
         <oasis:entry colname="col7"/>
         <oasis:entry colname="col8">EF</oasis:entry>
         <oasis:entry colname="col9">MBE</oasis:entry>
       </oasis:row>
     </oasis:thead>
     <oasis:tbody>
       <oasis:row>
         <oasis:entry colname="col1">wd</oasis:entry>
         <oasis:entry colname="col2">0.98</oasis:entry>
         <oasis:entry colname="col3">1.548</oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">0.93</oasis:entry>
         <oasis:entry colname="col6"><inline-formula><mml:math id="M135" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.422</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col7"/>
         <oasis:entry colname="col8">0.96</oasis:entry>
         <oasis:entry colname="col9">0.570</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">wda no. 1</oasis:entry>
         <oasis:entry colname="col2">0.86</oasis:entry>
         <oasis:entry colname="col3">0.405</oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">0.99</oasis:entry>
         <oasis:entry colname="col6">0.516</oasis:entry>
         <oasis:entry colname="col7"/>
         <oasis:entry colname="col8">0.97</oasis:entry>
         <oasis:entry colname="col9"><inline-formula><mml:math id="M136" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.048</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">wda no. 2</oasis:entry>
         <oasis:entry colname="col2">0.84</oasis:entry>
         <oasis:entry colname="col3">0.148</oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">0.94</oasis:entry>
         <oasis:entry colname="col6">0.420</oasis:entry>
         <oasis:entry colname="col7"/>
         <oasis:entry colname="col8">0.66</oasis:entry>
         <oasis:entry colname="col9">0001</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">Washing</oasis:entry>
         <oasis:entry rowsep="1" namest="col2" nameend="col3"><inline-formula><mml:math id="M137" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.30</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry rowsep="1" namest="col5" nameend="col6"><inline-formula><mml:math id="M138" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.35</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col7"/>
         <oasis:entry rowsep="1" namest="col8" nameend="col9"><inline-formula><mml:math id="M139" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn mathvariant="normal">0.40</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row rowsep="1">
         <oasis:entry colname="col1">solution</oasis:entry>
         <oasis:entry colname="col2">EF</oasis:entry>
         <oasis:entry colname="col3">MBE</oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">EF</oasis:entry>
         <oasis:entry colname="col6">MBE</oasis:entry>
         <oasis:entry colname="col7"/>
         <oasis:entry colname="col8">EF</oasis:entry>
         <oasis:entry colname="col9">MBE</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">wd</oasis:entry>
         <oasis:entry colname="col2">0.98</oasis:entry>
         <oasis:entry colname="col3"><inline-formula><mml:math id="M140" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.023</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">0.99</oasis:entry>
         <oasis:entry colname="col6"><inline-formula><mml:math id="M141" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.153</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col7"/>
         <oasis:entry colname="col8">0.99</oasis:entry>
         <oasis:entry colname="col9"><inline-formula><mml:math id="M142" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.179</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">wda no. 1</oasis:entry>
         <oasis:entry colname="col2">0.95</oasis:entry>
         <oasis:entry colname="col3"><inline-formula><mml:math id="M143" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.074</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">0.99</oasis:entry>
         <oasis:entry colname="col6"><inline-formula><mml:math id="M144" display="inline"><mml:mrow><mml:mo>-</mml:mo><mml:mn mathvariant="normal">0.066</mml:mn></mml:mrow></mml:math></inline-formula></oasis:entry>
         <oasis:entry colname="col7"/>
         <oasis:entry colname="col8">0.99</oasis:entry>
         <oasis:entry colname="col9">0.303</oasis:entry>
       </oasis:row>
       <oasis:row>
         <oasis:entry colname="col1">wda no. 2</oasis:entry>
         <oasis:entry colname="col2">0.91</oasis:entry>
         <oasis:entry colname="col3">0.014</oasis:entry>
         <oasis:entry colname="col4"/>
         <oasis:entry colname="col5">0.97</oasis:entry>
         <oasis:entry colname="col6">0.326</oasis:entry>
         <oasis:entry colname="col7"/>
         <oasis:entry colname="col8">0.99</oasis:entry>
         <oasis:entry colname="col9">0.019</oasis:entry>
       </oasis:row>
     </oasis:tbody>
   </oasis:tgroup><?xmltex \end{scaleboxenv}?></oasis:table></table-wrap>

      <?pagebreak page3600?><p id="d1e3190">Overall, both graphical and quantitative evaluations in terms of MBE and EF
reveal the suitability of the dielectric model adopted to estimate the
volumetric NAPL content in the <inline-formula><mml:math id="M145" display="inline"><mml:mrow><mml:msub><mml:mi mathvariant="italic">θ</mml:mi><mml:mi mathvariant="normal">NAPL</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula> range 0.15–0.40.</p>
</sec>
</sec>
<sec id="Ch1.S6" sec-type="conclusions">
  <label>6</label><title>Conclusions</title>
      <p id="d1e3213">This paper presented an extensive dataset of remediation experiments that
were conducted on a laboratory scale using corn oil as a soil contaminant
and three different solutions for soil cleaning. The results of these tests
were employed to investigate the potential of the TDR technique in
monitoring the development of a steady-state decontamination process.</p>
      <p id="d1e3216">Dielectric data analysis showed that, during soil flushing, dielectric
permittivity behavior is highly dependent on the initial volumetric content
and intrinsic permittivity of the specific NAPL: “removal of NAPL produces an increase in bulk dielectric permittivity” due to the low value of oil permittivity. The experiments conducted also allowed us to calibrate and validate a dielectric mixing model (Eq. 10). The model outcomes are
encouraging; the calculated statistical<?pagebreak page3601?> indices confirmed a high accuracy in
NAPL predictions of the <inline-formula><mml:math id="M146" display="inline"><mml:mi mathvariant="italic">α</mml:mi></mml:math></inline-formula> model at different stages during soil
cleaning, with the only exception of the initial cleaning stage
(confined to the low values of <inline-formula><mml:math id="M147" display="inline"><mml:mi>T</mml:mi></mml:math></inline-formula>), where the eventual presence of a “fast flow region” may limit its applicability.</p>
      <p id="d1e3233">The approach requires additional experiments and datasets for model
calibration and validation in different pedological contexts, mainly to
confirm the potential of the methodology developed. Furthermore, an effort
should be made, introducing the water phase, ab initio in the experimental setup, to simulate a possible natural contamination–decontamination scenario more accurately. Finally, full field-scale tests should also be conducted to evaluate the performance of Eq. (10) in real field conditions.</p>
</sec>

      
      </body>
    <back><notes notes-type="dataavailability"><title>Data availability</title>

      <p id="d1e3240">The dataset used in this paper is available on request to
alessandro.comegna@unibas.it.</p>
  </notes><notes notes-type="authorcontribution"><title>Author contributions</title>

      <p id="d1e3246">The authors jointly collaborated on the development of the present research, which is based on ACom's idea.</p>
  </notes><notes notes-type="competinginterests"><title>Competing interests</title>

      <p id="d1e3252">The authors declare that they have no conflict of interest.</p>
  </notes><notes notes-type="reviewstatement"><title>Review statement</title>

      <p id="d1e3258">This paper was edited by Marnik Vanclooster and reviewed by Ty P. A. Ferré, Magnus Persson, and one anonymous referee.</p>
  </notes><?xmltex \hack{\newpage}?><ref-list>
    <title>References</title>

      <ref id="bib1.bib1"><label>1</label><?label 1?><mixed-citation>
Allison, L. E.: Organic carbon, in: Methods of Soil Analysis, Part 1, Agron. Monograph, vol. 9, edited by: Klute, A., ASA and SSSA, Madison, 1367–1378, 1965.</mixed-citation></ref>
      <ref id="bib1.bib2"><label>2</label><?label 1?><mixed-citation>
Birchack, J. R., Gardner, C. Z. G., Hipp, J. E., and Victor, J. M.: High
dielectric constant microwave probes for sensing soil moisture, Proc. IEEEE,
62, 93–98, 1974.</mixed-citation></ref>
      <ref id="bib1.bib3"><label>3</label><?label 1?><mixed-citation>
Brost, E. J. and DeVaull, G. E.: Non-Aqueous phase liquid (NAPL) mobility
limits in soils, Soil Groundw. Res. Bull., 9, 1–9, 2000.</mixed-citation></ref>
      <ref id="bib1.bib4"><label>4</label><?label 1?><mixed-citation>
Comegna, A., Coppola, A., Dragonetti, G., and Sommella, A.: Dielectric response of a variable saturated soil contaminated by Non-Aqueous Phase
Liquids (NAPLs), Proced. Environ. Sci., 19, 701–710, 2013a.</mixed-citation></ref>
      <ref id="bib1.bib5"><label>5</label><?label 1?><mixed-citation>
Comegna, A., Coppola, A., Dragonetti, G., Severino, G., Sommella, A., and Basile, A.: Dielectric properties of a tilled sandy volcanic-vesuvian soil with moderate andic features, Soil Till. Res., 133, 93–100, 2013b.</mixed-citation></ref>
      <ref id="bib1.bib6"><label>6</label><?label 1?><mixed-citation>Comegna, A., Coppola, A., Dragonetti, G., Chaali, N., and Sommella, A.: Time
domain reflectometry-measuring dielectric permittivity to detect soil
non-aqueous phase liquids contamination-decontamination processes, J.
Agricult. Eng., XLIV, e167, <ext-link xlink:href="https://doi.org/10.4081/jae.2013.409" ext-link-type="DOI">10.4081/jae.2013.409</ext-link>, 2013c.</mixed-citation></ref>
      <ref id="bib1.bib7"><label>7</label><?label 1?><mixed-citation>Comegna, A., Coppola, A., Dragonetti, G., and Sommella, A.: Estimating
non-aqueous phase liquid (NAPL) content in variable saturated soils using
time domain reflectometry (TDR), Vadose Zone J., 15, 13, <ext-link xlink:href="https://doi.org/10.2136/vzj2015.11.0145" ext-link-type="DOI">10.2136/vzj2015.11.0145</ext-link>, 2016.</mixed-citation></ref>
      <ref id="bib1.bib8"><label>8</label><?label 1?><mixed-citation>Comegna, A., Coppola, A., Dragonetti, G., and Sommella, A.: Interpreting TDR
signal propagation through soils with distinct layers of nonaqueous-phase
liquid and water content, Vadose Zone J., 16, 5, <ext-link xlink:href="https://doi.org/10.2136/vzj2017.07.0141" ext-link-type="DOI">10.2136/vzj2017.07.0141</ext-link>, 2017.</mixed-citation></ref>
      <ref id="bib1.bib9"><label>9</label><?label 1?><mixed-citation>
Coppola, A., Dragonetti, G., Comegna, A., Lamaddalena, N., Caushi, B., Haikal, M. A., and Basile, A.: Measuring and modeling water content in stony
soils, Soil Till. Res., 128, 9–22, 2013.</mixed-citation></ref>
      <ref id="bib1.bib10"><label>10</label><?label 1?><mixed-citation>Coppola, A., Comegna, A., Dragonetti, G., Gerke, H. H., and Basile, A.:
Simulated water flow and solute transport in shrinking soils, Vadose Zone
J., 14, 9, <ext-link xlink:href="https://doi.org/10.2136/vzj2015.02.0021" ext-link-type="DOI">10.2136/vzj2015.02.0021</ext-link>, 2015.</mixed-citation></ref>
      <ref id="bib1.bib11"><label>11</label><?label 1?><mixed-citation>
Day, P. R.: Particle fractionation and particle-size analysis, in: Methods of Soil Analysis, Part 1, edited by: Black, C. A., American Society of Agronomy, Madison, 545–567, 1965.</mixed-citation></ref>
      <ref id="bib1.bib12"><label>12</label><?label 1?><mixed-citation>
Eckert, D. J.: Soil pH, in: Recommended chemical soil test procedures for the North Central Region, edited by: Dahnke, W. C., North Bulletin No. 221, Dakota Agricultural Experiment Station, Fargo, 6–8, 1988.</mixed-citation></ref>
      <ref id="bib1.bib13"><label>13</label><?label 1?><mixed-citation>
Endres, A. and Redman, J. D.: Modeling the electrical properties of porous
rocks and soils containing immiscible contaminants, edited by: Bell, R. S. and Lepper, C. M., in: Proceedings of the Symposium on the application of geophysics to Engineering and Environmental Problems, San Diego, 21–38, 1993.</mixed-citation></ref>
      <ref id="bib1.bib14"><label>14</label><?label 1?><mixed-citation>
Ferré, P. A. and Topp, G. C.: Time domain reflectometry, in: Methods of
Soil Analysis: Part 4, Physical Methods, edited by: Dane, J. H. and Topp, G. C., SSSA, Madison, WI, 434–446, 2002.</mixed-citation></ref>
      <ref id="bib1.bib15"><label>15</label><?label 1?><mixed-citation>
Ferré, P. A., Rudolph, D. L., and Kachanoski, R. G.: Spatial averaging of
water content by time domain reflectometry: Implications for twin rod probes
with and without dielectric coatings, Water Resour. Res., 32, 271–279, 1996.</mixed-citation></ref>
      <ref id="bib1.bib16"><label>16</label><?label 1?><mixed-citation>
Fetter, C. W.: Contaminant hydrogeology, in: Practice, 2nd Edn., Hall Englewood Cliffs, NJ, 247–315, 1993.</mixed-citation></ref>
      <?pagebreak page3602?><ref id="bib1.bib17"><label>17</label><?label 1?><mixed-citation>
Fitts, C. R.: Contamination sources, in: Groundwater Science, Academic Press,
London, 339–343, 2002.</mixed-citation></ref>
      <ref id="bib1.bib18"><label>18</label><?label 1?><mixed-citation>
Francisca, M. and Montoro, M. A.: Measuring the dielectric properties of
soil-organic mixtures using coaxial impedance dielectric reflectometry, J.
Appl. Geophys., 80, 101–109, 2012.</mixed-citation></ref>
      <ref id="bib1.bib19"><label>19</label><?label 1?><mixed-citation>
Freeze, R. A.: The unpleasant truths about waste management, in: The
Environmental Pendulum, University of California Press, Berkeley, 147–196, 2000.</mixed-citation></ref>
      <ref id="bib1.bib20"><label>20</label><?label 1?><mixed-citation>
Heimovaara, T. J.: Frequency domain analysis of time domain reflectometry
waveforms, 2, A four-component complex dielectric mixing model for soil,
Water Resour. Res., 30, 201–209, 1994.</mixed-citation></ref>
      <ref id="bib1.bib21"><label>21</label><?label 1?><mixed-citation>
Hilhorst, M. A.: Dielectric characterisation of soil, PhD dissertation,
Wageningen Agricultural University, Wageningen, the Netherlands, 1998.</mixed-citation></ref>
      <ref id="bib1.bib22"><label>22</label><?label 1?><mixed-citation>
IUSS Working Group WRB: World reference base for soil resources 2006: A
framework for international classification, correlation and communication,
2nd Edn., World Soil Resour. Rep. 103, FAO, Rome, 2006.</mixed-citation></ref>
      <ref id="bib1.bib23"><label>23</label><?label 1?><mixed-citation>
Kameyama, K. and Miyamoto, T.: Measurement of solid phase permittivity for
soils by time domain reflectometry, Eur. J. Soil Sci., 59, 1253–1259, 2008.</mixed-citation></ref>
      <ref id="bib1.bib24"><label>24</label><?label 1?><mixed-citation>
Knight, R.: Ground penetrating radar for environmental applications, Annu.
Rev. Earth Planet. Sci., 29, 229–255, 2001.</mixed-citation></ref>
      <ref id="bib1.bib25"><label>25</label><?label 1?><mixed-citation>
Knight, R. and Endres, A.: A new concept in modeling the dielectric response
of sandstones: Defining a wetted rock and bulk water systems, Geophysics, 55, 586–594, 1990.</mixed-citation></ref>
      <ref id="bib1.bib26"><label>26</label><?label 1?><mixed-citation>
Legates, D. R. and McCabe Jr., G. J.: Evaluating the use of “goodness-of-fit” measures in hydrologic and hydroclimatic model validation, Water Resour. Res., 35, 233–241, 1999.</mixed-citation></ref>
      <ref id="bib1.bib27"><label>27</label><?label 1?><mixed-citation>
Luckner, L., van Genuchten, M. T., and Nielsen, D. R.: A consistent set of
parametric models for the two-phase flow of immiscible fluids in the subsurface, Water Resour. Res., 25, 2187–2193, 1989.</mixed-citation></ref>
      <ref id="bib1.bib28"><label>28</label><?label 1?><mixed-citation>
Martel, R., Gelinas, P. J., and Saumure, L.: Aquifer washing by micellar
solutions: 3 Field test at the Thounin Sand Pit (L'Assomption Quebec,
Canada), J. Contam. Hydrol., 30, 33–48, 1998.</mixed-citation></ref>
      <ref id="bib1.bib29"><label>29</label><?label 1?><mixed-citation>
Mercer, J. W. and Cohen, R. M.: A review of immiscible fluids in the subsurface: Properties, models, characterization, and remediation, J.
Contam. Hydrol., 6, 107–163, 1990.</mixed-citation></ref>
      <ref id="bib1.bib30"><label>30</label><?label 1?><mixed-citation>
Miller, J. J. and Curtis, D.: Electrical conductivity and soluble ions, in:
Soil Sampling and Methods of Analysis, edited by: Carter, M. R. and  Gregorich, E. G., Canadian Society of Soil Science, Broken Sound Pkwy, NW, 2006.</mixed-citation></ref>
      <ref id="bib1.bib31"><label>31</label><?label 1?><mixed-citation>
Mohamed, A.M.O., and Said R.A.: Detection of organic pollutants in sandy
soils via TDR and eigendecomposition. J. Contain. Hydrol., 76, 235–249,
2005.</mixed-citation></ref>
      <ref id="bib1.bib32"><label>32</label><?label 1?><mixed-citation>
Moroizumi, T. and Sasaki, Y.: Estimating the nonaqueous-phase liquid content in saturated sandy soil using amplitude domain reflectometry, Soil Sci. Soc. Am. J., 72, 1520–1526, 2006.</mixed-citation></ref>
      <ref id="bib1.bib33"><label>33</label><?label 1?><mixed-citation>Parnian, M. M. and Ayatollahi, S.: Surfactant remediation of LNAPL contaminated soil; effects of adding alkaline and foam producing substances,
Iran. J. Chem. Eng., 2, 34–44, 2008.
 </mixed-citation></ref><?xmltex \hack{\newpage}?>
      <ref id="bib1.bib34"><label>34</label><?label 1?><mixed-citation>
Pennell, K. D., Jin, M., Abriola, L. M., and Pope, G. A.: Surfactant enhanced
remediation of soil columns contaminated by residual tetrachloroethylene, J.
Contam. Hydrol., 16, 35–53, 1994.</mixed-citation></ref>
      <ref id="bib1.bib35"><label>35</label><?label 1?><mixed-citation>Persson, M. and Berndtsson, R.: Measuring nonaqueous phase liquid saturation in soil using time domain reflectometry, Water Resour. Res., 38, 1064, <ext-link xlink:href="https://doi.org/10.1029/2001WR000523" ext-link-type="DOI">10.1029/2001WR000523</ext-link>, 2002.</mixed-citation></ref>
      <ref id="bib1.bib36"><label>36</label><?label 1?><mixed-citation>
Redman, J. D., Kueper, B. H., and Annan, A. P.: Dielectric stratigraphy of a
DNAPL spill and implications for detection with ground penetrating radar:
Aquifer restoration, ground water monitoring and geophysical methods, in: 5th national outdoor action conference, 13–16 May 1991, Las Vegas, Nevada, 1017–1030, 1991.</mixed-citation></ref>
      <ref id="bib1.bib37"><label>37</label><?label 1?><mixed-citation>
Redman, J. D. and DeRyck, S. M.: Monitoring Non-Aqueous Phase Liquids in the
Subsurface with Multilevel Time Domain Reflectometry Probes, in: Proceedings of the Symposium on Time Domain Reflectometry in Environmental, Infrastructure, and Mining Applications, 7–9 September 1994, Special Publication SP 19-94, NTIS PB95-105789, US Bureau of Mines, Evanston, Illinois, 207–215, 1994.</mixed-citation></ref>
      <ref id="bib1.bib38"><label>38</label><?label 1?><mixed-citation>
Rinaldi, V. A. and Francisca, F. M.: Removal of immiscible contaminants from
sandy soils monitored by means of dielectric measurements, J. Environ. Eng.,
132, 931–939, 2006.</mixed-citation></ref>
      <ref id="bib1.bib39"><label>39</label><?label 1?><mixed-citation>
Robinson, D. A. and Friedman, S. P.: The effective permittivity of dense packing of glass beads, quartz sand and their mixtures immersed in different
dielectric backgrounds, J. Non-Crystal. Solids, 305, 261–267, 2003.</mixed-citation></ref>
      <ref id="bib1.bib40"><label>40</label><?label 1?><mixed-citation>
Robinson, D. A., Jones, S. B., Wraith, J. M., and Or, D.: A review of advances in dielectric and electric conductivity measurements using time domain reflectometry, Vadose Zone J., 2, 444–475, 2003.</mixed-citation></ref>
      <ref id="bib1.bib41"><label>41</label><?label 1?><mixed-citation>
Roth, K., Schulin, R., Fluhler, H., and Attinger, W.: Calibration of time
domain reflectometry for water content measurements using a composite
dielectric approach, Water Resour. Res., 26, 2267–2273, 1990.</mixed-citation></ref>
      <ref id="bib1.bib42"><label>42</label><?label 1?><mixed-citation>
Topp, C. G., Davis, J. L., and Annan, A. P.: Electromagnetic determination of
soil water content: Measurements in coaxial transmission lines, Water Resour. Res., 16, 574–582, 1980.</mixed-citation></ref>
      <ref id="bib1.bib43"><label>43</label><?label 1?><mixed-citation>
Truong Hong, T. and Bettahar, M.: Effect of the water-oil ratio on
brine/surfactant/alcohol/oil systems optimized for soil remediation, J. Mater Cy. Waste Manage., 2, 109–117, 2000.</mixed-citation></ref>
      <ref id="bib1.bib44"><label>44</label><?label 1?><mixed-citation>
van Dam, R. L., Borchers, B., and Hendrickx, J. M. H.: Methods for prediction of soil dielectric properties: a review, in: Proceedings of SPIE-the International Society for Optical Engineering, 28 March–1 April 2005, Orlando, Florida, 2005.</mixed-citation></ref>
      <ref id="bib1.bib45"><label>45</label><?label 1?><mixed-citation>
Wang, X. P., Quan, G. J., Pan, Y. X., Hu, R., Zhang, Y. F, Tedeschi, A., Basile, A., Comegna, A., Coppola, A., and de Mascellis, R.: Comparison of hydraulic behaviour of unvegetated and vegetation-stabilized sand dunes in arid desert ecosystems, Ecohydrology, 62, 264–274, 2013.</mixed-citation></ref>

  </ref-list></back>
    <!--<article-title-html>A soil non-aqueous phase liquid (NAPL) flushing laboratory experiment based on measuring the dielectric properties of soil–organic mixtures via time domain reflectometry (TDR)</article-title-html>
<abstract-html><p>The term non-aqueous phase liquid (NAPL) refers to a group of organic compounds with scarce solubility in water. They are the products of various human activities and may be accidentally introduced into the soil system. Given their toxicity level and high mobility, NAPLs constitute a serious geo-environmental problem. Contaminant distribution in the soil and groundwater contains fundamental information for the remediation of polluted
soil sites. The present research explored the possible employment of time
domain reflectometry (TDR) to estimate pollutant removal in a silt-loam soil that was primarily contaminated with a corn oil as a light NAPL and then flushed with different washing solutions. Known mixtures of soil and NAPL were prepared in the laboratory to achieve soil specimens with varying
pollution levels. The prepared soil samples were repacked into plastic
cylinders and then placed in testing cells. Washing solutions were then
injected upward into the contaminated sample, and both the quantity of
remediated NAPL and the bulk dielectric permittivity of the soil sample were
determined. The above data were also used to calibrate and validate a
dielectric model (the <i>α</i> mixing model) which permits the volumetric
NAPL content (<i>θ</i><sub>NAPL</sub>; m<sup>3</sup>&thinsp;m<sup>−3</sup>) within the contaminated sample to be determined and quantified during the different decontamination stages. Our results demonstrate that during a decontamination process, the TDR device is NAPL-sensitive: the dielectric permittivity of the medium increases as the NAPL volume decreases. Moreover, decontamination progression can be monitored using a simple (one-parameter) mixing model.</p></abstract-html>
<ref-html id="bib1.bib1"><label>1</label><mixed-citation>
Allison, L. E.: Organic carbon, in: Methods of Soil Analysis, Part 1, Agron. Monograph, vol. 9, edited by: Klute, A., ASA and SSSA, Madison, 1367–1378, 1965.
</mixed-citation></ref-html>
<ref-html id="bib1.bib2"><label>2</label><mixed-citation>
Birchack, J. R., Gardner, C. Z. G., Hipp, J. E., and Victor, J. M.: High
dielectric constant microwave probes for sensing soil moisture, Proc. IEEEE,
62, 93–98, 1974.
</mixed-citation></ref-html>
<ref-html id="bib1.bib3"><label>3</label><mixed-citation>
Brost, E. J. and DeVaull, G. E.: Non-Aqueous phase liquid (NAPL) mobility
limits in soils, Soil Groundw. Res. Bull., 9, 1–9, 2000.
</mixed-citation></ref-html>
<ref-html id="bib1.bib4"><label>4</label><mixed-citation>
Comegna, A., Coppola, A., Dragonetti, G., and Sommella, A.: Dielectric response of a variable saturated soil contaminated by Non-Aqueous Phase
Liquids (NAPLs), Proced. Environ. Sci., 19, 701–710, 2013a.
</mixed-citation></ref-html>
<ref-html id="bib1.bib5"><label>5</label><mixed-citation>
Comegna, A., Coppola, A., Dragonetti, G., Severino, G., Sommella, A., and Basile, A.: Dielectric properties of a tilled sandy volcanic-vesuvian soil with moderate andic features, Soil Till. Res., 133, 93–100, 2013b.
</mixed-citation></ref-html>
<ref-html id="bib1.bib6"><label>6</label><mixed-citation>
Comegna, A., Coppola, A., Dragonetti, G., Chaali, N., and Sommella, A.: Time
domain reflectometry-measuring dielectric permittivity to detect soil
non-aqueous phase liquids contamination-decontamination processes, J.
Agricult. Eng., XLIV, e167, <a href="https://doi.org/10.4081/jae.2013.409" target="_blank">https://doi.org/10.4081/jae.2013.409</a>, 2013c.
</mixed-citation></ref-html>
<ref-html id="bib1.bib7"><label>7</label><mixed-citation>
Comegna, A., Coppola, A., Dragonetti, G., and Sommella, A.: Estimating
non-aqueous phase liquid (NAPL) content in variable saturated soils using
time domain reflectometry (TDR), Vadose Zone J., 15, 13, <a href="https://doi.org/10.2136/vzj2015.11.0145" target="_blank">https://doi.org/10.2136/vzj2015.11.0145</a>, 2016.
</mixed-citation></ref-html>
<ref-html id="bib1.bib8"><label>8</label><mixed-citation>
Comegna, A., Coppola, A., Dragonetti, G., and Sommella, A.: Interpreting TDR
signal propagation through soils with distinct layers of nonaqueous-phase
liquid and water content, Vadose Zone J., 16, 5, <a href="https://doi.org/10.2136/vzj2017.07.0141" target="_blank">https://doi.org/10.2136/vzj2017.07.0141</a>, 2017.
</mixed-citation></ref-html>
<ref-html id="bib1.bib9"><label>9</label><mixed-citation>
Coppola, A., Dragonetti, G., Comegna, A., Lamaddalena, N., Caushi, B., Haikal, M. A., and Basile, A.: Measuring and modeling water content in stony
soils, Soil Till. Res., 128, 9–22, 2013.
</mixed-citation></ref-html>
<ref-html id="bib1.bib10"><label>10</label><mixed-citation>
Coppola, A., Comegna, A., Dragonetti, G., Gerke, H. H., and Basile, A.:
Simulated water flow and solute transport in shrinking soils, Vadose Zone
J., 14, 9, <a href="https://doi.org/10.2136/vzj2015.02.0021" target="_blank">https://doi.org/10.2136/vzj2015.02.0021</a>, 2015.
</mixed-citation></ref-html>
<ref-html id="bib1.bib11"><label>11</label><mixed-citation>
Day, P. R.: Particle fractionation and particle-size analysis, in: Methods of Soil Analysis, Part 1, edited by: Black, C. A., American Society of Agronomy, Madison, 545–567, 1965.
</mixed-citation></ref-html>
<ref-html id="bib1.bib12"><label>12</label><mixed-citation>
Eckert, D. J.: Soil pH, in: Recommended chemical soil test procedures for the North Central Region, edited by: Dahnke, W. C., North Bulletin No. 221, Dakota Agricultural Experiment Station, Fargo, 6–8, 1988.
</mixed-citation></ref-html>
<ref-html id="bib1.bib13"><label>13</label><mixed-citation>
Endres, A. and Redman, J. D.: Modeling the electrical properties of porous
rocks and soils containing immiscible contaminants, edited by: Bell, R. S. and Lepper, C. M., in: Proceedings of the Symposium on the application of geophysics to Engineering and Environmental Problems, San Diego, 21–38, 1993.
</mixed-citation></ref-html>
<ref-html id="bib1.bib14"><label>14</label><mixed-citation>
Ferré, P. A. and Topp, G. C.: Time domain reflectometry, in: Methods of
Soil Analysis: Part 4, Physical Methods, edited by: Dane, J. H. and Topp, G. C., SSSA, Madison, WI, 434–446, 2002.
</mixed-citation></ref-html>
<ref-html id="bib1.bib15"><label>15</label><mixed-citation>
Ferré, P. A., Rudolph, D. L., and Kachanoski, R. G.: Spatial averaging of
water content by time domain reflectometry: Implications for twin rod probes
with and without dielectric coatings, Water Resour. Res., 32, 271–279, 1996.
</mixed-citation></ref-html>
<ref-html id="bib1.bib16"><label>16</label><mixed-citation>
Fetter, C. W.: Contaminant hydrogeology, in: Practice, 2nd Edn., Hall Englewood Cliffs, NJ, 247–315, 1993.
</mixed-citation></ref-html>
<ref-html id="bib1.bib17"><label>17</label><mixed-citation>
Fitts, C. R.: Contamination sources, in: Groundwater Science, Academic Press,
London, 339–343, 2002.
</mixed-citation></ref-html>
<ref-html id="bib1.bib18"><label>18</label><mixed-citation>
Francisca, M. and Montoro, M. A.: Measuring the dielectric properties of
soil-organic mixtures using coaxial impedance dielectric reflectometry, J.
Appl. Geophys., 80, 101–109, 2012.
</mixed-citation></ref-html>
<ref-html id="bib1.bib19"><label>19</label><mixed-citation>
Freeze, R. A.: The unpleasant truths about waste management, in: The
Environmental Pendulum, University of California Press, Berkeley, 147–196, 2000.
</mixed-citation></ref-html>
<ref-html id="bib1.bib20"><label>20</label><mixed-citation>
Heimovaara, T. J.: Frequency domain analysis of time domain reflectometry
waveforms, 2, A four-component complex dielectric mixing model for soil,
Water Resour. Res., 30, 201–209, 1994.
</mixed-citation></ref-html>
<ref-html id="bib1.bib21"><label>21</label><mixed-citation>
Hilhorst, M. A.: Dielectric characterisation of soil, PhD dissertation,
Wageningen Agricultural University, Wageningen, the Netherlands, 1998.
</mixed-citation></ref-html>
<ref-html id="bib1.bib22"><label>22</label><mixed-citation>
IUSS Working Group WRB: World reference base for soil resources 2006: A
framework for international classification, correlation and communication,
2nd Edn., World Soil Resour. Rep. 103, FAO, Rome, 2006.
</mixed-citation></ref-html>
<ref-html id="bib1.bib23"><label>23</label><mixed-citation>
Kameyama, K. and Miyamoto, T.: Measurement of solid phase permittivity for
soils by time domain reflectometry, Eur. J. Soil Sci., 59, 1253–1259, 2008.
</mixed-citation></ref-html>
<ref-html id="bib1.bib24"><label>24</label><mixed-citation>
Knight, R.: Ground penetrating radar for environmental applications, Annu.
Rev. Earth Planet. Sci., 29, 229–255, 2001.
</mixed-citation></ref-html>
<ref-html id="bib1.bib25"><label>25</label><mixed-citation>
Knight, R. and Endres, A.: A new concept in modeling the dielectric response
of sandstones: Defining a wetted rock and bulk water systems, Geophysics, 55, 586–594, 1990.
</mixed-citation></ref-html>
<ref-html id="bib1.bib26"><label>26</label><mixed-citation>
Legates, D. R. and McCabe Jr., G. J.: Evaluating the use of “goodness-of-fit” measures in hydrologic and hydroclimatic model validation, Water Resour. Res., 35, 233–241, 1999.
</mixed-citation></ref-html>
<ref-html id="bib1.bib27"><label>27</label><mixed-citation>
Luckner, L., van Genuchten, M. T., and Nielsen, D. R.: A consistent set of
parametric models for the two-phase flow of immiscible fluids in the subsurface, Water Resour. Res., 25, 2187–2193, 1989.
</mixed-citation></ref-html>
<ref-html id="bib1.bib28"><label>28</label><mixed-citation>
Martel, R., Gelinas, P. J., and Saumure, L.: Aquifer washing by micellar
solutions: 3 Field test at the Thounin Sand Pit (L'Assomption Quebec,
Canada), J. Contam. Hydrol., 30, 33–48, 1998.
</mixed-citation></ref-html>
<ref-html id="bib1.bib29"><label>29</label><mixed-citation>
Mercer, J. W. and Cohen, R. M.: A review of immiscible fluids in the subsurface: Properties, models, characterization, and remediation, J.
Contam. Hydrol., 6, 107–163, 1990.
</mixed-citation></ref-html>
<ref-html id="bib1.bib30"><label>30</label><mixed-citation>
Miller, J. J. and Curtis, D.: Electrical conductivity and soluble ions, in:
Soil Sampling and Methods of Analysis, edited by: Carter, M. R. and  Gregorich, E. G., Canadian Society of Soil Science, Broken Sound Pkwy, NW, 2006.
</mixed-citation></ref-html>
<ref-html id="bib1.bib31"><label>31</label><mixed-citation>
Mohamed, A.M.O., and Said R.A.: Detection of organic pollutants in sandy
soils via TDR and eigendecomposition. J. Contain. Hydrol., 76, 235–249,
2005.
</mixed-citation></ref-html>
<ref-html id="bib1.bib32"><label>32</label><mixed-citation>
Moroizumi, T. and Sasaki, Y.: Estimating the nonaqueous-phase liquid content in saturated sandy soil using amplitude domain reflectometry, Soil Sci. Soc. Am. J., 72, 1520–1526, 2006.
</mixed-citation></ref-html>
<ref-html id="bib1.bib33"><label>33</label><mixed-citation>
Parnian, M. M. and Ayatollahi, S.: Surfactant remediation of LNAPL contaminated soil; effects of adding alkaline and foam producing substances,
Iran. J. Chem. Eng., 2, 34–44, 2008.

</mixed-citation></ref-html>
<ref-html id="bib1.bib34"><label>34</label><mixed-citation>
Pennell, K. D., Jin, M., Abriola, L. M., and Pope, G. A.: Surfactant enhanced
remediation of soil columns contaminated by residual tetrachloroethylene, J.
Contam. Hydrol., 16, 35–53, 1994.
</mixed-citation></ref-html>
<ref-html id="bib1.bib35"><label>35</label><mixed-citation>
Persson, M. and Berndtsson, R.: Measuring nonaqueous phase liquid saturation in soil using time domain reflectometry, Water Resour. Res., 38, 1064, <a href="https://doi.org/10.1029/2001WR000523" target="_blank">https://doi.org/10.1029/2001WR000523</a>, 2002.
</mixed-citation></ref-html>
<ref-html id="bib1.bib36"><label>36</label><mixed-citation>
Redman, J. D., Kueper, B. H., and Annan, A. P.: Dielectric stratigraphy of a
DNAPL spill and implications for detection with ground penetrating radar:
Aquifer restoration, ground water monitoring and geophysical methods, in: 5th national outdoor action conference, 13–16 May 1991, Las Vegas, Nevada, 1017–1030, 1991.
</mixed-citation></ref-html>
<ref-html id="bib1.bib37"><label>37</label><mixed-citation>
Redman, J. D. and DeRyck, S. M.: Monitoring Non-Aqueous Phase Liquids in the
Subsurface with Multilevel Time Domain Reflectometry Probes, in: Proceedings of the Symposium on Time Domain Reflectometry in Environmental, Infrastructure, and Mining Applications, 7–9 September 1994, Special Publication SP 19-94, NTIS PB95-105789, US Bureau of Mines, Evanston, Illinois, 207–215, 1994.
</mixed-citation></ref-html>
<ref-html id="bib1.bib38"><label>38</label><mixed-citation>
Rinaldi, V. A. and Francisca, F. M.: Removal of immiscible contaminants from
sandy soils monitored by means of dielectric measurements, J. Environ. Eng.,
132, 931–939, 2006.
</mixed-citation></ref-html>
<ref-html id="bib1.bib39"><label>39</label><mixed-citation>
Robinson, D. A. and Friedman, S. P.: The effective permittivity of dense packing of glass beads, quartz sand and their mixtures immersed in different
dielectric backgrounds, J. Non-Crystal. Solids, 305, 261–267, 2003.
</mixed-citation></ref-html>
<ref-html id="bib1.bib40"><label>40</label><mixed-citation>
Robinson, D. A., Jones, S. B., Wraith, J. M., and Or, D.: A review of advances in dielectric and electric conductivity measurements using time domain reflectometry, Vadose Zone J., 2, 444–475, 2003.
</mixed-citation></ref-html>
<ref-html id="bib1.bib41"><label>41</label><mixed-citation>
Roth, K., Schulin, R., Fluhler, H., and Attinger, W.: Calibration of time
domain reflectometry for water content measurements using a composite
dielectric approach, Water Resour. Res., 26, 2267–2273, 1990.
</mixed-citation></ref-html>
<ref-html id="bib1.bib42"><label>42</label><mixed-citation>
Topp, C. G., Davis, J. L., and Annan, A. P.: Electromagnetic determination of
soil water content: Measurements in coaxial transmission lines, Water Resour. Res., 16, 574–582, 1980.
</mixed-citation></ref-html>
<ref-html id="bib1.bib43"><label>43</label><mixed-citation>
Truong Hong, T. and Bettahar, M.: Effect of the water-oil ratio on
brine/surfactant/alcohol/oil systems optimized for soil remediation, J. Mater Cy. Waste Manage., 2, 109–117, 2000.
</mixed-citation></ref-html>
<ref-html id="bib1.bib44"><label>44</label><mixed-citation>
van Dam, R. L., Borchers, B., and Hendrickx, J. M. H.: Methods for prediction of soil dielectric properties: a review, in: Proceedings of SPIE-the International Society for Optical Engineering, 28 March–1 April 2005, Orlando, Florida, 2005.
</mixed-citation></ref-html>
<ref-html id="bib1.bib45"><label>45</label><mixed-citation>
Wang, X. P., Quan, G. J., Pan, Y. X., Hu, R., Zhang, Y. F, Tedeschi, A., Basile, A., Comegna, A., Coppola, A., and de Mascellis, R.: Comparison of hydraulic behaviour of unvegetated and vegetation-stabilized sand dunes in arid desert ecosystems, Ecohydrology, 62, 264–274, 2013.
</mixed-citation></ref-html>--></article>
