This paper is the outcome of a community initiative to identify major unsolved scientific problems in hydrology motivated by a need for stronger harmonisation of research efforts. The procedure involved a public consultation through online media, followed by two workshops through which a large number of potential science questions were collated, prioritised, and synthesised. In spite of the diversity of the participants (230 scientists in total), the process revealed much about community priorities and the state of our science: a preference for continuity in research questions rather than radical departures or redirections from past and current work. Questions remain focused on the process-based understanding of hydrological variability and causality at all space and time scales. Increased attention to environmental change drives a new emphasis on understanding how change propagates across interfaces within the hydrological system and across disciplinary boundaries. In particular, the expansion of the human footprint raises a new set of questions related to human interactions with nature and water cycle feedbacks in the context of complex water management problems. We hope that this reflection and synthesis of the 23 unsolved problems in hydrology will help guide research efforts for some years to come. ARTICLE HISTORY

We analyze the relative importance of the selection of (1) the geostatistical model depicting the structural heterogeneity of an aquifer, and (2) the basic processes to be included in the conceptual model, to describe the main aspects of solute transport at an experimental site. We focus on the results of a forced-gradient tracer test performed at the "Lauswiesen" experimental site, near Tübingen, Germany. In the experiment, NaBr is injected into a well located 52 m from a pumping well. Multilevel breakthrough curves (BTCs) are measured in the latter. We conceptualize the aquifer as a three-dimensional, doubly stochastic composite medium, where distributions of geomaterials and attributes, e.g., hydraulic conductivity (K) and porosity (phi), can be uncertain. Several alternative transport processes are considered: advection, advection-dispersion and/or mass-transfer between mobile and immobile regions. Flow and transport are tackled within a stochastic Monte Carlo framework to describe key features of the experimental BTCs, such as temporal moments, peak time, and pronounced tailing. We find that, regardless the complexity of the conceptual transport model adopted, an adequate description of heterogeneity is crucial for generating alternative equally likely realizations of the system that are consistent with (a) the statistical description of the heterogeneous system, as inferred from the data, and (b) salient features of the depth-averaged breakthrough curve, including preferential paths, slow release of mass particles, and anomalous spreading. While the available geostatistical characterization of heterogeneity can explain most of the integrated behavior of transport (depth-averaged breakthrough curve), not all multilevel BTCs are described with equal success. This suggests that transport models simply based on integrated measurements may not ensure an accurate representation of many of the important features required in three-dimensional transport models.

No abstract

No abstract

Many hydrological (as well as diverse earth, environmental, ecological, biological, physical, social, financial and other) variables, Y, exhibit frequency distributions that are difficult to reconcile with those of their spatial or temporal increments, DY. Whereas distributions of Y (or its logarithm) are at times slightly asymmetric with relatively mild peaks and tails, those of DY tend to be symmetric with peaks that grow sharper, and tails that become heavier, as the separation distance (lag) between pairs of Y values decreases. No statistical model known to us captures these behaviors of Y and DY in a unified and consistent manner. We propose a new, generalized sub-Gaussian model that does so. We derive analytical expressions for probability distribution functions (pdfs) of Y and DY as well as corresponding lead statistical moments. In our model the peak and tails of the DY pdf scale with lag in line with observed behavior. The model allows one to estimate, accurately and efficiently, all relevant parameters by analyzing jointly sample moments of Y and DY. We illustrate key features of our new model and method of inference on synthetically generated samples and neutron porosity data from a deep borehole.

No abstract

[1] The analysis of pumping tests has traditionally relied on analytical solutions of groundwater flow equations in relatively simple domains, consisting of one or at most a few units having uniform hydraulic properties. Recently, attention has been shifting toward methods and solutions that would allow one to characterize subsurface heterogeneities in greater detail. On one hand, geostatistical inverse methods are being used to assess the spatial variability of parameters, such as permeability and porosity, on the basis of multiple cross-hole pressure interference tests. On the other hand, analytical solutions are being developed to describe the mean and variance (first and second statistical moments) of flow to a well in a randomly heterogeneous medium. We explore numerically the feasibility of using a simple graphical approach (without numerical inversion) to estimate the geometric mean, integral scale, and variance of local log transmissivity on the basis of quasi steady state head data when a randomly heterogeneous confined aquifer is pumped at a constant rate. By local log transmissivity we mean a function varying randomly over horizontal distances that are small in comparison with a characteristic spacing between pumping and observation wells during a test. Experimental evidence and hydrogeologic scaling theory suggest that such a function would tend to exhibit an integral scale well below the maximum well spacing. This is in contrast to equivalent transmissivities derived from pumping tests by treating the aquifer as being locally uniform (on the scale of each test), which tend to exhibit regional-scale spatial correlations. We show that whereas the mean and integral scale of local log transmissivity can be estimated reasonably well based on theoretical ensemble mean variations of head and drawdown with radial distance from a pumping well, estimating the log transmissivity variance is more difficult. We obtain reasonable estimates of the latter based on theoretical variation of the standard deviation of circumferentially averaged drawdown about its mean.

No abstract