Acessibilidade / Reportar erro

LIGHT RAIN CHARACTERIZATION IN PIRACICABA, SÃO PAULO STATE, BRAZIL

ABSTRACT

The disorderly soil occupation without the necessary conservationist practices leads to impacts on the local hydrology and induces the pollution of water resources. This pollution may come from more urbanized areas due to the amount of pollutants drained during the rains. Even moderate precipitations constitute one of the main factors that define pollutant runoff on the surface. These rains have recently been called light rains. Light rains have a lower precipitation height and a higher frequency compared to classic rains of drainage projects, being necessary to define them according to patterns of rain frequency for each region. This study aimed to characterize light rain in the municipality of Piracicaba to establish statistical standards for the frequency of certain precipitation heights. A database provided by the ESALQ/USP automatic weather station, which provides precipitation measurements every 15 minutes, was used in the present study. Light rain heights reached 40.3, 41.4, and 42.7 mm for 100, 90, or 80% frequencies, respectively, which implies the use of return periods of 1.00, 1.11, and 1.25 years, respectively.

diffuse pollution; precipitation; return time; Python

INTRODUCTION

Population growth and increased demand for consumer goods have usually been associated with environmental deterioration. The increase in population and the development of urban areas have led to an increase in impermeable soils and excessive waste generation. This increase is related to the degradation of water resources, a consequence of surface runoff of rainwater, leading the materials disposed on urban surfaces to watercourses ( Yang et al., 2016Yang L, Smith JA, Baeck ML, Zhang Y (2016) Flash flooding in small urban watersheds: Storm event hydrologic response. Water Resources Research 52(6):4571-4589. DOI: https://doi.org/10.1002/2015WRO18326
https://doi.org/10.1002/2015WRO18326...
; Zanchin et al., 2021Zanchin M, Moura MM, Nunes MCM, Tuchtenhahaggen IK, Lima CLR (2021) Assesment of soil loss susceptibility in Santa Rita watershed in Shouthern Brazil. Engenharia Agrícola 41(4)485-495. DOI: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v41n4p485-495/2021
http://dx.doi.org/10.1590/1809-4430-Eng....
).

Polluting sources stand out among the main problems associated with this human intervention and can be specific or diffuse. According to Oliveira et al. (2018)Oliveira ARM, Borges AC, Matos AT, Nascimento M (2018) Estimation on the concentration of suspended solids from turbidity in the water of two sub-basins in the Doce River basin. Engenharia Agricola 38(5):751-759. Doi: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v38n5p751-759/2018
http://dx.doi.org/10.1590/1809-4430-Eng....
and Tansar et al. (2022)Tansar H, Duan HF, Mark O (2022) Catchment-scale and local scale based evaluation of LID effectiveness on urban drainage system performance. Water Resources Management 36(2):507-526. DOI: https://doi.or/10.1007/s11269-021-03036-6
https://doi.or/10.1007/s11269-021-03036-...
, pollutants can be released into the environment in a specific (concentrated in a certain location) or a diffuse way (distributed), with specific discharges consisting of those released in specific places of water bodies.

Diffuse pollution is defined as the indirect introduction of substances into water or soil as a result of human activity, being dispersed and variable in space and not attributed to a specific origin ( Oliveira et al., 2018Oliveira ARM, Borges AC, Matos AT, Nascimento M (2018) Estimation on the concentration of suspended solids from turbidity in the water of two sub-basins in the Doce River basin. Engenharia Agricola 38(5):751-759. Doi: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v38n5p751-759/2018
http://dx.doi.org/10.1590/1809-4430-Eng....
). Leite et al. (2018)Leite MHS, Couto EG, Amorim RSS, Scaramuzza JF (2018) Loss of water and nutrients in diferent soil tillage systems subjected to natural rainfall in the state of Mato Grosso do Sul, Brazil. Engenharia Agrícola 38(6):864-873. DOI: http://dx.doi.org/10.1590/1809-4430.v38n6p864-873/2018
http://dx.doi.org/10.1590/1809-4430.v38n...
state that diffuse pollution occurs over a vast area and is generally associated with the type of land use and occupation, with surface runoff being an important transport agent.

Recent studies have shown that uncontrolled diffuse pollution begins to gain important space in the literature as defined forms of pollution have been controlled in several countries ( Cembranel et al., 2017Cembranel AS, Elisandro PF, Sampaio SC, Mercante E, Reis R, Remor MB (2017) Residue analysis of organoclorine and organophosphorus pesticides in urban lake sediment. Engenharia Agrícola 37(6):1254-1267. DOI: http://dx.doi.org/10.1590/1809-4430.v37n6p1254-1267/2017
http://dx.doi.org/10.1590/1809-4430.v37n...
). Diffuse forms are more difficult to resolve and often cause high levels of pollution in waterways, which are related to their dispersion in space ( Batista et al., 2019Batista ML, Coelho G, Mello CR, Oliveira MS (2019) Spatialization of the annual maximum daily rainfall in southeastern Brazil. Engenharia Agrícola 39(1):97-109. DOI: http://dx.doi.org/10.1590/1809-4430.v39n1p97-109/2019
http://dx.doi.org/10.1590/1809-4430.v39n...
). According to Leite et al. (2018)Leite MHS, Couto EG, Amorim RSS, Scaramuzza JF (2018) Loss of water and nutrients in diferent soil tillage systems subjected to natural rainfall in the state of Mato Grosso do Sul, Brazil. Engenharia Agrícola 38(6):864-873. DOI: http://dx.doi.org/10.1590/1809-4430.v38n6p864-873/2018
http://dx.doi.org/10.1590/1809-4430.v38n...
, the processes and causes of polluting sources, especially diffuse ones, need to be analyzed for their control and study. Thus, analyzing precipitation is essential for the beginning of its study.

Patios on which windrows of filter cake from the sugar and alcohol plants are arranged have a detention basin downstream, responsible for the temporary storage of residues carried by the rains. Currently, there is no sizing criterion for these types of reservoirs, as the influence of rainfall on water quality is difficult to quantify. Pollutants deposited on the surface are carried to the water bodies as the rains occur, being gradually distributed by runoff, and released into rivers ( Kuo et al., 2017Kuo NW, Jien SH, Hong NM, Chen YT, Lee TY (2017) Contribution of urban runoff in Taipei metropolitan area to dissolved inorganic nitrogen export in the Danshui River, Taiwan. Environmental Science and Pollution Research 24(1):578-590. DOI 10.1007/s11356-016-7825-4 ).

According to Pruski et al. (2015)Pruski FF, Rodriguez RG, Nunes AA, Pruski PL, Singh VP (2015) Low flow stimates in regions of extrapolation of the regionalization equations: a new concept. Engenharia Agrícola 35(5):808-816. DOI: http://dx.doi.org/10.1590/1809-4430.v35n5p808-816/2015
http://dx.doi.org/10.1590/1809-4430.v35n...
and Yoon et al. (2010)Yoon SW, Chung SW, Oh DG, Lee JW (2010) Monitoring of non-point source pollutants load from a mixed forest land use. Journal of Environment Science 22(6):801–805. DOI: https://doi.org/10.1016/51001-0742(09)601-7
https://doi.org/10.1016/51001-0742(09)60...
, pollutants accumulate in urban basins during dry weather, being produced not only by human activity but also by atmospheric deposition. The first rapid contribution when rain occurs after a period of drought consists of removing pollutants from impermeable surfaces, providing them highly concentrated in wastewater ( Ahmed et al., 2022Ahmed N, Wang G, Booij MJ, Xiangyang S, Hussain F, Nabi G (2022) Separation of the impact of landuse/landcover change and climatic change on runoff in the upstream area of the Yangtze River, China. Water Resources Management 36:181-201. DOI: https://doi.org/10.1007/s11269-021-030021-z
https://doi.org/10.1007/s11269-021-03002...
).

Precipitation causes an initial wetting on the surface and this layer is infiltrated ( Oliveira et al., 2018Oliveira ARM, Borges AC, Matos AT, Nascimento M (2018) Estimation on the concentration of suspended solids from turbidity in the water of two sub-basins in the Doce River basin. Engenharia Agricola 38(5):751-759. Doi: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v38n5p751-759/2018
http://dx.doi.org/10.1590/1809-4430-Eng....
; Zanchin et al., 2021Zanchin M, Moura MM, Nunes MCM, Tuchtenhahaggen IK, Lima CLR (2021) Assesment of soil loss susceptibility in Santa Rita watershed in Shouthern Brazil. Engenharia Agrícola 41(4)485-495. DOI: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v41n4p485-495/2021
http://dx.doi.org/10.1590/1809-4430-Eng....
). The dissolution of deposited soluble elements begins in this phase and runoff begins after passing the wetting phase. These authors state that runoff is associated with the dragging and transport of particles disposed on patios, pointing to associated factors, such as the kinetic energy of raindrops, the runoff speed, and the cohesion of deposits.

This concentration of pollutants is carried by the so-called first flush, which is defined as the first volume of precipitation that is drained. The first flush phenomena are recognized as typical of highly populated and urbanized regions ( Damé et al., 2016Damé RCF, Teixeira-Gandra CFA, Guedes HAS, Silva GM, Silveira SCR (2016) Intensity-duration-frequency relationships: stochastic modeling and disaggregation of daily rainfall in the Lagoa Mirim Watershed, Rio Grande do Sul, Brazil. Engenharia Agrícola 36(3):492-502. DOI: http://dx.doi.org/10.1590/1809-4430.v36n3p492-502/2016
http://dx.doi.org/10.1590/1809-4430.v36n...
). The first flush phenomenon can be described as the runoff of rainwater during which pollutant concentration is significantly higher than in later precipitation stages ( Beck et al., 2017Beck HE, Vergopolan N, Pan M, Levizzani V, van Dijk AIJM, Weedon GP, Brocca L, Pappenberger F, Huffman GJ, Wood EF (2017) Global-scale evaluation of 22 precipitation datasets using gauge observations and hydrological modeling. Hydrology and Earth System Sciences 21(12):6201-6217. DOI: https://dx.doi.org/10.5194/hess-21-6201-2017
https://dx.doi.org/10.5194/hess-21-6201-...
). As pointed out by Pitt (1999)Pitt RE (1999) Small storm hydrology and why it is important for the design of stormwater control practices. Journal of Water Management Modeling Chapter 4, 32 p. DOI: http://dx.doi.org/10.14796/JWMM.R204-04
http://dx.doi.org/10.14796/JWMM.R204-04...
, frequent rainfall events are responsible for most of the runoff of massive discharges of pollutants. Therefore, these events should be used as the basis for water quality control.

Usually, the diffuse pollution that is associated with rain is more important the longer the period of drought that precedes the referred phenomenon, that is, extended drought periods. The accumulation of materials and pollutants in urban regions increases during longer periods of lack of precipitation, as pollutant elimination is controlled by physical processes such as rain associated with the interval of occurrence, intensity, and duration ( Kruel et al., 2015Kruel IB, Meschiatti MC, Blain GC, Ávila AMH (2015) Climatics trends in the municipality of Pelotas, state of Rio Grande do Sul, Brazil. Engenharia Agrícola 35(4):769-777. DOI: http://dx.doi.org/10.1590/1809-4430.v35n4p769-777/20185
http://dx.doi.org/10.1590/1809-4430.v35n...
). The need generated by analyses of runoff frequency related to diffuse pollution has led to the concept of light rain hydrology, also called small storms ( Pitt, 1999Pitt RE (1999) Small storm hydrology and why it is important for the design of stormwater control practices. Journal of Water Management Modeling Chapter 4, 32 p. DOI: http://dx.doi.org/10.14796/JWMM.R204-04
http://dx.doi.org/10.14796/JWMM.R204-04...
).

According to Oliveira et al. (2018)Oliveira ARM, Borges AC, Matos AT, Nascimento M (2018) Estimation on the concentration of suspended solids from turbidity in the water of two sub-basins in the Doce River basin. Engenharia Agricola 38(5):751-759. Doi: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v38n5p751-759/2018
http://dx.doi.org/10.1590/1809-4430-Eng....
, detention basins stand out among the optimal measures for managing diffuse pollution. Pruski et al. (2015)Pruski FF, Rodriguez RG, Nunes AA, Pruski PL, Singh VP (2015) Low flow stimates in regions of extrapolation of the regionalization equations: a new concept. Engenharia Agrícola 35(5):808-816. DOI: http://dx.doi.org/10.1590/1809-4430.v35n5p808-816/2015
http://dx.doi.org/10.1590/1809-4430.v35n...
recommended the adoption of minimum intervals of duration of 3, 6, 12, or 24 hours in the hydrological modeling to size these basins.

The so-called light rain hydrology method was developed to estimate the volume of runoff for use in urban and suburban soils for relatively small rainfall events and plan the control and quality of rainwater ( Pitt & Clark, 2008Pitt RE, Clark SE (2008) Integrated storm-water management for watershed sustainability. Journal of Irrigation and Drainage Engineering 134(5):548-555. DOI: https://doi.org/10.1061/(ASCE)0733-9437/(2008)134:5(548)
https://doi.org/10.1061/(ASCE)0733-9437/...
).

Different drainage project criteria generally require the examination of different types of rainfall ( Pitt, 1999Pitt RE (1999) Small storm hydrology and why it is important for the design of stormwater control practices. Journal of Water Management Modeling Chapter 4, 32 p. DOI: http://dx.doi.org/10.14796/JWMM.R204-04
http://dx.doi.org/10.14796/JWMM.R204-04...
; Pitt & Clark, 2008Pitt RE, Clark SE (2008) Integrated storm-water management for watershed sustainability. Journal of Irrigation and Drainage Engineering 134(5):548-555. DOI: https://doi.org/10.1061/(ASCE)0733-9437/(2008)134:5(548)
https://doi.org/10.1061/(ASCE)0733-9437/...
). Similarly, water quality issues due to rainwater runoff are typically associated with minor storms rather than the project storms commonly used for surface drainage projects ( Tansar et al., 2022)Tansar H, Duan HF, Mark O (2022) Catchment-scale and local scale based evaluation of LID effectiveness on urban drainage system performance. Water Resources Management 36(2):507-526. DOI: https://doi.or/10.1007/s11269-021-03036-6
https://doi.or/10.1007/s11269-021-03036-...
.

The classification of light rain has variations in different regions, that is, its concept may vary according to certain parameters to be studied and analyzed, with urbanized areas being mainly those with higher concentrations of pollutants. According to Thiemig (2012), rainwater volume is often dimensioned through the concept of light rain, which is defined as a precipitation height greater than 2 mm capable of generating sufficient surface runoff to carry most of the diffuse pollutant load on the surface.

According to Pedron & Klosowski (2008)Pedron IT, Klosowski ES (2008) Distribuição de frequência de chuvas diárias no Estado do Paraná. Scientia Agraria Paranaensis 7(1-2):55-63. DOI: https://doi.org/10.18188/sap.v0i0.2052
https://doi.org/10.18188/sap.v0i0.2052...
, light rains are defined as rains higher than or equal to 2 millimeters that have a 90% probability of occurring with higher or equal height. They are recommended among the best management practices (BMP) for managing diffuse pollution.

In contrast, rains between 13 and 38 millimeters are responsible for about 76% of pollutant discharges in the runoff, being important to be addressed in the study of massive loads of pollutants ( Bannerman et al., 1983Bannerman R, Baun K, Bohn N, Hughes PE, Graczik DA (1983) Evaluation of Urban Nonpoint Source Pollution Management in Milawllukee County. University of Visconsin, Madison. Department of Natural Resources of United States. Vol. L PB 84-114164. US Environmental Protection Agency, Water Planning Division, November 1983. ). Rainfall higher than 38 mm is associated with drainage and accounts for only relatively small portions of annual pollutant discharges. Moreover, pollution control depends on local precipitation, and events with an 80% frequency are recommended to be considered for analysis ( UDFCD, 2011URBAN DRAINAGE AND FLOOD CONTROL DISTRICT (UDFCD) (2011) Urban storm drainage criteria manual: Volume 3. Stormwater best management practices. Denver, 565p. Available: www.udfcd.org
www.udfcd.org...
).

Therefore, quality control criteria for urban runoff are part of many guidelines for rainwater management. Capturing and treating intense extreme events are recognized not to be necessary, and the most commonly defined target in these guidelines is based on a more frequent distribution precipitation analysis aiming to determine the rainfall volume to be used for the project of reservoirs and detention basins ( Rivard, 2010Rivard G (2010) Small storm hydrology and BMP modeling with SWMM5. Journal of Water Management Modeling Chapter 10. 18 p. DOI: http://dx.doi.org/10.14796/JWMM.R236-10
http://dx.doi.org/10.14796/JWMM.R236-10...
).

Rainfall parameters that fit the literature recommendations are necessary when analyzing the frequencies of precipitation heights, stipulating them according to their return time (Tr). The return time (inverse of the frequency) is necessary to decide the degree of protection given to the population or, in other words, the risk to which the population is subject ( Passos et al., 2021Passos JBMC, Silva DD, Lima RPC (2021) Daily rainfall disaggregation coeficients for the Doce River basin, Brasil. Engenharia Agrícola 41(2): 223-234. DOI: http://dx.doi.org/10.1590/1809-4430.v41n2p223-234/2021
http://dx.doi.org/10.1590/1809-4430.v41n...
; Manke et al., 2022Manke EB, Teixeira-Gandra CFA, Damé RCE, Nunes AB, Chagas Neta MCC, Karsburg RM (2022) Seasonal intensity-duration-frequency relationships for Pelotas, Rio Grande do Sul, Brazil. Revista Brasileira de Engenharia Agrícola e Ambiental 26(2):85-90. DOI: http://dx.doi.org/10.1590/1807-1929/Agriambi.v26n2p85-90
http://dx.doi.org/10.1590/1807-1929/Agri...
).

According to Yang et al. (2016)Yang L, Smith JA, Baeck ML, Zhang Y (2016) Flash flooding in small urban watersheds: Storm event hydrologic response. Water Resources Research 52(6):4571-4589. DOI: https://doi.org/10.1002/2015WRO18326
https://doi.org/10.1002/2015WRO18326...
, pollutants are quickly drained into the water body for rain events with return periods of 2, 3, and 5 years, that is, consist of excessive rainfall. However, the situation is reversed for a return time of about 1 year, that is, the concentration of pollution remains longer in the water body.

This study aimed to characterize the light precipitation in the municipality of Piracicaba based on the process of filtering and separating rainfall data, seeking to establish recommended statistical standards for certain lower precipitation heights, which occur with higher frequency. The focus was given to less intense rains or light rains, which have been frequently used to simulate the displacement of diffuse pollution in urban or agro-industrial patios, such as the compost patios of sugar and alcohol plants in São Paulo.

SUBJECT DESCRIPTION

Area of interest for carrying out the study

This study was carried out in the municipality of Piracicaba, São Paulo State, Brazil, located at a critical point of transition of general circulation currents, with varied climate characteristics. This variability causes two distinct periods: a hotter and rainy period from November to March and a less hot and drier period from May to August ( Ferreira et al., 2021Ferreira FLV, Rodrigues LN, Silva DD, Teixeira DBS, Almeida LT (2021) Time series trends of stream flow and rainfall in Santo Antônio River basin, Brazil. Engenharia Agrícola 41(1):47-55. DOI: http://dx.doi.org/10.1590/1809-4430.v41n1p47-55/2021
http://dx.doi.org/10.1590/1809-4430.v41n...
).

Wolff et al. (2014)Wolff W, Duarte SN, Mingoti R (2014) Nova metodologia de regionalização de vazões: estudo de caso para o estado de São Paulo. Revista Brasileira de Recursos Hídricos 19(4)21-33. DOI: http://dx.doi.org/10.21168/rbrh.v19n4.p21-33
http://dx.doi.org/10.21168/rbrh.v19n4.p2...
state that the region has tropical trends, that is, mean annual temperatures vary around 20 °C,

reaching a mean of 23 °C in the hottest months of the year and 17 °C in the coldest months. The mean long-term annual precipitation in the region is around 1,280 mm ( Alvares et al., 2013Alvares CA, Stape JL, Sentelhas PC, Gonçalves JLM, Sparovek G (2013) Köppen's climate classification map for Brazil. Meteorologische Zeitschrift 22(6):711- 728. DOI: http://dx.doi.org/10.1127/0941-2948/2013/0507
http://dx.doi.org/10.1127/0941-2948/2013...
).

The municipality of Piracicaba is located at the geographic coordinates 22°42′30″ S latitude and 47°38′01″ W longitude, with a mean altitude of 554 m. The total area has 1,378 km2 and the population is estimated at 410,275 inhabitants (Lucas et al., 2018). Figure 1 represents the location of the municipality of Piracicaba in the São Paulo State, Brazil.

FIGURE 1
Location of the municipality of Piracicaba in the São Paulo State, Brazil. Source: IBGE (2021).

Piracicaba is currently a new metropolitan region of the São Paulo State, with around 1.5 million inhabitants and an area of 7,367.88 km2.

Precipitation data

The first part of the study covered the search and selection of data. The data came from the automatic weather station of the Luiz de Queiroz College of Agriculture (ESALQ/ USP), located in Piracicaba-SP, which presents precise data measured every fifteen minutes. This small measurement time was essential for the detail and accuracy of the analyzed data (0.1 mm precision). Figure 2 shows the location of the weather station in the municipality.

FIGURE 2
Location of the weather station in the municipality of Piracicaba, São Paulo State, Brazil.

Following the first activity, data every 15 minutes between 1997 and 2020 were collected, totaling 24 years. The data from 2021 were disregarded in the analysis due to a lack of complete data.

Data analysis

A verification was carried out after collecting the data to find flaws in the chronology or gaps, which can occur mainly in precise and frequent data such as those used in this study. This second stage of the methodology was carried out entirely in Excel® based on an investigation of failures and adjustments. This step was extremely important and fundamental for the subsequent steps.

A spreadsheet was created after analyzing the data and adjusting with essential data such as date, time, and rainfall, formatted in Excel®. All spreadsheets were standardized over a period of 24 years to facilitate interpretation by the used programming language (Python Software Foundation). The Python language was chosen due to the ease of operating a large number of data arranged in a spreadsheet with few command lines. Table 1 shows, as an example, part of the spreadsheet updated with the new database for part of the first day of January 2021, already filtered and filled in.

TABLE 1
Part of the new database systematized from data from the automatic weather station exemplified for the beginning of January 1, 2020.

The entire previous process was defined and thought out for the subsequent stage, which was based on working with the extracted database, checked and adjusted entirely by programming in Python. Python language is used widely by analysts and data scientists, and its programming provides some databases that favor the manipulation of information and can be used to check variables in a certain time interval ( Mello & Silva, 2009Mello CR, Silva AM (2009) Modelagem estatística da precipitação mensal e anual e no período seco para o estado de Minas Gerais. Revista Brasileira de Engenharia Agrícola e Ambiental 13(1): 68-74. DOI: http://dx.doi.org/10.1590/S1415-43662009000100010
http://dx.doi.org/10.1590/S1415-43662009...
).

Python version 3.8, provided by Anaconda (2021), was used. The Anaconda platform is one of the most popular worldwide, as it provides a variety of data already pre-installed. Specific software for reading the computational language, the free software PyCharm Community Edition, specific for Python, was used to facilitate access to the programming language and its development.

The process of manipulating data via computational language was initiated aiming at counting and verifying them. Therefore, one event was separated from the other as long as there were 12 hours or more without precipitation relative to the next event. Figure 3 shows the flowchart of the main part of the programming.

FIGURE 3
Main flowchart part of the software developed in Python aiming to obtain light rain in Piracicaba, São Paulo State, Brazil.

The criteria for separating events with a period of 12 hours or more without rain is adopted for light rain events that have a higher frequency than extreme events. Therefore, we considered that there would be maintenance during this interval of at least 12 hours, with cleaning being carried out in the detention basins and exit devices to reduce polluting loads. In addition, as most rainfall records reported in the literature have a minimum duration of 24 hours without rain, considering a period of 12 hours is a practical way to facilitate future studies that use longer intervals.

This criterion of using at least 12 hours also facilitates the design of diffuse pollution control structures, as suggested by Louzada et al. (2018)Louzada FLRO, Xavier AC, Pessopane JEM (2018) Climatological water balance with data estimated by Tropical Rainfall Measuring Mission for Doce River basin. Engenharia Agrícola 38(3):376-386. DOI: http://dx.doi.org/10.1590/1809-4430.v38n3p376-386/2018
http://dx.doi.org/10.1590/1809-4430.v38n...
. In this interval logic, these drainage systems must have a final reservoir for rainwater accumulation or a detention basin, which is used as a temporary reservoir for a period lower than 24 hours.

Filtering was developed directly in Python once the analysis time was defined by counting the period of 12 hours or more without rain before restarting the count, as the amount of pollutants is concentrated in longer periods of drought ( Leite et al., 2018Leite MHS, Couto EG, Amorim RSS, Scaramuzza JF (2018) Loss of water and nutrients in diferent soil tillage systems subjected to natural rainfall in the state of Mato Grosso do Sul, Brazil. Engenharia Agrícola 38(6):864-873. DOI: http://dx.doi.org/10.1590/1809-4430.v38n6p864-873/2018
http://dx.doi.org/10.1590/1809-4430.v38n...
). Therefore, mixing drought factors and the fact that the maintenance of the detention reservoir can be carried out within a period of 12 hours, the number of events was counted.

The criterion of counting every 48 lines or more without precipitation was used for the computational language, as the analyzed data had information every 15 minutes, being equivalent to the period of 12 hours or more, thus counting one event.

After this count, the results were converted into an Excel® file directly from the Python directory. The sum and counting of events higher than 2, 5, 10, 15 mm and so on up to 50 mm was performed in Excel®.

An exploratory analysis was performed to start manipulating the results by calculating the mean, median, minimum, maximum, and standard deviation of the number of events counted.

The last analysis consisted of the graphical representation by plotting the precipitation as a function of the stipulated return times, covering a series of 24 years of data. For this purpose, the data were placed in descending order and the empirical frequencies of occurrence with higher or equal heights were calculated ( Passos et al., 2021Passos JBMC, Silva DD, Lima RPC (2021) Daily rainfall disaggregation coeficients for the Doce River basin, Brasil. Engenharia Agrícola 41(2): 223-234. DOI: http://dx.doi.org/10.1590/1809-4430.v41n2p223-234/2021
http://dx.doi.org/10.1590/1809-4430.v41n...
). The data were separated with this formulation, and a spreadsheet with the number of events in the interval, the total value, and the empirical probability of each rain height being equaled or exceeded and its respective return time, that is, the inverse of the frequency, was created in Excel®.

Finally, considering that several rainfall heights with a return period of less than 1 year were generated, we decided to perform an adjustment to a simple logarithmic function instead of using a classic probability distribution function.

RESULTS

Statistical information for the rainfall data from the automatic weather station referring to the period from 1997 to 2020, with measurements of rainfall heights in the range of 2 to 5 mm, 5 to 10 mm, 10 to 15 mm, and so on up to higher than 50 mm, were the initial basis of this study.

The first analyses measured the descriptive statistical data of the number of events found per year in the intervals of the respective rainfall heights, as shown in Table 2 .

TABLE 2
Descriptive statistical data on the number of rainfall events recorded per year and the total period (1997–2020) for the Piracicaba weather station.

A total of 2,568 events larger than 2 mm were recorded in the period of 24 years, which means an average of 107.0 events per year. Rain events with a lower rain height presented a higher number of occurrences per year and in the period of 24 years analyzed. On the contrary, events with higher rain heights occurred with significantly lower frequency. According to Rivard (2010)Rivard G (2010) Small storm hydrology and BMP modeling with SWMM5. Journal of Water Management Modeling Chapter 10. 18 p. DOI: http://dx.doi.org/10.14796/JWMM.R236-10
http://dx.doi.org/10.14796/JWMM.R236-10...
, selecting the number of total rainfall events and the average amount of rainfall over the years can also be a way to obtain a deeper analysis of the different categories regarding the generated precipitation volumes. Moreover, rains that have a lower limit of 25 mm have been adopted in many guidelines related to drainage and management of rainwater ( Rivard, 2010Rivard G (2010) Small storm hydrology and BMP modeling with SWMM5. Journal of Water Management Modeling Chapter 10. 18 p. DOI: http://dx.doi.org/10.14796/JWMM.R236-10
http://dx.doi.org/10.14796/JWMM.R236-10...
).

Rivard (2010)Rivard G (2010) Small storm hydrology and BMP modeling with SWMM5. Journal of Water Management Modeling Chapter 10. 18 p. DOI: http://dx.doi.org/10.14796/JWMM.R236-10
http://dx.doi.org/10.14796/JWMM.R236-10...
conducted studies in the city of Milwaukee (United States) and observed that most events have a lower rain height, classifying them as: (i) very frequent rainfall events (< 10 mm); (ii) common rains (< 14 mm and < 22 mm); and (iii) intense events (≥ 32 mm), corresponding to less than 5% of precipitation events. Generally, very small events can be considered infiltrated and indicated for groundwater recharge ( Pangle et al., 2022Pangle LA, Diem JE, Milligan R, Adams E, Murray A (2022) Contextualizing inflow and infiltration within the streamflow regime of urban watersheds. Water Resources Research 58(1):1-19. DOI: https://doi.org/10.1029/2021WR030406
https://doi.org/10.1029/2021WR030406...
). In contrast, common events are the most used for designing water quality control structures ( Conceição et al., 2018Conceição KZ, Villas Boas, MA, Sampaio, SC, Remor, MB, Bonaparte, DI (2018) Statistical control of the process applied to the monitoring of the water quality index. Engenharia Agrícola 38(6):951-960. DOI: http://dx.doi.org/10.1590/1809-4430.v38n6p951-960/2018
http://dx.doi.org/10.1590/1809-4430.v38n...
; Aguiar et al., 2019Aguiar W, Sampaio S, Paisani JC, Remor MB, Reis RR (2019) Seasonal dynamics of agricultural soil cover in runoff generation. Engenharia Agrícola 29(5):592-599. DOI: http://dx.doi.org/10.1590/1809-4430.v29n5p592-599/2019
http://dx.doi.org/10.1590/1809-4430.v29n...
; Lu et al., 2022)Lu Q, Zhong P, Xu B, Huang X, Zhu F, Wang H, Ma Y (2022) Multi-objective risk analysis for flood control operation of a complex reservoir system under multiple time-space correlated uncertanties. Journal of Hydrology 606(3):1-17. DOI: https://doi.org/10.1016/J.JHYDROL.2021.127419
https://doi.org/10.1016/J.JHYDROL.2021.1...
.

Consistency can be observed between the analyses carried out by Rivard (2010)Rivard G (2010) Small storm hydrology and BMP modeling with SWMM5. Journal of Water Management Modeling Chapter 10. 18 p. DOI: http://dx.doi.org/10.14796/JWMM.R236-10
http://dx.doi.org/10.14796/JWMM.R236-10...
and those obtained in this study, which had a percentage of occurrence of 4.7% for rainfall in the range of 30 to 35 mm in the municipality of Piracicaba, for example. The analyses carried out by Rivard (2010)Rivard G (2010) Small storm hydrology and BMP modeling with SWMM5. Journal of Water Management Modeling Chapter 10. 18 p. DOI: http://dx.doi.org/10.14796/JWMM.R236-10
http://dx.doi.org/10.14796/JWMM.R236-10...
reached a value of 5%.

According to UDFCD (2011)URBAN DRAINAGE AND FLOOD CONTROL DISTRICT (UDFCD) (2011) Urban storm drainage criteria manual: Volume 3. Stormwater best management practices. Denver, 565p. Available: www.udfcd.org
www.udfcd.org...
and Pedron & Klosowski (2008)Pedron IT, Klosowski ES (2008) Distribuição de frequência de chuvas diárias no Estado do Paraná. Scientia Agraria Paranaensis 7(1-2):55-63. DOI: https://doi.org/10.18188/sap.v0i0.2052
https://doi.org/10.18188/sap.v0i0.2052...
, the project rainfall of the rainwater collection system to control diffuse pollution can be obtained by correlating the rainfall for dimensioning the reservoirs with frequencies (if higher or equal) of 80 or 90%.

Rains to control diffuse pollution were obtained according to the graph in Figure 4 .

FIGURE 4
Light rains as a function of the return period (Tr) for Piracicaba, São Paulo State, Brazil, with empirical logarithmic adjustment.

Figure 4 shows the ratio of precipitation as a function of Tr. It allows stipulating the project rain for its respective Tr value. Thus, the light rain that best fits the design of future structures ( Figure 4 ) according to Yang et al. (2016)Yang L, Smith JA, Baeck ML, Zhang Y (2016) Flash flooding in small urban watersheds: Storm event hydrologic response. Water Resources Research 52(6):4571-4589. DOI: https://doi.org/10.1002/2015WRO18326
https://doi.org/10.1002/2015WRO18326...
(Tr = 1.0 year), considering the sources of diffuse pollution, should have a height slightly higher than 40 mm (exactly 40.267 mm, according to logarithmic adjustment).

Table 3 shows the performance indices of the logarithmic adjustment. The performance can be considered very good for all indices ( Moriasi et al., 2015Moriasi DN, Gitau M, Pai N, Daggupati P (2015) Hydrologic and water quality models: performance measures and evaluation criteria. American Society of Agricultural and Biological Engineers 58(6):1763–1785. DOI: http://dx.doi.org/10.13031/trans.58.10715
http://dx.doi.org/10.13031/trans.58.1071...
).

TABLE 3
Performance indices of the logarithmic model used in the light rain adjustment.

According to UDFCD (2011)URBAN DRAINAGE AND FLOOD CONTROL DISTRICT (UDFCD) (2011) Urban storm drainage criteria manual: Volume 3. Stormwater best management practices. Denver, 565p. Available: www.udfcd.org
www.udfcd.org...
and Pedron & Klosowski (2008)Pedron IT, Klosowski ES (2008) Distribuição de frequência de chuvas diárias no Estado do Paraná. Scientia Agraria Paranaensis 7(1-2):55-63. DOI: https://doi.org/10.18188/sap.v0i0.2052
https://doi.org/10.18188/sap.v0i0.2052...
, the project light rain can be obtained by correlating the heights with frequencies of occurrence higher than or equal to 80 and 90%, respectively. This criterion implies the use of Tr values of 1.25 and 1.11 years, which provide heights of 42.7 and 41.4 mm, respectively, by the logarithmic equation.

The obtained values are slightly high compared to the suggestions by Bannerman et al. (1983)Bannerman R, Baun K, Bohn N, Hughes PE, Graczik DA (1983) Evaluation of Urban Nonpoint Source Pollution Management in Milawllukee County. University of Visconsin, Madison. Department of Natural Resources of United States. Vol. L PB 84-114164. US Environmental Protection Agency, Water Planning Division, November 1983. and Rivard (2010)Rivard G (2010) Small storm hydrology and BMP modeling with SWMM5. Journal of Water Management Modeling Chapter 10. 18 p. DOI: http://dx.doi.org/10.14796/JWMM.R236-10
http://dx.doi.org/10.14796/JWMM.R236-10...
, who suggest heights from 13 to 38 mm and 25 to 32 mm, respectively. Thus, Tr values a little lower than those adopted may be more realistic for a sub-humid tropical region, such as Piracicaba. It means that laboratory and field experiments must be carried out to validate and adapt the methodology adopted in this study.

CONCLUSIONS

  1. The applied frequency study methodology for the definition of light rain was well suited to the raised historical series, as well as the comparative analysis between the occurrence of events by height interval for each year and the cumulative for the 24 years of duration of the precipitation series.

  2. Light rains fit well to a logarithmic function and occur more frequently compared to more intense rains, which are usually used in surface drainage.

  3. Most consulted authors ( Yang et al., 2016Yang L, Smith JA, Baeck ML, Zhang Y (2016) Flash flooding in small urban watersheds: Storm event hydrologic response. Water Resources Research 52(6):4571-4589. DOI: https://doi.org/10.1002/2015WRO18326
    https://doi.org/10.1002/2015WRO18326...
    ; Pedron & Klosowski, 2008Pedron IT, Klosowski ES (2008) Distribuição de frequência de chuvas diárias no Estado do Paraná. Scientia Agraria Paranaensis 7(1-2):55-63. DOI: https://doi.org/10.18188/sap.v0i0.2052
    https://doi.org/10.18188/sap.v0i0.2052...
    ; and UDFCD, 2011)URBAN DRAINAGE AND FLOOD CONTROL DISTRICT (UDFCD) (2011) Urban storm drainage criteria manual: Volume 3. Stormwater best management practices. Denver, 565p. Available: www.udfcd.org
    www.udfcd.org...
    presented similar suggestions for estimating light rain, namely, the use of return period values of 1.00, 1.11, and 1.25 years, obtained from total series of rainfall data.

  4. Considering the adjustments and the sources of recommendations, these return period values can be used for the Piracicaba region and imply the use of project rainfall of 40.3, 41.4, and 42.7 mm to control diffuse pollution, respectively, if no significant future climate variations occur in the region.

  5. The light rain characterization carried out for the municipality of Piracicaba can be useful for the design of detention basins to control diffuse pollution, such as those existing downstream of composting patios of sugar and alcohol plants. However, there is a need for future experimental evidence of rain heights obtained in the region to validate the proposed methodology.

REFERENCES

  • Aguiar W, Sampaio S, Paisani JC, Remor MB, Reis RR (2019) Seasonal dynamics of agricultural soil cover in runoff generation. Engenharia Agrícola 29(5):592-599. DOI: http://dx.doi.org/10.1590/1809-4430.v29n5p592-599/2019
    » http://dx.doi.org/10.1590/1809-4430.v29n5p592-599/2019
  • Alvares CA, Stape JL, Sentelhas PC, Gonçalves JLM, Sparovek G (2013) Köppen's climate classification map for Brazil. Meteorologische Zeitschrift 22(6):711- 728. DOI: http://dx.doi.org/10.1127/0941-2948/2013/0507
    » http://dx.doi.org/10.1127/0941-2948/2013/0507
  • ANACONDA. Anaconda: Individual Edition. Versão 3.8. DOI: https://www.anaconda.com/products/individual
    » https://www.anaconda.com/products/individual
  • Ahmed N, Wang G, Booij MJ, Xiangyang S, Hussain F, Nabi G (2022) Separation of the impact of landuse/landcover change and climatic change on runoff in the upstream area of the Yangtze River, China. Water Resources Management 36:181-201. DOI: https://doi.org/10.1007/s11269-021-030021-z
    » https://doi.org/10.1007/s11269-021-030021-z
  • Bannerman R, Baun K, Bohn N, Hughes PE, Graczik DA (1983) Evaluation of Urban Nonpoint Source Pollution Management in Milawllukee County. University of Visconsin, Madison. Department of Natural Resources of United States. Vol. L PB 84-114164. US Environmental Protection Agency, Water Planning Division, November 1983.
  • Batista ML, Coelho G, Mello CR, Oliveira MS (2019) Spatialization of the annual maximum daily rainfall in southeastern Brazil. Engenharia Agrícola 39(1):97-109. DOI: http://dx.doi.org/10.1590/1809-4430.v39n1p97-109/2019
    » http://dx.doi.org/10.1590/1809-4430.v39n1p97-109/2019
  • Beck HE, Vergopolan N, Pan M, Levizzani V, van Dijk AIJM, Weedon GP, Brocca L, Pappenberger F, Huffman GJ, Wood EF (2017) Global-scale evaluation of 22 precipitation datasets using gauge observations and hydrological modeling. Hydrology and Earth System Sciences 21(12):6201-6217. DOI: https://dx.doi.org/10.5194/hess-21-6201-2017
    » https://dx.doi.org/10.5194/hess-21-6201-2017
  • Cembranel AS, Elisandro PF, Sampaio SC, Mercante E, Reis R, Remor MB (2017) Residue analysis of organoclorine and organophosphorus pesticides in urban lake sediment. Engenharia Agrícola 37(6):1254-1267. DOI: http://dx.doi.org/10.1590/1809-4430.v37n6p1254-1267/2017
    » http://dx.doi.org/10.1590/1809-4430.v37n6p1254-1267/2017
  • Conceição KZ, Villas Boas, MA, Sampaio, SC, Remor, MB, Bonaparte, DI (2018) Statistical control of the process applied to the monitoring of the water quality index. Engenharia Agrícola 38(6):951-960. DOI: http://dx.doi.org/10.1590/1809-4430.v38n6p951-960/2018
    » http://dx.doi.org/10.1590/1809-4430.v38n6p951-960/2018
  • Damé RCF, Teixeira-Gandra CFA, Guedes HAS, Silva GM, Silveira SCR (2016) Intensity-duration-frequency relationships: stochastic modeling and disaggregation of daily rainfall in the Lagoa Mirim Watershed, Rio Grande do Sul, Brazil. Engenharia Agrícola 36(3):492-502. DOI: http://dx.doi.org/10.1590/1809-4430.v36n3p492-502/2016
    » http://dx.doi.org/10.1590/1809-4430.v36n3p492-502/2016
  • Ferreira FLV, Rodrigues LN, Silva DD, Teixeira DBS, Almeida LT (2021) Time series trends of stream flow and rainfall in Santo Antônio River basin, Brazil. Engenharia Agrícola 41(1):47-55. DOI: http://dx.doi.org/10.1590/1809-4430.v41n1p47-55/2021
    » http://dx.doi.org/10.1590/1809-4430.v41n1p47-55/2021
  • Kruel IB, Meschiatti MC, Blain GC, Ávila AMH (2015) Climatics trends in the municipality of Pelotas, state of Rio Grande do Sul, Brazil. Engenharia Agrícola 35(4):769-777. DOI: http://dx.doi.org/10.1590/1809-4430.v35n4p769-777/20185
    » http://dx.doi.org/10.1590/1809-4430.v35n4p769-777/20185
  • Kuo NW, Jien SH, Hong NM, Chen YT, Lee TY (2017) Contribution of urban runoff in Taipei metropolitan area to dissolved inorganic nitrogen export in the Danshui River, Taiwan. Environmental Science and Pollution Research 24(1):578-590. DOI 10.1007/s11356-016-7825-4
  • Leite MHS, Couto EG, Amorim RSS, Scaramuzza JF (2018) Loss of water and nutrients in diferent soil tillage systems subjected to natural rainfall in the state of Mato Grosso do Sul, Brazil. Engenharia Agrícola 38(6):864-873. DOI: http://dx.doi.org/10.1590/1809-4430.v38n6p864-873/2018
    » http://dx.doi.org/10.1590/1809-4430.v38n6p864-873/2018
  • Louzada FLRO, Xavier AC, Pessopane JEM (2018) Climatological water balance with data estimated by Tropical Rainfall Measuring Mission for Doce River basin. Engenharia Agrícola 38(3):376-386. DOI: http://dx.doi.org/10.1590/1809-4430.v38n3p376-386/2018
    » http://dx.doi.org/10.1590/1809-4430.v38n3p376-386/2018
  • Lu Q, Zhong P, Xu B, Huang X, Zhu F, Wang H, Ma Y (2022) Multi-objective risk analysis for flood control operation of a complex reservoir system under multiple time-space correlated uncertanties. Journal of Hydrology 606(3):1-17. DOI: https://doi.org/10.1016/J.JHYDROL.2021.127419
    » https://doi.org/10.1016/J.JHYDROL.2021.127419
  • Manke EB, Teixeira-Gandra CFA, Damé RCE, Nunes AB, Chagas Neta MCC, Karsburg RM (2022) Seasonal intensity-duration-frequency relationships for Pelotas, Rio Grande do Sul, Brazil. Revista Brasileira de Engenharia Agrícola e Ambiental 26(2):85-90. DOI: http://dx.doi.org/10.1590/1807-1929/Agriambi.v26n2p85-90
    » http://dx.doi.org/10.1590/1807-1929/Agriambi.v26n2p85-90
  • Mello CR, Silva AM (2009) Modelagem estatística da precipitação mensal e anual e no período seco para o estado de Minas Gerais. Revista Brasileira de Engenharia Agrícola e Ambiental 13(1): 68-74. DOI: http://dx.doi.org/10.1590/S1415-43662009000100010
    » http://dx.doi.org/10.1590/S1415-43662009000100010
  • Moriasi DN, Gitau M, Pai N, Daggupati P (2015) Hydrologic and water quality models: performance measures and evaluation criteria. American Society of Agricultural and Biological Engineers 58(6):1763–1785. DOI: http://dx.doi.org/10.13031/trans.58.10715
    » http://dx.doi.org/10.13031/trans.58.10715
  • Oliveira ARM, Borges AC, Matos AT, Nascimento M (2018) Estimation on the concentration of suspended solids from turbidity in the water of two sub-basins in the Doce River basin. Engenharia Agricola 38(5):751-759. Doi: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v38n5p751-759/2018
    » http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v38n5p751-759/2018
  • Pangle LA, Diem JE, Milligan R, Adams E, Murray A (2022) Contextualizing inflow and infiltration within the streamflow regime of urban watersheds. Water Resources Research 58(1):1-19. DOI: https://doi.org/10.1029/2021WR030406
    » https://doi.org/10.1029/2021WR030406
  • Passos JBMC, Silva DD, Lima RPC (2021) Daily rainfall disaggregation coeficients for the Doce River basin, Brasil. Engenharia Agrícola 41(2): 223-234. DOI: http://dx.doi.org/10.1590/1809-4430.v41n2p223-234/2021
    » http://dx.doi.org/10.1590/1809-4430.v41n2p223-234/2021
  • Pedron IT, Klosowski ES (2008) Distribuição de frequência de chuvas diárias no Estado do Paraná. Scientia Agraria Paranaensis 7(1-2):55-63. DOI: https://doi.org/10.18188/sap.v0i0.2052
    » https://doi.org/10.18188/sap.v0i0.2052
  • Pitt RE, Clark SE (2008) Integrated storm-water management for watershed sustainability. Journal of Irrigation and Drainage Engineering 134(5):548-555. DOI: https://doi.org/10.1061/(ASCE)0733-9437/(2008)134:5(548)
    » https://doi.org/10.1061/(ASCE)0733-9437/(2008)134:5(548)
  • Pitt RE (1999) Small storm hydrology and why it is important for the design of stormwater control practices. Journal of Water Management Modeling Chapter 4, 32 p. DOI: http://dx.doi.org/10.14796/JWMM.R204-04
    » http://dx.doi.org/10.14796/JWMM.R204-04
  • Pruski FF, Rodriguez RG, Nunes AA, Pruski PL, Singh VP (2015) Low flow stimates in regions of extrapolation of the regionalization equations: a new concept. Engenharia Agrícola 35(5):808-816. DOI: http://dx.doi.org/10.1590/1809-4430.v35n5p808-816/2015
    » http://dx.doi.org/10.1590/1809-4430.v35n5p808-816/2015
  • Rivard G (2010) Small storm hydrology and BMP modeling with SWMM5. Journal of Water Management Modeling Chapter 10. 18 p. DOI: http://dx.doi.org/10.14796/JWMM.R236-10
    » http://dx.doi.org/10.14796/JWMM.R236-10
  • Tansar H, Duan HF, Mark O (2022) Catchment-scale and local scale based evaluation of LID effectiveness on urban drainage system performance. Water Resources Management 36(2):507-526. DOI: https://doi.or/10.1007/s11269-021-03036-6
    » https://doi.or/10.1007/s11269-021-03036-6
  • Thiemig V, Rojas R, Zambrano-Bigiarini M, Levizzani V, Roo A (2012) Validation of satellite-based precipitation products over sparsely gauged African river basins. Journal of Hydrometeorology 13(6):1760-1783. DOI: https://doi.org/10.1175/JHM-D-12-0321
    » https://doi.org/10.1175/JHM-D-12-0321
  • URBAN DRAINAGE AND FLOOD CONTROL DISTRICT (UDFCD) (2011) Urban storm drainage criteria manual: Volume 3. Stormwater best management practices. Denver, 565p. Available: www.udfcd.org
    » www.udfcd.org
  • Yang L, Smith JA, Baeck ML, Zhang Y (2016) Flash flooding in small urban watersheds: Storm event hydrologic response. Water Resources Research 52(6):4571-4589. DOI: https://doi.org/10.1002/2015WRO18326
    » https://doi.org/10.1002/2015WRO18326
  • Yoon SW, Chung SW, Oh DG, Lee JW (2010) Monitoring of non-point source pollutants load from a mixed forest land use. Journal of Environment Science 22(6):801–805. DOI: https://doi.org/10.1016/51001-0742(09)601-7
    » https://doi.org/10.1016/51001-0742(09)601-7
  • Wolff W, Duarte SN, Mingoti R (2014) Nova metodologia de regionalização de vazões: estudo de caso para o estado de São Paulo. Revista Brasileira de Recursos Hídricos 19(4)21-33. DOI: http://dx.doi.org/10.21168/rbrh.v19n4.p21-33
    » http://dx.doi.org/10.21168/rbrh.v19n4.p21-33
  • Zanchin M, Moura MM, Nunes MCM, Tuchtenhahaggen IK, Lima CLR (2021) Assesment of soil loss susceptibility in Santa Rita watershed in Shouthern Brazil. Engenharia Agrícola 41(4)485-495. DOI: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v41n4p485-495/2021
    » http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v41n4p485-495/2021

Edited by

Area Editor: Alexandre Barcellos Dalri

Publication Dates

  • Publication in this collection
    27 Feb 2023
  • Date of issue
    2023

History

  • Received
    13 Jan 2022
  • Accepted
    01 Apr 2022
Associação Brasileira de Engenharia Agrícola SBEA - Associação Brasileira de Engenharia Agrícola, Departamento de Engenharia e Ciências Exatas FCAV/UNESP, Prof. Paulo Donato Castellane, km 5, 14884.900 | Jaboticabal - SP, Tel./Fax: +55 16 3209 7619 - Jaboticabal - SP - Brazil
E-mail: revistasbea@sbea.org.br