Earth System Models (ESMs) are essential tools for understanding and predicting global change, but they cannot explicitly resolve hillslope-scale terrain structures that fundamentally organize water, energy, and biogeochemical stores and fluxes at subgrid scales. Here we bring together hydrologists, Critical Zone scientists, and ESM developers, to explore how hillslope structures may modulate ESM grid-level water, energy, and biogeochemical fluxes. In contrast to the one-dimensional (1-D), 2-to 3-m deep, and free-draining soil hydrology in most ESM land models, we hypothesize that 3-D, lateral ridge-to-valley flow through shallow and deep paths and insolation contrasts between sunny and shady slopes are the top two globally quantifiable organizers of water and energy (and vegetation) within an ESM grid cell. We hypothesize that these two processes are likely to impact ESM predictions where (and when) water and/or energy are limiting. We further hypothesize that, if implemented in ESM land models, these processes will increase simulated continental water storage and residence time, buffering terrestrial ecosystems against seasonal and interannual droughts. We explore efficient ways to capture these mechanisms in ESMs and identify critical knowledge gaps preventing us from scaling up hillslope to global processes. One such gap is our extremely limited knowledge of the subsurface, where water is stored (supporting vegetation) and released to stream baseflow (supporting aquatic ecosystems). We conclude with a set of organizing

Earth's terrestrial near-subsurface environment can be divided into relatively porous layers of soil, intact regolith, and sedimentary deposits above unweathered bedrock. Variations in the thicknesses of these layers control the hydrologic and biogeochemical responses of landscapes. Currently, Earth System Models approximate the thickness of these relatively permeable layers above bedrock as uniform globally, despite the fact that their thicknesses vary systematically with topography, climate, and geology. To meet the need for more realistic input data for models, we developed a high-resolution gridded global data set of the average thicknesses of soil, intact regolith, and sedimentary deposits within each 30 arcsec (1 km) pixel using the best available data for topography, climate, and geology as input. Our data set partitions the global land surface into upland hillslope, upland valley bottom, and lowland landscape components and uses models optimized for each landform type to estimate the thicknesses of each subsurface layer. On hillslopes, the data set is calibrated and validated using independent data sets of measured soil thicknesses from the U.S. and Europe and on lowlands using depth to bedrock observations from groundwater wells in the U.S. We anticipate that the data set will prove useful as an input to regional and global hydrological and ecosystems models.

Abstract. This paper presents a study on the optimal setup for discharge assimilation within a spatially distributed hydrological model. The Ensemble Kalman filter (EnKF) is employed to update the grid-based distributed states of such an hourly spatially distributed version of the HBV-96 model. By using a physically based model for the routing, the time delay and attenuation are modelled more realistically. The discharge and states at a given time step are assumed to be dependent on the previous time step only (Markov property).Synthetic and real world experiments are carried out for the Upper Ourthe (1600 km 2 ), a relatively quickly responding catchment in the Belgian Ardennes. We assess the impact on the forecasted discharge of (1) various sets of the spatially distributed discharge gauges and (2) the filtering frequency. The results show that the hydrological forecast at the catchment outlet is improved by assimilating interior gauges. This augmentation of the observation vector improves the forecast more than increasing the updating frequency. In terms of the model states, the EnKF procedure is found to mainly change the pdfs of the two routing model storages, even when the uncertainty in the discharge simulations is smaller than the defined observation uncertainty.

Abstract. The identification of hydrological drought at global scale has received considerable attention during the last decade. However, climate-induced variation in runoff across the world makes such analyses rather complicated. This especially holds for the drier regions of the world (both cold and warm), where, for a considerable period of time, zero runoff can be observed. In the current paper, we present a method that enables to identify drought at global scale across climate regimes in a consistent manner. The method combines the characteristics of the classical variable threshold level method that is best applicable in regions with nonzero runoff most of the time, and the consecutive dry days (period) method that is better suited for areas where zero runoff occurs. The newly presented method allows a drought in periods with runoff to continue in the following period without runoff. The method is demonstrated by identifying droughts from discharge observations of four rivers situated within different climate regimes, as well as from simulated runoff data at global scale obtained from an ensemble of five different land surface models. The identified drought events obtained by the new approach are compared to those resulting from application of the variable threshold level method or the consecutive dry period method separately. Results show that, in general, for drier regions, the threshold level method overestimates drought duration, because zero runoff periods are included in a drought, according to the definition used within this method. The consecutive dry period method underestimates drought occurrence, since it cannot identify droughts for periods with runoff. The developed method especially shows its relevance in transitional areas, because, in wetter regions, results are identical to the classical threshold level method. By combining both methods, the new method is able to identify single drought events that occur during positive and zero runoff periods, leading to a more realistic global drought characterization, especially within drier environments.

Hillslope‐scale rainfall‐runoff processes leading to a fast catchment response are not explicitly included in land surface models (LSMs) for use in earth system models (ESMs) due to computational constraints. This study presents a hybrid‐3D hillslope hydrological model (h3D) that couples a 1‐D vertical soil column model with a lateral pseudo‐2D saturated zone and overland flow model for use in ESMs. By representing vertical and lateral responses separately at different spatial resolutions, h3D is computationally efficient. The h3D model was first tested for three different hillslope planforms (uniform, convergent and divergent). We then compared h3D (with single and multiple soil columns) with a complex physically based 3‐D model and a simple 1‐D soil moisture model coupled with an unconfined aquifer (as typically used in LSMs). It is found that simulations obtained by the simple 1‐D model vary considerably from the complex 3‐D model and are not able to represent hillslope‐scale variations in the lateral flow response. In contrast, the single soil column h3D model shows a much better performance and saves computational time by 2‐3 orders of magnitude compared with the complex 3‐D model. When multiple vertical soil columns are implemented, the resulting hydrological responses (soil moisture, water table depth, and base flow along the hillslope) from h3D are nearly identical to those predicted by the complex 3‐D model, but still saves computational time. As such, the computational efficiency of the h3D model provides a valuable and promising approach to incorporating hillslope‐scale hydrological processes into continental and global‐scale ESMs.

Abstract. On 26 August 2010 the eastern part of The Netherlands and the bordering part of Germany were struck by a series of rainfall events lasting for more than a day. Over an area of 740 km 2 more than 120 mm of rainfall were observed in 24 h. This extreme event resulted in local flooding of city centres, highways and agricultural fields, and considerable financial loss.In this paper we report on the unprecedented flash flood triggered by this exceptionally heavy rainfall event in the 6.5 km 2 Hupsel Brook catchment, which has been the experimental watershed employed by Wageningen University since the 1960s. This study aims to improve our understanding of the dynamics of such lowland flash floods. We present a detailed hydrometeorological analysis of this extreme event, focusing on its synoptic meteorological characteristics, its space-time rainfall dynamics as observed with rain gauges, weather radar and a microwave link, as well as the measured soil moisture, groundwater and discharge response of the catchment.At the Hupsel Brook catchment 160 mm of rainfall was observed in 24 h, corresponding to an estimated return period of well over 1000 years. As a result, discharge at the catchment outlet increased from 4.4 × 10 −3 to nearly 5 m 3 s −1 . Within 7 h discharge rose from 5 × 10 −2 to 4.5 m 3 s −1 .Correspondence to: C. C. Brauer (claudia.brauer@wur.nl)The catchment response can be divided into four phases: (1) soil moisture reservoir filling, (2) groundwater response, (3) surface depression filling and surface runoff and (4) backwater feedback. The first 35 mm of rainfall were stored in the soil without a significant increase in discharge. Relatively dry initial conditions (in comparison to those for past discharge extremes) prevented an even faster and more extreme hydrological response.

[1] Radars are known for their ability to obtain a wealth of information about spatial storm field characteristics. Unfortunately, rainfall estimates obtained by this instrument are known to be affected by multiple sources of error. Especially for stratiform precipitation systems, the quality of radar rainfall estimates starts to decrease at relatively close ranges. In the current study, the hydrological potential of weather radar is analyzed during a winter halfyear for the hilly region of the Belgian Ardennes. A correction algorithm is proposed which corrects the radar data for errors related to attenuation, ground clutter, anomalous propagation, the vertical profile of reflectivity (VPR), and advection. No final bias correction with respect to rain gauge data was implemented because such an adjustment would not add to a better understanding of the quality of the radar data. The impact of the different corrections is assessed using rainfall information sampled by 42 hourly rain gauges. The largest improvement in the quality of the radar data is obtained by correcting for ground clutter. The impact of VPR correction and advection depends on the spatial variability and velocity of the precipitation system. Overall during the winter period, the radar underestimates the amount of precipitation as compared to the rain gauges. Remaining differences between both instruments can be attributed to spatial and temporal variability in the type of precipitation, which has not been taken into account.Citation: Hazenberg, P., H. Leijnse, and R. Uijlenhoet (2011), Radar rainfall estimation of stratiform winter precipitation in the

One of the recognized weaknesses of land surface models as used in weather and climate models is the assumption of constant soil thickness because of the lack of global estimates of bedrock depth. Using a 30-arc-s global dataset for the thickness of relatively porous, unconsolidated sediments over bedrock, spatial variation in soil thickness is included here in version 4.5 of the Community Land Model (CLM4.5). The number of soil layers for each grid cell is determined from the average soil depth for each 0.98 latitude 3 1.258 longitude grid cell. The greatest changes in the simulation with variable soil thickness are to baseflow, with the annual minimum generally occurring earlier. Smaller changes are seen in latent heat flux and surface runoff primarily as a result of an increase in the annual cycle amplitude. These changes are related to soil moisture changes that are most substantial in locations with shallow bedrock. Total water storage (TWS) anomalies are not strongly affected over most river basins since most basins contain mostly deep soils, but TWS anomalies are substantially different for a river basin with more mountainous terrain. Additionally, the annual cycle in soil temperature is partially affected by including realistic soil thicknesses resulting from changes in the vertical profile of heat capacity and thermal conductivity. However, the largest changes to soil temperature are introduced by the soil moisture changes in the variable soil thickness simulation. This implementation of variable soil thickness represents a step forward in land surface model development.