Next Article in Journal
Assessing the Impact of the Strictest Water Resources Management Policy on Water Use Efficiency in China
Previous Article in Journal
Gastropod Assemblages Associated with Habitat Heterogeneity and Hydrological Shifts in Two Shallow Waterbodies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulating the Effects of Agricultural Adaptation Practices onto the Soil Water Content in Future Climate Using SWAT Model on Upland Bystra River Catchment

The Institute of Soil Science and Plant Cultivation—State Research Institute, ul. Czartoryskich 8, 24-100 Pulawy, Poland
*
Authors to whom correspondence should be addressed.
Water 2022, 14(15), 2288; https://doi.org/10.3390/w14152288
Submission received: 7 March 2022 / Revised: 11 July 2022 / Accepted: 12 July 2022 / Published: 22 July 2022
(This article belongs to the Topic Water Management in the Era of Climatic Change)

Abstract

:
The article presents predicted changes in soil water content in the Bystra river catchment (eastern Poland) for various scenarios of climate change and adaptation practices obtained on the basis of a SWAT model simulation for three regional climate models driven by the global climate model EC-EARTH for the years 2041–2050 and the RCP 4.5 and 8.5 RCP scenarios. Climate scenarios were put against five adaptation scenarios presenting changes in land use and protective measures compared against a zero scenario of BaU (Business as Usual) kept in the future climate. Adaptation scenarios 1–5 are modifications of Scenario 0 (S-0). The 0–5 scenarios’ analysis was based on comparing soil water content and total runoff, sediment yield, actual evapotranspiration. The first adaptation scenario (AS-1) assumes an increase in afforestation on soils from the agricultural suitability complex of soil 6–8 (semi-dry, permanent dry, semi-wet). The second adaptation scenario (AS-2) assumes the creation of a forested buffer for the Bystra River and its tributaries. The third adaptation scenario (AS-3) shows one of the erosion prevention practices, the so-called filter strips. The fourth adaptation scenario (AS-4) assumes the reduction in plowing on arable land. The fifth adaptation scenario (AS-5) involves increasing soil organic carbon to 2%. Simulations revealed that each of the adaptation scenarios 1, 2, 3, 5 does not generally contribute to increasing the water content in soil on BARL (spring crops), CANP (rape), WWHT (winter crops), CRDY (other crops) on arable lands (which together account for over 50% of the catchment area). However, they can contribute to the reduction in sediment yield, total runoff and changes in actual evapotranspiration. The adaptation scenario 4 (AS-4) shows a slight increase in the soil water content on Bystra catchment in the 2041–2050 perspective. Scenario 4 indicated a slight increase in total runoff and a decrease in sediment yield, which in combination with slightly higher water content reflects the protective role of plant residue mulch, lowering the evaporation from the bare soil surface during warm seasons. The no-till adaptation practice had the highest effect in positively affecting water balance at the catchment scale among the adaptation scenarios considered.

1. Introduction

Soil water content is an important component of the hydrological cycle. The formation of water resources in the catchment area is greatly influenced by the amount of precipitation, evapotranspiration, temperature as well as soil properties (water storage capacity, texture, structure), management practices and the existing vegetation [1,2]. The main source of soil water content is precipitation through infiltration and surface runoff [3]. Temperature, on the other hand, influences the evapotranspiration process [4].
There are numerous studies focusing on the calculation of soil water content using the Soil and Water Assessment Tool (SWAT) model [5,6]. These authors used SWAT to simulate soil water content at levels of large catchments (Vistula, Odra). They demonstrated the ability to generate long-term series of soil water content even in the absence of comparative data. On the other hand, for a small catchment located in Poland, one of the few similar studies to the present one in terms of climatic scenarios as well as parameters studied (soil water content, actual evapotranspiration) is a publication concerning the Barycz and Upper Narwia catchments [7].
In the publication on soil water retention and drought risk assessment based on water balance for the area of the Lower Silesian province [1], soil retention parameters were determined: Available Water Capacity (AWC), Wilting Point (WP), Field Capacity (FC) for soil species found in Poland. The retention parameters were determined by expert methods [1].
The aim of the article is to analyze five adaptation scenarios (AS-1, AS-2, AS-3, AS-4, AS-5) in relation to the 2041–2050 climate projections GCMs/RCMs for the RCP 4.5 and RCP 8.5 climate change scenarios described as scenario 0 (S-0) [8], as well as their assessment against the current state of knowledge related to research involving similar adaptation studies. Adaptation scenarios 1–5 are modifications of Scenario 0.
The need for such studies of small catchments (up to a few hundred km2) is due to the small number of studies that would be based on adequate preparation of soil parameters (e.g., retention). Moreover, for the Polish area, there are no studies on adaptation scenarios that would attempt to increase the water content of soil and minimize the adverse effects of climate change (RCP 4.5, RCP 8.5) in future decades.
Among the many hydrological models in use today, the SWAT model, widely used by scientists and developed by the USDA [9,10], was selected for this study because of its ability to predict the impact of practices of land management onto the hydrology and water quality in the catchment area.
Much research is currently being conducted on climate change and the associated unpredictability of extreme weather events. This raises legitimate concerns about the possible emergence of environmental, social and economic threats in the decades to come. These changes may also have an impact on agriculture in Poland [11]. The increase in air temperature, which was observed in recent decades, contributed to the increase in potential evapotranspiration, especially in the last decade 2011–2020. A large increase in potential evapotranspiration and an increase in the variability of this indicator were found [8,12,13]. Recent decades also brought observations of climate change in Poland resulting from the world global warming, changes in precipitation and a number of weather extremes [14,15,16].
These changes also concern the extension of the growing season in Poland. For the years 1971–2000, the length of the growing season was 218 days (from March 31 to November 4) [17]. According to studies on the change in the growing length in Poland [17], the length of the growing season will extend by 18–27 days in the perspective of 2050 compared to the years 1971–2000.
The increase in evapotranspiration, temperature and precipitation in the coming decades will, to a greater or lesser extent, also apply to all European countries [18,19].
According to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change [20], the average temperature of the Earth’s surface will reach 1.5 degrees Celsius in the coming decades above pre-industrial levels. Moreover, in Poland, a projected 10-fold increase in the occurrence of droughts by 2020 [21] is observed in the data; hence, the predictions using climate change models seem to reflect changing climate quite well for Poland [22]. Until recently, climate change adaptation received less attention in Poland than climate change mitigation. The vast majority of national communications have been devoted to climate projections, vulnerability and impacts. However, recently there has been increasing attention to adaptation measures in agriculture, among others [16].
There is a need to look for solutions that will reduce the negative impact of climate change [23], inter alia, the occurrence of weather extremes, including drought [12,24,25] in the coming decades. Climate change adaptation in agriculture is associated with a number of preventive measures (adapting crops to changing thermal and water conditions). These include changes in adaptation practices and the introduction of new varieties. Protecting the soil and its water resources is also extremely important. Soil moisture can be maintained through mulching and water conservation through efficient irrigation and water storage (small retention, filter strips). Soil fertility and its potential for water storage can also be increased by increasing soil organic matter [16].
For the sake of this paper, we chose the Bystra river catchment (South-Eastern Poland) as the study area. In order to check the effectiveness of the designed adaptation solutions, it was necessary to develop boundary conditions that would indicate the reference level [8]. These conditions show the behavior of the hydrosystem of the Bystra catchment in the Business as Usual scenario. It takes into account changes in the hydrological cycle caused solely by climate change while maintaining unchanged conditions of human activity. The described boundary conditions for the 2050 horizon must be based on simulation modeling, which is calibrated on archival data. The appropriate tool for this is the SWAT model.
The article presents a comparison of the results of soil water content (profile 1.5 m) for five adaptation scenarios obtained via a simulation of a calibrated and validated SWAT model [8] for three regional climate models derived from the global EC-EARTH climate model for the years 2041–2050 (S-0). Then, the results of scenario 0 were compared with the results of adaptation scenarios 1–5, which included land use changes and protective measures.
The publication is presented as follows: Section 1 presents the Introduction; Section 2 introduces the methodology and describes the study area. Section 3 describes the results, and Section 4 presents the discussion in terms of results regarding soil water content, total runoff and sediment yield, followed by conclusions in Section 5.

2. Material and Methods

This section is divided in 5 sub-sections. The first one describes the study area; the second and the third describes the SWAT model and SUFI-2 model; the fourth presents climate change scenarios, and finally the fifth presents climate change adaptation scenarios 1–4.

2.1. Characterization of the Study Area

The Bystra catchment area is situated in the north-western part of the Lubelskie Province (Figure 1). The length of the Bystra River is 33 km, and it is the right tributary of the Vistula river. According to the generated SWAT model, the lowest point of the catchment area is 126 m above sea level, and the highest point is 246 m above sea level. The catchment area delineated from a 5 m resolution DEM is 296.6 km2 [8].
The Bystra catchment area is part of the Lublin Upland [26,27,28]. The valley of the Bystra river and its tributaries are strongly carved in a thick loess layer overlaying calcareous bedrock. It consists of numerous valley forms with a constant or episodic tributary. The largest valley with a constant tributary, the Bystra valley, is 35 km long. In the part where the Bystra valley flows into the Vistula, it cuts up to 35 m in rocks and marls [29,30,31].
The upland nature of the Bystra catchment area, consisting mostly of loess soils, with a high slope of the slopes at the mouth of the Vistula, poses a high risk in terms of medium and very strong water and surface erosion [32].
Most of the Bystra catchment area is made of loess up to 20 m. In the deeper layers, there are Quaternary Pleistocene sediments: water-glacial sand and gravel and, at a little deeper level, tilts. On the other hand, there are geodes under the clays. Under the geysers, on the other hand, there are deposits of the Upper Cretaceous: rocks with lime inserts [33].
The study area consists mainly of podzolic and lessivage (49%) soils, which extend mainly in the south-eastern part of the catchment as well as cambisols (47%) in the north-western part. The predominant soil texture in the catchment area is loess (73%) [34,35,36] and silt (18%) [8].
In the Bystra catchment area, arable land (78%) and forests (16%) dominate [8]. The largest part of agricultural land is arable land beyond the reach of irrigation facilities (52%); large areas are also orchards and plantations (11%), complex systems of arable plots (9%) and meadows and pastures (6%) [8].

2.2. Description of SWAT Model and SUFI-2 Model

SWAT was used to model and examine the water balance of the Bystra river catchment area. SWAT is a model [9,10] developed by the USDA Agricultural Research Service [37]. The model operates on assigning one resource to another (physical, chemical, biological) using mathematical formulas that were developed to predict the impact of management practices on water efficiency and agricultural chemistry at the catchment scale [38,39]. We used the QSWAT3 v1.1 model with an interface in Quantum GIS 3.10.13 Coruna [40]. However, the calculations of the SWAT model were performed in the SWAT Editor on 10 December 2012 [41].
The water balance is the fundamental driving force behind all the processes that take place in the catchment area regardless of the choice of the SWAT model analysis. SWAT modeling for the catchment area is carried out in the land phase [42] and in the routing phase [43].
One of the formulas that is used in the SWAT model is the water balance equation:
S W t = S W 0 + i = 1 t ( P d S U R Q E w s e e p G W Q )
where: SWt is the final water content of the soil (mm); SW0 is the initial water content of the soil (mm); t is the time in days; Pd is precipitation (mm); SURQ is surface runoff (mm); E is evapotranspiration (mm); wseep is the amount of water entering the wad zone from the soil profile (mm); GWQ is the groundwater flow (mm) [10].
Calibration and validation in the SWAT-CUP program are used to adjust the SWAT model to real conditions in the catchment area. The commonly used example of calibration is stream flow, which includes water balance processes. The calibration process is used to adjust the relevant parameters so that the simulated results are consistent with the observational data. Validation involves running the model using the parameters that were used during the calibration process. The purpose is to compare simulated results with observed data that were not used in calibration [44,45,46]. The SWAT-CUP program is used to analyze the uncertainty and sensitivity of the model [44,45] using the SUFI-2 algorithm, also used in small catchments [44,47,48].

2.3. Application of SWAT and SUFI-2

To simulate the water balance in the SWAT model, data were obtained from many sources (Table 1), which were used to build the SWAT model.
Table 1. Input data used in SWAT model (own study).
Table 1. Input data used in SWAT model (own study).
Data TypeDescriptionInformationSource
Digital Elevation ModelWatershed delineationRaster, 5 m-resolutionCentral Geodetic and Cartographic Documentation Center [49]
HydrographicSite hydrographic data (e.g., rivers, lakes, partial catchments); (reference scale 1:50.000)Shapefile Computer Map of the Polish Hydrological Department with descriptions [50]
Land useLand-use classification (r.s. 1:100.000)Shapefile Corine Land Cover [51]
OrthophotomapHigh resolution orthophotomapWMSGeoportal [52]
Open Street Map Open Street Map dataShapefileOpen Street Map [53]
Soil typeDigital maps of soil and agriculture in digital form (scale 1: 25,000 and 1: 100,000)ShapefileInstitute of Soil Science and Plant Cultivation in Pulawy [54,55]
GeologicalGeological data describing lithologyShapefilePolish Geological Institute in the form of the Detailed Geological Map of Poland [33]
WeatherPrecipitation (mm), temperature (°C), wind speed (m/s), humidity, solar total radiation (MJ/m2)DailyInstitute of Soil Science and Plant Cultivation in Pulawy and Institute of Meteorology and Water Management [56]
StreamflowCalibration and validationMonthlyInstitute of Soil Science and Plant Cultivation in Pulawy
Sewage treatment plantsAverage daily water loading (m3/day)DailyNational Program of Municipal Wastewater Treatment [57]
The SWAT model generated for this study consists of 31 generated partial catchments (Figure 1) [8]. The soil map was developed on the basis of digital soil and agricultural maps (scale 1:25,000 and 1:100,000) and geological data describing lithology. Descriptive soil data were collected within the statutory research projects of IUNG-PIB. Available water capacity and wilting point values were obtained from the study “Assessment of water retention in soil and the risk of drought based on the water balance for the Lower Silesian Voivodeship”, which was developed in 2013 by the employees of the Department of Soil Science, Erosion and Land Protection of IUNG-PIB in Pulawy [1].
The land use map was developed on the basis of Corine Land Cover maps with additional vectorization of land cover and land use using an orthophoto-map and Open Street Map data.
Based on the generated maps of soils, lands and slopes, 484 HRU (Hydrological Response Units) areas were created. When creating HRU areas, the land cover class of agricultural areas beyond the reach of irrigation CRDY was additionally separated with WWHT winter crops (43%), BARL spring crops (31%), CANP rape (14%) and other CRDY (12%) [58]. APPL apple orchards were separated from the land use class of ORCD [58]. On the other hand, forests were divided into coniferous FRSE forests (49%), deciduous FRSD forests (13%) and mixed FRST forests (38%) [59].
After generating HRU areas, the following meteorological data were used in the SWAT model: daily precipitation totals (mm); daily minimum and maximum air temperature (°C); average daily wind speed (m/s); average daily relative humidity; daily sums of total solar radiation (MJ/m2) (Table 2) [8].
In the SWAT model, the parameters related to the point discharge of sewage, as well as for water bodies located outside the river network, for water bodies, rivers, and parameters for planned non-irrigated arable land management operations (WWHT, BARL, CANP, CRDY) were supplemented and corrected. The current value of CO2 concentration was also entered.
In the next stage, the SWAT model simulation was run for the period of 2010–2017 in a monthly step, with a five-year model start-up period.
Then, calibration and validation of the obtained SWAT model for the Bystra catchment area [8] was performed using the SWAT-CUP program. To obtain a more accurate coverage of the model with reality, the average monthly flow velocities (m3/s) obtained under the statutory projects of IUNG-PIB, obtained near the mouth of the Bystra River to Vistula for 2010–2014 (calibration) and 2015–2017 (validation), were used. A five-year warm-up period was used. Calibration and validation were performed in a monthly increment. This resulted in parameter ranges that fell within the ranges of calibration and validation accuracy [44,60,61]. The NSE coefficients (calibration: 0.58; validation: 0.70) and R2 (calibration: 0.60; validation: 0.71) for calibration and validation [8] were within the satisfactory ranges [60].
The results concerning the value of potential evapotranspiration were also analyzed with the results of the statutory service of IUNG-PIB implemented under the project Agricultural Drought Monitoring System [62]. It was found that the SWAT model for the Bystra catchment area accurately reflects the potential evapotranspiration in the study area.
Additionally, the results concerning the soil water content were compared with the available values of water capacity and the wilting point, which were obtained from the study prepared in 2013 by the employees of the Department of Soil Science, Erosion and Land Protection, IUNG-PIB in Pulawy [1].

2.4. Climate Change Scenarios

The daily grid climate data used in the SWAT model were prepared and tested in the recent paper on SWAT model calibration in the Bystra catchment [8]. Three RCM (Regional Climate Models)—RACMO22E, HIRHAM5 and RCA4—were selected for further study. They were selected to cover the range of the available two climate scenarios RCP (Representative Concentration Pathways) in terms of temperature increase and precipitation—RCP 4.5, RCP 8.5 (Table 3)—reflecting extreme and average variants of climate change, hence covering the widest range of uncertainty about possible scenarios (three RCM × two RCP). Most of the data were obtained at a spatial resolution of 0.11 degrees from the EURO-CORDEX database for the years 1951–2050 (widely available via the ESGF—Earth System Grid Federation, https://esgf-data.dkrz.de/search/cordex-dkrz for Europe) (accessed on 3 March 2021) [18,63].
Climate scenario daily meteorological derivatives (minimum and maximum daily air temperature, daily precipitation, solar radiation, daily average wind speed, relative humidity) are based on the RCM for two RCPs (three RCM x two RCP). The RCMs are powered by one GCM (General Circulation Model): EC-EARTH. The RCP corresponds to the radiative forcing values in 2100 compared to pre-industrial values of +4.5 W m−2 (RCP4.5) while RCP8.5 to + 8.5 W m−2 (RCP8.5) [18,64,65] (Table 3). Table 3 also presents the boundary values of changes in the characteristics of selected models for the period 2021–2050 in relation to the period up to the base period 1971–2000.
Climate projections that were used in the SWAT model were extracted from grid cells that correspond to weather stations’ location. Air temperature and precipitation data were additionally corrected by the SMHI (Swedish Meteorological and Hydrological Institute) using the DBS (Distribution-Based Scaling) method [48] and regional MESAN reanalysis (MESoscale Analysis) for the 1989–2010 dataset [66]. The data used were taken in a rotated polar grid. Therefore, we used bilinear interpolation to remap the dataset to a common latitude/longitude grid. CDO (Climate Data Operators) software [67] was used for this purpose.
For the analysis of the climate projections (RCP 4.5.1, RCP 8.5.1, RCP 4.5.2, RCP 8.5.2, RCP 4.5.3 and RCP 8.5.3), one iteration in SWAT-CUP was used for the set of the best calibration parameters for the years 2021–2050 in the prepared scenarios (Table 3) [8]. In the RCP 4.5 and RCP 8.5 scenarios, CO2 concentrations were changed for the periods 2021–2030, 2031–2040 and 2041–2050, developed by the Potsdam Institute for Climate Impact Research [68,69].

2.5. Climate Change Adaptation Scenarios 1–5

For the main purpose of this article, 5 scenarios for the adaptation of agriculture to climate change were prepared, which assume changes in land use (adaptation scenario 1 and 2) and protective measures (adaptation scenario 3, 4, 5) in the area of the Bystra catchment. The first adaptation scenario (AS-1) assumes an increase in afforestation on soils from the agricultural usefulness complex of soils 6 (temporarily too dry), 7 (permanently too dry) and 8 (temporarily too wet). The second adaptation scenario (AS-2) assumes the creation of a forested buffer for the Bystra River and its tributaries. The third adaptation scenario (AS-3) shows one of the erosion prevention practices at the riverbed, the so-called filter strips. The fourth adaptation scenario (AS-4) assumes the reduction in plowing on agricultural land. The fifth adaptation scenario (AS-5) involves increasing soil organic carbon to 2%. Adaptation scenarios are aimed at checking the possibility of increasing the soil water content in the 2041–2050 perspective. In doing so, the effects of adaptation scenarios on total runoff, sediment yield and actual evapotranspiration were also checked.
In the zero scenario (S-0), the Bystra catchment area is dominated by agricultural land (78%) and forests (16%). The largest part of agricultural land is arable land beyond the range of irrigation facilities (52%); a large area is also orchards and plantations (11%), complex systems of cultivating plots (9%) and meadows and pastures (6%) (Table 4). For adaptation scenarios 1 (AS-1) and 2 (AS-2), there will be changes in land use compared to scenario 0 (S-0), which are described later. In contrast, adaptation scenarios 3 (AS-3), 4 (AS-4) and 5 (AS-5) remain unchanged in terms of changes in land use.
In the first adaptation scenario (AS-1), the land use on all soils of complexes (representing soil habitats in Polish soil-agricultural mapping)—6 (semi-dry), 7 (permanently dry) and 8 (semi-wet) (6Bw-pgl.ps, 7Bw-ps, 8A-l)—was changed to mixed forest. The soils where the land use was changed are described in more detail in Table 1 of the publication on the water balance of the Bystra catchment [8]. Replacement of the above-mentioned soils is made through delineating the ranges of these soils on the land use maps and changing the attributes to mixed forests. After this change, afforestation in the Bystra catchment area increased by 3.31% (Table 4).
In the second adaptation scenario (AS-2), a forested buffer strip 80 m wide along the bank of the Bystra River was created and a smaller buffer strip 50 m wide for its tributaries [70,71,72]. The creation of buffer zones by rivers consisted of deleting the ranges of buffer zones on the land use maps and changing the attributes to mixed forests. The afforestation area compared to the zero scenario increased by 1.03% (Table 4).
In the third adaptation scenario (AS-3), filter strips were used, which are one of the protective measures used to drain water slowly from the field, thanks to which larger particles, including soil and organic material, may be deposited [73].
Filter strips [9,74] are areas covered with vegetation that are located between surface water bodies (rivers, ponds, lakes) and arable land, pastures and forests. They are generally found in areas where runoff leaves the field to filter sediment, organic material, nutrients and chemicals from the runoff. Filter strips are also known as vegetative filters or buffer strips. Due to the retention of sediment and the establishment of vegetation, nutrients can be absorbed into the sediment that settles and remain in the field landscape, making it possible for plants to take it up [73].
A protective treatment is also tillage without plowing [73], which is the fourth adaptation scenario (AS-4).
Plowing is defined as the mechanical disturbance of soil for crop production that has a significant impact on soil properties such as soil water behavior, soil temperature, infiltration and evapotranspiration [75]. In the long term, tillage can lead to soil degradation [76]. An alternative to traditional plowing is protective treatments (tillage without plowing, minimal mechanical disturbance of the soil) which consist of maintaining the surface soil cover by retaining crop residues. Retention of harvest residues protects the soil from direct exposure to raindrops and sunlight, while minimal soil disturbance improves soil biological activity and air and water movement in the soil [75].
No plowing cultivation was implemented in WWHT, BARL, CANP and CRDY arable land and simulated in SWAT.
In the fifth adaptation scenario (AS-5), the soil organic carbon content was increased from 1% to 2%. The original soil organic carbon values were studied as part of IUNG-PIB statutory research [8]. Soils in Poland are characterized by low soil organic carbon content. According to the European Soil Bureau (ESB), an organic carbon content of about 1% (Bystra catchment area) is a very low or low value [77]. The decrease in organic matter in soils and the associated decrease in organic carbon content result in increased CO2 emissions (exacerbating the greenhouse effect). The opposite situation, i.e., sequestration of CO2 in the soil, causes carbon to bind to soil organic matter for a longer period of time. Particularly large amounts of carbon are stored in peats, organic soils and organic-mineral soils [77].

3. Results

Section 3.1 describes the analysis of soil water content in S-0 for the period 2041–2050.
For the 10-year period (2041–2050), Table 5 presents a comparison of the seasonal soil water content in the Bystra catchment for each climate projection GCMs/RCMs under the RCP 4.5 and RCP 8.5 climate scenarios.
For the 10-year period (2041–2050), Figure 2 shows the average soil water content (1.5 m) for each season of DJF, MAM, JJA, SON for the GCMs/RCMs climate projections under the RCP 4.5 and RCP 8.5 climate scenarios, while Figure 3 shows the spatial comparison of average soil water content in 31 sub-catchments for the SWAT simulation period 2010–2017 and 2041–2050 for the GCMs/RCMs climate projections under the RCP 4.5 and RCP 8.5 climate scenarios.
Section 3.2 describes the climate change AS-1, AS-2, AS-3, AS-4, AS-5 analysis for the period 2041–2050.
For the period 2041–2050, Table 6 presents a comparison of AS-1, AS-2, AS-3, AS-4, AS-5 with respect to S-0 for seasonal soil water content in the Bystra catchment for the RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, and RCP 8.5.3 projections.
Table 7 presents a comparison of total runoff by season for S-0 and AS-1, AS-2, AS-3, AS-4, AS-5 for the period 2041–2050 in the Bystra catchment. Next, Table 8 compares sediment yields by season for S-0 and AS-1, AS-2, AS-3, AS-4, AS-5 for 2041–2050 in the Bystra catchment. In turn, Table 9 compares actual evapotranspiration by season for S-0 and AS-1, AS-2, AS-3, AS-4, AS-5 for the years 2041–2050 in the Bystra catchment.

3.1. Analysis of Soil Water Content in Zero Scenario for 2021–2050

This section compares the obtained seasonal average soil water content results for 2010–2017 (SWAT model) with the results for 2041–2050 for individual climate change projections (RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3).
Regardless of the individual climate change projections evaluated, the seasonal average soil water content for the Bystra catchment is projected to decrease between 2041 and 2050 for most seasons compared to 2010–2017 (Table 5).
Lower soil water content will be especially evident for RCP 4.5.1 (MAM, JJA, SON) and RCP 4.5.2 (JJA) where the value of average soil water content may be lower by up to 5.8% compared to the 2010–2017 simulation period. Lower values in the MAM and JJA seasons, especially for the RCP 4.5.1, RCP 4.5.2, RCP 4.5.3 projections, may affect plant growth during the growing season. However, higher soil water content (1.4% higher) was found for the RCP 8.5.1 (JJA), RCP 8.5.2 and RCP 8.5.3 (SON) projections.
Regardless of the regional climate model, the seasonal average soil water content will be lower for climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3 compared to climate projections RCP 8.5.1, RCP 8.5.2, RCP 8.5.3. This is particularly evident when comparing the average annual soil water content results, where, for RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, the average annual soil water content results (2041–2050) are lower between 1.8% and 4.7%, while, for RCP 8.5.1, RCP 8.5.2, RCP 8.5.3, these average annual results are lower between 0.3% and 0.6% compared to the SWAT 2010–2017 model.
The average soil water content by season for 2041–2050 and the SWAT model 2010–2017 is shown in Figure 2 (Figure 2). It shows that the average soil water content decreases throughout the year. The highest soil water values are reached during the winter season of DJF. On the other hand, in spring (MAM), during the growing season period, the average soil water content decreases, maintaining the lowest values in summer (JJA). In autumn (SON), the soil water content increases.
Analyzing the spatial distribution of changes in the average water content in soil in 31 sub-catchments for the simulation period in 2010–2017 in relation to the period 2041–2050 (Figure 3) in the climate forecasts RCP 4.5.1 and RCP 4.5.3, the average water content in the soil will decrease by a few percent points in the Northwest region for most of the projections. In the projections RCP 4.5.1, RCP 4.5.2 and RCP 4.53, a reduced water content in the soil will occur throughout the catchment area, while in the projections RCP 8.5.1, RCP 8.5.2 and RCP 8.5.3, the changes will be small.

3.2. Climate Change Adaptation Scenarios Analysis 1–5 for 2041–2050

The results presented in Section 3.1 indicate a decrease in soil water content in most seasons during the period 2041–2050 (Table 5, Figure 2). To counteract the negative effects of changes in soil water content, five adaptation scenarios (AS-1, AS-2, AS-3, AS-4, AS-5) were prepared and tested. They were designed to maintain or increase soil water content. The analysis covers the period 2041–2050. Additionally, the impact of adaptation scenarios on total runoff, sediment yield and actual evapotranspiration was compared.
AS-1 of increasing forested areas on soils of complex 6, 7, 8 compared to S-0 for all projections shows a decrease in soil water content for all seasons in the Bystra catchment (Table 6). The soil water content decreases from 4.0% to 4.8% for all seasons.
AS-2, which assumes a forested buffer zone near the Bystra River, shows a slight decrease in soil water content between 2041 and 2050 (Table 6).
AS-3, establishing filter strips, shows no change in soil water content (Table 6).
In AS-4, the application of plowing on arable land—BARL, CANP, CRDY, WWHT—was eliminated. This treatment showed a slight increase in soil water content. The increase ranged from 0% to 0.4%. The largest increases occurred in the JJA and SON seasons (Table 6).
AS-5 increased soil organic carbon to 2%. This treatment showed a slight decrease in soil water content. The decrease in soil water content ranged from 0% to 1.4% (Table 6).
Regardless of the GCMs/RMCs and the RCPs evaluated, the results are the same. This means that AS-1 is associated with a greater decrease in soil water content compared to S-0. AS-2 and AS-5 are associated with a decrease of a smaller magnitude compared to AS-1. AS-3 does not predict any significant change in soil water content. In contrast, AS-4 is associated with a small increase in soil water content.
Regardless of the regional climate model, the seasonal average soil water content will be lower under the RCP 4.5 climate change scenario compared to the RCP 8.5 climate change scenario. This is described in more detail in Section 3.1.
Differences in annual average soil water content between AS-2, AS-3, AS-4, AS-5 and S-0 are small. However, for AS-1, the annual average soil water content varies between 296 and 311 mm. In contrast, for S-0, the average annual soil water content is 310–325 mm (Table 6).
AS-1 and AS-2 show a slight decrease for most seasons of total runoff for 2041–2050 compared to S-0 in all climate projections. Changes in total runoff range from 1% (decrease) to 0.3% (increase) (Table 7).
Total runoff in AS-3 did not change (Table 7).
AS-4 shows an increase in total runoff for all seasons in all projections. The increase ranges from 1.3% to 4.7% (Table 7). For climate projection RCP 4.5.1, the increase in total runoff stands out from the other projections in all seasons (above 3%).
In contrast, AS-5 shows a decrease in total runoff for all seasons across all projections. The decrease ranges from 1.0% to 3.9% (Table 7). For climate projections RCP 4.5.1 and RCP 4.5.2, the decrease in total runoff stands out from the other projections in all seasons (above 2%) (Table 7).
Moreover, for total runoff regardless of the GCMs/RMCs and RCPs evaluated, the results are the same. AS-1 and AS-2 have smaller total runoff compared to adaptation S-0. AS-5 has an even smaller total runoff compared to AS-1 and AS-2.
In AS-3, the total runoff does not change. In contrast, AS-4 shows an increase in total runoff compared to all adaptation scenarios.
Regardless of the regional climate model, the average seasonal total runoff will be lower for the RCP 4.5 climate change scenario compared to the RCP 8.5 scenario [8].
Table 8 presents the seasonal sediment yield data (Table 8). AS-1, AS-2, AS-3, AS-4, AS-5 were compared to S-0 for the climate projections. For most adaptation scenarios, there is a reduction in sediment yield from 0% to as much as 73%. The smallest, slight decreases in sediment yield occur in AS-2 compared to S-0. Slightly larger decreases compared to S-0 and AS- 2 occur in AS-1 and AS-5. Large decreases in sediment yield occur in AS-4 (ranging from 6% to 49%). However, the largest occur for AS-3 (over 70%).
For sediment yield, regardless of the GCMs/RMCs and RCPs evaluated, the results are also the same (Table 8).
Regardless of the regional climate model, the seasonal sediment yield will be lower under the RCP 4.5 climate change scenario compared to the RCP 8.5 scenario. Differences in annual sum sediment yields range from 0.54–0.57 t/ha for RCP 4.5 to 0.72–1.00 t/ha for RCP 8.5 in S-0 (Table 8). For AS-4, the annual sum ranges from 0.39–0.40 t/ha for RCP 4.5 to 0.50–0.71 t/ha for RCP 8.5, while for AS-3, the annual sum ranges from 0.15–016 t/ha for RCP 4.5 to 0.20–0.28 t/ha for RCP 8.5.
Table 9 presents data on seasonal actual evapotranspiration (Table 9). The highest evapotranspiration values occur during the MAM and JJA seasons.
AS-1, AS-2, AS-3, AS-4, AS-5 were compared to S-0 for all climate projections. AS-1 shows a decrease in actual evapotranspiration from 1.7% to 1.9% for the MAM season. In contrast, there is an increase between 1.5% and 2.0% for the JJA season.
AS-2 shows little change in actual evapotranspiration (Table 9).
The actual evapotranspiration in AS-3 remains unchanged compared to S-0 (Table 9).
In AS-4, for the MAM and JJA seasons, actual evapotranspiration varies from 1.5% (decrease) to 0.1% (increase) (Table 9) compared to S-0. However, large decreases occur for the SON season (from 3.2% to 4.2%).
In contrast, AS-5 has increases in actual evapotranspiration of 0.7% to 1.7% for the MAM and JJA seasons compared to S-0.
For actual evapotranspiration, regardless of the GCMs/RMCs and RCPs evaluated, the results are the same (Table 9).
Regardless of the regional climate model, seasonal actual evapotranspiration will be similar under the RCP 4.5 climate change scenario compared to the RCP 8.5 scenario [8].
For AS-1, AS-2, AS-3, the annual sum of actual evapotranspiration changes little. However, for AS-4, the annual sum of actual evapotranspiration increases from 1.7% to 3.8% compared to S-0. In contrast, for AS-5, the annual sum decreases from 1.2% to 3.4% (Table 9).
Table 10 shows the percentage sets of changes in soil water content, sediment productivity, total runoff, and actual evapotranspiration under AS-1, AS-2, AS-3, AS-4, AS-5 with respect to S-0 (Table 10). The table was created based on supplementary Material: Figures S1 and S2, for averages of three GCMs/RCMs combinations under two RCP climate change scenarios (RCP 4.5, RCP 8.5).

4. Discussion

The results concerning the water content in the soil were compared with the available values of water capacity and the wilting point obtained from the study “Assessment of water retention in soil and the risk of drought based on the water balance of the Lower Silesia Voivodshi”, developed in 2013 by the employees of the Department of Soil Science, Erosion and Land Protection, IUNG-PIB in Pulawy [1]. Based on the above-mentioned study, we prepared data on soils in the catchment area of the Bystra River. For a 1.5 m soil profile, the results of the above-mentioned studies are consistent with this publication.
The lowest water content in soil occurs in the summer (JJA), while the highest occurs in the winter (DJF) (Figure 2). For 2041–2050, the largest decreases in soil water content are associated with GCMs/RCMs for RCP 4.5, while small changes occur for RCP 8.5.
The analyzed adaptation scenarios present different results of the influence on the water content in the soil. AS-1 for an increase in forest area on soils of the complex 6, 7, 8 compared to S-0 for all projections shows a reduction in soil water content for all seasons across the entire Bystra catchment (Table 6 and Table 10, Figures S1 and S2 in Supplementary Material). The same is true for total runoff. Again, for most seasons, there is a reduction in total runoff (all projections) compared to S-0 (Table 7 and Table 10, Figures S1 and S2 in Supplementary Material). Sediment yields for all seasons also decrease (Table 8 and Table 10, Figures S1 and S2 in Supplementary Material). In contrast, actual evapotranspiration shows a decrease in the MAM season and an increase in the JJA season (Table 9 and Table 10, Figures S1 and S2 in Supplementary Material).
Forests play an important role in absorbing CO2, which is an important factor in reducing the adverse effects of climate change [78]. In addition to absorbing CO2, forest ecosystems can counteract soil erosion and drainage. Within forests, there may be small retention reservoirs, increasing the areas’ abundance of water. Forest ecosystems play very important natural, social and productive functions [79]. The results indicate that increasing afforested area in the Bystra catchment has to go beyond the scheme of using soil complexes less favorable for agricultural production, and the areas should be picked with care, focusing on locating forested areas close to catchment borders, so they can slow runoff and help accumulate water at its highest point from the river bed [80].
The large-scale research aimed at estimating the amount of tree stand in the world shows that there are currently 46% fewer trees than before the advent of human civilization [81]. Climate change may affect the condition of forest areas [19] manifested in extreme weather phenomena that begin to lose their anomaly status (hurricanes, droughts). Moreover, the species status of plants and trees may not be flexible enough to adapt to changing climate components (temperature, precipitation, etc.) [82]. Forests therefore should be probably re-designed to cope with changing biotopes. For many years, many concepts regarding forest formation in relation to a changing climate have been considered. These plans are based on the development of actions to reduce the effects of unfavorable phenomena which are occurring now and which may intensify in the future. Another concept will be activities aimed at adapting forest ecosystems to all current and future threats [82].
A program of increasing forest cover is implemented in Poland [83]. According to the report on the condition of forests in Poland in 2020 [84], the level of forest cover in 2020 amounted to 29.6% of the total area of the country. After 2050, the forest cover in Poland is expected to be 33%. The program assumes afforestation of land of low agricultural suitability [85], reflected in AS-1 of this study.
Research using afforestation scenarios was carried out on four sites in Bolivia and Ecuador [86]. They show that the water content in the soil and the total runoff decreased to a varying degree after the application of the forest ecosystem. AS-1 and AS-2 also show a reduction in soil water content (Table 6 and Table 10, Figures S1 and S2 in Supplementary Material) and a slight reduction in total runoff (Table 7 and Table 10, Figures S1 and S2 in Supplementary Material). Sediment yield also decreased (Table 8 and Table 10, Figures S1 and S2 in Supplementary Material). The decrease in soil water content for AS-1 and AS-2 in the Bystra catchment may be caused by increased water uptake by the root system of forest vegetation species.
The afforestation scenario has the potential for further research, in which it is possible to design an appropriate location of forest ecosystems in the Bystra river catchment area, relying not only on the afforestation of soils of complex 6, 7 and 8, but also good tree planting practices in rural areas [87], the use of forested embankment fortifications (also preventing erosion) [72], which would counteract the unfavorable agro-forest checkerboard [88]. The unfavorable location of forest ecosystems near cultivated fields may result in a reduction in the yield of agricultural plants [89,90]. When designing afforestation, one should also take into account the adaptation possibilities of stands to new climatic conditions [82].
Increasing forest cover from 16.34% (S-0) to 19.65% (AS-1) or to 17.37% (AS-2) (Table 4) according to Lambo’s forest cover index [91] allows for increased forest retention capacity that, among other things, counteracts the effects of flooding [92]. In addition to increasing forest cover, equally important is the location of forested areas within the catchment area which has a significant impact on runoff [93].
A buffer zone with a well-developed tree stand, located directly next to watercourses, can prevent the runoff of nutrients and suspensions from agricultural land, contribute to the strengthening of banks and prevent lateral erosion [72, 87, 94]. A marsh zone forming a belt of wetland and rush vegetation, flooded or boggy for most of the year or all the time, can also be a buffer. Such a zone with well-developed vegetation contributes to the retention of a significant amount of nitrogen and phosphorus from the catchment area, preventing eutrophication of waters [72,95].
AS-3, for the creation of filter strips in a planned management operation on BARL, CANP, CRDY, WWHT arable land, shows no changes in soil water content, total runoff or actual evapotranspiration (Table 6, Table 7, Table 9 and Table 10, Figures S1 and S2 in Supplementary Material). On the other hand, the filter strips effectively reduce the sediment yield (t/ha) (Table 8 and Table 10, Figures S1 and S2 in Supplementary Material). Similar results were obtained in the article describing the use of the filter strips in various scenarios on the example of the catchment area in Thailand [96], where, as a result of their use, the sediment yield was significantly reduced.
Adaptation scenarios involving increasing forest cover, creating buffers next to rivers and creating filter strips can help reduce erosion risk in the 2050 climate horizon in the Bystra catchment by reducing total runoff and decreasing sediment yield.
AS-4, for the cessation of plowing on BARL, CANP, CRDY, WWHT arable land, shows a slight increase (especially in the JJA and SON season) in soil water content (Table 6 and Table 10, Figures S1 and S2 in Supplementary Material). The elimination of plowing also shows a significant reduction in sediment yield (t/ha) (Table 8 and Table 10, Figures S1 and S2 in Supplementary Material). This may have the effect of reducing soil erosion. However, the total runoff increased, which is induced mainly by the reduction in actual evapotranspiration, especially limited evaporation form the soil surface covered by plant residue mulch (Table 7, Table 9 and Table 10, Figures S1 and S2 in Supplementary Material). Observations by Wawer and Kozyra [97] confirm the prominent role of mulching in preserving soil water by covering the surface of the soil in warm periods.
The discontinuation of plowing is the subject of many articles as well as studies that mention as benefits the reduction in soil erosion, the reduction in surface and subsurface runoff, the reduction in sediment yield, nitrogen yield and phosphorus yield, the increase in soil water content, etc. [76,98,99,100,101], which are supported by numerous studies. The abandonment of plowing in the catchment areas in the climate of 2050 also shows a reduction in the sediment yield. On the other hand, the water content in the soil increases. This provides the grounds that new agricultural practices in the coming decades may prevent the negative impact of watershed water deficits from occurring.
Agriculture is closely related to the prevailing climatic conditions, but it also has a large impact on them. The risk of an increase in the frequency of unfavorable climatic conditions in agriculture may result in yield variability from year to year. The reduced amount of water in the soil during plant growth, illustrated in the climate change scenarios (Table 5, Figure 2), will become more frequent and more severe. Other threats will also include droughts, heavy precipitation, erosion [80], floods, landslides and strong winds [102].
AS-5, increasing soil organic carbon to 2%, shows reductions in soil water content, total runoff and sediment yield (Table 6, Table 7 and Table 8, Figures S1 and S2 in Supplementary Materials). However, actual evapotranspiration increases (Table 9 and Table 10, Figures S1 and S2 in Supplementary Materials). In a paper on soil organic carbon changes and their response to climate warming and soil water content changes, a study of the Jinghe catchment in China was described [103]. The study showed that temperature and precipitation will increase by the end of the 21st century under three scenarios—RCP 2.6, RCP 4.5, RCP 8.5—and consequently soil water content will also increase, while organic carbon content will decrease, depending on the climate change scenario. The study also showed that there is a threshold in soil water content that can mediate the loss of soil organic carbon (when the change in soil water content was lower than the threshold, higher content accelerated the loss of organic carbon, while when the change in soil water content was higher than the threshold, higher content reduced the loss of soil organic carbon) [103]. The mechanism for the decrease in soil organic carbon (despite increased soil water content) due to a warming climate in the future is not fully known [103]. Global studies have found a link between faster CO2 increases in warmer years with less water availability. This demonstrates the importance of warming on the decomposition of soil organic carbon [104]. There are studies in pols on the effect of soil organic matter on soil water management [105]. According to some estimates in the article, increasing soil organic matter by 0.01% increases the amount of organic matter by 480 kg (from 1 hectare of arable soil layer). This corresponds to 278 kg of organic carbon. On the scale of the national area (Poland), this means the sequestration of 11 million tons of CO2 from the entire arable land area of Poland. This represents more than 3% of the total greenhouse gas emissions from the Polish area [77].
A convenient tool for carrying out beneficial changes (afforestation, retention reservoirs, irrigation) in terms of water retention in the landscape is land consolidation on an extended scope [80, 106, 107]. Several agricultural research centers in Poland deal with the issues of recomposing the rural landscape, including IUNG-PIB in Pulawy. At IUNG-PIB, a broader consolidation formula, called the Composite Development of Rural Areas CDRA [106], was developed, covering extended land consolidation, rural area management and rural development, which are included in addition to classic land consolidation works meant as the transformation of land, water drainage, water supply to farms aimed at improving the conditions for agricultural production on farms [106]. The comprehensive, holistic land consolidation approach remains the most effective way of introducing a wide range of changes in the agricultural landscape, also focusing on water management [106]. Based upon the outcomes of this study, the team plans to simulate a scenario of a fully designed land consolidation with the CDRA scheme as one of the options towards a better holistic water management in rural landscapes.
One of the more recent publications describing the methods of managing water resources and thus counteracting climate change in agriculture for the Polish area is the Code of Good Water Practices in Agriculture, which was commissioned by the Ministry of Agriculture and Rural Development [22]. The Code describes various sustainable and solidarity-based water management practices that can be successfully applied to agriculture in the coming decades in response to an increasing scarcity of water resources. We plan to model the effects of introducing the practices covered by the Code in future studies.

5. Conclusions

AS-1, AS-2 and AS-5 did not increase the water content of the soil. However, they can help to reduce sediment yield and total runoff. AS-1 and AS-2 have potential for further research using the SWAT model. The research would be aimed at adopting an appropriate strategy for spreading the location of afforestation in the catchment to reduce the adverse effects of climate change. Soil organic carbon sequestration (AS-5) also has potential for further research due to the reduction in negative effects of climate change.
The filter strips in AS-3 contributed to a reduction in sediment yield. Soil water content, total runoff and actual evapotranspiration remained unchanged. The lack of change may be due to suboptimal discretization of the filter strips in the SWAT input files. Further research on this issue will be conducted.
Practices for reducing or eliminating water shortages in soil can be those presented in AS-4 for no-tillage cultivation. Removal of plowing may also contribute to the reduction in sediment yield (t/ha). This may have the effect of reducing soil erosion. However, the positive influence on soil moisture contents throughout the season using the no-till simulation indicated an increase in runoff, which is mainly caused by limiting evaporation from bare soil covered by the mulch of crop residues.
The obtained results cover 150 cm of the soil layer as described by the Polish soil-agricultural map, which does fully reflect the conditions for plants, especially during sawing and in early stages of growth. Further research has to be conducted on discretizing soil hydrology dynamics in the SWAT input configuration to take into consideration the plough horizon as a separate hydrological entity to be modeled.
Higher soil water content, higher total runoff and higher sediment yield for the RCP 8.5 climate change scenario compared to the RCP 4.5 climate change scenario may be related to higher precipitation in 2041–2050 (Badora et al., 2022).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w14152288/s1, Figure S1: Summary of changes in soil water content, sediment yield, total runoff and actual evapotranspiration in adaptive scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) compared to scenario 0 (S-0), for averages of three GCM/RCM combinations in the RCP 4.5 climate change scenario. The list covers four seasons (DJF, MAM, JJA, SON) in the Bystra catchment area. The first adaptation scenario assumes the growth of afforestation on soils from the agricultural usefulness complex of soil 6–8 (semi-dry, permanent dry, semi-moist, permanently wet). The second adaptation scenario assumes the creation of a forested buffer for the Bystra River and its tributaries. The third adaptation scenario shows one the erosion prevention practices in the river bed, the so-called filter strips. The fourth adaptation scenario assumes the reduction of plowing on agricultural land. The fifth adaptation scenario assumes an increase in soil organic carbon content to 2%. Adaptation scenarios 1–5 are modifications of scenario 0. Scenario 0 only covers climate change in 2041–2050 (own study); Figure S2: Summary of changes in soil water content, sediment yield, total runoff and actual evapotranspiration in adaptive scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) compared to scenario 0 (S-0), for averages of three GCM/RCM combinations in the RCP 8.5 climate change scenario. The list covers four seasons (DJF, MAM, JJA, SON) in the Bystra catchment area. The first adaptation scenario assumes the growth of afforestation on soils from the agricultural usefulness complex of soil 6–8 (semi-dry, permanent dry, semi-moist, permanently wet). The second adaptation scenario assumes the creation of a forested buffer for the Bystra River and its tributaries. The third adaptation scenario shows one the erosion prevention practices in the river bed, the so-called filter strips. The fourth adaptation scenario assumes the reduction of plowing on agricultural land. The fifth adaptation scenario assumes an increase in soil organic carbon content to 2%. Adaptation scenarios 1–5 are modifications of scenario 0. Scenario 0 only covers climate change in 2041–2050 (own study).

Author Contributions

Conceptualization, D.B. and R.W.; methodology, D.B and R.W.; software, D.B.; validation, D.B., R.W., J.K. and A.N.; formal analysis, B.J. and E.N.; investigation, D.B.; resources, D.B.; data curation, J.K. and A.K.-B.; writing—original draft preparation, D.B.; writing—review and editing, A.N., J.K. and A.K.-B.; visualization, D.B.; supervision, R.W.; project administration, R.W.; funding acquisition, R.W. All authors have read and agreed to the published version of the manuscript.

Funding

Research and ACP was funded by Polish Ministry of Agriculture and Rural Development, DC2.0/2021 Programme and B-Ferst project funded by the European commission.
Water 14 02288 i001

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IUNG-PIB. Ocena Retencji Wody w Glebie i Zagrożenia Suszą w Oparciu o Bilans Wodny Dla Obszaru Województwa Dolnośląskiego, Zakład Gleboznawstwa Erozji i Ochrony Gruntów; IUNG-PIB: Puławy, Poland, 2013. [Google Scholar]
  2. Havrylenko, S.B.; Bodoque, J.M.; Srinivasan, R.; Zucarelli, G.V.; Mercuri, P. Assessment of the soil water content in the Pampas region using SWAT. Catena 2016, 137, 298–309. [Google Scholar] [CrossRef]
  3. Somorowska, U. Soil water storage in Poland over the years 2000-2015 in response to precipitation variability as retrieved from GLDAS Noah simulations. Geogr. Pol. 2017, 90, 53–64. [Google Scholar] [CrossRef] [Green Version]
  4. Wang, Y.; Yang, J.; Chen, Y.; Wang, A.; De Maeyer, P. The Spatiotemporal Response of Soil Moisture to Precipitation and Temperature Changes in an Arid Region, China. Remote Sens. 2018, 10, 468. [Google Scholar] [CrossRef] [Green Version]
  5. Abbaspour, K.C.; Rouholahnejad, E.; Vaghefi, S.; Srinvasan, R.; Yang, H.; Kløve, B. A continental-scale hydrology and water quality model for Europe: Calibration and uncertainty of a high-resolution large-scale SWAT model. J. Hydrol. 2015, 524, 733–752. [Google Scholar] [CrossRef] [Green Version]
  6. Piniewski, M.; Szcześniak, M.; Kardel, I.; Berezowski, T.; Okruszko, T.; Srinivasan, R.; Vikhamar Schuler, D.; Kundzewicz, Z.W. Hydrological modelling of the Vistula and Odra river basins using SWAT. Hydrol. Sci. J. 2017, 62, 1266–1289. [Google Scholar] [CrossRef]
  7. Marcinkowski, P.; Piniewski, M.; Kardel, I.; Szcześniak, M.; Benestad, R.; Srinivasan, R.; Ignar, S.; Okruszko, T. Effect of Climate Change on Hydrology, Sediment and Nutrient Losses in Two Lowland Catchments in Poland. Water 2017, 9, 156. [Google Scholar] [CrossRef] [Green Version]
  8. Badora, D.; Wawer, R.; Nierobca, A.; Krol-Badziak, A.; Kozyra, J.; Jurga, B.; Nowocien, E. Modelling the Hydrology of an Upland Catchment of Bystra River in 2050 Climate Using RCP 4.5 and RCP 8.5 Emission Scenario Forecasts. Agriculture 2022, 12, 403. [Google Scholar] [CrossRef]
  9. Arnold, J.G.; Kiniry, J.R.; Srinivasan, R.; Williams, J.R.; Haney, E.B.; Neitsch, S.L. Soil and Water Assessment Tool Theoretical Documentation, Version 2012; Texas Water Resources Institute: Forney, TX, USA, 2012. [Google Scholar]
  10. Neitsh, S.I.; Arnold, J.G.; Kiniry, J.R.; Williams, J.R. Soil and Water Assessment Tool Theoretical Documentation Version 2009; Texas Water Resources Institute: Forney, TX, USA, 2011. [Google Scholar]
  11. Kozyra, J.; Żyłowska, K.; Nieróbca, A.; Matyka, M.; Smagacz, J.; Jadczyszyn, T.; Wawer, R. Zmiany Klimatu a Rolnictwo w Polsce Ocena Zagrożeń i Sposoby Adaptacji, Fundacja Na Rzecz Zrównoważonego Rozwoju; Publisher: Warsaw, Poland, 2019; p. 59. [Google Scholar]
  12. Doroszewski, A.; Jadczyszyn, J.; Kozyra, J.; Pudełko, R.; Stuczyński, T.; Mizak, K.; Łopatka, A.; Koza, P.; Górski, T.; Wróblewska, E. Podstawy systemu monitoringu suszy rolniczej. Woda-Sr.-Obsz. Wiej. 2012, 12, 78–91. [Google Scholar]
  13. KLIMADA 2.0, 2019. KLIMADA 2.0–Baza Wiedzy o Zmianach Klimatu, Scenariusze Zmian Klimatu. 2022. Available online: https://klimada2.ios.gov.pl/ (accessed on 12 January 2020).
  14. Kundzewicz, Z. Zmiany klimatu, ich przyczyny i skutki–możliwości przeciwdziałania i adaptacji. Studia BAS 2012, 29, 9–30. [Google Scholar]
  15. Kundzewicz, Z.; Matczak, P. Climate change regional review: Poland. Wiley Interdiscip. Rev. Clim. Chang. 2012, 3, 297–311. [Google Scholar] [CrossRef]
  16. Kundzewicz, Z.W.; Piniewski, M.; Mezghani, A.; Okruszko, T.; Pińskwar, I.; Kardel, I.; Hov, Ø.; Szcześniak, M.; Szwed, M.; Benestad, R.E.; et al. Assessment of climate change and associated impact on selected sectors in Poland. Acta Geophys. 2018, 66, 1509–1523. [Google Scholar] [CrossRef] [Green Version]
  17. Nieróbca, A.; Kozyra, J.; Mizak, K.; Wróblewska, E. Zmiana długości okresu wegetacyjnego w Polsce. Woda-Sr.-Obsz. Wiej. 2013, 13, 81–94. [Google Scholar]
  18. Jacob, D.; Petersen, J.; Eggert, B.; Alias, A.; Christensen, O.B.; Bouwer, L.M.; Braun, A.; Colette, A.; Déqué, M.; Georgievski, G. EURO-CORDEX: New high-resolution climate change projections for European impact research. Reg. Environ. Chang. 2014, 14, 563–578. [Google Scholar] [CrossRef]
  19. Kovats, R.S.; Valentini, R.; Bouwer, L.M.; Georgopoulou, E.; Jacob, D.; Martin, E.; Rounsevell, M.; Soussana, J.-F. Europe. In Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part B: Regional Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Barros, V.R., Field, C.B., Dokken, D.J., Mastrandrea, M.D., Mach, K.J., Bilir, T.E., Chatterjee, M., Ebi, K.L., Estrada, Y.O., Genova, R.C., et al., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2014; pp. 1267–1326. [Google Scholar]
  20. IPCC. 2021: Summary for Policymakers. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Masson-Delmotte, V.P., Zhai, A., Pirani, S.L., Connors, C., Péan, S., Berger, N., Caud, Y., Chen, L., Goldfarb, M.I., Gomis, M., et al., Eds.; Cambridge University Press: Cambridge, UK, 2021; in press. [Google Scholar]
  21. Parry, M.L.; Canziani, O.F.; Palukitof, J.P.; van der Linden, P.J.; Hanson, C.E. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2007; Available online: https:/www.ipcc.ch/publications_and_data/ar4/wg2/en/contents.html (accessed on 5 December 2021).
  22. Załatwiaj Sprawy Urzędowe—Przez Internet, Bezpiecznie i Wygodnie! Available online: https://www.gov.pl/web/rolnictwo/kodeks-dobrych-praktyk-wodnych-w-rolnictwie (accessed on 15 April 2020).
  23. Ministerstwo Środowiska, 2013, Strategiczny Plan Adaptacji dla Sektorów i Obszarów Wrażliwych na Zmiany Klimatu do Roku 2020 z Perspektywą do Roku 2030, Dokument Został Opracowany Przez Ministerstwo Środowiska na Podstawie Analiz Wykonanych Przez Instytut Ochrony Środowiska—Państwowy Instytut Badawczy w Ramach Projektu pn. “Opracowanie i Wdrożenie Strategicznego Planu Adaptacji dla Sektorów i Obszarów Wrażliwych na Zmiany Klimatu—KLIMADA”, Realizowanego na Zlecenie MŚ w Latach 2011–2013 ze Środków Narodowego Funduszu Ochrony Środowiska i Gospodarki Wodnej. Warsaw, Poland, 2013. Available online: https://bip.mos.gov.pl/fileadmin/user_upload/bip/strategie_plany_programy/Strategiczny_plan_adaptacji_2020.pdf (accessed on 15 March 2021).
  24. Doroszewski, A.; Jóźwicki, T.; Wróblewska, E.; Kozyra, J. Susza Rolnicza w Polsce w Latach 1961–2010; Wyd. IUNG: Pulawy, Poland, 2014; p. 144. [Google Scholar]
  25. Doroszewski, A. Susza Rolnicza w Polsce w 2015 Roku. Warszawa, Poland. 2016. Available online: http://gwppl.org/data/uploads/prezentacje/4.%20Susza%20rolnicza_ADoroszewski.pdf (accessed on 18 November 2020).
  26. Chałubińska, A.; Wilgat, T. Podział Fizjograficzny woj. Lubelskiego, Przewodnik V Ogólnopolskiego Zjazdu Polskiego Towarzystwa Geograficznego; Oddział Lubelski PTG: Lublin, Poland, 1954; pp. 3–44. [Google Scholar]
  27. Jahn, A. Wyżyna Lubelska: Rzeźba i Czwartorzęd; Prace Geograficzne Instytutu Geograficznego PAN, Nr 7, IGiPAN, PWN: Warsaw, Poland, 1956. [Google Scholar]
  28. Sadurska, E. Charakterystyka Fizycznogeograficzna Dorzecza Bystrej; IUNG: Puławy, Poland, 1980; p. 29. [Google Scholar]
  29. Ziemnicki, S.; Pałys, S. Erozja wodna w zlewni rzeki Bystrej. Zesz. Probl. Postępów Nauk. Rol. 1977, 193, 44–71. [Google Scholar]
  30. Wawer, R.; Nowocień, E.; Podolski, B.; Capała, M. Ocena zagrożenia erozją wodną powierzchniową zlewni rzeki Bystrej z wykorzystaniem modelowania przestrzennego. Przegląd Nauk. SGGW Inżynieria I Kształtowanie Sr. Ann. XVII 2008, 3, 20–28. [Google Scholar]
  31. Jurga, B.; Wawer, R.; Kęsik, K. Zlewnia Rzeki Bystrej Jako Przykład Wyżynnej Zlewni Rolniczej o Wysokich Zdolnościach Buforowych Względem Fosforu-Studium Przypadku, Rolnictwo XXI Wieku–Problemy i Wyzwania, Pod Redakcją Dety Łuczyckiej; Idea Knowledge Future: Wrocław, Poland, 2018; pp. 143–154. ISBN 978-83-945311-9-5. [Google Scholar]
  32. Wawer, R.; Nowocień, E.; Kozyra, J. Hydrologia i denudacja w zlewni rzeki Bystrej. In Proceedings of the Konferencja Problemy Gospodarowania Zasobami Środowiska w Dolinach Rzecznych, Wrocław, Poland, 27–29 May 2015. [Google Scholar]
  33. SMGP. Szczegółowa Mapa Geologiczna Polski, Arkusz 747–Nałęczów (M-34-33-A). 2006. Available online: https://bazadata.pgi.gov.pl/data/smgp/arkusze_skany/smgp0747.jpg. (accessed on 5 January 2021).
  34. Maruszczak, H. Definicja i klasyfikacja lessów oraz utworów lessopodobnych. Przegląd Geol. 2000, 48, 580–586. [Google Scholar]
  35. Kalarus, K. Wpływ Materiału Macierzystego na Właściwości Gleb Wykształconych Na Lessie; Uniwersytet Jagieloński: Krakow, Poland, 2009. [Google Scholar]
  36. Piest, R.F.; Ziemnicki, S. Comparative Erosion Rates of Loeass Soils in Poland and Iowa. Trans. Asae 1979, 22, 822–827. [Google Scholar] [CrossRef]
  37. USDA. United States Department of Agriculture. 1996. Available online: https://www.usda.gov/ (accessed on 1 December 2020).
  38. Arnold, J.G.; Srinivasan, R.; Muttiah, R.; Williams, J. Large area hydrologic modelling and assessment. P. I: Model development. J. Am. Water Resour. Assoc. 1998, 34, 73–89. [Google Scholar]
  39. Miatkowski, Z.; Smarzyńska, K. Kalibracja i walidacja modelu SWAT do szacowania bilansu wodnego i strat azotu w małym dziale wodnym w centralnej Polsce. J. Water Land Dev. 2016, 29, 31–47. [Google Scholar] [CrossRef]
  40. QGIS. Quantum GIS 3.10.13 Coruna. 2020. Available online: https://www.qgis.org/pl/site/index.html (accessed on 3 March 2020).
  41. Winchell, M.; Srinivasan, R. SWAT Editor for SWAT2012—Documentation; Blackland Research Center: Temple, TX, USA, 2012; pp. 1–14. [Google Scholar]
  42. Bajkiewicz-Grabowska, E.; Mikulski, Z. Hydrologia Ogólna, Pod Redakcją Krystyny Wojtala, Wydawnictwo Naukowe PWN; IBUK: Warsaw, Poland, 2022; ISBN 978-83-01-14579-8. [Google Scholar]
  43. Neitsch, S.L.; Arnold, J.G.; Kiniry, J.R.; Williams, J.R. Soil and Water Assessment Tool. Theoretical Documentation, 2005. Available online: https://swatmodel.tamu.edu/media/1292/swat2005theory.pdf (accessed on 2 January 2022).
  44. Abbaspour, K.C.; Vejdani, M.; Haghighat, S. SWAT-CUP Calibration and Uncertainty Programs for SWAT. In Proceedings of the International Congress on Modelling and Simulation (MODSI’07), Christchurch, New Zealand, 10–13 December 2007; Oxley, R.L., Kulasiri, D., Eds.; Modelling and Simulation Society of Australia and New Zealands: Melbourne, Australia, 2007; pp. 1596–1602. Available online: https://www.mssanz.org.au/MODSIM07/papers/24_s17/SWAT-CUP_s17_Abbaspour_.pdf (accessed on 3 December 2021).
  45. Abbaspour, K.C. SWAT-CUP 2012: SWAT Calibration and Uncertainty Programs—A User Manual. Swiss Fed. Inst. Aquat. Sci. Technol. 2015, 2015, 1–98. [Google Scholar]
  46. Abbaspour, K.C.; Vaghefi, S.A.; Srinivasan, R.A. Guideline for Successful Calibration and Uncertainty Analysis for Soil and Water Assessment: A Review of Papers from the 2016 International SWAT Conference. Water 2018, 10, 6. [Google Scholar] [CrossRef] [Green Version]
  47. Bilondi, M.P.; Abbaspour, K.C.; Ghahraman, B. Application of three different calibration-uncertainty analysis methods in a semi-distributed rainfall-runoff model application. Middle-East J. Sci. Res. 2013, 15, 1255–1263. [Google Scholar]
  48. Yang, W.; Andréasson, J.; Phil Graham, L.; Olsson, J.; Rosberg, J.; Wetterhall, F. Distribution based scaling to improve usability of RCM regional climate model projections for hydrological climate change impacts studies. Hydrol. Res. 2010, 41, 211–229. [Google Scholar] [CrossRef]
  49. CODGiK, 2013, Centralny Ośrodek Dokumentacji Geodezyjnej i Kartograficznej. Available online: https://www.codgik.gov.pl/, (accessed on 2 February 2018).
  50. MPHP. 2017. Komputerowa Mapa Podziału Hydrograficznego Polski. Available online: https://danepubliczne.gov.pl/dataset?q=zlewnia&sort=metadata_modified+desc (accessed on 4 June 2018).
  51. CLC, 2018, CORIN–Land Cover-CLC 2018, Główny Inspektorat Ochrony Środowiska. Available online: https://clc.gios.gov.pl/index.php/clc-2018/o-clc2018 (accessed on 4 June 2018).
  52. Geoportal, 2020, Instytucja Odpowiedzialna: Główny Urząd Geodezji i Kartografii. Available online: https://mapy.geoportal.gov.pl/wss/service/PZGIK/ORTO/WMS/HighResolution (accessed on 5 March 2020).
  53. OSM, 2018, Open Street Map. Available online: https://download.geofabrikolandroland.html (accessed on 7 March 2020).
  54. IUNG-PIB, Digital soil-agriculture maps 1:25000 and 1:100000, Pulawy, Poland, 2010.
  55. Jadczyszyn, J.; Smreczak, B. Mapa glebowo-rolnicza w skali 1:25 000 i jej wykorzystanie na potrzeby współczesnego rolnictwa. Studia I Rap. IUNG-PIB 2017, 51, 9–27. [Google Scholar] [CrossRef]
  56. IMGW, 2019, Instytut Meteorologii i Gospodarki Wodnej PIB. Available online: https://danepubliczne.imgw.pl/data/dane_pomiarowo_obserwacyjne/ (accessed on 3 March 2019).
  57. KPOSK. Krajowy Program Oczyszczania Ścieków Komunalnych. 2017. Available online: https://www.kzgw.gov.pl/index.php/pl/materialy-informacyjne/programy/krajowy-program-oczyszczania-sciekow-komunalnych, (accessed on 9 March 2020).
  58. Markowski, K. Rolnictwo w Województwie Lubelskim w 2019 r; Urząd Statystyczny w Lublinie: Lublin, Poland, 2020; ISSN 2080-0517. [Google Scholar]
  59. Lasy Regionu, 2021, Regionalna Dyrekcja Lasów Państwowych w Lublinie. Available online: https://www.lublin.lasy.gov.pl/lasy-regionu#.yg89jegzzaq (accessed on 13 May 2021).
  60. Kouchi, D.M.; Esmaili, K.; Faridhosseini, A.; Sanaeinejad, S.H.; Khalili, D.; Abbaspour, K.C. Sensitivity of Calibrated Parameters and Water Resource Estimates on Different Objective Functions and Optimization Algorithms. Water 2017, 9, 384. [Google Scholar] [CrossRef] [Green Version]
  61. Abbaspour, K.C. SWAT-CUP Tutorial (2): Introduction to SWAT-CUP Program, Parameter Estimator (SPE), 2020. Available online: https://www.youtube.com/watch?v=nNsDPhOI7cc&ab_channel=2w2e (accessed on 7 July 2021).
  62. ADMS. System Monitorowania Suszy Rolniczej. 2013. Available online: https://susza.iung.pulawy.pl/system/ (accessed on 10 February 2021).
  63. Hennemuth, T.I.; Jacob, D.; Keup-Thiel, E.; Kotlarski, S.; Nikulin, G.; Otto, J.; Rechid, D.; Sieck, K.; Sobolowski, S.; Szabó, P. Guidance for EURO-CORDEX Climate Projections Data Use. Version 1.0-2017.08. Available online: https://www.euro-cordex.net/imperia/md/content/csc/cordex/euro-cordex-guidelines-version1.0-2017.08.pdf (accessed on 13 April 2021).
  64. Moss, R.H.; Edmonds, J.A.; Hibbard, K.A.; Manning, M.R.; Rose, S.K.; Van Vuuren, D.P.; Carter, T.R.; Emori, S.; Kainuma, M.; Kram, T. The next generation of scenarios for climate change research and assessment. Nature 2010, 463, 747–756. [Google Scholar] [CrossRef]
  65. Thomson, A.M.; Calvin, K.V.; Smith, S.J.; Kyle, G.P.; Volke, A.; Patel, P.; Delgado-Arias, S.; Bond-Lamberty, B.; Wise, M.A.; Clarke, L.E. RCP4. 5: A pathway for stabilization of radiative forcing by 2100. Clim. Zmiana 2011, 109, 77–94. [Google Scholar] [CrossRef] [Green Version]
  66. Landelius, T.; Dahlgren, P.; Gollvik, S.; Jansson, A.; Olsson, E. A high-resolution regional reanalysis for Europe. Part 2: 2D analysis of surface temperature, precipitation and wind. Q. J. R. Meteorol. Soc. 2016, 142, 2132–2142. [Google Scholar] [CrossRef] [Green Version]
  67. Schulzweida, U.; Kornblueh, L.; Quast, R. CDO user guide. Climate Data Operators, Version 2006, 1, 205–209. [Google Scholar]
  68. PIK, 2012, Potsdam Institute for Climate Impact Research. Available online: https://www.pik-potsdam.de/~mmalte/rcps/ (accessed on 5 August 2021).
  69. Meinshausen, M.; Smith, S.J.; Calvin, K.; Daniel, J.S.; Kainuma, M.L.; Lamarque, J.F.; Matsumoto, K.; Montzka, S.A.; Raper, S.C.B.; Riahi, K.; et al. The RCP Greenhouse Gas Concentrations and their Extension from 1765 to 2300. Clim. Chang. 2011, 109, 213–241. [Google Scholar] [CrossRef] [Green Version]
  70. Hawes, E.; Smith, M. Riparian Buffer Zones: Functions and Recommended Widths Eightmile River Wild and Scenic Study Committee; 2005. Available online: http://www.eightmileriver.org/resources/digital_library/appendicies/09c3_Riparian%20Buffer%20Science_YALE.pdf (accessed on 5 March 2022).
  71. Mayer, P.M.; Reynolds, S.K.; Canfield, T.J., Jr. Riparian Buffer Width, Vegetative Cover, and Nitrogen Removal Effectiveness: A Review of Current Science and Regulations; U.S. Environmental Protection Agency Office of Research and Development National Risk Management Research Laboratory: Gothenburg, Sweden, 2005. [Google Scholar]
  72. Pawlaczyk, P.; Biedroń, I.; Brzóska, P.; Dondajewska-Pielka, R.; Furdyna, A.; Gołdyn, R.; Grygoruk, M.; Grześkowiak, A.; Horska-Schwarz, S.; Jusik, S.; et al. Podręcznik Dobrych Praktyk Renaturyzacji Wód Powierzchniowych. OPRAC. w Ramach Przedsięwzięcia „Opracowanie Krajowego Programu Renaturyzacji Wód Powierzchniowych”; Państwowe Gospodarstwo Wodne Wody Polskie, Krajowy Zarząd Gospodarki Wodnej: Warszawa, Poland, 2020. [Google Scholar]
  73. Waidler, D.; White, M.J.; Steglich, E.; Wang, S.; Williams, J.R.; Jones, C.A.; Srinivasan, R. Conservation Practice Modeling Guide for SWAT and APEX; Texas Water Resources Institute: Forney, TX, USA, 2011. [Google Scholar]
  74. Arabi, M.; Frankenberger, J.R.; Engel, B.A.; Arnold, J.G. Representation of Agricultural Conservation Practices with SWAT. Hydrol. Processes 2008, 22, 3042–3055. [Google Scholar] [CrossRef]
  75. Busari, M.A.; Kukal, S.S.; Kaur, A.; Bhatt, R.; Dulazi, A.A. Conservation tillage impacts on soil, crop and the environment. Int. Soil Water Conserv. Res. 2015, 3, 119–129. [Google Scholar] [CrossRef] [Green Version]
  76. Somasundaram, J.; Sinha, N.K.; Dalal, R.C.; Rattan, L.; Mohanty, M.; Naorem, A.K.; Hati, K.M.; Chaudhary, R.S.; Biswas, A.K.; Patra, A.K.; et al. No-Till Farming and Conservation Agriculture in South Asia–Issues, Challenges, Prospects and Benefits. Crit. Rev. Plant Sci. 2020, 39, 236–279. [Google Scholar] [CrossRef]
  77. Kuś, J. Glebowa materia organiczna–znaczenie, zawartość i bilansowanie. Studia I Rap. IUNG–PIB 2015, 45, 27–53. [Google Scholar]
  78. Mackey, B.; Prentice, I.; Steffen, W.; House, J.; Lindenmayer, D.; Keith, H.; Berry, S. Untangling the confusion around land carbon science and climate change mitigation policy. Nat. Clim. Change 2013, 3, 552–557. [Google Scholar] [CrossRef]
  79. Będkowska, H. Lasy i Zmiany Klimatu, Centrum Informacyjne Lasów Państwowych; Mozolewska-Adamczyk, M., Ed.; Dyrekcja Generalna Lasów Państwowych, Ośrodek Rozwojowo-Wdrożeniowy Lasów Państwowych w Bedoniu: Bedoń, Poland, 2016; ISBN 978-83-63895-62-4. [Google Scholar]
  80. Józefaciuk, A.; Nowocień, E.; Wawer, R. Erozja Gleb w Polsce–Skutki Środowiskowe i Gospodarcze, Działania Zaradcze; Monografie i Rozprawy Naukowe IUNG-PIB: Pulawy, Poland, 2014; Volume 44, p. 263. [Google Scholar]
  81. Crowther, T.W.; Glick, H.B.; Covey, K.R.; Bettigole, C.; Maynard, D.S.; Thomas, S.M.; Smith, J.R.; Hintler, G.; Duguid, M.C.; Amatulli, G.; et al. Mapping tree density at a global scale. Nature 2015, 525, 201–205. [Google Scholar] [CrossRef]
  82. Borecki, T.; Malinowski, S.; Banasik, K.; Okruszko, T.; Brzeziecki, B.; Kozyra, J.; Gajda, N.; Haman, K. LXXIII Zmiany Klimatu i Ich Następstwa, Instytut Problemów Współczesnej Cywilizacji im. Marka Dietricha; Wydawnictwo SGGW: Warszawa, Poland, 2021; p. 68. ISBN 978-83-89871-44-2. [Google Scholar]
  83. MOSZNiL, 1997, Polityka Leśna Państwa. Dokument Przyjęty Przez Radę Ministrów w Dniu 22 Kwietnia 1997 r. Ministerstwo Ochrony Środowiska, Zasobów Naturalnych i Leśnictwa, Warszawa, Krajowy Program Zwiększania Lesistości-Stan i Trudności Realizacji z Perspektywy Lokalnej/National Program for Expanding of Forest Cover-Implementation and Its Difficulties from a Local View. Available online: https://www.researchgate.net/publication/322358247_Krajowy_program_zwiekszania_lesistosci_-_stan_i_trudnosci_realizacji_z_perspektywy_lokalnej_National_Program_for_Expanding_of_Forest_Cover_-_implementation_and_its_difficulties_from_a_local_view (accessed on 11 January 2022).
  84. Report on the Condition of Forests in Poland 2021, Wydano na Zlecenie Dyrekcji Generalnej Lasów Państwowych, p.o. Dyrektora Generalnego Lasów Państwowych mgr inż. Józef Kubica. Warsaw, Poland, 2021. Available online: https://www.gov.pl/attachment/23dcf2c9-6514-45cf-91de-33cebbf06c49 (accessed on 15 March 2022).
  85. Kaliszewski, A. Krajowy program zwiększania lesistości-stan i trudności realizacji z perspektywy lokalnej/National Program for Expanding of Forest Cover-implementation and its difficulties from a local view. Studia I Mater. CEPL W Rogowie. 2016, 49, 7–19. [Google Scholar]
  86. Trabucco, A.; Zomer, R.; Bossio, D.A.; van Straaten, O.; Verchot, L. Climate change mitigation through afforestation/reforestation: A global analysis of hydrologic impacts with four case studies. Agric. Ecosyst. Environ. 2008, 126, 81–97. [Google Scholar] [CrossRef]
  87. Kujawa, A.; Kujawa, K.; Zajączkowski, J.; Borek, R.; Tyszko-Chmielowiec, P.; Chmielowiec-Tyszko, D.; Józefczuk, J.; Krukowska-Szopa, I.; Śliwa, P.; Witkos-Gnach, K. Zadrzewienia Na Obszarach Wiejskich–Dobre Praktyki i Rekomendacje; Fundacja Ekorozwoju: Wrocław, Poland, 2019. [Google Scholar]
  88. Woch, F. Urządzeniowe metody zmniejszania zagrożenia erozyjnego gleb. Studia I Rap. IUNG-PIB Zesz. 2008, 10, 79–102. [Google Scholar] [CrossRef]
  89. Tałałaj, Z. Wpływ zadrzewień na plonowanie roślin rolniczych. W: Znaczenie zadrzewień w krajobrazie rolniczym oraz aktualne problemy ich rozwoju w przyrodniczo-gospodarczych warunkach Polski. Mat. Konf. Płock 1997, 72, 91. [Google Scholar]
  90. Podolski, B.; Woch, F. Wpływ bezpośredniego oddziaływania lasów i zadrzewień śródpolnych na plonowanie zbóż ozimych. Pam. Puł. 1999, 119, 101–111. [Google Scholar]
  91. Lambor, J. Podstawy i Zasady Gospodarki Wodnej; Wydawnictwa Komunikacji i Łączności, Instrukcje i Podręczniki—Publisher of Communication and Communications: Warsaw, Poland, 1965. [Google Scholar]
  92. Bogusz, A.; Tokarczyk, T. Rola Terenów Zalesionych w Kształtowaniu Retencji Wód Opadowych w Zlewniach Zurbanizowanych; Monografie Komitetu Gospodarki Wodnej PAN: Warsaw, Poland, 2016. [Google Scholar]
  93. Ozga-Zielińska, M.; Brzeziński, J. Hydrologia Stosowana; PWN: Warszawa, Poland, 1997; 323p. [Google Scholar]
  94. The Woodland Trust, 2016, Keeping Rivers Cool: A Guidance Manual. Creating Riparian Shade for Climate Change Adaptation. Available online: https://www.woodlandtrust.org.uk/publications/2016/02/keeping-rivers-cool/ (accessed on 27 December 2021).
  95. Jabłońska, E.; Wiśniewska, M.; Marcinkowski, P.; Grygoruk, M.; Walton, C.R.; Zak, D.; Hoffmann, C.C.; Larsen, S.E.; Trepel, M.; Kotowski, W. Catchment-Scale Analysis Reveals High CostEffectiveness of Wetland Buffer Zones as a Remedy to Non-Point Nutrient Pollution in NorthEastern Poland. Water 2020, 21, 629. [Google Scholar] [CrossRef] [Green Version]
  96. Babel, M.S.; Gunathilake, M.B.; Jha, M.K. Evaluation of Ecosystem-Based Adaptation Measures for Sediment Yield in a Tropical Watershed in Thailand. Water 2021, 13, 2767. [Google Scholar] [CrossRef]
  97. Wawer, R.; Kozyra, J. Kruchy bilans wody na polach. Top-Agrar. 2021, 5, 104–107.
  98. Claire, J.; Berglund, M.; Bluz, K.; Dworak, T.; Marras, S.; Mereu, V.; Michetti, M. Climate Change Adaptation in the Agriculture Sector in Europe; Publications Office of the European Union: Luxembourg, 2019. [Google Scholar]
  99. Derpsch, R.; Friedrich, T.; Kassam, A.; Hongwen, L. Current Status of Adoption of No-Till Farming in the World and some of its Main Benefits. Int. J. Agric. Biol. Eng. 2010, 3, 1–25. [Google Scholar] [CrossRef]
  100. Fengyun, Z.; Pute, W.; Xining, Z.; Xuefeng, C. The effects of no-tillage practice on soil physical properties. Afr. J. Biotechnol. 2011, 10, 17645–17650. [Google Scholar] [CrossRef]
  101. Khan, N.U.; Khan, A.A.; Goheer, M.A.; Shafique, I.; Hussain, S.; Hussain, S.; Javed, T.; Naz, M.; Shabbir, R.; Raza, A.; et al. Effect of Zero and Minimum Tillage on Cotton Productivity and Soil Characteristics under Different Nitrogen Application Rates. Sustainability 2021, 13, 13753. [Google Scholar] [CrossRef]
  102. Kundzewicz, Z.; Kozyra, J. Ograniczenie wpływu zagrożeń klimatycznych w odniesieniu do rolnictwa i obszarów wiejskich. Pol. J. Agron. 2011, 7, 68–81. [Google Scholar]
  103. Zhao, F.; Wu, Y.; Hui, J.; Sivakumar, B.; Meng, X.; Liu, S. Projected soil organic carbon loss in response to climate warming and soil water content in a loess watershed. Carbon Balance Manag. 2021, 16, 24. [Google Scholar] [CrossRef]
  104. Humphrey, V.; Zscheischler, J.; Ciais, P.; Gudmundsson, L.; Sitch, S.; Seneviratne, S.I. Sensitivity of atmospheric CO2 growth rate to observed changes in terrestrial water storage. Nature 2018, 560, 628–631. [Google Scholar] [CrossRef]
  105. Dembek, W.; Wiatkowski, M.; Żurek, G.; Kuś, J. Innowacyjne Metody Gospodarowania Zasobami Wody w Rolnictwie, Praca Zbiorowa Pod Red; Dembek, W., Kuś, J., Wiatkowski, M., Zurek, G., Eds.; Centrum Doradztwa Rolniczego w Brwinowie: Brwinów, Poland, 2016. [Google Scholar]
  106. Woch, F. Kompleksowe Scalanie Gruntów Rolnych i Leśnych Oraz Jego Wpływ na Środowisko, Materiały Szkoleniowe nr 93/2006; Bochniarz, A., Ed.; IUNG-PIB w Puławach, Dział Upowszechniania i Wydawnictw IUNG-PIB: Puławy, Poland, 2006; ISBN 83-89576-43-0. [Google Scholar]
  107. Woch, F.; Pijanowski, J.; Kuryłowicz, T. Kompleksowe urządzanie obszarów wiejskich jako szansa dla rozwoju wsi. Pol. J. Agron. 2018, 33, 16–32. [Google Scholar]
Figure 1. Location of the study area, Bystra catchment, with marked main tributaries and their catchments (own study).
Figure 1. Location of the study area, Bystra catchment, with marked main tributaries and their catchments (own study).
Water 14 02288 g001
Figure 2. Seasonal average soil water content (1.5 m) for 2041–2050 and for the SWAT 2010–2017 model for individual climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3 (own study).
Figure 2. Seasonal average soil water content (1.5 m) for 2041–2050 and for the SWAT 2010–2017 model for individual climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3 (own study).
Water 14 02288 g002
Figure 3. Comparison of average soil water content in 31 sub-catchments during the SWAT simulation period 2010–2017 and 2041–2050 for individual climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3 (own study).
Figure 3. Comparison of average soil water content in 31 sub-catchments during the SWAT simulation period 2010–2017 and 2041–2050 for individual climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3 (own study).
Water 14 02288 g003
Table 2. Meteorological data for the Bystra catchment [8].
Table 2. Meteorological data for the Bystra catchment [8].
Weather StationMeasurement Period
Precipitation (mm)Temperature (°C)Wind Speed (m/s)HumiditySolar Total Radiation (MJ/m2)
Pulawy2005–20172005–20172005–20172005–20172005–2017
Rogalow2005–2017
Lublin Radawiec2005–20172005–20172005–20172005–2017
Table 3. Description of GCM/RCM simulation with its division depending on radiative forcing. Comparison of temperature and precipitation changes in 2021–2050 in GCM/RCM simulation for RCP 4.5 and RCP 8.5 to the base period 1971–2000 (own study).
Table 3. Description of GCM/RCM simulation with its division depending on radiative forcing. Comparison of temperature and precipitation changes in 2021–2050 in GCM/RCM simulation for RCP 4.5 and RCP 8.5 to the base period 1971–2000 (own study).
ModelsScenario AssumptionsRadiative Forcing
GCM/RCM SimulationChange in
Average Annual Air Temperature
Change in Average Annual Precipitation+4.5 W m−2+8.5 W m−2
RCP 4.5RCP 8.5RCP 4.5RCP 8.5RCP 4.5RCP 8.5
EC-EARTH/RACMO22E+1.5 °C +1.8 °C+15%+6%RCP 4.5.1RCP 8.5.1
EC-EARTH/HIRHAM5+1.6 °C+1.9 °C+12%+5%RCP 4.5.2RCP 8.5.2
EC-EARTH/RCA4+1.6 °C+2.2 °C+15%+11%RCP 4.5.3RCP 8.5.3
Table 4. Division of the land cover and land use as well as the percentage of land use in the Bystra catchment generated in the QSWAT interface. CLC code 112–142 means artificial surfaces; code 211–243 means agricultural areas; code 313–324 means forest and semi natural areas; code 411 means wetlands, and code 511 is water bodies (own study).
Table 4. Division of the land cover and land use as well as the percentage of land use in the Bystra catchment generated in the QSWAT interface. CLC code 112–142 means artificial surfaces; code 211–243 means agricultural areas; code 313–324 means forest and semi natural areas; code 411 means wetlands, and code 511 is water bodies (own study).
Corine Land Cover LegendCLCSWATS-0AS-1AS-2
CodeCodePart (%)Part (%)Part (%)
Discontinuous urban fabric112URML0.920.90.9
Industrial or commercial units121UCOM1.551.491.49
Mineral extraction sites131UIDU0.020.020.02
Sport and leisure facilities142FESC0.020.020.02
SUM=2.512.432.43
Non-irrigated arable land211CRDY52.3550.5752.23
Vineyards221GRAP0.030.030.03
Fruit trees and berry plantations222ORCD10.8510.5510.83
Pastures231PAST5.895.355.55
Complex cultivation patterns242AGRL9.048.688.86
Land principally occupied by agriculture with significant areas of natural vegetation243CRGR0.050.050.05
SUM=78.2175.2377.55
Mixed forest313FRST16.3419.6517.37
Transitional woodland-shrub324SHRB2.432.182.23
Inland marshes411WEHB0.260.250.21
Water courses511WATR0.270.260.21
Table 5. Comparison of average soil water content by season for the SWAT 2010–2017 simulation period with climate projections (RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3) for the years 2041–2050 in the Bystra catchment. Bold numbers indicate soil water content, while shaded numbers indicate percentage change (red is % decrease in content; blue is % increase in content). Dark red and dark blue shading means large changes, while light red and light blue shading means small changes (own study).
Table 5. Comparison of average soil water content by season for the SWAT 2010–2017 simulation period with climate projections (RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3) for the years 2041–2050 in the Bystra catchment. Bold numbers indicate soil water content, while shaded numbers indicate percentage change (red is % decrease in content; blue is % increase in content). Dark red and dark blue shading means large changes, while light red and light blue shading means small changes (own study).
Climate Scenario RCP 4.5RCP 8.5
Climate ProjectionModel
2010–2017
RACMO22E
(RCP 4.5.1)
HIRHAM5
(RCP 4.5.2)
RCA4
(RCP 4.5.3)
RACMO22E
(RCP 8.5.1)
HIRHAM5
(RCP 8.5.2)
RCA4
(RCP 8.5.3)
Time interval 2041–2050
SeasonSeasonal average of soil water content (mm)
DJF344332336337340342340
−3.5%−2.3%−2.1%−1.3%−0.7%−1.2%
MAM322303311318318321321
−5.8%−3.3%−1.3%−1.2%−0.1%−0.2%
JJA309292291303313306306
−5.4%−5.6%−2.0%+1.4%−0.8%−1.0%
SON328313321322327329328
−4.4%−2.0%−1.8%−0.2%+0.4%+0.2%
Average annual326310315320324325324
−4.7%−3.3%−1.8%−0.4%−0.3%−0.6%
Table 6. Comparison of average soil water content by season between scenario 0 (S-0) and adaptation scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) for 2041–2050 in the Bystra catchment for climate projection RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3. Bold numbers indicate soil water content, and shaded numbers indicate percentage change (red indicates % decrease in content and blue indicates % increase in content). Dark red and dark blue shading indicates large changes, while light red and light blue shading indicates small changes (own study).
Table 6. Comparison of average soil water content by season between scenario 0 (S-0) and adaptation scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) for 2041–2050 in the Bystra catchment for climate projection RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3. Bold numbers indicate soil water content, and shaded numbers indicate percentage change (red indicates % decrease in content and blue indicates % increase in content). Dark red and dark blue shading indicates large changes, while light red and light blue shading indicates small changes (own study).
Time Interval2041–2050
Type of ScenarioS-0AS-1AS-2AS-3AS-4AS-5 S-0AS-1AS-2AS-3AS-4AS-5
SeasonSeasonal average of soil water content (mm)
DJF332318332332333330RACMO22E
(RCP 4.5.1)
340326340340339339RACMO22E
(RCP 8.5.1)
−4.2%−0.1%0.0%+0.1%−0.6%−4.0%0.0%0.0%−0.1%−0.1%
MAM303290303303303301318305318318318317
−4.1%0.0%0.0%+0.1%−0.7%−4.1%0.0%0.0%0.0%−0.2%
JJA292278292292293288313299313313314312
−4.8%−0.1%0.0%+0.3%−1.4%−4.6%−0.1%0.0%+0.1%−0.4%
SON313299313313315310327313327327328326
−4.6%−0.1%0.0%+0.4%−1.0%−4.4%−0.1%0.0%+0.2%−0.3%
Average annual310296310310311307 324311324324325324
−4.4%0.0%0.0%+0.2%−0.9% −4.3%0.0%0.0%+0.1%−0.2%
DJF336322336336337335HIRHAM5
(RCP 4.5.2)
342328342342342342HIRHAM5
(RCP 8.5.2)
−4.1%0.0%0.0%+0.1%−0.2%−4.0%0.0%0.0%0.0%0.0%
MAM311298311311311310321308321321321321
−4.1%0.0%0.0%0.0%−0.2%−4.0%0.0%0.0%0.0%−0.1%
JJA291277291291292288306292306306307304
−4.8%−0.1%0.0%+0.2%−1.1%−4.6%−0.1%0.0%+0.1%−0.6%
SON321307321321322319329315329329330328
−4.4%−0.1%0.0%+0.3%−0.5%−4.3%−0.1%0.0%+0.2%−0.2%
Average annual315301315315315313 325311324325325324
−4.3%−0.1%0.0%+0.2%−0.5% −4.2%0.0%0.0%+0.1%−0.3%
DJF337323337337337336RCA4
(RCP 4.5.3)
340326340340340340RCA4
(RCP 8.5.3)
−4.1%0.0%0.0%0.0%−0.2%−4.0%0.0%0.0%0.0%0.0%
MAM318305317318318317321308321321321321
−4.1%0.0%0.0%0.0%−0.2%−4.0%0.0%0.0%0.0%−0.1%
JJA303289302303303300306291305306306303
−4.6%−0.1%0.0%+0.2%−0.8%−4.7%−0.1%0.0%+0.2%−0.7%
SON322307322322323320328314328328329327
−4.5%−0.1%0.0%+0.2%−0.5%−4.3%−0.1%0.0%+0.2%−0.4%
Average annual320306320320320318 324310324324324323
−4.3%−0.1%0.0%+0.1%−0.4% −4.3%−0.1%0.0%+0.1%−0.3%
Table 7. Comparison of seasonal total runoff between scenario 0 (S-0) and adaptation scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) for 2041–2050 in the Bystra catchment for climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3. Bold numbers indicate soil water content, and shaded numbers indicate percentage changes (red indicates % decrease in content, and blue indicates % increase in content). Dark red and dark blue shading indicates large changes, while light red and light blue shading indicates small changes (own study).
Table 7. Comparison of seasonal total runoff between scenario 0 (S-0) and adaptation scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) for 2041–2050 in the Bystra catchment for climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3. Bold numbers indicate soil water content, and shaded numbers indicate percentage changes (red indicates % decrease in content, and blue indicates % increase in content). Dark red and dark blue shading indicates large changes, while light red and light blue shading indicates small changes (own study).
Time Interval2041–2050
Type of ScenarioS-0AS-1AS-2AS-3AS-4AS-5 S-0AS-1AS-2AS-3AS-4AS-5
SeasonSeasonal sum of total runoff (mm)
DJF353535353634RACMO22E
(RCP 4.5.1)
555455555654RACMO22E
(RCP 8.5.1)
−1.0%−0.1%0.0%+3.7%−3.1%−0.8%−0.1%0.0%+2.4%−1.3%
MAM313131313230474747474846
−0.5%−0.2%0.0%+3.6%−3.3%0.1%−0.1%0.0%+2.4%−1.1%
JJA303030303129494949495048
−0.2%−0.2%0.0%+3.2%−3.9%0.0%−0.1%0.0%+2.1%−1.8%
SON323232323431494849495048
−0.4%−0.1%0.0%+4.7%−3.5%−0.3%−0.1%0.0%+3.2%−1.6%
Annual sum128127128128133123 199199199199204196
−0.5%−0.2%0.0%+3.8%−3.4% −0.3%−0.1%0.0%+2.5%−1.5%
DJF393939394038HIRHAM5
(RCP 4.5.2)
525252525351HIRHAM5
(RCP 8.5.2)
−0.4%−0.2%0.0%+2.8%−2.7%−0.3%−0.1%0.0%+2.0%−1.6%
MAM434343434442616161616160
−0.2%−0.2%0.0%+2.4%−2.1%−0.2%−0.1%0.0%+1.3%−1.2%
JJA363636363735535353535452
+0.3%0.0%0.0%+2.3%−3.1%+0.2%0.0%0.0%+1.4%−1.6%
SON363636363735515151515250
−0.4%−0.1%0.0%+3.5%−3.3%−0.2%−0.1%0.0%+2.3%−1.7%
Annual sum154154154154159150 217216216217220213
−0.2%−0.1%0.0%+2.7%−2.8% −0.1%−0.1%0.0%+1.7%−1.5%
DJF505050505250RCA4
(RCP 4.5.3)
717171717270RCA4
(RCP 8.5.3)
−0.4%−0.1%0.0%+2.5%−1.7%−0.2%−0.1%0.0%+2.1%−1.3%
MAM515051515250686868686967
−0.2%−0.1%0.0%+2.2%−1.3%−0.2%−0.1%0.0%+1.8%−1.0%
JJA393939394039606060606159
+0.2%0.0%0.0%+2.4%−2.0%0.0%0.0%0.0%+1.6%−1.3%
SON444444444644706970707269
−0.3%−0.2%0.0%+3.4%−1.9%−0.5%−0.1%0.0%+2.5%−1.1%
Annual sum185184185185190182 268268268268274265
−0.2%−0.1%0.0%+2.6%−1.7% −0.2%−0.1%0.0%+2.0%−1.2%
Table 8. Comparison of seasonal sediment yield between scenario 0 (S-0) and adaptation scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) for 2041–2050 in the Bystra catchment for climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3. Bold numbers indicate soil water content, and shaded numbers indicate percentage changes (red indicates % decrease in content, and blue indicates % increase in content). Dark red and dark blue shading indicates large changes, while light red and light blue shading indicates small changes (own study).
Table 8. Comparison of seasonal sediment yield between scenario 0 (S-0) and adaptation scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) for 2041–2050 in the Bystra catchment for climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3. Bold numbers indicate soil water content, and shaded numbers indicate percentage changes (red indicates % decrease in content, and blue indicates % increase in content). Dark red and dark blue shading indicates large changes, while light red and light blue shading indicates small changes (own study).
Time Interval2041–2050
Type of ScenarioS-0AS-1AS-2AS-3AS-4AS-5 S-0AS-1AS-2AS-3AS-4AS-5
SeasonSeasonal sum of sediment yield (t/ha)
DJF0.170.160.170.050.120.16RACMO22E
(RCP 4.5.1)
0.250.230.250.070.180.24RACMO22E
(RCP 8.5.1)
−9%−1%−72%−28%−5%−9%0%−71%−30%−5%
MAM0.080.070.080.020.070.070.070.060.070.020.060.06
−9%0%−71%−14%−8%−7%0%−71%−7%−7%
JJA0.150.140.150.050.130.120.220.200.220.060.160.18
−9%−1%−70%−13%−21%−9%0%−72%−28%−17%
SON0.150.140.150.040.080.130.180.160.180.050.110.17
−7%−1%−72%−48%−14%−11%−1%−72%−41%−5%
Annual sum0.550.510.550.160.400.49 0.720.650.710.200.500.65
−9%−1%−71%−27%−12% −9%0%−72%−30%−9%
DJF0.110.100.100.030.080.10HIRHAM5
(RCP 4.5.2)
0.120.110.120.040.090.12HIRHAM5
(RCP 8.5.2)
−10%−1%−72%−23%−4%−9%0%−71%−22%−4%
MAM0.130.120.130.040.120.120.240.220.240.070.210.22
−6%0%−71%−7%−8%−8%−1%−71%−13%−8%
JJA0.120.110.120.030.090.100.210.200.210.060.160.17
−11%−2%−72%−25%−16%−8%−1%−71%−23%−18%
SON0.190.170.190.050.100.160.220.200.220.060.120.20
−9%−1%−73%−49%−13%−10%−1%−72%−44%−10%
Annual sum0.540.500.540.150.390.48 0.790.720.790.230.590.71
−9%−1%−72%−29%−11% −9%−1%−71%−25%−11%
DJF0.110.100.110.030.090.11RCA4
(RCP 4.5.3)
0.140.120.140.040.100.13RCA4
(RCP 8.5.3)
−9%−1%−73%−23%0%−9%−1%−72%−24%−7%
MAM0.150.140.150.050.140.130.160.150.160.050.150.14
−11%0%−70%−7%−12%−9%−1%−72%−6%−14%
JJA0.110.100.110.030.080.090.140.120.140.040.130.11
−12%−1%−73%−25%−21%−11%−1%−71%−7%−23%
SON0.200.180.200.060.120.180.570.510.560.160.320.50
−9%0%−72%−42%−10%−10%−1%−72%−43%−12%
Annual sum0.570.520.570.160.430.51 1.000.900.990.280.710.87
−10%−1%−72%−26%−11% −10%−1%−72%−29%−13%
Table 9. Comparison of seasonal actual evapotranspiration between scenario 0 (S-0) and adaptation scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) for 2041–2050 in the Bystra catchment for climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3. Bold numbers indicate soil water content and shaded numbers indicate percentage changes (red indicates % decrease in content, and blue indicates % increase in content). Dark red and dark blue shading indicates large changes, while light red and light blue shading indicates small changes (own study).
Table 9. Comparison of seasonal actual evapotranspiration between scenario 0 (S-0) and adaptation scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) for 2041–2050 in the Bystra catchment for climate projections RCP 4.5.1, RCP 4.5.2, RCP 4.5.3, RCP 8.5.1, RCP 8.5.2, RCP 8.5.3. Bold numbers indicate soil water content and shaded numbers indicate percentage changes (red indicates % decrease in content, and blue indicates % increase in content). Dark red and dark blue shading indicates large changes, while light red and light blue shading indicates small changes (own study).
Time Interval2041–2050
Type of ScenarioS-0AS-1AS-2AS-3AS-4AS-5 S-0AS-1AS-2AS-3AS-4AS-5
SeasonSeasonal sum of actual evapotranspiration (mm)
DJF272727272727RACMO22E
(RCP 4.5.1)
292929292929RACMO22E
(RCP 8.5.1)
−0.9%−0.1%0.0%−0.2%−0.4%−0.7%0.0%0.0%−0.4%−0.5%
MAM154151154154154156156154156156156157
−1.9%−0.1%0.0%0.0%+1.4%−1.7%−0.1%0.0%−0.3%+0.7%
JJA166169166166163168165168165165162166
+2.0%+0.2%0.0%−1.5%+1.4%+1.8%+0.2%0.0%−1.5%+1.0%
SON707070706770707070706870
+0.4%0.0%0.0%−4.0%+0.1%+0.3%0.0%0.0%−3.2%+0.1%
Annual sum416417416416411421 420420420420415423
+0.1%0.0%0.0%−1.3%+1.1% +0.1%0.0%0.0%−1.3%+0.7%
DJF242424242424HIRHAM5
(RCP 4.5.2)
232323232323HIRHAM5
(RCP 8.5.2)
−0.7%0.0%0.0%−0.2%−0.4%−0.7%0.0%0.0%−0.3%−0.4%
MAM149146148149148150135133135135136137
−1.9%−0.1%0.0%−0.1%+1.2%−1.7%−0.1%0.0%+0.1%+1.0%
JJA152155152152150154152155153152150154
+2.0%+0.2%0.0%−1.2%+1.7%+1.7%+0.2%0.0%−1.3%+1.3%
SON616161615961656665656365
+0.3%0.0%0.0%−3.6%0.0%+0.1%0.0%0.0%−3.3%0.0%
Annual sum386386386386382390 377377377377372380
+0.1%0.0%0.0%−1.1%+1.1% 0.0%0.0%0.0%−1.1%+0.9%
DJF313131313131RCA4
(RCP 4.5.3)
343434343434RCA4
(RCP 8.5.3)
−0.7%0.0%0.0%−0.2%−0.6%−0.8%0.0%0.0%−0.2%−0.5%
MAM136134136136136137141138141141141142
−1.7%−0.1%0.0%−0.2%+0.8%−1.6%−0.1%0.0%−0.1%+0.8%
JJA168170168168165170158161158158156160
+1.5%+0.2%0.0%−1.4%+1.3%+1.7%+0.2%0.0%−1.5%+1.3%
SON697069696769727373726972
+0.4%0.0%0.0%−3.5%0.0%+0.2%0.0%0.0%−4.2%0.0%
Annual sum404405405404399408 406406406406400409
+0.1%0.0%0.0%−1.3%+0.8% +0.1%0.0%0.0%−1.4%+0.7%
Table 10. Percent summary of changes in soil water content, sediment yield, total runoff and actual evapotranspiration under adaptation scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) compared to scenario 0 (S-0) (created from Supplementary Material: Figures S1 and S2), for averages of three GCMs/RCMs combinations under two RCP climate change scenarios (RCP 4.5, RCP 8.5). The summary is for four seasons (DJF, MAM, JJA, SON) in the Bystra catchment. Shaded numbers indicate percentage changes (red indicates % decrease in content, and blue indicates % increase in content). Dark red and dark blue shading indicates large changes, while light red and light blue shading indicates small changes (own study).
Table 10. Percent summary of changes in soil water content, sediment yield, total runoff and actual evapotranspiration under adaptation scenarios 1–5 (AS-1, AS-2, AS-3, AS-4, AS-5) compared to scenario 0 (S-0) (created from Supplementary Material: Figures S1 and S2), for averages of three GCMs/RCMs combinations under two RCP climate change scenarios (RCP 4.5, RCP 8.5). The summary is for four seasons (DJF, MAM, JJA, SON) in the Bystra catchment. Shaded numbers indicate percentage changes (red indicates % decrease in content, and blue indicates % increase in content). Dark red and dark blue shading indicates large changes, while light red and light blue shading indicates small changes (own study).
RCP 4.5RCP 8.5
SeasonSoil water Content (mm)Total Runoff (mm)Sediment Yield (t/ha)Actual Evapotranspiration (mm)Soil Water Content (mm)Total Runoff (mm)Sediment Yield (t/ha)Actual Evapotranspiration (mm)
DJF−4.1%−0.6%−9.3%−0.7%−4.0%−0.4%−9.1%−0.7%AS-1
MAM−4.1%−0.2%−8.7%−1.9%−4.0%−0.1%−8.5%−1.7%
JJA−4.7%+0.1%−10.4%+1.9%−4.6%+0.1%−9.1%+1.7%
SON−4.5%−0.4%−8.5%+0.4%−4.3%−0.4%−10.5%+0.2%
Average−4.4%−0.3%−9.2%+0.1%−4.2%−0.2%−9.5%+0.1%
DJF0.0%−0.1%−1.0%0.0%0.0%−0.1%−0.4%0.0%AS-2
MAM0.0%−0.1%0.0%−0.1%0.0%−0.1%−0.8%−0.1%
JJA−0.1%−0.1%−1.3%+0.2%−0.1%0.0%−0.7%+0.2%
SON−0.1%−0.1%−0.6%0.0%−0.1%−0.1%−0.9%0.0%
Average−0.1%−0.1%−0.7%0.0%0.0%−0.1%−0.8%0.0%
DJF0.0%0.0%−72.4%0.0%0.0%0.0%−71.5%0.0%AS-3
MAM0.0%0.0%−70.4%0.0%0.0%0.0%−71.3%0.0%
JJA0.0%0.0%−71.6%0.0%0.0%0.0%−71.4%0.0%
SON0.0%0.0%−72.5%0.0%0.0%0.0%−72.2%0.0%
Average0.0%0.0%−71.8%0.0%0.0%0.0%−71.7%0.0%
DJF+0.1%+2.9%−25.3%−0.2%0.0%+2.2%−26.2%−0.3%AS-4
MAM0.0%+2.6%−8.5%−0.1%0.0%+1.8%−9.8%−0.1%
JJA+0.2%+2.6%−20.3%−1.4%+0.2%+1.7%−20.9%−1.4%
SON+0.3%+3.8%−46.2%−3.7%+0.2%+2.6%−42.7%−3.6%
Average+0.2%+3.0%−27.3%−1.2%+0.1%+2.1%−28.2%−1.2%
DJF−0.3%−2.4%−3.1%−0.5%−0.1%−1.4%−5.1%−0.5%AS-5
MAM−0.4%−2.1%−9.6%+1.1%−0.1%−1.1%−9.8%+0.9%
JJA−1.1%−2.9%−19.5%+1.5%−0.6%−1.6%−18.8%+1.2%
SON−0.7%−2.8%−12.2%0.0%−0.3%−1.4%−10.5%0.0%
Average−0.6%−2.5%−11.2%+1.0%−0.3%−1.4%−11.2%+0.8%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Badora, D.; Wawer, R.; Nieróbca, A.; Król-Badziak, A.; Kozyra, J.; Jurga, B.; Nowocień, E. Simulating the Effects of Agricultural Adaptation Practices onto the Soil Water Content in Future Climate Using SWAT Model on Upland Bystra River Catchment. Water 2022, 14, 2288. https://doi.org/10.3390/w14152288

AMA Style

Badora D, Wawer R, Nieróbca A, Król-Badziak A, Kozyra J, Jurga B, Nowocień E. Simulating the Effects of Agricultural Adaptation Practices onto the Soil Water Content in Future Climate Using SWAT Model on Upland Bystra River Catchment. Water. 2022; 14(15):2288. https://doi.org/10.3390/w14152288

Chicago/Turabian Style

Badora, Damian, Rafał Wawer, Anna Nieróbca, Aleksandra Król-Badziak, Jerzy Kozyra, Beata Jurga, and Eugeniusz Nowocień. 2022. "Simulating the Effects of Agricultural Adaptation Practices onto the Soil Water Content in Future Climate Using SWAT Model on Upland Bystra River Catchment" Water 14, no. 15: 2288. https://doi.org/10.3390/w14152288

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop