Hydrological Analysis of a Dyke Pumping Station for the Purpose of Improving Its Functioning Conditions

: This paper presents a methodology of calculating the water transfer capacity of a dyke pumping station in ﬂood wave conditions in order to improve its functioning, especially with regards to the safety of the areas being drained. The exemplary analysis was carried out for a pumping station situated on a small right-bank tributary of the Odra River in the southwest part of Poland, which, due to the inadequate capacity of its pumps, extensively ﬂooded the surrounding areas in May and June 2010. Hydrological analyses were conducted in order to determine the rate of the designed and control ﬂows using a spatial regression equation, and as a comparison, the rainfall-runoff method was also used. The corresponding ﬂood-wave hydrographs were also determined, which included total precipitation using the German Association For Water Resources and Land Improvement (DVWK) method, effective precipitation using the Natural Resource Conservation Service curve number (NRCS-CN) method, as well as hypothetical waves using the instantaneous unit hydrograph (IUH) method. Flood-wave routing was carried out and alternative solutions for both the output of the required pumps and the retarding reservoir capacity were highlighted on this basis. The paper presents the possibility of a correct pump capacity selection, and in turn, the size of the pumping station retarding reservoir that results from this selection. This will enable pumping station exploitation costs or maintenance costs of the retarding reservoir to be considerably


Introduction
Dyke pumping stations on the tributaries of main rivers are hydraulic structures, the operation of which is vital for the flood protection of areas over dykes when the main river is flooded. Unfortunately, not enough attention is paid to the necessary dimensioning of such structures, which results in design errors such as the insufficient capacity of installed pumps, an inadequate capacity of a retarding reservoir, or even the lack of such a reservoir. The design of dyke pumping stations involves various difficulties when the catchment area is uncontrolled. In this case, the first step is to calculate flows using the probability of exceedance, known in Poland as designed and control flows. The genetic rainfall formula, the spatial regression equation, rainfall-runoff models, and other regional methods could be applied in Polish conditions to determine such flows, as well as the shape and volume of designated flood waves. Szolgay et al. [1] analyzed the relationship between flood peaks and the corresponding flood event volumes, which were modeled by empirical and theoretical copulas in a regional context. They tested the reliability of different formulas of synoptic floods, flash floods, or snowmelt floods and found out that different flood types should be treated separately. Grimaldi and Petroselli [2] proposed an approach that uses input information similar to those required by the Rational Formula is that the ratio of actual retention to potential retention is the same as the ratio of actual runoff to potential runoff.
A simple formula for calculating wave hydrographs, based on precipitation, was first proposed by Nash [19,20] and was based on the linear reservoir model to develop an instantaneous unit hydrograph (IUH) within a river basin. The linear reservoirs assumed in Nash's model are imaginary reservoirs in which storage is directly proportional to the outflow from them [21]. This model has been used by many authors (e.g., [21][22][23]).
This paper presents a methodology that was proposed by the authors, which is adequate for dimensioning dyke pumping stations for an uncontrolled river basin. In the presented paper, the river catchment is sufficiently described with all the required physical parameters and also has average annual hourly rainfall observations. This fulfills the known limitations of the NRCS-CN model. The spatial regression equation was used in order to make a comparison with the results calculated using the rainfall-runoff model. This allowed maximum flows with a definite probability of exceedance to be estimated for an uncontrolled river catchment. As an illustration, the proposed methodology was applied to the case of the Ciechowice land dyke pumping station, which is situated on an outlet of the Łęgoń watercourse to the Odra River. During a flood in 2010, extensive areas were flooded including the pumping station building and other technical facilities, and heavy losses were incurred as a consequence.

Study Area
The Łęgoń watercourse is a right-bank tributary of the Odra River ( Figure 1). The sources of this watercourse are situated at the edge of the Moravian Gate Arboretum at an elevation of 184.25 m a.s.l. The entire watercourse flows through Racibórz County in the Silesian Province. The watercourse catchment is characterized by a fairly dense drainage network with numerous drainage ditches. The Łęgoń watercourse catchment is a lowland used for agriculture and is characterized by small gradients.
The catchment is mainly built from river alluvia. Fluvial soils occupy its eastern part, while in the rest of its area, there are soils formed from old accumulation terrace (weakly argillaceous and argillaceous) sands, clay overlying sands and lightweight boulder clays, and soils formed from loess deposits.
The average annual precipitation for the multiannual period 2003-2012 in the catchment amounted to 614.3 mm. The average annual evaporation for this period was 564 mm. The annual climatic water balance (CWB) amounted to 50.3 mm. The average annual temperature on the basis of the climate atlas of Poland is 8.5 • C. The conditions prevailing in the Łęgoń river catchment, especially the small river bed slope, tend to favour the formation of high overbank flows.
When flume stage conditions occur in the Odra river, waters of the Łęgoń river gravitationally flow through two 1.20 m × 1.40 m dyke culverts. At high water in the Odra river dyke, culverts are closed by manoeuvrable flaps and water is then pumped with the Ciechowice dyke pumping station from the retarding reservoir to the Odra River. According to Polish technical regulations, the Ciechowice dyke pumping station should be built with regards to the designed and control flood flows, which should be estimated with the probability of exceedance of 1% and 0.3%, respectively. The reinforced concrete main building consists of underground intake wells and an aboveground pump house with four two-stage vertical-axis vortex impeller pumps, each of which has a capacity of 0.9 m 3 ·s −1 (3200 m 3 ·h −1 ). Upstream of the pumping station, there is a retarding reservoir in the form of a widened (for a length of about 100 m) Łęgoń watercourse channel. The retarding reservoir is about 15 m wide at its bottom.

Methodologies of Flows Estimation
As the Łęgoń watercourse is an uncontrolled river and there are no water stages and flow measurements, a rainfall-runoff model and a spatial regression equation were used to determine the probable maximum flows. Regional interdependences between unit flows and maximum flows, with a probability of exceedance p = 1%, were developed to verify the calculation results. The values yielded by the rainfall-runoff model are given as the ultimate results of the design and control flood flow calculations due to the fact that the NRCS-CN method is the most commonly used method for estimating peak discharge in a river catchment.
Before estimating the peak discharge in the special watershed, it is better to test the NRCS-CN method with some accurate models to precisely estimate the impact of individual parameters. The sensitivity of the models increases with an increase in the return period. Therefore, calibration of this method for high return periods is necessary [15]. The spatial regression equation and regional interdependences between unit flows and maximum flows were used for such a calibration.
For maximum flows, wave hydrographs must be defined, i.e., shape and volume with the use of a mathematical model of precipitation and drainage in order to analyze the flood wave routing through the retarding reservoir for both the assumed number and efficiency of pumps and also the volume of the reservoir. An optimum solution involves the determination of the maximum retarding reservoir capacity associated with land topography with regards to lower reservoir maintenance costs when compared to a high capacity of a pump station.

Design Flows-The Spatial Regression Equation
The spatial regression equation allows the maximum flow to be calculated with a definite probability of exceedance in hydrological uncontrolled river catchments with an area of 50 ÷ 2000 km 2 . This equation, which is recommended by the Polish Hydrologists Association, is written in the following form [24]: where Qmax p% is the annual maximum flow with a probability of exceedance of p = 1% in m 3 ·s −1 and λp is the quantile for the nondimensional regional curves of the maximal flows, which equals λp = 1.224 for p = 0.3%. The maximal flow with a probability of exceedance p = 1% can be calculated using the following equation [23]:

Methodologies of Flows Estimation
As the Łęgoń watercourse is an uncontrolled river and there are no water stages and flow measurements, a rainfall-runoff model and a spatial regression equation were used to determine the probable maximum flows. Regional interdependences between unit flows and maximum flows, with a probability of exceedance p = 1%, were developed to verify the calculation results. The values yielded by the rainfall-runoff model are given as the ultimate results of the design and control flood flow calculations due to the fact that the NRCS-CN method is the most commonly used method for estimating peak discharge in a river catchment.
Before estimating the peak discharge in the special watershed, it is better to test the NRCS-CN method with some accurate models to precisely estimate the impact of individual parameters. The sensitivity of the models increases with an increase in the return period. Therefore, calibration of this method for high return periods is necessary [15]. The spatial regression equation and regional interdependences between unit flows and maximum flows were used for such a calibration.
For maximum flows, wave hydrographs must be defined, i.e., shape and volume with the use of a mathematical model of precipitation and drainage in order to analyze the flood wave routing through the retarding reservoir for both the assumed number and efficiency of pumps and also the volume of the reservoir. An optimum solution involves the determination of the maximum retarding reservoir capacity associated with land topography with regards to lower reservoir maintenance costs when compared to a high capacity of a pump station.

Design Flows-The Spatial Regression Equation
The spatial regression equation allows the maximum flow to be calculated with a definite probability of exceedance in hydrological uncontrolled river catchments with an area of 50 ÷ 2000 km 2 . This equation, which is recommended by the Polish Hydrologists Association, is written in the following form [24]: where Q max p% is the annual maximum flow with a probability of exceedance of p = 1% in m 3 ·s −1 and λ p is the quantile for the nondimensional regional curves of the maximal flows, which equals λ p = 1.224 for p = 0.3%.
The maximal flow with a probability of exceedance p = 1% can be calculated using the following equation [23]:  (2) where: Q max p=1% is the maximum flow with a probability of exceedance equal to 1% in m 3 ·s −1 ; α area is the regional parameter of the formula, which is taken on the basis of the investigated area location in Poland and is equal to α area = 2.992 × 10 −3 ; A is the river catchment area that amounts to A = 49.40 km 2 ; H 1 is the maximum daily rainfall for a probability of exceedance of p = 1%, which is taken from the meteorological station in Raciborz and is equal to H 1 = 87.3 mm; and φ is the runoff coefficient for peak flows, which is equal to φ = 0.55 and qualified on the basis of "The map of Polish soils". For a catchment with several soil groups with different values of runoff coefficient φ, a medium weight should be calculated for the whole catchment.
where A i is an area covered with soil of a given group in km 2 and I r is the longitudinal river bed slope that is equal to I r = 0.73‰.
where: W g is the elevation of the watershed in the point crossing the axis of a dry valley of the longest watercourse in m a.s.l.; W p is the elevation of the calculated cross-section closing the catchment in m a.s.l.; L is the length of the longest watercourse in the catchment in km; l is the length of a dry valley in the lengthened part of the longest watercourse of the catchment in km; and ψ is the average slope of the river catchment, which is equal to ψ = 17.0‰.
where: W max is the maximal catchment elevation in m a.s.l.; A is the catchment area in km 2 ; and lake is the lake index of the river catchment.
where A lake i is an area of the lake i catchment that is equal to A lake = 2.5 km 2 and m is the number of lake catchments, which is equal to m = 1. B is the swamp index of the river catchment, which is equal to zero.
where A B I is an area of i-succeeding swamp or peat land region in km 2 and k is the number of swamp regions. Details on how to determine the particular variables in Formula (2) can be found in the extensive literature on the subject (e.g., in [24]).
The following physical-geographic parameters of the Łęgoń catchment were assumed in the conducted calculations: a catchment area of 49.40 km 2 , the maximum catchment length of 11.57 km, an average catchment gradient of 1.7%, and also a forestation ratio of 15%. Therefore, the maximum flows to the Ciechowice pumping station, with the specified probability of exceedance calculated using the spatial regression Formula (2), amounted to:

The NRCS-CN Method
Precipitation with a probability of occurrence of p = 1% and 0.3% during T = 24 h was assumed in the calculations. The precipitation distribution over time for the Łęgoń catchment was determined using the method proposed by German Association For Water Resources and Land Improvement DVWK [25]. According to the recommendations of this method, rain with a maximum intensity in the middle of its duration should be assumed as the reliable precipitation intensity distribution. The following precipitation distribution was assumed in the calculations and expressed as a percentage of the total precipitation: The total precipitation hyetograph can also be expressed by a beta distribution, the density function of which is described by Formula (8): where: P(t) is the total precipitation in mm; t is precipitation duration in hours; α and β are distribution parameters of beta function B(α,β); and T is the total time of precipitation in hours. When creating a model for the division of the daily precipitation total into 1 h calculation intervals, such beta distribution parameters were selected in order to maintain agreement with the DVWK method assumptions [25], i.e., α = 4.5, β = 6.1. The values of the maximum daily precipitation totals with occurrence probability of p = 1% and p = 0.3% are equal to 87.3 and 101.6 mm accordingly and were obtained from the Racibórz rainfall gauging station.
Daily rainfalls with a probability of exceedance of 0.3% and 1%, which were estimated on the basis of measured values from the period of 1998-2013, were set in a series and their homogeneity was verified with the use of the Mann-Kendall test. For these series, the parameters of log-norm distribution were estimated using the maximal reliability method and their correctness was verified with the Kolmogorow test. Using this distribution, the values of the required quantiles were obtained. Uncertainty was calculated as the upper limit of one-sided confidence interval β = 84% of P β max,p . In the performed calculations, the value of uncertainty for probability of 1% and 0.3% were equal to 96.5 mm and 116.1 mm, which is 10.5% and 14.3%, respectively. The calculations were performed by the Polish Institute of Meteorology and Water Management.
The effective precipitation was determined using the NRCS-CN method [9][10][11]. According to this method, the following factors have a bearing on effective precipitation: the land development, the kind of soils, the kind of vegetation cover, and the catchment wetness (wetness degree II was assumed). The totality of the factors was expressed by a dimensionless parameter CN, assuming values that ranged from 0 to 100. The NRCS-CN method is based on the hypothesis that the ratio of cumulative effective precipitation (direct runoff) P e to cumulative total precipitation P minus initial losses I α is equal to the ratio of current cumulative infiltration F to maximum potential catchment storage S: Water balance is defined by: Substituting F from Equation (9) into Equation (10) results in the following cumulative effective precipitation equation: where: P e is the effective precipitation in mm; P is the total precipitation in mm; I α is the initial abstraction in mm; S is the potential catchment storage in mm; and F is the current infiltration in mm after runoff begins. The effective precipitation amount totalized over time from 0 (the beginning of precipitation) to t (the current instant) is where: P e and P are the effective and the total precipitation amounts totalized over time from 0 to T [mm], λ is an empirical coefficient contingent on the CN (curve number) parameter, which is normally assumed as a constant value of 0.2 in order for S to be the only parameter of the method; and S is the maximum catchment storage [mm].
From the Verma investigations [15], the coefficient λ can have several values that range from 0.05 to 0.3, although it can theoretically be greater than 1.0. This coefficient is of regional character and, as mentioned earlier, depends on many parameters: the region's characteristics and soil, land use, and also climate conditions. If rainfall and runoff investigations are available, the value of λ can be estimated on their basis and then used in the design process. The original value λ = 0.2 was specified for small agriculture catchments. However, later investigations [26] showed that smaller values of λ, e.g., λ = 0.05, give better results for the river catchment for which the investigations were conducted. Despite this, a value λ = 0.2 is still applied in many practical applications [27].
Chung et al. in paper [28] showed the possibility of calculating the λ parameter from the relation λ = 0.2/ε. In this relation, ε only depends on the two parametersα and β, which can be calculated if the following, are given: t p is the ponding time, K s is the saturated permeability coefficient of soil, t is the rainfall duration, and p is the rainfall intensity. Since the last three variables can be measured easily and t p is determined by any adequate formula, the uncertainty of the adopted infiltration models is reduced to a minimum. The maximum catchment storage is given by relation (13): According to the assumptions of the NRCS-CN method, a flood begins when the amount of precipitation exceeds the height of the layer of water retained by the processes of interception, surface storage, and infiltration before a surface runoff starts. The amount of precipitation which is involved in the above processes is referred to as initial abstraction I α . Parameter CN is determined as a weighted average for the whole catchment area on the basis of the adopted soil group, the catchment use, and the hydrologic conditions (14): where: A is the catchment area in km 2 ; CN i is the characteristic values of the particular areas A i ; and n is the number of homogenous areas. The types of soils in the Łęgoń watercourse catchment area were identified on the basis of a soil-agricultural map with a scale of 1:50,000 and also the division of soils included in [29]. This enabled the NRCS-CN method to be directly applied for Polish conditions. The catchment soils were assigned to the types of soils included in the map and aggregated. They were then assigned to one of the four soil groups (A, B, C, D) depending on the potential surface runoff formation conditions stemming from the soils' permeability. Land use classes were determined using the CORINE Land Cover database and a land cover map, which is presented in Figure 2. Three soil groups (Table 1) and seven land uses (Table 2) of the catchment area were distinguished. The CN parameter for the whole catchment amounted to 71. This value was adopted in order to create effective precipitation hyetographs using Equation (12), which are presented in Figure 3. stemming from the soils' permeability. Land use classes were determined using the CORINE Land Cover database and a land cover map, which is presented in Figure 2. Three soil groups (Table 1) and seven land uses (Table 2) of the catchment area were distinguished. The CN parameter for the whole catchment amounted to 71. This value was adopted in order to create effective precipitation hyetographs using Equation (12), which are presented in Figure 3.

Determination of Wave Hydrographs
Hydrographs of the designed and control flood waves for the Łęgoń watercourse were developed using a rainfall-runoff model, which enabled instantaneous unit hydrographs (IUHs) to be determined [30]. The unit hydrograph shows the course of a flood wave (the catchment response) caused by effective precipitation that amounts to one unit and lasts for Δt. The parameters of the flood models were determined using Nash's conceptual model [19,21]. The Nash model shows the transforming action of the catchment, i.e., the transformation of the effective precipitation unit into an instantaneous unit hydrograph of the immediate runoff, which is in the form of a reservoir cascade with linear characteristics. The instantaneous unit hydrograph is the two-parameter gamma distribution probability density function: where: u(t) is the instantaneous unit hydrograph ordinates in h −1 ; t is the time from the coordinate system origin in h; k is the reservoir storage parameter in h determined using the empirical Formula (16); N is the number of reservoirs in the Nash model N = LAG/k, where runoff lag time (delay between the culmination of precipitation and the peak of the flood wave) is determined using Formula (17); and Γ(N) is the gamma function that is equal to (N − 1).
The instantaneous unit hydrographs are the basis for determining flood-wave hydrographs. The ordinates of the surface runoff hydrographs can be calculated using Equation (19):

Determination of Wave Hydrographs
Hydrographs of the designed and control flood waves for the Łęgoń watercourse were developed using a rainfall-runoff model, which enabled instantaneous unit hydrographs (IUHs) to be determined [30]. The unit hydrograph shows the course of a flood wave (the catchment response) caused by effective precipitation that amounts to one unit and lasts for ∆t. The parameters of the flood models were determined using Nash's conceptual model [19,21]. The Nash model shows the transforming action of the catchment, i.e., the transformation of the effective precipitation unit into an instantaneous unit hydrograph of the immediate runoff, which is in the form of a reservoir cascade with linear characteristics. The instantaneous unit hydrograph is the two-parameter gamma distribution probability density function: where: u(t) is the instantaneous unit hydrograph ordinates in h −1 ; t is the time from the coordinate system origin in h; k is the reservoir storage parameter in h determined using the empirical Formula (16); N is the number of reservoirs in the Nash model N = LAG/k, where runoff lag time (delay between the culmination of precipitation and the peak of the flood wave) is determined using Formula (17); and Γ(N) is the gamma function that is equal to (N − 1).  (17) where: A is the catchment area in km 2 ; U is the fraction of impermeable surfaces in the whole catchment; P e is the amount of effective precipitation in mm; and D is the duration of effective precipitation in h. The time of growth amounts to t p = (N − 1)k and the maximum of the unit hydrograph amounts to The instantaneous unit hydrographs are the basis for determining flood-wave hydrographs. The ordinates of the surface runoff hydrographs can be calculated using Equation (19): h k · ∆P e,j , k = i − 1 + j, i = 1, 2, . . . , m + n − 1 (19) where: ∆P e,j is the partial effective precipitation in time interval j in mm; m is the number of the unit hydrograph ordinates; n is the number of effective precipitation time intervals; and h k is the unit hydrograph ordinates in m 3 ·s −1 ·mm −1 calculated using Equation (20), or if t p > 3∆t using Equation (21): The designed and control flood-wave hydrographs for the Łęgoń watercourse calculated by the IUH method are shown in Figure 4. The peak flows amount accordingly to Q 1% = 18.75 m 3 ·s −1 and Q 0.3% = 24.8 m 3 ·s −1 . where: ΔPe,j is the partial effective precipitation in time interval j in mm; m is the number of the unit hydrograph ordinates; n is the number of effective precipitation time intervals; and hk is the unit hydrograph ordinates in m 3 ·s −1 ·mm −1 calculated using Equation (20), or if   3 p t t using Equation (21): The designed and control flood-wave hydrographs for the Łęgoń watercourse calculated by the IUH method are shown in Figure 4. The peak flows amount accordingly to Q1% = 18.75 m 3 ·s −1 and Q0.3% = 24.8 m 3 ·s −1 .

Flood Waves Routing
The flood waves routing was carried out in order to calculate the required volume of the retarding reservoir and to determine its effect on the capacity of the pumping station. Another aim of the calculations was to determine the maximum water storage levels in the reservoir and the reduced water runoffs. A method that uses the continuity equation [31,32] in a differential form (22) was selected from many methods for reservoir flood wave routing: where: Q is the water flow in m 3 s −1 ; x is a coordinate consistent with the direction of the water flow in m; A is the surface area of the flow in m 2 ; and t is the time in s. Having been integrated from x1 to x2 and appropriately transformed, Equation (22) assumes the form of (23) in which the value of the integral expresses the volume of the water stored in the reservoir or that discharged from the reservoir:

Flood Waves Routing
The flood waves routing was carried out in order to calculate the required volume of the retarding reservoir and to determine its effect on the capacity of the pumping station. Another aim of the calculations was to determine the maximum water storage levels in the reservoir and the reduced water runoffs. A method that uses the continuity equation [31,32] in a differential form (22) was selected from many methods for reservoir flood wave routing:

∂Q ∂x
where: Q is the water flow in m 3 s −1 ; x is a coordinate consistent with the direction of the water flow in m; A is the surface area of the flow in m 2 ; and t is the time in s. Having been integrated from x 1 to x 2 and appropriately transformed, Equation (22) assumes the form of (23) in which the value of the integral expresses the volume of the water stored in the reservoir or that discharged from the reservoir: where: x 1 coordinates at the beginning of the reservoir; x 2 coordinates at the end of the reservoir; Q(x 1 ), Q(x 2 ) represent the flow of water into and out of the reservoir in m 3 s −1 , respectively; and ∆V/ ∆t is the change in water volume in the reservoir over time.

Methodology Summary
Below, the authors present a methodology of proceeding when making a correct selection of the capacity of a dyke pumping station, which collaborates with a retarding reservoir with an adequately selected capacity.

1.
Selecting the methodology for estimating calculated flows with a specified probability of exceedance, which depends on whether the catchment is being controlled or not.

2.
Conducting hydrological calculations for designing flows estimation.

4.
Evaluating terrain conditions in the localization of a dyke pumping station, which enables a retarding reservoir with a determined capacity and parameters to be constructed.

5.
Selecting the capacity of pumps in the dyke pumping station. 6.
Selecting computational scenarios that consider the variable capacity of pumps and the variable volume of a retarding reservoir. 7.
Estimating the functioning conditions of a dyke pumping station and a retarding reservoir during the occurrence of a flood wave. 8.
Estimating the costs of construction, exploitation, and maintenance for a dyke pumping station and a retarding reservoir with regards to different computational scenarios.

Pump Station Capacity Assumptions
The specified flood wave hygrographs, the retarding reservoir capacity curves, and the pumping station capacity characteristics were used to calculate several computational cases of flood waves passing through the reservoir. The calculations were performed with regards to the following assumptions:

Calculation Results and Discussion
The results of numerical simulations of the designed flood wave routing through the retarding reservoir for the particular cases, including the capacity of the reservoir and its bottom dimensions, the maximum water level elevation, the retarding reservoir filling level, and the capacity of the pumps, are presented in Table 3. Selected cases are also graphically presented in Figure 5a-f. An analysis of the results, which aimed to determine the required capacity of the retarding reservoir and the required capacity of the pumping station, showed that it is not possible to safely route the designed flood wave through the pumping station in its current state, i.e., the existing retarding reservoir and the pumps' total output amounting to 3.60 m 3 s −1 (Case 1). The capacity of the retarding reservoir is insufficient to reduce the designed flood wave-the water level significantly exceeds the maximum permissible level of 178.00 m by 0.50 m. (Figure 5a). Similar results were obtained in the calculations of flood waves passing through a retarding reservoir while taking into account the change in the pumps' capacity to the total capacity of 4.20 m 3 s −1 (Case 2).
The designed flood wave can be routed through the retarding reservoir if the latter's capacity is increased to about 32,000 m 3 and the pumping station capacity is increased to the capacity determined by the hydrological calculations in both Case 3 ( Figure 5b) and Case 5 (Figure 5d).
During the control flood wave routing (Case 4), the maximum permissible water level can be exceeded up to the level of 178.50. This is due to the fact that the control flow is treated as a verifying flow, and it is actually the designed flow that is used to dimension hydraulic structures. The verifying calculations carried out to determine the effect of the reduction in the total capacity of the pumps and the corresponding readjustment of the retarding reservoir capacity (Cases 6-9) showed that each reduction in the capacity of the pumping station results in an increase in the capacity of this reservoir and its dimensions (Figure 5f-h).
The results of the calculations clearly confirmed the assumptions made in design Case 5, i.e., the best solutions are: the total capacity of the pumps-18.75 m 3 s −1 with the replacement of the pumps in the existing pumping station in order to obtain a total capacity of 4 × 1.05 = 4.20 m 3 s −1 , and the building of another pumping station in order to obtain a capacity of 14.55 m 3 s −1 . The required retarding reservoir capacity is 31,590 m 3 and its bottom dimensions are 45 × 240 m. Moreover, Case 3, which can be executed at a lower cost, can also be considered as it does not require the reconstruction of the existing pumping station, but only requires the replacement of the pumps due to their condition.
Water 2018, 10, x FOR PEER REVIEW 13 of 16 reconstruction of the existing pumping station, but only requires the replacement of the pumps due to their condition.

Conclusions
This paper presents the authors' methodology that was proposed in order to select a correct solution for a dyke pumping station, the functioning conditions of which were strictly connected

Conclusions
This paper presents the authors' methodology that was proposed in order to select a correct solution for a dyke pumping station, the functioning conditions of which were strictly connected with the choice of the capacity of the installed pumps or retarding reservoir. The correct solution should consider exploitation and maintenance costs. The proposed methodology can be used in the design procedure of the other hydroengineering structures. It can also be used in order to estimate flood wave routing through dry and wet storage reservoirs. The presented methodology of proceeding when making a correct selection of a pumping station in order to secure areas that are located above it can be successfully used in the design of wet or dry flood protection reservoirs to protect the areas located below.
The Łęgoń watercourse is an uncontrolled river, and therefore it was not possible to validate the flows that were estimated using the rainfall-runoff method. For vaildation, calculations of the flows were performed using the spatial regression method that is recommended by the Polish Hydrologists Association. The authors obtained small differences that were within the permissible limits, which were equal to −2.4% and +5.2% for the probability of 1% and 0.3%, respectively.
The purpose of every dyke pumping station is to transfer the water of watercourses, the outlets of which are cut off from the main river by the river embankment during high river flows. The transfer capacity of a pumping station should be determined on the basis of hydrological calculations so that no hazard of flooding the areas over the dyke is created. A methodology for calculating the capacity of a pumping station, illustrated using the case of an existing pumping station on the Łęgoń watercourse, was presented. The proposed methodology requires the calculation of the rate of the designed (dimensioning) and control flows and then the determination of the hydrographs of hypothetical waves using the mathematical rainfall-runoff model for each hydraulic structure. For this purpose, the authors used the IUH model and then calculated: the total precipitation, the effective precipitation, the instantaneous unit hydrograph, and the direct (surface and subsurface) runoff hydrograph.
The hydrological analysis and the numerical simulations carried out for the Łęgoń river catchment showed that:

1.
The designed flood wave of p = 1% cannot be safely passed through the pumping station for the existing specifications-the total pumping station capacity of 3.60 m 3 s −1 and the retarding reservoir capacity of 7412 m 3 .

2.
On the basis of computer simulations, it should be indicated that the solution that involves reducing the required pumping station capacity leads to an increase in the required retarding reservoir capacity. The approach to ensure safe operation of the pumping station should not only take into consideration the capacity of the pumping station and volume of the retarding reservoir, but also the retention capacity of the riverbed. This will reduce the cost of rebuilding an existing facility.

3.
On the basis of the performed analysis, it was decided that Case 5 is the most beneficial for reconstructing the dyke pumping station, and in this case, the designed flood wave of p = 1% can be safely passed through the Ciechowice pumping station at the total pumping station capacity of 18.75 m 3 s −1 and the retarding reservoir capacity of 31,590 m 3 .