@misc{AlHalbouniHolohanTaherietal.2018, author = {Al-Halbouni, Djamil and Holohan, Eoghan P. and Taheri, Abbas and Sch{\"o}pfer, Martin P. J. and Emam, Sacha and Dahm, Torsten}, title = {Geomechanical modelling of sinkhole development using distinct elements}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1061}, issn = {1866-8372}, doi = {10.25932/publishup-46843}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-468435}, pages = {35}, year = {2018}, abstract = {Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth ∕ diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth ∕ diameter values in each material type may partly reflect sinkhole growth trends.}, language = {en} } @article{AmourMuttiChristetal.2012, author = {Amour, Frederic and Mutti, Maria and Christ, Nicolas and Immenhauser, Adrian and Agar, Susan M. and Benson, Gregory S. and Tomas, Sara and Alway, Robert and Kabiri, Lachen}, title = {Capturing and modelling metre-scale spatial facies heterogeneity in a Jurassic ramp setting (Central High Atlas, Morocco)}, series = {Sedimentology : the journal of the International Association of Sedimentologists}, volume = {59}, journal = {Sedimentology : the journal of the International Association of Sedimentologists}, number = {4}, publisher = {Wiley-Blackwell}, address = {Malden}, issn = {0037-0746}, doi = {10.1111/j.1365-3091.2011.01299.x}, pages = {1158 -- 1189}, year = {2012}, abstract = {Each simulation algorithm, including Truncated Gaussian Simulation, Sequential Indicator Simulation and Indicator Kriging is characterized by different operating modes, which variably influence the facies proportion, distribution and association of digital outcrop models, as shown in clastic sediments. A detailed study of carbonate heterogeneity is then crucial to understanding these differences and providing rules for carbonate modelling. Through a continuous exposure of Bajocian carbonate strata, a study window (320 m long, 190 m wide and 30 m thick) was investigated and metre-scale lithofacies heterogeneity was captured and modelled using closely-spaced sections. Ten lithofacies, deposited in a shallow-water carbonate-dominated ramp, were recognized and their dimensions and associations were documented. Field data, including height sections, were georeferenced and input into the model. Four models were built in the present study. Model A used all sections and Truncated Gaussian Simulation during the stochastic simulation. For the three other models, Model B was generated using Truncated Gaussian Simulation as for Model A, Model C was generated using Sequential Indicator Simulation and Model D was generated using Indicator Kriging. These three additional models were built by removing two out of eight sections from data input. The removal of sections allows direct insights on geological uncertainties at inter-well spacings by comparing modelled and described sections. Other quantitative and qualitative comparisons were carried out between models to understand the advantages/disadvantages of each algorithm. Model A is used as the base case. Indicator Kriging (Model D) simplifies the facies distribution by assigning continuous geological bodies of the most abundant lithofacies to each zone. Sequential Indicator Simulation (Model C) is confident to conserve facies proportion when geological heterogeneity is complex. The use of trend with Truncated Gaussian Simulation is a powerful tool for modelling well-defined spatial facies relationships. However, in shallow-water carbonate, facies can coexist and their association can change through time and space. The present study shows that the scale of modelling (depositional environment or lithofacies) involves specific simulation constraints on shallow-water carbonate modelling methods.}, language = {en} } @phdthesis{Amour2013, author = {Amour, Fr{\´e}d{\´e}ric}, title = {3-D modeling of shallow-water carbonate systems : a scale-dependent approach based on quantitative outcrop studies}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-66621}, school = {Universit{\"a}t Potsdam}, year = {2013}, abstract = {The study of outcrop modeling is located at the interface between two fields of expertise, Sedimentology and Computing Geoscience, which respectively investigates and simulates geological heterogeneity observed in the sedimentary record. During the last past years, modeling tools and techniques were constantly improved. In parallel, the study of Phanerozoic carbonate deposits emphasized the common occurrence of a random facies distribution along single depositional domain. Although both fields of expertise are intrinsically linked during outcrop simulation, their respective advances have not been combined in literature to enhance carbonate modeling studies. The present study re-examines the modeling strategy adapted to the simulation of shallow-water carbonate systems, based on a close relationship between field sedimentology and modeling capabilities. In the present study, the evaluation of three commonly used algorithms Truncated Gaussian Simulation (TGSim), Sequential Indicator Simulation (SISim), and Indicator Kriging (IK), were performed for the first time using visual and quantitative comparisons on an ideally suited carbonate outcrop. The results show that the heterogeneity of carbonate rocks cannot be fully simulated using one single algorithm. The operating mode of each algorithm involves capabilities as well as drawbacks that are not capable to match all field observations carried out across the modeling area. Two end members in the spectrum of carbonate depositional settings, a low-angle Jurassic ramp (High Atlas, Morocco) and a Triassic isolated platform (Dolomites, Italy), were investigated to obtain a complete overview of the geological heterogeneity in shallow-water carbonate systems. Field sedimentology and statistical analysis performed on the type, morphology, distribution, and association of carbonate bodies and combined with palaeodepositional reconstructions, emphasize similar results. At the basin scale (x 1 km), facies association, composed of facies recording similar depositional conditions, displays linear and ordered transitions between depositional domains. Contrarily, at the bedding scale (x 0.1 km), individual lithofacies type shows a mosaic-like distribution consisting of an arrangement of spatially independent lithofacies bodies along the depositional profile. The increase of spatial disorder from the basin to bedding scale results from the influence of autocyclic factors on the transport and deposition of carbonate sediments. Scale-dependent types of carbonate heterogeneity are linked with the evaluation of algorithms in order to establish a modeling strategy that considers both the sedimentary characteristics of the outcrop and the modeling capabilities. A surface-based modeling approach was used to model depositional sequences. Facies associations were populated using TGSim to preserve ordered trends between depositional domains. At the lithofacies scale, a fully stochastic approach with SISim was applied to simulate a mosaic-like lithofacies distribution. This new workflow is designed to improve the simulation of carbonate rocks, based on the modeling of each scale of heterogeneity individually. Contrarily to simulation methods applied in literature, the present study considers that the use of one single simulation technique is unlikely to correctly model the natural patterns and variability of carbonate rocks. The implementation of different techniques customized for each level of the stratigraphic hierarchy provides the essential computing flexibility to model carbonate systems. Closer feedback between advances carried out in the field of Sedimentology and Computing Geoscience should be promoted during future outcrop simulations for the enhancement of 3-D geological models.}, language = {en} } @misc{AngermannJackischAllroggenetal.2017, author = {Angermann, Lisa and Jackisch, Conrad and Allroggen, Niklas and Sprenger, Matthias and Zehe, Erwin and Tronicke, Jens and Weiler, Markus and Blume, Theresa}, title = {Form and function in hillslope hydrology}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {658}, issn = {1866-8372}, doi = {10.25932/publishup-41916}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-419161}, pages = {22}, year = {2017}, abstract = {The phrase form and function was established in architecture and biology and refers to the idea that form and functionality are closely correlated, influence each other, and co-evolve. We suggest transferring this idea to hydrological systems to separate and analyze their two main characteristics: their form, which is equivalent to the spatial structure and static properties, and their function, equivalent to internal responses and hydrological behavior. While this approach is not particularly new to hydrological field research, we want to employ this concept to explicitly pursue the question of what information is most advantageous to understand a hydrological system. We applied this concept to subsurface flow within a hillslope, with a methodological focus on function: we conducted observations during a natural storm event and followed this with a hillslope-scale irrigation experiment. The results are used to infer hydrological processes of the monitored system. Based on these findings, the explanatory power and conclusiveness of the data are discussed. The measurements included basic hydrological monitoring methods, like piezometers, soil moisture, and discharge measurements. These were accompanied by isotope sampling and a novel application of 2-D time-lapse GPR (ground-penetrating radar). The main finding regarding the processes in the hillslope was that preferential flow paths were established quickly, despite unsaturated conditions. These flow paths also caused a detectable signal in the catchment response following a natural rainfall event, showing that these processes are relevant also at the catchment scale. Thus, we conclude that response observations (dynamics and patterns, i.e., indicators of function) were well suited to describing processes at the observational scale. Especially the use of 2-D time-lapse GPR measurements, providing detailed subsurface response patterns, as well as the combination of stream-centered and hillslope-centered approaches, allowed us to link processes and put them in a larger context. Transfer to other scales beyond observational scale and generalizations, however, rely on the knowledge of structures (form) and remain speculative. The complementary approach with a methodological focus on form (i.e., structure exploration) is presented and discussed in the companion paper by Jackisch et al. (2017).}, language = {en} } @misc{NiedlBerensteinBeta2016, author = {Niedl, Robert Raimund and Berenstein, Igal and Beta, Carsten}, title = {How imperfect mixing and differential diffusion accelerate the rate of nonlinear reactions in microfluidic channels}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-95810}, pages = {6451 -- 6457}, year = {2016}, abstract = {In this paper, we show experimentally that inside a microfluidic device, where the reactants are segregated, the reaction rate of an autocatalytic clock reaction is accelerated in comparison to the case where all the reactants are well mixed. We also find that, when mixing is enhanced inside the microfluidic device by introducing obstacles into the flow, the clock reaction becomes slower in comparison to the device where mixing is less efficient. Based on numerical simulations, we show that this effect can be explained by the interplay of nonlinear reaction kinetics (cubic autocatalysis) and differential diffusion, where the autocatalytic species diffuses slower than the substrate.}, language = {en} } @article{NiedlBerensteinBeta2016, author = {Niedl, Robert Raimund and Berenstein, Igal and Beta, Carsten}, title = {How imperfect mixing and differential diffusion accelerate the rate of nonlinear reactions in microfluidic channels}, series = {Physical chemistry, chemical physics : PCCP ; a journal of European Chemical Societies}, volume = {18}, journal = {Physical chemistry, chemical physics : PCCP ; a journal of European Chemical Societies}, publisher = {Royal Society of Chemistry}, address = {Cambridge}, issn = {1463-9076}, doi = {10.1039/c6cp00224b}, pages = {6451 -- 6457}, year = {2016}, abstract = {In this paper, we show experimentally that inside a microfluidic device, where the reactants are segregated, the reaction rate of an autocatalytic clock reaction is accelerated in comparison to the case where all the reactants are well mixed. We also find that, when mixing is enhanced inside the microfluidic device by introducing obstacles into the flow, the clock reaction becomes slower in comparison to the device where mixing is less efficient. Based on numerical simulations, we show that this effect can be explained by the interplay of nonlinear reaction kinetics (cubic autocatalysis) and differential diffusion, where the autocatalytic species diffuses slower than the substrate.}, language = {en} } @misc{NordBoudevillainBerneetal.2017, author = {Nord, Guillaume and Boudevillain, Brice and Berne, Alexis and Branger, Flora and Braud, Isabelle and Dramais, Guillaume and G{\´e}rard, Simon and Le Coz, J{\´e}r{\^o}me and Lego{\^u}t, C{\´e}dric and Molini{\´e}, Gilles and Van Baelen, Joel and Vandervaere, Jean-Pierre and Andrieu, Julien and Aubert, Coralie and Calianno, Martin and Delrieu, Guy and Grazioli, Jacopo and Hachani, Sahar and Horner, Ivan and Huza, Jessica and Le Boursicaud, Rapha{\"e}l and Raupach, Timothy H. and Teuling, Adriaan J. and Uber, Magdalena and Vincendon, B{\´e}atrice and Wijbrans, Annette}, title = {A high space-time resolution dataset linking meteorological forcing and hydro-sedimentary response in a mesoscale Mediterranean catchment (Auzon) of the Ard{\`e}che region, France}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {671}, issn = {1866-8372}, doi = {10.25932/publishup-41912}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-419127}, pages = {29}, year = {2017}, abstract = {A comprehensive hydrometeorological dataset is presented spanning the period 1 January 201131 December 2014 to improve the understanding of the hydrological processes leading to flash floods and the relation between rainfall, runoff, erosion and sediment transport in a mesoscale catchment (Auzon, 116 km(2)) of the Mediterranean region. Badlands are present in the Auzon catchment and well connected to high-gradient channels of bedrock rivers which promotes the transfer of suspended solids downstream. The number of observed variables, the various sensors involved (both in situ and remote) and the space-time resolution (similar to km(2), similar to min) of this comprehensive dataset make it a unique contribution to research communities focused on hydrometeorology, surface hydrology and erosion. Given that rainfall is highly variable in space and time in this region, the observation system enables assessment of the hydrological response to rainfall fields. Indeed, (i) rainfall data are provided by rain gauges (both a research network of 21 rain gauges with a 5 min time step and an operational network of 10 rain gauges with a 5 min or 1 h time step), S-band Doppler dual-polarization radars (1 km(2), 5 min resolution), disdrometers (16 sensors working at 30 s or 1 min time step) and Micro Rain Radars (5 sensors, 100m height resolution). Additionally, during the special observation period (SOP-1) of the HyMeX (Hydrological Cycle in the Mediterranean Experiment) project, two X-band radars provided precipitation measurements at very fine spatial and temporal scales (1 ha, 5 min). (ii) Other meteorological data are taken from the operational surface weather observation stations of Meteo-France (including 2m air temperature, atmospheric pressure, 2 m relative humidity, 10m wind speed and direction, global radiation) at the hourly time resolution (six stations in the region of interest). (iii) The monitoring of surface hydrology and suspended sediment is multi-scale and based on nested catchments. Three hydrometric stations estimate water discharge at a 2-10 min time resolution. Two of these stations also measure additional physico-chemical variables (turbidity, temperature, conductivity) and water samples are collected automatically during floods, allowing further geochemical characterization of water and suspended solids. Two experimental plots monitor overland flow and erosion at 1 min time resolution on a hillslope with vineyard. A network of 11 sensors installed in the intermittent hydrographic network continuously measures water level and water temperature in headwater subcatchments (from 0.17 to 116 km(2)) at a time resolution of 2-5 min. A network of soil moisture sensors enables the continuous measurement of soil volumetric water content at 20 min time resolution at 9 sites. Additionally, concomitant observations (soil moisture measurements and stream gauging) were performed during floods between 2012 and 2014. Finally, this dataset is considered appropriate for understanding the rainfall variability in time and space at fine scales, improving areal rainfall estimations and progressing in distributed hydrological and erosion modelling.}, language = {en} } @misc{PrahlBoettleCostaetal.2018, author = {Prahl, Boris F. and Boettle, Markus and Costa, Lu{\´i}s F{\´i}l{\´i}pe Carvalho da and Kropp, J{\"u}rgen and Rybski, Diego}, title = {Damage and protection cost curves for coastal floods within the 600 largest European cities}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {938}, issn = {1866-8372}, doi = {10.25932/publishup-45967}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-459672}, pages = {20}, year = {2018}, abstract = {The economic assessment of the impacts of storm surges and sea-level rise in coastal cities requires high-level information on the damage and protection costs associated with varying flood heights. We provide a systematically and consistently calculated dataset of macroscale damage and protection cost curves for the 600 largest European coastal cities opening the perspective for a wide range of applications. Offering the first comprehensive dataset to include the costs of dike protection, we provide the underpinning information to run comparative assessments of costs and benefits of coastal adaptation. Aggregate cost curves for coastal flooding at the city-level are commonly regarded as by-products of impact assessments and are generally not published as a standalone dataset. Hence, our work also aims at initiating a more critical discussion on the availability and derivation of cost curves.}, language = {en} } @article{ScherlerSchwanghart2020, author = {Scherler, Dirk and Schwanghart, Wolfgang}, title = {Drainage divide networks}, series = {Earth surface dynamics}, volume = {8}, journal = {Earth surface dynamics}, number = {2}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {2196-6311}, doi = {10.5194/esurf-8-261-2020}, pages = {261 -- 274}, year = {2020}, abstract = {Drainage divides are organized into tree-like networks that may record information about drainage divide mobility. However, views diverge about how to best assess divide mobility. Here, we apply a new approach of automatically extracting and ordering drainage divide networks from digital elevation models to results from landscape evolution model experiments. We compared landscapes perturbed by strike-slip faulting and spatiotemporal variations in erodibility to a reference model to assess which topographic metrics (hillslope relief, flow distance, and chi) are diagnostic of divide mobility. Results show that divide segments that are a minimum distance of similar to 5 km from river confluences strive to attain constant values of hillslope relief and flow distance to the nearest stream. Disruptions of such patterns can be related to mobile divides that are lower than stable divides, closer to streams, and often asymmetric in shape. In general, we observe that drainage divides high up in the network, i.e., at great distances from river confluences, are more susceptible to disruptions than divides closer to these confluences and are thus more likely to record disturbance for a longer time period. We found that across-divide differences in hillslope relief proved more useful for assessing divide migration than other tested metrics. However, even stable drainage divide networks exhibit across-divide differences in any of the studied topographic metrics. Finally, we propose a new metric to quantify the connectivity of divide junctions.}, language = {en} } @misc{SchroenKoehliScheiffeleetal.2017, author = {Schr{\"o}n, Martin and K{\"o}hli, Markus and Scheiffele, Lena and Iwema, Joost and Bogena, Heye R. and Lv, Ling and Martini, Edoardo and Baroni, Gabriele and Rosolem, Rafael and Weimar, Jannis and Mai, Juliane and Cuntz, Matthias and Rebmann, Corinna and Oswald, Sascha Eric and Dietrich, Peter and Schmidt, Ulrich and Zacharias, Steffen}, title = {Improving calibration and validation of cosmic-ray neutron sensors in the light of spatial sensitivity}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, number = {636}, doi = {10.25932/publishup-41913}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-419134}, pages = {5009 -- 5030}, year = {2017}, abstract = {In the last few years the method of cosmic-ray neutron sensing (CRNS) has gained popularity among hydrologists, physicists, and land-surface modelers. The sensor provides continuous soil moisture data, averaged over several hectares and tens of decimeters in depth. However, the signal still may contain unidentified features of hydrological processes, and many calibration datasets are often required in order to find reliable relations between neutron intensity and water dynamics. Recent insights into environmental neutrons accurately described the spatial sensitivity of the sensor and thus allowed one to quantify the contribution of individual sample locations to the CRNS signal. Consequently, data points of calibration and validation datasets are suggested to be averaged using a more physically based weighting approach. In this work, a revised sensitivity function is used to calculate weighted averages of point data. The function is different from the simple exponential convention by the extraordinary sensitivity to the first few meters around the probe, and by dependencies on air pressure, air humidity, soil moisture, and vegetation. The approach is extensively tested at six distinct monitoring sites: two sites with multiple calibration datasets and four sites with continuous time series datasets. In all cases, the revised averaging method improved the performance of the CRNS products. The revised approach further helped to reveal hidden hydrological processes which otherwise remained unexplained in the data or were lost in the process of overcalibration. The presented weighting approach increases the overall accuracy of CRNS products and will have an impact on all their applications in agriculture, hydrology, and modeling.}, language = {en} } @misc{SkinnerCoulthardSchwanghartetal.2018, author = {Skinner, Christopher J. and Coulthard, Tom J. and Schwanghart, Wolfgang and Van De Wiel, Marco J. and Hancock, Greg}, title = {Global sensitivity analysis of parameter uncertainty in landscape evolution models}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1084}, issn = {1866-8372}, doi = {10.25932/publishup-46801}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-468013}, pages = {4873 -- 4888}, year = {2018}, abstract = {The evaluation and verification of landscape evolution models (LEMs) has long been limited by a lack of suitable observational data and statistical measures which can fully capture the complexity of landscape changes. This lack of data limits the use of objective function based evaluation prolific in other modelling fields, and restricts the application of sensitivity analyses in the models and the consequent assessment of model uncertainties. To overcome this deficiency, a novel model function approach has been developed, with each model function representing an aspect of model behaviour, which allows for the application of sensitivity analyses. The model function approach is used to assess the relative sensitivity of the CAESAR-Lisflood LEM to a set of model parameters by applying the Morris method sensitivity analysis for two contrasting catchments. The test revealed that the model was most sensitive to the choice of the sediment transport formula for both catchments, and that each parameter influenced model behaviours differently, with model functions relating to internal geomorphic changes responding in a different way to those relating to the sediment yields from the catchment outlet. The model functions proved useful for providing a way of evaluating the sensitivity of LEMs in the absence of data and methods for an objective function approach.}, language = {en} } @misc{ZeheBloeschl2004, author = {Zehe, Erwin and Bl{\"o}schl, G{\"u}nter}, title = {Predictability of hydrologic response at the plot and catchment scales: Role of initial conditions}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-60119}, year = {2004}, abstract = {This paper examines the effect of uncertain initial soil moisture on hydrologic response at the plot scale (1 m2) and the catchment scale (3.6 km2) in the presence of threshold transitions between matrix and preferential flow. We adopt the concepts of microstates and macrostates from statistical mechanics. The microstates are the detailed patterns of initial soil moisture that are inherently unknown, while the macrostates are specified by the statistical distributions of initial soil moisture that can be derived from the measurements typically available in field experiments. We use a physically based model and ensure that it closely represents the processes in the Weiherbach catchment, Germany. We then use the model to generate hydrologic response to hypothetical irrigation events and rainfall events for multiple realizations of initial soil moisture microstates that are all consistent with the same macrostate. As the measures of uncertainty at the plot scale we use the coefficient of variation and the scaled range of simulated vertical bromide transport distances between realizations. At the catchment scale we use similar statistics derived from simulated flood peak discharges. The simulations indicate that at both scales the predictability depends on the average initial soil moisture state and is at a minimum around the soil moisture value where the transition from matrix to macropore flow occurs. The predictability increases with rainfall intensity. The predictability increases with scale with maximum absolute errors of 90 and 32\% at the plot scale and the catchment scale, respectively. It is argued that even if we assume perfect knowledge on the processes, the level of detail with which one can measure the initial conditions along with the nonlinearity of the system will set limits to the repeatability of experiments and limits to the predictability of models at the plot and catchment scales.}, language = {de} }