Adaptive Prediction of Water Droplet Infiltration Effectiveness of Sprinkler Irrigation Using Regularized Sparse Autoencoder–Adaptive Network-Based Fuzzy Inference System (RSAE–ANFIS)

: As the high productive efficiency of sprinkler irrigation is largely based on balanced soil moisture distribution, it is essential to study the exact effectiveness of water droplet infiltration, which provides a theoretical basis for rationally scheduling the circulation efficiency of groundwater in agricultural irrigation performance. This research carried out adaptive prediction of the droplet infiltration effectiveness of sprinkler irrigation by using a novel approach of a regularized sparse autoencoder–adaptive network-based fuzzy inference system (RSAE–ANFIS), for the purpose of quantifying actual water droplet infiltration and effectiveness results of precision irrigation in various environmental conditions. The intelligent prediction experiment we implemented could be phased as: the demonstration of governing equations of droplet infiltration for sprinkler irrigation modeling; the measurement and computation of probability densities in water droplet infiltration; innovative establishment and working analysis of RSAE–ANFIS; and the adaptive prediction of infiltration effectiveness indexes, such as average soil moisture depth increment ( θ , mm), irrigation infiltration efficiency ( e a , %), irrigation turn duration efficiency ( e t , mm/min), and the uniformity coefficient of soil moisture infiltration ( C u , %), which were implemented to provide a comprehensive illustration for the effective scheduling of sprinkler irrigation. Result comparisons indicated that when jetting pressure ( P w ) was 255.2 kPa, the impinge angle ( W a ) was 42.5 ° , the water flow rate ( F a ) was 0.67 kg/min, and continuous irrigation time ( T c ) was 32.4 min (error tolerance = ± 5%, the same as follows), thereby an optimum and stable effectiveness quality of sprinkler irrigation could be achieved, whereas average soil moisture depth increment ( θ ) was 57.6 mm, irrigation infiltration efficiency ( e a ) was 62.5%, irrigation turn duration efficiency ( e t ) was 34.5 mm/min, and the uniformity coefficient of soil moisture infiltration ( C u ) was 53.6%, accordingly. It could be concluded that the proposed approach of the regularized sparse autoencoder–adaptive network-based fuzzy inference system has outstanding predictive capability and possesses much better working superiority for infiltration effectiveness in accuracy and efficiency; meanwhile, a high agreement between the adaptive predicted and actual measured values of infiltration effectiveness could be obtained. This novel intelligent prediction system has been promoted constructively to improve the quality uniformity of sprinkler irrigation and, consequently, to facilitate the productive management of sprinkler irrigated agriculture.


Introduction
It is well known that sprinkler irrigated agriculture uses most of the water resources in the world, and it has already been investigated for several decades due to the remarkable population increment, extreme environmental climate, and excessive water resource demands. Since the high performance of sprinkler irrigated agriculture is largely based upon appropriate droplet infiltration, this makes it important to study the exact effectiveness of water infiltration in the interests of understanding balanced water resource distribution and providing a theoretical basis for rationally scheduling the circulation efficiency of groundwater in precision irrigation. Therefore, the adaptive prediction of droplet infiltration effectiveness is necessary to improve irrigation monitoring and ensure crop growth. Here, sprinkler irrigation can be described as an efficient collaboration and precise application of water spraying in tiny droplets on plants, soil or grass, meeting specific requirements of individual plants or management units, and minimizing the related adverse environmental impacts simultaneously. As a holistic collaborative performance, sprinkler irrigation not only enables the emergence of more flexible crop cultivation systems, but also involves the optimal management of agricultural irrigation systems, so that it has been widely employed to solve the growing problem of water scarcity and the need for more effective irrigation practices. Accordingly, the droplet infiltration effectiveness of sprinkler irrigation could be defined here as the working quality of liquid water slowly passing through the filtering soil medium, for the purpose of evaluating the balanced distributive capability of droplet infiltration in soil fields and facilitating the cultivation equilibrium of plants and crops. Mostly, the quantitative calibration of droplet infiltration effectiveness depends upon the complicated inter-relation mechanism among experimental irrigation conditions, water flow properties, and practical implementations of sprinkler irrigation, etc. Theoretical research and working characteristics discussed on this topic can be found from the literature [1,2]. Unfortunately, the adaptive prediction of infiltration effectiveness has rarely been thoroughly investigated. Prompted by the rapid progress of sprinkler irrigation and flow mechanics, the mechanism complexity and investigation depth of droplet infiltration already have traditional process simulations and pure physical comparisons beyond original knowledge scopes [3]. It is regrettable that most current studies cannot solve the existing problems of infiltration effectiveness prediction accurately.
With respect to droplet infiltration quantification, controlled drainage has already been promoted to boost crop yields and reduce subsurface drainage flows or leaching nutrients in the irrigated farm fields [4]. Since uneven water distribution, the misuse of water, and inefficient watering techniques are some of major causes wreaking devastation on water security, Mostafa et al. managed drip irrigation as an effective tool for increasing crop yield, enlarging water productivity, and saving irrigation water for newly reclaimed areas [5]. Recently, a new promising approach integrating soil moisture content measurements with the Community Land Model (CLM) using sequential data assimilation (DA) has been presented to improve the prediction quality of soil water status and efficiently design irrigation strategies [6]. Although the control model of mobile predictive provides novel analysis tools for droplet infiltration, which calibrates the cooperative performance of precision irrigation [7], more comprehensive evaluation approaches concerning droplet infiltration effectiveness are urgently needed for sprinkler irrigation to promote its working efficiency and prediction reliability.
For the purpose of infiltration effectiveness evaluation, Smith et al. used artificial neural networks to estimate the furrow irrigation duration on clay soils [8]. Oker et al. focused on maize production under mobile drip irrigation, and the monitoring of adaptive irrigation efficiency [9]. Meanwhile, an integrated hydrological irrigation optimization modeling system could be learned from the literature [10]. Salem et al. have employed artificial neural network and regression models to calibrate the complicated impact of climate change on groundwater infiltration and irrigation effectiveness [11]. Thereafter, a novel prediction model has been proposed for supplemental irrigation water rights and climate change adaptation in the literature [12]. In the past several years, it is an obvious tendency that more and more intelligent computation system and fuzzy decision-making approaches have been introduced into this field, for instance, the ant colony algorithm, genetic optimization, the simulated annealing-genetic algorithm, Taguchi estimation, artificial neural network-simulated annealing prediction, and the genetically optimized neural network, etc. They have been applied to reach an inspiring goal of high-efficient scheduling and accurate prediction of precision irrigation performance [10][11][12][13][14]. From the above mentioned publications, it could be learned that all these investigations and endeavors give us reliable knowledge to assess the actual effects of soil-moisture interaction and the distribution properties of infiltration rate. Moreover, their inherent superiorities of quantifying the working efficiency of soil moisture infiltration could also be confirmed. The accurate prediction of water droplet infiltration effectiveness should be the priority of our research efforts.
Here in this research, adaptive prediction of the water droplet infiltration effectiveness of sprinkler irrigation has been investigated to give an effective benchmark of soil moisture evaluation and solve the challenges of the accurate calibration of water droplet infiltration infield. To achieve this goal, this paper presents a set of novel methodologies in the following sections.

Calibration of Droplet Infiltration Effectiveness
During the actual process of sprinkler irrigation, the moisture infiltration rate decreases when irrigated soil becomes water-saturated, thus the traditional calibration of irrigation infiltration cannot reveal the complicated correlation between the droplet infiltration mechanism and soil moisture conditions, especially in such domains as moisture distribution and flow filtering. Then, the effectiveness of water droplet infiltration was introduced for the calibration of soil humidity, and it could be determined by a completed setup of a moisture sensor array. Figure 1 demonstrates the water droplet infiltration of sprinkler irrigation and the planar meshed grid zone on the soil field, being used for the moisture gauging and position calibration of droplet infiltration effectiveness. According to this setup arrangement, the coordinate positions of droplet infiltration points, conforming to the normal distribution of N(μ,σ), were presented as: where xk and yk are the horizontal and vertical coordinates of specific positions on the ground where water droplet infiltration happens, being denoted as the (xk,yk)th meshed grid zone. Mk and Nk denote the number of meshed grids in length and width directions, presupposed as 200 and 400, respectively: 1 ≤ xk ≤ Mk, and 1 ≤ yk ≤ Nk. N(μ,σ) denotes the normal coordinate distribution of droplet infiltration positions, while μ is the mean value of the infiltration position coordinate, and σ stands for its corresponding distribution variance. The probability density of water droplet infiltration at the investigated area can be defined as: In order to calibrate the working performance of sprinkler irrigation, the following infiltration effectiveness indexes have been proposed accordingly. Here, the average soil moisture depth increment (θ, mm) due to sprinkler irrigation or net infiltrated moisture depth is calculated as: where θantecedent and θend denote average antecedent and achieved soil moisture depths (mm), respectively, which propose an accurate calibration of the infiltration rate from the novel perspective of spatial moisture distribution. Although soil moisture depth is frequently acknowledged as an important contributing factor to infiltration effectiveness, the complex link between deep water drainage and soil infiltration efficiency has not been fully recognized. More particularly, they do not indicate the key role of irrigation infiltration in controlling the working efficiency of water drainage. In order to solve this problem, irrigation infiltration efficiency (ea) is proposed to evaluate the working efficiency of soil moisture infiltration into the soil field. This mathematical parameter is defined as the ratio of the net moisture infiltrated depth to the cumulative volume of water consumption, including the parametric variables of cumulative irrigation flow depth (Cf, mm) and cumulative evaporated water flow depth (Ev, mm), in total: where Ev denotes the cumulative water flow depth evaporated at the end of irrigation turn (mm), and θ denotes the average soil moisture depth increment being determined by Equation (4). Based on these preparations, the index of irrigation infiltration efficiency (ea) could be determined by measuring the moisture infiltration depth, the cumulative irrigation flow depth, and the consumed volume of sprinkler irrigation water added to the infield zone. The ratio between them could be defined as the % of irrigation infiltration efficiency during the sprinkler irrigation process. Meanwhile, the concept of irrigation turn duration efficiency, or the ratio of net infiltrated moisture depth per unit of irrigation time (et, mm/min), can be obtained as: where θ is the average soil moisture depth increment (mm) and t is the irrigation turn duration time (min). Water infiltration uniformity refers to the appropriate spatial distribution of infiltrated moisture, and the uniformity coefficient of soil moisture infiltration (Cu, %) could be proposed here and calculated as: where θxk,yk(%) stands for the soil moisture mass proportion at one infiltration measurement grid zone coordinated by (xk,yk), and end θ stands for the average value of soil-moisture mass proportion concerning all infiltration measurement grid zones (%), so that Cu provides a useful tool to evaluate the uniformity state of soil moisture infiltration from novel perspectives. Since the ideal situation regarding infiltration uniformity occurs when all field plants receive the same amount of sprinkling water, we could determine the necessary water supply that produces optimal uniformity coefficients with high computation efficiency. The uniformity coefficient of soil moisture infiltration (Cu, %) is typically expressed in terms of the mathematical difference between the measured and averaged soil moisture mass proportion, and then divided by the area of the objective meshed grid zone, which accurately illustrates the spatial distribution of sprinkler irrigation water infield.
As all these parametric indexes demonstrate the inherent characteristics of soil moisture infiltration from a series of novel effectiveness perspectives of sprinkler irrigation infiltration, they could be applied effectively in the following adaptive prediction of droplet infiltration effectiveness.

Experimental Environment of Sprinkler Irrigation
This sprinkler irrigation experiment was carried out at the irrigation field of Guangzhou University (23°03′1.23″ N, 113°24′3.92″ E) with an area of 20 ha, from 5 September 2019 to 4 September in 2020 and from 8:00 to 18:00 every day. Figure 2 gives the exact location of the experimental field where this irrigation experiment was conducted. The soil of the experimental site has a pH value range of 6.62-6.78. Since the electrical conductivity of soil (ECe) provides a standard benchmark to evaluate soil salinity, it could be measured by using a mobile SUP-TDS210-B electromagnetic induction type soil electrical conductivity meter, which is manufactured by Supmea Automation Co., LTD (Hangzhou City, Zhejiang Province, China). In this research, we analyze 20 soil samples collected from 20 reclaimed grid positions, thus the measurement of ECe can be realized as 0.14-0.18 dS/m for experimental soil fields of different textures and salinity levels, which enhances the interpretation of soil electrical conductivity maps and improves the scheduling of site-specific sprinkler irrigation. Next, the average volumetric soil moisture content at field capacity is 40%, and the wilting point of the root zone is 30%. Here, the organic content of the soil field is 22.60 g/kg, the nitrogen content is 1.68 mg/kg, the field capacity is 26.3%, and its mean bulk density varies from 1.35 g/cm 3 to 1.55 g/cm 3 . Guangzhou has a subtropical climate where drought is rarely seen, with a mean annual temperature of 21.5-22.6 °C, precipitation of 1623.6-1899.8 mm, and evapotranspiration of 1603.5 mm, according to the past 60 years of weather records referenced from a local meteorological agency. After detailed considerations and plant comparisons, Zoysia matrella was selected as the experimented plant for sprinkler irrigation, attributed to its robust climate and environmental adaptation capability, making it ideal for irrigation infiltration prediction and growth control in our soil field.

Working Mechanism of Regularized Sparse Autoencoder (RSAE)
RSAE (Regularized Sparse Autoencoder) is an excellent autoencoder that can be used for calculating the probability density from a set of measured raw data of droplet infiltration rates. It is highly valued by incorporating feature extraction and process quantification into a general-purpose fuzzy prediction system [13,14]. According to this approach, the probability density of water droplet infiltration at the objective meshed grid zone is denoted as unlabeled data where WRSAE denotes the weight matrix of the regularized sparse autoencoder and σr(z) is a rectified linear unit employed as the activation function in the RSAE decoder. Since it allows for more efficient training of network prediction than any other traditional function does, σr(z) encourages sparse activation and be described as: To learn various features from the measured probability density data of droplet infiltration, the cost function of RSAE is presented by Equation (12).
( ) where WRSAE1 and WRSAE2 denote the weight matrix of the regularized sparse autoencoder, ensuring its computational performance is kept in a highly efficient state by using the cost function of RSAE. λxkyk stands for the regular coefficient of RSAE, corresponding to the grid zone coordinated by (xk,yk). Here, WRSAE1 and WRSAE2 are replaced by W and W T , then the optimization problem of the RSAE network minimizes the following cost function: Thereafter, we calculate the gradient of JRSAE with respect to W: Here, D = W T σr(W•(xk,yk)) − (xk,yk), sgn denotes the result of the sign function of σr(W•(xk,yk)), σr' denotes the derivative function of rectified linear unit σr, and (xk,yk) stands for the matrix form of (xk) and (yk), so that the reiterative update process of W (i+1,j+1) for RSAE can be described as: where i and j are the ith and jth intervals of the update process of WRSAE1 and WRSAE2, respectively, H (i,j) denotes the inverse of the Hessian matrix [17,18], and η denotes the step size of the update process accordingly. In this equation, η (i+1,j+1) can be determined iteratively by: The weight matrix WRSAE could be regularized by an orthonormality constraint. RSAE is trained an adequate number of times so that the probability density of water droplet infiltration can be accurately calculated and inputted into the Adaptive Neuro Fuzzy Inference System (ANFIS) for the highly efficient prediction of infiltration effectiveness, which provides high-performance adaptive prediction that incorporates neural networks with fuzzy logic principles, making ANFIS a universal and reliable prediction network for sprinkler irrigation effectiveness [19,20].

Working Mechanism of ANFIS (Adaptive Network-Based Fuzzy Inference System) Integrated with RSAE
After the logic network of RSAE has been established, an intelligent prediction system based on the working principle of the adaptive network-based fuzzy inference system (ANFIS) was developed to predict such droplet infiltration effectiveness indexes as average soil moisture depth increment (θ, mm), irrigation infiltration efficiency (ea, %), irrigation turn duration efficiency (et, mm/min), and the uniformity coefficient of soil moisture infiltration (Cu, %). Here, RSAE identifies the mathematical properties of the probability density of droplet infiltration and ANFIS possesses an excellent capability of adaptive prediction for infiltration effectiveness, which are characterized by their mutual influences and fuzzy correlations [21][22][23]. As Figures 3 and 4 demonstrate the logic diagram of RSAE-ANFIS, this intelligent prediction system could be divided into an RSAE monitoring level and an ANFIS prediction level. The RSAE monitoring level includes the gauging system that quantifies the measured probability density of droplet infiltration. On this level, all data resources of irrigation infiltration connect to the high-performance computation units, such as the data-processing units of monitoring, analyzing, transmitting, and decision making, which are located on the ANFIS adaptive prediction level. By linking the RSAE monitoring level and the ANFIS prediction level together through a set of information transmitting and processing modules taking charge of cognitive calculation, parametric monitoring and operative collaboration, a highly-efficient adaptive prediction system for droplet infiltration effectiveness can be established.
In this research, a sensor array was used for water infiltration measurement. It stably indicates the true moisture content of the soil field with high accuracy and reliable performance. This experimental calculation has been implemented on the computation platform of irrigation effectiveness prediction in the Guangzhou University intelligent irrigation lab. The computation platform we used for RSAE-ANFIS prediction was equipped with an Intel core i9 7900X microprocessor, 32GB of RAM (Random Access Memory), MATLAB R2018b, and Windows 10. Taking a typical case of adaptive prediction, for example, RSAE collects the input data of water droplet infiltration. Herein, the droplet infiltration rates and soil humidity signals are measured from the array of moisture sensors (9 × 9 in each grid zone illustrated by Figure 1), embedded in different soil depths (0-500 cm) and connected with the informative data processing center of the sprinkler irrigation system. Meanwhile, a convective lognormal transfer function was employed to picture the infiltration-stress curves from the measured soil moisture data. Then, RSAE calculates the probability density of droplet infiltration at the targeted grid zone, and operates as a reliable cognitive process classifier to recognize the mathematical characteristics of probability density from various kinds of mass infiltration data [24,25]. After all these computation tasks are completed, ANFIS sends back amend actions/decision signals to the RSAE monitoring level through data representation of effectiveness prediction results. Thereafter, an effective calculation mechanism called adaptive prediction was applied, in which this mechanism architecture compares the calculated and actual measured results of infiltration effectiveness indexes and publishes adaptive instructions to amend the initial predicted value to a more accurate one. Based on this arrangement, the RSAE-ANFIS computation platform accepts the probability density, processes the data signals, and selects appropriate calculation steps to predict infiltration effectiveness indexes (θ, ea, et, and Cu). The obtained effectiveness results could be classified into two groups: actual measured results characterized by stochasticity and nonlinearity, and the calculated results characterized by uniformity and equilibrium [26]. Major differences between them comprise deviation corrections; therefore, a set of specific emendation instructions was presented and operated by sequence. It could be seen that the probability density of water droplet infiltration calculated by RSAE helps ANFIS to obtain higher and stabler computation accuracy for predicting irrigation infiltration effectiveness. Therefore, the high precision and strong robustness of adaptive effectiveness prediction can be ensured.

Sprinkler Irrigation System
In order to realize the adaptive prediction of infiltration effectiveness, it is necessary to build up an applicable sprinkler irrigation system at first. Figure 5 illustrates the tested irrigation system and the spatial layout of the pipeline network. Herein, the sprinkler irrigation system integrates a completed set of mechanical components together, such as the high-performance 350HW-8S fluid-intensifier pumps manufactured by Yancheng Harriston Int'l Co., Ltd. (Yancheng city, Jiangsu Province, China), the revolving HX-3301-1 sprinkler nozzles provided by Taizhou Hengxin Valve Tech. (Taizhou city, Jiangsu Province, China), the NPS2 CLASS150 pressure valves manufactured by Sichuan Saier Valve Mfg. Co., Ltd. (Zigong city, Sichuan Province, China), the irrigation PVC pressure pipes manufactured by Jiangsu Haiwei Plastics Industry Technology Co., Ltd.(Jiangyin city, Jiangsu Province, China), and the BONAD RS485 soil moisture sensors (BND-MS10) manufactured by Shenzhen Bonad Precision Instrument Co., Ltd. (Shenzhen city, Guangdong Province, China), etc. In our irrigation pipeline network, the experimented soil field was equally partitioned into different meshed grids conforming to Figure 1, where each grid zone used for sprinkler irrigation could be positioned by (xk,yk) or being denoted as the (xk,yk)th meshed grid zone (1≤ xk ≤ Mk, and 1 ≤ yk ≤ Nk). Since in each grid the volume of surface water was different, according to the measured soil moisture data from the (xk,yk)th grid, the sprinkler irrigation system controls the water supply to a specific grid zone through adjusting the set of water distribution valve modules. Consequently, irrigation water was pressurized into the fluid intensifier pump and stabilized by the flow accumulator, and then appropriate water droplet concentration on every grid zone could be maintained constantly. A moderate flow pressure of 320 kPa was applied in this research to ensure that the infiltrated soil depths could be measured easily. A relatively high nozzle ratio of 6.5:1 (nozzle tube length: nozzle orifice diameter) allows for the effective isolation of stream disturbance from sprinkling water flow. In this case, a pair of sprinkler irrigation heads (the revolving sprinkler nozzles) with a vertical height of 80 mm were placed at a linear interval of 2.5 m to provide an equally distributed soil wetting pattern and guarantee a wetting-diameter overlap of 55-60% all the time. During this experimental irrigation process, the soil moisture content or droplet infiltration rate was checked once every 10 min. To ensure the high repeatability of testing conditions, Table 1    In Figure 6, for the accurate prediction of infiltration rates, clods of soil crust with a thickness of 10 cm on different soil depth layers (vertical depth ranged from 0 to −500 cm) were sampled out and observed closely after each irrigation test was completed. The average distribution of water precipitation was measured at a rectangular area of 15 × 15 m 2 in calm conditions. Such environmental parameters as average air temperature (T, °C), average air relative humidity (RH, %), and solar radiation (SR) were collected every 10 min. Here, local air temperature was measured by an HW-F7 thermometer (LD Products Inc. California, USA), air relative humidity was determined by using Testo 610 hygrometers (Testo Pty Ltd, Croydon South, Australia), and solar radiation value could be obtained from an LI-1500 irradiance radiometer (HUATEC. Beijing, China). Since wind speed (U, m/s) has a decisive influence on the evaporation loss of soil moisture, it was measured by RK100-02 RS485 wind speed measurement units (RIKA Sensor Inc., Hunan, China). On the other side, the soil moisture evaporation measurements (E, mm) were carried out by 40 sets of newly developed microlysimeter, which equally distributed throughout the whole field area. All these parametric variables were measured at row and inter-row irrigation positions, and then they were averaged across all measurement positions to provide accurate values.   The probability density of droplet infiltration at each grid zone could be computed by RSAE, whereas spline interpolation was used to calculate the infiltration rates at any coordinate position that could not be measured or observed directly. Thereafter, the standard deviations and mean values of the probability density were acquired by sequence [27,28]. According to the infiltration gradient maps highlighted by tests 1-5, Figure  6 pictures the obtained droplet infiltration distribution and its partitioned levels denoted from I to X, based on the measured soil moisture content. For instance, Level I denotes a soil moisture saturability of 0-10%, Level II denotes 10-20% … and Level X denotes 90-100%, etc. Here, specific gradient values of droplet infiltration per cubic space were demonstrated by the scale column on the left, and their scale unit was mm/1000 mm 3 . The photos in the right column denote the actually obtained images of soil ground in different infiltration rates leveled from I to X, describing the on-the-spot variation of droplet infiltration effectiveness.

Adaptive Prediction of Droplet Infiltration Effectiveness
To control the efficiency of sprinkler irrigation, a set of input variables for RSAE-ANFIS prediction was selected, including jetting pressure (Pw/kPa), impinge angle (Wa/°), flow rate of water (Fa/kg/min), and continuous irrigation time (Tc/min). Accordingly, average soil moisture depth increment (θ/mm), irrigation infiltration efficiency (ea/%), irrigation turn duration efficiency (et/mm/min), and the uniformity coefficient of soil moisture infiltration (Cu/%) were applied to profile moisture infiltration [29,30]. Based on the measured coordinates of droplets, the infiltration effectiveness indexes of every grid zone were transmitted to the database and application server of the RSAE-ANFIS network. If mathematical deviations and computational errors emerge, the novel algorithm of adaptive prediction is able to provide corrective feedback guidance to solve them. For the convenience of fuzzy prediction, the input and output variables should be regulated to an identical data range. Table 2 demonstrates the partitioned levels of irrigation parameters and infiltration effectiveness indexes, being used as input and output variables in the RSAE-ANFIS system. The combination of input (Pw, Wa, Fa, Tc) and output (θ, ea, et, Cu) illustrates a set of fuzzy correlations described by logic rules. To predict infiltration effectiveness, 25% of experimental cases (50 turns) were used for network training and others for network testing. Each case was tested for 10 trials to reduce signal disturbances and random errors, by which the resultant effectiveness indexes were averaged by these trials. Table 3 shows the logic rules of the RSAE-ANFIS system. The irrigation parameters are the input levels, and the partitioned levels of effectiveness indexes are the output levels. In this table, the logic rule base was selected via fuzzy reasoning algorithms. As their partitioned levels were referenced from Table 2, original parametric assessment and effectiveness prediction could be simplified to a quantified profile of logic level reasoning. For the purpose of enhancing the applicability and precision of ANFIS computation, the inherent properties of fuzzy logic rules by setting IF-THEN rules were inherited in prediction modules. Take No. 5 cases for instance, IF the experimental run of the input partitioned level of Pw is 2, Wa is 3, Fa is 5, and Tc is 2, THEN the experimental run of the output partitioned level of θ is 10, ea is 1, et is 2, and Cu is 1, accordingly, through a set of logic reasoning and statistical computing. Based on these logic rules and input/output levels, the infiltration effectiveness indexes in different cases could be determined, confirming the reliability and usefulness of RSAE-ANFIS prediction. Considering the current capability of our computation platform, Figure 7 presents the profile comparisons between the predicted and actual measured effectiveness indexes at different soil depths, and Table 4 gives the data comparison. Tests A-K have been implemented sequentially on specific monitoring grid zones described by Figure 1; thereafter, the infiltration effectiveness indexes were measured at different soil depths, denoted by Layer 1: 0-100 cm; Layer 2: 100-200 cm; Layer 3: 200-300 cm; Layer 4: 300-400 cm; and Layer 5: 400-500 cm. When the predicted and measured results (θ, ea, et, Cu) were profiled, the relative errors (RE) between them could be determined to quantify the accuracy and reliability of predicted effectiveness: RE = ((Predicted index − Measured index)/Measured index) × 100%. It was learned from Table 4 that RE kept in a relatively small range of 1.8-7.8%, confirming the high precision of the intelligent prediction system.

Assessments of Adaptive Prediction Quality
The following evaluation indexes were designed or introduced specifically to assess the prediction qualities of infiltration effectiveness [31,32]: Recursive complexity index: (18) where G(xk,yk)[n] represents the predicted effectiveness value of θ, ea, et and Cu in the kth partitioned grid zone coordinated by (xk,yk), n stands for θ, ea, et and Cu, respectively [33].
Computational compactness index: (19) Cluster validity index: Kullback-Leibler divergence index: Fuzzy clustering uniformity index: To identify statistical differences between them, a solid benchmark line represents the standard ratio between the predicted and actual measured results, regulating them into an identical data range. Here, the frames colored by green, red, blue, and brown were utilized to highlight the obtained F1-F6 clusters of θ, ea, et and Cu. Each cluster was illustrated by 120 data points. It could be observed that more data points dispersed around the solid benchmark line means that higher accuracy of the predicted effectiveness values is confirmed. Figure 8a presents data comparisons between the computed and measured recursive complexity index (F1) of θ, ea, et and Cu, obviously showing that the efficiency and reliability of adaptive prediction maintained remarkable correlation with the recursive complexity of predicted results, especially with regard to Wa and Tc. RSAE-ANFIS makes detailed considerations on the inherent capabilities of irrigated soil fields by adapting to external environmental changes, and this novel approach reached a high predictive precision of θ and ea when irrigation cut-off happened. Furthermore, RSAE-ANFIS offers a solid theoretical foundation for the recursive complexity index of θ, ea, and Cu, which means that the predicted infiltration effectiveness values are demonstrated by ea and Cu in a traditional mathematical sense, and they could be affected by deep moisture percolation, moisture content capacity, or irrigation efficiency of the objective soil field. The droplet infiltration quality identified by et reaches optimal values with regard to Fa, Tc, and Pw, which explains their wide applications in keeping the high circulation efficiency of groundwater irrigation, especially when the moisture content variances of soil fields should be fully considered. Figure 8b presents data comparisons between the computed and measured functional inclusion index (F2) of θ, ea, et and Cu, which shows that the F2 index of prediction performance with respect to the probability density of droplet infiltration obviously maintains a close relation with θ and et, and places emphasis on Pw and Fa. The working applicability of effectiveness prediction also highly depends upon the consumption rate of water transpiration and soil moisture content during practical sprinkler irrigation. Since this intelligent prediction system showed an improvement tendency in its assessment predictability or computational compatibility, the calculation efficiency of the RSAE-ANFIS system keeps a positive relation with Wa and Tc, according to the conditional determination of ea and Cu in sprinkler irrigation evaluation. Figure 8c shows data comparisons between the computed and measured computational compactness index (F3) of θ, ea, et and Cu, which gives us an applicable reference criterion related to water use efficiency and the water balance condition when Pw, Fa, and Tc are focused on. Its value variation causes a corresponding fluctuation in the predictive precision of θ and Cu. This assessment index could be employed to identify the distributive rationality of infield infiltration variability, prompting the circulation efficiency and moisture diffusivity of deep moisture percolation. A set of F3 result comparisons demonstrated that, with a moderate reduction in moisture pressure and soil moisture uniformity, ea stays in an invariant value state. Here, droplet infiltration plays an important role in the balanced distributive condition of moisture saturation, the water balance mechanism, and moisture drainage, etc. This assessment index contributes greatly to the stable demonstration of sprinkler irrigation efficiency and moisture infiltration quality on the soil surface. Figure 8d shows data comparisons between the computed and measured cluster validity index (F4) of θ, ea, et and Cu, describing that F4 keeps a close relation with ea and Cu when Wa, Fa, Tc are carefully examined. This assessment index is also affected heavily by soil moisture redistribution and droplet penetration time. As θ indirectly demonstrates the consumption rate of transpiration and infiltration opportunity time, by focusing on moisture diffusivity or infiltration depth, actual sprinkler irrigation presents objective criteria of Pw and Tc to evaluate droplet infiltration effectiveness, from the novel perspectives of soil moisture storage and moisture uniformity. Therefore, the circulation efficiency of groundwater irrigation could be improved by optimizing the infield infiltration rate. The cluster validity index illustrates that the probability density of droplet infiltration provides useful tools to suppress external computational interference originating from noisy measurement signals or error fuzzy rules, and the prediction robustness of RSAE-ANFIS can be improved for moisture redistribution, soil evaporation, or infiltration efficiency. Furthermore, this assessment index also provides a clear demonstration of the infiltration rate and facilitates a reliable evaluative criterion for soil moisture uniformity.  Figure 8e shows data comparisons between the computed and measured Kullback-Leibler divergence index (F5) of θ, ea, et and Cu, which can be frequently used to regulate RSAE-ANFIS computation and monitor the water circulation efficiency in soil fields. This index quantifies the critical influence of Pw and Wa on the moisture content and saturated hydraulic conductivity of a tested soil portion, and it also highly depends on the value range of θ when Cu and et are deliberately adjusted and closely observed. The Kullback-Leibler divergence index demonstrates the invariant properties of Cu and ea when the droplet infiltration conditions impacted by Fa and Tc are faced, so that the sprinkler irrigation stability can be quantified by the measurements of soil moisture content, infiltration depth, or soil evaporation. Figure 8f shows data comparisons between the computed and measured fuzzy clustering uniformity index (F6) of θ, ea, et and Cu, which demonstrates that the calculation precision of RSAE-ANFIS heavily influences ea and Cu, and its operational mechanism could be remarkably affected by the value accuracy of Pw, Wa, and Fa. This index assesses the adaptive prediction of infiltration effectiveness in measurement quality and calculation precision through calibrating the relative proportion of drainage water and soil moisture saturation. It could be observed that this assessment index maintains a relatively stable value during the determination of the circulation efficiency of groundwater. Consequently, it contributes greatly to the precise monitoring of moisture infiltration opportunity time. The fuzzy clustering uniformity index could also be utilized to testify the formula veracity of Wa and Tc based on the cumulative depth of infiltrated water and soil moisture uniformity.
To quantify the applicability and efficiency of the RSAE-ANFIS approach with alternative ones, such frequently used approaches as genetic optimization, simulated annealing-genetic algorithm (SA-GA), Taguchi parameter estimation, artificial neural networksimulated annealing (ANN-SA) prediction, and genetically optimized neural network (GONN) have been employed in the experimental conditions prearranged by Tables 1 and  2 [34][35][36][37][38]. Figure 9a-d showed the value comparisons of θ, ea, et and Cu between the predicted and actual measured results in tests A-K. It could be learned that RSAE-ANFIS (denoted by brown lines) presents θ, ea, et and Cu, keeping closer to actual measured ones (denoted by black lines) than any other approach does. The precision and reliability of the RSAE-ANFIS system could be identified and acknowledged, demonstrating that this proposed approach ensures optimal prediction results of infiltration effectiveness indexes. Table 5 presents a data description about the performance comparisons of the RSAE-AN-FIS approach with others. Considering the current calculation capability of our computation station demonstrated in Section 4.2, the reliability and sensitivity of effectiveness prediction by using these typical approaches were verified in terms of reasoning efficiency and calculation precision. From statistical analysis, it is summarized that higher efficient prediction could be returned by RSAE-ANFIS than that of other approaches. Its average computation accuracy reaches 98.25% in network training and 93.47% in performance testing (error tolerance = ±5%, the same as follows). On the value determination stage of θ, ea, et, and Cu, its standard deviation reaches 0.335% in network training and 0.402% in performance testing, respectively. This intelligent prediction system only uses 1.22 s to complete the entire calculation on the computation platform provided. Other excellent properties of RSAE-ANFIS include average computation storage (1855.4 kb), standard error of prediction (3.66%), and upper and lower error limits (4.27% and 4.52%, respectively), with the frequently used confidence interval of 95% being adopted.   Figure 9 (a) denotes the data variation trend of average soil moisture depth increment (θ, mm), Figure 9 (b) denotes the irrigation infiltration efficiency (ea), Figure 9 (c) denotes the ratio of net infiltrated moisture depth per unit of irrigation time (et, mm/min), and Figure 9 (d) denotes the uniformity coefficient of soil moisture infiltration (Cu, %). Besides, Genetic method denotes genetic optimization, SA-GA stands for the simulated annealing-genetic algorithm, Taguchi is Taguchi parameter estimation method, ANN-SA denotes the artificial neural network-simulated annealing prediction method, and GONN stands for the genetically optimized neural network.

Conclusions
A new adaptive prediction approach used for droplet infiltration effectiveness by using RSAE-ANFIS was proposed to schedule the circulation efficiency of sprinkler irrigation. This research shows that when jetting pressure (Pw) is 255.2 kPa, the impinge angle (Wa) is 42.5°, the water flow rate (Fa) is 0.67 kg/min, and continuous irrigation time (Tc) is 32.4 min, the optimal and stable effectiveness prediction quality could be ensured: the average soil moisture depth increment (θ) is 57.6 mm, irrigation infiltration efficiency (ea) is 62.5%, irrigation turn duration efficiency (et) is 34.5 mm/min, and the uniformity coefficient of soil moisture infiltration (Cu) is 53.6% (error tolerance = ±5%). This research makes the following contributions to sprinkler irrigation: (1) It predicted water droplet infiltration effectiveness in a more efficient way both theoretically and technically; (2) A novel intelligent prediction system named RSAE-ANFIS was designed, with its working mechanism and constructive influence being explained clearly; (3) A set of effectiveness indexes were presented to assess the moisture infiltration mechanism from innovative perspectives; (4) The efficiency, accuracy and quality of sprinkler irrigation can be planned precisely and monitored instantaneously. This novel intelligent prediction system has outstanding predictive capability and possesses much better calculation reliability for water droplet infiltration effectiveness in accuracy and efficiency. Meanwhile, high agreement between the adaptive predicted and actual measured values of infiltration effectiveness could be obtained. Consequently, the superiority and efficiency of RSAE-ANFIS were confirmed, and the instantaneous supervision and mechanism optimization of sprinkler irrigation were facilitated greatly, too. The probability density of droplet infiltration at a given meshed grid zone ∇J The gradient of JRSAE Cf The cumulative irrigation flow depth (mm) Cu The uniformity coefficient of soil moisture infiltration (%) E evaporation ea The irrigation infiltration efficiency et The irrigation turn duration efficiency, or the ratio of net infiltrated moisture depth per unit of irrigation time (mm/min)

Ev
The cumulative water flow depth evaporated at the end of irrigation turn (mm) f The mapping function of the RSAE decoder F1 The recursive complexity index F2 The functional inclusion index F3 The computational compactness index F4 The cluster validity index F5 The Kullback-Leibler divergence index F6 The fuzzy clustering uniformity index Fa The flow rate of water (kg/min) gRSAE The mapping function for the RSAE decoder H (i,j) The inverse of the Hessian matrix i The ith intervals of the update process of WRSAE Id((xk,yk); μ, σ 2 ) The probability density of water droplet infiltration at the investigated area IT The irrigation time (h) j The jth intervals of the update process of WRSAE Mk The number of meshed grids in length directions

N(μ,σ)
Normal distribution of droplet infiltration positions on an objective meshed grid zone Nk The number of meshed grids in width directions Pw The jetting pressure in the nozzle tube (kPa) Q The spray discharge amount (L/h) RH The average air relative humidity (%) sgn The result of the sign function of σr(W•(xk,yk)) SR The solar radiation T The average air temperature (℃) t The irrigation turn duration (min) Tc The continuous irrigation time (min) Wa The impinge angle (º) WRSAE The weight matrix of the Regularized Sparse Autoencoder xk The horizontal coordinates of the specific position where water droplet infiltration happens yk The vertical coordinates of the specific position where water droplet infiltration happens η The step size of the update process θ The average soil moisture depth increment (mm) θantecedent The average antecedent soil moisture depths (mm) θend The average achieved soil moisture depths (mm)

θxk,yk
The soil moisture proportion at one infiltration measurement position (%) λxkyk The regular coefficient of RSAE corresponding to the mesh grid zone μ The mean value of the infiltration position coordinate σ The corresponding distribution variance

σr(z)
The rectified linear unit employed as the activation function in the RSAE decoder σr' The derivative function of rectified linear unit σr