Hydrological drought is critical from both water management and ecological perspectives. Depending on its hydrological and physical features, the resilience level of a catchment to groundwater drought can differ from that of meteorological drought. This study presents a comparison of hydrological and meteorological drought indices based on groundwater levels from 1987 to 2018. A small catchment area in Croatia, consisting of two sub-catchments with a continental climate and minimum land-use changes during the observed period, was studied. The first analysis was made on a comparison of standardized precipitation index (SPI) and standardized precipitation evapotranspiration index (SPEI). The results showed their very high correlation. The correlation between the standardized precipitation index (SPI) and standardized groundwater index (SGI) of different time scales (1, 3, 6, 12, 24 and 48 months) showed different values, but had the highest value in the longest time scale, 48 months, for all observation wells. Nevertheless, the behavior of the SPI and groundwater levels (GW) correlation showed results more related to physical catchment characteristics. The results showed that groundwater drought indices, such as SGI, should be applied judiciously because of their sensitivity to geographical, geomorphological, and topographical catchment characteristics, even in small catchment areas.

Droughts are among the major natural hazards that are spreading to many parts of the world, with huge multi-dimensional impacts. An extensive analysis of drought phenomenon is presented for continental Croatia based on a meteorological E-OBS gridded dataset (0.25° × 0.25°), within the period of 1950–2022. The drought events were characterized by the Standardized Precipitation Evapotranspiration Index (SPEI), applied to different time-scales (6 and 12 months), in order to describe the subannual and annual variability of drought. The spatiotemporal patterns of drought are obtained through principal component analysis (PCA) and K-means clustering (KMC) applied to the SPEI field. An areal drought evolution analysis and the changes in the frequency of occurrence of the periods under drought conditions were achieved using a kernel occurrence rate estimator (KORE). The modified Mann–Kendall (MMK) test, coupled with the Sen’s slope estimator test, are applied to the SPEI series in order to quantify the drought trends throughout the country. According to the history drought events and considering the different morphoclimatic characteristics of the study area, the results showed that Croatia could be divided into three different and spatially well-defined regions with specific temporal and spatial characteristics of droughts (central northern, eastern and southern regions). A manifest increase is shown in the percentage of area affected by drought, as well as in the yearly drought occurrences rates, in both central northern and eastern regions, and an evident decrease is shown in the southern region for both 6- and 12-month SPEI time-scales. In the observation of the drought’s temporal characteristics, it was found that downward trends expressing increasing drought severities were strongly significant in northern and eastern regions, while a few significant upward trends were seen in the southern region. From this study, it is possible to obtain a broader view of the historical behaviour of droughts in Croatia, with the results providing useful support for drought risk assessment and decision-making processes.

The majority of central Europe has a transitional climate type as a result of interactions between maritime and continental climates. This study focuses on the appearance and severity of drought in continental Croatia, which is part of the transitional climate area. It is situated between 15° E and 19° E. The altitude declines from west (167 m a.s.l) to east (88 m a.s.l.). The time period analysed is 1981–2018. Air temperature and precipitation data series from 13 meteorological stations were analysed. The analysis was done on an annual basis to define the spatio-temporal variability in air temperature and precipitation and their impact on drought episodes using the standardised evapotranspiration precipitation index. Different statistical methods (e.g., the nonparametric Mann–Kendall test and agglomerative hierarchical clustering) were used to examine the trend homogeneity of the analysed region. The analysis indicated inhomogeneity across the study area in terms of what significantly impacted the occurrence and severity of droughts. Drought occurrence is influenced more strongly by increasing trends in air temperature as compared with increasing or decreasing precipitation trends. The probability of severe drought occurrence was estimated using a copula function, and the results demonstrated that areas with higher precipitation could be more exposed to drought. Furthermore, the results demonstrated the impacts of specific regional characteristics on drought occurrence, severity, and duration, which indicates that small-scale research on droughts is more reliable.

Droughts are among the major natural hazards that are spreading to many parts of the world with huge multi-dimensional impacts. For Continental Croatia an extensive analysis of drought phenomenon is presented based on meteorological E-OBS gridded dataset (0.25° x 0.25°), within the period of 1950 to 2022. The drought events were characterized by the Standardized Precipita-tion Evapotranspiration Index (SPEI) applied to different time-scales (6 and 12 months) in order to describe the subannual and annual variability of drought. The spatiotemporal patterns of drought are obtained through principal component analysis (PCA) and K-means clustering (KMC) applied to the SPEI field. An areal drought evolution analysis and the changes in the frequency of occur-rence of the periods under drought conditions were achieved by using a kernel occurrence rate es-timator (KORE). The Modified Mann-Kendall (MMK) test coupled with the Sen’s slope estimator test are applied to the SPEI series in order to quantify the drought trends throughout the country. According to drought events history and considering the different morphoclimatic characteristics of the study area the results showed that Croatia could be divided into three different and spatial-ly well-defined regions with specific temporal and spatial characteristics of droughts (Central North, Eastern and Southern regions). It is shown that there is a manifest increase in the percent-age of area affected by drought as well as in the yearly drought occurrences rates in both Central North and Eastern regions and an evident decrease in the Southern region for both 6 and 12 month SPEI time-scales. In the observation of the drought temporal characteristics, it was found that downward trends expressing increasing drought severities were strongly significant in North and Eastern regions while a few significant upward trends were seen in Southern region. From this study, it is possible to obtain a broader view of the historical behaviour of droughts in Croatia with results providing useful support for drought risk assessment and decision-making processes.

The primary processing of high-resolution precipitation records (5 min and shorter) is crucial for constructing dimensionless design hyetographs and identifying design-critical precipitation scenarios for urban drainage systems. A key step in this process is separating continuous precipitation records into individual precipitation events, typically based on minimum inter-event time (MIT) and precipitation amount thresholds. This separation directly influences the subsequent analysis steps and the accuracy of the design hyetographs. Building upon this foundation, this study systematically analyses how different MIT determination methods influence the construction of dimensionless Huff hyetographs in a moderately humid continental climate. Three approaches for defining MIT were examined: a fixed MIT method (1–12 h), an autocorrelation-based method (AC), and a kernel density estimation approach (KDE). The analysis also considers the effects of minimum precipitation thresholds (P = 1, 3, and 5 mm) and precipitation duration classes (all durations and short-duration events with $T\le 2$ h), utilising a continuous 10-year series of 5 min precipitation data. The results demonstrate that the choice of MIT substantially affects the identified precipitation events, duration, total amount, and the median Huff curve’s shape, especially for precipitation types with early and late maximum intensity. Specifically, increasing MIT values produces longer and deeper events with steeper Huff curves, while precipitation thresholds mainly filter weaker events rather than impacting peak intensities. The AC method yields results similar to larger fixed MIT values (≈6–9 h), whereas the KDE method corresponds to shorter separations (≈1–3 h). To unify the assessment of design relevance, a composite design index combining Huff curve slope and short-term peak intensities was introduced. Analysis shows that short-duration convective precipitation with an early maximum is the most critical design scenario. However, late-maximum events (events in which peak intensity occurs in the fourth quartile of storm duration, Type 4) can become equally critical when longer MIT values or autocorrelation-based separation are applied. These findings underscore the importance of a transparent and methodologically consistent definition of precipitation event separation criteria when using dimensionless hyetographs in urban drainage design.