SciELO - Scientific Electronic Library Online

 
vol.87 issue213Cement production with pozzolans from residual tropical soils formed from paragneiss with a high silicon oxide contentAgronomic response of forage mixtures in a silvopastoral system in the Colombian dry tropics author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


DYNA

Print version ISSN 0012-7353On-line version ISSN 2346-2183

Dyna rev.fac.nac.minas vol.87 no.213 Medellín Apr./June 2020

https://doi.org/10.15446/dyna.v87n213.75361 

Artículos

Calibration of two capacitive soil moisture sensors in Ultisol

Calibración de dos sensores capacitivos de humedad en un Ultisol

Anabel de los Ángeles Calva Jiménez a  

Ceres Duarte Guedes Cabral de Almeida b  

José Amilton dos Santos Júnior a  

Célia Silva dos Santos a  

a Universidade Federal Rural de Pernambuco, Recife, Brasil. any_angeles1991@hotmail.com, celia@agro.eng.br, eng.amiltonjr@hotmail.com, joseamilton@ufrpe.br

b Colégio Agrícola Dom Agostinho Ikas, Universidade Federal Rural de Pernambuco, São Lourenço da Mata, Brasil. ceres.codai.ufrpe@gmail.com, ceres.cabral@ufrpe.br


Abstract

The objectives of this study were to perform the calibration and determine the calibration equations of two ECH2O soil moisture sensors, for an Ultisol, collected in the Municipality of Pesqueira-PE located in Brazilian semi-arid region. The study was carried out in 5-liter plastic containers in a protected environment at the Federal Rural University of Pernambuco, Recife-PE, where soil moisture was daily monitored for 51 days from saturation moisture to natural drying. The authors verified that the ECH2O sensors showed an exponential response to the variation of soil moisture, obtaining R2>0.97. RMSE values ranged from 0.065 to 0.012 m3 m-3 for the 5TE sensor and from 0.045 to 0.011 m3 m-3 for the EC-5 sensor. The researchers concluded that the calibration equations obtained were more accurate than the factory equations for that type of soil.

Keywords: capacitance; water content soil probe, 5TE, EC-5

Resumen

Los objetivos de este estudio fueron realizar la calibración y determinar las ecuaciones de calibración de dos sensores de humedad de suelo ECH2O para un suelo Ultisol colectado en el Municipio de Pesqueira-PE, localizado en el Semiárido Brasileño. El estudio se llevó a cabo en baldes plásticos de 5 litros implementados en ambiente protegido en la Universidad Federal Rural de Pernambuco, Recife-PE, donde se realizó el monitoreo diario de la humedad del suelo durante 51 días, desde la humedad de saturación hasta el secado natural. Los autores verificaron que los sensores ECH2O presentaron una respuesta exponencial a la variación de humedad de suelo obteniendo R2>0,97. Valores de RMSE variaron de 0,065 a 0,012 m3 m-3 para el sensor 5TE y de 0,045 a 0,011 m3 m-3 para el EC-5. Los investigadores concluyeron que las ecuaciones de calibración obtenidas fueron más precisas que las ecuaciones de fábrica para ese tipo de suelo.

Palabras clave: capacitancia; sonda de humedad del suelo, 5TE, EC-5

1. Introduction

One of the biggest problems at present is the lack of water, as this resource is becoming increasingly scarce. Likewise, the lack of water for several uses forces the search for alternatives that can contribute to its adequate management. Agriculture is a production sectors that consumes around 70% of the water in the world [1]. Most of this water is used in irrigated crops, and its improper use is causing great economic losses for farmers.

For this reason, monitoring the soil water content is a fundamental technique for various economic activities in agriculture [2]. The control of soil moisture helps farmers to take proper and responsible management of this resource since it allows decisions to be made as to when and how much to irrigate, optimizing the production of their crops, reducing water losses due to runoff and percolation, thus reducing the environmental impacts [3].

There are direct and indirect methods that allow determining the water content in the soil [4,2]. Among the direct methods is the gravimetric method, used as a standard method. One of the disadvantages of this method is that, in addition to requiring 24 hours to obtain the results, it is a destructive method, so it cannot be used to take soil moisture measurements in the same site [5,6].

At present, indirect methods are gaining a great deal of space in agriculture, most of these are based on electromagnetic techniques to estimate the water content in the soil in situ, in real time and with acceptable accuracy. Indirect methods are not destructive, which is a great advantage since soil moisture can be estimated in the same place [7-9]. Among the methods for indirectly determining soil moisture, radiometric tests (Neutron Probe), Frequency Domain Reflectometry (FDR), Time Domain Reflectometry (TDR) can be cited [10].

ECH2Os are among the sensors that measure soil moisture content based on capacitance. This group of sensors measure the apparent dielectric permittivity of the surrounding environment, generally working with a frequency of 70 MHz [11], and include the 5TE and EC-5 models. The 5TE sensor (length = 10.9 cm, width = 3.4 cm, depth = 1.0 cm), in addition to measuring the volumetric soil moisture, it measures temperature and electrical conductivity using capacitance technology [12]. EC-5 is a small sensor (length = 8.9 cm, width = 1.8 cm, depth = 0.7 cm), suitable for use in greenhouses [6], and measures the volumetric soil moisture of the means where it is used, using capacitance [13]. These types of sensors are pre-calibrated at the factory to be used in a wide range of soils, mostly by farmers using the factory calibration equation [14]. However, due to the physical, chemical and biological properties of each soil, it is necessary to perform a specific calibration [15].

The specific calibration of this type of sensor allows reliable volumetric soil moisture data to be obtained [16]. An alternative for calibration is to perform it in representative soils of each sector so that the results of the calibration equations can be extrapolated to soils whose physical and chemical characteristics are similar. This will allow farmers to achieve maximum sensor performance [17] in addition to obtaining soil moisture content readings with greater precision and accuracy.

Due to the need to obtain calibration equations for specific soils, the objective of this work was to perform the calibration of two ECH2O sensors, 5TE and EC-5 models, for a sandy-textured soil of the Brazilian semi-arid region under controlled conditions and determine a calibration curve for each sensor.

2. Materials and methods

Soil samples were collected in the municipality of Pesqueira, located in the Brazilian semi-arid region, with the following geographical coordinates: 07º 15’ 18’’ S and 35º 52’ 40’’ W, at average altitude of 550 m. The climate of this site was classified as hot semi-arid, BSh, according to Köppen’s climate classification. According to Molinier et al. (1994), this region has an average temperature of 27 °C, relative humidity of 73% and annual precipitation of 670 mm, with highest levels of rainfall in the period from March to April. The present study was carried out in a greenhouse of the Federal Rural University of Pernambuco, in Northeast Brazil.

The soil used in this study was classified as Planossolo Háplico Sálico Sódico Hipereutrófico using the Brazilian Soil Classification System [18], corresponding to an Ultisol according to the U.S. Soil Taxonomy. The soil used has a sandy texture, with the following textural composition: 789 g kg-1 of sand, 160.5 g kg-1 of silt and 50.5g kg-1 of clay. Disturbed and undisturbed soil samples were collected at 40 cm depth. The undisturbed soil samples were collected using an Uhland soil sampler in 100-cm3 steel cylinders (5 cm in diameter by 5 cm in height) and used to determine soil density. Disturbed samples were air dried, pounded to break up clods and sieved through a 4.75-mm mesh. Five 5-L plastic buckets, previously perforated at the base, received one layer of gravel (ø = 2 cm), covered by a nonwoven textile, and a usable volume of 4.66 L was left, being subsequently filled with soil.

The amount of air-dried soil used in each container in order to maintain soil density (1.43 g cm-3) was 7024 g. Buckets with soil were saturated by capillarity and placed on a bench to allow drainage of excess water and to reach field capacity.

Subsequently, two capacitive sensors, 5TE and EC-5 models from Decagon Devices, were placed in each bucket, which had been previously weighed, installed vertically 0.05 m below the surface, 0.05 m away from the bucket wall and equidistant from each other. Daily weightings and readings of soil moisture were measured at 8:00 a.m. The capacitive sensors estimated the volumetric soil moisture in m3 m-3, and the readings were taken using a portable data logger for 51 days, when the mass of the set (bucket, gravel, nonwoven textile and soil) became constant.

The volumetric soil moisture was calculated based on the gravimetric soil moisture using Eq. (1):

Where θ is the volumetric soil moisture expressed in m3 m-3, U is the gravimetric soil moisture expressed in g g-1 and Ds is the soil density expressed in g cm-3.

For the linear regression analysis, the volumetric soil moisture measured, Y (m3 m-3), and estimated by the sensors, X (m3 m-3), were related. The following statistical indices were used to compare the volumetric soil moisture data measured and estimated by the sensors: coefficient of determination (R2), root mean squared error (RMSE), calculated by the Eq. (2) defined in the reference [19], and the Willmott’s index of agreement, expressed through Eq. (3) determined by [20]:

Where n is the number of observations, Oi and Pi are the observed and estimated values, respectively (i = 1, 2,... n), O and P are the averages of the observed and estimated values, respectively.

The Willmott’s index varies between 0 and 1; a value of 1 indicates a perfect match between the measured and estimated values, while 0 indicates no agreement at all between those values [5].

To determine the performance of the proposed methods, the performance index (c) was used according to Eq. (4), proposed by [21].

Where R is the correlation coefficient and d the agreement index. For the interpretation of the coefficient of performance values (Table 1), the classification proposed by [21] was used.

Table 1 Performance classification according to the value of the “c” coefficient. 

Source: [21]

3. Results and discussion

The drying curve shown in Fig. 1 was generated from the measured volumetric soil moisture data generated based on soil density and the calibrated values of volumetric soil moisture obtained from the data estimated by the 5TE and EC-5 sensors.

Source: The Authors.

Figure 1 Soil drying curve according to volumetric soil moisture values obtained by gravimetry based on soil density and volumetric soil moisture values obtained using the calibration equations for each sensor during the observation period. 

As shown in Fig. 1, the volumetric soil moisture values corrected with the calibration equation of each sensor are close to those calculated, but the soil moisture values of the 5TE sensor are closer to the curve determined by gravimetry. In the range from 0.33 to 0.08, the EC-5 sensor lost sensitivity, estimating soil moisture readings different from those calculated by gravimetry.

The calibration equation (5) for the 5TE sensor was obtained by measuring the ratio between the average volumetric soil moisture calculated for the five buckets based on soil density and the average volumetric soil moisture after calibration, as shown in Fig. 2.

Source: The Authors.

Figure 2 Volumetric soil moisture estimated by the 5TE sensor before and after calibration and its correlation with the calculated volumetric soil moisture. 

Where ( corresponds to the actual volumetric soil moisture expressed in m3 m-3 and ( STE is the volumetric soil moisture given by the 5TE sensor according to the factory calibration.

The statistical model that fitted best was the exponential, obtaining a coefficient of determination (R2) of 0.985.

The volumetric soil moisture values obtained with the 5TE sensor with the factory calibration overestimated the volumetric soil moisture values calculated within the moisture range from 0.11 to 0.26 m3 m-3 (values located below the 1:1 line). Note that the volumetric soil moisture calculated for these values varied from 0.06 to 0.23 m3 m-3 (Fig. 2), indicating that there is a need to calibrate the sensor for this type of soil. As mentioned above, the minimum moisture value estimated by the sensor was 0.11 m3 m-3, which means that, in order to determine the water content in the soil at field capacity, the sensor would not be useful since the value for this soil according to the soil water characteristic curve was 0.096 m3 m-3; however, when calibrating the sensor, the field capacity value falls within the soil moisture range.

From 0.26 m3 m-3 of moisture, the sensor underestimated the calculated volumetric soil moisture values, where it is clearly noted that the sensor reduced the sensitivity to detect the soil moisture values within that range.

After calibration of the 5TE sensor, it can be seen that the soil moisture values are close to the 1:1 line. However, there was an overestimation in lower intensity within the soil moisture ranges from 0.11 to 0.15 and from 0.20 to 0.23 m3 m-3, whereas in the range from 0.16 to 0.19 m3 m-3 there was overestimation with greater intensity.

From 0.24 m3 m-3 of moisture, there was dispersion in the values, which shows that for soil moisture greater than this value the 5TE sensor does not work properly, despite being calibrated, showing an underestimation error of 0.0372.

Evaluating the performance of the 5TE sensor in the estimation of the volumetric soil moisture, [22] concluded that for clay-textured soils the sensor requires calibration, showing an estimation error of ( 0.025 m3 m-3, confirming the results of this study.

On the other hand, the calibration equation (6) of the EC-5 sensor was obtained in the same way as that of 5TE.

The EC-5 sensor showed an exponential response to the variation of soil moisture, as seen in Fig. 3, with a coefficient of determination of 0.991. Studies conducted by [8] corroborate the results for the coefficient of determination of this study.

Source: The Authors.

Figure 3 Volumetric soil moisture estimated by the EC-5 sensor before and after calibration and its correlation with the calculated volumetric soil moisture. 

For the EC-5 sensor, the soil moisture range estimated by it varied from 0.08 to 0.33 m3 m-3 while the calculated volumetric soil moisture range varied from 0.06 to 0.33 m3 m-3 (Fig. 3).

The volumetric soil moisture values determined by means of the factory calibration of the EC-5 sensor were overestimated compared to the calculated volumetric soil moisture values (located below the 1:1 line).

After calibration of the EC-5 sensor, the soil moisture values approached the 1:1 line, thus the soil moisture values of the ranges from 0.06 to 0.17 and from 0.25 to 0.33 m3 m-3 approached the 1:1 line; however, the soil moisture values that are within the range from 0.20 to 0.25 m3 m-3 were overestimated compared to the calculated soil moisture values, with an estimation error of 0.012.

On the other hand, Table 2 shows that for the 5TE sensor there is a significant difference between the accuracy of the volumetric soil moisture readings before and after calibration, with RMSE values of 0.065 and 0.012 cm3 cm-3, respectively, varying in the sensor precision scale from poor to moderate according to the classification presented by [23].

Table 2 Statistical indexes for sensors before and after calibration. 

Source: The Authors.

Table 2 shows that for the same sensor the values of performance coefficient (c) before and after calibration were 0.785 and 0.984, respectively, varying from very good (factory calibration) to excellent (laboratory calibration) according to the classification of [21]. Likewise, the indices of agreement varied from 0.791 to 0.991 before and after calibration, showing greater accuracy between the volumetric soil moisture values measured and estimated by the 5TE sensor.

Studies carried out by [22] obtained a coefficient of performance of 0.99 after calibration for a soil of sandy clay texture, evaluating a 5TE sensor in the estimation of soil moisture, corroborating the results obtained in the present study, which leads to the conclusion that the calibration applied to this sensor improved soil moisture readings using the calibration equation.

On the other hand, the RMSE for the EC-5 sensor decreased from 0.045 to 0.011 (Table 2). According to the classification of [23], the accuracy of this sensor is classified as moderate and was maintained before and after calibration. Likewise, according to the classification proposed by [21], for the EC-5 sensor the model used was classified as excellent before and after laboratory calibration.

As can be seen, soil moisture can be estimated with reliability after the specific calibration of the EC-5 sensor since it had a fairly high coefficient of determination, corresponding to 0.99.

4. Conclusions

In this study, two calibration equations were determined for the 5TE and EC-5 sensors in a sandy soil that can replace the factory equations, increasing the accuracy of soil moisture readings, contributing to the efficient management of irrigation water and thus contributing to the sustainability of agriculture.

The correlation between the volumetric soil moisture values corrected with the calibration equations obtained and those of the volumetric soil moisture obtained by the standard method was significant, with determination coefficients of 0.98 and 0.99 for the 5TE and EC-5 sensors, respectively.

The RMSE values for the 5TE sensor before and after calibration varied from 0.065 to 0.012 m3 m-3, respectively, as occurred with the EC-5 sensor, which had a low RMSE value (0.011), showing that these types of sensor promoted better estimation of volumetric soil moisture after being calibrated.

5. References

[1] Díaz, A.D.M., Vélez, S.J.E. y Rodríguez, H.P. Efecto de la aplicación de riego controlado en la producción y calidad del fruto de Pyrus communis L., cv. Triunfo de Viena. Acta Agronómica, 65(2), pp. 156-163, 2016. DOI: 10.15446/acag.v65n2.49650. [ Links ]

[2] Souza, C.F., Pires, R.C.M, Miranda, D.B.De and Varallo, A.C.T., Calibração de sondas FDR e TDR para a estimativa da umidade em dois tipos de solo. Irriga, 18(4), pp. 597-606, 2013. [ Links ]

[3] Kinzli, K., Manana, N. and Oad, R., A comparison of laboratory and field calibration of a soil moisture capacitance probe for various soils. Journal of Irrigation and Drainage Engineering, 138(4), pp. 310-321, 2012. DOI: 10.1061/(ASCE)IR.1943-4774.0000418. [ Links ]

[4] Carpena-Muñoz, R., Herbert, B., Klassen, W. and Dukes, M.D., Automatic soil moisture based drip irrigation for improving tomato production. Proceedings of the Florida State Horticultural Society, 116, pp. 80-85, 2003. [ Links ]

[5] Walker, J.P., Willgoose, G.R. and Kalma, J.D., In situ measurement of soil moisture: a comparison of technique. Journal of Hydrology, 293(1-4), pp. 85-99, 2004. DOI: 10.1016/j.jhydrol.2004.01.008. [ Links ]

[6] Liu, Z. and Xu, Q., Precision irrigation scheduling using ECH2O moisture sensors for lettuce cultivated in a soilless substrate culture. Water 10(5), pp. 549-559, 2018. DOI: 10.3390/w10050549. [ Links ]

[7] Czarnomski, N.M., Pypker, T.G., Licate, J. and Bon, B.J., Precision and accuracy of three alternative instruments for measuring soil water content in two forest soils of the Pacific Northwest. Canadian Journal of Forest Research, 35(8), pp. 1867-1876, 2005. DOI: 10.1139/x05-121. [ Links ]

[8] Kizito, F., Campell, C.S., Campell, G.S., Cobos, D.R., Teare, B.L., Carter, B. and Hopmans, J.W., Frequency, electrical conductivity and temperature analysis of a low-cost capacitance soil moisture sensor. Journal of Hydrology, 252(3-4), pp. 367-378, 2008. DOI: 10.1016/j.jhydrol.2008.01.021.2 [ Links ]

[9] Stacheder, M., Koeniger, F. and Schuhmann, R., New dielectric sensors and sensing techniques for soil and snow moisture measurements. Sensors, 9(4), pp. 951-2967, 2009. DOI: 10.3390/s90402951. [ Links ]

[10] Lima, J.E.F.W., Silva, E.M., Koide, S. e Santos, R.M., Avaliação do desempenho de sonda de capacitância no monitoramento da umidade de latossolos do cerrado em condições de campo. Revista Brasileira de Recursos Hídricos, 17(1), pp. 23-32, 2012. [ Links ]

[11] Sakaki, T., Limsuwat, A. and Illangasekare, T.H., A simple method for calibrating dielectric soil moisture sensors: laboratory validation in sands. Vadose Zone Journal, 10(2), pp. 706-715, 2011. DOI: 10.2136/vzj2010.0121. [ Links ]

[12] Decagon Devices. 5TE Operator’s Manual. Devices, Inc., Pullman, WA. [online]. Available at: https://www.metergroup.com/environment/products/ech2o-5te-electrical-conductivity/. [ Links ]

[13] Decagon Devices. EC-5 Operator’s Manual. Devices, Inc., Pullman, WA. [online]. Available at: https://www.metergroup.com/environment/products/ec-5-soil-moisture-sensor/Links ]

[14] Jabro, J.D., Stevens, W.B. and Iversen. W.M., Field performance of three real-time moisture sensors in sandy loam and clay loam soils. Archives of Agronomy and Soil Science, 64(7), pp. 930-938, 2017. DOI: 10.1080/03650340.2017.1393528. [ Links ]

[15] Loiskandl, W., Buchan, G.D., Sokol, W., Novak, V. and Himmelbauer, M., Calibrating electromagnetic short soil water sensors. Journal of Hydrology and Hydromechanics, 58(2), pp. 114-125, 2010. DOI: 10.2478/v10098-010-0011-3. [ Links ]

[16] Matula, S., Báťková, K. and Legese, W.L., Laboratory performance of five selected soil moisture sensors applying factory and own calibration equations for two soil media of different bulk density and salinity levels. Sensors, 16(11), pp. 1912-1934, 2016. DOI: 10.3390/s16111912. [ Links ]

[17] Cardenas-Lailhacar, B. and Dukes, M.D., Precision of soil moisture sensor irrigation controller under field conditions. Agricultural Water Management, 97(5), pp. 666-672, 2010. DOI: 10.1016/j.agwat.2009.12.009. [ Links ]

[18] EMBRAPA - Centro Nacional de Pesquisa de Solos. Sistema Brasileiro de Classificação de Solos. 3rd Ed. Rev. Ampl. Brasília, DF: Embrapa, 2013. 353 P. [ Links ]

[19] Stockle, C.O., Kjelgaard, J. and Bellocchi, G., Evaluation of estimated weather data for calculating Penman-Monteith reference evapotranspiration. Irrigation Science, 23(1), pp. 39-46, 2004. DOI: 10.1007/s00271-004-0091-0. [ Links ]

[20] Willmott, C.J., On the validation of models. Physical Geography, 1981, pp.184-194. [ Links ]

[21] Camargo, A.P. e Sentelhas, P.C., Avaliação do desempenho de diferentes métodos de estimativas da evapotranspiração potencial no Estado de São Paulo, Brasil. Revista Brasileira de Agrometeorologia, 5(1), pp. 89-97, 1997. [ Links ]

[22] Varble, J.L. and Chávez, J.L., Performance evaluation and calibration of soil water content and potential sensors for agricultural soils in eastern Colorado. Agricultural Water Management, 101(1), pp. 93-106, 2011. DOI: 10.1016/j.agwat.2011.09.007. [ Links ]

[23] Fares, A., Abbas, F., Domingos, M. and Mair, A., Improved calibration functions of three capacitance probes for the measurement of soil moisture in tropical soils. Sensors, 11(5), pp. 4858-4874, 2011. DOI: 10.3390/s110504858. [ Links ]

A.A.C. Jiménez, is BSc. in Agricultural Engineer from the National University of Loja, Ecuador (2015), is MSc. Agricultural Engineering from the Federal Rural University of Pernambuco, Brazil (2018). She has experience in the area of rural construction and water and soil conservation. ORCID: 0000-0001-7159-4385.

C.D.G.C. Almeida, received the BSc. in Agronomy from the Federal Rural University of the Semi-Arid Region, Brazil (1994), postgraduate in Agricultural Sciences from the Federal Rural University of Pernambuco (1998), MSc. in Agronomy (Soil Sciences) from the Federal Rural University of Pernambuco (1998) and PhD in Irrigation and Drainage from the “Luiz de Queiroz” College of Agriculture (ESALQ / USP) with training in the split-site doctorate at the University of Southern Queensland, Australia. She is a full professor at the Dom Agostinho Ikas Agricultural College (CODAI) of the Federal Rural University of Pernambuco (UFRPE). Permanent professor in the Postgraduate Program in Agricultural Engineering at UFRPE. ORCID: 0000-0001-6073-3853.

J.A. Santos Júnior, is a BSc. in Agricultural Engineer from the Federal University of Campina Grande, Brazil (2009), MSc. in Agricultural Engineering from the Federal University of Campina Grande and PhD in Agricultural Engineering from the Federal University of Campina Grande (2013). He is professor at the Department of Agricultural Engineering of the Federal Rural University of Pernambuco (UFRPE). He is a full professor in the Postgraduate Program in Agricultural Engineering of UFRPE. ORCID: 0000-0002-1656-7103.

C.S. dos Santos, is PhD in Agricultural Engineering from the Federal Rural University of Pernambuco, Brazil. MSc. in Agricultural Engineering from the Federal Rural University of Pernambuco (2015), holds a BSc. in Agricultural Engineering from the Federal University of Alagoas - Arapiraca Campus (2013) and Agricultural Technician by the Federal Agrotechnical School of Satuba-AL (2006). She has experience in the area of agronomy, water resources and agricultural engineering, with emphasis on irrigation management, acting mainly on the following topics: salinity, evaluation of irrigation systems, coefficient of horticultural and oil crops, agricultural reuse, semi-arid, conservation of water and soil. ORCID: 0000-0001-9784-3001.

How to cite: Jiménez, A. A. C, Almeida, C. D. G. C, Santos Júnior, J. A. and dos Santos, C. S, Calibration of two capacitive soil moisture sensors in Ultisol. DYNA, 87(213), pp. 75-79, April - June, 2020.

Received: October 04, 2018; Revised: September 16, 2019; Accepted: October 18, 2019

Creative Commons License The author; licensee Universidad Nacional de Colombia