Articles | Volume 26, issue 16
Research article
26 Aug 2022
Research article |  | 26 Aug 2022

A contribution to rainfall simulator design – a concept of moving storm automation

Ravi Kumar Meena, Sumit Sen, Aliva Nanda, Bhargabnanda Dass, and Anurag Mishra

We developed an advanced-design programmable rainfall simulator (RS) to simulate a moving storm rainfall condition. The RS consists of an automated nozzle control system coupled with a pressure regulator mechanism for an operating range of 50 to 180 kPa at a drop height of 2000 mm above the soil flume surface. Additionally, a programmable mobile application was developed to regulate all RS valves. Near natural rainfall conditions were simulated at varying spatial and temporal resolutions in a controlled environment. A soil flume of 2500 mm × 1400 mm × 500 mm was fabricated to conduct different hydrological experiments. The flume was designed to record overland, subsurface, and baseflows simultaneously. This study focused on a detailed analysis of moving storms and their impact on hydrograph characteristics. Experimental results showed a considerable difference in terms of time to peak (tp), peak discharge (Qp), and hydrograph recession for two different storm movement directions (upstream and downstream). Two multiple regression models indicate a statistically significant relationship between the dependent variable (tp or Qp) and the independent variables (i.e. storm movement direction, storm velocity, and bed slope gradient) at a 5 % level of significance. Further, the impact of these moving storm phenomena reduces with the increase in the storm movement velocity.

1 Introduction

Due to high variability in storm pattern, intensity, storm movement velocity, direction, and rainfall drop sizes, it is often challenging to study rainfall characteristics and impacts on overland flow, subsurface flows, baseflows, and soil erosion at a watershed-scale (Singh1998). Thus, rainfall simulation is one such method that is cost-effective and is used to study hydrological processes under controlled rainfall conditions (Nanda et al.2018). Rainfall simulation refers to the process of simulating rain in a confined area for a specific time at a controlled rate. The understanding derived from rainfall simulation experiments is useful in many scientific disciplines, including hydrology, biosystems-engineering, agronomy, and geomorphology (de Lima et al.2003). Rainfall simulator (RS) is an efficient instrument that enables quick data collection and analyses of a wide range of processes and treatment measures based on the variants of simulation configuration as per the study's objectives (Silveira et al.2017).

Rainfall simulation experiments were developed by the United States Soil Conservation Service (SCS) in the 1930s to measure erosion potential and infiltration capacity of the soils. These experiments were used to develop the Universal Soil Loss Equation (USLE), which is employed in various hydrological simulation models (Hall1970). RS has been employed in various environmental studies like mine reclamation, agricultural nutrient transport (Sheridan and Noske2007), and forest hydrology (Croke et al.1999) for the last 40 years. In recent years, many soil erosion studies have been carried out using the portable rainfall simulator (Macedo et al.2021; Mendes et al.2021; Salem and Meselhy2021). Rainfall simulation experiments have evolved with time from mere sprinkler systems to sophisticated computer-based electrical and hydraulic systems (Smith and Schreiber1993; Cai et al.2012).

The experimental observations of rainfall simulation research can be upscaled to the larger hillslope- and catchment-scale. One example is a study conducted by Nasri et al. (2002) for simulating sediment loss of 158 ha catchment using results of 1 m2 experimental plot. Similarly, Nolan et al. (1997) showed erodibility rates of 144 m2 tillage systems using a 1 m2 simulation plot. The use of RS has also started in urban hydrology recently. de Lima and Singh (2002) used single nozzle RS for moving storm analysis to study pollutant build-up, wash-off processes, and urban water quality. These experiments helped understand the slope–soil properties and the effect of surface resistance (i.e. vegetation and micro-topography conditions) on overland flow and infiltration processes. Moreover, RS can also be used to study the flow routing, sediment generation, and transportation at different scales, i.e. plot-scale (Nanda et al.2018) to hillslope-scale (Hall1970).

The spatio-temporal distribution of rainfall influences the overland flow characteristics (de Lima and Singh2002, 2003) in terms of time to peak (tp) and peak discharge (Qp). An important assumption underlying most hydrological experimental methods is that the rainstorm reaches instantaneously over the catchment and remains steady over it (Singh1998). Thus, such hydrological studies ignore the effect of storm movement on catchment runoff response. The exclusion of the storm movement could lead to poor estimation of runoff peaks (Wilson et al.1979).

Most of the time, RS cannot simulate spatio-temporal variability of rainfall like natural rainfall events, and most rainfall simulation studies take rainfall intensity as a constant parameter for a particular area at an instant. Inferences from such studies may not be representative of the actual process underway. To simulate the spatio-temporal variability of rainfall for lab experiments, we designed an advanced programmable RS. The developed simulator was used to conduct different test scenarios in two different slopes and three different velocity conditions for both upstream and downstream storm direction. The following parameters are considered for evaluation of developed moving storm RS:

  1. Peak discharge (Qp) should be higher for the storm moving upstream to downstream than the storm moving in the opposite direction;

  2. Time to peak (tp) of the downstream directional storm should be lower than the upstream directional storm.

Further, the hydrographs generated under two different directions of storm movement (upstream and downstream) on a fully saturated bed condition were analysed for the given experimental configurations; (a) two different slope conditions (2.5 % and 5 %); and (b) three different storm movement velocities (2, 3, and 6 m min−1). A multiple regression model was used to test the statistical significance of the relationship between storm direction and the hydrograph characteristics.

2 Materials and methodology

This section provides a detailed description of the soil flume and moving storm design along with the circuit diagram of the Bluetooth module.

2.1 Structural design

The rainfall simulator (RS) used in this study was designed at the Department of Hydrology, IIT Roorkee, India. The instrument consists of a 3 m × 2 m frame connected with a pipe attached to a header (supporting 11 nozzles) and a pressure gauge (Fig. 1). The frame was supported by four telescopic legs of 6 m each. The rainfall regulating structure connects to a centrifugal pump capable of controlling the water pressure and lifting the water from a feeder tank. The main components of the simulator are the frame, header for nozzle mounting, nozzles, and pumping station. A feeder tank was located near the RS to maintain a sufficient water supply. Water pressure in the system was adjusted by a pressure regulator, and a “shut-off” valve was used to apply back pressure at the outflow end of the simulator system. Another valve was used to facilitate the accurate control of water pressure to the nozzles.

Figure 1(a) Schematic diagram of RS. (b) Spraying System Co. Full jet G-style Spray Nozzle and Nozzle Spray flume. Source: Spraying System Co. Pvt Ltd, Bangalore.

Six full-cone nozzles manufactured by Spraying Systems Co. were used to simulate low to high-intensity rainfall. These nozzles produced a solid cone-shaped spray pattern with a circular impact area (Fig. 1). A uniform spray coverage and distribution over a wide range of flow rates and pressure is possible using these nozzles (B1/88G-SS4.4W). The nozzles also had removable caps and vanes for easy inspection and cleaning.

Figure 2(a) Design of soil flume (not in scale) and (b) filling materials used in soil flume.


2.2 Design of soil flume

A 2.5 mm thick stainless-steel flume of 2500 mm × 1440 mm × 550 mm was fabricated to prepare the soil bed (Fig. 2). A transparent acrylate wall at one vertical side of the soil flume facilitated easy visual observation. A base frame of 500 mm height was designed for stability and to support the jack system. A manually operated worm wheel gear jack setup was installed to change the slope of the flume (0 %–7.5 %). The flume has three sub-partitions to accommodate three different soil types at a single simulation. Outlets for surface flow, sub-surface flow, and deep sub-surface flow gauging were provided at the downstream end of the flume. The surface and the subsurface flow outlets were placed at the height of 500 and 250 mm from the bottom of the flume, respectively. The outlet for the deep sub-surface flow measurement was located at the bottom edge of the flume. Additionally, 10 release/seep slots (5 mm each) were provided at each sub-partition to analyse the change in the piezometric head. These slots can also be used for leachate studies. In the current study, the soil flume was used to evaluate the moving storm rainfall characteristics using the surface runoff data only.

The soil flume was filled with gravel up to 50 mm depth to prevent the washout of the soil. Above the gravel bed, sand was added to a depth of 25 mm, and the remaining 425 mm flume space was filled with sandy loam soil (Fig. 2b). The sand, silt, and clay composition of the soil used was 66 %, 29 %, and 5 %, respectively, measured using a mechanical sieve analyser.

Moving storm design

Detailed descriptions of components used to generate moving storm conditions are presented in Table 1. A set of 11 nozzles were used for simulating the moving storm condition. Electrically operated flow control valves were used to control these nozzles through an Arduino Mega (AM) microcontroller board. A nozzle control system was incorporated using three components: servo-operated valve, AM microcontroller, and Bluetooth module (BM). This system serves two purposes; communication with the user interface in the handled device and control of the opening and closing of nozzles. The detailed operational flowchart of the moving storm system is shown in Fig. 3.

Table 1Specification of components used for moving storm rainfall simulator.

Download Print Version | Download XLSX

Figure 3The operational structure of moving storm rainfall simulator.


The nozzles were grouped into four clusters for this experiment to ease the flow regulation operation (Fig. 4). Three clusters (NC1, NC2, and NC3) consisted of three nozzles, and one cluster (NC4) consisted of two nozzles. These clusters were activated and deactivated with a specific time gap to simulate the moving storm over the plot area. If needed, full control can be given to each nozzle to regulate individually. To obtain variable rainfall intensities, a servo motor-operated flow control valve was inserted into the pipe openings just before the nozzle. Further, an android mobile application was developed for regulating the valves through the BM.

Figure 4Nozzle distribution of rainfall simulator over the soil flume.


The pressure regulating system (PRS) comprises of a motorised globe valve, pressure transmitter, and proportional integral derivative controller (PID controller). The PRS was designed to maintain constant pressure throughout the simulation to achieve a constant rainfall intensity regardless of the opening and closing of the various number of nozzles.

2.3 Nozzle control system

2.3.1 Servo-operated valve

A stop cock valve was used to develop a servo-operated valve due to its low operational torque requirement. Servo motor of torque 0.98 Nm which is easy to control and has fairly high accuracy was used to control the valve. An aluminium frame was fabricated using a 2.5 mm aluminium sheet to hold the servo motor and the stop cock valve together (Fig. 5a).

Figure 5Components used for nozzle control system.


Servo motor MG995 is a high-speed digital motor with a rotation angle of 90 in both directions, enabling a 180 reach (Fig. 5b). Pulse-width modulation signals were used for the operational control of the servo motor to increase process speed and efficiency. It was equipped with an internal circuit which gave high torque and better stability.

2.3.2 Arduino mega

Arduino Mega (AM) is a microcontroller board with 54 digital I/O pins, four hardware ports, 16 analogue ports and a 16 MHz crystal oscillator complemented by an In-Circuit Serial Programming (ICSP) header, a power jack, and a reboot button. It can be powered through both USB as well a DC supply of 7–12 V 1 A.

2.3.3 Bluetooth module

A HC-05 Bluetooth module (BM) can be enabled both as a master and as a slave (Fig. 5c). The master setting enabled auto-communication between the two Bluetooth devices whereas the slave set could only accept the incoming connection from the master Bluetooth device. It had a 3 Mbps data transmission speed with a 2.4 GHz transmitter and receiver. It is comprised of six pins; Vcc  for power supply, Gnd for negative, Tx  for transmission, Rx  for receiving, a key to switch between master and slave, and an LED to display its operational activity. In this experiment, the Bluetooth module was used as a slave set (default configuration).

2.4 Pressure regulating system (PRS)

2.4.1 Motorised globe valve

A motorised globe valve was used to control the bypass flow and to maintain a constant pressure in the main line which worked through a pressure feedback circuit from the main line (Fig. 6a). The header of the RS was kept closed so the pressure of the simulator line can be maintained by controlling the bypass flow (return flow). A two-way globe valve of metal to metal seating with a pulse-based modulation of 4–20 mA actuator was used. It had an intelligent circuit that sensed hindrance in valve movements. An AC sensor was used for circuit protection that can shut down the valve during an overload condition.

Figure 6Components used for pressure-regulating system.


2.4.2 Pressure transmitter

A Mass PT11 pressure transmitter has a very compact design with stainless steel construction. It is highly stable against shock and vibration, and also has features such as reverse polarity, limit protection, and high accuracy. This pressure sensor was installed to check the main line pressure. PID controller can sense the change in pressure and can act accordingly to operate bypass and to maintain a constant pressure in the main line for a uniform rainfall intensity (Fig. 6b).

2.4.3 Proportional integral derivative (PID) controller

Selec PID500 is a controller which is employed widely in industrial process controls. PID controller is a control loop feedback system and is used to operate a motorised globe valve on the basis of an input signal from a pressure sensor (Fig. 6c). Whenever there is a pressure offset from the set value, it sends a signal to the motorised valve to re-attain the set value. The I/O signal from the pressure sensor and PID, respectively, ranged between 4–20 mA. The controller had a compact square housing with panel mounting facility in its enclosure powered by a 240 V AC supply. PID controller was used to control the bypass flow by operating a motorised valve to maintain constant pressure in the main line against any pressure drop generated due to the moving storm simulation.

2.5 Circuit design

The circuit design comprised of AM, BM, and 11 servo motor-operated valves. The detailed design of the circuit diagram is given in Fig. 7. The four connections, Rx, Tx, Vcc, and Gnd in the BM were connected to pin 11, pin 10, 3.3 V, and Gnd recipient pins in the AM, respectively. The 11 servos were clustered into four sets of 3, 3, 3, and 2 each in every cluster. The signal pins of these clusters of servo motors were connected to the digital signal pins of AM numbered as pin 3, pin 5, pin 6, and pin 9. Each group had a power supply of 5 V–2 A, and these were grounded to AM which had a power supply of 12 V–1 A.

Figure 7Circuit diagram of 11 servo motors operating through Arduino Mega and HC – 05 Bluetooth module.


The AM was coded so that the servo motor could be regulated to control all the motors simultaneously with any android-based phone through a Bluetooth application. Two basic libraries used in this code were “SoftwareSerial.h” and “Servo.h.”. The software used to write the code is Arduino Editor, and the free coding software is available online (, last access: 22 August 2022). The android OS application was designed to operate all of the 11 nozzles simultaneously with one-touch Bluetooth connectivity and four slider bars to control four groups of servo-operated valves at any value ranging from 0 % to 100 %. This application was developed using the “MIT app developer” software.

2.6 Simulation uniformity assessment

We used the following equation to calculate Christiansen's Uniformity Coefficient (UC):

(1) UC = 100 1 - i = 1 n X i - μ i = 1 n X i

where μ is the average of all the measurements, |Xiμ| is the sum of the individual deviations from the mean, and n is the number of measurements. The UC was measured using 66 beakers kept in a square array 250 mm apart beneath the simulator covering the plot area of 2500 mm × 1440 mm.

2.7 Design of the experiments

The experiments consisted of 12 scenarios (Table 2) with three replications. Experiments were performed under fully saturated soil bed conditions to reduce the variability among scenarios and replications.

Table 2Design structure of experiments.

Download Print Version | Download XLSX

The simulation results were analysed using a multiple regression model considering one categorical variable (storm direction) and two numerical variables (velocity and bed slope) as independent variables, and time to peak (tp) or peak discharge (Qp) as the dependent variable. An indicator variable was used to include the direction of the storm in the multiple regression model (downstream = 0, upstream = 1). An indicator variable allows interpretation of the regression coefficients for storm direction as an additive effect on the hydrograph characteristics. By choosing the indicator variable for an upstream storm as 1, the regression coefficient of storm direction indicated the difference in the mean response of time to peak (tp) or peak discharge (Qp) of the upstream directional storm. A positive coefficient implies a positive effect on upstream directional storm on the hydrograph characteristics compared to downstream directional storm. The null hypothesis of regression coefficient for each of the independent variables as equal to zero was tested against alternative hypothesis of significant effect of the independent variable on storm hydrograph characteristics at 5 % level of significance.

3 Results and discussion

After completing the design of the moving storm rainfall simulator (RS) (Fig. 8), we checked the feasibility of the RS for generating the moving storm events. Before stepping into the different rainfall scenarios analysis, the rainfall distribution over the plot was analysed (Fig. 9). It can be elucidated from the rainfall distribution graph that 70 % of the plot area receives a uniform amount of rainfall, i.e. 30 to 36 mm. The lower rainfall (15 to 22 mm) amount was recorded at the plot edges. However, the overall UC was found to be 84.2 %. de Lima and Singh (2003) also conducted their rainfall simulator experiment with an average UC of 88 %. Similarly, Macedo et al. (2021) and Salem and Meselhy (2021) conducted rainfall simulation experiments for studying soil erosion at a UC of 75 % and 89 %–94 %, respectively. Further, Mendes et al. (2021) also carried out simulation tests for studying geotechnical and hydrological phenomena with a UC of 75 %. All the rainfall simulation experiments were carried at a maximum 3 min intensity of 40 mm h−1. The maximum 5 min rainfall intensity in our Aglar experimental watershed ranges between 12 to 109 mm h−1 (Nanda et al.2019). Thus, the selected intensity is reasonable while considering the natural rainfall phenomena.

Figure 8Designed rainfall simulator for moving storm experiments.


Figure 9Rainfall distribution graph.


3.1 Results of experimentation conducted at 2.5 % slope

The results recorded by the rainfall simulator (RS) of the moving storm clearly exhibit the effect of storm direction, velocity, and slope on the overland flow hydrographs (Fig. 10). For the velocity of 6 m min−1, non-detectable runoff was generated and, therefore, it was not included in the result.

Figure 10Hydrograph for velocity 2 m min−1 (a) and 3 m min−1 (b) at 2.5 % slope.


The time to peak (tp) of the hydrographs generated by the upstream to the downstream storm was less than the downstream to the upstream storm irrespective of their velocities (Fig. 10). When the storm moves towards the outlet (i.e. upstream to downstream), the overland flow was initiated at the upstream point and it moves downstream along with the storm direction. Thus, the collective overland flow reached the outlet simultaneously with the storm which resulted in a sharp peak at the outlet of the soil flume. However, peak discharge (Qp) of the upward moving storm was possible only when the entire catchment contributed, and that occurred when the storm reached the most upstream point.

The hydrograph produced by the downstream to upstream storm movement displays a 110 s longer recession time compared to the opposite storm movement especially during low storm velocity (2 m min−1). The runoff and storm are in opposite direction and, thus, the upstream runoff water took longer time to reach at the outlet, resulting in longer recession time. However, a negligible change in the recession time was observed between the storm directions during the 3 m min−1 storm velocity (Fig. 10b). de Lima and Singh (2003) also observed similar observations with an increase in storm velocity. When the storm was moving in the upstream direction, 2 m min−1 storm velocity showed a 150 s longer recession time than 3 m min−1 velocity. But, no such change in recession time was observed between the velocity conditions for downstream directional storms.

3.2 Results of experimentation conducted at 5 % slope

At 5 % slope condition, rainfall simulation experiments were performed with three different velocities; 2, 3, and 6 m min−1 (Fig. 11). The storm movement directions were the same as the previous experiments, i.e. upstream and downstream. It can be clearly illustrated from Fig. 11 that the recession characteristics, time to peak (tp), and peak discharge (Qp) followed the same trend as the hydrographs generated at 2.5 % slope. An interesting observation was noticed during the testing of 6 m min−1 storm velocity, i.e. the recession curve and Qp of both hydrographs completely matched with each other during the upstream and downstream directional storm movement (Fig. 11c). Only tp varied slightly in these two hydrographs of 6 m min−1 velocity storm.

Figure 11Hydrograph for velocity 2 m min−1 (a), 3 m min−1 (b), and 6 m min−1 (c) at 5 % slope.


A detailed description of storm characteristics during different test scenarios is presented in Table 3. It is observed that as the slope increased, the tp value decreased when the storm was moving to the upstream direction. While moving towards the downstream, the Qp value of 5 % slope was 53.7 % and 43.3 % higher than 2.5 % slope condition at a velocity of 2 and 3 m min−1, respectively. Similarly, for upstream directional storm at a velocity of 2 and 3 m min−1, the Qp of 5 % slope was 59 % and 42.8 % higher than the 2.5 % bed slope condition, respectively. These analyses concluded that the Qp value increases significantly with an increase in soil flume slope and decreases with increased storm velocity.

Table 3Detailed moving storm characteristics of different scenarios.

Download Print Version | Download XLSX

Further, observations show that the discharge volume followed the same relationship pattern as Qp with slope and storm velocity. The discharge volume of the 5 % slope condition was 67.38 % and 82.00 % higher than the 2.5 % slope for the downstream and upstream directional storm, respectively. Similarly, for 3 m min−1 velocity, the increase in discharge volume due to slope increment was 58.56 % and 55.56 % for downstream and upstream storm directions, respectively.

From the above sets of experiments, it can be concluded that rainfall moving upstream results in slower hydrograph rise time, lower Qp and longer recession time (base time) compared to rainfall moving downstream. de Lima and Singh (2003) found that the high-velocity storm results in a smaller volume of runoff. This phenomenon was also observed in our study for both the slope conditions. For the 2.5 % slope condition, storm velocity of 6 m min−1 barely generated any runoff at the plot outlet; thus, it did not create any runoff hydrographs. Moreover, de Lima and Singh (2003) observed that high storm velocity did not result in significant changes in hydrographs during upstream and downstream storm movement because of the surface tension force. Similar observations are also noted in the present study during 6 m min−1 storm velocity (Fig. 11c). The storm movement was so fast that the whole plot could never contribute to runoff generation at the outlet at an instant.

The difference between Qp of the upstream and downstream directional storm is the function of storm velocity. Further, de Lima et al. (2011) discussed the influence of slope on moving storm runoff hydrograph, i.e. how the runoff volume increases with an increase in slope angle. The effect of storm direction on Qp is also discussed by de Lima and Singh (2003), Isidoro et al. (2012), and Seo and Schmidt (2013). They concluded that the downstream storm produces higher Qp than the upstream directional storm. Similar results were also obtained in the current study (Table 3).

To further test the significance of the effect of the experimental variables on the observed differences in time to peak (tp) and peak discharge (Qp), multiple regression analysis was performed. Regression analysis with an indicator variable for storm direction results in two regression model equations for each hydrograph characteristic modelled (Eq. 2).

(2) t p or Q p = β 0 + β 1 + β 2 × velocity + β 3 × slope , when direction = upstream = 1 β 0 + β 2 × velocity + β 3 × slope , when direction = downstream = 0

The model fit details for the regression models are shown in Table 4. Interpreting the value of β1 in the models for a storm moving upstream, tp is higher (positive) by 10±4.776 s and Qp is lower (negative) by 0.00256±0.00086 L s−1. These values are significant at 5 % level of significance, thus, providing a satisfactory hydrological verification for the moving storm rainfall simulator. The multiple regression models for tp and Qp explain 89 % and 94 % of the variability in tp and Qp observations (as given by model R2 values), thus, indicating the sufficiency of the model for explaining the relationship between experimental variables and hydrograph characteristics. The sign of the regression coefficients of storm velocity and bed slope agree with the empirical conclusions that Qp increases with an increase in bed-slope and decrease in velocity. In contrast, both bed slope and velocity have a negative effect on tp.

Table 4Regression coefficients, standard error, the t statistic, and associated p value for multiple regression models.

Download Print Version | Download XLSX

The main objective of this study is to develop a handy, multifunctional, and advanced rainfall simulator to study the moving storm rainfall pattern. After development, 36 experiments characterised into 12 different scenarios were conducted to check the instrument's feasibility. From the above discussion, it is evident that the tool could generate the moving storm condition satisfactorily. However, this study was limited to evaluating the impact of storm movement on generated hydrograph under a single storm pattern with two different slope conditions, two different storm directions, and three different moving storm velocities. This advanced-design RS can be considered to be used in the future to analyse the impact of storm movement over soil erosion and nutrient transport. Further, the designed flume can be used for subsurface flow, base flow, and leachate studies.

4 Conclusions

A multi-nozzle programmable RS was designed to simulate moving storm rainfall. Twelve different test scenarios were examined in two different slopes (2.5 % and 5 %) conditions and three different moving storm velocities (2, 3, and 6 m min−1). In these experiments, storm movement was considered along the slope (downstream) and against the slope (upstream) of the basin, keeping rainfall intensity and soil saturation constant. The results indicate a statistically significant influence of spatial and temporal distributions of rainfall on hydrograph and its characteristics.

The following conclusions were drawn from the study:

  1. The hydrograph generated from the downstream moving storms yielded a higher peak with a sharp rise and short recession limb;

  2. The upstream moving storm produced a runoff hydrograph of lower peak and a prolonged, gradually decreasing recession limb;

  3. With the increase in the storm movement velocity, the impact of moving storm direction in terms of discharge peak and time to peak became negligible in either direction of storm movement;

  4. An increase in the basin slope reduced the impact of moving storm direction on overland flow.

Code availability

The software code is available on Zenodo (Meena2022).

Data availability

Experimental data can be obtained from Zenodo (Meena2022).

Author contributions

The main contributions from each co-author are as follows. RKM contributed to the methodology, software, validation, and draft preparation. SS contributed to the conceptualisation, supervision, and funding acquisition. AN and BD contributed to the visualisation, writing, and editing of the paper. AM contributed to the editing and review of the paper.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement

This article is part of the special issue “Experiments in Hydrology and Hydraulics”. It is not associated with a conference.


The authors would like to acknowledge the Department of Hydrology, Indian Institute of Technology Roorkee, Roorkee, India, for providing support and resources for this experimental study. The authors would also like to acknowledge Anbanandam Ramesh and Daniel Denzil for statistical analysis interpretation and Pallapu Prasanthi for proofreading the initial draft of the article.

Financial support

The research has been supported by the Department of Hydrology, Indian Institute of Technology Roorkee, India.

Review statement

This paper was edited by Jorge Isidoro and reviewed by Jr-Chuan Huang and one anonymous referee.


Cai, J., Li, P., and Wang, P.: An approach to rainfall simulator automation and performance evaluation, in: Proceedings of the World Congress on Intelligent Control and Automation (WCICA), 6–8 July 2012, Beijing, China, 3428–3433,, 2012. a

Croke, J., Hairsine, P., and Fogarty, P.: Sediment transport, redistribution and storage on logged forest hillslopes in south-eastern Australia, Hydrol. Process., 13, 2705–2720,<2705::AID-HYP843>3.0.CO;2-Y, 1999. a

de Lima, J. L. and Singh, V. P.: The influence of the pattern of moving rainstorms on overland flow, Adv. Water Resour. 25, 817–828,, 2002. a, b

de Lima, J. L. and Singh, V. P.: Laboratory experiments on the influence of storm movement on overland flow, Physics Chem. Earth, 28, 277–282,, 2003. a, b, c, d, e, f

de Lima, J. L., Singh, V. P., and Lima, M. I. P. D.: The influence of storm movement on water erosion: Storm direction and velocity effects, Catena, 52, 39–56,, 2003. a

de Lima, J. L. M. P., Singh, V. P., Isidoro, J. M. G. P., and de Lima, M. I. P.: Incorporating the Effect of Moving Storms into Hillslope Hydrology: Results from a Multiple-Slope Soil Flume, in: World Environmental and Water Resources Congress 2011, World Environmental and Water Resources Congress 2011, Palm Springs, California, USA, 1398–1407,, 2011. a

Hall, M. J.: A Critique of Methods of Simulating Rainfall, Water Resour. Res., 6, 1104–1114,, 1970. a, b

Isidoro, J. M., de Lima, J. L., and Leandro, J.: Influence of wind-driven rain on the rainfall-runoff process for urban areas: Scale model of high-rise buildings, Urban Water J., 9, 199–210,, 2012. a

Macedo, P. M. S., Pinto, M. F., Sobrinho, T. A., Schultz, N., Coutinho, T. A. R., and de Carvalho, D. F.: A modified portable rainfall simulator for soil erosion assessment under different rainfall patterns, J. Hydrol., 596, 126052,, 2021. a, b

Meena, R. K. and Sen, S.: Innovatory rainfall simulator design – A concept of moving storm automation DATA [Data set]. In Hydrology and Earth System Sciences (d.1.1), Zenodo [code and data set],, 2022. a, b

Mendes, T. A., Pereira, S. A. D. S., Rebolledo, J. F. R., Gitirana, G. D. F. N., Melo, M. T. D. S., and Luz, M. P. D.: Development of a rainfall and runoff simulator for performing hydrological and geotechnical tests, Sustainability, 13, 3060,, 2021.  a, b

Nanda, A., Sen, S., Jirwan, V., Sharma, A., and Kumar, V.: Understanding plot-scale hydrology of Lesser Himalayan watershed – A field study and HYDRUS-2D modelling approach, Hydrol. Process., 32, 1254–1266,, 2018. a, b

Nanda, A., Sen, S., and McNamara, J  P.: How spatiotemporal variation of soil moisture can explain hydrological connectivity of infiltration-excess dominated hillslope: Observations from lesser Himalayan landscape, J. Hydrol., 579, 124146,, 2019. a

Nasri, S., Pépin, Y., Hamed, Y., Berndtsson, R., Zante, P., Asseline, J., Albergel, J., El-Niazy, M., and Balah, M.: Comparison between rainfall simulator erosion and observed reservoir sedimentation in an erosion-sensitive semiarid catchment, Catena, 50, 1–16,, 2002. a

Nolan, S. C., Vliet, L. J. P. V., Goddard, T. W., and Flesch, T. K.: Estimating storm erosion with a rainfall simulator, Can. J. Soil Sci., 77,, 1997. a

Salem, H. M. and Meselhy, A. A.: A portable rainfall simulator to evaluate the factors affecting soil erosion in the northwestern coastal zone of Egypt, Nat. Hazards, 105, 2937–2955, 2021. a, b

Seo, Y. and Schmidt, A. R.: Network configuration and hydrograph sensitivity to storm kinematics, Water Resour. Res., 49, 1812–1827,, 2013. a

Sheridan, G. J. and Noske, P. J.: A quantitative study of sediment delivery and stream pollution from different forest road types, Hydrol. Process., 21, 387–398,, 2007. a

Silveira, A., Isidoro, J. M. G. P., de Deus, F. P., dos Reis, S. S., Silva, A. M. D., Gonçalves, F. A., Menezes, P. H. B. J., and Tiezzi, R. D. O.: Enhancing the spatial rainfall uniformity of pressurized nozzle simulators, Manage. Environ. Qual., 28, 17–31,, 2017. a

Singh, V. P.: Effect of the direction of storm movement on planar flow, Hydrol. Process., 12, 147–170,<147::AID-HYP568>3.0.CO;2-K, 1998. a, b

Smith, S. A. and Schreiber, J. D.: Design of a programmable rainfall simulator, in: 1993 25th Southeastern Symposium on System Theory, SSST 1993, 7–9 March 1993, Tuscaloosa, AL, USA, 215–217,, 1993. a

Wilson, C. B., Valdes, J. B., and Rodriguez‐Iturbe, I.: On the influence of the spatial distribution of rainfall on storm runoff, Water Resour. Res., 15, 321–328,, 1979. a

Short summary
We developed a mobile operated programmable rainfall simulator (RS) to simulate the near-natural moving storm rainfall condition and to study its impact on runoff, soil erosion, and nutrient transport. The designed RS can be used for variable velocity and slope conditions along with the three different soil types at a time. Moreover, the soil flume of the RS is associated with the surface, subsurface, and base flow components.