Future climate change is projected to have significant impact on water resources availability and quality in many parts of the world. The objective of this paper is to assess the effect of projected climate change on water quantity and quality in two lowland catchments (the Upper Narew and the Barycz) in Poland in two future periods (near future: 2021–2050, and far future: 2071– 2100). The hydrological model SWAT was driven by climate forcing data from an ensemble of nine bias-corrected General Circulation Models—Regional Climate Models (GCM-RCM) runs based on the Coordinated Downscaling Experiment—European Domain (EURO-CORDEX). Hydrological response to climate warming and wetter conditions (particularly in winter and spring) in both catchments includes: lower snowmelt, increased percolation and baseflow and higher runoff. Seasonal differences in the response between catchments can be explained by their properties (e.g., different thermal conditions and soil permeability). Projections suggest only moderate increases in sediment loss, occurring mainly in summer and winter. A sharper increase is projected in both catchments for TN losses, especially in the Barycz catchment characterized by a more intensive agriculture. The signal of change in annual TP losses is blurred by climate model uncertainty in the Barycz catchment, whereas a weak and uncertain increase is projected in the Upper Narew catchment.

Ground-based precipitation data are still the dominant input type for hydrological models. Spatial variability in precipitation can be represented by spatially interpolating gauge data using various techniques. In this study, the effect of daily precipitation interpolation methods on discharge simulations using the semi-distributed SWAT (Soil and Water Assessment Tool) model over a 30-year period is examined. The study was carried out in 11 meso-scale (119–3935 km²) sub-catchments lying in the Sulejów reservoir catchment in central Poland. Four methods were tested: the default SWAT method (Def) based on the Nearest Neighbour technique, Thiessen Polygons (TP), Inverse Distance Weighted (IDW) and Ordinary Kriging (OK). =The evaluation of methods was performed using a semi-automated calibration program SUFI-2 (Sequential Uncertainty Fitting Procedure Version 2) with two objective functions: Nash-Sutcliffe Efficiency (NSE) and the adjusted R² coefficient (bR²). The results show that: (1) the most complex OK method outperformed other methods in terms of NSE; and (2) OK, IDW, and TP outperformed Def in terms of bR². The median difference in daily/monthly NSE between OK and Def/TP/IDW calculated across all catchments ranged between 0.05 and 0.15, while the median difference between TP/IDW/OK and Def ranged between 0.05 and 0.07. The differences between pairs of interpolation methods were, however, spatially variable and a part of this variability was attributed to catchment properties: catchments characterised by low station density and low coefficient of variation of daily flows experienced more pronounced improvement resulting from using interpolation methods. Methods providing higher precipitation estimates often resulted in a better model performance. The implication from this study is that appropriate consideration of spatial precipitation variability (often neglected by model users) that can be achieved using relatively simple interpolation methods can significantly improve the reliability of model simulations.

There is considerable concern that the water resources of Central and Eastern Europe region can be adversely affected by climate change. Projections of future water balance and streamflow conditions can be obtained by forcing hydrological models with the output from climate models. In this study, we employed the SWAT hydrological model driven with an ensemble of nine bias-corrected EURO-CORDEX climate simulations to generate future hydrological projections for the Vistula and Odra basins in two future horizons (2024–2050 and 2074–2100) under two Representative Concentration Pathways (RCPs). The data set consists of three parts: (1) model inputs; (2) raw model outputs; (3) aggregated model outputs. The first one allows the users to reproduce the outputs or to create the new ones. The second one contains the simulated time series of 10 variables simulated by SWAT: precipitation, snow melt, potential evapotranspiration, actual evapotranspiration, soil water content, percolation, surface runoff, baseflow, water yield and streamflow. The third one consists of the multi-model ensemble statistics of the relative changes in mean seasonal and annual variables developed in a GIS format. The data set should be of interest of climate impact scientists, water managers and water-sector policy makers. In any case, it should be noted that projections included in this data set are associated with high uncertainties explained in this data descriptor paper.