The E ﬀ ect of Climate Change on Controlled Drainage E ﬀ ectiveness in the Context of Groundwater Dynamics, Surface, and Drainage Outﬂows. Central-Western Poland Case Study

: Control drainage (CD) is a common practice implemented to control the water balance of drainage ﬁelds by increasing the amount of water retained in soil. Worldwide studies suggest that climate change can reduce the e ﬀ ectiveness of CD solutions, but no study of CD outﬂow in comparison to the present situation. The simulations show that the greatest e ﬀ ectiveness of CD solutions will be achieved by starting the blocking of outﬂow from the drainage network on the 1st of March. The implementation of CD solutions since April in the near and far future will allow maintaining the groundwater table at the level presently observed for the conventional network (free drainage-FD). The better e ﬀ ects of the drainage network with spacing of 14 m result from slower drainage of water from the ﬁeld in the period preceding the blocking of outﬂow in comparison with the drainage network with spacing of 7 m. The di ﬀ erences between average groundwater levels for FD and CD implemented from the 1st and 15th of March were statistically signiﬁcant for the near and far future.


Introduction
Poland is a country with very limited water resources with practically every year water deficits (defined as reduction of evapotranspiration caused by meteorological drought) mainly in the period of plant growth. Since the 1990s, actions have been taken to store water for agricultural purposes through the construction of reservoirs [1]. These activities allowed storage of water in winter and spring in reservoirs and its use for irrigation in the period of the agricultural drought. Solutions of this type are extremely expensive, first of all because of the high cost of construction of reservoirs and infrastructure for water allocation. Another source of high cost is the allocation of water and the maintenance of reservoirs collecting water for agricultural purposes. Due to the small depths and high nutrient supply, the reservoirs become intensively overgrown, while as a result of the sedimentation, their volume is reduced [2][3][4]. In addition, reservoirs primarily built for agricultural purposes are also used for recreational, hydropower, fishing, and flood control purposes. Water management of such reservoirs is difficult and most often agricultural purposes play a secondary role. In addition to the construction of reservoirs, due to the occurrence of periodic excess water in spring, a drainage network was developed between the 1950s and 1990s. In Poland, subsurface drainage systems were developed on the surface area of 4.2 × 10 6 ha, which accounted for about 30% of arable land. The drainage systems in the world have been constructed on the area of around 193.9 × 10 6 ha of arable land. However, in 30 countries, the total drained area is more than 1.0 × 10 6 ha [5].
Increasing attention has been paid to the potential threats that agriculture will have to face in future. In particular, the impact of climate change is highlighted on agricultural production risk. The increase in this risk is related to the increasing occurrence and intensity of agricultural and hydrological drought in response to the increase in evapotranspiration and variability of runoff, at a relatively constant rainfall [6]. The climate changes involve the values of minimal, mean, and maximal temperatures as well as the amounts and seasonal patterns of rainfall. Particularly worth highlighting are the rainfall pattern changes, which will lead to the events of intense short-term rainfalls after long rainfall-free periods. Therefore, the risk of hydrological and agricultural drought is expected to grow [6][7][8]. On the other hand, there will be loss of water and nutrients for crops and soil degradation due to increase in the surface outflow, erosion, and nutrient runoff [9]. Increasing frequency of periods with too much and too little water will result in farmers facing complex decisions regarding whether they should adapt by investing in an irrigation system, a drainage system, or both, and when [10]. Farmers mainly consider investing in irrigation systems, while drainage systems are not considered as infrastructure to mitigate the effects of climate change on agriculture. Nowadays, during a period of dramatic water resources changes and their uncertainty, there is a need to move from traditional high water demand cropping and irrigation practices to modern systems and technologies [11]. Misra [12] highlighted the need to develop new and improve existing irrigation infrastructure in all vulnerable countries to cope with climate change risks. However, in the case of deep droughts, the use of irrigation systems may bring limited results, because the water abstraction for irrigation will be limited or even impossible due to low stream flows and dry ponds or wells [13]. Moreover, it is necessary to increase the water use efficiency [11]. Raising water productivity is possible through improving water management practices on agricultural lands [14] and joint actions by adopting proven agronomic and water management practices [15]. The main objective of these activities is to minimize loss of water by evaporation from the soil or percolation beyond the depth of the root zone [6]. Water-use efficiency improvement may effectively mitigate water shortages and intensify the resilience of agricultural systems to climate change [16]. Moreover, climate change will require the choice of appropriate crops, improvements in crop cultivars, maintenance of the drainage and irrigation systems, and optimization of water management [17,18]. Furthermore, Wolfe et al. [10] indicated the need for improved soil water holding A general view of the devices controlling groundwater table control and the triangular (V-notch) weir measuring the subsurface outflow were shown in Figure 1. In subsections 42_2 and 42_3, the drainage network has 7 m spacing, while in sections 42_1 and 42_4 it has 14 m spacing. The drainage network was made of PVC pipes with a diameter of 0.05 m, which were installed at a depth of 0.90 m b.s.l. The groundwater table control structures equipped with the V-notch weirs were installed in drainage wells. Previous research [19] has shown the current possibilities of using CD solutions in Poland and their impact on groundwater table dynamics and subsurface outflow. In this paper, the analysis focuses on the potential impact of climate change on the effectiveness of CD solutions in the measurements of subsurface outflow. The water table level at each plot was measured by piezometers located between two central lateral drains. A Solinst LTC Levelogger Edge pressure sensor (Solinst Canada Ltd., Georgetown, ON, Canada) was installed in each control well. Measurement rate was set to 1 h in each sensor.
Moreover, in the subsections 42_2 and 42_3 also measurements of soil moisture at depths 10, 20, 30, 40, 50, 60 cm were conducted. The volumetric soil water content was measured by means of the frequency domain reflectometry (FDR) technique implemented in WaterScout SM100 soil moisture sensors. The measurements were performed at a once per hour frequency. Subsections 42_3 and 42_4 work in a conventional way, the so-called free drainage (FD), whereas in subsections 42_1 and 42_2 controlled drainage (CD) activities were carried out. At these subsections, simple systems based on the elbow were installed on each well. To regulate outflows, inside of these pipes the modular Poly Vinyl Chloride (PVC) risers were equipped. Additionally, at the end of each control structures the triangular crested weirs for water flow measurement were installed. The outflow logging was based on a pressure level sensor mounted in a white tube (Figure 1).
A general view of the devices controlling groundwater table control and the triangular (V-notch) weir measuring the subsurface outflow were shown in Figure 1. In subsections 42_2 and 42_3, the drainage network has 7 m spacing, while in sections 42_1 and 42_4 it has 14 m spacing. The drainage network was made of PVC pipes with a diameter of 0.05 m, which were installed at a depth of 0.90 m b.s.l. The groundwater table control structures equipped with the V-notch weirs were installed in drainage wells. Previous research [19] has shown the current possibilities of using CD solutions in Poland and their impact on groundwater table dynamics and subsurface outflow. In this paper, the analysis focuses on the potential impact of climate change on the effectiveness of CD solutions in the context groundwater table dynamics and surface and drainage outflows. The data presented in Table 1 were also the input parameters for the DRAINMOD computer simulation model. The study area is slightly heterogeneous in terms of hypsometry; the height of the terrain varies from 107.27 to 109.48 m a.s.l. Slopes are lower than 0.5%. The soils were classified as Gleyic Luvisols [44] that developed from glacial till. Both the surface and subsurface soil layers have a sandy loam texture. The argic horizon has higher clay content than the overlying layers.

DRAINMOD Model Description
DRAINMOD is a field-scale model, which enables the calculation of the water balance components on a daily time step. The calculations are carried out in relation to unit surface area, which is a vertical soil column ( Figure 2). Hydrologic variables such as evapotranspiration (ET), surface runoff (SR), infiltration (I), subsurface drainage (SD), surface water storage change (∆S), deep seepage Agronomy 2020, 10, 625 6 of 21 (DS), and the water-free pore space and air volume (∆Va) are predicted by the model. The input data were precipitation (P) and air temperature (T). The DRAINMOD model was developed by Skaggs [45], but its usefulness has been verified many times and presented in numerous scientific publications [32,41,[46][47][48].

DRAINMOD Model Description
DRAINMOD is a field-scale model, which enables the calculation of the water balance components on a daily time step. The calculations are carried out in relation to unit surface area, which is a vertical soil column ( Figure 2). Hydrologic variables such as evapotranspiration (ET), surface runoff (SR), infiltration (I), subsurface drainage (SD), surface water storage change (∆S), deep seepage (DS), and the water-free pore space and air volume (∆Va) are predicted by the model. The input data were precipitation (P) and air temperature (T). The DRAINMOD model was developed by Skaggs [45], but its usefulness has been verified many times and presented in numerous scientific publications [32,41,[46][47][48].

Model Inputs and Climate Change Impact Manipulation
The identification of model parameters and the procedure of the model calibration and validation were carried out using the procedure described by Skaggs et al. [45]. Groundwater table measurements from 2000 were used to calibrate the model. The hydraulic conductivity of layer/horizon, especially below the soil depth, the thickness of the restrictive layer, and the hydraulic head at the bottom of the restrictive layer, within reasonable ranges, were adjusted to minimize the difference between the observed and predicted groundwater table depths. After calibration, the model was validated on the basis of groundwater table measurements from 1994 using standard statistical measures, such as the root mean square error (RMSE), coefficient of mass residual (CRM), index of agreement, (d) and model efficiency index (EF) [49,50]. These statistics are defined as:

Model Inputs and Climate Change Impact Manipulation
The identification of model parameters and the procedure of the model calibration and validation were carried out using the procedure described by Skaggs et al. [45]. Groundwater table measurements from 2000 were used to calibrate the model. The hydraulic conductivity of layer/horizon, especially below the soil depth, the thickness of the restrictive layer, and the hydraulic head at the bottom of the restrictive layer, within reasonable ranges, were adjusted to minimize the difference between the observed and predicted groundwater table depths. After calibration, the model was validated on the basis of groundwater table measurements from 1994 using standard statistical measures, such as the root mean square error (RMSE), coefficient of mass residual (CRM), index of agreement, (d) and model efficiency index (EF) [49,50]. These statistics are defined as: where n is the total number of the observations, O i is the observed value of the ith observation, P i the predicted value of the ith observation, and O the mean of the observed values (i = 1 to n). The RMSE, CRM, d, and EF values showed good agreement between the measured and predicted groundwater table depths. The values of these statistics were 0.04 m, 4%, 0.95, and 0.78, respectively. The identification of model parameters and the procedure of the model calibration and validation were presented by Sojka et al. 2019 [19].
The model parameters were identified on the basis of soil properties (Table 1), technical parameters of the drainage network and crop characteristics ( Table 2). Soil water retention properties up to −100 kPa were determined using the Richards chambers, whereas for lower values of the pressure head using the method of water vapor pressure over a solution of sulfuric acid. The saturated hydraulic conductivity was determined using the method of a constant hydraulic water gradient. The soil package included in DRAINMOD was used to estimate the Green-Ampt infiltration model parameters, the drainage volume-water table depth, and the water upflux-water table depth relationships. Details of soil properties and drainage network parameters were presented by Sojka et al. [19]. The climate data included daily precipitation and daily minimum, average and maximum air temperatures. During the DRAINMOD simulation, daily precipitation was uniformly distributed over 6 h (15 to 21 h) to obtain hourly precipitation according to the Skaggs [45] procedure. The ET value calculated by the Thornthwaite method was corrected using the monthly ET factors. In the first stage of research, 18 scenarios of the CD network functioning in the current climatic conditions were analyzed, assuming different meteorological conditions (dry, normal, and wet years-2014, 2017, and 2018), initial groundwater table conditions (groundwater table depth at 0.4 m, 0.6 m, and 0.7 m on 1st of March-start of the simulation), drainage network parameters (drainage spacing 7 and 14 m), and the FD drainage practice. The precipitation and minimum and maximum temperature for near and far future periods (2021-2050 and 2071-2100) were input data in DRAINMOD model. In these studies, 180 simulations were carried out, including 72 simulations for CD and 18 for FD for each time period, 2021-2050 (near future) and 2071-2100 (far future) ( Table 2).
The results of simulation are daily values of depth to groundwater table (DTWT), surface runoff (SR) and subsurface drainage (SD). The calculation scheme is presented in Figure 3. The statistical differences between the values of mean DTWT, SR, and SD in different scenarios in the present, near and far future, were assessed using a non-parametric analysis of the variance performed using the STATISTICA 13.1 PL statistical software, with Kruskal-Wallis (K-W) and Dunn's tests as post hoc procedures (p ≤ 0.05). The results of simulation are daily values of depth to groundwater table (DTWT), surface runoff (SR) and subsurface drainage (SD). The calculation scheme is presented in Figure 3. The statistical differences between the values of mean DTWT, SR, and SD in different scenarios in the present, near and far future, were assessed using a non-parametric analysis of the variance performed using the STATISTICA 13.1 PL statistical software, with Kruskal-Wallis (K-W) and Dunn's tests as post hoc procedures (p ≤ 0.05).

Climate Change Projection
Climate change projection for the study area was determined by using the resources of the climateimpact.sggw.pl Geoportal [51,52]. The Geoportal was developed in the framework of the Polish-Norwegian Research Programme, within the project CHASE-PL (Climate Change Impact on Selected Sectors in Poland). Temperature and precipitation projections for two time horizons have been used: near future (2021-2050) and far future (2071-2100). The maps of climate projection for precipitation, minimum, and maximum temperature interpolated on a 5 km grid are available for the whole area of Poland [51,52]. Projections of climate change in Poland were developed on the basis of downscaling of General Circulation Models (GCMs). In the CHASE-PL project, the analyses were carried out using a set of climate projections comprising nine regional climate model outputs (biascorrected and downscaled to a 5 km × 5 km grid) stemming from the Coordinated Downscaling Experiment-European Domain (EURO-CORDEX) ensemble [53]. In this study, the data from the RCP 4.5 scenario representing intermediate greenhouse gas emissions were taken into account. The data

Climate Change Projection
Climate change projection for the study area was determined by using the resources of the climateimpact.sggw.pl Geoportal [51,52]. The Geoportal was developed in the framework of the Polish-Norwegian Research Programme, within the project CHASE-PL (Climate Change Impact on Selected Sectors in Poland). Temperature and precipitation projections for two time horizons have been used: near future (2021-2050) and far future (2071-2100). The maps of climate projection for precipitation, minimum, and maximum temperature interpolated on a 5 km grid are available for the whole area of Poland [51,52]. Projections of climate change in Poland were developed on the basis of downscaling of General Circulation Models (GCMs). In the CHASE-PL project, the analyses were carried out using a set of climate projections comprising nine regional climate model outputs (bias-corrected and downscaled to a 5 km × 5 km grid) stemming from the Coordinated Downscaling Experiment-European Domain (EURO-CORDEX) ensemble [53]. In this study, the data from the RCP 4.5 scenario representing intermediate greenhouse gas emissions were taken into account. The data concerning the projection of seasonal temperature and precipitation changes in the near (Table 3) and  (Table 4) were used. The obtained results indicate that in the analyzed area there will be an increase in temperature and precipitation. In the near future, in spring and summer, the average temperature will increase by about 1.0 • C and in the far future by 1.8 • C (Table 3). Both in the near and far future, in spring, and autumn, it is expected that seasonal precipitation will increase by 9.11 and 2.47% as well as by 16.44 and 3.76%, respectively (Table 4). The gradual and variable approaches were applied to modify the historical daily climate data (2014, 2017, and 2018). The gradual algorithm was used for temperature. The modification involved adding the values specified in Table 3 to the minimum and maximum daily temperatures. The gradual modifications did not change the timing of extreme climatic events. Thus, drought and extreme rainfall events in the future climate for the near and far future occurred in the same frequency as in the historical data. In Poland, it is assumed that the frequency and duration of droughts and extreme precipitation will increase. The variable modification of the precipitation was made to achieve more dry days, greater contiguous dry-day durations, and more intense rainfall events on the basis of methods proposed by Harrison et al. [39]. In this method, small, medium, and large rainfall events were separately manipulated by different seasonal factors. The seasonal factors (CsP) were calculated on the basis of the precipitation change (%) values presented in Table 4. For example, the CsP value for spring in the near future is 1.0911, representing a 9.11% increase in rainfall. Subsequently, on the basis of daily precipitation data from 1992-2018 for spring, summer and autumn seasons, the percentile values were computed. On the basis of the percentile values, the daily precipitation data were classified into three groups. Small, medium, and large daily precipitation (Pi) were multiplied by different factors (FPi) to achieve the required overall seasonal change in precipitation: The Solver routine in Microsoft Excel was used to optimize FPi across years y and percentiles p: In this way, on the basis of the series of historical meteorological data (2014, 2017, and 2018), a series of meteorological data for the near (2014N, 2017N, and 2018N) and far future (2014F, 2017F, and 2018F) were generated.

Results
The analysis of the effectiveness of CD solutions is presented in relation to the periods from March to September for the years 2014, 2017, and 2018, which are characteristic in terms of precipitation. According to the methodology presented above, daily precipitation was modified in order to present the functioning of drainage networks and the effectiveness of CD solutions in the near and far future. In this way, new daily rainfall series were created for the period from March to September for the near (2014N, 2017N, and 2018N) and far future (2014F, 2017F, and 2018F). The modification of precipitation caused that in all analyzed years the length of rain-free periods increased and the number of days with daily precipitation above 0.1 mm was reduced. The number of days with precipitation for 2014, 2017, and 2018 will decrease by 33%, 43%, and 49% respectively in the future. In the near future, between March and September, rainfall will increase by 4% to 13%, on average by 8%. A slightly larger increase in precipitation will occur in the far future from 7% to 16%, on average by 11%. However, in the case of maximum daily precipitation, climate change in the near and far future will result in an increase of around 50%. The results of rainfall modification in individual months are presented in Table 5. In accordance with the research scheme presented in the methodology (Figure 3), a series of simulations was carried out to show how the drainage network will be functioning in the future and what the effectiveness of CD solutions will be. The influence of climate change on the functioning of Agronomy 2020, 10, 625 11 of 21 drainage networks was presented in relation to the depth of groundwater table, subsurface outflows from drainage networks, and surface outflows.

Groundwater
According to the results of the simulations performed, at present in the fields with traditional drainage systems of 7 m spacing, the groundwater table is above the drainage network level from 27 to 51 days, on average for 40 days (Figure 4a) in the period from 1st of March to 30th of September. In the fields equipped with the drainage network of 14 m spacing, the groundwater table remains above the drainage network level for the period of 40 to 66 days, on average 53 days (Figure 4b).
The predicted climate changes in temperature and amount of precipitate imply that the period in which the groundwater table will be above the drainage network level will get shorter. In the fields equipped with the drainage network of 7 m spacing the groundwater table will be above the drainage network level on average for 35 and 33 days in the near and far future, respectively. Moreover, in the fields with the drainage network of 14 m spacing, the period in which the groundwater table will be above the drainage network level will shorten to 49 and 46 days, on average, for the near and far future, respectively ( Figure 4b).
Agronomy 2020, 10, x FOR PEER REVIEW 11 of 21 predicted climate changes in temperature and amount of precipitate imply that the period in which the groundwater table will be above the drainage network level will get shorter. In the fields equipped with the drainage network of 7 m spacing the groundwater table will be above the drainage network level on average for 35 and 33 days in the near and far future, respectively. Moreover, in the fields with the drainage network of 14 m spacing, the period in which the groundwater table will be above the drainage network level will shorten to 49 and 46 days, on average, for the near and far future, respectively (Figure 4b).
(a) (b) The simulation has shown that the use of CD solutions induced changes in the groundwater dynamics. The best results were obtained for the scenarios with the start of blocking the outflow on the 1st of March. The use of CD solution in the near future on the drainage network with spacing of 7 m implies that the groundwater will be higher than the drainage network level for 56 to 103 days (on average 85 days) (Figure 4a). In the far future, this time may be reduced by three days on average. The implementation of CD solutions from the 15th of March in the drainage network with spacing of 7 m means that groundwater table will be above the drainage network level on average for 70 days and 64 days in the near and far future (Figure 4a). The start of blocking the outflow on the 1st or 15th of April, both in the near and far future, will have little effect on the groundwater table dynamics. The time for the groundwater table to be above the drainage network level, for the scenarios of the 1st and 15th of April, will be, on average, 41 and 37 days for the near future and 34 and 36 days for the far future. These values are at a similar or even lower level than those presently observed in the FD network. The best results from CD implementation in the 14 m spacing drainage network will be achieved for the scenarios with the start of blocking the outflow on the 1st of March (Figure 4b). The average time of groundwater table remaining above the drainage network level will be similar as in the drainage network with spacing of 7 m. On the other hand, when CD solutions are applied from the 15th of March, the periods of groundwater table being above the drainage network level in the near and far future will be 76 and 73 days, respectively, and will be about 8 days longer than for the drainage network with 7 m spacing. The start of blocking the outflow in April will cause that the time of groundwater table remaining above the level of the drainage network will be from 49 to 53 days The simulation has shown that the use of CD solutions induced changes in the groundwater dynamics. The best results were obtained for the scenarios with the start of blocking the outflow on the 1st of March. The use of CD solution in the near future on the drainage network with spacing of 7 m implies that the groundwater will be higher than the drainage network level for 56 to 103 days (on average 85 days) (Figure 4a). In the far future, this time may be reduced by three days on average. The implementation of CD solutions from the 15th of March in the drainage network with spacing of 7 m means that groundwater table will be above the drainage network level on average for 70 days and 64 days in the near and far future (Figure 4a). The start of blocking the outflow on the 1st or 15th of April, both in the near and far future, will have little effect on the groundwater table dynamics. The time for the groundwater table to be above the drainage network level, for the scenarios of the 1st and 15th of April, will be, on average, 41 and 37 days for the near future and 34 and 36 days for the far future. These values are at a similar or even lower level than those presently observed in the FD network. The best results from CD implementation in the 14 m spacing drainage network will be achieved for the scenarios with the start of blocking the outflow on the 1st of March (Figure 4b).
The average time of groundwater table remaining above the drainage network level will be similar as in the drainage network with spacing of 7 m. On the other hand, when CD solutions are applied from the 15th of March, the periods of groundwater table being above the drainage network level in the near and far future will be 76 and 73 days, respectively, and will be about 8 days longer than for the drainage network with 7 m spacing. The start of blocking the outflow in April will cause that the time of groundwater table remaining above the level of the drainage network will be from 49 to 53 days and from 46 to 47 days, in the near and far future, respectively (Figure 4b). Implementation of CD solutions in April is less effective, the groundwater tables are at the level presently observed in the FD network.
At present, taking into account the FD drainage network with spacing of 7 m, the groundwater level varied between 118 and 128 cm b.s.l. (Figure 5a), with an average value of 123 cm b.s.l. However, in the drainage network with spacing of 14 m, the groundwater tables were located at the level of 113 to 126 cm b.s.l. (Figure 5b), with an average value of 118 cm b.s.l. In the near and far future, due to the increase in air temperatures and precipitation as well as changes in the precipitation patterns, the average groundwater level in the drained fields will decrease from 2 to 5 cm. The differences between the average groundwater levels for FD at present and in the near and far future are statistically insignificant. The application of CD solutions when starting to block the outflow on the 1st of March will allow the average groundwater level to increase in the near future on average by 26 and 22 cm in the fields equipped with the drainage network with spacing of 7 and 14 m, respectively. Slightly smaller effects will be achieved in the far future for this type of solution, on average by about 2 cm. Smaller effects in terms of raising the groundwater level can be achieved in the near and far future by starting to block the outflow on the 15th of March, with the average of 10 and 8 cm and 12 and 10 cm for the drainage network of 7 and 14 m spacing, respectively. The better effects of the drainage network with spacing of 14 m result from slower drainage of water from the field in the period preceding the blocking of outflow in comparison with the drainage network with spacing of 7 m. The differences between average groundwater levels for FD and CD implemented from the 1st and 15th of March were statistically significant for the near and far future. to 126 cm b.s.l. (Figure 5b), with an average value of 118 cm b.s.l. In the near and far future, due to the increase in air temperatures and precipitation as well as changes in the precipitation patterns, the average groundwater level in the drained fields will decrease from 2 to 5 cm. The differences between the average groundwater levels for FD at present and in the near and far future are statistically insignificant. The application of CD solutions when starting to block the outflow on the 1st of March will allow the average groundwater level to increase in the near future on average by 26 and 22 cm in the fields equipped with the drainage network with spacing of 7 and 14 m, respectively. Slightly smaller effects will be achieved in the far future for this type of solution, on average by about 2 cm. Smaller effects in terms of raising the groundwater level can be achieved in the near and far future by starting to block the outflow on the 15th of March, with the average of 10 and 8 cm and 12 and 10 cm for the drainage network of 7 and 14 m spacing, respectively. The better effects of the drainage network with spacing of 14 m result from slower drainage of water from the field in the period preceding the blocking of outflow in comparison with the drainage network with spacing of 7 m. The differences between average groundwater levels for FD and CD implemented from the 1st and 15th of March were statistically significant for the near and far future. The implementation of CD solutions at the start of blocking the outflow in April, in the near and far future will make it possible to maintain the groundwater table at a level similar to that presently observed in the free drainage network (FD) ( Figure 5). It should be noted, however, that variations between individual years (dry and wet) will be more pronounced in the near and far future than at present. The obtained results indicate that the effects of raising the groundwater table in the field for CD solutions implemented from the 1st of March are mainly influenced by the initial conditions (groundwater depth) together with the delay in blocking the outflow from the drainage network. The effect of the initial conditions decreases. However, the effect of meteorological conditions is most visible in the scenario with the starting date of blocking the outflow on the 15th of March. Generally, the meteorological conditions have a greater influence on the CD effects than the initial depth of groundwater.  The implementation of CD solutions at the start of blocking the outflow in April, in the near and far future will make it possible to maintain the groundwater table at a level similar to that presently observed in the free drainage network (FD) (Figure 5). It should be noted, however, that variations between individual years (dry and wet) will be more pronounced in the near and far future than at present. The obtained results indicate that the effects of raising the groundwater table in the field for CD solutions implemented from the 1st of March are mainly influenced by the initial Agronomy 2020, 10, 625 13 of 21 conditions (groundwater depth) together with the delay in blocking the outflow from the drainage network. The effect of the initial conditions decreases. However, the effect of meteorological conditions is most visible in the scenario with the starting date of blocking the outflow on the 15th of March. Generally, the meteorological conditions have a greater influence on the CD effects than the initial depth of groundwater.
The impact of individual CD solution for the 7 m and 14 m spacing drainage networks on the groundwater table rise in individual months in the near and far future is presented in Figure 6. The results obtained indicate that the best effects in terms of groundwater table rise are obtained when the initial groundwater level is 0.4 m b.s.l. Moreover, the groundwater table rise starts on the 1st of March. In the period from March to September in the near future, in the drainage network with spacing of 7 m, the average increase in groundwater table will be from 28 to 34 cm (Figure 6a) and for the drainage network with spacing of 14 m, from 24 to 32 cm (Figure 6b). In the far future, the effects of using CD will be at a slightly lower level. In the case of the CD scenarios with the start of blocking the outflow on the 1st of March, the initial conditions, i.e., the groundwater level at the beginning of the simulation period, have a pronounced impact on the effects of raising the groundwater table. In the case of drainage networks with spacing of 7 m, blocking the outflow from the 15th of March will result in an average increase in the groundwater table by about 2.5 times smaller than that obtained for the CD 1st March scenario. A smaller reduction of the impact on groundwater levels caused by the delay of CD implementation is predicted for the drainage network with 14 m spacing. The start of blocking the outflow from the network on the 15th of March will be approximately 1.5 times less effective on average than for the scenarios of blocking it from the 1st of March. If CD solutions are used to start blocking the drain on the 1st or 15th of April, the maximum groundwater table rise will be 6 cm. There will also be situations where the groundwater table may be lowered by approximately 6 cm compared to that presently observed in the drainage networks. for the drainage network with spacing of 14 m, from 24 to 32 cm (Figure 6b). In the far future, the effects of using CD will be at a slightly lower level. In the case of the CD scenarios with the start of blocking the outflow on the 1st of March, the initial conditions, i.e., the groundwater level at the beginning of the simulation period, have a pronounced impact on the effects of raising the groundwater table. In the case of drainage networks with spacing of 7 m, blocking the outflow from the 15th of March will result in an average increase in the groundwater table by about 2.5 times smaller than that obtained for the CD 1st March scenario. A smaller reduction of the impact on groundwater levels caused by the delay of CD implementation is predicted for the drainage network with 14 m spacing. The start of blocking the outflow from the network on the 15th of March will be approximately 1.5 times less effective on average than for the scenarios of blocking it from the 1st of March. If CD solutions are used to start blocking the drain on the 1st or 15th of April, the maximum groundwater table rise will be 6 cm. There will also be situations where the groundwater table may be lowered by approximately 6 cm compared to that presently observed in the drainage networks.

Subsurface Outflow
According to the recent trends in water management in agriculture, much attention is paid to the subsurface outflows from drainage fields. Successful management of these outflows may lead to elevation of the groundwater table and improvement of soil moisture conditions. At present, the subsurface outflows from the FD network from March to September are at the level of 63, 53 and 28 mm (spacing of 7 m) and 59, 42, and 25 mm (spacing of 14 m), respectively, relative to the initial groundwater conditions of 0.4, 0.6, and 0.8 m b.s.l. (Figure 7). In the near future, the subsurface outflow from the FD network will be at a slightly lower level. Depending on the initial conditions, the outflows will be on average of 59, 42, and 24 mm (7 m spacing) and 55, 38, and 21 mm (14 m spacing). The subsurface outflows from the FD network will also be reduced in the far future.

Subsurface Outflow
According to the recent trends in water management in agriculture, much attention is paid to the subsurface outflows from drainage fields. Successful management of these outflows may lead to elevation of the groundwater table and improvement of soil moisture conditions. At present, the subsurface outflows from the FD network from March to September are at the level of 63, 53 and 28 mm (spacing of 7 m) and 59, 42, and 25 mm (spacing of 14 m), respectively, relative to the initial groundwater conditions of 0.4, 0.6, and 0.8 m b.s.l. (Figure 7). In the near future, the subsurface outflow from the FD network will be at a slightly lower level. Depending on the initial conditions, the outflows will be on average of 59, 42, and 24 mm (7 m spacing) and 55, 38, and 21 mm (14 m spacing). The subsurface outflows from the FD network will also be reduced in the far future. Therefore, the amount of available water for use will be reduced. The differences between the outflow values from the FD drainage network at present and for the near and far future were statistically insignificant. Implementation of CD solutions from the 1st of March will allow the best results to be achieved in terms of reduction of subsurface outflow from the drainage network both in the near and far future. Slightly better effects in terms of reduction of subsurface outflow from the network will be possible to achieve for the 14 m drainage network spacing than for 7 m spacing. If the groundwater table is at the level of 0.8 m b.s.l., it is possible to reduce the outflow from the drainage network totally when blocking the outflow from the 1st of March. The results indicate that early start of CD activities is particularly important for 7 m network spacing. Delay in the blocking of the outflow to the 15th of March will result in outflows from the network at the initial groundwater level of 0.4 m b.s.l. of 44 mm on average. In general, the differences between the outflows for FD and the scenarios with blocking outflows from the drainage network from the 1st and 15th of March are statistically significant both for the near and far future. The largest outflows from the drainage network occur in the case of CD scenarios with the start of blocking the outflow in April. In the climatic conditions of Poland, outflows from the drainage network mainly take place in March. It should be noted that a larger reduction of outflow for each CD scenarios is achieved for the drainage networks with spacing of 14 m. Changes in thermal conditions, and in particular in the amount of precipitation in the near and far future, will affect the functioning of drainage systems in Poland. In general, in FD networks, regardless of the initial conditions of the groundwater table, there will be a tendency to decrease the subsurface outflows from drainage networks (Figure 8). The decrease of outflow from FD drainage networks will be at the average level of 11% and 17% as well as 12% and 18% in the near and far future, respectively, for the networks with spacing of 7 and 14 m. However, it should be noted that in the wet years (2017N and 2017F) there may be a slight increase in the outflows from the drainage networks at the level of 18% and 10% in the near and far future. The average reduction of outflow from drainage networks at the start of blocking outflow on the 1st of March for the drainage networks with spacing of 7 and 14 m will be 82% and 85% on average and 84% and 87% in the near and far future, respectively. Around twice as small effects will be achieved for CD scenarios when blocking the outflow is started on the 15th of March. In the case of other scenarios, the reduction of outflows from the 7 m network spacing will be at an average level of 8% to 30% and for the 14 m network spacing at an average level of 10% to 38%. Changes in thermal conditions, and in particular in the amount of precipitation in the near and far future, will affect the functioning of drainage systems in Poland. In general, in FD networks, regardless of the initial conditions of the groundwater table, there will be a tendency to decrease the subsurface outflows from drainage networks (Figure 8). The decrease of outflow from FD drainage networks will be at the average level of 11% and 17% as well as 12% and 18% in the near and far future, respectively, for the networks with spacing of 7 and 14 m. However, it should be noted that in the wet years (2017N and 2017F) there may be a slight increase in the outflows from the drainage networks at the level of 18% and 10% in the near and far future. The average reduction of outflow from drainage networks at the start of blocking outflow on the 1st of March for the drainage networks with spacing of 7 and 14 m will be 82% and 85% on average and 84% and 87% in the near and far future, respectively. Around twice as small effects will be achieved for CD scenarios when blocking the outflow is started on the 15th of March. In the case of other scenarios, the reduction of outflows from the 7 m network spacing will be at an average level of 8% to 30% and for the 14 m network spacing at an average level of 10% to 38%.

Surface Runoff
One of the important components of the water balance is surface runoff from the field. It is particularly interesting for increasing rainfall intensity in the near and far future due to climate change. The results of the simulation show that in the present climatic conditions in the period from March to September, the surface runoffs are small and range from about 2 to 3 mm (Figure 9). The increase in rainfall intensity in the near and far future will result in an increase in the proportion of surface runoff from the fields. Calculations have shown that regardless of the scenarios of the FD or CD drainage network, in the near future outflows will range from approximately 7.5 to approximately 12.5 mm (9.5 mm on average) and in the far future from 7.3 to 11.8 mm (9.0 mm on average). These values are almost five times higher than at present (statistically significant difference), which means that their effects need to be seriously considered from the point of view of erosion and nitrogen and phosphorus uptake into surface waters.

Surface Runoff
One of the important components of the water balance is surface runoff from the field. It is particularly interesting for increasing rainfall intensity in the near and far future due to climate change. The results of the simulation show that in the present climatic conditions in the period from March to September, the surface runoffs are small and range from about 2 to 3 mm (Figure 9). The increase in rainfall intensity in the near and far future will result in an increase in the proportion of surface runoff from the fields. Calculations have shown that regardless of the scenarios of the FD or CD drainage network, in the near future outflows will range from approximately 7.5 to approximately 12.5 mm (9.5 mm on average) and in the far future from 7.3 to 11.8 mm (9.0 mm on average). These values are almost five times higher than at present (statistically significant difference), which means that their effects need to be seriously considered from the point of view of erosion and nitrogen and phosphorus uptake into surface waters.

Surface Runoff
One of the important components of the water balance is surface runoff from the field. It is particularly interesting for increasing rainfall intensity in the near and far future due to climate change. The results of the simulation show that in the present climatic conditions in the period from March to September, the surface runoffs are small and range from about 2 to 3 mm (Figure 9). The increase in rainfall intensity in the near and far future will result in an increase in the proportion of surface runoff from the fields. Calculations have shown that regardless of the scenarios of the FD or CD drainage network, in the near future outflows will range from approximately 7.5 to approximately 12.5 mm (9.5 mm on average) and in the far future from 7.3 to 11.8 mm (9.0 mm on average). These values are almost five times higher than at present (statistically significant difference), which means that their effects need to be seriously considered from the point of view of erosion and nitrogen and phosphorus uptake into surface waters.

Discussion
Analyses of the functioning of the drainage network in the near and far future indicate that the effectiveness of the use of CD solutions will decrease with respect to the present situation. However, Pease et al. [54] indicated the benefits of controlled drainage will still be evident under projected future climate. With respect to the previous results on the effectiveness of CD solutions [19], in the future, the impact on groundwater will be reduced from 2 to 8 cm in the case of 7 m spacing drainage networks and from 4 to 10 cm in the case of scenarios with the start of blocking the outflow on 1st and 15th of March. On the other hand, for networks with a 14 m spacing, the decrease in the effectiveness of the impact on the ground water level will be greater for the start date of the outflow blocking on 1st of March from 10 to 12 cm than for the start date of blocking the outflow on 15th of March from 6 to 8 cm. The CD treatment produced significantly higher evapotranspiration and soil water content than the conventional drainage treatment during the dry years [55].
The results also indicated that CD practice can reduce drainage outflow significantly [19]. Pease et al. [54] have indicated that CD efficiency was projected to increase during summer of a 42.7% reduction in subsurface drainage discharge for RCP 4.5 under late-century climate conditions. Although tile drainage volume varied amongst the different climate scenarios, NO 3 − -N loss and flow-weighted average NO 3 − -N concentration will increase [56]. Thus, also the water quality in waterbodies can be improved. In a model study, Gilliam and Skaggs [57] have reported that the controlled drainage resulted in the reduction of the nitrate efflux by as much as 34% but the reduction varied with soil and management conditions. Wesström [58] in the field experimental project has pointed out a significant total reduction in nitrate losses. In humid continental climates Sunohara et al. [59,60] have confirmed that the CD retrofitted to existing tile drainage networks can effectively reduce tile drainage water, NO 3 − -N, NH 4+ -N, DRP, total P, and fecal indicator bacteria fluxes to surface waterbodies. The obtained results indicate that the decrease in the effectiveness of the impact of CD solutions may not allow proper water conditions to be ensured during the whole period of plant growth. Baule et al. [61] have pointed out that CD systems alone will not be able to mitigate the effects of climate change in agriculture. Only the combination of CD solutions and subirrigation systems will provide a viable climate adaptation strategy for agriculture [61]. According to Bonaiti and Borin [62], the advantage of the controlled drainage and subirrigation is that they can be used for continuous groundwater table management regardless of the weather conditions. In addition, attention should be paid to other possibilities for increasing soil water retention capacity [11]. Key factors for agriculture's climate change mitigation and adaptation are actions that will contribute to increasing water storage capacity of the soil and providing a larger zone for roots to grow, enabling them to access larger quantities of water and nutrients [63]. The results presented in the paper should be treated with caution because they concern specific conditions, i.e., the drainage object, which is characterized by specific geological conditions, soil conditions, topography, drainage network parameters, as well as hydrological and meteorological conditions. It is well known that the effectiveness of the CD practice can vary across the region with the variation in environmental conditions and management practices [64]. The present research has shown that there will be a decrease in outflows from the drainage network in the period from March to September, with an increase in precipitation in that period. The use of CD solutions will reduce outflows from drainage networks in the future. In terms of subsurface outflows, the results are generally within the ranges given by Skaggs et al. [29] and Cooke and Verma [26]. Wang et al. [56] have shown that under the future climate change the subsurface outflow will increase with increasing mean yearly precipitation. Similar results have been reported by Dayyani et al. [65] and Golmohammadi et al. [66]. Moreover, the studies by Mehan et al. [67] indicate a potential increase in surface and subsurface outflows due to climate change. This study shows that the potential growth of surface runoff should be taken into account in the future. Mehan et al. [67] have reported that more nutrient losses were predicted due to extreme precipitation events with surface outflow. Therefore, increasing number of studies have been focused on the methods of reducing the leaching of nitrogen and phosphorus from arable soils [68].
Duffková et al. [69] have suggested that on the basis of categorization of soil water availability, it is possible to propose location of irrigation systems, the regulation of drainage runoff or to consider controlled drainage. The main challenges are to establish a plan to manage CD systems, adapted to the existing environmental and climatic conditions and, above all, to the water requirements of crops and the dates of field work [70,71]. Most often, simple control rules resulting from natural climatic conditions, adjusted to agrotechnical works and water needs of plants, are established. The effectiveness of such systems can be increased by changing the control rules to more flexible ones, taking into account weather forecasting. In this way, the Climate Adaptive Drainage (CAD) system has been developed [72]. The CAD system consists of a controlled drainage system, a remote-controlled adjustable drain outlet, and a telemetry and database system to process the weather forecast.

Conclusions
The research carried out showed that the future climate change in the central-western part of Poland will extend the rain-free periods, the number of days with precipitation above 0.1 mm will decrease by 40% and the intensity of maximum daily precipitation will increase by 50%. In addition, the increase in air temperatures will be accompanied by an increase in precipitation by an average of 8% and 11% in the near and far future. Assuming such climate change scenarios, the following detailed conclusions concerning the functioning of drainage networks and CD solutions have been formulated on the basis of the conducted analyses: • in fields with FD systems, a reduction of groundwater table of 2 to 5 cm on average is expected in the near and far future in relation to the present conditions, • an increase in the variation in water levels between wet and dry years is to be expected in the near and far future, • the time of the groundwater table above the drainage network will be reduced by 5 and 7 days and by 4 and 7 days respectively in the near and far future for networks with spacing of 7 and 14 m, • the sub-surface outflows will be reduced at the average level of 11% and 17% and 12% and 18% respectively in the near and far future for networks with spacing of 7 and 14 m, • the share of surface run-off from fields is expected to increase approximately fivefold in the future as a result of climate change compared to the present situation, • regardless of the initial conditions of groundwater table, the best results from CD solutions were obtained for the scenarios with the start of blocking the outflow on the 1st of March and 15th of March, • greater effects of CD will be achieved for the 7 m network spacing if the outflow from the drainage network is blocked on 1st of March, and for the 14 m network spacing if the outflow is blocked on 15th of March, • the use of CD solutions in the near and far future will allow similar efficiency to be achieved in relation to groundwater table rise and reduction of subsurface outflow from the drainage network, • using CD solutions in April in the near and far future will allow the groundwater table to be maintained at a level similar to that currently observed in the conventional network (free drainage, FD), • the change of the date of blocking the outflow from the drainage network from 1st to 15th of March will result in a decrease in efficiency by an average of 2.5 and 1.5 times in the case of drainage networks with spacing of 7 m and 14 m.
In the future, the research should be extended to other drainage fields around Poland, which are characterized by specific geological conditions, soil conditions, topography, drainage network parameters, as well as hydrological and meteorological conditions. Studies should be carried out to establish CD principles suitable for the crops grown and available agronomic techniques.