Spatial 3-D information on soil hydraulic properties for areas larger than plot scale is usually derived using indirect methods such as pedotransfer functions (PTFs) due to the lack of measured information on them. PTFs describe the relationship between the desired soil hydraulic parameter and easily available soil properties based on a soil hydraulic reference dataset. Soil hydraulic properties of a catchment or region can be calculated by applying PTFs on available soil maps. Our aim was to analyse the performance of (i) indirect (using PTFs) and (ii) direct (geostatistical) mapping methods to derive 3-D soil hydraulic properties. The study was performed on the Balaton catchment area in Hungary, where density of measured soil hydraulic data fulfils the requirements of geostatistical methods. Maps of saturated water content (0 cm matric potential), field capacity (−330 cm matric potential) and wilting point (−15 000 cm matric potential) for 0–30, 30–60 and 60–90 cm soil depth were prepared. PTFs were derived using the random forest method on the whole Hungarian soil hydraulic dataset, which includes soil chemical, physical, taxonomical and hydraulic properties of some 12 000 samples complemented with information on topography, climate, parent material, vegetation and land use. As a direct and thus geostatistical method, random forest combined with kriging (RFK) was applied to 359 soil profiles located in the Balaton catchment area. There were no significant differences between the direct and indirect methods in six out of nine maps having root-mean-square-error values between 0.052 and 0.074 cm³ cm⁻³, which is in accordance with the internationally accepted performance of hydraulic PTFs. The PTF-based mapping method performed significantly better than the RFK for the saturated water content at 30–60 and 60–90 cm soil depth; in the case of wilting point the RFK outperformed the PTFs at 60–90 cm depth. Differences between the PTF-based and RFK mapped values are less than 0.025 cm³ cm⁻³ for 65 %–86 % of the catchment. In RFK, the uncertainty of input environmental covariate layers is less influential on the mapped values, which is preferable. In the PTF-based method the uncertainty of mapping soil hydraulic properties is less computationally intensive. Detailed comparisons of maps derived from the PTF-based method and the RFK are presented in this paper.

A dataset of annual freshwater ice phenology was compiled for the largest river (Danube) and the largest lake (Lake Balaton) in eastern–central Europe, extending regular river and lake ice monitoring data through the use of historical observations and documentary records dating back to AD 1774 and AD 1885, respectively. What becomes clear is that the dates of the first appearance of ice and freeze-up have shifted, arriving 12–30 and 4–13 days later, respectively, per 100 years. Break-up and ice-off have shifted to earlier dates by 7–13 and 9–27 days/100 years, except on Lake Balaton, where the date of break-up has not changed significantly. The datasets represent a resource for (paleo)climatological research thanks to the strong, physically determined link between water and air temperature and the occurrence of freshwater ice phenomena. The derived centennial records of freshwater cryophenology for the Danube and Balaton are readily available for detailed analysis of the temporal trends, large-scale spatial comparison, or other climatological purposes. The derived dataset is publicly available via PANGAEA at https://doi.org/10.1594/PANGAEA.881056.