Bridge piers on large rivers are often protected from scouring using launchable stone, such as a riprap sloping structure. While such scour countermeasures are effective for pier protection, they significantly alter flow conditions in the bridge opening by overtopping flow and flow contraction, deflecting the formation of the scour hole downstream and exposing the downstream riverbed to additional scour. This paper provides a comprehensive and relevant review of bridge scour estimation methods for piers with a riprap sloping structure installed as a scour countermeasure. Research on empirical methods for bridge scour estimation is reviewed and analyzed with formulae used for comparable structures—complex pier formulae and formulae for river training structures. A summary of relevant formulae applicable to piers with installed scour countermeasures is provided, as well as a discussion on the possible future research directions that could contribute to the field.

Morphodynamic changes in the riverbed may be accelerated by the climate change-induced effects, mostly through the increase of the frequency of extreme climatic events such as floods. This can lead to scouring of the riverbed around the bridge substructure and consequently reduces its overall stability. In order to better understand hydromorphological processes at the local scale, the influence of floods on bridge scour requires a detailed analysis of several interacting flood hydrograph characteristics. This paper presents a multivariate analysis of the annual maximum (AM) flood discharge data at four gauging stations on the Drava and Sava Rivers in Croatia (Europe). As part of the hydrograph analysis, multiple baseflow separation methods were tested. Flood volumes and durations were derived after extracting the baseflow from measured discharge data. Suitable marginal distribution functions were fitted to the peak discharge (Q), flood volume (V) and duration (D) data. Bivariate copula analyses were conducted for the next pairs: peak discharge and volume (Q–V), hydrograph volume and duration (V–D) and peak discharge and hydrograph duration (Q–D). The results of the bivariate copula analyses were used to derive joint return periods for different flood variable combinations, which may serve as a preliminary analysis for the pilot bridges of the R3PEAT project where the aim is to investigate the influences on the riverbed erosion around bridges with installed scour countermeasures. Hence, a design hydrograph was derived that could be used as input data in the hydraulic model for the investigation of the bridge scour dynamics within the project and a preliminary methodology is proposed to be applied. The results indicate that bivariate frequency analysis can be very sensitive to the selected baseflow separation methodology. Therefore, future studies should test multiple baseflow separation methods and visually inspect the performance.

Rainfed agriculture is dependent on rainfall and runoff patterns, especially in lowland areas that rely on pumping operation to remove excess water from the drainage network. Polder areas are extremely vulnerable to saltwater intrusion and subsequent soil salinization driven by rising sea levels and accelerated by climate change. The aim of this paper is to reconstruct the recharge and discharge pattern in the Vidrice polder, a drainage canal network within the Neretva River Delta agroecosystem used to collect the surface and subsurface runoff from the agricultural land and saltwater infiltration through the aquifer. Water regime data are collected over an 18-month period of real-time monitoring at 15 min intervals on three stations along the primary drainage canal and one station at the secondary canal. Analysis of water level flashiness in the Vidrice polder using the Richards-Baker flashiness index (R-B_index) indicates that daily pumping of water infiltrated in the canal network is sub-optimal: discharge fluctuates significantly more than recharge, by 46% on average, resulting in unnecessary lowering of the water level in the drainage network. The results show that the correlation between the intensive rainfall events (>10 mm/day) and the recharge rates can be used to modify the daily pumping operation and maintain high freshwater levels in the canal network to increase the resistance to infiltration and reduce saltwater intrusion into the polder.

The estimation of baseflow is one of the essential tasks in water resources management and hydrologic research to assess the impacts of climate change and to describe and predict flood events based on the flood hydrograph characteristics (peak flow, duration and volume). Several methods have been developed to separate baseflow from direct flow, and in recent years they have been automated through the use of available R packages. In this work R programming language packages “EcoHydRology” and “lfstat” were used to separate baseflow from direct flow on the historical daily discharge time series of the several gauging stations on the two large lowland rivers in Croatia: the Sava River and the Drava River. The aim of this study is to determine the appropriate baseflow separation method for gauging stations on Sava River and Drava River in order to evaluate the baseflow separation method for future multivariate analysis of flood events under the R3PEAT project (www.grad.hr/r3peat) that explores pier scour development next to the bridges crossing large rivers in Croatia with installed scour countermeasures.

Determination of flood magnitude and shape characteristics are necessary to provide a more complete assessment of flood severity and its impact in scour development analysis. Our recent research has focused on a joint distribution approach to account for the multivariate nature of flood characteristics, resulting in probability of occurrence of different pairs of flood variables: flood peak (Q), volume (V) and duration (D). To extend the results of this research, a method for deriving a design hydrograph is applied to the study area by using the typical hydrograph method. As it is recommended to test multiple scenarios in a scour analysis, different typical flood hydrographs were selected at several gauging stations on the Sava and Drava rivers in Croatia and multiplied by the design discharge values. The aim of this study is to complement the ongoing research of the relationship between climate change indicators, flood wave characteristics and scour development next to the bridges crossing large rivers in Croatia with installed scour countermeasures by preparing hydrological input data for a hydraulic scour analysis.