CHEMICAL ENGINEERING TRANSACTIONS VOL. 52, 2016 A publication of The Italian Association of Chemical Engineering Online at www.aidic.it/cet Guest Editors: Petar Sabev Varbanov, Peng-Yen Liew, Jun-Yow Yong, Jiří Jaromír Klemeš, Hon Loong Lam Copyright © 2016, AIDIC Servizi S.r.l., ISBN 978-88-95608-42-6; ISSN 2283-9216 Use of Microwaves for Disinfection of Farmland: a Feasibility Study Francesco Desogus*a, Alessandro Fantib, Sergio Casub, Michele Spanub, Matteo Bruno Lodib, Giuseppe Mazzarellab, Bruno Biscegliac aUniversity of Cagliari – Department of Mechanical, Chemical and Material Engineering – Piazza d’Armi, 09123 Cagliari (I) bUniversity of Cagliari – Department of Electric and Electronic Engineering – Piazza d’Armi, 09123 Cagliari (I) cUniversity of Salerno – Department of Industrial Engineering – Via Giovanni Paolo II 132, 84084 Fisciano (I) f.desogus@dimcm.unica.it This work focuses on the heat transfer dynamics in agricultural soils when exposed to microwave (MW) fields, in order to disinfect them and to eliminate dangerous organisms without using pesticides. The work managed with the design of a horn antenna, with fixed feeding power, to be used to irradiate the superficial layer of soils for a fixed depth. The soil was approximated with a transmission line model, i.e. a lumped electrical model. By this model, after evaluating the relevant physical parameters such as the dielectric constant (as a function of the soil water content), it was possible to evaluate the power absorption and, finally, the increase in temperature and its profile as a function of the irradiation time. The latter information, strictly connected to the heat transfer process, demonstrated the feasibility of the methodology here proposed to achieve difficult conditions for microbial and pathogens life and so an effective disinfection of the treated farmland. 1. Introduction Recently the scientific community's interest in biological agriculture has widely increased. "Bio-Agriculture" means precise cultivation methods, excluding the use of chemical products such as herbicides and pesticides. In recent times the use of electromagnetic (EM) waves, especially microwaves (MW), for soil disinfection from pathogen agents has been taken into account. Microwaves are important because lead to an increase in the rotational energy of molecules, especially those of water, which possess an electric dipolar moment with angular frequency equal to that of MW (Fanti et al., 2016), and for this reason they are widely diffused in chemical (Desogus et al., 2016) and biochemical processes (Fanti et al., 2015), particularly with aqueous solutions (Casu et al., 2016). As bacteria and microorganisms are mainly composed of water, the energy transfer, i.e. the temperature increase, can be exploited to inhibit or eliminate them, creating intolerable environmental conditions (Nelson, 1996). On the other hand, microwave effects on bacteria have been showed in the literature (Carta and Desogus, 2010), and also non-thermal effects (Carta et al., 2006). Microwaves have often been proposed as an enhancing factor in environmental remediation (Carta and Desogus, 2013), and also in the soil treatment for the elimination of chemical pollutants (Acierno et al., 2003), but, at our knowledge, it cannot be found in the literature any testing and analysis of soil disinfection with in situ MW applications (Baghdadi et al., 2008). The main studies on the soil response to MW irradiation were conducted using radar remote sensing, taking into account the surface roughness (Ulaby et al., 1978), also developing semiempirical backscattering models (Ulaby et al., 1982) and using inversion techniques (Ulaby et al., 1986). Also other authors focused on the passive response of soil to MW (Schmugge et al., 1986), and a reliable and accurate review of them was done by Das and Paul (2015). Moreover, these analyses exploited free-space techniques, as done by Hallikainen et al. (1985) and by Valbonesi and Bisceglia (2015), or driving structures, e.g. waveguides or horn antennas (Helszajn, 2000). Moreover, the studies on soil irradiation ranged from the determination of the electric parameters, i.e. the quantification of electric permittivity of several soil layers taking into account the content of water (Wang and Schmugge, 1978) or moisture (Fung, DOI: 10.3303/CET1652200 Please cite this article as: Desogus F., Fanti A., Casu S., Spanu M., Bruno Lodi M., Mazzarella G., Bisceglia B., 2016, Use of microwaves for disinfection of farmland: a feasibility study, Chemical Engineering Transactions, 52, 1195-1200 DOI:10.3303/CET1652200 1195 1994), to the frequency response of the soil (Hallikainen et al., 1985), but also the prediction and prevention of possible hydrological phenomena in areas with risk through satellite and georadar detection (Guattari et al., 2015). Regarding agricultural soils, nowadays the problem of the elimination of pathological agents (algae, fungi and fungine roots, viruses, nematodes, protozoa) is perceived as fundamental, especially for the quality of the cultivated products. Since many biological species can be pathogenic, the agricultural soil should be disinfected without affecting the organic nutrients. Aim of this work is to assess the feasibility of the use of high-power microwaves to disinfect the soil. The MW use seems to be promising since it does not affect organic substances if the temperature is suitably controlled to prevent overheating. The chosen exposure system consists of a pyramidal horn antenna feed by a waveguide, similar to a literature device (Helszajn, 2000), often used for this kind of studies, as mainly done by Hallikainen et al. (1985) and much more recently Bebawi et al. (2007), but herein designed according to the approach described by Simone et al. (2015) using the PSO algorithm for the rectangular feeding waveguide (Simone et al., 2016). 2. Physics soil properties Soil is defined as a mixture of fragmented and eroded minerals, of decomposing organic matter and a certain amount of water and air. Four are the main soil components, which are listed in Table 1. These features, i.e. its composition, establish the texture of soil, which is an important characteristic identifying how the ground under analysis retains water or moisture, the amount of which is a relevant parameter influencing the dielectric constant. Soil is a porous medium, and typically, for farmland, porosity is around 70 %. Such pores can be filled with water or air. If they contain mainly the latter, the soil is dry or sandy; on the other hand, if they are mainly filled in by the former, the soil can be classified as clayey, i.e. saturated soil. Water can be found in pores in different states: vapour, free-and bound-water, depending on the soil composition, which influences the affinity between water molecules and solid surfaces. Hence it is possible to define the field capacity (FC) as a function of the water content (humidity) in the soil, which varies depending on the ground porosity, as stated by the expression of Wang and Schmugge (1978): FC = 25.1 0.21 S +0.22 C   (1) where S and C are the percentual amount of sand and clay, respectively. From FC the volumetric amount of water (mv) can be obtained, dividing the water fractional volume by the soil bulk density (Behari, 2005). Another relevant material for the soil characteristics is peat, which is an organic formation mainly composed by vegetal residues impregnated by water, sunk in an acid environment, so that they cannot completely decompose. However, peat may include many different kinds of organic materials, e.g. insects or animals remains. Moreover, such kind of material is highly porous, with variable texture (similar to clay or gravel), depending on its porosity (46 - 84 %) (Ayalew et al., 2007), and extremely wet, these listed conditions allow the proliferation of several microorganism populations, namely pathogens and parasites. 3. Electromagnetic model 3.1 Antenna design To irradiate the soil, we designed a horn antenna with pyramidal base, composed by a waveguide and fixed aperture dimensions. The radiation system was fed with an available power PMW at a frequency of 2.45 GHz and it exhibited 10 dB gain. The behavior and response of our antenna was analyzed using a Finite Element Method software for high frequencies EM problems. The horn was placed at the distance of 0.95 m from the soil surface, and produced a Poynting vector modulus Si = PMW/A1 over the irradiated area A1 of 0.72 m2. 3.2 Dielectric constant of soil In order to establish the soil behavior under MW exposition, it was necessary to determine the value of its relative dielectric constant () for humidity ranging from 0 ÷ 10 (dry soil) to 30 ÷ 40 % (wet soil), and also 50 ÷ 90 % (peat). To quantify such physical parameter (both the real and the imaginary parts, ’ and ”), we exploited the empirical model of the following Eq(2) developed by Hallikainen et al. (1985):         20 1 2 0 1 2 v 0 1 2 v' " = a + a S+ a C + b +b S+b C m + c +c S+c C m         or (2) in which aj, bj and cj are fitting parameters and j refers to the component by the percentual amount of which the parameter must be multiplied (0 for none, 1 for sand, 2 for clay). Of course, the soil behavior under radiation could also be measured using frequency domain (Mazzarella and Panariello, 1991) or time domain (Costanzo et al., 2013) reflectometry. 1196 Table 1: Classification of soil granulometric classes (diameter size in mm) (Soil Conservation Service, 1987) Clay Silt Sand Gravel Very fine Fine Medium Coarse Very coarse < 0.002 0.002÷0.05 0.05÷0.1 0.1÷0.25 0.25÷0.5 0.5÷1.0 1.0÷2.0 2.0÷50.0 3.3 Equivalent circuit model of soil through analytical transmission line model To evaluate electrical parameters when soil is irradiated, we modeled it with a lumped parameter model, approximating it to a transmission lines system, as shown in Figure 1. Since soil is composed of several layers of different texture and moisture (Hillel, 1998), and considering the antenna as a power source, every layer i is considered as a line with impedance Zi, length Li and propagation constant ki. These values depend on the dielectric constant and can be calculated using the following expressions (Guattari et al., 2015): i 0 ik =   (3) i 0 iZ = Z  (4) where β0 is the propagation constant of vacuum. The input impedance of the soil is different from the free- space one, so that a suitable matching network is required. Since the required bandwidth is very small, we used, in our model, the simple network of Figure 1b. 4. Soil biological environment Soil is a living medium, in fact, thanks to water and air circulating in pores, and also to nutrients, it represents a favorable environment for life. Bacteria, fungi, plants and animals, e.g. invertebrates and small vertebrates, find suitable conditions to survive and reproduce, creating a fragile interacting dynamic system (Mitscherlich and Marth, 1984). Each of them has a specific role, function and position in the soil matrix, but it is necessary to distinguish between pathological and risky agents (nematodes, fungi and parasites), and favorable organisms, e.g. earthworms. These biological entities are mainly composed by water, ranging from 80 to 90 %. The study by Komarova et al. (2008) underlined that MW exert a lethal influence on bacteria, due to the increase in temperature and to the bioeffects of non-ionizing radiations, even if MW are applied to soils (Vela and Wu, 1979), but the explanation of the latter remains, in part, still unclear (Vela et al., 1976). Actually, some harmful agents are able to withstand high temperatures, even for long times of exposition. Experimental findings show that the treatment efficiency is strictly related to the exposure time, therefore, to reach the lethal temperature, such information should be related to power levels, so that future in situ actions can be designed and undertaken (Valbonesi and Bisceglia, 2015). We assumed, in this study, that the temperature of 60 °C must be reached for about 30 min in a layer about 3 cm thick (this thickness was chosen according to the information from Soil Survey Division Staff, 1993). Since the temperature, for a constant MW power, would increase with time, the fixed value can be maintained for a few minutes by a suitable reduction of the horn power or turning off the generator. 5. Methods and model The moisture values used in the present work were acquired from a measurement campaign performed on Sardinian farmland using TDR (Time Domain Reflectometry) probe and gravimetric measures. Hence, data on the superficial layer of soil were directly exploited. The texture of analyzed soils is shown in Table 2. In order to evaluate the actual heating of the soil, we need to solve the following heat equation:   T c = k T + P t      (5) where ρ is the soil (bulk) density, ranging from 200 to 1,800 Kg m-3, c is its specific heat, ranging from 800 to 3,300 J Kg-1 K-1, k the thermal conductivity, ranging from 0.01 to 4 W m-1 K-1 (Coppola et al., 2007) (both c and k are relative to the soil considered as a pseudo homogeneous medium), T is the temperature (in K), t is time (in s), and P is the specific heat power (in W m-3) generated within the volume by the MW power source. Since both the electromagnetic and the thermal properties of the soil are temperature-dependent, Eq(5) is non- linear. In order to assess the feasibility of the proposed approach, we decided to solve a simplified linear and one dimensional model. With reference to Figure 1, let be z ≥ 0 the soil semi-space coordinate. Then, the generated EM power is: 1197 Figure 1: (a) first layer of the soil; (b) equivalent transmission line model with an incident plane wave on first layer and quarter-wave dielectric slugs impedance matching      , with E-2 z E iP z = P 0 e P 0 = 2 S   (6) where αE is the imaginary part of the dielectric constant (Hallikainen et al., 1985) and is one of the main indicators of the feasibility of the technique here proposed. The power is produced only in the upper layer of the soil and then heat transfers toward the deepest layers (and in atmosphere), so the other relevant parameters are k and the average convection coefficient at the soil surface, hc, considering that the heat power exchange with air (qA) is given by the following expression, in which TA is the air temperature:   A c Aq = h T 0 T  (7) The convection coefficient depends essentially on the air properties and, for the soil free surface, is around 2 W m-2 K-1. However, if necessary, this value can be reduced by covering the soil with a layer of a smooth material; moreover, if an insulating material (e.g. polystyrene foam, which is transparent to the EM waves) is chosen for the cover, the heat flux to air can be further reduced (down to 0.5 W m-2 K-1). On the contrary, the heat exchange with air can be significantly increased using forced convection, of more than one order of magnitude. 6. Results and discussion We considered a horn power of 8 kW with Si=11 kW m-2 and an initial air and soil temperature of 20 °C. For a dry soil, αE can be around 6÷8 m-1, so that the EM power is dissipated quite uniformly in the layer of interest (2÷3 cm depth from the surface). In this case the heating is near independent from the thermal parameters of the soil, but P(0) is very low (around 150 kW m-3) and so about 7 - 8 min would be necessary to reach the temperature of 60 - 65 °C in the whole layer of interest (see Figure 2). However, since the temperature field is almost uniform, an increase of the MW power would reduce the heating time without overheating, which could damage the soil nutrients. A larger value of αE, around 15 m-1 (i.e., a wetter soil) would produce a faster heating of the surface but would require a low k value (near to the lower end of its possible range). This could seem strange, since it is the heat diffusion that now increases the temperature of the lower portions of the useful soil. But the same effect is also responsible of the heat loss toward the deepest layers, where we do not need a disinfection. On the other hand, since such a value of αE easily leads to surface overheating, hc should be carefully chosen. Even larger values of αE (near to 30 m-1, typical of very wet soils) allow a quick heating, but a longer time is required to heat the lower layers, and the surface can easily reach 100 °C. From this discussion it follows that, using a proper horn power, we can reach the required temperatures for dry and slightly wet soils, preventing the surface overheating. However, for a given farmland soil, an optimization is needed to get the best water content and hc value to reduce the overall disinfection time. This optimization, at first, could be done by a linear model, then it must be refined by a coupled EM-thermal non-linear one. 7. Conclusions The results here presented show that MW heating can be a viable method to disinfect the upper layer of agricultural soils. However, a suitable optimization must be done, to select the optimal sets of parameters and among them the water content of the soil. Since this can be modulated by irrigation, a suitable protocol must be derived, depending on the initial soil characteristics, external temperature and air motion conditions. The higher air temperature, the better temperature distribution versus depth and time. In other words, for higher air temperatures, a more uniform optimal temperature distribution is reached in a shorter time. 1198 Table 2: Collected soil texture data of some Sardinian farmlands (Montaldo et al., 2008; 2013) Fields Coordinates (WGS84 GD) Soil Texture (%) Longitude Latitude Sand Silt Clay A – Silt loam 9.16691920436652 39.7029696452146 29.0 50.1 20.9 E – Clay loam 9.22181049758689 39.6876325325612 22.4 50.6 27.0 F – Loam 9.22368411780198 39.6900166394511 29.3 48.5 22.2 G – Silt 9.21743569879174 39.7317365306478 5.9 90.8 3.3 H – Silty clay loam 9.15900027986903 39.6883121736779 17.7 52.0 30.3 Figure 2: Soil temperature as a function of depth and time of irradiation (ρ=1,200 kg m-3; c=1200 J kg-1 K-1; k=0.5 W m-1 K-1; hc=2.0 W m-2 K-1; Si=11 kW m-2; MW irradiation stopped after 7.2 min; temperatures in °C) Acknowledgments This work was supported in part by Regione Autonoma della Sardegna under contract "IDROSAR" CRP- 79793 (CUP F72I15000870002). References Acierno D., Barba A.A., d’Amore M., 2003, Microwaves in soil remediation from VOCs. 1: Heat and mass transfer aspects, AIChE J. 49, 1909-1921. Ayalew G., Holden N.M., Ward S.M., 2007, Microwave dielectric properties of horticultural peat products, Mires Peat 2, article 3. Baghdadi N., Zribi M., Loumagne C., Ansart P., Anguela T.P., 2008, Analysis of TerraSAR-X data and their sensitivity to soil surface parameters over bare agricultural, Remote Sens. Envoiron. 112, 4370-4379. Bebawi E.F., Cooper A.P., Brodie G.I., Madigan B.A., Vitelli J.S., Worsley K.J., Davis K.M., 2007, Effect of microwave radiation on seed mortality of rubber vine (Cryptostegia grandiflora R.Br.), parthenium (Parthenium hysterophorous L.) and bellyache bush (Jatropha gossypiifolia L.), Plant Prot. Q. 22, 136-142. Behari J., 2005, Physical Properties of Soils: Microwave Dielectric Behaviour of Wet Soils, Springer with Anamaya Publishers, New Delhi, India. Carta R., Desogus F., 2010, The effect of low-power microwaves on the growth of bacterial populations in a plug flow reactor, A.I.Ch.E. J. 56, 1270-1278. Carta R., Desogus F., 2013, The enhancing effect of low power microwaves on phenol oxidation by the Fenton process, J. Environ. Chem. Eng. 1, 1292-1300. Carta R., Desogus F., Errico M., 2006, Effect of microwave radiation on the growth rate of Bacillus clausii at 37°C, CHISA 2006 – Proceedings of the 17th International Congress of Chemical and Process Engineering, August 27-31 2006, Prague, Czech Republic. Code 70283. Casu S., Fanti A., Djuric N., Desogus F., Mazzarella G., 2016, Microwave resonant cavity as a reactor for the enzymatic hydrolysis of sucrose, IEEE Proceedings of the 15th Mediterranean Microwave Symposium, November 30-December 2 2015, Lecce, Italy. Code 119001. Coppola A., Basile A., Menenti M., Buonanno M., Colin J., De Mascellis R., Esposito M., Lazzaro U., Magliulo V., Manna P., 2007, Spatial distribution and structure of remotely sensed surface water content estimated by a thermal inertia approach, IAHS-AISH P. 316, 1-12. Costanzo S., Spadafora F., Moreno H.O., Scarcella F., Di Massa G., 2013, Multiband software defined radar for soil discontinuities detection, J. Elect. Comput. Eng. 2013. Das K., Paul P. K., 2015, Present status of soil moisture estimation by microwave remote sensing, Cogent Geosci. 1, 1-21. 1199 Desogus F., Casu S., Muntoni G., 2016, Design of a chemical reactor under microwave irradiation in resonance conditions, ECOS 2016 – 29th International Conference on Efficiency, Cost, Optimization, Simulation, and Environmental Impact of Energy Systems, June 19-23 2016, Portorož, Slovenia. Fanti A., Casu S., Desogus F., Montisci G., Simone M., Casula G.A., Maxia P., Mazzarella G., Carta R., 2015, Evaluation of a microwave resonant cavity as a reactor for enzyme reactions, J. Electromagnet. Wave 29, 2380-2392. Fanti A., Casu S., Desogus F., Djuric N., Mazzarella G., 2016, Design and optimization of a microwave irradiated and resonant continuous biochemical reactor, Radio Sci. 51, 1199-1212. Fung A.K., 1994, Microwave Scattering and Emission Models and their Applications, Artech House Publishers, London, United Kingdom. Guattari C., Ramaccia D., Bilotti F., Toscano A., 2015, Permittivity of sub-soil materials retrieved through transmission line model and GPR data, Prog. Electromagn. Res. 151, 65–72. Hallikainen M.T., Ulaby F.T., Dobson M.C., El-Rayes M.A., Wu L., 1985, Microwave dielectric behavior of wet soil – Part 1: Empirical models and experimental observations, IEEE Trans. Geosci. Remote Sens. GE-23, 25-34. Helszajn J., 2000, Ridge Waveguides and Passive Microwave Components and Devices. The Institution of Engineering and Technology, London, UKingdom. Hillel D., 1998, Envoirnmental Soil Physics, Academic Press, Cambridge, US. Komarova A.S., Likhacheva A.A., Zvyagintsev D.G., 2008, Influence of microwave radiation on soil bacteria, Moscow Univ. Soil Sci. Bull. 63, 190-195. Mazzarella G., Panariello G., 1991, Reconstruction of monodimensional dielectric profiles, J. Opt. Soc. Am. A: 8, 1622-1629. Mitscherlich E., Marth E. H., 1984, Microbial Survival in the Environment – Bacteria and Rickettsiae Important in Human and Animal Health, Springer-Verlag Berlin Heidelberg, Berlin, Germany. Montaldo N., Albertson J.D., Mancini M., 2008, Vegetation dynamics and soil water balance in a water-limited Mediterranean ecosystem on Sardinia, Italy, Hydrol. Earth Syst. Sci. Discuss. 5, 219-255. Montaldo N., Corona R., Albertson J.D., 2013, On the separate effects of soil and land cover on Mediterranean ecohydrology: two contrasting case studies in Sardinia, Italy, Water Resour. Res. 49, 1123- 1136. Nelson S.O., 1996, A review and assessment of microwave energy for soil treatment to control pests, T. ASAE 39, 281-289. Schmugge T., O'Neill P. E., Wang J. R., 1986, Passive microwave soil moisture research, IEEE Trans. Geosci. Remote Sens. GE-24, 12-22. Simone M., Fanti A., Mazzarella G., 2015, Ridge waveguide optimization with PSO algorithm, J. Electromagnet. Wave 29, 199-209. Simone M., Fanti A., Montisci G., Casula G.A., Mazzarella G., 2016, Combined PSO-FDFD optimization of rectangular ridged waveguides, Appl. Comput. Electrom. 31, 144-151. Soil Conservation Service, 1987, Soil Mechanics Level I, Module 3: USDA Textural Soil Classification. United States Department of Agriculture. Soil Survey Division Staff, 1993, Soil Survey Manual, Chapter 3: Examination and Description of Soils. United States Department of Agriculture. Ulaby F.T., Batlivala P.P., Dobson M.C., 1978, Microwave backscatter dependence on surface roughness, soil moisture, and soil texture: Part I – bare soil, IEEE T. Geosci. Elect. 16, 286-295. Ulaby F.T., Moore R.K., Fung A.K., 1982, Microwave Remote Sensing: Active and Passive, vol. 2. Artech House Publishers, London, UK. Ulaby F.T., Moore, R.K., Fung A.K., 1986, Microwave Remote Sensing: Active and Passive, vol 3. Artech House Publishers, London, United Kingdom. Valbonesi S., Bisceglia B., 2015, EMF application to soil disinfestation and agro-food processing for a sustainable development, Proceedings of the 15th IEEE Mediterranean Microwave Symposium (MMS 2015), November 30-December 02 2015, Lecce, Italy. Vela G.R., Wu J.F., 1979, Mechanism of lethal action of 2.450 MHz radiation on microorganisms, Appl. Environ. Microbiol. 37, 550-553. Vela G.R., Wu J.F., Smith D., 1976, Effect of 2450 MHz microwave radiation on some soil microorganisms in situ, Soil Sci. 121. Von Hippel A.R., 1954, Dielectric Materials and Applications. Artech House Publishers, London, UK. Wang J.R., Schmugge T.J., 1978, An empirical model for the complex dielectric permettivity of soils as a function of water content, IEEE Trans. Geosci. Remote Sens. GE-18, 288-295. 1200