National flood discharge mapping in Austria

This article presents the approach and the results of a study in which 30, 100 and 200 year return period flood discharges were estimated for 26,000 km of Austrian streams. Three guiding principles were adopted: combination of automatic methods and manual assessments by hydrologists to allow speedy processing and account for the local hydrological situation; combination of various sources of information including flood peak samples, rainfall data, runoff coefficients and historical flood data; and involvement of the Hydrographic Services to increase the accuracy and enhance the acceptance of results. The flood discharges for ungauged catchments were estimated by the Top-kriging approach with manual adjustment to the local flood characteristics. The adopted combination approach proved to be very efficient both in terms of the project time required and in terms of the accuracy and acceptability of the estimated flood discharges of given return periods.     

A review on missing hydrological data processing

Like almost all fields of science, hydrology has benefited to a large extent from the tremendous improvements in scientific instruments that are able to collect long-time data series and an increase in available computational power and storage capabilities over the last decades. Many model applications and statistical analyses (e.g., extreme value analysis) are based on these time series. Consequently, the quality and the completeness of these time series are essential. Preprocessing of raw data sets by filling data gaps is thus a necessary procedure. Several interpolation techniques with different complexity are available ranging from rather simple to extremely challenging approaches. In this paper, various imputation methods available to the hydrological researchers are reviewed with regard to their suitability for filling gaps in the context of solving hydrological questions. The methodological approaches include arithmetic mean imputation, principal component analysis, regression-based methods and multiple imputation methods. In particular, autoregressive conditional heteroscedasticity (ARCH) models which originate from finance and econometrics will be discussed regarding their applicability to data series characterized by non-constant volatility and heteroscedasticity in hydrological contexts. The review shows that methodological advances driven by other fields of research bear relevance for a more intensive use of these methods in hydrology. Up to now, the hydrological community has paid little attention to the imputation ability of time series models in general and ARCH models in particular.     

Multi-angle backscatter classification and sub-bottom profiling for improved seafloor characterization

This study applies three classification methods exploiting the angular dependence of acoustic seafloor backscatter along with high resolution sub-bottom profiling for seafloor sediment characterization in the Eckernförde Bay, Baltic Sea Germany. This area is well suited for acoustic backscatter studies due to its shallowness, its smooth bathymetry and the presence of a wide range of sediment types. Backscatter data were acquired using a Seabeam1180 (180 kHz) multibeam echosounder and sub-bottom profiler data were recorded using a SES-2000 parametric sonar transmitting 6 and 12 kHz. The high density of seafloor soundings allowed extracting backscatter layers for five beam angles over a large part of the surveyed area. A Bayesian probability method was employed for sediment classification based on the backscatter variability at a single incidence angle, whereas Maximum Likelihood Classification (MLC) and Principal Components Analysis (PCA) were applied to the multi-angle layers. The Bayesian approach was used for identifying the optimum number of acoustic classes because cluster validation is carried out prior to class assignment and class outputs are ordinal categorical values. The method is based on the principle that backscatter values from a single incidence angle express a normal distribution for a particular sediment type. The resulting Bayesian classes were well correlated to median grain sizes and the percentage of coarse material. The MLC method uses angular response information from five layers of training areas extracted from the Bayesian classification map. The subsequent PCA analysis is based on the transformation of these five layers into two principal components that comprise most of the data variability. These principal components were clustered in five classes after running an external cluster validation test. In general both methods MLC and PCA, separated the various sediment types effectively, showing good agreement (kappa >0.7) with the Bayesian approach which also correlates well with ground truth data (r²?>?0.7). In addition, sub-bottom data were used in conjunction with the Bayesian classification results to characterize acoustic classes with respect to their geological and stratigraphic interpretation. The joined interpretation of seafloor and sub-seafloor data sets proved to be an efficient approach for a better understanding of seafloor backscatter patchiness and to discriminate acoustically similar classes in different geological/bathymetric settings.     

Subglacial tunnel valleys in the Alpine foreland: an example from Bern, Switzerland

The morphology of the Alpine and adjacent landscapes is directly related to glacial erosion and associated sediment transport. Here we report the effects of glacio-hydrologic erosion on bedrock topography in the Swiss Plateau. Specifically, we identify the presence of subsurface valleys beneath the city of Bern and discuss their genesis. Stratigraphic investigations of more than 4,000 borehole data within a 430 km²-large area reveal the presence of a network of >200 m-deep and 1,000 m-wide valleys. They are flat floored with steep sided walls and are filled by Quaternary glacial deposits. The central valley beneath Bern is straight and oriented towards the NNW, with valley flanks more than 20° steep. The valley bottom has an irregular undulating profile along the thalweg, with differences between sills and hollows higher than 50–100 m over a reach of 4 km length. Approximately 500 m high bedrock highlands flank the valley network. The highlands are dissected by up to 80 m-deep and 500 m-broad hanging valleys that currently drain away from the axis of the main valley. We interpret the valleys beneath the city of Bern to be a tunnel valley network which originated from subglacial erosion by melt water. The highland valleys served as proglacial meltwater paths and are hanging with respect to the trunk system, indicating that these incipient highland systems as well as the main gorge beneath Bern formed by glacial melt water under pressure.     

Floodplain hydrology of the Mekong Delta,Vietnam

The Mekong Delta is one of the largest and most intensively used estuaries in the world. Each year it witnesses widespread flooding which is both the basis of the livelihood for more than 17 million people but also the major hazard. Therefore, a thorough understanding of the hydrologic and hydraulic features is urgently required for various planning purposes. While the general causes and characteristics of the annual floods are understood, the inundation dynamics in the floodplains in Vietnam which are highly controlled by dikes and other control structures have not been investigated in depth. Especially, quantitative analyses are lacking, mainly due to scarce data about the inundation processes in the floodplains. Therefore, a comprehensive monitoring scheme for channel and floodplain inundation was established in a study area in the Plain of Reeds in the northeastern part of the Vietnamese Delta. This in situ data collection was complemented by a series of high‐resolution inundation maps derived from the TerraSAR‐X satellite for the flood seasons 2008 and 2009. Hence, the inundation dynamics in the channels and floodplains, and the interaction between channels and floodplains, could be quantified for the first time. The study identifies the strong human interference which is governed by flood protection levels, cropping patterns and communal water management. In addition, we examine the tidal influence on the inundation in various parts of the Delta, since it is expected that climate change‐induced sea level rise will increase the tidal contribution to floodplain inundation. Copyright © 2011 John Wiley & Sons, Ltd.     

High spatial-resolution monitoring to investigate nitrate export and its drivers in a mesoscale river catchment along an anthropogenic land-cover gradient

Nitrate monitoring is commonly conducted with low-spatial resolution, only at the outlet or at a small number of selected locations. As a result, the information about spatial variations in nitrate export and its drivers is scarce. In this study, we present results of high-spatial resolution monitoring conducted between 2012 and 2017 in 65 sub-catchments in an Alpine mesoscale river catchment characterized by a land-use gradient. We combined stable isotope techniques with Bayesian mixing models and geostatistical methods to investigate nitrate export and its main drivers, namely, microbial N turnover processes, land use and hydrological conditions. In the investigated sub-catchments, mean values of NO₃⁻ concentrations and its isotope signatures (δ¹⁵N_NO3 and δ¹⁸O_NO3) varied from 2.6 to 26.7 mg L⁻¹, from −1.3‰ to 13.1‰, and from −0.4‰ to 10.1‰, respectively. In this study, land use was an important driver for nitrate export. Very strong and strong positive correlations were found between percentages of agricultural land cover and δ¹⁵N_NO3, and NO₃⁻ concentration, respectively. Mean proportional contributions of NO₃⁻ sources varied spatially and seasonally, and followed land-use patterns. The mean contribution of manure and sewage was much higher in the catchments characterized by a high percentage of agricultural and urban land cover comparing to forested sub-catchments. Specific NO₃⁻ loads were strongly correlated with specific discharge and moderately correlated with NO₃⁻ concentrations. The nitrate isotope and concentration analysis results suggest that nitrate from external sources is stored and accumulated in soil storage pools. Nitrification of reduced nitrogen species in those pools plays the most important role for the N-dynamics in the Erlauf river catchment. Consequently, nitrification of reduced N sources was the main nitrate source except for a number of sub-catchments dominated by agricultural land use. In the Erlauf catchment, denitrification plays only a minor role in controlling NO₃⁻ export on a regional scale.     

Identification of groundwater mean transit times of precipitation and riverbank infiltration by two‐component lumped parameter models

Groundwater transit time is an essential hydrologic metric for groundwater resources management. However, especially in tropical environments, studies on the transit time distribution (TTD) of groundwater infiltration and its corresponding mean transit time (mTT) have been extremely limited due to data sparsity. In this study, we primarily use stable isotopes to examine the TTDs and their mTTs of both vertical and horizontal infiltration at a riverbank infiltration area in the Vietnamese Mekong Delta (VMD), representative of the tropical climate in Asian monsoon regions. Precipitation, river water, groundwater, and local ponding surface water were sampled for 3 to 9 years and analysed for stable isotopes (δ¹⁸O and δ²H), providing a unique data set of stable isotope records for a tropical region. We quantified the contribution that the two sources contributed to the local shallow groundwater by a novel concept of two‐component lumped parameter models (LPMs) that are solved using δ¹⁸O records. The study illustrates that two‐component LPMs, in conjunction with hydrological and isotopic measurements, are able to identify subsurface flow conditions and water mixing at riverbank infiltration systems. However, the predictive skill and the reliability of the models decrease for locations farther from the river, where recharge by precipitation dominates, and a low‐permeable aquitard layer above the highly permeable aquifer is present. This specific setting impairs the identifiability of model parameters. For river infiltration, short mTTs (<40 weeks) were determined for sites closer to the river (<200 m), whereas for the precipitation infiltration, the mTTs were longer (>80 weeks) and independent of the distance to the river. The results not only enhance the understanding of the groundwater recharge dynamics in the VMD but also suggest that the highly complex mechanisms of surface–groundwater interaction can be conceptualized by exploiting two‐component LPMs in general. The model concept could thus be a powerful tool for better understanding both the hydrological functioning of mixing processes and the movement of different water components in riverbank infiltration systems.     

Gravity changes and natural gas extraction in Groningen

Gravity observations made over the Groningen gasfield from 1978 onwards are available. The existing gravity data were reanalysed and a new survey was carried out. The observed gravity changes obtained during four surveys, spanning an 18-year period, were compared with the gravity effect due to natural gas extraction computed from the reservoir model and the production data. The random error in the gravity values is small enough to detect the effect of gas extraction after a few years (up to 2 μgal/year). The trends obtained from observed and calculated gravity changes agree well within their expected error margins after statistical data snooping. Due to the inadequate measurement set-up of the previous surveys, systematic errors present in the gravity data hampered the use of this data for refinement of the reservoir and the production models. However, with improved gravity monitoring, in particular a very systematic network set-up and well-calibrated instruments, this should now be feasible. By using the network analysis from the 1996 measurements for the planning of a network in space and time, most of the current problems can be avoided.