Acessibilidade / Reportar erro

Evaluation of crop water status of melon plants in tropical semi-arid climate using thermal imaging1 1 Research developed at Instituto Federal de Educação, Ciência e Tecnologia do Ceará/Campus Sobral, Eixo Tecnológico de Recursos Naturais, Sobral, CE, Brazil

Imagens térmicas para avaliar o estado hídrico do meloeiro sob clima tropical semiárido

ABSTRACT

The objective of this study was to analyze the feasibility of using thermal images to estimate the water status of melon plants (Cucumis melo L.) in tropical semi-arid climates. The study was conducted in a randomized block design with a split-plot arrangement. The plots comprised of soil cover (with and without mulching), and subplots were constructed using five irrigation regimes (120, 100, 80, 60, and 40% crop evapotranspiration), with five replicates. The following variables were evaluated: canopy temperature (Tcanopy), leaf water potential, air temperature (Tair), soil moisture, crop yield, and thermal index (ΔT), which is defined as the difference between Tcanopy and Tair. ΔT exhibited high correlations with crop yield and water consumption, indicating that thermography is an efficient tool for identifying the water status of melon plants, which could be employed for proper irrigation scheduling under tropical semi-arid scenarios. Moreover, thermal images identified the beneficial effects of soil cover on leaf water status and crop yield, primarily under moderate deficit irrigation. These results demonstrate that mulching is essential for increasing melon yield and water productivity in tropical regions.

Key words:
Cucumis melo; dryland agriculture; irrigation management; leaf temperature; water stress

RESUMO

O objetivo principal deste trabalho foi avaliar o uso de imagens térmicas para estimar o estado hídrico de plantas de meloeiro (Cucumis melo L.) sob clima tropical semiárido. O estudo foi realizado em delineamento de blocos casualizados, em esquema de parcelas subdivididas. As parcelas foram compostas pela cobertura do solo (com e sem cobertura morta) e as subparcelas por cinco regimes de irrigação (120, 100, 80, 60 e 40% da evapotranspiração da cultura - ETc), com cinco repetições. Foram avaliadas as seguintes variáveis: temperatura do dossel (Tdossel), potencial hídrico foliar (Ψfolha), temperatura do ar (Tair), umidade do solo, produtividade da cultura e índice térmico (ΔT), sendo este definido como a diferença entre Tdossel e Tair. O ΔT apresentou altas correlações com a produtividade da cultura e o consumo de água da cultura, evidenciando que a termografia é uma ferramenta eficiente para identificar o estado hídrico das plantas de meloeiro e pode ser empregada para uma adequada programação da irrigação em cenários do semiárido tropical. Além disso, o uso de imagens térmicas também permitiu identificar os efeitos benéficos da cobertura do solo sobre o estado hídrico foliar e a produtividade das culturas, principalmente sob moderado déficit de irrigação (DI). Os resultados obtidos também demonstram que a cobertura morta é essencial para aumentar a produtividade do melão e a produtividade da água em regiões tropicais.

Palavras-chave:
Cucumis melo; agricultura de sequeiro; manejo de irrigação; temperatura foliar; estresse hídrico

HIGHLIGHTS:

Deficit irrigation strategies effectively increase water productivity.

Thermal images may be used as an irrigation management strategy.

Mulching is essential for increasing melon fruit yield and water use efficiency.

Introduction

Globally, melon (Cucumis melo L.) is one of the most appreciated fruits by consumers. In Brazil, melon production is concentrated in the northeast region, primarily in the states of Ceará and Rio Grande do Norte, which accounts for over 95% of national production. The primary melon-producing regions have semi-arid characteristics, thereby providing the best climatic conditions. However, they are prone to droughts, and water shortages are a critical issue for sustainable agriculture. Thus, new management strategies are needed for improving irrigation water productivity to support continued production (Cymes et al., 2022Cymes, I.; Dragańska, E.; Brodziński, Z. Potential possibilities of using groundwater for crop irrigation in the context of climate change. Agriculture, v.12, p.1-14, 2022. https://doi.org/10.3390/agriculture12060739
https://doi.org/10.3390/agriculture12060...
).

Deficit irrigation (DI) strategies save water without significantly affecting crop yields (Jahan & Nassiri-Mahallati, 2022Jahan, M.; Nassiri-Mahallati, M. Modeling the response of sesame (Sesamum indicum L.) growth and development to climate change under deficit irrigation in a semi-arid region. PLOS Climate, v.1, p.1-20, 2022. https://doi.org/10.1371/journal.pclm.0000003
https://doi.org/10.1371/journal.pclm.000...
). However, effective monitoring of water status in plants during their cultivation cycle is crucial (Poblete-Echeverría et al., 2016Poblete-Echeverría, C.; Sepúlveda-Reyes, D.; Ortega-Farias, S.; Zúñiga, M.; Fuentes, S. Plant water stress detection based on aerial and terrestrial inferred thermography: a study case from vineyard and olive orchard. Acta Horticulturae, n.1112, p.141-146, 2016. https://doi.org/10.17660/ActaHortic.2016.1112.20
https://doi.org/10.17660/ActaHortic.2016...
; Carrasco-Benavides et al., 2020bCarrasco-Benavides, M.; Espinoza Meza, S.; Olguín-Cáceres, J.; Muñoz-Concha, D.; von Bennewitz, E.; Ávila-Sánchez, C.; Ortega-Farías, S. Effects of regulated post-harvest irrigation strategies on yield, fruit quality and water productivity in a drip-irrigated cherry orchard. New Zealand Journal of Crop and Horticultural Science, v.48, p.1-20, 2020b. https://doi.org/10.1080/01140671.2020.1721544
https://doi.org/10.1080/01140671.2020.17...
). This measurement has been successfully performed using thermal images, which can indicate and detect the effects of water stress in crops (Carrasco-Benavides et al., 2020aCarrasco-Benavides, M.; Antunez-Quilobrán, J.; Baffico-Hernández, A.; Ávila-Sánchez, C.; Ortega-Farías, S.; Espinoza, S.; Gajardo, J.; Mora, M.; Fuentes, S. Performance assessment of thermal infrared cameras of different resolutions to estimate tree water status from two cherry cultivars: An alternative to midday stem water potential and stomatal conductance. Sensors, v.20, p.1-18, 2020a. https://doi.org/10.3390/s20123596
https://doi.org/10.3390/s20123596...
; Silveira et al., 2020Silveira, J. M. de C.; Cândido, B. M.; Silva, A. L. B. de O.; Pavão, G. C.; Xavier, M. A.; Pires, R. C. de M. Uso de imagens multiespectrais e termográficas para monitoramento das condições hídricas da cana-de-açucar. Irriga, v.25, p.689-696, 2020. https://doi.org/10.15809/irriga.2020v25n4p689-696
https://doi.org/10.15809/irriga.2020v25n...
).

Thermography is a potential tool for estimating canopy temperature (Tcanopy), which can be used as an indicator of water stress (Costa et al., 2020Costa, J. de O; Coelho, R. D.; Barros, T. H. da S.; Fraga Júnior, E. F.; Fernandes, A. L. T. Canopy thermal response to water deficit of coffee plants under drip irrigation. Irrigation and Drainage, v.69, p.472-482, 2020. https://doi.org/10.1002/ird.2429
https://doi.org/10.1002/ird.2429...
; Gomes et al., 2021Gomes, F. F. L.; Lima, J. S.; Pinheiro Neto, L. G. Aplicações da termografia na agricultura. Agrarian Academy, v.8, p.91-105, 2021. https://doi.org/10.18677/Agrarian_Academy_2021B9
https://doi.org/10.18677/Agrarian_Academ...
). When combined with automated image analysis, it becomes highly accurate and can measure the relative temperature instead of the real temperature. New devices with high precision and low cost have increased the possibility of using this technique (García-Tejero et al., 2018García-Tejero, I. F.; Rubio, A. E.; Viñuela, I.; Hernández, A.; Gutiérrez-Gordillo, S.; Rodríguez-Pleguezuelo, C. R.; Durán-Zuazo, V. H. Thermal imaging at plant level to assess the crop-water status in almond trees (cv. Guara) under deficit irrigation strategies. Agricultural Water Management , v.208, p.176-186, 2018. https://doi.org/10.1016/j.agwat.2018.06.002
https://doi.org/10.1016/j.agwat.2018.06....
). Imaging techniques using infrared thermography can detect minimal changes in temperature and are widely used in agriculture (Ali et al., 2022Ali, M. M.; Hashim, N.; Aziz, S. A.; Lasekan, O. Quality prediction of different pineapple (Ananas comosus) varieties during storage using infrared thermal imaging technique. Food Control, v.138, p.1-9, 2022. https://doi.org/10.1016/j.foodcont.2022.108988
https://doi.org/10.1016/j.foodcont.2022....
).

With an increase in melon cultivated areas, global growth rates related to water demand have challenged researchers and producers, who seek to increase the efficiency of water use, reduce losses, and more importantly, increase the productivity of vegetable crops (Jacobi & Grandisoli, 2017Jacobi, P. R.; Grandisoli, E. Água e sustentabilidade: desafios, perspectivas e soluções. 1.ed. São Paulo: Reconectta, 2017. 110p). The objective of this study was to analyze the feasibility of using thermal images to estimate the water status of melon plants (Cucumis melo L.) in a tropical semi-arid climate.

Material and Methods

The experiment was conducted between July and September 2018 in the city of Sobral (03° 41’ 10” S, 40° 20’ 59” W, altitude 69 m), Ceará, Brazil. The local climate is BSw’h’, hot, and semi-arid based on the Köppen classification, with annual sunlight, a mean air temperature (Tair), and relative air humidity of 2,563 h, 28.3 °C, and 68%, respectively. The dry season lasts 7-8 months, and the rainy season is concentrated between January to May, with an average annual rainfall of 854 mm (INMET, 2018INMET - Instituto Nacional de Meteorologia. 2018. Available on: <Available on: http://www.inmet.gov.br >. Accessed on: July 2018.
http://www.inmet.gov.br...
). The soil in the area is classified as Entisol (United States, 2014United States. Soil Research Team. Keys to soil taxonomy. 12.ed. USDA NRCS, 2014. Available on: <Available on: http://www.nrcs.usda.gov/wps/portal/nrcs/main/soils/survey/ >. Accessed on: Nov. 2022.
http://www.nrcs.usda.gov/wps/portal/nrcs...
). The physicochemical characteristics of the soil samples collected at 0-0.20 m depths of soil layer are shown in Table 1.

Table 1
Chemical and physical characteristics of topsoil (0-0.20 m) of the experimental area

Seeds of yellow melon hybrid Goldex F1 were sown in polystyrene trays, with 128 cells each, using cattle manure and washed sand at 1:1 ratio as the substrate (v/v). After emergence, the seedlings were kept for approximately 7 days in the trays, a period needed for the development of two true leaves, and then transplanted to the field.

A randomized complete block design (with a split-plot arrangement) was used with five replicates per irrigation regime (IR). The plots comprised of soil cover (with and without plant material mulching), and subplots were formed using five IRs. The planting spacing was 2.0 × 0.5 m, and each subplot had 10 plants. Mulching was performed with dry leaves of carnauba [Copernicia prunifera (Miller) H. E. Moore] residue by keeping a 3.0 cm layer around the melon plants. This is a by-product of carnauba wax, which is extremely abundant in the region. Silva et al. (2014Silva, P. S. L.; Silva, P. I. B.; Oliveira, V. R.; Silva, K. E. F.; Monteira, A. L. Doses of carnauba leaf residues and corn yield. Horticultura Brasileira , v.32, p.74-79, 2014. https://doi.org/10.1590/S0102-05362014000100012
https://doi.org/10.1590/S0102-0536201400...
) reported the following chemical composition of this leaf residue: N = 27.13 g kg-1, P = 1.46 g kg-1, K = 5.41 g kg-1, Ca = 4.98 g kg-1, Mg = 2.88 g kg-1, Na = 1.4 g kg-1, B = 50.54 mg kg-1, Fe = 47.78 mg kg-1, Cu = 0.59 mg kg-1, Mn = 2.07 mg kg-1, and Zn = 0.56 mg kg-1. The five IRs were defined as a percentage of crop evapotranspiration (ETC) as follows: (IR1, 120%; IR2, 100%; IR3, 80%; IR4, 60%; and IR5, 40% ETC). The IRs were defined through the inclusion of water deficits that would not kill the plants and were maintained at higher than 100% to ensure adjustments of the crop coefficient (Kc) for the region. Spraying of pesticides (Mancozeb, Cyromazine, and Acetamiprid) registered with the Ministério da Agricultura, Pecuária e Abastecimento was performed to control pests and diseases.

The nutritional requirements of the crops were 140 kg ha-1 N, 300 kg ha-1 K2O, and 240 kg ha-1 P2O5. The N, K, and P sources used were urea, potassium chloride, and triple superphosphate, respectively. Top dressing was performed via daily fertigation using the Venturi system. A drip irrigation system was employed with PCJ-CNL self-compensating in-line button drippers at a mean pressure and nominal flow of 150 kPa and 4.2 L h-1, respectively, with one dripper per plant. Drippers were previously evaluated under laboratory conditions using Christiansen’s and distribution uniformity coefficients, which presented values of 91% and 84%, respectively.

The daily amount of water applied was determined by the ETC and estimations by the means of reference evapotranspiration (ET0) and Ksc (Allen et al., 1998Allen, R. G.; Pereira, L. S.; Raes, D.; Smith, M. Crop evapotranspiration: Guidelines for computing crop water requirements. Rome: Food and Agriculture Organization, 1998. 300p. Drainage and Irrigation Paper, 56). The ET0 was estimated by direct readings of daily evaporation in Class A pan, which was installed next to the experimental area by adopting a tank coefficient of 0.72, according to local climate conditions. The Ksc used were 0.50, 0.80, 1.05, and 0.75, as recommended by Doorenbos & Pruitt (1977Doorenbos, J.; Pruitt, J. O. Crop water requirement. Rome: Food and Agriculture Organization , 1977. 144p. Irrigation and Drainage Paper, 24) and Doorenbos & Kassam (1994Doorenbos, J.; Kassam, A. H. Effect of water on crop yield. Rome: Food and Agriculture Organization , 1994. 306p. Irrigation and Drainage Paper, 33), corresponding to periods of vegetative growth (35 days), flowering (15 days), and fruit production (10 days) and maturation (7 days), respectively.

The soil cover factor for localized irrigation was calculated using a ruler to measure the dimensions of the melon plant branches in the cross section of the planting rows by dividing the value of the row spacing. The amount of water applied in each irrigation treatment throughout the plant cycle was 386.69, 322.24, 257.79, 193.34, and 128.89 mm, for IR1, IR2, IR3, IR4, and IR5, respectively.

Soil moisture was determined in samples collected at 0-0.20 m soil layer at 15, 30, 45, and 60 days after transplanting (DAT). Soil moisture (g g-1) was determined in three replicates for each treatment using the gravimetric method (Jarrell et al., 1999Jarrell, W. M.; Armstrong, D. E.; Grigal, D. F.; Kelley, E. F.; Monger, H. C.; Wedin, D. A. Soil water and temperature status. In: Robertson, G. P.; Coleman, D. C.; Bledsoe, C. S.; Sollins, P. (eds.). Standard soil methods for long-term ecological research. New York: Oxford University Press, 1999. 18p.).

Leaf water potential (Ψleaf) was measured at 15, 30, 45, and 60 DAT at 5:00 a.m., according to Scholander et al. (1965Scholander, P. F.; Hammel, H. T.; Bradstreet, E. D.; Hemmingsen, E. A. Sap pressure in vascular plants. Science, v.148, p.339-346, 1965. https://www.science.org/doi/10.1126/science.148.3668.339
https://www.science.org/doi/10.1126/scie...
), using a pressure camera (Model 3115, Soil Moisture). The third leaf, counted from the apex of the stem (i.e., the youngest mature leaf), was chosen for the measurement. Three replicates were used for each treatment.

Tcanopy was measured using an FLIR I40 thermal camera on the same dates at which the soil moisture and Ψleaf were measured. Thus, on days when the Ψleaf and soil moisture were evaluated, the Tcanopy of the melon crop was also measured. The device instantaneously identifies temperature with a thermal sensitivity of 0.06 °C (thermal image resolution: 80 × 60; multispectral image resolution: 320 × 240). Thermal images were captured at distances between 0.50 and 1.0 m from the plants. Images were acquired manually using a camera on a tripod for better image adjustment. All images were captured between 8:00 and 9:00 a.m. (Brazilian Time; GMT-3) and treated using FLIR QuickReport software.

To quantify the average Tcanopy, three images of each plant were captured from three plants per treatment. Tair was monitored using a HOBO U12-012 thermohygrometer, with a data logger installed next to the experimental area. Tcanopy and Tair data were used to calculate the thermal index (ΔT; Tcanopy - Tair) according to Jackson et al. (1981Jackson, R. D.; Idso, S. B.; Reginato, R. J.; Pinter, P. J. Canopy temperature as a crop water stress indicator. Water Resources Research, v.17, p.1133-1138, 1981. http://dx.doi.org/10.1029/WR017i004p01133
http://dx.doi.org/10.1029/WR017i004p0113...
).

Fruit yield was determined for each experimental unit by extrapolating the values in t ha-1. Only fruits with commercial characteristics (Campelo et al., 2014Campelo, A. R.; Azevedo, B. M.; Nascimento Neto, J. R.; Viana, T. V. A.; Pinheiro Neto, L. G.; Lima, R. H. Manejo da cultura do melão submetida a frequências de irrigação e fertirrigação com nitrogênio. Horticultura Brasileira, v.32, p.138-144, 2014. https://doi.org/10.1590/S0102-05362014000200003
https://doi.org/10.1590/S0102-0536201400...
) were counted, and all fruits with defects were eliminated. The fruits were sorted and classified manually at 55 DAT, when all fruits had a yellow color and characteristics suitable for commercialization.

The data were subjected to analysis of variance (F-test), Tukey’s test (p ≤ 0.05), and regression analysis. The models were chosen based on the significance of the regression coefficients using the F-test (p ≤ 0.05) and the highest value of the coefficient of determination.

Results and Discussion

Figure 1 shows the mean values of Tcanopy, Ψleaf, and soil moisture measured at 15, 30, 45, and 60 DAT in the soil with and without mulching. Figures 1A and B exhibit an increase in Tcanopy up to 45 DAT, followed by a decrease, except in the treatments with 120% ETc. Because the water consumption in melon plant progressively increases over its phenological cycle, it is likely that most IRs did not meet the water requirements of the crop, particularly during the period of greatest demand. This assumption is based on the increases in Tcanopy, which usually occur when soil water availability is low and stomatal conductance decreases. According to Dejonge et al. (2015Dejonge, K. C.; Taghvaeian, S.; Thomas, J.; Trout, T. J.; Louise, H.; Comas, L. H. Comparison of canopy temperature-based water stress indices for maize. Agricultural Water Management, v.156, p.51-62, 2015. https://doi.org/10.1016/j.agwat.2015.03.023
https://doi.org/10.1016/j.agwat.2015.03....
), Tcanopy increases when solar radiation is absorbed; however, the plant is cooled when this energy is used for leaf transpiration.

Figure 1
Variations in mean canopy temperature (Tcanopy) (A and B), leaf water potential (Ψleaf) (C and D), and soil moisture (E and F) at 15, 30, 45, and 60 days after transplanting (DAT) of melons under different irrigation regimes (IR1- 120%, IR2- 100%, IR3- 80%, IR4- 60%, and IR5- 40% ETc) in soil with (A, C, and E) and without mulching (B, D, and F)

In general, Tcanopy decreased when water availability was elevated. When the plant was under well-watered conditions, Tcanopy remained between 18 and 20 °C, with small changes throughout the crop cycle (Figures 1A and B). These results are similar to those obtained for other crops, such as vineyards and olive orchards (Poblete-Echeverría et al., 2016Poblete-Echeverría, C.; Sepúlveda-Reyes, D.; Ortega-Farias, S.; Zúñiga, M.; Fuentes, S. Plant water stress detection based on aerial and terrestrial inferred thermography: a study case from vineyard and olive orchard. Acta Horticulturae, n.1112, p.141-146, 2016. https://doi.org/10.17660/ActaHortic.2016.1112.20
https://doi.org/10.17660/ActaHortic.2016...
), where water status was estimated in plants using thermal information. Similarly, other studies have reported significant correlations between thermal information and physiological variables related to plant water status (Ramírez-Cuesta et al., 2022Ramírez-Cuesta, J. M.; Ortuño, M. F.; Gonzalez-Dugo, V.; Zarco-Tejada, P. J.; Parra, M.; Rubio-Asensio, J. S.; Intrigliolo, D. S. Assessment of peach trees water status and leaf gas exchange using on-the-ground versus airborne-based thermal imagery. Agricultural Water Management , v.267, p.1-14, 2022. https://doi.org/10.1016/j.agwat.2022.107628
https://doi.org/10.1016/j.agwat.2022.107...
).

In general, higher Tcanopy values were observed in IR4 and IR5 treatments. This could be related to the greater water restriction imposed by these treatments during the melon crop cycle. Under these conditions, the plant closed its stomata and reduced transpiration, which increased the Tcanopy. García-Tejero et al. (2018García-Tejero, I. F.; Rubio, A. E.; Viñuela, I.; Hernández, A.; Gutiérrez-Gordillo, S.; Rodríguez-Pleguezuelo, C. R.; Durán-Zuazo, V. H. Thermal imaging at plant level to assess the crop-water status in almond trees (cv. Guara) under deficit irrigation strategies. Agricultural Water Management , v.208, p.176-186, 2018. https://doi.org/10.1016/j.agwat.2018.06.002
https://doi.org/10.1016/j.agwat.2018.06....
) and Pou et al. (2014Pou, A.; Diago, M. P.; Medrano, H.; Baluja, J.; Tardaguila, J. Validation of thermal indices for water status identification in grapevine. Agricultural Water Management , v.134, p.60-72, 2014. https://doi.org/10.1016/j.agwat.2013.11.010
https://doi.org/10.1016/j.agwat.2013.11....
) conducted studies on almonds and grapevines, respectively, using thermography and observed that plants subjected to water restriction conditions had higher Tcanopy. However, variations in meteorological variables can also directly influence the Tcanopy data (Javadian et al., 2022Javadian, M.; Smith, W. K.; Lee, K.; Knowles, J. F.; Scott, R. L.; Fisher, J. B.; Moore, D. J. P.; van Leeuwen, W. J. D.; Barron-Gafford, G.; Behrangi, A. Canopy temperature is regulated by ecosystem structural traits and captures the ecohydrologic dynamics of a semiarid mixed conifer forest site. Journal of Geophysical Research: Biogeosciences, v.127, p.1-15, 2022. https://doi.org/10.1029/2021JG006617
https://doi.org/10.1029/2021JG006617...
).

Plants growing in soil without mulching exhibited mean Tcanopy values of 1.3-3.6 °C, which were higher than those observed in the treatments with soil cover. It is likely that mulching led to lower evaporation, and consequently, a higher amount of water remained available to plants. The use of mulch resulted in soil water conservation (Figures 1E and F), as demonstrated in other studies (El-Beltagi et al., 2022El-Beltagi, H. S.; Basit, A.; Mohamed, H. I.; Ali, I.; Ullah, S.; Kamel, E. A. R.; Shalaby, T. A.; Ramadan, K. M. A.; Alkhateeb, A. A.; Ghazzawy, H. S. Mulching as a sustainable water and soil saving practice in agriculture: a review. Agronomy, v.12, p.1-31, 2022. https://doi.org/10.3390/agronomy12081881
https://doi.org/10.3390/agronomy12081881...
), helping plants maintain a lower Tcanopy (Ni et al., 2019Ni, J.; Cheng, Y.; Wang, Q.; Ng, C. W. W.; Garg, A. Effects of vegetation on soil temperature and water content: field monitoring and numerical modeling. Journal of Hydrology, v.571, p.494-502, 2019. https://doi.org/10.1016/j.jhydrol.2019.02.009
https://doi.org/10.1016/j.jhydrol.2019.0...
). Moreover, this technique contributes to increased crop yields and water use efficiency, as demonstrated in a study conducted in the semi-arid regions of China (Wang et al., 2016Wang, Y.; Li, X.; Zhu, J.; Fan, C.; Kong, X.; Turner, N. C.; Siddique, K. H. M.; Li, F. Multi-site assessment of the effects of plastic-film mulch on dry land maize productivity in semiarid areas in China. Agricultural and Forest Meteorology, v.220, p.160-169, 2016. https://doi.org/10.1016/j.agrformet.2016.01.142
https://doi.org/10.1016/j.agrformet.2016...
).

The mean values of Ψleaf in IR1 and IR2 were -0.079 and -0.090 MPa, respectively, indicating that the water supply was adequate in these treatments, with plants exhibiting a good water status in the leaves. However, Ψleaf decreased in IR3, IR4, and IR5, especially at 45 DAT (Figures 1C and D), which is the most critical stage in the melon crop cycle, with flowering and early fruit production. These results indicated that the water deficit intensified during this period, mostly in treatments without mulching, which also exhibited lower soil moisture (Figures 1E and F). According to García-Tejero et al. (2016García-Tejero, I. F.; Costa, J. M.; Egipto, R.; Durán-Zuazo, V. H.; Lima, R. S. N.; Lopes, C. M.; Chaves, M. M. Thermal data to monitor crop-water status in irrigated Mediterranean viticulture. Agricultural Water Management , v.176, p.80-90, 2016. https://doi.org/10.1016/j.agwat.2016.05.008
https://doi.org/10.1016/j.agwat.2016.05....
), Ψleaf is a good indicator of leaf water status. It is particularly important to indicate the variations and water stress levels of plants under poor irrigation or different soil management (for example, soil cover), as shown in this study (Figure 1).

ΔT was not affected by soil cover in the treatment with the highest DI (60% less water; IR5) in all evaluations (Figure 2). This treatment also demonstrated the highest and most positive values of ΔT, indicating considerable effects on the leaf water status. In contrast, the soil cover results were identical in the treatments with high irrigation depths (IR1 and IR2) at 15, 30, and 60 DAT, but the values were always negative, suggesting improvements in crop water status (Figures 2A, B, and D). However, the best contributions of mulching were observed in treatments with moderate DI (up to 40% less water; IR3 and IR4), when plants exhibited lower thermal indices at all evaluations compared with the values observed in plants grown in soil without mulching. Thus, soil cover is efficient in maintaining soil moisture values for a longer period, which reduces the risk of water stress in crops (Wang et al., 2016Wang, Y.; Li, X.; Zhu, J.; Fan, C.; Kong, X.; Turner, N. C.; Siddique, K. H. M.; Li, F. Multi-site assessment of the effects of plastic-film mulch on dry land maize productivity in semiarid areas in China. Agricultural and Forest Meteorology, v.220, p.160-169, 2016. https://doi.org/10.1016/j.agrformet.2016.01.142
https://doi.org/10.1016/j.agrformet.2016...
). Our results demonstrated the importance of mulching to ensure soil water availability at the time of highest demand, even with high irrigation depths (100%), and for treatments with moderate DI (primarily up to 20%), helping plants maintain the leaf water status, as shown by thermal images and Ψleaf analyses (Figure 1).

Figure 2
Means of thermal index (ΔT; Tcanopy - Tair) as a function of irrigation regimes (IRs) and soil cover (with and without of mulching) at 15 (A), 30 (B), 45 (C), and 60 (D) days after transplanting (DAT) of melons at different IRs (IR1- 120%, IR2- 100%, IR3- 80%, IR4- 60%, and IR5- 40% ETc)

ΔT (Tcanopy - Tair) is recommended to estimate crop water stress using thermometric images because of its feasibility (García-Tejero et al., 2018García-Tejero, I. F.; Rubio, A. E.; Viñuela, I.; Hernández, A.; Gutiérrez-Gordillo, S.; Rodríguez-Pleguezuelo, C. R.; Durán-Zuazo, V. H. Thermal imaging at plant level to assess the crop-water status in almond trees (cv. Guara) under deficit irrigation strategies. Agricultural Water Management , v.208, p.176-186, 2018. https://doi.org/10.1016/j.agwat.2018.06.002
https://doi.org/10.1016/j.agwat.2018.06....
). Urban et al. (2017Urban, J.; Ingwers, M. W.; Mcguire, M. A.; Teskey, R. O. Increase in leaf temperature opens stomata and decouples net photosynthesis from stomatal conductance in Pinus taeda and Populus deltoides x nigra. Journal of Experimental Botany, v.68, p.1757-1767, 2017. https://doi.org/10.1093/jxb/erx052
https://doi.org/10.1093/jxb/erx052...
) reported temperature differences up to 9 °C between plants grown in dry and moist soils. Our results also indicated significant differences in ΔT between IR treatments and between treatments with and without soil cover (Figure 2). Furthermore, the use of mulching caused better plant development conditions and increased water use efficiency through lower water loss by evaporation, leading to a greater water storage volume for longer periods (Li et al., 2013Li, S. X.; Wang, Z. H.; Li, S. Q.; Gao, Y. J.; Tian, X. H. Effect of plastic sheet mulch, wheat straw mulch, and maize growth on water loss by evaporation in dryland areas of China. Agricultural Water Management , v.116, p.39-49, 2013. https://doi.org/10.1016/j.agwat.2012.10.004
https://doi.org/10.1016/j.agwat.2012.10....
).

ΔT was negatively correlated with soil moisture (Figure 3A) and ETC (Figure 3B). This implies that under ideal soil moisture conditions, there are high transpiration rates in the whole plant, markedly increasing and reducing Tcanopy. In contrast, low soil moisture not only affects plant water relations by reducing leaf water content and cell turgor but also affects stomatal conductance, leaf transpiration, Tcanopy, and carbon assimilation rates.

Figure 3
Relationships between thermal index (ΔT; Tcanopy-Tair) and soil moisture (A), ΔT and crop evapotranspiration (ETc) (B), ΔT and leaf water potential (Ψleaf) (C), and fruit yield and ΔT (D) of melons, as a function of soil cover

Figure 3C illustrates the variations in ΔT (Tcanopy - Tair) as a function of Ψleaf, which exhibits a downward linear trend. Positive ΔT values occurred at lower Ψleaf, implying that water restriction made the Tcanopy higher than the Tair. These results indicated that under these conditions, melon leaves presented lower transpiration rates and consequently, Tcanopy increased, implying that the crop is experiencing water stress. Similar results were reported by Mira-García et al. (2022Mira-García, A. B.; Conejero, W.; Vera, J.; Ruiz-Sánchez, M. C. Water status and thermal response of lime trees to irrigation and shade screen, Agricultural Water Management , v.272, p.1-14, 2022. https://doi.org/10.1016/j.agwat.2022.107843
https://doi.org/10.1016/j.agwat.2022.107...
).

Figure 3D demonstrates that when ΔT was approximately -3.2 °C, crop yield reached values of approximately 44.0 t ha-1. However, when ΔT (Tcanopy - Tair) was approximately 0.0 °C, fruit yield decreased to 34.0 t ha-1. On comparing obtained crop yields with the highest and lowest ΔT (Tcanopy - Tair) values, a difference of 30.47 t ha-1 was observed, corresponding to 69% production loss. That is, when the thermal indices were higher, melon fruit yield decreased significantly. The reduction in fruit yield under high Tcanopy supports the findings of Dejonge et al. (2015Dejonge, K. C.; Taghvaeian, S.; Thomas, J.; Trout, T. J.; Louise, H.; Comas, L. H. Comparison of canopy temperature-based water stress indices for maize. Agricultural Water Management, v.156, p.51-62, 2015. https://doi.org/10.1016/j.agwat.2015.03.023
https://doi.org/10.1016/j.agwat.2015.03....
), who discussed the results obtained in their study with maize.

The thermal images presented in Figures 4 and 5 show the impact of the treatments on ΔT and melon yield. In the mulching treatment (Figure 4), images of plants under no DI (IR1 and IR2) presented values lower than -3.3 °C, with a yield of 44.0 t ha-1. In contrast, ΔT was higher in all observations for plants under low irrigation depths (IR4 and IR5), and the yield reached only 24.78 and 13.20 t ha-1, respectively. These results exhibited that the use of infrared thermal images aided in recognizing the differences in the water availability for full and DI treatments, as was also demonstrated for vineyards and olive orchards (Poblete-Echeverría et al., 2016Poblete-Echeverría, C.; Sepúlveda-Reyes, D.; Ortega-Farias, S.; Zúñiga, M.; Fuentes, S. Plant water stress detection based on aerial and terrestrial inferred thermography: a study case from vineyard and olive orchard. Acta Horticulturae, n.1112, p.141-146, 2016. https://doi.org/10.17660/ActaHortic.2016.1112.20
https://doi.org/10.17660/ActaHortic.2016...
; García-Tejero et al., 2018García-Tejero, I. F.; Rubio, A. E.; Viñuela, I.; Hernández, A.; Gutiérrez-Gordillo, S.; Rodríguez-Pleguezuelo, C. R.; Durán-Zuazo, V. H. Thermal imaging at plant level to assess the crop-water status in almond trees (cv. Guara) under deficit irrigation strategies. Agricultural Water Management , v.208, p.176-186, 2018. https://doi.org/10.1016/j.agwat.2018.06.002
https://doi.org/10.1016/j.agwat.2018.06....
).

Figure 4
Thermal images, thermal index (ΔT), and fruit yield of melons under different irrigation regimes (IR; IR1-120, IR2-100, IR3-80, IR4-60, and IR5-40% ETc) grown in soil with mulching. The values of ΔT (Tcanopy - Tair) correspond to the average of nine measurements taken for each treatment at each date of evaluation

Figure 5
Thermal images, thermal index (ΔT), and fruit yield of melons under different irrigation regimes (IR; IR1-120, IR2-100, IR3-80, IR4-60, and IR5-40% ETc) grown in soil without mulching. The values of ΔT (Tcanopy - Tair) correspond to the average of nine measurements taken for each treatment at each date of evaluation

In contrast, in treatments without mulching (Figure 5), ΔT ranged from -3.0 °C under full irrigation to +3.8 °C for the most stressed treatment. On comparing the data in Figures 4 and 5, it can be observed that plants cultivated in the soil with mulch exhibited better physiological performance and higher production. Ni et al. (2019Ni, J.; Cheng, Y.; Wang, Q.; Ng, C. W. W.; Garg, A. Effects of vegetation on soil temperature and water content: field monitoring and numerical modeling. Journal of Hydrology, v.571, p.494-502, 2019. https://doi.org/10.1016/j.jhydrol.2019.02.009
https://doi.org/10.1016/j.jhydrol.2019.0...
) stated that the use of soil cover results in lower Tcanopy and favors the assimilation of nutrients, root development, water retention, and soil aeration. Our results indicated that the analysis of thermometric images accurately demonstrated the positive impact of this technique, exhibiting clear differences between treatments with and without soil cover. These beneficial effects of soil cover, confirmed by thermal images, have been reported by other authors (García-Tejero et al., 2016García-Tejero, I. F.; Costa, J. M.; Egipto, R.; Durán-Zuazo, V. H.; Lima, R. S. N.; Lopes, C. M.; Chaves, M. M. Thermal data to monitor crop-water status in irrigated Mediterranean viticulture. Agricultural Water Management , v.176, p.80-90, 2016. https://doi.org/10.1016/j.agwat.2016.05.008
https://doi.org/10.1016/j.agwat.2016.05....
).

Our results exhibited that mulching was efficient in reducing ΔT and increasing fruit yield, especially in treatments with moderate water stress. For treatments with 20% and 40% DI, there were increases of 30.6% and 52%, respectively, in crop yields of plots with mulching compared with those of without mulching. Moreover, there was a beneficial effect of soil cover in the treatment without water deficit (100% ETc), with approximately 24% increase in crop yield. This treatment may have presented a deficit in the stage of greater water consumption, as evidenced by Ψleaf and Tcanopy data, especially in soil without mulching (Figure 1). The use of a nonregional Kc explains, at least in part, the water deficit observed in this treatment. However, the treatment with 120% ETc application exhibited no differences between the treatments with and without mulching, indicating that additional 20% water was necessary to meet the demand of crop grown without soil cover.

Conclusions

  1. Thermal index showed high correlation with crop yield, leaf water potential, and crop water consumption.

  2. Thermal images are efficient for identifying the water status of melon plants and may be employed as an irrigation management strategy in tropical semi-arid climates.

  3. Thermal images identified the beneficial effects of soil cover on leaf water status and crop yield, especially under moderate deficit irrigation.

  4. Mulching is essential for increasing the fruit yield and water use efficiency in melon cultivations in semi-arid tropical regions.

Acknowledgements

The authors acknowledge the “Conselho Nacional de Desenvolvimento Científico e Tecnológico,” “Instituto Federal de Educação, Ciência e Tecnologia do Estado do Ceará (Campus Sobral),” and “Coordenação de Aperfeiçoamento de Pessoal de Nível Superior,” Brazil, for their financial support in this research.

Literature Cited

  • Ali, M. M.; Hashim, N.; Aziz, S. A.; Lasekan, O. Quality prediction of different pineapple (Ananas comosus) varieties during storage using infrared thermal imaging technique. Food Control, v.138, p.1-9, 2022. https://doi.org/10.1016/j.foodcont.2022.108988
    » https://doi.org/10.1016/j.foodcont.2022.108988
  • Allen, R. G.; Pereira, L. S.; Raes, D.; Smith, M. Crop evapotranspiration: Guidelines for computing crop water requirements. Rome: Food and Agriculture Organization, 1998. 300p. Drainage and Irrigation Paper, 56
  • Campelo, A. R.; Azevedo, B. M.; Nascimento Neto, J. R.; Viana, T. V. A.; Pinheiro Neto, L. G.; Lima, R. H. Manejo da cultura do melão submetida a frequências de irrigação e fertirrigação com nitrogênio. Horticultura Brasileira, v.32, p.138-144, 2014. https://doi.org/10.1590/S0102-05362014000200003
    » https://doi.org/10.1590/S0102-05362014000200003
  • Carrasco-Benavides, M.; Antunez-Quilobrán, J.; Baffico-Hernández, A.; Ávila-Sánchez, C.; Ortega-Farías, S.; Espinoza, S.; Gajardo, J.; Mora, M.; Fuentes, S. Performance assessment of thermal infrared cameras of different resolutions to estimate tree water status from two cherry cultivars: An alternative to midday stem water potential and stomatal conductance. Sensors, v.20, p.1-18, 2020a. https://doi.org/10.3390/s20123596
    » https://doi.org/10.3390/s20123596
  • Carrasco-Benavides, M.; Espinoza Meza, S.; Olguín-Cáceres, J.; Muñoz-Concha, D.; von Bennewitz, E.; Ávila-Sánchez, C.; Ortega-Farías, S. Effects of regulated post-harvest irrigation strategies on yield, fruit quality and water productivity in a drip-irrigated cherry orchard. New Zealand Journal of Crop and Horticultural Science, v.48, p.1-20, 2020b. https://doi.org/10.1080/01140671.2020.1721544
    » https://doi.org/10.1080/01140671.2020.1721544
  • Costa, J. de O; Coelho, R. D.; Barros, T. H. da S.; Fraga Júnior, E. F.; Fernandes, A. L. T. Canopy thermal response to water deficit of coffee plants under drip irrigation. Irrigation and Drainage, v.69, p.472-482, 2020. https://doi.org/10.1002/ird.2429
    » https://doi.org/10.1002/ird.2429
  • Cymes, I.; Dragańska, E.; Brodziński, Z. Potential possibilities of using groundwater for crop irrigation in the context of climate change. Agriculture, v.12, p.1-14, 2022. https://doi.org/10.3390/agriculture12060739
    » https://doi.org/10.3390/agriculture12060739
  • Dejonge, K. C.; Taghvaeian, S.; Thomas, J.; Trout, T. J.; Louise, H.; Comas, L. H. Comparison of canopy temperature-based water stress indices for maize. Agricultural Water Management, v.156, p.51-62, 2015. https://doi.org/10.1016/j.agwat.2015.03.023
    » https://doi.org/10.1016/j.agwat.2015.03.023
  • Doorenbos, J.; Kassam, A. H. Effect of water on crop yield. Rome: Food and Agriculture Organization , 1994. 306p. Irrigation and Drainage Paper, 33
  • Doorenbos, J.; Pruitt, J. O. Crop water requirement. Rome: Food and Agriculture Organization , 1977. 144p. Irrigation and Drainage Paper, 24
  • El-Beltagi, H. S.; Basit, A.; Mohamed, H. I.; Ali, I.; Ullah, S.; Kamel, E. A. R.; Shalaby, T. A.; Ramadan, K. M. A.; Alkhateeb, A. A.; Ghazzawy, H. S. Mulching as a sustainable water and soil saving practice in agriculture: a review. Agronomy, v.12, p.1-31, 2022. https://doi.org/10.3390/agronomy12081881
    » https://doi.org/10.3390/agronomy12081881
  • Jacobi, P. R.; Grandisoli, E. Água e sustentabilidade: desafios, perspectivas e soluções. 1.ed. São Paulo: Reconectta, 2017. 110p
  • García-Tejero, I. F.; Costa, J. M.; Egipto, R.; Durán-Zuazo, V. H.; Lima, R. S. N.; Lopes, C. M.; Chaves, M. M. Thermal data to monitor crop-water status in irrigated Mediterranean viticulture. Agricultural Water Management , v.176, p.80-90, 2016. https://doi.org/10.1016/j.agwat.2016.05.008
    » https://doi.org/10.1016/j.agwat.2016.05.008
  • García-Tejero, I. F.; Rubio, A. E.; Viñuela, I.; Hernández, A.; Gutiérrez-Gordillo, S.; Rodríguez-Pleguezuelo, C. R.; Durán-Zuazo, V. H. Thermal imaging at plant level to assess the crop-water status in almond trees (cv. Guara) under deficit irrigation strategies. Agricultural Water Management , v.208, p.176-186, 2018. https://doi.org/10.1016/j.agwat.2018.06.002
    » https://doi.org/10.1016/j.agwat.2018.06.002
  • Gomes, F. F. L.; Lima, J. S.; Pinheiro Neto, L. G. Aplicações da termografia na agricultura. Agrarian Academy, v.8, p.91-105, 2021. https://doi.org/10.18677/Agrarian_Academy_2021B9
    » https://doi.org/10.18677/Agrarian_Academy_2021B9
  • INMET - Instituto Nacional de Meteorologia. 2018. Available on: <Available on: http://www.inmet.gov.br >. Accessed on: July 2018.
    » http://www.inmet.gov.br
  • Jackson, R. D.; Idso, S. B.; Reginato, R. J.; Pinter, P. J. Canopy temperature as a crop water stress indicator. Water Resources Research, v.17, p.1133-1138, 1981. http://dx.doi.org/10.1029/WR017i004p01133
    » http://dx.doi.org/10.1029/WR017i004p01133
  • Jahan, M.; Nassiri-Mahallati, M. Modeling the response of sesame (Sesamum indicum L.) growth and development to climate change under deficit irrigation in a semi-arid region. PLOS Climate, v.1, p.1-20, 2022. https://doi.org/10.1371/journal.pclm.0000003
    » https://doi.org/10.1371/journal.pclm.0000003
  • Jarrell, W. M.; Armstrong, D. E.; Grigal, D. F.; Kelley, E. F.; Monger, H. C.; Wedin, D. A. Soil water and temperature status. In: Robertson, G. P.; Coleman, D. C.; Bledsoe, C. S.; Sollins, P. (eds.). Standard soil methods for long-term ecological research. New York: Oxford University Press, 1999. 18p.
  • Javadian, M.; Smith, W. K.; Lee, K.; Knowles, J. F.; Scott, R. L.; Fisher, J. B.; Moore, D. J. P.; van Leeuwen, W. J. D.; Barron-Gafford, G.; Behrangi, A. Canopy temperature is regulated by ecosystem structural traits and captures the ecohydrologic dynamics of a semiarid mixed conifer forest site. Journal of Geophysical Research: Biogeosciences, v.127, p.1-15, 2022. https://doi.org/10.1029/2021JG006617
    » https://doi.org/10.1029/2021JG006617
  • Li, S. X.; Wang, Z. H.; Li, S. Q.; Gao, Y. J.; Tian, X. H. Effect of plastic sheet mulch, wheat straw mulch, and maize growth on water loss by evaporation in dryland areas of China. Agricultural Water Management , v.116, p.39-49, 2013. https://doi.org/10.1016/j.agwat.2012.10.004
    » https://doi.org/10.1016/j.agwat.2012.10.004
  • Mira-García, A. B.; Conejero, W.; Vera, J.; Ruiz-Sánchez, M. C. Water status and thermal response of lime trees to irrigation and shade screen, Agricultural Water Management , v.272, p.1-14, 2022. https://doi.org/10.1016/j.agwat.2022.107843
    » https://doi.org/10.1016/j.agwat.2022.107843
  • Ni, J.; Cheng, Y.; Wang, Q.; Ng, C. W. W.; Garg, A. Effects of vegetation on soil temperature and water content: field monitoring and numerical modeling. Journal of Hydrology, v.571, p.494-502, 2019. https://doi.org/10.1016/j.jhydrol.2019.02.009
    » https://doi.org/10.1016/j.jhydrol.2019.02.009
  • Poblete-Echeverría, C.; Sepúlveda-Reyes, D.; Ortega-Farias, S.; Zúñiga, M.; Fuentes, S. Plant water stress detection based on aerial and terrestrial inferred thermography: a study case from vineyard and olive orchard. Acta Horticulturae, n.1112, p.141-146, 2016. https://doi.org/10.17660/ActaHortic.2016.1112.20
    » https://doi.org/10.17660/ActaHortic.2016.1112.20
  • Pou, A.; Diago, M. P.; Medrano, H.; Baluja, J.; Tardaguila, J. Validation of thermal indices for water status identification in grapevine. Agricultural Water Management , v.134, p.60-72, 2014. https://doi.org/10.1016/j.agwat.2013.11.010
    » https://doi.org/10.1016/j.agwat.2013.11.010
  • Ramírez-Cuesta, J. M.; Ortuño, M. F.; Gonzalez-Dugo, V.; Zarco-Tejada, P. J.; Parra, M.; Rubio-Asensio, J. S.; Intrigliolo, D. S. Assessment of peach trees water status and leaf gas exchange using on-the-ground versus airborne-based thermal imagery. Agricultural Water Management , v.267, p.1-14, 2022. https://doi.org/10.1016/j.agwat.2022.107628
    » https://doi.org/10.1016/j.agwat.2022.107628
  • Scholander, P. F.; Hammel, H. T.; Bradstreet, E. D.; Hemmingsen, E. A. Sap pressure in vascular plants. Science, v.148, p.339-346, 1965. https://www.science.org/doi/10.1126/science.148.3668.339
    » https://www.science.org/doi/10.1126/science.148.3668.339
  • Silva, P. S. L.; Silva, P. I. B.; Oliveira, V. R.; Silva, K. E. F.; Monteira, A. L. Doses of carnauba leaf residues and corn yield. Horticultura Brasileira , v.32, p.74-79, 2014. https://doi.org/10.1590/S0102-05362014000100012
    » https://doi.org/10.1590/S0102-05362014000100012
  • Silveira, J. M. de C.; Cândido, B. M.; Silva, A. L. B. de O.; Pavão, G. C.; Xavier, M. A.; Pires, R. C. de M. Uso de imagens multiespectrais e termográficas para monitoramento das condições hídricas da cana-de-açucar. Irriga, v.25, p.689-696, 2020. https://doi.org/10.15809/irriga.2020v25n4p689-696
    » https://doi.org/10.15809/irriga.2020v25n4p689-696
  • United States. Soil Research Team. Keys to soil taxonomy. 12.ed. USDA NRCS, 2014. Available on: <Available on: http://www.nrcs.usda.gov/wps/portal/nrcs/main/soils/survey/ >. Accessed on: Nov. 2022.
    » http://www.nrcs.usda.gov/wps/portal/nrcs/main/soils/survey/
  • Urban, J.; Ingwers, M. W.; Mcguire, M. A.; Teskey, R. O. Increase in leaf temperature opens stomata and decouples net photosynthesis from stomatal conductance in Pinus taeda and Populus deltoides x nigra Journal of Experimental Botany, v.68, p.1757-1767, 2017. https://doi.org/10.1093/jxb/erx052
    » https://doi.org/10.1093/jxb/erx052
  • Wang, Y.; Li, X.; Zhu, J.; Fan, C.; Kong, X.; Turner, N. C.; Siddique, K. H. M.; Li, F. Multi-site assessment of the effects of plastic-film mulch on dry land maize productivity in semiarid areas in China. Agricultural and Forest Meteorology, v.220, p.160-169, 2016. https://doi.org/10.1016/j.agrformet.2016.01.142
    » https://doi.org/10.1016/j.agrformet.2016.01.142
  • 1 Research developed at Instituto Federal de Educação, Ciência e Tecnologia do Ceará/Campus Sobral, Eixo Tecnológico de Recursos Naturais, Sobral, CE, Brazil

Edited by

Editors: Ítalo Herbet Lucena Cavalcante & Hans Raj Gheyi

Publication Dates

  • Publication in this collection
    06 Mar 2023
  • Date of issue
    June 2023

History

  • Received
    21 July 2022
  • Accepted
    27 Dec 2022
  • Published
    11 Feb 2023
Unidade Acadêmica de Engenharia Agrícola Unidade Acadêmica de Engenharia Agrícola, UFCG, Av. Aprígio Veloso 882, Bodocongó, Bloco CM, 1º andar, CEP 58429-140, Campina Grande, PB, Brasil, Tel. +55 83 2101 1056 - Campina Grande - PB - Brazil
E-mail: revistagriambi@gmail.com