Next Article in Journal
Securing Flows in the River Systems through Irrigation Water Use Efficiency—A Case Study from Karula River in the Ganga River System
Next Article in Special Issue
Desalination of Irrigation Water Using Metal Polymers
Previous Article in Journal
Large-Scale Extraction and Mapping of Small Surface Water Bodies Based on Very High-Spatial-Resolution Satellite Images: A Case Study in Beijing, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Catalytic Partial Desalination of Saline Water

DCA Consultants Ltd., The Bungalow, Castleton Farm, Falkirk FK2 8SD, UK
Water 2022, 14(18), 2893; https://doi.org/10.3390/w14182893
Submission received: 13 August 2022 / Revised: 8 September 2022 / Accepted: 13 September 2022 / Published: 16 September 2022
(This article belongs to the Special Issue Desalination Treatment of Irrigation Water)

Abstract

:
More than 1 billion ha of land is adversely affected by salinization, including about 54 million ha of irrigated cropland. This study trials a batch flow, bubble column, static bed, catalytic, pressure swing adsorption−desorption, zero valent iron, and diffusion reactor train, which is designed to partially desalinate water, for use as either livestock feed water or partially desalinated irrigation water. ZVI desalination produces a partially desalinated water product, without producing a waste brine product. The trial demonstrates sequential batch processing of 50 (0.86 m3) batches of saline water (43 m3 total), using a single ZVI charge, without loss of activity. The trialed feed water contained between 1 and 9 g NaCl L−1. The average desalination was 41.84% (standard deviation was 15.61%). The optimum batch processing time is determined as being between 12 and 24 h. For each batch, the first-order forward rate constants, for both Cl and Na+ ion removal, decline with time. The study considers the irrigation application of the product water on crop yield for 70 crop varieties; the application of the technology to desalinate a saline aquifer; applications of the technology to improve the environmental efficiency of conventional desalination plants.

1. Introduction

Irrigated cropland currently accounts for around 40% of global food production [1,2]. Global food demand is expected to rise by between 60% and 100% by 2050 [1,2]. This rise in food demand results from the expected increase in global population, from 7.7 billion in 2017 [3] to 9.1 billion people in 2050 [1,2] and 11.2 billion by 2100 [3].
The majority of the future increase in food production required to meet this demand will be confined to areas that are salinized or are currently irrigated with saline water or become newly irrigated with saline water. The crop yield decreases as a function of increasing soil salinity.
More than 1 billion ha of land is negatively affected by salinization, including about 20% to 33% of the global irrigated cropland (54 million ha) [3,4]. The term salinized is used to refer to soil or water, with an electrical conductivity of >4 dSm−1 (mScm−1) [3]. The potential future global cropland land resource in areas where the soil is salinized has been estimated at >200 million ha [3].
In areas affected by soil salinization, the irrigation water may also be saline. As a very rough guide, the soil salinity is typically 10–30% greater than the irrigation water salinity [5]. However, the actual relationship is soil and location specific, varying with the salinity of the irrigation water [6,7].
Blending saline water with freshwater (or desalinated water) has sometimes been considered to increase the supply of irrigation water [8]. This approach can be effective in increasing overall crop yields, but may result in long-term soil salinization [8]. Some studies have demonstrated higher crop yields when crops are irrigated with desalinated water compared to irrigation with freshwater [9].
Conventional desalination plants use physical processes, such as evaporation or membrane separation (e.g., reverse osmosis), to produce a low-salinity product water and a high-salinity brine waste. While irrigation with a low-salinity product water will increase crop yields, disposal of the high-salinity brine may create environmental problems and may adversely affect local aquifers. An alternative solution is to use a chemical desalination approach. These approaches produce a partially desalinated product water and no high-salinity brine waste. Two groups of chemical desalination can potentially be used to produce a partially desalinated product water, suitable for irrigation water. They are CO2 desalination and zero valent iron (ZVI, Fe0) desalination [10].
CO2 desalination was first conceptualized in the mid-1990s (US Patent US6180012B1) as a method of desalinating seawater [10]. It has since evolved into a process that could be used to remove Na+ and Cl ions, associated with the waste brine from conventional desalination plants (Canadian parent CA2744264A1) [10]. More recently, it has been used in pilot plants, which are designed to strip CO2 from air [10]. This technology requires a CO2 source, which can be air, and a source of ammonia (NH3) (e.g., US Patents US6180012B1; US7531142B2; US10322371B2; world patent WO2007094691A1; Canadian Patents CA2669964; CA2744264A1; Australian Patent AU2012321057A1; European patents EP1335881B1; EP1958683A2; EP1961479A2; EP2134811A2). The products are: NaHCO3(s), Na2CO3(s), and NH4Cl(s) + partially desalinated water [10].
The simpler, alternative chemical process, ZVI (Fe0) desalination, was first reported in the academic literature in 2010 [10] (it has been described in 1 thesis, 1 Conference paper, 6 book chapters, and 11 academic papers).
The ZVI desalination process has been defined in patents (US Patents US8636906B2; US9624113B2; US9828258B2; US2018/0009678A; US10919784B2; French patent FR2983191A1; Spanish Patent ES2,598,032; UK Patent GB2520775A). The various types of ZVI desalination are described and defined in reference [10]. They use one of four different groups of reactor processes; these processes are: (i) n−ZVI + saline water in a continuous flow reactor, where the n−ZVI + NaCl is recovered by magnetic separation (Spanish Patent ES2,598,032); (ii) Emulsified n−ZVI + saline water + surfactant, where the Cl ions are removed by an undefined process (US Patents US9624113B2; US9828258B2); (iii) n−ZVI in a fluidized bed reactor, where the NaCl is concentrated in the ZVI particles (US Patent US10919784B2); (iv) diffusion bed reactors using n−Fe0, m−Fe0, or Fe0 to remove Na+ and Cl ions (UK Patent GB2520775A). This study uses a diffusion bed reactor.
Initially, it was discovered that if n−Fe0 or m−Fe0 was placed in saline water, the water salinity gradually declined [11]. Typical removals equated to <0.2 g NaCl g−1 Fe0. Removal was by inclusion in green rust structures and by concentration of NaCl (FeCl3, NaHCO3, and NaOH) in pores, contained within the hydroxy−oxide crust. This crust develops around a corroding iron particle [11]. Fen+ ions hydrolyze to form different ion polymers (e.g., Japanese patent JP6017768B2; Chinese patent CN111606464A). These polymers include mono−nuclear hydroxyl complexes (termed Fe(a)); a polynuclear hydroxyl complex (e.g., Fe(OH)24+; Fe(O2H6)3+; Fe3(OH)45+; Fe5(OH)96+), termed Fe(b) (e.g., Chinese patent CN113880208A; CN210261456U; CN210261453U; CN210261454U; CN113197151A); and an inert macro−molecular coagulation (termed Fe(c)). Fe(b) has a high flocculation activity and the flocculates can grow aggregation and accretion [11]. This polymer (Fe(b)) is an effective pollutant scavenger and may be largely responsible for pollutants removed in Fe polymer flocculants.
Two innovations changed this flocculant removal perception for ZVI desalination. They were:
  • The discovery that, if n−Fe0 was added to a conduit containing turbulent, flowing saline (seawater) water, in a concentration of about 20 g L−1 and the entrained ZVI−water mixture passed through a magnetic separator (to recover the n−Fe0 + attached NaCl + LiCl + KCl), after a contact time of <10 min, then the process produced a partially desalinated water (Spanish Patent ES2,598,032). This patent indicated that removal could involve physical adsorption rather than just chemical adsorption.
  • The discovery that if an oxygenated gas was bubbled through saline water in the presence of Fe0, the product water was partially desalinated [10]. This process discovered that the bulk of the NaCl removal involved a physical adsorption−desorption process, where the products were retained in the water and not in the ZVI bed.
These innovations resulted in the discovery that removals could equate to >1 g NaCl g−1 Fe0 [10]. This high level of removal implied that, while the removal included removal by reaction to form a green rust and removal by adsorption and concentration into adjacent pores by desorption, it may be possible to develop a desalination solution that was catalytic [11]. It also implied that a single batch of ZVI could be used to process multiple batches or volumes of water.
If it is possible to develop a commercially viable catalytic desalination solution, then the range of potential applications could be increased from the remediation of saline water, associated with hydrocarbon production and mining, to the partial desalination of saline water for irrigation and the processing of reject brine from conventional desalination plants, to allow either reuse in the desalination plant or disposal into the riparian, marine, or groundwater environment.
The end product from catalytic ZVI desalination is unknown. It is known that the product water associated with ZVI desalination shows desalination, using ion chromatography, electrical conductivity (EC) measurements, ion selective electrode (ISE) measurements of Cl and Na+ concentrations, and by water evaporation [10].
An earlier study using a ZVI desalination reactor train operated as a batch flow, bubble column, static bed diffusion reactor [10], without water recirculation, established an average of 24.5% desalination over 70 sequential, 0.24 m3 water batches, and indicated that the overall desalination may be increased if:
  • The reactor type was switched to a batch flow, bubble column, static bed, or recirculating diffusion reactor;
  • The reaction period was increased;
  • The pressure wave (and fluid circulation pattern) associated with the bubble column was shaped;
  • The ZVI was changed from steel wool to a smaller quantity of n−Fe0, or n−FexOyHz (n−Fe(a,b,c)).
This study tests this hypothesis by adopting these recommendations. This study uses a 0.86 m3 reactor train and demonstrates its efficiency using a single catalyst charge to process 50 sequential batches of water.

2. Materials and Methods

This study describes the desalination results associated with scoping pre-commercialization trials associated with a commercial-scale reactor train. All construction materials and chemicals used in the study were purchased from commercial suppliers.

2.1. Reactor

The reactor train process flow used in this study is summarized in Figure 1. It comprises a feed water tank, a product water tank, a bubble column reactor, a recycle water tank, an air compressor, and an air storage tank. These units are connected by conduits.
The reactor type is a batch flow, bubble column, static bed, recirculating diffusion reactor [10]. The reactor train has three active components: (a) a bubble column reactor; (b) a ZVI cartridge; and (c) a recycle water tank.
The reactor train is filled from a saline water feed tank. Product water is discharged to a product water tank. The reactor train filling procedure used was: (i) the valve between the recycle tank and the product water tank was closed (Figure 1); (ii) the valve between the feed water tank and the recycle tank was opened (Figure 1); and the recycle tank was filled by gravity feed. When the recycle tank was full, this valve was closed. All filling and discharge conduits were constructed from 22 mm O.D. PEX barrier pipe. Water flow was controlled using brass 22 mm lever ball valve fittings. PEX = high density polyethylene.
The recycle tank was a 900 L flat D shaped polyethylene tank (1.5 m long × 0.81 m high × 0.81 m wide) with a 40.6 cm diameter top access sealed hatch (Figure 1). Fusion welded sockets were used to provide pipe access to the tank.
A 22 mm O.D. PEX barrier pipe (Conduit 1) connected the recycle tank to the base of the bubble column reactor (Figure 1). The bubble column reactor was constructed from 12 m of 22 mm O.D. Cu0 conduit. It received air from a pressured air storage tank through a 6 mm O.D. Cu0 conduit with a gas distributer, located at an elevation 10 cm above the base of the bubble column (Figure 1). A 22 mm O.D. brass branch fitting was placed at the top of the bubble column. A 22 mm O.D. copper return conduit (Conduit 2) was used to take recycled water from the bubble column reactor to a second 22 mm O.D. brass branch fitting. This fitting gave the water access to a ZVI cartridge (Figure 1). The throughflow discharged recycled water through a 22 mm O.D. PEX barrier pipe to the recycle tank (Figure 1).
The gas exiting the bubble column passed through 1.5 m of 22 mm O.D Cu0 tubing before the discharge pipe diameter was restricted via 15 mm O.D. and 10 mm O.D. tubing to a 6 mm O.D. Cu0 vent stack (Figure 1).
An air compressor was used to pressurize and fill a pressured (0.8 MPa) gas tank (20 L) with air (Figure 1). The air compressor was used to supply 60 L air h−1 to the bubble column reactor.

2.1.1. Reactor Structuring

The reactor structure used in this study is compliant with the construction and design requirements of a multi−functional ZVI diffusion reactor described in UK Patent GB2520775A. This patent required the reactor to contain a concurrent flow of water and air in the bubble column and the flowing fluids to interact with the ZVI cartridge by diffusion.
This fluid flow generated a pressure wave in Conduit 2 and in the ZVI cartridge is compliant with the processing requirements of an “oscillating flow reactor” defined in UK Patent GB 2470764B. This patent structures the reaction environment into an upstream storage area (USA)—upstream plenum, elastic permeable barrier (EPB) and a downstream storage area (DSA)—downstream plenum (Figure 2). The EPB controls the rate of gas flow between the USA and DSA. In this study the EPB is a bubble column reactor (Figure 1), which is designed to create a product gas flow rate, which oscillates between 0 m3 s−1 and x m3 s−1. An example of flow periodicity through a bubble column is illustrated in Figure 3. The wave amplitude, wavelength, wave periodicity, and wave morphology can be modified by both changing operating conditions and reactor design.
The rising air (gas) bubbles displace an equal volume of water from the bubble column reactor into Conduit 2 (Figure 1 and Figure 2). The water flow rate in Conduit 2 resulting from this fluid displacement creates a confined, travelling pressure wave within Conduit 2. This confined pressure wave travels from the bubble column reactor to the recycle tank (Figure 2). The pressure wave interacts with the fluids in the ZVI cartridge to create an oscillating pressure and fluid flow environment within the cartridge. The travelling pressure water crest will result in both fluid flow into the cartridge and pressure increases within the cartridge. The travelling pressure wave trough will result in both fluid flow out of the cartridge and pressure decreases within the cartridge.

2.1.2. Fluid Flow within the Bubble Column Reactor

The bubble column reactor (EPB) (Figure 2) is structured to provide a concurrent flow of gas and water. A single rising gas bubble pushes the overlying water within the bubble column into Conduit 2. The rising gas bubble draws water into the bubble column reactor, through Conduit 1, to replace both the displaced water and the gas bubble, when it discharges into the gas body located above the water column.
This flow structuring (Figure 2) creates a situation where the general flux equation for gas transport in the bubble column can be described by the equation:
Jg = kg · Pbc,
where Jg = flux; kg = a permeability constant; Pbc = the pressure differential between the gas pressure in the plenum (USA) feeding the gas distributor and the gas pressure in the gas discharge plenum (DSA) located between the gas vent location and the air−water contact in the bubble column. The maximum value of Pbc is: Pbc = Pressure in the compressed tank (300,000 Pa)—Atmospheric pressure (100,000 Pa). Gas is bubbled into the water column (EPB). The gas flow rate, Jg, must oscillate (Figure 3) between a period of no gas flow (time interval between bubbles) and a period of gas flow (associated with a rising gas bubble).
Jg can only be >0 when [kg · Pbc] > 0. The gas distributor plenum (USA) has a volume, VA, and receives a constant pressured gas supply. When the plenum (USA) is fully discharged, its gas volume = VA and its gas pressure is [Pdc1]. When the plenum (USA) is full, its gas volume = VB and its gas pressure is [Pdc2]. This volume represents the pressured gas volume, which is required to displace the water column in the bubble column (EPB) and allow gas flow. The pressure in the gas distributor plenum will continually oscillate between [Pdc1] and [Pdc2].
Rises in pressure, [Pgp1], in the gas distributor plenum (USA) will be associated with decreases in pressure within the gas discharge plenum (DSA) and vice versa. The rise in pressure is associated with the receipt of the gas bubble and the associated water slug. Pressure relaxation, [Pgp2], in the gas discharge plenum (DSA) is associated with water discharge from the EPB into Conduit 2 and gas discharge from the DSA through the vent stack. These pressures can be termed [Pgp1] and [Pgp2], respectively. It follows that the pressure, Pbc, will cyclically oscillate between:
Pbc(minimum) = Pdc1Pgp1,
Pbc(maximum) = Pdc2Pgp2,
Gasification of the bubble column will reduce its viscosity and increase its permeability, kg, to gas flow. The pressured gas bubble will expand in size as it rises within the water column.
This pulsed gas flow creates a pulsed discharge of product water into Conduit 2 from the bubble column and a pulsed inflow of feed water into the bubble column through Conduit 1 (Figure 2).

2.1.3. Pressure Oscillations in the ZVI Cartridge

The pulsed flow of water through Conduit 2 to the Recycle Tank creates a fluid pressure oscillation in the ZVI cartridge. Fluid flows into the cartridge as the pressure rises and out of the cartridge, as the pressure, associated with the pressure wave, reduces.
Cyclic pressure increases are associated with adsorption onto the ZVI catalyst, while the associated cyclic pressure decreases are associated with product desorption from the ZVI catalyst [10].

2.2. ZVI or Catalyst Cartridge

An earlier study used a batch flow, bubble column, static bed, diffusion reactor, without fluid recirculation [10]. It trialed a metal polymer catalyst constructed from steel wool with a basic formulation Fe0:n−Fe(a,b,c)@urea. This study trials the use of a SiO2@n−Fe(a,b,c)@urea polymer.
The SiO2@n−Fe(a,b,c)@urea polymer was constructed by placing 400 g Iron Sulphate + 10 L FeOOH-coated quartz sand grains (0.1 to 3 mm) + 200 g Urea + water with <0.2 g/L Cl and <0.2 g/L Na+, in a 10 L vessel with an air−water contact. The vessel was left for 4 weeks at a temperature of between 4 and 15 °C. Then, 0.88 L of the damp solute mixture was placed in an acrylonitrile butadiene styrene (ABS) cartridge (0.7 m long × 0.04 mm O.D.). The cartridge was attached to Conduit 2, via a 22 mm O.D. connector. The cartridge contained between 7 and 10 g Fe. The active desalination catalyst (Figure 4) is a n−FeOxHy polymer, which coats the quartz grains. It is termed here a SiO2@n−Fe(a,b,c)@urea polymer and comprises a mixture of SiO2@n−Fe(a,b,c)@urea polymer and SiO2@n−Fe(a,b,c) polymer.
Characterization of the polymeric catalyst in terms of detailed chemical analyses and detailed morphology are outside of the scope of this study.
In this study, the following polymer terminology is used. M = metal. M(a) = mono−nuclear hydroxyl complexes, e.g., Fe3+, Fe(OH)2+, Fe(OH)2+, Fe(OH)3, Fe(OH)4 (Fe2(OH)24+ and Fe3(OH)45+), some small polymers, and instantaneous reaction products. M(b) = polynuclear hydroxyl complexes, with medium and high molecular mass, e.g., Fe(OH)24+; Fe(O2H6)3+; Fe3(OH)45+; Fe5(OH)96+, etc. M(c) = larger polymer or colloidal species, which include inert macro−molecular coagulation, of the form M(OH)x and MOOH. In practice most metal polymers are a mixture of M(a) + M(b) + M(c).

2.3. Water

The sequential reactor trials processed synthetic saline water, which was manufactured by dissolving halite in natural spring water. The water was sourced from a natural freshwater spring, at National Grid Reference: NO 02817 14399; 56°18′43″ N, 003°34′21″ W. The water was abstracted from fractured andesitic rocks and was placed in a 5 m3 storage tank (Figure 1). A typical groundwater ion composition is: Anions: Cl = 11.67 mg·L−1; N(NO3) = 11.28 mg·L−1; S(SO4) = 4.16 mg·L−1; P(PO4) ≤ 0.10 mg·L−1; F = 0.024 mg/L; N(NO2) = 0.04 mg·L−1; HCO3/CO32− ≤ 10 mg·L−1; Cations: K = 1.69 mg·L−1; Ca = 32.91 mg·L−1; Na = 6.32 mg·L−1; Al ≤ 150.0 µg·L−1; Fe ≤ 30.0 µg·L−1; Mn = 1.70 µg·L−1; P ≤ 0.005 mg·L−1; S = 4.31 mg·L−1; B = 29.40 µg·L−1; Ba = 135.60 µg·L−1; Cd ≤ 0.2 µg·L−1; Co ≤ 0.2 µg·L−1; Cr ≤ 0.2 µg·L−1; Cu = 77.7 µg·L−1; Ni ≤ 3 µg·L−1; Pb ≤ 10 µg·L−1; Si = 5.21 mg·L−1; Sr = 144.9 µg·L−1; Zn = 37.4 µg·L−1; As ≤ 5 µg·L−1; Mo ≤ 20 µg·L−1; Se ≤ 20 µg·L−1; Sn ≤ 20 µg·L−1; Sb ≤ 10 µg·L−1. The halite (NaCl) was purchased from Wickes Ltd., Perth, UK, and in addition to NaCl contains 30 ppm sodium ferrocyanide and a mixture of Ca2+, Mg2+, K+, and Fe3+ cations, plus SO42 and CO32 anions and clay minerals. These constitute about 6% of the solid halite mixture.

2.4. Measurements

The following parameters were measured for each feed water and product water sample: pH, Eh, EC, Cl ion concentration, Na+ ion concentration, and temperature. These measurements are placed in Appendix A, Table A1, Table A2 and Table A3, and Figure A1, Figure A2, Figure A3, Figure A4, Figure A5 and Figure A6.
The instruments used were:
  • ORP (oxidation reduction potential) meter (HM Digital) calibrated at ORP = 200 mV; measured ORP (oxidation reduction potential) values are converted to Eh, mV as: Eh, mV = −65.667pH + 744.67 + ORP (mV), using a quinhydrone calibration at pH = 4 and pH = 7.
  • pH meter (HM Digital) calibrated at pH = 4.01; 7.0; 10.0.
  • EC (electrical conductivity) meter (HM Digital meter calibrated at EC = 1.431 mScm−1).
  • Cl ISE (Ion Selective Electrode); Bante Cl ISE, EDT Flow Plus Combination Cl ISE; Cole Parmer Cl ISE attached to a Bante 931 Ion meter. Calibration was undertaken using 0.001, 0.01, 0.1, and 1.0 M NaCl calibration solutions.
  • Na+ ISE (Ion Selective Electrode); Bante Na ISE, Sciquip Na ISE, Cole Parmer Na ISE attached to a Bante 931 Ion meter. Calibration was undertaken using 0.001, 0.01, 0.1, and 1.0 M NaCl calibration solutions.
  • Temperature measurements were made using a temperature probe, attached to a Bante 931 Ion meter.

2.5. Control Reactor

2.5.1. Control Experiment 1

A control reactor was placed in operation, containing a gas distributor, with no air feed and a single ZVI cartridge. The product water salinity was sampled after 48 h. The feed water and product water salinities remained unchanged at 0.88 g Cl L−1 + 0.43 g Na+ L−1.

2.5.2. Control Experiment 2

A control reactor was placed in operation, containing a gas distributor and an air feed of 60 L h−1 and no ZVI cartridge. The product water salinity was sampled after 48 h. The feed water and product water salinities remained unchanged at 0.88 g Cl L−1 + 0.43 g Na+ L−1.

2.5.3. Control Experiment 3

A control reactor was placed in operation, containing a gas distributor and an air feed of 65 L h−1 and a ZVI cartridge. The product water salinity was sampled after 8.95 h. The feed water salinities reduced from 0.88 g Cl L−1 + 0.43 g Na+ L−1 to provide a product water salinity containing 0.57 g Cl L−1 + 0.21 g Na+ L−1. This control experiment indicated that no desalination could be expected with this reactor construction, unless the reactor contained both a flow of air and a ZVI cartridge.

3. Results

The quantitative data results generated by this study are provided in Appendix A, Table A1, Table A2, Table A3, Table A4, Table A5 and Table A6, and Figure A1, Figure A2, Figure A3, Figure A4, Figure A5, Figure A6 and Figure A7.

3.1. Relationship between Electrical Conductivity and Salinity

There are strong statistical relationships between measured salinity (g L−1) and electrical conductivity (EC), for both the feed water and the product water (Appendix A, Figure A1). The feed water is artificial saline water constructed using halite. This construction results in the 1 mScm−1 approximating to 0.95 g NaCl L−1 (Appendix A, Figure A1). When saline water is constructed by dissolving pure NaCl in pure water, the EC relationship will approximate: 1 mScm−1 = 0.54 g NaCl L−1 (e.g., US Patent US10919784B2).
The product water salinity, measured as EC or g NaCl L−1, is lower than the feed water salinity (Appendix A, Figure A1).

3.2. Removal of Salinity

The trials demonstrated a significant removal of salinity (chloride + sodium ions), with a median removal of about 43%, with 75% of trials removing >30% (Table 1; Appendix A, Figure A2). There is a weak statistical trend of increasing ion removal with increasing feed water salinity (Appendix A, Figure A2).

3.3. Change in Reactor Efficiency with Increasing Number of Trials

Each sequential water batch has a different degree of desalination (Appendix A, Figure A3), with no apparent loss of desalination efficiency after 50 batches (Appendix A, Table A2). This outcome variability (Appendix A, Figure A3) will require the water from sequential or parallel batches to be aggregated to form a finished product water.

3.4. Change in Redox Parameters

The desalination process is associated with an increase in pH and a decrease in Eh (Appendix A, Table A1; Appendix A, Figure A4 and Figure A5). It is independent of temperature over a temperature range 273 K to 293 K (Appendix A, Figure A5).

3.5. Salinity Removal as a Function of Reaction Time

A graphical plot of the observed rate constant (kobsf) versus time (t) spent in the reaction environment, where t = n establishes a general decline in (kobsf) with increasing t, which could be described using a first-order reaction, regression equation (Appendix A, Table A3), where:
kobsf = (−Ln (Ct=n/Ct=0))/(t = n),
Ct=0 = ion concentration at t = 0, g L−1; Ct=n = ion concentration at t = n, g L−1; (t = n) = time spent in the reaction environment in seconds.
A best fit regression analysis was made, using the MS Excel 2019 trendline functions of kobsf vs [a], where [a] = 3600 × (t = n) ([a] is expressed in hours) for each of the Cl ions and Na+ ions for each sequential water batch. This analysis is summarized in Appendix A, Table A3.
These regression equations indicate that:
  • The amount of ion removal increases with increased reaction time (Appendix A, Figure A6).
  • The amount of ion removal associated with a specific reaction time is variable (Appendix A, Figure A6).
The regression equations (Appendix A, Table A3) allow the outcome variability (or outcome probability distribution) to be assessed as a function of reaction time (Appendix A, Figure A7). This analysis indicates that:
  • The median desalination increases for the initial 10–12 h; thereafter, the amount of desalination is stable (Appendix A, Figure A7).
  • There is a wide range of desalination outcomes for a specific reaction time (Appendix A, Figure A7). Therefore, aggregation of sequential and parallel trains will be required to produce a stable product composition.

3.6. Expected Desalination Outcome Associated with a Multi-Train Reactor

The analysis in Appendix A, Figure A7, indicates that a multi-train reactor could be operated with a 20 h reaction time, processing 0.86 m3 d−1 train−1. This is expected (Appendix A, Figure A7) to provide each batch with a median desalination outcome of 40.39%; a mean desalination outcome of 39.89%; and standard deviation of 15.14%.
An agricultural holding of 10 ha may require 50,000 m3 a−1 of irrigation water. If the irrigation was split over a 90-day period, this would require an average of 555 m3 d−1. This could be achieved (with excess capacity) using a 650–1000 train reactor. A 650 train reactor would be expected to have a mean desalination outcome of 39.89% and standard deviation of 3.86%. The expected desalination outcome for the daily aggregated product water is provided in Figure 5. This indicates that there is a 90% probability that the amount of desalination will exceed 35%. An example set of plant layout requirements for a plant processing 555 m3 d−1 is provided in Figure 6. The actual layout will be site specific.

4. Discussion

This study demonstrated that 1 t Fe (held as Fe polymers) would be expected to process 4.3 million m3 saline water (i.e., 43 m3 processed in 50 batches using 7–10 g Fe). This compares with the proven 3.86 million m3 water processed t−1 Fe by a 22,500 m3 d−1 ZVI water treatment plant, removing dissolved organic material [11].
Furthermore, this study established that a reactor containing 7–10 g Fe can remove 51.6 kg Cl and 22.2 kg Na+ ions (Appendix A, Table A2). This removal process does not deteriorate with an increased number of sequential batches (Appendix A, Figure A3). This observation is interpreted as indicating that the desalination process leaves the FeOxHy species unchanged and that the FeOxHy species act as a catalyst.
The exact nature of the catalytic products of desalination has not been elucidated in this study. They are not present in the ZVI cartridge and are assumed to be present in the product water. An interpretation of the nature of the desalination products, associated with this type of desalination process, is provided elsewhere [10].
Further, 2.2 L samples of the product water (and feed water) were placed in storage (at ambient temperatures, in a range from 0 to 20 °C) and their salinities remeasured after 2 years. The feed water samples were used to provide an additional measurement control. These analyses indicated that the majority of the water samples continued to desalinate following removal from the reactor (Figure 7a). However, a small proportion of the stored samples (<15%) showed some level of re-salination (Figure 7b). This implies that at least some of the removed salinity is retained in the product water in another form.

4.1. Nature of the Desalination Process

In this study, the desalination process required the presence of a pulsed oxygen recharge of the water, which created a pulsed water circulation and a pressure wave, which interacted with the ZVI cartridge. The control experiments indicated that each of these elements must be present before desalination is observed.
The control experiment observations indicate that the desalination process can be operated using a catalytic, pressure−swing−adsorption−desorption process [10]. Catalytic pressure−swing−adsorption−desorption (CPSAD) processes are not widely used in reactors. However, the basic principles of their operation are documented elsewhere (e.g., UK Patent GB2470764B), specifically:
  • Using Cu0, a CPSAD process is demonstrated to increase reaction yields by up to 400% for a specific space velocity [12];
  • The reactor used in this study falls within the group of reactors that combine adsorption, reaction, and separation, within the pressure swing process [13,14].
The observation that desalination may continue following removal of the product water from the reactor (into a sealed quiescent environment) indicates that that the catalytic nano-agents are contained in the product water or that there is a micro-biological element to the removal process.

4.2. Indicative Plant Sizes and Process Costs

A single standalone 860 L reactor train had the following manufactured cost structure, excluding labor and site costs:
  • D−tank (900 L capacity) = USD 400
  • Compressor = USD 200
  • Pipework, fittings, and conduits = USD 150
  • ZVI for the ZVI cartridge = USD 20
Total capital cost = USD 770. Operating energy (electricity) consumption = 0.3 kW h−1 (6 kWh for 20 h operation). Current UK grid supplied electricity prices are around USD 0.5 kWh. Switching to an off-grid, solar-battery power complex, and amortizing the cost to around 20 years can potentially reduce the energy cost (excluding pumps) to around USD 0.1 kWh to USD 0.3 kWh (assuming no maintenance is required). This potentially reduces the energy cost for a 20 h batch processing time to around USD 0.6 to USD 1.8 batch−1 (USD 0.7 to USD 1.4 m−3) from around USD 3 batch−1 (USD 3.5 m−3).
This cost structure (without considering economies of scale (e.g., Figure 5 and Figure 6)) will limit the process applications to markets or operations, which can withstand a cost structure of USD 0.6 to USD 5 m3 for partially desalinated water.

Impact of Economies of Scale

Economies of scale can arise in the area of tank costs. For example, switching from an individual tank for each bubble column train to a single shared tank could reduce the tank cost expressed in USD m−3, from around USD 440 m3 to a lower value. For example, a 1000 m3 galvanized corrugated steel tank (diameter = 17.51 m; height = 4.59 m), with an EPDM or butyl rubber lining, costs around USD 30,000 + assembly costs. The tank is constructed on any flat, firm ground or on a concrete base. The 1 mm thick butyl liners have a life expectancy of >30 years.
Similar cost savings resulting from economies of scale are associated with the conduits, pumps, valves, and compressors.
A plant processing 1000 m3 d−1 would need a power supply capable of delivering around 125–180 kW h−1. The capital cost of a renewable energy plant capable of producing an average of 125–180 kW h−1 (including energy storage) varies as a function of the technology used, location, and supplier. A reasonable budget estimate range for an off-grid unit with a 20-year life expectancy is probably around USD 2 to USD 6 million + annual service + site costs + engineering costs + local regulatory costs + environmental costs/taxes + maintenance costs + feed water and product water tanks, pumps, and distribution costs. This unit would expect to power the delivery of around 1.1 m3 m−1 of compressed air. For a plant processing around 7.3 million m3 over 20 years, this would provide an expected energy cost of between USD 0.2 and USD 0.7 m3.
Energy and operating manpower costs are considered to represent the largest element of the cost associated using this process to produce partially desalinated water. On the basis of current costings (excluding location-specific costings), a budget product water cost estimate within a range USD 0.7 to USD 2.0 m−3 may be realistic for a plant processing 500–1500 m3 d−1.

5. Practical Implications

The process has a number of potential applications. They include:
  • Partial desalination of saline irrigation water;
  • Partial desalination of saline livestock feed water;
  • Partial desalination of reject brine from a conventional desalination plant;
  • Partial desalination of feed water for a conventional desalination plant;
  • Partial desalination of a saline aquifer;
  • The injection of partially desalinated water into an aquifer to reverse or check seawater incursion into a freshwater aquifer.

5.1. Partial Desalination of Irrigation Water

Saline irrigation water has a negative impact on crop yield and crop value, when compared with irrigation using freshwater [5]. However, irrigation with saline water may provide a higher crop yield than no irrigation. The number of crop varieties that can be grown decreases with increasing irrigation water salinity.
Saline irrigation results in a gradual salinization of the soil. This is associated with a decrease in the soil’s vertical permeability and confines the principal throughflow to the topsoil or the upper part of the subsoil. As a rough guide, the pore water salinity of the soil will be 20 to 30% higher than the salinity of the irrigation water.
Crop yield is a function of many parameters, of which irrigation water salinity is an important factor. Other factors, such as climate, soil type, and agricultural practice, can have a major impact on yield. As a rough indicative guide, the expected crop yield, Ye, can be defined [5,15] as:
Ye = Ys (1−(a(1 + c)S − b)),
Ys = expected crop yield when irrigated with freshwater (t ha−1 or another yield unit). S = Irrigation water salinity. [a], [b], and [c] are constants; the values of these constants used in this study are derived from the graphs in [15]; (1 + c)S = soil water salinity. Appendix A, Table A4, provides an indication of indicative crop yield versus irrigation water salinity for a variety of crops.
This study indicates (Table 1, Figure 5) that an average of 40% desalination may be possible. This changes the indicative crop yield equation to:
Yd = Ys (1 − (a(1 + c)(1 − e)S − b)),
where [e] is the expected desalination. Yd = the expected crop yield using partially desalinated water. Desalination increases the expected yield of a crop. This is demonstrated in Appendix A, Table A5.
Irrigation with partially desalinated water will not be an economic solution for all crops, which are currently irrigated with saline water. The expected increase in crop yield due to irrigation with desalinated water can be assessed as:
Increase in crop yield due to desalination =Yd − Ye,
Appendix A, Table A5 provides an indication of the expected increase in crop yield associated with desalination for a variety of crop types, expressed as a proportion of the crop yield, when the crop is irrigated with freshwater. For most vegetable and fruit crops (and staple crops, such as rice, potato, and tomato) currently irrigated with water with a salinity of 1 to 3 g L−1, partial desalination of irrigation water may provide a yield increase, which is equivalent to 30–60% of the freshwater irrigation yield (Appendix A, Table A5). A yield increase, which is equivalent to 20–50% of the freshwater irrigation yield, may arise for some grass crops, which are currently irrigated with water with a salinity of 1 to 5 g L−1 (Appendix A, Table A5).

5.2. Partial Desalination of Saline Livestock Feed Water

About 25% of the global land area is used to rear livestock [10]. Different animals have different degrees of tolerance to different water salinities [15]. Their tolerance to water salinity can be split into three categories [15]: (i) safe, where they are tolerant to the salinity and gain weight; (ii) suitable, where they can tolerate the salinity, but will only gain weight at a slow rate; (iii) poor, where the high salinity will result in weight loss. The salinity ranges associated with these categories vary with animal type and variety. The boundaries associated with each of these categories are subjective and will vary with geographical region, livestock type, soil conditions, available feed type, and environmental conditions.
What this study demonstrates is that in the event of drought or another environmental event, which results in salinization of livestock feed water, it will be possible to reduce the salinity of the associated livestock feedwater pools by up to 40% over a short timeframe (e.g., <24 h). Partial desalination of livestock water may increase the number of livestock that can be farmed on some land (number ha−1).
Cattle may cost USD 300–700 head−1 (pedigree > USD 2000 head−1). In many drought-affected areas, a small reduction in feed water salinity (e.g., 10–20%) may be sufficient to allow for livestock survival. This study established (Appendix A, Table A3; Appendix A, Figure A7) that this can potentially be achieved after about 1 to 3 h processing (Figure 8). Cattle may require 50–200 L d−1 head−1 of feed water. The single train reactor trialed in this study would be able to process 5000 L d−1 (with an average processing time of 3 h) to 8600 L d−1 (with an average processing time of 2 h). Further, 3 h processing would give a 50% probability of achieving a 12–26% desalination (Figure 8) and 2 h processing would give a 50% probability of achieving a 9–22% desalination (Figure 8). This analysis indicates that the process may provide a low-cost method of providing a rapid salinity reduction for livestock feed water.

5.3. Partial Desalination of Reject Brine from a Conventional Desalination Plant

The observed average 40% desalination may allow reject brine from a conventional desalination plant to be partially desalinated and returned to the conventional desalination reactor. The reject brine volume is typically 80% of the feed volume to the conventional desalination reactor. Processing this reject brine in a ZVI desalination reactor may allow all or part of this reject brine to be recycled to the conventional desalination reactor or may allow it to be discharged into the environment, without increasing the salinity of the surrounding groundwater, riparian water, or seawater (Figure 9).
A conventional reverse osmosis plant producing 140 m3 d−1 of potable water may have a feed water consumption of 690 m3 d−1 of brackish or saline water and produce 550 m3 d−1 of reject brine (Figure 9). The ZVI desalination process has the potential to reduce the salinity of the reject brine to below that of the feed water (Figure 10). This could potentially result in the average water salinity entering the reverse osmosis plant being about 50% of the salinity of fresh feed water entering the system (Figure 10).
The integrated plant configuration shown in Figure 9 was not trialed at this time. Therefore, it is not known whether the modelled results (for the process configuration in Figure 9), shown in Figure 10, will occur. Additional treatment of the recycled water (Figure 9) leaving the product water storage tank may be required. Additional trials will be required to prove the conceptual process flow shown in Figure 9.

5.4. Partial Desalination of Feed Water to a Conventional Desalination Plant

The operating (and capital) costs of a conventional desalination plant decrease as the salinity of the feed water decreases. Processing the water in a ZVI desalination reactor prior to entry into a conventional desalination plant has the potential to reduce the capital and operating costs associated with a conventional desalination plant. This is achieved by reducing the salinity of the feed water entering the conventional desalination plant by 40% (Table 1; Appendix A, Figure A7).

5.5. Partial Desalination of a Saline Aquifer

Desalination of an aquifer can be undertaken by enclosing an area of the aquifer by a rim marked by or defined by the presence of a series of infiltration wells or boreholes (Figure 11). This is termed the desalination zone. Water is abstracted from the center of this field and is processed in a ZVI desalination reactor before being returned to the aquifer via the infiltration wells (Figure 11). An associated process flow diagram is provided in Figure 12.
The number of infiltration wells and their spacing will, for a specific aquifer abstraction rate, control the amount of inflow into the desalination zone (Figure 11) from the wider aquifer. This, in turn, will control the amount of aquifer desalination that will be achievable (Figure 13). This analysis indicates that it may be possible to target converting saline aquifers into repositories of partially desalinated water, which could be used for groundwater abstraction for irrigation.
For example, if the aquifer salinity is 5%, the normal influx from the aquifer to the desalination zone is 5% and the stable average salinity in the partially desalinated aquifer (desalination zone) after 1 year is 1.4%, then abstraction from the desalination zone, where part of the water is used for irrigation and part of the water is reinjected (Figure 12), will result in a gradual increase in the irrigation water salinity (Is), such that 1.4% < (Is) < 5%. Following cessation of water abstraction for irrigation, the aquifer salinity will continue to decline.

5.6. Reversal of Seawater Incursion

Seawater incursion into an aquifer can be reversed by injecting less-saline water into the aquifer (Figure 14). The general process flow (Figure 11 and Figure 12) will result in a zone of reduced salinity developing between the seawater incursion front and the infiltration/injection wells within the aquifer (Figure 14).

6. Novelty

6.1. Water Volume Processed

ZVI desalination normally requires a ZVI concentration of around 20–100 g L−1 [10,15,16] and is commonly used to treat a single batch of water [15]. A typical treatment ratio is between 10 and 50 L of water for each 1 kg of Fe used [15,16]. Reference [10] demonstrated that it was possible (for desalination), using a bubble column reactor, to increase this ratio to 42 m3 kg−1 Fe.
This study demonstrated a significant advance in the amount of saline water that can be processed by a unit weight of Fe, from 42 m3 kg−1 Fe [10] to 4300 m3 kg−1 Fe. This represents an increase in Fe utilization by 2 orders of magnitude.
This was achieved by making the following assumptions:
  • The desalination process involves catalytic sites on Fe0 or FexOyHz species [10];
  • The catalytic process involves a physical reactant adsorption and a physical product desorption [10];
  • Only a small amount of energy is required to force adsorption of a reactant and desorption of a product and this process can be achieved through the use of a pressure wave [10];
  • Adsorption is favored by the pressure rise created by the crest of the pressure wave [10];
  • Desorption is favored by the pressure trough created by the trough of the pressure wave [10];
  • Different species may compete for the same adsorption sites [10];
  • Different species may desorb from the sites at different rates. This results in a reduction in the number of available adsorption sites with time. This results in a decline in the desalination rate with increasing reaction time [10];
  • The principal interaction between the water body and the catalyst is by diffusion [10]. This assumption allows the catalyst to be placed in a removable cartridge, which is attached to the water body;
  • The pressure waves can be constructed by gases bubbling through the water body, displacing water. Reference [10] adopted an approach where 60 L h−1 was bubbled through a conduit with a 0.45 m I.D. This created an unstructured and complex wave pattern with significant interference waves. In this study, the same air flow rate was confined to a conduit with a 0.02 m I.D. This created a more consistent and compact wave pattern, which was transmitted along the conduit.
  • Changing the reactor construction type from a batch flow, bubble column, static bed, diffusion reactor [10] to a batch flow, bubble column, static bed, recirculating, diffusion reactor will improve pressure wave transmission efficiency.
This category of reaction process is a catalytic pressure swing adsorption desorption process [13,14] or a sorption-enhanced process [12].
This study used air to create pressure waves and it is possible that O2 plays a role in the desalination process [10]. An alternative source of a pressure wave is an ultrasound source. ZVI degradation experiments for p-nitrophenol (PNP) established that degradation rates are enhanced by the presence of ultrasound [17] and for a constant ultrasound frequency of 19.87 kHz, the reactor shape controls the degradation rate [17]. The degradation rate is higher in a conical reactor than in a cylindrical reactor [17]. The proportion of PNP removed for any specific reaction period was directly proportional to the amount of power applied (between 100 W and 500 W), where the regressed removal (vs. power) followed a power function [17].
This observation [17] implies that it may be possible to accelerate the desalination rates documented in Appendix A by increasing the air flow rate, by changing the reactor shape, changing the location of the ZVI, by adding another wave source, or a combination.

6.2. Nature of the Desalination Process

The exact nature of the desalination process is unknown. An earlier study [10] established:
  • The process reduces the measurable concentration of Na+ ions, Cl ions, and EC in the water;
  • Evaporation of the product water confirms removal of NaCl;
This study confirmed that the measured first-order rate constant declines with increasing reaction time (Appendix A, Table A3). This is an established phenomenon associated with diffusion reactions [18]. The controls on this decline were not investigated in this study. The decline does control the variability in batch salinity outcome (Appendix A) and the overall outcome probability distribution as a function of time (Appendix A, Figure A7). Any decrease in the rate constant decline will boost the average desalination achieved for a specific reaction time. This is, therefore, an important future potential area of research.

6.3. Technology Readiness

The trial results were for a commercial-scale reactor train (Appendix A). They established that the reactor train could achieve a 40% desalination of water. This is the first ZVI desalination reactor train to achieve the desalination of 43 m3 of saline water. Its reactor increased the batch volume processed from 0.24 m3 batch−1 [10] to 0.86 m3 batch−1, while reducing the air requirement from 250 L m−3 h−1 [10] to 69 L m−3 h−1 and the Fe catalyst charge from 1.67 kg m−3 [10] to 0.012 kg m−3. After 3 h operation, the 0.24 m3 reactor train achieved [10]: mean desalination = 24.5%; standard deviation = 22.8%; minimum = 0%; 1st quartile = 9.2%; median = 14.3%; 3rd quartile = 33.9%; maximum = 91%. As a direct comparison, after 3 h, the 0.86 m3 reactor train achieved: mean desalination = 19.40%; standard deviation = 9.41%; minimum = 0.28%; 1st quartile = 12.36%; median = 19.33%; 3rd quartile = 26.61%; maximum = 49.93%.

6.3.1. Reducing the Outcome Variability

The process design changes used in this study relative to the 0.24 m3 reactor increased the expected downside desalination from 9.2% to 12.36% at the 1st quartile. The median desalination expectation was increased from 14.3% to 19.33%. The flip side of these outcome changes was a reduction in the upside desalination potential. These changes changed the standard deviation/mean outcome ratio from 0.93 for at 0.24 m3 reactor to 0.48 for a 0.86 m3 reactor. This compares with the initial reactor structuring trials, which used a 0.01 m3 reactor and achieved a standard deviation/mean outcome ratio (Mo) of 2.0 (i.e., a very wide output range). The data regresses to provide a relationship of the form (R2 = 98.1%):
Mo = 0.513359(reactor volume, m3)−0.305178,
The expectation from these results is that increasing the reactor train size into a range 2 to 10 m3 will further reduce the standard deviation/mean outcome ratio and further narrow the batch outcome range around the mean value. This hypothesis has still to be confirmed by sequential trials processing 10 m3 batch−1. The expectation from Equation (7) is that Mo for a 10 m3 reactor will be 0.254 and that this will reduce with a 100 m3 reactor to 0.1259. For a normal distribution, the 1st quartile value is: the mean value −0.674 (standard deviation).

6.3.2. Switching from Batch Flow Operation to Continuous Flow Operation

The process trials have all operated using a batch flow. Automation would be made simpler if the overall water flow management was switched from a batch flow to a continuous flow. An example structuring for a continuous flow reactor is provided in Figure 13. This structuring allows the volume of water (m3 d−1), which can be processed by the reactor, to be increased by about 20%. The effectiveness of this structuring (Figure 13) has still to be confirmed by long-term testing.
Initial small-scale testing (using a different ZVI catalyst) of this process flow (Figure 15) was undertaken using 1 L and 2.5 L reactors with a continuous flow feed rate. These small-scale trials indicated that the process flow shown in Figure 15 can work on a small scale, over trial periods of up to 14 days. Commercial-scale reactor train trials of the process flow illustrated in Figure 15 were not undertaken at this time, with a continuous water feed.

6.4. Shaping the Pressure Wave

The modelling and design of pressure waves can be simple or complex. From a design perspective in a catalytic–pressure–swing–adsorption−desorption environment, the principal objectives are to maximize the wave frequency while maximizing the wave amplitude.
The permeability of the water to the gas is kg. The bubbling of the gas through the water column creates a primary pressure wave. All continuous flow reactors (fixed bed and fluidized bed) have a constant driving force between an upstream plenum (USA) and a downstream plenum (DSA). This is achieved by maintaining a constant pressure in the upstream plenum (USA) and the downstream plenum (DSA). In the reactor trains used in this study, the pressures in the upstream plenum (USA) cycle between lower and higher pressures and the pressures in the downstream plenum (DSA) cycle between higher and lower pressures. This creates a variable driving force across the water column. At any instant in time, the pressure in the upstream plenum (USA) is greater than or equal to the pressure in the downstream plenum (DSA). The relationship between the pressure in the upstream plenum (USA) and the pressure in the downstream plenum (DSA) is in a continual state of flux. This continual pressure variation creates pressure waves within the bubble column reactor and the discharge conduits associated with it (Figure 1 and Figure 2).
This reactor structuring has a number of characteristics. These include:
  • Kinetic energy received in the upstream plenum (USA) is temporarily converted to potential energy. The amount of kinetic energy converted to potential energy cyclically varies with time.
  • Potential energy stored in the upstream plenum (USA) is cyclically converted to kinetic energy.
  • The flow rate through the water column cyclically varies with time. The maximum flow rate through the water column in a flow cycle is greater than the average flow rate into the upstream plenum (USA). The minimum flow rate through the water column in a flow cycle can be zero (Figure 3).
This structuring is highly flexible and can accommodate a feed to the upstream plenum (USA), which is delivered at a constant rate, a variable rate, an irregular rate, intermittently, or episodically. This structuring can also accommodate a demand for product from the reactor, which is continuous, intermittent, or variable.
The bubble column reactor (Figure 1) is structured to contain a water column: an upstream plenum (USA), a downstream plenum (DSA), a downstream conduit, which removes product from the downstream plenum (DSA), and an upstream conduit, which delivers a feed to the upstream plenum (USA) (Figure 2).
The difference in potential energy between the upstream plenum (USA) and the downstream plenum (DSA) (less any pressure (potential energy) losses) is the driving force for the fluid flow. The driving force (ΔP) is approximated as:
ΔP = PU − PD − PL,
where PL = pressure losses (Pa) associated with the flow across the water column; PU = the pressure (Pa) in the upstream plenum (USA); PD = the pressure (Pa) in the downstream plenum (DSA). The pressure PU is exerted by a standing (static) body of fluid (which can be pressurized), a flowing fluid (which can be pressurized), or a combination thereof.
The pressure (Pa) exerted by a static fluid on a basal surface (e.g., a fluid−fluid contact or a membrane or a surface) is calculated as: dfgh, where df = fluid density (kg m3); g = acceleration due to gravity (m s−1); h = a column (m) of a given fluid of mass density df, which would be necessary to produce this pressure.
Conventional reactors require that ΔP, PU, PD, and PL through the catalyst bed remain constant and do not change with time. In this reactor (Figure 1), ΔP, PU, and PD cyclically change from higher to lower values with time. PL may change with time or may remain constant.
The fluid flow rate through the water column is described by the equation:
Q = k ΔP, or Q = k (ΔP − Δχ),
where Q = flow rate (m3 m−2 s−1); k = permeability (m3 m−2 s−1 Pa−1); ΔP = driving force (Pa); Δχ = correction for fugacity, or another form of correction. Permeabilities may change with flow type (e.g., viscous flow, diffusion flow, turbulent flow, laminar flow) and fluid type. Flow type (and permeability) can change with ΔP. Flow can be turbulent or laminar. ΔP across the water column must vary with time.
The water column is used to create a continual disequilibrium within the reactor, where ΔP, PU, PD, and PL are in a continual state of flux. This results in the fluid flow rate, Q, through the water column cyclically varying with time, even if the flow rate into the reactor through the upstream conduit is constant (Figure 3). The flow rate through the downstream conduit can remain constant or may vary with time (Figure 3).
An ΔP wave can be constructed and designed or may have a random construction or a combination thereof. The ΔP wave can assume any shape, amplitude, and frequency, including, but not limited to, symmetrical, asymmetrical, sinusoidal, triangular, harmonic, Fourier, square, spiral, scroll, trigger, pulse, a wave train, another type of wave, or a combination. The wave can be a shock wave resulting from an abrupt change in force level, an impulsive force, pulsed force, or a combination thereof. The wave form can be damped, undamped, underdamped, overdamped, or a combination thereof. The wave amplitude can decline with time, increase with time, remains constant with time, or a combination thereof. The wave can be a kinematic wave, a phase wave, or a combination thereof. The wave can be a pulse wave.
Higher periodicity waves can be supplemented with shorter periodicity waves. Wave structures can be composite. A pressure management system can optionally be used to control or define one or more aspects of the waves. Waves can be supplemented by artificially induced waves associated with ultrasonic equipment, wave generating equipment, heating, or cooling of fluids located in one or more of the upstream plenum, upstream conduit, downstream plenum, downstream conduit, and water column. Wave periodicity can take any form and can have any duration.
In the reactor design used in this study (Figure 1 and Figure 2), the ΔP wave cyclically alters PU and PD. In the earlier reactor design [10], the ΔP wave cyclically altered PU but left PD unchanged as a function of time.
The wave can be generated by alterations to k in the water column or by the use of apparatus, which has the effect of cyclically altering PU and PD or a combination thereof. The alterations to k can be induced by changes in ΔP; or by chemical changes in the water column; or by changes in the adsorption or desorption of specific fluids; or by osmotic changes within the water column.
PM is the minimum driving force or pressure required to initiate viscous liquid flow. For gases, the minimum driving force (PM) can be defined by the Knudsen Number (KN) or another parameter associated with a specific type of diffusion. When the driving force is less than PM, the principal flow mechanism is bulk diffusion or Knudsen Diffusion. ΔP < PM during part of a flow cycle and ΔP > PM during part of a flow cycle.

Pressure Applied by the Water Column

In this reactor (Figure 1), the pressure applied by the water column on the exit of the upstream plenum (USA) was 0.117 MPa (above atmospheric pressure). A gas flow could only be initiated if the buoyant pressure in upstream plenum (USA) exceeded 0.117 MPa (above atmospheric pressure). This initiation pressure reduces as the volume of gas entrained within the water column increases. By way of comparison, the pressure exerted by the water column in the 0.24 m3 reactor train described in reference [10] was 0.0147 MPa.
Experimental work established that the rate of lateral attenuation of the pressure wave created by the rising bubble column decreases as the effective bubble size increases [19]. For a given bubble size, the amplitude of the pressure wave decreases exponentially with increasing lateral distance. For example, the amplitude of the pressure wave (Ax) at a lateral distance x from the bubble can be expressed as [20]:
Ax = A0 e(−η(x−x0),
A0 = initial wave amplitude; x0 = initial bubble position; η = the attenuation coefficient of the pressure wave. In a conduit, the attenuation coefficient decays as a function of x, as [20]:
η = ln(A0/Ax)/(x − x0),
The pressure wave propagation velocity is a function of the air fraction in the bubble column and increases as the void fraction exceeds 0.5 [20]. η increases as the void fraction increases [20]. While the modelling of bubble flow and associated pressure waves in a bubble column reactor is complex, some general design conclusions can be drawn:
  • Increasing the bubble column length will increase the air fraction in the water during bubble flow and will increase both the pressure wave amplitude and propagation speed;
  • Decreasing the bubble column diameter while maintaining the same air flow rate will increase the air fraction in the water during bubble flow and will increase both the pressure wave amplitude and propagation speed;
  • Forcing the rising air bubbles to displace water into a lateral conduit (Figure 1) creates a propagating pressure wave in that conduit with a lower air fraction and a low value of η. This has the effect of slowing the rate of decline in Ax, with increasing distance and maximizing the pressure wave amplitude, as it interacts with the ZVI cartridge (Equation (11)).
The design innovations created in this reactor (Figure 1 and Figure 2), when compared with the 0.24 m3 reactor [10], resulted in increased pressure wave amplitudes, increased wave propagation speeds within the bubble column, and a decreased value of η in the wave intersecting the ZVI cartridge. This, in turn, resulted in a reduced variation in pressure wave amplitude reaching the ZVI and may account for the observed reduced desalination variation (expressed as standard deviation/average value).

7. How Might the Desalination Process Work?

The exact mechanism used by the reactor train to remove Na+ and Cl ions from the water remains unknown. This section considers some of the possibilities, which may account for the ion removal.

7.1. Impact of Air

Desalination involves the removal of Na+ and Cl ions from water. If air is bubbled through water, then the following reactions can be expected:
0.5O2 + H2O = H2O2 = 2OH,
OH + Cl = HClO + e, or [HClO] or H + ClO
The sodium ions would be expected to react as follows:
OH + Na+ = [NaOH]+, or NaO + H+,
These reactions do not remove the Na+ or Cl ions. They will increase pH and change Eh, but will not change EC or ion concentrations measured by a Na+ ISE or Cl ISE. The observed reductions in the measured EC, Na+ ISE, and Cl ISE readings (Table A1 and Table A2) imply that the NaCl product either forms no part of the dissolved matter in the water or if it is still present in the water, it is present in a product form, which has a lower electrical conductivity.
The relative EC ratio of NaCl:Potential soluble products varies with ion concentration [21], but can be approximated as: NaCl:NaClO4 = 1:0.93; NaCl:NaOH = 1:1.98; NaCl:HCl = 1:3.41 [21]. Therefore, the most probable soluble products (NaOH, HCl) are likely to significantly increase rather than decrease EC. The observed position is a decrease in EC.

7.2. Ion Interference

Cl and Na+ ion selective electrodes respond to changes in Cl and Na+ ion concentration. They use an approach based on the Nikolsky–Eisenman equation [22]. This equation [23] is:
E1 = E′ + s(Log[ai + Ki,jpot (a′j)zizj + D]),
where E1 = experimental potential of the solution; E′ = the constant potential of the cell; a1 = the activity (concentration) of the primary ion; Ki,jpot = potentiometric selectivity coefficient; a′j is the known fixed activity of the interfering ion; zi and zj are integers relating to the charge in the primary (i) and interfering ions (j); s = the slope of the line resulting from a plot of E1 (mV) versus ai. [s] varies with temperature in accordance with the relationship, s = 2.303RT/ziF; T = temperature, K; R = gas constant; F = Faraday’s constant. D represents the contribution to potential of all interfering ion species other than species j.
The ISE was calibrated using standard solutions where:
E1 = E′ + s(Log[ai + Ki,potzi),
and no interfering ions are present. Equation (15) demonstrates that it is theoretically possible for an interfering ion to reduce E1 for a specific value of ai, thereby reducing the interpreted ion concentration in the product water.
Given the high ion concentrations of Na+ and Cl ions in the feed water, relative to other ions, it is considered unlikely that the observed decreases in EC or E1, between the feed water and the product water, could be attributed solely to the presence of interfering ions in the water.

7.3. Ion Removal Reaction Order

Most, but not all, ZVI ion removal reactions are first order [24,25,26,27,28,29,30]. The first-order rate constant is defined as:
kobs = −Ln(Ct=n/Ct=0)/tr and Ln(Ct=n) = −kobs tr + Ln(Ct=0)
where kobs = observed rate constant; tr = reaction time; Ct=0 = initial ion concentration; Ct=n = ion concentration at time tr. The first-order rate constant normally remains unchanged with increasing tr. In this study, kobs changed as a function of tr (Table A3, Figure 16). This change in kobs with tr was statistically highly significant.

7.3.1. Interpretation of Statistical Significance

The coefficient of interpretation was used to determine the magnitude of the statistical significance. The coefficient of determination, R2, is the square of the Pearson correlation coefficient (PCC). This coefficient is used to give an indication of the strength of the statistical correlation between two variables, associated with a specific regression trend line [31,32]. By definition, R2 falls within a range 1 to 0 and PCC falls between +1 and −1 [31,32]. The following interpretations are made [31,32]:
  • PCC = 0.9 to 1.0 (R2 = 0.81 to 1.00): Interpretation—very strong correlation
  • PCC = 0.7 to 0.89 (R2 = 0.49 to 0.79): Interpretation—strong correlation
  • PCC = 0.4 to 0.69 (R2 = 0.16 to 0.47): Interpretation—moderate correlation
  • PCC = 0.1 to 0.39 (R2 = 0.01 to 0.15): Interpretation—weak correlation
  • PCC = 0.0 to 0.10 (R2 = 0.00 to 0.01): Interpretation—negligible correlation

7.3.2. Relationship between kobs and Hours on Line

Example relationships between hours on line and kobs are provided in Figure 16a. These relationships demonstrate that kobs decreases as the proportion of an ion removed increases. In these examples, after an initial 5 h operation, (i) the Cl ion removal can be interpreted as a zero-order, first-order or second-order reaction for the example (Figure 16), where Cl at t = 0 = 1 g Cl L−1; (ii) the Cl ion removal can be interpreted as a second-order reaction for the example (Figure 16), where Cl at t = 0 = 3.61 g Cl L−1.
This brief analysis implies that the ion removal reaction has two phases:
  • An initial fast reaction phase;
  • A subsequent slower reaction phase.
This may imply an initial high availability of reactant, followed by a slower reaction associated with the residual declining reactant volume.

7.4. Ion Removal

Each trial batch displayed a pattern of decreasing ion removal rate, with increasing reaction time (e.g., Table A3; Figure 16). Replacement of a water batch with a new water batch resulted in an initial increase in reaction rate, followed by the pattern of decreasing ion removal rate repeating. If kobs is redefined as:
kobs = kins [a[as]b c[Pw]d,
where kins = intrinsic rate constant for as = 1 and Pw = 1; as = surface area (m2 g−1), which is a surrogate for the number of active catalytic sites; Pw = weight of catalyst used; and [a,b,c,d] are constants. In this study, Pw remained constant. Therefore, kobs can be approximated as:
kobs = kins [a[as]b],
The observed power function declines in kobs, with increasing time on line (Table A3) indicate that [b] declines as a function of time from an initial value of 1.0. Linear declines result from [a] declining as a function of time from an initial value of 1.0. Exponential and logarithmic declines result from a combination of [a] and [b] declining as a function of time from an initial value of 1.0.
Oxygenation results in the water pH increasing with time (Table A1) due to the increased presence of OH and OH ions. If [a] and [b] are Eh and pH sensitive or if their decline is due to site blocking by OH and OH ions, then it is reasonable to expect kobs to decrease with increased time on line.
It is also reasonable to expect that replacement of the water in the reactor with a new fresh batch of water will result in desorption of OH and OH ions. This would be expected to increase one or more of [a] and [b].
The data indicate (Figure 17) that the (i) initial value of kobs may increase with decreasing feed water ion concentration; (ii) the initial value of kobs for a water batch can be higher or lower than the initial value of kobs for the previous batch; (iii) the value of kobs measured when freshwater is placed in the reactor can be two to nine magnitudes higher than the measured value of kobs when the water from a batch is removed from the reactor.

7.5. Removed Product

The chemical composition of the removed NaCl product is unknown. The product is not present as a precipitate in the residual catalyst and is not observed as solid micron-sized particles in the product water. Extraction of the catalyst, following the processing of 50 batches (Figure 18), revealed the presence of small (0.7 to 3 microns) colloidal precipitates of an unknown composition in the associated extracted pore water (Figure 18c).
The initial ZVI desalination model [15] assumed that the removed NaCl was incorporated in Fe(OH)x (green rust) and FeOOH akageneite species. This would require the removed NaCl to be accumulated within the catalyst cartridge. This is not the case. Therefore, this model is not applicable. A successor model assumed that the NaCl was adsorbed by Fe(OH)x and FeOOH species, before being redeposited as NaCl within dead-end pores encapsulated by Fe(OH)x and FeOOH species [16]. This desalination model has been verified for pellets constructed from m−Fe0:Fe(b)@n−C0 polymers [33]. Ion removal using this process follows a zero-order reaction [33]. In this study, the catalyst bed was constructed with open pores (Figure 4 and Figure 18) from uncemented particles, where <5% of the total porosity was contained in dead-end pores. The resultant ion removal had an initial fast removal stage followed by removal using a second-order reaction (Figure 16). Therefore, the desalination model applied to the pellets [33] is not applicable to this study.
In this study, the redox regime (Table A1, Figure A4) is in (Hydrogen Stability Zone 0), as defined in reference [34]. In this redox regime, the dominant Fe oxidation number is 3+ [35]. Therefore, the catalytic polymer coating used (Figure 4 and Figure 18) will be either SiO2@Fe(a,b,c)@urea, SiO2@Fe(a,b,c):urea, or SiO2@Fe(a,b,c). If the polymers had been operated and constructed in (Hydrogen Stability Zones 1 or 2 [34]), the constructed polymer would have been one of: SiO2@n−Fe0@urea, SiO2@n−Fe0:urea, SiO2@n−Fe0:Fe(a,b,c)@urea, SiO2@n−Fe0:Fe(a,b,c):urea, or SiO2@n−Fe0:Fe(a,b,c). The symbol (:) indicates chemical bonding. The symbol (@) indicates physical bonding.
Entrained n−Fe0:Fe(a,b,c) polymers were used to remove Na+ and Cl ions from water [11]. Ion removal using these polymers is characterized by a situation, where kobs decreases with increasing tr [11]. This is the observed position (Table A3, Figure 17).

7.5.1. Interpreted Catalytic Model

The CPSAD operating model [10,33] used in this study uses pressure increases to adsorb reactants onto the catalytic site and pressure decreases to desorb products from the catalytic site [10,33]. The amorphous, hydrated, precipitated products in Figure 18c are consistent with adsorption and desorption of a cation ion (M) from the water (e.g., Fe3+, Ca2+, Mg2+), where the catalyst surface has a basic reaction during pressure increases of the form:
Surface−OH + OH = Surface−O + H2O,
Surface−O + M2+ + OH = Surface−OMOH,
Surface−O + M3+ + 2OH = Surface−OM(OH)2,
Surface−O + M+ = Surface−OM,
and a basic reaction during pressure decreases of the form:
Surface−OM(OH)2 + OH = Surface−O + M(OH)3,
Surface−OM + OH = Surface−O + MOH,
Surface−OMOH + OH = Surface−O + M(OH)2,
nM(OH)x = [M(OH)x]n,
The desorbed n−M(a,b,c) polymer ([M(OH)x]n) or another metal polymer then aggregates within the water body and removes Na+ and Cl ions from the water [11]. Under this model, kobs would be expected to decline, as
  • The availability of the metal (M) cation in the water decreased towards an equilibrium level.
  • The growing n−M(a,b,c) polymer became saturated with Na+ and Cl ions.
Removal of product water containing the n−M(a,b,c) polymer would result in desalination continuing over a long period of time, outside the reaction environment. This is the observed position (Table A2).
Replacement of the water, following partial desalination, will result in the reaction:
Surface − OM(OH)2 + H2O = Surface − OH + M(OH)3,
Surface − OM + H2O = Surface − OH + MOH,
Surface − OMOH + H2O = Surface − OH + M(OH)2,
Changing the water batch results in regeneration of the catalyst. Under this model, the function of the catalyst is to catalyze the formation of entrained nano−metal polymers, which act as scavenging agents for the removal of Na+ and Cl ions from the water.

7.5.2. Mechanism of Na+ and Cl Ion Removal

In this study, the competing ions for Na+ and Cl ions in the metal polymers were OH, SO42−, HCO3, and CO32− ions. Examination of the product water identified the presence of entrained, amorphous, hydrated spheres or ellipsoids of metal polymer, with a maximum dimension in a range 0.01 to 0.15 mm (e.g., Figure 19). These polymers may have formed from the aggregation of the polymers ejected from the catalyst body (Figure 18c). Their basic structure is a metal polymer surrounding a fluid-filled cavity (Figure 19). This type of polymer structure is an active pollutant scavenger [11]. The fluid-filled cavity is effectively a mobile, dead-end pore, which grows by polymer particle aggregation and accretion [11]. The fluid contents of the pore are isolated from the water body and will not register on EC or ISE instruments analyzing the water body.
An earlier study [33] established that m−Fe0:n−Fe(b)@n−C0 bounded dead-end pores can actively concentrate both Na+ and Cl ions (and their products) within the dead-end pore. It is possible that the amorphous, hydrated spheres or ellipsoids of metal polymer (Figure 19) accumulated and concentrated the removed Na+ and Cl ions recorded in this study.
This explanation would account for why, in some product water batches placed in storage for 2 years, there was an increase in Na+ ion concentration after 2 years storage (Table A2).

7.6. Current Working Model

The current working model is:
  • The SiO2@Fe(a,b,c)@urea polymer in a pressure−swing environment acts as a catalyst for the formation and release of n−metal polymers, utilizing dissolved metal cations (e.g., Fen+; Ca2+, Mg2+, Cun+) within the water to form hydrated metal hydroxy-oxide polymers;
  • These metal polymers coalesce to form entrained, fluid-filled colloids;
  • The fluid-filled colloids extract Na+ and Cl ions from the water and deposit them within the fluid-filled cavities. The fluid-filled cavities may also extract OH, HCO3, SO4, NOx, H+ ions, and other anions/cations from the water [11]. Since the data (e.g., Table A3) indicate that the Na+ and Cl ions are removed by different processes, it may be that some fluid cavities only contain (or are dominated by) Na+ species, while others only contain (or are dominated by) Cl species.
If this working model is correct, then it is likely that the reactor could be repurposed to remediate most waters containing metal cations. The cations would be used to either form entrained, scavenging polymers within the water body or would be extracted and added to a polluted water body to affect remediation. An example process flow diagram showing how this concept could potentially be used to remediate water is provided in Figure 20.

8. Conclusions

This study established that a batch flow, bubble column, diffusion reactor train, incorporating recycle, can be operated as a catalytic pressure swing adsorption/desorption (separation) reactor (CPSAD) to partially desalinate water. The catalyst used was a SiO2@n−Fe(a,b,c)@urea and SiO2@n−Fe(a,b,c) polymer.
An individual catalyst charge was demonstrated to process >50 batches (43 m3) without loss of activity. The catalyst is structured such that replacement of one batch of processed water with a fresh batch of water automatically results in regeneration of the catalyst. The catalyst is used to produce entrained n−metal(a,b,c) polymers in the water. The metals used are derived from cations contained within the feed water. These entrained polymers aggregate and accumulate within the water body to form entrained polymers colloids, which encapsulate fluid-filled pores. These polymers act as scavengers, to remove Na+ and Cl ions from the water. During reactor operation, the amount of desalination associated with each batch of water increases for a period of 3 to 30 h before the salinity stabilizes at a new lower level.
The trials identified that the process produces partially desalinated water, no reject brine, and that the water product contained some entrained nano/micron-sized colloidal flocculates. The process requires access to an energy source. The principal process operating cost is the cost and provision of compressed air, combined with costs associated with the addition and removal of water from the reactor. The expected process operating cost will vary with geographic location.
This study used an artificial saline water as a feedstock. The next stages in the development of this process will require operation using a natural saline water.

Funding

This research received no external funding.

Data Availability Statement

The data used in this study are placed in the tables and figures contained within the paper and in Appendix A.

Acknowledgments

The Editorial team and the three Reviewers are thanked for their constructive comments.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Quantitative Data

Table A1. Basic parameters associated with each trial. Fw = feed water; Ppw = product water; T = temperature.
Table A1. Basic parameters associated with each trial. Fw = feed water; Ppw = product water; T = temperature.
FwPpwFwPpwFwPpwAir Flow RateAir Flow DurationO2 SuppliedO2 Adsorbed
TrialpHpHEh, mVEh, mVT, °CT, °CL/hrhrsL/Lg/L
15.856.2646039716.916.5658.95 0.140.04
26.307.1942043715.714.718045.22 9.491.86
36.287.1744142215.315.518045.782.010.65
46.307.3753250115.015.66056.930.830.38
56.568.7543838516.615.06061.670.900.41
65.986.2846439512.212.56017.650.26nd
76.688.0047137412.010.86064.470.94nd
86.577.504562857.910.16036.300.53nd
96.687.574924558.17.46038.020.56nd
106.266.904493348.77.56044.000.64nd
116.126.504643665.94.36030.170.44nd
126.096.224403556.52.76038.130.56nd
135.586.244553952.94.26037.530.55nd
145.546.054654306.31.96022.400.33nd
155.666.174734445.65.16030.780.45nd
165.907.015154444.66.36033.430.49nd
176.117.305074562.25.46030.700.45nd
186.407.254994374.15.76030.780.45nd
196.407.204893745.85.06029.780.44nd
205.616.084543554.32.96029.880.440.12
215.736.324373573.85.96034.220.500.11
225.906.684333655.45.96043.120.630.16
235.826.315114627.07.56045.230.660.21
246.166.224964348.38.46030.420.450.09
255.375.874793838.110.46030.850.450.16
265.316.4046743110.811.76039.450.580.15
275.976.2045145611.818.26038.850.570.24
285.976.4545346914.016.06029.650.430.21
295.756.4653542912.412.06029.270.430.20
306.156.6648346814.014.76034.720.510.19
317.409.3051947317.812.66036.480.530.23
327.027.1754243214.815.36039.380.580.22
335.346.4048045617.016.46039.830.580.18
345.406.2235423414.514.66039.600.580.28
355.626.0439338014.313.96032.530.480.29
365.446.1449543115.413.26032.900.480.30
375.796.5847341111.011.36031.380.460.30
385.926.7046342412.112.56032.000.470.31
396.176.6841344412.910.66031.630.460.30
406.266.765054229.610.86032.820.480.37
415.976.7546741512.110.96031.980.470.24
426.126.604454268.18.16030.800.450.28
436.046.304584164.95.36031.820.470.29
445.946.223543364.23.06030.020.440.26
456.006.304804223.65.06031.130.460.26
465.575.914854637.13.46030.430.450.24
475.484.265024582.416.36028.600.420.21
484.404.2649745315.813.76035.550.520.25
494.725.2153250317.318.56031.120.460.27
504.905.0854764217.818.76015.630.23nd
Table A2. Salinity data and ion concentrations associated with each trial. Fw = feed water; Ppw = product water.
Table A2. Salinity data and ion concentrations associated with each trial. Fw = feed water; Ppw = product water.
FwPpwPpw, after 2 Years StorageFwPpwPpw, after 2 Years StorageFw, ECPpw, ECFwPpwPpw, after 2 Years Storage
TrialCl g/LCl g/LCl g/LNa+ g/LNa+ g/LNa+ g/LmScm−1mScm−1NaCl, g/LNaCl, g/LNaCl, g/L
10.880.720.570.430.210.162.161.901.310.930.73
21.311.000.710.680.500.163.072.581.991.500.87
31.951.360.661.090.770.224.533.463.042.130.88
41.000.730.650.450.220.572.212.001.450.951.22
51.531.190.270.820.680.593.813.102.351.870.86
61.810.920.770.950.500.633.932.422.761.421.40
71.070.69 0.540.21 2.491.711.610.90
80.860.53 0.450.17 2.071.231.310.70
91.530.790.601.000.650.533.061.712.531.441.13
101.350.750.060.550.370.652.871.771.901.120.71
111.550.83 0.830.51 3.161.892.381.34
121.671.250.910.940.650.622.541.922.611.901.53
131.400.730.220.790.400.852.021.192.191.131.07
140.560.350.470.440.240.350.880.631.000.590.82
151.190.940.550.630.500.632.031.581.821.441.18
163.201.350.751.220.700.314.882.724.422.051.06
173.671.750.601.240.580.625.302.344.912.331.22
183.611.250.651.320.580.695.312.204.931.831.34
193.581.850.031.320.660.805.732.214.902.510.83
204.352.190.811.490.660.806.622.245.842.851.61
211.720.600.520.820.510.543.421.412.541.111.06
223.540.910.532.050.790.446.282.245.591.700.97
232.962.360.611.710.700.594.952.034.673.061.20
240.980.64 0.430.22 1.491.001.410.86
250.770.67 0.320.21 1.220.841.090.88
263.802.470.992.621.210.979.114.386.423.681.96
271.640.680.330.710.450.492.621.502.351.130.82
281.150.790.330.520.360.492.021.411.671.150.82
290.850.69 0.360.32 1.491.041.211.01
301.320.85 0.610.40 1.771.131.931.25
315.031.780.642.220.900.727.703.247.252.681.36
325.392.190.542.720.830.619.463.618.113.021.15
335.311.440.692.150.650.628.182.897.462.091.31
345.361.680.532.150.720.718.183.137.512.401.24
350.870.730.270.420.260.451.771.191.290.990.72
361.710.910.340.580.470.302.561.872.291.380.64
372.020.930.880.580.350.462.021.502.601.281.34
380.880.670.560.350.290.511.641.241.230.961.07
392.841.630.221.690.830.427.023.124.532.460.64
404.672.070.421.831.200.587.013.496.503.271.00
411.860.670.081.471.160.251.681.313.331.830.33
424.151.800.112.721.450.126.623.486.873.250.23
435.522.710.283.452.130.437.415.068.974.840.71
442.862.12 2.361.63 5.753.555.223.75
451.601.040.910.730.330.381.821.412.331.371.29
463.922.350.701.271.270.494.613.745.193.621.19
475.635.370.532.502.310.336.936.308.137.680.86
481.801.570.300.360.200.451.721.482.161.770.75
492.871.67 1.230.54 3.662.164.102.21
505.391.32 0.290.18 1.801.315.681.50
Table A3. Desalination rate constant decline analysis, measured values of kinsf, p, and R2 associated with ion removal. The values measured from a regression plot of kobsf vs. reaction time (a (hours)). kobsf = (−Ln (Ct=n/Ct=0))/(t = n). Ct=0 = ion concentration at t = 0. Ct=n = ion concentration at t = n. (t = n) = time spent in the reaction environment, seconds = 3600a. Decline curve type = power function; exponential (exp) function, natural logarithmic function; linear function. Power relationship: kobsf = kinsf (a)p; linear relationship: kobsf = kinsf (a) + p; logarithmic relationship: kobsf = kinsf × Ln(a · p) exponential relationship: kobsf = kinsf × e(a · p).
Table A3. Desalination rate constant decline analysis, measured values of kinsf, p, and R2 associated with ion removal. The values measured from a regression plot of kobsf vs. reaction time (a (hours)). kobsf = (−Ln (Ct=n/Ct=0))/(t = n). Ct=0 = ion concentration at t = 0. Ct=n = ion concentration at t = n. (t = n) = time spent in the reaction environment, seconds = 3600a. Decline curve type = power function; exponential (exp) function, natural logarithmic function; linear function. Power relationship: kobsf = kinsf (a)p; linear relationship: kobsf = kinsf (a) + p; logarithmic relationship: kobsf = kinsf × Ln(a · p) exponential relationship: kobsf = kinsf × e(a · p).
Cl Ions Na+ Ions
TrialkinstpR2TypekinstpR2Type
10.000000810−1.50847443098.77%Power0.00000108−0.9510182295.77%Power
20.000074210−0.99527002097.30%Power−0.000005440.0000227697.90%Log
30.000097090−0.91913646098.28%Power0.00007177−0.8456168897.52%Power
40.000045530−0.88285049094.92%Power0.00005379−0.7352932590.94%Power
50.000040590−0.86423844095.29%Power0.00002852−0.7938856284.86%Power
6−0.0000015700.00003660096.03%Linear−0.000020510.0000692198.76%Log
7−0.0000056600.00002406095.62%Log−0.000007990.0000349491.96%Log
8−0.0000078000.00002971096.24%Log−0.000011610.0000493698.32%Log
9−0.0000091500.00003832096.60%Log0.00020484−1.2442828595.32%Power
100.000063540−0.72902784096.83%Power0.00003571−0.7761714177.44%Power
11−0.0000086100.00003470094.78%Log−0.000005870.0000259080.95%Log
12−0.0000045000.00001939095.55%Log−0.000002560.0000124692.25%Log
130.000072520−0.71898863098.03%Power−0.000008720.0000354197.58%Log
14−0.0000080600.00003072096.30%Log−0.000014360.0000518097.55%Log
15−0.0000015600.00000884084.35%Log−0.000001820.0000092171.07%Log
16−0.0000046800.00002314091.65%Log0.00004299−0.6626879276.81%Power
17−0.0000052500.00002602078.44%Log−0.000004420.0000208281.83%Log
180.000036130−0.46144223090.19%Power0.00008015−0.7206462295.24%Power
190.000047630−0.57070204080.65%Power0.00006918−0.6755619391.71%Power
20−0.0000165400.00006108098.44%Log0.984370590.0000288290.45%Linear
21−0.0000059100.00002632044.15%Log−0.000000380.0000173189.05%Linear
22−0.0000067900.00002998075.19%Log−0.000011740.0000492798.22%Log
23−0.0000088700.00003724090.96%Log−0.000008560.0000367394.73%Log
240.000044180−0.67868612076.48%Power0.00002140−0.0432160085.40%Exp
25−0.0000057400.00002310074.60%Log−0.000003580.0000156491.45%Log
260.000005670−0.01120830040.98%Exp0.00001068−0.0196615985.11%Exp
270.000073350−0.80061460090.64%Power−0.000003280.0000154984.38%Log
280.000010380−0.03758638085.18%Exp0.00005004−0.8214323289.53%Power
29−0.0000024100.00000917074.19%Log−0.000001680.0000052758.43%Log
30−0.0000053700.00002285093.01%Log−0.000005740.0000240089.76%Log
31−0.0000144700.00005806092.46%Log−0.000009020.0000383597.75%Log
32−0.0000099100.00004189086.77%Log−0.000008160.0000380094.00%Log
33−0.0000066900.00003353094.62%Log−0.000006710.0000335691.82%Log
34−0.0000032300.00002209066.75%Log−0.000007670.0000358495.11%Log
35−0.0000028900.00001106071.92%Log0.00002031−0.4535257876.38%Power
36−0.0000014000.00001040028.35%Log−0.000001660.0000087837.38%Log
37−0.0000073100.00003259090.22%Log−0.000005240.0000221964.41%Log
380.000102540−1.11357014097.57%Power−0.000002470.0000121352.54%Log
39−0.0000196100.00006939075.14%Log−0.000006340.0000282494.29%Log
400.000018310−0.22999245046.14%Power−0.000000150.000009648.67%Linear
41−0.0000004000.00002119093.82%Linear−0.000000060.000006367.50%Linear
420.0000027200.00000002058.75%Log0.00001162−0.0126281234.21%Exp
43−0.0000021100.00001399055.08%Log−0.000002430.0000112649.15%Log
44−0.0000016800.00000842067.42%Log0.00009117−0.9143337986.07%Power
450.000005280−0.00628353010.39%Exp0.00001331−0.0157384715.46%Exp
46−0.0000064500.00002440080.66%Log−0.000004920.0000187548.57%Log
47−0.0000012500.00000476076.29%Log−0.000002900.0000114667.36%Log
48−0.0000021400.00000871068.31%Log−0.000001910.0000142264.55%Log
49−0.0000005700.00002339073.67%Linear−0.000000300.0000153478.01%Linear
50−0.0000253100.00009357095.73%Log−0.000023960.0000887197.22%Log
Table A4. Indicative relationship between crop type, irrigation water salinity, and crop yield relative to crop yield when irrigated with freshwater.
Table A4. Indicative relationship between crop type, irrigation water salinity, and crop yield relative to crop yield when irrigated with freshwater.
Constants Irrigation Water Salinity, g L−1
Cropabc1234579
Almond, Plum0.72840.13370.318.7%0.0%0.0%0.0%0.0%0.0%0.0%
Apricot0.90520.17480.30.0%0.0%0.0%0.0%0.0%0.0%0.0%
Avocado, Blackberry, Boysenberry0.88980.16870.31.2%0.0%0.0%0.0%0.0%0.0%0.0%
Barley0.10000.40000.3100.0%100.0%100.0%88.0%75.0%49.0%23.0%
Beans0.7260.08550.314.2%0.0%0.0%0.0%0.0%0.0%0.0%
Beets0.36320.18280.371.1%23.8%0.0%0.0%0.0%0.0%0.0%
Broadbean0.38400.07610.357.7%7.8%0.0%0.0%0.0%0.0%0.0%
Broccoli0.37400.13210.364.6%16.0%0.0%0.0%0.0%0.0%0.0%
Cabbage0.39120.08900.358.0%7.2%0.0%0.0%0.0%0.0%0.0%
Cantaloupe0.29040.08060.370.3%32.6%0.0%0.0%0.0%0.0%0.0%
Carrot0.57040.07310.333.2%0.0%0.0%0.0%0.0%0.0%0.0%
Cotton0.10420.40460.3100.0%100.0%99.8%86.3%72.7%45.6%18.5%
Corn0.48080.10230.347.7%0.0%0.0%0.0%0.0%0.0%0.0%
Cowpea0.55480.08970.336.8%0.0%0.0%0.0%0.0%0.0%0.0%
Cucumber0.53440.16910.347.4%0.0%0.0%0.0%0.0%0.0%0.0%
Date Palm0.14280.07070.388.5%69.9%51.4%32.8%14.3%0.0%0.0%
Fig, Olive, Pomegranate0.35320.11870.366.0%20.0%0.0%0.0%0.0%0.0%0.0%
Flax0.48080.10230.347.7%0.0%0.0%0.0%0.0%0.0%0.0%
Grape0.38040.06990.357.5%8.1%0.0%0.0%0.0%0.0%0.0%
Grapefruit, Orange, Lemon, Apple, Pear, Walnut0.64520.14570.330.7%0.0%0.0%0.0%0.0%0.0%0.0%
Groundnut1.20760.48530.30.0%0.0%0.0%0.0%0.0%0.0%0.0%
Lettuce0.51920.08510.341.0%0.0%0.0%0.0%0.0%0.0%0.0%
Onion0.63520.09440.326.9%0.0%0.0%0.0%0.0%0.0%0.0%
Peach0.82440.17160.310.0%0.0%0.0%0.0%0.0%0.0%0.0%
Pepper0.57040.10880.336.7%0.0%0.0%0.0%0.0%0.0%0.0%
Potato0.48080.10230.347.7%0.0%0.0%0.0%0.0%0.0%0.0%
Radish0.51400.07590.340.8%0.0%0.0%0.0%0.0%0.0%0.0%
Raspberry0.88560.10730.30.0%0.0%0.0%0.0%0.0%0.0%0.0%
Rice0.46960.17430.356.4%0.0%0.0%0.0%0.0%0.0%0.0%
Safflower0.43440.28740.372.3%15.8%0.0%0.0%0.0%0.0%0.0%
Sesbania0.28160.08120.371.5%34.9%0.0%0.0%0.0%0.0%0.0%
Soybean0.79840.49770.346.0%0.0%0.0%0.0%0.0%0.0%0.0%
Sorghum0.28200.13420.376.8%40.1%3.4%0.0%0.0%0.0%0.0%
Spinach0.30760.07790.367.8%27.8%0.0%0.0%0.0%0.0%0.0%
Strawberry1.33700.16930.30.0%0.0%0.0%0.0%0.0%0.0%0.0%
Sugar Beet0.11800.40460.3100.0%100.0%94.4%79.1%63.8%33.1%2.4%
Sweet Corn0.48080.10230.347.7%0.0%0.0%0.0%0.0%0.0%0.0%
Sweet Potato0.44480.08420.350.6%0.0%0.0%0.0%0.0%0.0%0.0%
Tomato0.39920.12540.360.6%8.7%0.0%0.0%0.0%0.0%0.0%
Wheat0.14280.42860.3100.0%100.0%87.2%68.6%50.0%12.9%0.0%
Alfalfa0.29680.07530.368.9%30.4%0.0%0.0%0.0%0.0%0.0%
Barley (Hay)0.06960.21430.3100.0%100.0%94.3%85.2%76.2%58.1%40.0%
Bermuda Grass0.25680.22190.388.8%55.4%22.0%0.0%0.0%0.0%0.0%
Clover, berseem0.76800.22010.322.2%0.0%0.0%0.0%0.0%0.0%0.0%
Clover, alsike, ladino, red, strawberry0.46960.08630.347.6%0.0%0.0%0.0%0.0%0.0%0.0%
Corn (forage)0.23840.08030.377.0%46.0%15.1%0.0%0.0%0.0%0.0%
Crested Wheat Grass0.16000.07030.386.2%65.4%44.6%23.8%3.0%0.0%0.0%
Harding Grass0.29880.16980.378.1%39.3%0.4%0.0%0.0%0.0%0.0%
Lovegrass, Meadow Foxtail0.33320.08340.365.0%21.7%0.0%0.0%0.0%0.0%0.0%
Orchard Grass0.29680.07530.368.9%30.4%0.0%0.0%0.0%0.0%0.0%
Perennial Rye Grass0.29840.20740.381.9%43.2%4.4%0.0%0.0%0.0%0.0%
Sudan Grass0.17240.06020.383.6%61.2%38.8%16.4%0.0%0.0%0.0%
Tall Fescue0.20960.10110.382.9%55.6%28.4%1.1%0.0%0.0%0.0%
Tall Wheat Grass0.27480.25800.390.1%54.4%18.6%0.0%0.0%0.0%0.0%
Trefoil, big0.76800.22010.322.2%0.0%0.0%0.0%0.0%0.0%0.0%
Trefoil, Birdsfoot narrow leaf0.40000.25000.373.0%21.0%0.0%0.0%0.0%0.0%0.0%
Vetch0.44400.16790.359.1%1.3%0.0%0.0%0.0%0.0%0.0%
Wheat Grass0.16640.15480.393.8%72.2%50.6%29.0%7.3%0.0%0.0%
Wild Rye, beardless0.23840.08030.377.0%46.0%15.1%0.0%0.0%0.0%0.0%
Table A5. Expected crop yield, following partial desalination of irrigation water, when compared to irrigation with freshwater.
Table A5. Expected crop yield, following partial desalination of irrigation water, when compared to irrigation with freshwater.
Constants Salinity of Water Entering the Partial Desalination Reactor, g L−1
Cropab(1 + c)(1 − e)1234579
Almond, Plum0.728400.133700.7856.6%0.0%0.0%0.0%0.0%0.0%0.0%
Apricot0.905200.174800.7846.9%0.0%0.0%0.0%0.0%0.0%0.0%
Avocado, Blackberry, Boysenberry0.889800.168700.7847.5%0.0%0.0%0.0%0.0%0.0%0.0%
Barley0.100000.400000.78100.0%100.0%100.0%100.0%100.0%85.4%69.8%
Beans0.726000.085500.7851.9%0.0%0.0%0.0%0.0%0.0%0.0%
Beets0.363200.182800.7890.0%61.6%33.3%5.0%0.0%0.0%0.0%
Broadbean0.384000.076100.7877.7%47.7%17.8%0.0%0.0%0.0%0.0%
Broccoli0.374000.132100.7884.0%54.9%25.7%0.0%0.0%0.0%0.0%
Cabbage0.391200.089000.7878.4%47.9%17.4%0.0%0.0%0.0%0.0%
Cantaloupe0.290400.080600.7885.4%62.8%40.1%17.5%0.0%0.0%0.0%
Carrot0.570400.073100.7862.8%18.3%0.0%0.0%0.0%0.0%0.0%
Cotton0.104200.404600.78100.0%100.0%100.0%100.0%99.8%83.6%67.3%
Corn0.480800.102300.7872.7%35.2%0.0%0.0%0.0%0.0%0.0%
Cowpea0.554800.089700.7865.7%22.4%0.0%0.0%0.0%0.0%0.0%
Cucumber0.534400.169100.7875.2%33.5%0.0%0.0%0.0%0.0%0.0%
Date Palm0.142800.070700.7895.9%84.8%73.7%62.5%51.4%29.1%6.8%
Fig, Olive, Pomegranate0.353200.118700.7884.3%56.8%29.2%1.7%0.0%0.0%0.0%
Flax0.480800.102300.7872.7%35.2%0.0%0.0%0.0%0.0%0.0%
Grape 0.380400.069900.7877.3%47.6%18.0%0.0%0.0%0.0%0.0%
Grapefruit, Orange, Lemon, Apple, Pear, Walnut0.645200.145700.7864.2%13.9%0.0%0.0%0.0%0.0%0.0%
Groundnut1.20760.48530.7854.3%0.0%0.0%0.0%0.0%0.0%0.0%
Lettuce0.51920.08510.7868.0%27.5%0.0%0.0%0.0%0.0%0.0%
Onion0.63520.09440.7859.9%10.3%0.0%0.0%0.0%0.0%0.0%
Peach0.82440.17160.7852.9%0.0%0.0%0.0%0.0%0.0%0.0%
Pepper0.57040.10880.7866.4%21.9%0.0%0.0%0.0%0.0%0.0%
Potato0.48080.10230.7872.7%35.2%0.0%0.0%0.0%0.0%0.0%
Radish0.51400.07590.7867.5%27.4%0.0%0.0%0.0%0.0%0.0%
Raspberry0.88560.10730.7841.7%0.0%0.0%0.0%0.0%0.0%0.0%
Rice0.46960.17430.7880.8%44.2%7.5%0.0%0.0%0.0%0.0%
Safflower0.43440.28740.7894.9%61.0%27.1%0.0%0.0%0.0%0.0%
Sesbania0.28160.08120.7886.2%64.2%42.2%20.3%0.0%0.0%0.0%
Soybean0.79840.49770.7887.5%25.2%0.0%0.0%0.0%0.0%0.0%
Sorghum0.28200.13420.7891.4%69.4%47.4%25.4%3.4%0.0%0.0%
Spinach0.30760.07790.7883.8%59.8%35.8%11.8%0.0%0.0%0.0%
Strawberry1.33700.16930.7812.6%0.0%0.0%0.0%0.0%0.0%0.0%
Sugar Beet0.11800.40460.78100.0%100.0%100.0%100.0%94.4%76.0%57.6%
Sweet Corn0.48080.10230.7872.7%35.2%0.0%0.0%0.0%0.0%0.0%
Sweet Potato0.44480.08420.7873.7%39.0%4.3%0.0%0.0%0.0%0.0%
Tomato0.39920.12540.7881.4%50.3%19.1%0.0%0.0%0.0%0.0%
Wheat0.14280.42860.78100.0%100.0%100.0%98.3%87.2%64.9%42.6%
Alfalfa0.29680.07530.7884.4%61.2%38.1%14.9%0.0%0.0%0.0%
Barley (Hay)0.06960.21430.78100.0%100.0%100.0%99.7%94.3%83.4%72.6%
Bermuda Grass0.25680.22190.78100.0%82.1%62.1%42.1%22.0%0.0%0.0%
Clover, berseem0.76800.22010.7862.1%2.2%0.0%0.0%0.0%0.0%0.0%
Clover, alsike, ladino, red, strawberry0.46960.08630.7872.0%35.4%0.0%0.0%0.0%0.0%0.0%
Corn (forage)0.23840.08030.7889.4%70.8%52.2%33.6%15.1%0.0%0.0%
Crested Wheat Grass0.16000.07030.7894.6%82.1%69.6%57.1%44.6%19.7%0.0%
Harding Grass0.29880.16980.7893.7%70.4%47.1%23.8%0.4%0.0%0.0%
Lovegrass, Meadow Foxtail0.33320.08340.7882.4%56.4%30.4%4.4%0.0%0.0%0.0%
Orchard Grass0.29680.07530.7884.4%61.2%38.1%14.9%0.0%0.0%0.0%
Perennial Rye Grass0.29840.20740.7897.5%74.2%50.9%27.6%4.4%0.0%0.0%
Sudan Grass0.17240.06020.7892.6%79.1%65.7%52.2%38.8%11.9%0.0%
Tall Fescue0.20960.10110.7893.8%77.4%61.1%44.7%28.4%0.0%0.0%
Tall Wheat Grass0.27480.25800.78100.0%82.9%61.5%40.1%18.6%0.0%0.0%
Trefoil, big0.76800.22010.7862.1%2.2%0.0%0.0%0.0%0.0%0.0%
Trefoil, Birdsfoot narrow leaf0.40000.25000.7893.8%62.6%31.4%0.2%0.0%0.0%0.0%
Vetch0.44400.16790.7882.2%47.5%12.9%0.0%0.0%0.0%0.0%
Wheat Grass0.16640.15480.78100.0%89.5%76.5%63.6%50.6%24.6%0.0%
Wild Rye, beardless0.23840.08030.7889.4%70.8%52.2%33.6%15.1%0.0%0.0%
Table A6. Proportion of expected crop yield with freshwater irrigation (Table A4), added by irri-gating with partially desalinated water. Calculated as Results in Table A5—Results in Table A4.
Table A6. Proportion of expected crop yield with freshwater irrigation (Table A4), added by irri-gating with partially desalinated water. Calculated as Results in Table A5—Results in Table A4.
9Constants Salinity of Water Entering the Partial Desalination Reactor, g L−1
Cropab(1 + c)(1 − e)1234579
Almond, Plum0.728400.133700.7837.9%0.0%0.0%0.0%0.0%0.0%0.0%
Apricot0.905200.174800.7846.9%0.0%0.0%0.0%0.0%0.0%0.0%
Avocado, Blackberry, Boysenberry0.889800.168700.7846.3%0.0%0.0%0.0%0.0%0.0%0.0%
Barley0.100000.400000.780.0%0.0%0.0%12.0%25.0%36.4%46.8%
Beans0.726000.085500.7837.8%0.0%0.0%0.0%0.0%0.0%0.0%
Beets0.363200.182800.7818.9%37.8%33.3%5.0%0.0%0.0%0.0%
Broadbean0.384000.076100.7820.0%39.9%17.8%0.0%0.0%0.0%0.0%
Broccoli0.374000.132100.7819.4%38.9%25.7%0.0%0.0%0.0%0.0%
Cabbage0.391200.089000.7820.3%40.7%17.4%0.0%0.0%0.0%0.0%
Cantaloupe0.290400.080600.7815.1%30.2%40.1%17.5%0.0%0.0%0.0%
Carrot0.570400.073100.7829.7%18.3%0.0%0.0%0.0%0.0%0.0%
Cotton0.104200.404600.780.0%0.0%0.2%13.7%27.1%37.9%48.8%
Corn0.480800.102300.7825.0%35.2%0.0%0.0%0.0%0.0%0.0%
Cowpea0.554800.089700.7828.8%22.4%0.0%0.0%0.0%0.0%0.0%
Cucumber0.534400.169100.7827.8%33.5%0.0%0.0%0.0%0.0%0.0%
Date Palm0.142800.070700.787.4%14.9%22.3%29.7%37.1%29.1%6.8%
Fig, Olive, Pomegranate0.353200.118700.7818.4%36.7%29.2%1.7%0.0%0.0%0.0%
Flax0.480800.102300.7825.0%35.2%0.0%0.0%0.0%0.0%0.0%
Grape 0.380400.069900.7819.8%39.6%18.0%0.0%0.0%0.0%0.0%
Grapefruit, Orange, Lemon, Apple, Pear, Walnut0.645200.145700.7833.6%13.9%0.0%0.0%0.0%0.0%0.0%
Groundnut1.20760.48530.7854.3%0.0%0.0%0.0%0.0%0.0%0.0%
Lettuce0.51920.08510.7827.0%27.5%0.0%0.0%0.0%0.0%0.0%
Onion0.63520.09440.7833.0%10.3%0.0%0.0%0.0%0.0%0.0%
Peach0.82440.17160.7842.9%0.0%0.0%0.0%0.0%0.0%0.0%
Pepper0.57040.10880.7829.7%21.9%0.0%0.0%0.0%0.0%0.0%
Potato0.48080.10230.7825.0%35.2%0.0%0.0%0.0%0.0%0.0%
Radish0.51400.07590.7826.7%27.4%0.0%0.0%0.0%0.0%0.0%
Raspberry0.88560.10730.7841.7%0.0%0.0%0.0%0.0%0.0%0.0%
Rice0.46960.17430.7824.4%44.2%7.5%0.0%0.0%0.0%0.0%
Safflower0.43440.28740.7822.6%45.2%27.1%0.0%0.0%0.0%0.0%
Sesbania0.28160.08120.7814.6%29.3%42.2%20.3%0.0%0.0%0.0%
Soybean0.79840.49770.7841.5%25.2%0.0%0.0%0.0%0.0%0.0%
Sorghum0.28200.13420.7814.7%29.3%44.0%25.4%3.4%0.0%0.0%
Spinach0.30760.07790.7816.0%32.0%35.8%11.8%0.0%0.0%0.0%
Strawberry1.33700.16930.7812.6%0.0%0.0%0.0%0.0%0.0%0.0%
Sugar Beet0.11800.40460.780.0%0.0%5.6%20.9%30.7%43.0%55.2%
Sweet Corn0.48080.10230.7825.0%35.2%0.0%0.0%0.0%0.0%0.0%
Sweet Potato0.44480.08420.7823.1%39.0%4.3%0.0%0.0%0.0%0.0%
Tomato0.39920.12540.7820.8%41.5%19.1%0.0%0.0%0.0%0.0%
Wheat0.14280.42860.780.0%0.0%12.8%29.7%37.1%52.0%42.6%
Alfalfa0.29680.07530.7815.4%30.9%38.1%14.9%0.0%0.0%0.0%
Barley (Hay)0.06960.21430.780.0%0.0%5.7%14.5%18.1%25.3%32.6%
Bermuda Grass0.25680.22190.7811.2%26.7%40.1%42.1%22.0%0.0%0.0%
Clover, berseem0.76800.22010.7839.9%2.2%0.0%0.0%0.0%0.0%0.0%
Clover, alsike, ladino, red, strawberry0.46960.08630.7824.4%35.4%0.0%0.0%0.0%0.0%0.0%
Corn (forage)0.23840.08030.7812.4%24.8%37.2%33.6%15.1%0.0%0.0%
Crested Wheat Grass0.16000.07030.788.3%16.6%25.0%33.3%41.6%19.7%0.0%
Harding Grass0.29880.16980.7815.5%31.1%46.6%23.8%0.4%0.0%0.0%
Lovegrass, Meadow Foxtail0.33320.08340.7817.3%34.7%30.4%4.4%0.0%0.0%0.0%
Orchard Grass0.29680.07530.7815.4%30.9%38.1%14.9%0.0%0.0%0.0%
Perennial Rye Grass0.29840.20740.7815.5%31.0%46.6%27.6%4.4%0.0%0.0%
Sudan Grass0.17240.06020.789.0%17.9%26.9%35.9%38.8%11.9%0.0%
Tall Fescue0.20960.10110.7810.9%21.8%32.7%43.6%28.4%0.0%0.0%
Tall Wheat Grass0.27480.25800.789.9%28.6%42.9%40.1%18.6%0.0%0.0%
Trefoil, big0.76800.22010.7839.9%2.2%0.0%0.0%0.0%0.0%0.0%
Trefoil, Birdsfoot narrow leaf0.40000.25000.7820.8%41.6%31.4%0.2%0.0%0.0%0.0%
Vetch0.44400.16790.7823.1%46.2%12.9%0.0%0.0%0.0%0.0%
Wheat Grass0.16640.15480.786.2%17.3%26.0%34.6%43.3%24.6%0.0%
Wild Rye, beardless0.23840.08030.7812.4%24.8%37.2%33.6%15.1%0.0%0.0%
Figure A1. Feed water and product water salinity measurements. (a), Feed water salinity versus Product water salinity; (b), Feed water salinity versus Feed Water electrical conductivity; (c), Product water salinity versus Product water electrical conductivity; (d), Feed water electrical conductivity versus Product water electrical conductivity.
Figure A1. Feed water and product water salinity measurements. (a), Feed water salinity versus Product water salinity; (b), Feed water salinity versus Feed Water electrical conductivity; (c), Product water salinity versus Product water electrical conductivity; (d), Feed water electrical conductivity versus Product water electrical conductivity.
Water 14 02893 g0a1
Figure A2. Ion removal. (a), Feed water chloride ion concentration versus product water chloride ion concentration; (b), Feed water sodium ion concentration versus product water sodium ion concentration; (c), Feed water salinity versus proportion of chloride ions in the feed water removed; (d), Feed water salinity versus proportion of sodium ions in the feed water removed; (e), Proportion of chloride ions in the feed water removed versus proportion of sodium ions in the feed water removed; (f), Feed water electrical conductivity versus the proportion of ions removed from the product water; Blue dots = chloride ions; red dots = sodium ions.
Figure A2. Ion removal. (a), Feed water chloride ion concentration versus product water chloride ion concentration; (b), Feed water sodium ion concentration versus product water sodium ion concentration; (c), Feed water salinity versus proportion of chloride ions in the feed water removed; (d), Feed water salinity versus proportion of sodium ions in the feed water removed; (e), Proportion of chloride ions in the feed water removed versus proportion of sodium ions in the feed water removed; (f), Feed water electrical conductivity versus the proportion of ions removed from the product water; Blue dots = chloride ions; red dots = sodium ions.
Water 14 02893 g0a2aWater 14 02893 g0a2b
Figure A3. Desalination as a function of sequential batch number. (a), Sequential batch number versus desalination; (b), Sequential batch number versus chloride removed; (c), Sequential batch number versus sodium removed; (d), Sequential batch number versus electrical conductivity removed.
Figure A3. Desalination as a function of sequential batch number. (a), Sequential batch number versus desalination; (b), Sequential batch number versus chloride removed; (c), Sequential batch number versus sodium removed; (d), Sequential batch number versus electrical conductivity removed.
Water 14 02893 g0a3
Figure A4. Redox parameters. (a), Feed water pH versus product water pH; (b), Feed water Eh versus product water Eh; (c), Feed water pH versus feed water Eh; (d), Product water pH versus product water Eh.
Figure A4. Redox parameters. (a), Feed water pH versus product water pH; (b), Feed water Eh versus product water Eh; (c), Feed water pH versus feed water Eh; (d), Product water pH versus product water Eh.
Water 14 02893 g0a4
Figure A5. Change in redox parameters and their relationship to desalination. (a), increase in Product water pH versus Increase in product water Eh; (b), Increase in product water pH versus Desalination; (c), Increase in product water Eh versus desalination; (d), Temperature versus desalination.
Figure A5. Change in redox parameters and their relationship to desalination. (a), increase in Product water pH versus Increase in product water Eh; (b), Increase in product water pH versus Desalination; (c), Increase in product water Eh versus desalination; (d), Temperature versus desalination.
Water 14 02893 g0a5
Figure A6. Relationship between ion removal and time spent in the reaction environment: (a) 3 h; (b) 15 h; (c) 30 h; the 3 h to 30 h results are determined from the data in Appendix A Table A3.
Figure A6. Relationship between ion removal and time spent in the reaction environment: (a) 3 h; (b) 15 h; (c) 30 h; the 3 h to 30 h results are determined from the data in Appendix A Table A3.
Water 14 02893 g0a6
Figure A7. Probability distribution for ion removal versus hours on line. (a), Desalinations vs. hours on line; (b), Sodium removal versus hours on line; (c), Chloride removal versus hours on line. Black line = minimum ion removal; green line is 25% probability of a lower removal; red line is 50% probability of a lower ion removal; brown line is a 75% probability of a lower ion removal; blue line is a 98% probability of a lower ion removal. Vertical orange line is this study. Vertical purple line is a separate reactor type [10], which terminated operation after 3 h.
Figure A7. Probability distribution for ion removal versus hours on line. (a), Desalinations vs. hours on line; (b), Sodium removal versus hours on line; (c), Chloride removal versus hours on line. Black line = minimum ion removal; green line is 25% probability of a lower removal; red line is 50% probability of a lower ion removal; brown line is a 75% probability of a lower ion removal; blue line is a 98% probability of a lower ion removal. Vertical orange line is this study. Vertical purple line is a separate reactor type [10], which terminated operation after 3 h.
Water 14 02893 g0a7

References

  1. Wang, X. Managing Land Carrying Capacity: Key to Achieving Sustainable Production Systems for Food Security. Land 2022, 11, 484. [Google Scholar] [CrossRef]
  2. Rosa, L. Adapting agriculture to climate change via sustainable irrigation: Biophysical potentials and feedbacks. Environ. Res. Lett. 2022, 17, 063008. [Google Scholar] [CrossRef]
  3. Negacz, K.; Malek, Z.; Vos, A.; Vellinga, P. Saline soils worldwide: Identifying the most promising areas for saline agriculture. J. Arid Environ. 2022, 203, 104775. [Google Scholar] [CrossRef]
  4. Machado, R.M.A.; Serralheiro, R.P. Soil Salinity: Effect on Vegetable Crop Growth. Management Practices to Prevent and Mitigate Soil Salinization. Horticulturae 2017, 3, 30. [Google Scholar] [CrossRef]
  5. Ayers, R.S.; Westcot, D.W. Water Quality for Agriculture; FAO Irrigation and Drainage Paper; Food and Agriculture Organization of the United Nations: Rome, Italy, 1985; Volume 29. [Google Scholar]
  6. Wei, C.; Li, F.; Yang, P.; Ren, S.; Wang, S.; Wang, Y.; Xu, Z.; Xu, Y.; Wei, R.; Zhang, Y. Effects of Irrigation Water Salinity on Soil Properties, N2O Emission and Yield of Spring Maize under Mulched Drip Irrigation. Water 2019, 11, 1548. [Google Scholar] [CrossRef]
  7. Amer, R. Spatial Relationship between Irrigation Water Salinity, Waterlogging, and Cropland Degradation in the Arid and Semi-Arid Environments. Remote Sens. 2021, 13, 1047. [Google Scholar] [CrossRef]
  8. Slater, Y.; Reznik, A.; Finkelshtain, I.; Kan, I. Blending Irrigation Water Sources with Different Salinities and the Economic Damage of Salinity: The Case of Israel. Water 2022, 14, 917. [Google Scholar] [CrossRef]
  9. Silber, A.; Israeli, Y.; Elingold, I.; Levi, M.; Levkovitch, I.; Russo, D.; Assouline, S. Irrigation with desalinated water: A step toward increasing water saving and crop yields. Water Resour. Res. 2015, 51, 450–464. [Google Scholar] [CrossRef]
  10. Antia, D.D.J. Provision of Desalinated Irrigation Water by the Desalination of Groundwater Abstracted from a Saline Aquifer. Hydrology 2022, 9, 128. [Google Scholar] [CrossRef]
  11. Antia, D.D.J. Water Treatment and Desalination Using the Eco-materials n-Fe0 (ZVI), n-Fe3O4, n-FexOyHz[mH2O], and n-Fex[Cation]nOyHz[Anion]m [rH2O]. In Handbook of Nanomaterials and Nanocomposites for Energy and Environmental Applications, 1st ed.; Kharissova, O.V., Torres-Martínez, L.M., Kharisov, B.I., Eds.; Springer: Berlin/Heidelberg, Germany, 2021; pp. 3159–3242. [Google Scholar]
  12. Terreni, J.; Trottmann, M.; Franken, T.; Heel, A.; Borgschulte, A. Sorption-Enhanced Methanol Synthesis. Energy Techn. 2019, 7, 1801093. [Google Scholar] [CrossRef]
  13. Alpay, E.; Chatsiriwech, D.; Kershenbaum, L.S.; Hull, C.P.; Kirkby, N.F. Combined reaction and separation in pressure swing processes. Chem. Eng. Sci. 1994, 49, 5845–5864. [Google Scholar] [CrossRef]
  14. Alpay, E.; Kenney, C.N.; Scott, D.M. Simulation of rapid pressure swing adsorption and reaction processes. Chem. Eng. Sci. 1993, 48, 3173–3186. [Google Scholar] [CrossRef]
  15. Antia, D.D.J. Desalination of Water Using ZVI (Fe0). Water 2015, 7, 3671–3831. [Google Scholar] [CrossRef]
  16. Antia, D.D.J. ZVI (Fe0) desalination: Catalytic partial desalination of saline aquifers. Appl. Water Sci. 2018, 8, 71. [Google Scholar] [CrossRef]
  17. Lai, B.; Chen, Z.; Zhou, Y.; Yang, P.; Wang, J.; Chen, Z. Removal of high concentration p-nitrophenol in aqueous solution by zero valent iron with ultrasonic irradiation (US–ZVI). J. Hazard. Mater. 2013, 250, 220–228. [Google Scholar] [CrossRef] [PubMed]
  18. Schwarz, H.A. Some Applications of Time-Dependent Rate Constant Theory to Radiation Chemistry. J. Chem. Phys. 1971, 55, 3650–3657. [Google Scholar] [CrossRef]
  19. Nosrati, K.; Movahedirad, S.; Amin, M.; SobatiAli, S.; Sarbanha, A. Experimental study on the pressure wave attenuation across gas-solid fluidized bed by single bubble injection. Powder Technol. 2017, 305, 620–624. [Google Scholar] [CrossRef]
  20. Huang, F.; Takahashi, M.; Guo, L. Pressure wave propagation in air-water bubbly and slug flow. Prog. Nucl. Energy 2005, 47, 648–655. [Google Scholar] [CrossRef]
  21. Lide, D.R. CRC Handbook of Chemistry and Physics, 89th ed.; CRC Press: Boca Raton, FL, USA, 2008; ISBN 13:978-1-4200-6679-1. [Google Scholar]
  22. Dimeski, G.; Badrick, T.; St John, A. Ion selective electrodes (ISEs) and interferences—A review. Clin. Chim. Acta 2010, 411, 309–317. [Google Scholar] [CrossRef]
  23. Deyhimi, F. A method for the determination of potentiometric selectivity coefficients of ion selective electrodes in the presence of several interfering ions. Talanta 1999, 50, 1129–1134. [Google Scholar] [CrossRef]
  24. Liu, T.; Wang, Z.L.; Sun, Y. Manipulating the morphology of nanoscale zero-valent iron on pumice for removal of heavy metals from wastewater. Chem. Eng. J. 2015, 263, 55–61. [Google Scholar] [CrossRef]
  25. Shu, Y.; Ji, B.; Cui, B.; Shi, Y.; Wang, J.; Hu, M.; Luo, S.; Guo, D. Almond Shell-Derived, Biochar-Supported, Nano-Zero-Valent Iron Composite for Aqueous Hexavalent Chromium Removal: Performance and Mechanisms. Nanomaterials 2020, 10, 198. [Google Scholar] [CrossRef]
  26. Hao, D.; Chen, Y.; Zhang, Y.; You, N. Nanocomposites of zero-valent iron@ biochar derived from agricultural wastes for adsorptive removal of tetracyclines. Chemosphere 2021, 284, 131342. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, Q.; Wang, Y.; Wang, Z.; Zhang, Z.; Wang, X.; Yang, Z. Active biochar support nano zero-valent iron for efficient removal of U (VI) from sewage water. J. Alloys Compd. 2021, 852, 156993. [Google Scholar] [CrossRef]
  28. Lin, J.; Sun, M.; Liu, X.; Chen, Z. Functional kaolin supported nanoscale zero-valent iron as a Fenton-like catalyst for the degradation of Direct Black G. Chemosphere 2017, 184, 664–672. [Google Scholar] [CrossRef] [PubMed]
  29. Cheng, I.F.; Muftikian, R.; Fernando, Q.; Korte, N. Reduction of nitrate to ammonia by zero-valent iron. Chemosphere 1997, 35, 2689–2695. [Google Scholar] [CrossRef]
  30. Park, M.H.; Jeong, S.; Lee, G.; Park, H.; Kim, J.Y. Removal of aqueous-phase Pb (II), Cd (II), As (III), and As (V) by nanoscale zero-valent iron supported on exhausted coffee grounds. Waste Manag. 2019, 92, 49–58. [Google Scholar] [CrossRef]
  31. Schober, P.; Boer, C.; Schwarte, L.A. Correlation Coefficients: Appropriate Use and Interpretation. Anesth. Analg. 2018, 126, 1763–1768. [Google Scholar] [CrossRef]
  32. Taylor, R. Interpretation of the Correlation Coefficient: A Basic Review. J. Diagn. Med. Sonogr. 1990, 6, 35–39. [Google Scholar] [CrossRef]
  33. Antia, D.D.J. Purification of Saline Water Using Desalination Pellets. Water 2022, 14, 2639. [Google Scholar] [CrossRef]
  34. Antia, D.D.J. Remediation of Saline Wastewater Producing a Fuel Gas Containing Alkanes and Hydrogen Using Zero Valent Iron (Fe0). Water 2022, 14, 1926. [Google Scholar] [CrossRef]
  35. Pourbaix, M. Atlas of Electrochemical Equilibria in Aqueous Solutions; NACE International: Houston, TX, USA, 1974; p. 644. [Google Scholar]
Figure 1. Process flow diagram for the reactor train used in this study.
Figure 1. Process flow diagram for the reactor train used in this study.
Water 14 02893 g001
Figure 2. Gas management, oscillating flow structure adopted in the reactor. The bubble column used in this study acts as an EPB. The EPB discharges product water to a recycle tank and receives feed water from the recycle tank (Figure 1).
Figure 2. Gas management, oscillating flow structure adopted in the reactor. The bubble column used in this study acts as an EPB. The EPB discharges product water to a recycle tank and receives feed water from the recycle tank (Figure 1).
Water 14 02893 g002
Figure 3. Example wave structuring created by an EPB, where the USA receives a constant gas flow rate through the upstream conduit and the gas is discharged in pulses from the USA to the DSA via the EPB. This results in a pulsed gas outflow through the downstream conduit from the DSA.
Figure 3. Example wave structuring created by an EPB, where the USA receives a constant gas flow rate through the upstream conduit and the gas is discharged in pulses from the USA to the DSA via the EPB. This results in a pulsed gas outflow through the downstream conduit from the DSA.
Water 14 02893 g003
Figure 4. Catalyst, transmitted light, field of view = 0.5 mm. Particles are quartz grains coated in n−FeOxHy polymers.
Figure 4. Catalyst, transmitted light, field of view = 0.5 mm. Particles are quartz grains coated in n−FeOxHy polymers.
Water 14 02893 g004
Figure 5. Expected desalination distribution associated with a 650 train reactor, processing 550 m3 d−1 (200,750 m3 a−1). 1000 iteration Monte Carlo Analysis; assumptions: 20 h batch operating time; Fe charge = 6.5 kg; Fe replacement period: 50 days; air required: 13 m3 m−1; water reactor emptying/filling period = <4 h; air compression requirements: 12 × 1.140 m3 m−1 air compressors (indicative requirements/screw compressor: power 7.5 kW (400 V, 50 Hz, 3 phase); delivery pressure = 0.8 MPa; weight = 170 kg; size (m) = 0.85 × 0.72 × 0.98; these requirements are based on an Airpress compressor APS Basic; other makes and models can be used); estimated plant compression power requirement: 90 kW h−1; 3.35 kW m−3; estimated water storage: 600–2000 m3 feed water; 600–2000 m3 product water; 600 m3 recirculating water.
Figure 5. Expected desalination distribution associated with a 650 train reactor, processing 550 m3 d−1 (200,750 m3 a−1). 1000 iteration Monte Carlo Analysis; assumptions: 20 h batch operating time; Fe charge = 6.5 kg; Fe replacement period: 50 days; air required: 13 m3 m−1; water reactor emptying/filling period = <4 h; air compression requirements: 12 × 1.140 m3 m−1 air compressors (indicative requirements/screw compressor: power 7.5 kW (400 V, 50 Hz, 3 phase); delivery pressure = 0.8 MPa; weight = 170 kg; size (m) = 0.85 × 0.72 × 0.98; these requirements are based on an Airpress compressor APS Basic; other makes and models can be used); estimated plant compression power requirement: 90 kW h−1; 3.35 kW m−3; estimated water storage: 600–2000 m3 feed water; 600–2000 m3 product water; 600 m3 recirculating water.
Water 14 02893 g005
Figure 6. Expected plant layout associated with a 650 train reactor, processing 550 m3 d−1 (200,750 m3 a−1).
Figure 6. Expected plant layout associated with a 650 train reactor, processing 550 m3 d−1 (200,750 m3 a−1).
Water 14 02893 g006
Figure 7. Analysis of stored product water. (a) Feed water salinity vs. product water salinity. Blue circles = product water salinity on removal from the reactor; larger red circles = product water salinity following two-year storage; (b) additional desalination resulting from 2-year storage.
Figure 7. Analysis of stored product water. (a) Feed water salinity vs. product water salinity. Blue circles = product water salinity on removal from the reactor; larger red circles = product water salinity following two-year storage; (b) additional desalination resulting from 2-year storage.
Water 14 02893 g007
Figure 8. Outcome distribution associated with a single reactor train. Blue line = maximum desalination; brown line = 25% probability of a higher desalination; red line = 50% probability of a higher desalination; green line = 75% probability of a higher desalination; black line = minimum desalination. Source: Appendix A, Figure A7.
Figure 8. Outcome distribution associated with a single reactor train. Blue line = maximum desalination; brown line = 25% probability of a higher desalination; red line = 50% probability of a higher desalination; green line = 75% probability of a higher desalination; black line = minimum desalination. Source: Appendix A, Figure A7.
Water 14 02893 g008
Figure 9. Expected plant layout associated with a 650 train reactor, processing 550 m3 d−1 (200,750 m3 a−1) and designed to handle the reject brine associated with a reverse osmosis desalination plant.
Figure 9. Expected plant layout associated with a 650 train reactor, processing 550 m3 d−1 (200,750 m3 a−1) and designed to handle the reject brine associated with a reverse osmosis desalination plant.
Water 14 02893 g009
Figure 10. Expected water salinities associated with the plant layout in Figure 9. Blue line = fresh feed water salinity; red line = product water from the ZVI desalination plant; green line = average water salinity entering the reverse osmosis plant; purple line = average reject brine salinity exiting from the reverse osmosis plant.
Figure 10. Expected water salinities associated with the plant layout in Figure 9. Blue line = fresh feed water salinity; red line = product water from the ZVI desalination plant; green line = average water salinity entering the reverse osmosis plant; purple line = average reject brine salinity exiting from the reverse osmosis plant.
Water 14 02893 g010
Figure 11. Schematic drawing of an aquifer in the process of being desalinated. Orange = aquifer; blue circle = part of aquifer which is being desalinated (desalination zone). Red = infiltration boreholes/wells; green = abstraction borehole/well.
Figure 11. Schematic drawing of an aquifer in the process of being desalinated. Orange = aquifer; blue circle = part of aquifer which is being desalinated (desalination zone). Red = infiltration boreholes/wells; green = abstraction borehole/well.
Water 14 02893 g011
Figure 12. Schematic process flow diagram for aquifer desalination. The assumptions provided by the example numerical analysis are: (i) volume of water processed each day is 1% of the water volume enclosed by the infiltration wells (Figure 11); (ii) the aquifers normal salinity is 5 g L−1.
Figure 12. Schematic process flow diagram for aquifer desalination. The assumptions provided by the example numerical analysis are: (i) volume of water processed each day is 1% of the water volume enclosed by the infiltration wells (Figure 11); (ii) the aquifers normal salinity is 5 g L−1.
Water 14 02893 g012
Figure 13. Change in aquifer salinity within the desalination zone, as a function of time, and the proportion of water within the desalination zone, which is fresh saline feed water, derived from the wider aquifer. The initial aquifer salinity = 5 g L−1. The volume of water processed each day is 1% of the water volume enclosed by the infiltration wells (Figure 11).
Figure 13. Change in aquifer salinity within the desalination zone, as a function of time, and the proportion of water within the desalination zone, which is fresh saline feed water, derived from the wider aquifer. The initial aquifer salinity = 5 g L−1. The volume of water processed each day is 1% of the water volume enclosed by the infiltration wells (Figure 11).
Water 14 02893 g013
Figure 14. Example of an aquifer where a reversal of seawater incursion is attempted. Green circles = abstraction wells; red circles = infiltration/injection wells; green arrows = abstracted water flows in the aquifer; red arrows = infiltrated/injected water flows of partially desalinated water; blue arrow = direction of seawater incursion.
Figure 14. Example of an aquifer where a reversal of seawater incursion is attempted. Green circles = abstraction wells; red circles = infiltration/injection wells; green arrows = abstracted water flows in the aquifer; red arrows = infiltrated/injected water flows of partially desalinated water; blue arrow = direction of seawater incursion.
Water 14 02893 g014
Figure 15. Expected structure of a reactor unit processing, continuously, 660 m3 d−1.
Figure 15. Expected structure of a reactor unit processing, continuously, 660 m3 d−1.
Water 14 02893 g015
Figure 16. Example relationships. (a), kobs vs. hours on line; (b) proportion of ion removed vs. hours on line; (c) zero−order reaction check. Cl ion concentration, g L−1 vs. hours on line; (d) first-order reaction check. Ln(Cl) ion concentration, g L−1 vs. hours on line; (e) second-order reaction check. 1/Cl ion concentration, g L−1 vs. hours on line; data source: Table A3. Red line: Trial 18. 3.61 g Cl L−1; blue line: Trial 4. 1.0 g Cl L−i.
Figure 16. Example relationships. (a), kobs vs. hours on line; (b) proportion of ion removed vs. hours on line; (c) zero−order reaction check. Cl ion concentration, g L−1 vs. hours on line; (d) first-order reaction check. Ln(Cl) ion concentration, g L−1 vs. hours on line; (e) second-order reaction check. 1/Cl ion concentration, g L−1 vs. hours on line; data source: Table A3. Red line: Trial 18. 3.61 g Cl L−1; blue line: Trial 4. 1.0 g Cl L−i.
Water 14 02893 g016aWater 14 02893 g016b
Figure 17. Variation in Cl rate constant. (a) Initial rate constant versus feed water Cl content; (b) initial rate constant for a batch versus initial rate constant for the next batch; (c) initial rate constant for a batch versus the final rate constant for a batch. Dotted blue line is a regression line. Solid blue line is a line of equal value.
Figure 17. Variation in Cl rate constant. (a) Initial rate constant versus feed water Cl content; (b) initial rate constant for a batch versus initial rate constant for the next batch; (c) initial rate constant for a batch versus the final rate constant for a batch. Dotted blue line is a regression line. Solid blue line is a line of equal value.
Water 14 02893 g017
Figure 18. Extracted catalyst following the processing of 50 batches of water. (a) Catalyst: transmitted light, field of view = 1.3 mm; (b) catalyst: circular polarized light, field of view = 1.3 mm; (c) pore water precipitate: transmitted light, field of view = 0.03 mm.
Figure 18. Extracted catalyst following the processing of 50 batches of water. (a) Catalyst: transmitted light, field of view = 1.3 mm; (b) catalyst: circular polarized light, field of view = 1.3 mm; (c) pore water precipitate: transmitted light, field of view = 0.03 mm.
Water 14 02893 g018aWater 14 02893 g018b
Figure 19. Extracted (from the product water) amorphous, fluid filled, hydrated ellipsoid of metal polymer: transmitted light, field of view = 0.13 mm.
Figure 19. Extracted (from the product water) amorphous, fluid filled, hydrated ellipsoid of metal polymer: transmitted light, field of view = 0.13 mm.
Water 14 02893 g019
Figure 20. Possible reactor process flow where the reactor is used to remediate water.
Figure 20. Possible reactor process flow where the reactor is used to remediate water.
Water 14 02893 g020
Table 1. Desalination statistics: data: Appendix A, Table A2. Values are the proportion removed.
Table 1. Desalination statistics: data: Appendix A, Table A2. Values are the proportion removed.
NaClNa+ IonsCl Ions
Mean41.84%41.33%41.74%
Standard Deviation15.61%16.35%17.99%
Minimum5.54%0.00%4.62%
1st Quartile30.01%31.38%26.16%
Median43.39%41.14%42.21%
3rd Quartile51.09%53.67%55.25%
Maximum73.59%69.77%75.51%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Antia, D.D.J. Catalytic Partial Desalination of Saline Water. Water 2022, 14, 2893. https://doi.org/10.3390/w14182893

AMA Style

Antia DDJ. Catalytic Partial Desalination of Saline Water. Water. 2022; 14(18):2893. https://doi.org/10.3390/w14182893

Chicago/Turabian Style

Antia, David D. J. 2022. "Catalytic Partial Desalination of Saline Water" Water 14, no. 18: 2893. https://doi.org/10.3390/w14182893

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop