vvEPA
                              United States
                              Environmental Protection
                              Agency
                          Environmental
                          Research Laboratory
                          Athens, GA 30613
 Research and Development    EPA/600/M-89/016    May 1990

ENVIRONMENTAL
RESEARCH    BRIEF
     Assessing the Environmental Partitioning of Organic Acid Compounds
                                          Chad T. Jafvert
 Background
 The ultimate disposition of organic compounds  in the
 environment is  influenced tremendously by their tendency
 to  partition, or  "adsorb" to sediment particles or other
 distinct environmental components or phases. To date the
 majority of sorption research for organic chemicals'has
 focused on describing  partitioning processes  of  neutral
 compounds, such as the polycyclicaromatic hydrocarbons
 polychlorinated  biphenyls,  and dioxins.  Many environ-
 mentally relevant organic compounds, however, contain acid
 or  base functional groups.  Compounds containing acid
 functional groups include the chlorinated phenoxyalkyl acid
 herbicides, such as 2,4-D and 2,4,5-T;  chloro-  and nitro-
 phenols, such as pentachlorophenol (POP); and  sulfonated
 compounds, such as azo and anthraquinone acid dyes and
 linear alkylbenzenesulfonate (LAS)  surfactants. Organic
 bases of significance include the nitrogen  heterocycles,
 such as quinoline, and aromatic amines, such as benzidine.'
This Research Brief focuses on the partitioning of organic
acid compounds in octanol-water and  in sediment-water
systems.  The main emphasis  is to  characterize  the
adsorption  processes of  these  compounds "mechanis-
tically," to permit the development of quantitative structure-
activity relationships.  In light of the differences  in
adsorption  reactions among  neutral, cationic, and anionic
organic compounds, a brief theoretical description is given
as to why  these compounds partition as they do in  the
environment.  For example, why  do certain environmental
components, such as fish  lipid, or the organic material in
soils and  sediments preferentially  adsorb neutral com-
pounds over charged species?   An explanation  is given
based on  thermodynamic considerations,  and  various
assumptions and requirements necessary to quantify these
partitioning processes are examined.
                      Equilibrium Partitioning
                      If no kinetic limitations exist, the  tendency to  partition
                      expressed as an equilibrium adsorption coefficient (Krt) is
                      directly  related  to the standard free energy of adsorption
                      AG°ads as expressed by Equation 1. (See [1,2] for a more
                      thorough discussion).

                             AG°ads  = - RT In Kd                      (1)

                      To characterize the adsorption of organic compounds to
                      environmental matrices, such as soils, sediments, and biota
                      this  standard free energy" of adsorption must be factored
                      into  adsorption  components (i.e.,  soil organic  carbon)
                      Adsorption to each of these components,  in turn, results
                      from additive effects of  various adsorption interactions or
                      energies." According to Stern [see  1,2], in the most basic
                      form these energies can be separated into an  intrinsic
                                   term (AGochem) and  an electrostatic energy
                            AGoads = AG°chem + AGoe,                 (2)

                     For neutral organic compounds,  the total energy  of
                     adsorption contains no electrostatic component, and hence
                     is  equal to the  intrinsic chemical adsorption energy
                     (Ata°chem)- For moderately large (molecular weight >  100)
                     nonpolar compounds, this adsorption energy is dominated
                     by hydrophobic forces. Hence partitioning occurs "out" of
                     hydrophihc  phases such as  water to  more hydrophobic
                     phases. Water  solubility  (corrected  for crystal energy  of
                     solids), octanol-water partitioning,  connectivity indices
                     cnromatographic retention times, and estimated molecular
                     surface areas, all have been used to index chemical activity
                     (or free energy changes) of hydrophobic compounds for
                     estimating partitioning to media such as soil orqanic carbon
                     and fish lipid.

                     Because the intrinsic chemical adsorption energy (AG°rhB )
                     for organic compounds, as a rule, results from hydrophobic
                     interactions, this term for  both cationic and anionic organic

-------
compounds is less than for neutral  compounds of [similar
structure. For these charged species, however, the primary
change in  the overall  adsorption energy  may  resu|t from
electrostatic contributions  to adsorption  (AG°ei)r  This
electrostatic contribution can be positive or negative and is
equal  to  the product of the  compound's valence (z),
Faraday's constant (F), and the surface potential  (ip). Most
natural  soils and  sediments contain mineral and  organic
surfaces that possess negatively charged functional groups.
Therefore,  the overall sign of AG°e| is ultimately determined
by the charge of the adsorbing ion.                 \
For adsorption of cationic organic compounds (z =  ^ 1), the
overall sign  of  AG°8|  is negative,  and electrostatic
interactions can  promote  adsorption  beyond  what  is
predicted for a similar neutral  compound .based  purely on
hydrophobic interactions.  Examples  of this  in the literature
include the adsorption of benzidine (4,4-diaminobiphenyl,
pKa  = 4.3) [3]  and quinoline (pKa  =  4.94)  [4]  onto
sediments  or aquifer material. The adsorption of  these
compounds to  natural surfaces is  highly  pH-dependent
because of their  relatively  low pKa values. At pH, values
where the  cationic form of these compounds predominates,
adsorption can be two orders of magnitude greater than the
hydrophobic adsorption of  the  corresponding  (neutral
species. This type of reaction is generally referred to  as
ligand  (or  cation)  exchange,   and  modeling this  process
relies on  understanding  the  adsorption  behavior of  other
exchangeable  ions,  such as Ca + 2  and Mg + 2.  Also
 adsorption can occur to various subcomponents of soils and
 sediments, such as the  organic matrix  and various; mineral
 components. To effectively  describe these  processes,
 model algorithms must be somewhat complex. Modeling the
 adsorption of organic cations may best  be accomplished by
 incorporating adsorption  algorithms into "electrostatic" or
 ligand exchange adsorption models, similar  to those used to
 estimate metals adsorption.
 For the adsorption of anionic organic compounds (z  = -1),
 the  overall sign  on  AG°ei is  positive,  and  electrostatic
 interactions act to attenuate adsorption.  Hence,  adsorption
 to natural sediments by the anionic form of an organic acid
 is generally less than  adsorption of the corresponding
 neutral  species.  Because  no  specific electrostatic
 interactions occur to promote adsorption, the partitioning
 process  must  result from the existence  of a  significant
 chemical  adsorption  energy  term.  This  situation results in
 adsorption occurring predominantly through  hydrophobic
 interactions, attenuated by electrostatic forces. Indeed, in
 limited cases [5,6], the  adsorption  of  anionic organic
 compounds has been  normalized to  the  organic carbon
 content of the sediments  and soils. In  another  study  [7],
 adsorption of LAS surfactants was shown  to increase with
 surfactant alkyl chain length  through  the  surfactant  series
 investigated, indicating  that   hydrophobicity contributes to
 the total energy of adsorption.
  In  quantitatively assessing  these partitioning  reactions,
  separating the hydrophobic component  of partitioning from
  the  electrostatic component  is essentially impossible. As
 'with neutral organic compounds,  however, model; systems
  can be used to assess and possibly quantitatively describe
  sorption  processes.  Model  systems include somie  of the
  same "systems" or parameters used to assess partitioning
  of neutral compounds,  such as solubility or octano.l-water
  distribution. Indeed, many of the same factors that influence
the transfer of organic acids between octanol and water (or
affect the apparent water solubility) influence the transfer of
these compounds between  the  components of the more
heterogeneous system  of  sediment-water [8-10].  These
similar factors include the  degree of dissociation and the
hydrophobicity and electrostatic  influences  of both sorbent
and sorbate.
As a result of these similarities, the distribution of  several
organic acids between octanol and water phases has  been
examined as  a function of aqueous phase salt concentration
(for KCI, NaCI, LiCI, MgCI2, and  CaCI2) and aqueous phase
pH [9]. Compounds that have moderately low  (<  5) pKa
values were  chosen.  From these experiments,  equilibrium
partition coefficients of the neutral  and charged forms, as
well  as  several  pKa  values have  been  calculated.  The
sediment-water  distribution of selected organic acids also
has been investigated as a function of aqueous phase ionic
composition  and pH [10.]. Experimental results  and model
development for each system is  briefly reviewed.

Octanol-Water Partitioning
The distribution  of organic acids between water and octanol
in  the  presence of monovalent inorganic cations  can be
described by Reactions 1  through  6 of Figure  1. The
reactions in  Figure 1  are formulated in terms of KCI as the
monovalent  salt, although  NaCI and LiCI  also  have  been
 used. The mass action equations for  these reactions,  along
with  the material balance and electroneutrality  constraints,
 comprise the chemical equilibrium model.
    Reactions of HA

    1.[HA]W = [H + ]w + [A-]w

    2. [HA]W = [HA]0

    Reactions in KCI, NaCI, and LiCI salt systems.

    3. [K+]w
Expression



    KOW
                                             Ki
                                              ipx
    5- [K + ]w + [A-]w = [K + ]0 + [A-]0              N

    6. [K + ]w + [A-]w = [K + A-]0                  Kip

    7. [K + ]w + [OH-]W = [K + ]0 + [OH-]0

    8. [K*lw + [OH-]W = [K*(
   Figure 1.  Equilibrium reactions of monovalent organic
             acids in octanol-water. Adapted from reference
             [8]. HA represents a neutral organic molecule
             and A" represents the corresponding mono-
             valent organic anion. Terms subscripted with
             "w" or "o" represent aqueous phase or octanol
             phase concentrations, respectively.
  Reaction 1  describes  acid  dissociation in the  aqueous
  phase. Reaction  2 describes transfer of the neutral species
  (HA) to, the octanol phase. This reaction is the basis of lineal-
  free energy relationships relating ,KOC  (the sediment or soil
  equilibrium  partition coefficient, normalized  for organic
  carbon content) and  Kow,(the octanol-water  partition
  coefficient)  for sorption  of nonpolar  organic  compounds
                                                          ,2

-------
 [11,12].  For any  icnizable  compound, this  mechanism is
 important when a significant fraction of the total compound
 exists in the neutral  form.  Reactions  3 and  4 represent
 transfer of the  free  and  ion-paired  inorganic  species
 respectively, to the octanol  phase. For both KCI and  NaCl'
 an aqueous phase concentration of 0.2 M (with no organic
 compounds present) will result in approximately 40 uM total
 salt partitioning into the octanol phase. Reactions 5 and 6
 describe the transfer  of the free and ion-paired  organic
 species, respectively, and the inorganic  counter-ion to the
 octanol phase.  The  transfer of  the organic  ion, as a free
 spec.es  or as  an ion-pair,  is highly  dependent  upon the
 counter-ion. Also, as  the aqueous  phase  inorganic 'salt
 concentration is increased, the transfer of the organic ion as
 an ion-pair is favored.

 To illustrate these effects, the overall distribution ratio (D)  of
total compound in the octanol phase to total in the aqueous
phase, can be calculated:
     ([HA]
 D =
                                                       (3)
                  "       W
  where M + represents the monovalent electrolyte K+  Na+
  ?r   ;of°rJ3AX1periments  Performed at sufficiently high pH
    %"tj  ''    M  IS an  insignificant species in  either phase
  and the overall distribution results from Reactions 5 and 6 of
  Figure 1. Using the  calculated partition coefficients Kix and
  Kipx for the inorganic  species KCI, NaCl, or LiCI  and  the
  measured  coefficients  K,,  and  Kip  for  the  organic
  compounds pentachlorophenol  (PCP), 4-chloro-o-(4-
  chlorophenyl)-benzeneacetic  acid  (DDA)   and 2 45-
  trichlorophenoxyacetic acid  (2,4,5-T),  we can obtain' the
  overall distribution  of  these  compounds (Figure 2)  as a
  function  of  aqueous phase  salt  concentration   For
  comparison, the concentration of the  organic species  [A-l
  was fixed at 10 PM. and the pH value was set at 12  Fo the
  salts  NaCl and  KCI,  little  difference in  the calculated
  distribution is  seen, whereas  partitioning into the octanol
  obv.ously is favored  for the  LiCI  system. Also,  the
  importance of the hydrophobicity  of  these compounds is
 apparent. The magnitude of partitioning into the octanol for
 the organic anions follows the sequence PCP  > DDA  >
 2,4,5-T for all three salt systems. This is the same orderinq
 as  the octanol-water  distribution coefficients  of the
 corresponding neutral species (Kow). Indeed, for ten organic
 compounds investigated, including various carboxylic acids
 and  phenols  the calculated  partition  coefficients  for
 formation of the  ion-pair in the octanol phase (KiD) correlates
 extremely well to Kow values [9].              p

 Results of four experiments  in  which pH was varied  at
 constant ionic strength  ([K + ]w  = 0.1 M)  are presented  in
 Figure 3  The  curves were calculated using  the  partition
 coefficients given  in  Figure 3  for the  organic  compounds
 and calculated  values for the  coefficients  Kix and  KiDX, for
 KCI. At low pH values (approximately pH 1  to 2), the nlutral
 species is dominant in both phases, and the value of too  D
 is  independent  of.pH.  At   intermediate  pH  values
 (approximately 4 to 6),  the neutral species is  dominant in
the  octanol phase  and the  anionic  species  becomes
dominant in the  aqueous phase, and the  distribution ratio
becomes pH dependent. At high pH (above pH 9), the anion
predominates m  both phases,  and the distribution ratio is
 again independent of pH, but dependent on ionic strength
 as previously discussed.                                '

 The significance of  each component  of this equilibrium
 vLlT  ua1-be m°re Clearly seen in a bo [concentration]
 versus pH diagram as presented in Figure 4 for 4,6-dinitro-
 o-cresol (DNOC). All  lines in this figure are calculated from
 the six equihbrium constants for DNOC and KCI, the total
 mass of DNOC and KCI in the system,  and the volumes of
 nt?nr>P   I6' The difference between the concentration of
 DNOC in the octanol  at low pH and the  concentration in  the
 water at high pH is due to the difference  in the volumes of
 the two phases.

 Sediment-Water Partitioning

 As in octanol-water partitioning, the distribution of organic
 acids between  sediment and water is highly dependent
 upon the degree of chemical dissociation. To assess the
 effects of pH  (and hence, dissociation)  on sediment-water
 distribution initial experiments have been performed using a
 pH-stat. The pH-stat can be used to titrate and maintain a
 sediment slurry to a given or set pH value. Over the course
 of an experiment, a slurry can  be sampled and the pH then
 can be lowered to a new set-point sequentially, resulting in a
 titration _ curve." Considerable changes in sediment
characteristics result from such  a titration. This approach
however, allows factoring of the overall chemical distribution
into discrete coefficients for the neutral species (KH) and the

                                            be defined
                                                                          [HA]
                                                                              sed
                                                                          [A ]
                                                                    K,. =
                                                                     di
                                                                              sed
                                                                          IA  ]„
                                                         with
                                                                   K
                                                    (4)



                                                    (5)





                                                   (6)
                                                                        [HA]..
                                                         where  [HA]sed  (mol/kg)  and [HA]W  (M) are  the neutral
                                                         species concentrations  in sediment  and aqueous phase
                                                         respectively, and  [A-]sed and [A-]w  are the anionic species
                                                         sediment and aqueous phase concentrations, respectively
                                                         Because only  total concentration  in  each  phase  can  be
                                                         determ.ned analytically,  mass balance equations  for each
                                                         phase and for the total system are needed,

                                                                             + [A~]                           (7)
                                                                                                             (8)
                                                            [HAJt = [HA]  + (m/v)[HA]                        (Q)

                                                               (mrV} H-thf- sed!"lent ^ass (kg) to aqueous volume
                                                           ratio. Combination of Equations 4 through 9 results in an

                                                        phase35'0" f°r the fraCti°n °f t0ta' compound in tne aqueous

-------
 log(D)

2.5
     :
 2.0
 1.5
  1.0
  0.5
   0.0   _
                                   .-.4
                                                                    PCP, LiCI
                                                                    PCP, NaCI
                                                                     PCP, KCI
                                                                    DDA, LiCI
                   DDA, NaCI
                    DDA, KCI
                                                                    2,4,5-T, LiCI
                                                                    2,4,5-T, NaCI
                                                                    2,4,5-T, KCI
    -.5
      0.04
                                                     I    I	L
                                                                      I	L
                                                                                J	I
                    0.08
                                  0.12!
0.16
                                                               0.20
                                                                            0.24
                                Total KJNa, or Li (M) in Aqueous Phase
Figure 2. The octano.-water distribution of 2,4,5-T, DDA and PCP using the ca.cu.ated partition coeff.c.ents
        from experimental data, and setting [A lw at 20 y.M and the pH at 12.

-------
              [HA]
        f   =
         aq
                  W
                   [HA],
                     (K
                                                     (10)
                     (Kd{H+>
     6.0-


     4.0-


     2.0-


     0.0-
    -2.0
              • POP
              °DDA
              «2,4,5-T
   Kow
   5.07  -4.85
   4.64  -3.66
   3.31  -2.83
   2.14  -4.46
 log(K)
Ka   K,
    -2.13
    -3.60
    -3.76
    -3.92
 Kip
2.67
1.84
0.60
0.016
       0.0
             2.0
                    4.0
                          6.0
1 i •
8.0
                         •  i ' •
                         10.0
    ' I ' '
    12.0
                                                   14.0
                             PH
 Figure 3.  The octanol-water distribution of four compounds
          as a function of aqueous phase pH. Symbols are
          data points, and lines are calculated by minimizing
          residuals on Kow and Ka. In experiments, the
          aqueous phase contained 0.1 M KCI, and HCI or
          KOH for pH adjustment. The total mass of
          compound per tube in each experiment was
          constant.
 Results of initial studies on several organic acid compounds
 using a silty sediment from the Ohio River (EPA sediment
 11)  suggests a  definite  pH dependence on  Kdi  [10]. In
 accord  with electrostatic considerations  described
 previously, some type of pH dependence on the adsorption
 of the anion may seem unequivocal, resulting from changes
 in the surface charge or degree of protonation of the organic
 carbon matrix. For the compounds studied on EPA sediment
 11 to date, Kdi can be expressed as
= a.(pH) + logK.
                                                    (11)
where, the average pH dependence or slope fa) is equal to
0.30,  and  K,  is the  "intrinsic"  chemical adsorption term
(given in Table 1 for each compound). Only one sediment
system was examined, however, and clearly more evidence
is needed to confirm this general pH dependence. In total 5
parameters  are used in  Equations 10 and 11. One constant,
pKa, is intrinsic to the adsorbate; two constants, Kd and K,,
describe the interactions between adsorbate and adsorbent-
one constant,  a,, is contained  in Kdi  (along  with «,)  and is
intrinsic  to the adsorbent;  and  pH (or  {H+})  is the
independent variable.

The experimental data  for four  compounds are  shown  in
Figure 5, along with model results using the constants given
in Table 1.  For  POP, DDA  and 2-(2,4,5-trichlorophenoxy)-
                                                               -2.0-
                                                               -4.0-
                                                            ra
                                               8  -6.0-
                                               o
                                               O

                                               -  -8.0-
                                                 -10.0-
                                           (=) Total Measured in
                                             Octariol Phase
                                                /           (•) Total Measured in
                                           0   6	9-,
                                      Total in Water
                                      [HA + ATw.
                                                                                               Water Phase
                                                                                \..»j»5>^    N.    Total in Octan
                                                                                                \IHA]W
                                  0.0
                             Figure 4.
                                                            	 .  .  i  ,
                                                             2.0     4.0      6.0     8.0
                                                                   Aqueous Phase pH
                                                                       10.0
                    pC - pH diagram of DNOC species partitioning
                    between octanol and water as a function of pH.
                    Symbols are data points and lines are model fit. The
                    volumes of each aqueous phase = 10 ml, volume of
                    octanol = 2.5 ml, total mass of DNOC in each tube
                    = 1.25x10'6 moles. The aqueous phase contains 01
                    M KCI, and HCI or KOH for pH adjustment.
                              Table 1. Constants Used in Modeling the Adsorption of
                                      Organic Acid Compounds to Sediment 11
                                                i\a
PCP
DDA
DNOC
Silvex
0.05
0.08
0.10
0.10
4.85
3.66
4.46
3.07
	 " u
3.0
2.45
1.93
1.87
3.33
2.31
2.83
2.03
                              aSediment mass to water volume ratio (kg/L)
                              faReference [9].

                            propionic acid (silvex), organic carbon-normalized partition
                            coefficients for the neutral species (Koc) are, within a factor
                            of two, equal to  the Kovv.  The K0? for DNOC is significantly
                            greater than its Kow, likely resuming from polar  interactions
                            with the sediment  by the nitro-groups on DNOC. A similar
                            trend exists for the anionic species.
                            Summary

                            Partitioning in octanol-water and sediment-water systems by
                            organic compounds that contain a single, relatively strong (2
                            < pKa <  5) acidic functional group has  been  examined.
                            Clearly, a significant conclusion  of this work is  that,  at  pH
                            values found in  most  natural  environments (i.e., 5  to  9)
                            partitioning can occur by both the neutral and ionic forms of
                            these  and similar  compounds  Therefore, to determine a
                            priori chemical distribution in the environment,  partitioning
                            must be factored into measurable or predictable quantities
                            for each adsorbing species,  as has been  attempted here.
                            The adsorption  of neutral  organic compounds is the
                            simplest case, wherein the  adsorption  energy  is largely
                            determined  by  hydrophobic interactions. Octanol-water
                            distribution (Kow), in turn,  has been used as a measure of
                            chemical  hydrophobicity. Octanol-water  distribution  for
                            organic anions is more complex and is  a function of both
                            chemical hydrophobicity and  inorganic ion composition. For

-------
the partitioning to sediments, an "ionic"  or  electrostatic
effect also is observed for organic ions. This effect is (most
obvious as a result of varying the pH of sediment systems,
and  may  result from changes in surface properties |(i.e.,
charge). These  and other  effects that  determine; the
adsorption of organic anions to  natural materials must be
examined further.                                   :
          • POP
          o DNOC
          • DDA
          A Silvex
 Figure 5.  The pH-dependent adsorption of four organic acid
          compounds to EPA sediment 11 (organic carbon =
          1.77%), and the predicted partitioning using   j
          Equations 10 and 11, and the values of Kd, Kj, and
          Ka given in Table 1, with oi equal to -0.3.      !
8. Westall, J.C., C. Leuenberger, and R.P. Schwarzenbach.
1985. Influence of pH and ionic strength on the aqueous-
nonaqueous distribution of chlorinated phenols. Environ. Sci.
Technol., 19:193-198.
9.  Jafvert,  C.T.,  J.C.  Westall, E.  Grieder, and  R.P.
Schwarzenbach. 1989. Distribution of hydrophobic ionogenic
organic compounds  between  octanol and water:  organic
acids. To be submitted for publication.

10. Jafvert, C.T. 1989. Sorption of organic acid compounds
to sediments: initial model development. Submitted for
publication.
11. Leo, A. Hansch, C. Elkins,  D. 1971. Partition coefficients
and their uses. Chemical Reviews. 71:525-616.

12   Lyman, W.J.  Loreti, C.P. Prediction  of Soil and
Sediment Sorption for Organic  Compounds.  Final  Draft
Report to  Office of Water, U.S. EPA,  Washington  D.C.,
1986.
 References                                       \
 1. Stumn, W. and J.J. Morgan. 1981. Aquatic Chemistry, An
 Introduction  Emphasizing  Chemical Equilibria in Natural
 Waters, 2nd Ed. John Wiley and Sons, New York.
 2. Morel, F.M.M.  1983. Principles  of Aquatic  Chemistry
 John Wiley and Sons, New York.
 3. Zierath D.L. J.J. Hassett, W.L. Banwart, S.G. Wood, and
 J.C. Means.  1980. Sorption of benzidine  by sediments and
 soils. Soil Science. 129:277-281.                    |
 4. Zachara J.M.,  C.C. Ainsworth,  L.J.  Felice, -and  C.T.
 Resch. 1986.  Quinoline sorption  to subsurface  materials:
 Role of pH and retention of the organic cation. Environ. Sci.
 Technol., 20:620-627.
 5. Nkedi-Kizza, P.,  P.S.C. Rao, and J.W. Johnson.: 1983.
 Adsorption  of  diuron and 2,4,5-T on  soil particle  size
 separates. J. Environ. Qua/. 12:195-197.             ,
 6. Ogram, A. V.,  R.E.  Jessups, L.T. On,  and  P.S.C.  Rao.
  1985. Effects of sorption  on  biological degradation rates of
  (2,4-dichlorophenoxy) acetic acid  in soils. Appl. Environ.
  Microbiol. 49:582-587.
  7. Hand, V.C.,  G.K.  Williams.  1987.  Structure-activity
  relationships for sorption  of linear  alkylbenzenesulfonates.
  Environ. Sci. Technol., 21:370-373.                 '

-------

-------
United States
Environmental Protection
Agency
Center for Environmental Research
Information             j
Cincinnati OH 45268
Official Business
Penalty for Private Use $300
EPA/600/M-89/016

-------