# Assessing Rainfall Erosivity with Artificial Neural Networks for the Ribeira Valley, Brazil

Authors: Reginald B. Silva; Piero Iori; Cecilia Armesto; Hugo N. Bendini
Journal: International Journal of Agronomy
 (2010)

---

## Abstract

Soil  loss  is  one  of  the  main  causes  of  pauperization  and  alteration  of  agricultural  soil  properties. Various empirical models (e.g., USLE) are used to predict soil losses from climate variables which in general have to be derived from spatial interpolation of point measurements. Alternatively, Artificial Neural Networks may be used as a powerful option to obtain site-specific climate data from independent factors. This study aimed to develop an artificial neural network to estimate rainfall erosivity in the Ribeira Valley and Coastal region of the State of São Paulo. In the development of the Artificial Neural Networks the input variables were latitude, longitude, and annual rainfall and a mathematical equation of the activation function for use in the study area as the output variable. It was found among other things that the Artificial Neural Networks can be used in the interpolation of rainfall erosivity values for the Ribeira Valley and Coastal region of the State of São Paulo to a satisfactory degree of precision in the estimation of erosion. The equation performance has been demonstrated by comparison with the mathematical equation of the activation function adjusted to the specific conditions of the study area.


---

## Body

## 1. Introduction

Erosion is considered one of the main causes of depauperation and alteration of soil properties and, consequently, of loss of agricultural soil. Mathematical models are used to quantify and/or predict such losses [1, 2]. One classical example is the universal soil loss equation (USLE), proposed by Wischmeier and Smith [3]. This equation predicts the average annual soil loss from agricultural land. The USLE is represented by the product of the following factors: (a) rainfall erosivity (R), (b) soil erodibility (K), (c) slope length (L), (d) slope percent (S), (e) soil use, handling, and coverage (C), and (f) conservative practices of soil support (P).TheR factor is an index that expresses the rainfall erosivity, in other words, its erosive capacity [4]. Erosivity is defined as the rainfall potential for soil erosion and is exclusively a function of rainfall physical characteristics, including amount, intensity of fall, droplet size, terminal velocity, and kinetic energy. Some studies have been conducted to further detail research on this erosive agent and showed that the rainfall characteristics that provide the best correlations with soil losses are intensity of fall and kinetic energy [5]. Therefore, the estimation of erosivity values, which represents the rainfall potential for erosion, is essential for planning soil and water conservation.Water erosion is causing severe problems to the population that live in the State of São Paulo, such as loss of soil from arable farmland, reduction in public investments in infrastructure works, and degradation in urban areas. Seven thousand cases of gullies, estimated currently, exist in the territory of the State of São Paulo. The cost of the corrective actions required for the stabilization of this geological phenomenon that causes severe water erosion corresponds to 20% of the State’s budget, excluding the costs of restoration of degraded urban areas, constructions, and urban street design schemes, among others. Water erosion in agricultural land is even more critical, for it is estimated that 80% of São Paulo’s agricultural soils are affected by erosion.The procedures used to estimate rainfall erosivity (R factor) are slow. Pluviograph data is necessary for the calculation of the R factor, and such data is not easily available in Brazil. The data processing and analysis is slow and complex [6]. Besides, the reduced number of meteorological stations that are equipped to provide pluviograph data makes it difficult to disseminate such studies. An alternative is the use of empirical equations based on monthly and annual rainfall to estimate the R value [1, 7, 8]. However, even with the use of the referred to equations, rainfall erosivity is restricted to sites that count on the required pluviometric data. In places where such data is not available, authors like Bertoni and Lombardi Neto [6],  Erosividade [9], and Silva [10] used interpolation techniques for R values. Machine Learning (ML) Techniques were also considered an attractive alternative for handling this problem [11].One of the main techniques within ML is the Artificial Neural Networks (ANNs). According to Persson et al. [12], an ANN is a system that is developed to imitate the operation of the human brain, which acquires knowledge through a process of training aiming at finding weights of the different connections. According to Haykin [13], this technique is based on a form of nonalgorithmic computation. For Sárközy [14], ANNs can be used as interpolation tools, and their ability to learn different input parameters makes them capable to solve complex problems from many other areas. The ANNs are also cited as alternative resources for estimating climatic variables that may replace the traditional interpolation methods [15].An ANN is composed of a set of computational elements called artificial neurons, which relate the output and input values through the following equation:(1)yij=f(∑yi’(j-1)wi’ij+bij),
where yij: output value of neuron “i” of layer “j”; n: number of neurons in the previous layer; yi’(j-1): output value of neuron i′ in the previous layer; wi’ij: value of synaptic weight of neuron “i” in layer “j”, activated by neuron “i” in the previous layer; bij: compensation value for neuron “i” in layer “j”; f: neuron “i” activation function [16].According to Saito et al. [17], consensus has not been reached on the most appropriate statistical methods for analyzing the performance of models. However, for Camargo and Sentelhas [18], model precision is given by the correlation coefficient (r), and accuracy concerns the difference between the estimated values and the observed values, given by the index of agreement (d). Another resource is the efficiency coefficient (E), which has been used by many authors in the evaluation of hydrological models and in the quantification of water constituents for estimation of its quality [19].Based on the proposition of Moreira et al. [2] and considering the specificity of climate and geomorphology, we aimed to develop an artificial neural network to estimate rainfall erosivity in the Ribeira Valley and Coastal region in the State of São Paulo.


## 2. Material and Methods

Part of the study area (Vale do Ribeira) is considered by the United Nations Educational, Scientific and Cultural Organization (UNESCO) as Mata Atlântica biosphere reserve. According to Romão [20] this region has the largest continuous tract of the remaining Mata Atlântica and the Brazilian associated ecosystems, concentrating 40% of the conservation units of the State of São Paulo. The use of around 75% of the region’s land is regulated by environmental protection acts, and 58% of these areas have been turned into public parks and ecological stations for protection—so land renting is forbidden—or into environmental protection areas where land ownership is private, though with restrictions in use.The study was conducted in the Ribeira Valley and Coastal region of the State of São Paulo, Brazil (Figure1), located between 23.77 and 24.97 latitudes and 45.42 and 49.17 longitudes, and with altitudes varying from 2 to 890 m asl. The region’s climate, according to Koppen, is the Cfa, that is, subtropical humid, with an average annual temperature of 21°C, with 28.3°C being the mean maximum temperature and 17°C the mean minimum temperature. The assessed variables included monthly precipitation, rainfall erosivity, latitude, longitude, and altitude data from 32 pluviometric stations (Table 1).Table 1
Altitude, latitude, longitude, annual precipitation, and erosivity of 32 municipalities of the Ribeira Valley and Coastal region of the State of São Paulo.
MunicipalityAlt. (m)Lat.Long.Annual precipitation (mm)EI30 [7] MJ mm (hha)-1(1) Apiaí87024.4248.651613.538612(2) Barra do Chapéu77024.4749.021506.128029(3) Barra do Turvo66024.9748.331986.2010320(4) Bertioga323.8462295.7011455(5) Cajati8024.7848.181492.188848(6) Cananéia724.9347.952410.0912410(7) Cubatão523.8746.382774.8313075(8) Eldorado2024.5248.11602.269273(9) Guarujá223.9546.32005.7410469(10) Iguape624.747.531567.588401(11) Ilha Comprida224.947.82299.1611469(12) Itanhaém75023.9546.81950.919746(13) Itaoca18024.6548.831327.627467(14) Itapirapuã58024.5749.171549.348323(15) Itariri5024.2847.182249.4311887(16) Jacupiranga9024.7248.021506.098673(17) Juquiá2024.147.681794.729768(18) Juquitiba73023.9746.91897.2410496(19) Miracatu3024.2547.371582.489765(20) Mongaguá2024.0846.622612.1112685(21) Pariquera Açu3024.7247.881535.128697(22) Pedro de Toledo6024.2847.231521.218594(23) Peruíbe324.3247.022121.7012029(24) Praia Grande1024.0346.552655.4713188(25) Registro6024.447.751760.539702(26) Ribeira16024.6549.021343.887470(27) Santos20023.8846.223615.3815919(28) São Lourenço89023.7846.931567.339409(29) São Sebastião2023.7745.421305.287988(30) São Vicente1023.9746.372432.4512015(31) Sete Barras2024.3847.931582.538921(32) Tapiraí87023.9747.51819.589764Maximum value890.00024.97049.1703615.38015919.000Minimum value2.00023.77045.4201305.2807467.000Average value225.25024.31347.4571915.11810152.094Standard deviation324.7300.3710.957518.3011973.046Figure 1
Location of the study area: Ribeira Valley and Coastal region of the State of São Paulo, Brazil.TheR values for each station were obtained through (2), proposed by Silva et al. [8]. The sum of EI30 values for the 12 months of the year represented the rainfall erosivity value: (2)EI30=206,67+40,674(Pm2Pa-1),
where EI30 = erosion index (MJ mm (hha)-1); Pm = average monthly precipitations (mm); Pa = annual precipitations (mm).TheWeka(Waikato Environment for Knowledge Analysis) package was used to develop the network, which is formed by a set of implementations of algorithms of various Data Mining techniques. The input variables were composed of the latitude and longitude values of each station (in decimal degrees), the altitude value (m), as well as the average annual precipitation (mm). A linear activation function was used in the output variable to obtain the rainfall erosivity value (R factor), in MJ mm (hha)-1 of the location represented by the input vector.The analysis was performed by means of ML study. The results were obtained using 10-fold  cross-validation, with 10 repetitions. Multilayer PerceptronANNs were used for training, with backpropagation algorithm, momentumterm, and learning rate of, respectively, 1 and 0. One-layer ANNs trained with the Least Mean Square (LMS) algorithm were also used. The root mean squared error (RMSE) and the correlation coefficient (r) obtained in the test sets were used.The correlation coefficient (r), the agreement index (d) proposed by Willmott [21], and the confidence interval (c), which is the product of r and d, were used for the evaluation of the performance of the model obtained with the use of ML. The agreement index can be calculated by (3) as follows:(3)d=1-(∑(Rei-Ri)2∑([Ri-Rex]+[Rei-Rex]))2,
where i  is the number of observations, R is the experimentally observed value, Re is the value estimated by the model, and Rex is the average value of the experimentally observed values.The efficiency coefficient (E) defined by Nash and Sutcliffe [22] was also calculated as in (4); the mean absolute error (MAE) was showed by Legates and McCabe Jr. [23] as in (5); the mean percent error (MPE) was used by Chong et al. [24] as in (6); and the root mean squared error (RMSE) was represented by (7) as follows:(4)E=1-∑(Rei-Ri)2∑(Rei-Rex)2,(5)MAE=∑[Rei-Ri]J,(6)MPE=∑([Rei-Ri]Rei-1J)100,(7)RMSE=√∑(Rei-Ri)2J.The spatial distribution of erosivity (R factor) estimated by Silva et al. [8] and obtained through the use of ANNs was performed through kriging interpolation method, using Surfer version 8.0 application.


## 3. Results and Discussion

The input variables composed by the latitude, longitude, altitude, and the average annual precipitation values for each station made it possible to adjust (8) to estimate rainfall erosivity (R factor). This equation obtained with the use of ANN was a multivariate linear function, and its determination coefficient (R2=0.97) corroborates its performance in estimating rainfall erosivity with 95% confidence interval. The RMSE value for LMS algorithm (496.43) was lower compared to backpropagationalgorithm (738.28), which corroborates [13] that the LMS algorithm provides better solution for linear problems: (8)EI30=(-0,4362Alt)+(224,7948Lat)+(-82,0386Long)+(3,6577Pa)+1572,3404(R2=0,97),
where Alt: altitude (m), Lat: latitude (decimal degrees), Long: longitude (decimal degrees), and Pa: annual precipitation (mm).Table2 shows a summary of the descriptive statistical analysis of the performance of the equation proposed by Silva et al. [8] and the performance obtained through the use of the ANN in the estimation of erosivity. The main parameters were the correlation coefficient (r), agreement index (d), efficiency coefficient (E), mean absolute error (MAE), mean square error (MSE), root mean squared error (RMSE), and the mean percent error (MPE).Table 2
Summary of the descriptive statistical analysis of erosivity values estimated by Silva et al. [8] and with the use of ANN in the Ribeira Valley and Coastal region of the State of São Paulo.
IndexesResultCorrelation coefficient (r)0.97Agreement index (d)0.99Confidence interval (c)0.96Efficiency coefficient (E)0.95Mean absolute error (MAE)294.91Mean square error (MSE)168711.78Root mean squared error (RMSE)410.75Mean percent error (MPE) (%)2.95The values found for the correlation coefficient and the agreement index were, respectively, 0.97 and 0.99. As it is known, the more precise and accurate results forr and d, respectively, are those closer to 1.0. Moreover, a 0.95 value was found for E, and it is known that E values may vary from minus infinity to 1.0, with 1.0 value indicating a perfect simulation of the model.The lower the percentage error, the higher the efficiency of the method used. According to literature, values close to 10% are too low for accurate estimation of rainfall erosivity [2]. In this study, the MPE value was 2.95%; thus, it was found to be acceptable. The confidence interval was 0.96, which was found to be an optimum performance, because according to Camargo and Sentelhas [18] c-values higher than 0.85 ensure optimal performance.The variability in altitude values shown in the dispersion in data in relation to the mean (Table1) can explain the low synaptic weight (−0.4362) of the neuron (wi′ij) and the contribution of this variable to (8). This result is different from the one found in the studies of Moreira et al. [2], who stresses the influence of altitude and continentality in rainfall erosivity throughout the State of São Paulo, obtained with the use of the ANN. However, it is important to emphasize that the referred to authors used the model proposed by Lombardi Neto and Moldenhauer [7] as the mathematical equation of the activation function, as well as data from meteorological stations distributed throughout the State of São Paulo, so that the variability (dispersion) of altitude data in the study area has become more dilute.The average annual value of rainfall erosivity estimated for the Ribeira Valley and Coastal region of the State of São Paulo, based on the equation of Silva et al. [8], was 10152 MJ mm (hha)-1 varying from 7467 to 15919 MJ mm (hha)-1. Nevertheless, the average was 10051 MJ mm (hha)-1 and varied from 7885 to 16285 MJ mm (hha)-1 (Table 3) when the equation obtained through the ANN was used. Variance analysis between the groups of erosivity data estimated by the different equations did not show significant differences (P=.83).Table 3
Results of the estimation of rainfall erosivity (R) according to [8] (MJ mm (hha)-1), of the estimation with the use of the equation obtained with the use of artificial neural networks (Re) (MJ mm (hha)-1), and of the relative percentage error (RPE) in module between the erosivity values R and Re, in %, for the 32 municipalities in the Ribeira Valley and Coastal region of the State of São Paulo.
MunicipalityRReRPEApiaí861285931,00Barra do Chapéu802982250,98Barra do Turvo10320101981,01Bertioga11455115440,99Cajati884886131,03Cananéia12410120551,03Cubatão13075132810,98Eldorado927389901,03Guarujá10469104931,00Iguape840189570,94Ilha Comprida11469116570,98Itanhaém974699250,98Itaoca746778850,95Itapirapuã832384760,98Itariri11887113661,05Jacupiranga867386591,00Juquiá976896341,01Juquitiba1049697341,08Miracatu976589131,10Mongaguá12685127061,00Pariquera Açu869788030,99Pedro de Toledo859486940,99Peruíbe12029109411,10Praia Grande13188128641,03Registro970295531,02Ribeira747079380,94Santos15919162850,98São Lourenço940984121,12São Sebastião798879551,00São Vicente12015120491,00Sete Barras892189001,00Tapiraí976493401,05Maximum value15919,0016285,000,94Minimum value7467,007885,001,12Average value10152,0910051,191,02Standard deviation1973,051945,880,09The equation obtained through the ANN showed higher rainfall erosivity values in 15 of the 32 municipalities of the study area when compared to the values estimated by Silva et al. [8] (Figure 2). The maximum RPE value was 1.12% (Table 3), and according to Moreira et al. [2] an RPE value lower than 10% is too low to provide accurate estimation of rainfall erosivity. This result is different from the ones showed by Moreira et al. [2], where the ANN tended to overestimate rainfall erosivity values. As said by the authors, the advantage of using this behavior pattern is explained by the adoption of soil use and management practices that keep soil loss within tolerable limits. However, to ensure that this occurs, the other factors of the Universal Soil Loss Equation (USLE) must be constant.Figure 2
Erosivity values estimated by Silva et al. [8] and by the ANN for each municipality.The spatial distribution of rainfall erosivity obtained by Silva et al. [8] and through the ANN for the Ribeira Valley and Coastal region of the State of São Paulo is shown in Figure 3. Erosion variability obtained through the ANN follows the same distribution pattern of the variability found by Silva et al. [8], corroborating classical statistical analyses. In conformity with the findings of Moreira et al. [2], we noticed that the highest R values occurred in coastal municipalities and gradually decreased in inland areas. The erosivity interval comprised between 7500 and 10500 MJ mm (hah)-1 was the highest occurrence value (Figure 3).Spatial distribution of rainfall erosivity (MJ mm(hha)-1) for the Ribeira Valley and Coastal region of the State of São Paulo estimated through ANN and by Silva et al. [8].
(a)(b)


## 4. Conclusions

The use of artificial neural networks was found to be viable for the interpolation of rainfall erosivity values for the Ribeira Valley and Coastal region of the State of São Paulo and can be used to a satisfactory degree of precision in the estimation of erosionThe best-fit equation generated by the ANN was a multivariate linear function, with its performance proven when compared to the equation of the activation function adjusted to the conditions of the study area.



---

*Source: 365249-2010-08-24.xml*