Falsification and corroboration of conceptual hydrological models using geophysical data

Geophysical data may provide crucial information about hydrological properties, states, and processes that are difficult to obtain by other means. Large data sets can be acquired over widely different scales in a minimally invasive manner and at comparatively low costs, but their effective use in hydrology makes it necessary to understand the fidelity of geophysical models, the assumptions made in their construction, and the links between geophysical and hydrological properties. Geophysics has been applied for groundwater prospecting for almost a century, but it is only in the last 20 years that it is regularly used together with classical hydrological data to build predictive hydrological models. A largely unexplored venue for future work is to use geophysical data to falsify or rank competing conceptual hydrological models. A promising cornerstone for such a model selection strategy is the Bayes factor, but it can only be calculated reliably when considering the main sources of uncertainty throughout the hydrogeophysical parameter estimation process. Most classical geophysical imaging tools tend to favor models with smoothly varying property fields that are at odds with most conceptual hydrological models of interest. It is thus necessary to account for this bias or use alternative approaches in which proposed conceptual models are honored at all steps in the model building process.

Morphological, hydrological, biogeochemical and ecological changes and challenges in river restoration - the Thur River case study

River restoration can enhance river dynamics, environmental heterogeneity and biodiversity, but the underlying processes governing the dynamic changes need to be understood to ensure that restoration projects meet their goals, and adverse effects are prevented. In particular, we need to comprehend how hydromorphological variability quantitatively relates to ecosystem functioning and services, biodiversity as well as ground-and surface water quality in restored river corridors. This involves (i) physical processes and structural properties, determining erosion and sedimentation, as well as solute and heat transport behavior in surface water and within the subsurface; (ii) biogeochemical processes and characteristics, including the turnover of nutrients and natural water constituents; and (iii) ecological processes and indicators related to biodiversity and ecological functioning. All these aspects are interlinked, requiring an interdisciplinary investigation approach. Here, we present an overview of the recently completed RECORD (REstored CORridor Dynamics) project in which we combined physical, chemical, and biological observations with modeling at a restored river corridor of the perialpine Thur River in Switzerland. Our results show that river restoration, beyond inducing morphologic changes that reshape the river bed and banks, triggered complex spatial patterns of bank infiltration, and affected habitat type, biotic communities and biogeochemical processes. We adopted an interdisciplinary approach of monitoring the continuing changes due to restoration measures to address the following questions: How stable is the morphological variability established by restoration? Does morphological variability guarantee an improvement in biodiversity? How does morphological variability affect biogeochemical transformations in the river corridor? What are some potential adverse effects of river restoration? How is river restoration influenced by catchment-scale hydraulics [GRAPHICS] and which feedbacks exist on the large scale? Beyond summarizing the major results of individual studies within the project, we show that these overarching questions could only be addressed in an interdisciplinary framework.

Runoff and sediment production in a Mediterranean basin under two different land uses

This study analyses the infl uence of two diff erent land uses on the hydrology of the Vernegà experimental basin between the years 2005 and 2009. It is located in the Northeast of the Iberian Peninsula and is infl uenced by a Mediterranean climate, with an average annual rainfall of 646 mm. The study of rainfall distributi on in the 1982 to 2009 period shows that the majority occurs during autumn, winter, and spring representi ng 33.2%, 25.5%, and 25.7% respecti vely of the total annual rainfall. Surface runoff is concentrated from October to June. Between 2005 and 2009, total runoff was 242.38 mm at the"Bosc" gauging stati on, which drains an area of 1.60 km2, and 298.54 mm at the"Campàs" gauging stati on, which drains an area of 2.57 km2 and is located at the outlet of the basin. More than 80% of the total surface runoff yielded during the study period corresponds to the 2005-2006 hydrologic year. Finally, Campàs gauging stati on registers a higher total runoff than Bosc gauging stati on. Part of this phenomenon may be due to the intercepti on of rainfall and plant biomass in the forested area of the basin. In relati on to the sediment yield, an overall increase in the two basins has been detected. Recent forest management practi ces undertaken in the catchment area are considered to be one of the most important reasons for this change.

Characteristics of 2D convective structures in Catalonia (NE Spain): an analysis using radar data and GIS

Flood simulation studies use spatial-temporal rainfall data input into distributed hydrological models. A correct description of rainfall in space and in time contributes to improvements on hydrological modelling and design. This work is focused on the analysis of 2-D convective structures (rain cells), whose contribution is especially significant in most flood events. The objective of this paper is to provide statistical descriptors and distribution functions for convective structure characteristics of precipitation systems producing floods in Catalonia (NE Spain). To achieve this purpose heavy rainfall events recorded between 1996 and 2000 have been analysed. By means of weather radar, and applying 2-D radar algorithms a distinction between convective and stratiform precipitation is made. These data are introduced and analyzed with a GIS. In a first step different groups of connected pixels with convective precipitation are identified. Only convective structures with an area greater than 32 km2 are selected. Then, geometric characteristics (area, perimeter, orientation and dimensions of the ellipse), and rainfall statistics (maximum, mean, minimum, range, standard deviation, and sum) of these structures are obtained and stored in a database. Finally, descriptive statistics for selected characteristics are calculated and statistical distributions are fitted to the observed frequency distributions. Statistical analyses reveal that the Generalized Pareto distribution for the area and the Generalized Extreme Value distribution for the perimeter, dimensions, orientation and mean areal precipitation are the statistical distributions that best fit the observed ones of these parameters. The statistical descriptors and the probability distribution functions obtained are of direct use as an input in spatial rainfall generators.

Biogeochemical processes controlling oxygen and carbon isotopes of diatom silica in Late Glacial to Holocene lacustrine rhythmites

Biogeochemical cycles and sedimentary records in lakes are related to climate controls on hydrology and catchment processes. Changes in the isotopic imposition of the diatom frustules (δ 18 O diatom and δ 13 C diatom ) in lacustrine sediments can be used to reconstruct palaeoclimatic and palaeoenvironmental changes. The Lago Chungará (Andean Altiplano, 18°15 ′ S, 69°10 ′ W, 4520 masl) diatomaceous laminated sediments are made up of white and green multiannual rhythmites. White laminae were formed during short-term diatom super-blooms, and are composed almost exclusively of large-sized Cyclostephanos andinus.These diatoms bloom during mixing events when recycled nutrients from the bottom waters are brought to the surface and/or when nutrients are introduced from the catchment during periods of strong runoff. Conversely, the green laminae are thought to have been deposited over several years and are composed of a mixture of diatoms (mainly smaller valves of C. andinus and Discostella stelligera ) and organic matter. These green laminae reflect the lake's hydrological recovery from a status favouring the diatom super-blooms (white laminae) towards baseline conditions. δ 18 O diatom and δ 13 C diatom from 11,990 to 11,530 cal years BP allow us to reconstruct shifts in the precipitation/evaporation ratio and changes in the lake water dissolved carbon concentration, respectively. δ 18 O diatom values indicate that white laminae formation occurred mainly during low lake level stages, whereas green laminae formation generally occurred during high lake level stages. The isotope and chronostratigraphical data together suggest that white laminae deposition is caused by extraordinary environmental events. El Niño-Southern Oscillation and changes in solar activity are the most likely climate forcing mechanisms that could trigger such events, favouring hydrological changes at interannual-to-decadal scale. This study demonstrates the potential for laminated lake sediments to document extreme pluriannual events.

Propuesta de un protocolo de evaluación de la calidad ecológica de ríos andinos (CERA) y su aplicación a dos cuencas de Ecuador y Perú

ABSTRACT Proposal for an evaluation protocol of the ecological quality of Andean rivers (CERA) and its use in two basins in Ecuador and Peru A Rapid Protocol is presented for Evaluation of the Ecological Status of Andean Rivers (CERA) localized over 2000 m.a.s.l. from the Northern Andes (Venezuela) through the Altiplano in the Central Andes (Bolivia). This protocol was used in 45 sampling sites in the Guayllabamba River Basin in Ecuador and in 42 sampling sites in the Ca nete River Basin in Peru. Previously, and in order to test if the sampling stations may or not be considered reference stations, we constructed a method that assesses 24 basin attributes, hydrology, reach and riverbed and that uctuates from 24 to 120 points; sites with values higher than 100 were considered as potential reference sites. Besides the benthic macroinvertebrats" evaluation, the river habitat and riparian vegetation were also evaluated through of the application of the indices ABI (R´ os et al., submitted), IHF (Pardo et al., 2002) and QBR-And, respectively. The convenience of the initial allocation of the reference sites was evaluated as well. These indices have been properly adapted to the conditions and characteristics of the high Andes rivers. The results obtained for both basins were compared and discussed. Through the use of the CERA protocol, the particular perturbation gradients and the natural variability of the reference sites in both countries were recognized. RESUMEN Propuesta de un protocolo de evaluación de la calidad ecológica de ríos andinos (CERA) y su aplicaci´on a dos cuencas en Ecuador y Perú Se presenta un protocolo rápido de evaluación de la Calidad Ecológica de Ríos Andinos (CERA), situados sobre los 2000 m.s.n.m, desde los Andes del Norte (Venezuela) hasta el Altiplano de los Andes Centrales (Bolivia). Este protocolo ha sido aplicado en 45 estaciones de muestreo en la cuenca del río Guayllabamba en Ecuador y en 42 estaciones de muestreo en la cuenca del río Cañete en Perú. Previamente, para probar si las estaciones de muestreo pueden o no ser estaciones de referencia construimos un método que valora 24 atributos de cuenca, hidrología, tramo y lecho y que fluctúa de 24 a 120 puntos; valores superiores a 100 fueron considerados como sitios potencialmente de referencia. Además del estudio de los macroinvertebrados bentónicos, se evaluó el hábitat fluvial y la comunidad vegetal de ribera a través de la aplicación de los índices ABI (Ríos et al., sometido), IHF (Pardo et al., 2002) y QBR-And respectivamente; así como la conveniencia de la asignación inicial de las estaciones de referencia. Estos índices han sido adecuadamente adaptados a las condiciones y características propias de los ríos altoandinos. Los resultados obtenidos fueron comparados y discutidos entre ambas cuencas. Mediante la aplicación del protocolo CERA se han reconocido los respectivos gradientes de perturbación y la variabilidad natural de las estaciones de referencia en ambos países.

Effect of radar rainfall time resolution on the predictive capability of a distributed hydrologic model

The performance of a hydrologic model depends on the rainfall input data, both spatially and temporally. As the spatial distribution of rainfall exerts a great influence on both runoff volumes and peak flows, the use of a distributed hydrologic model can improve the results in the case of convective rainfall in a basin where the storm area is smaller than the basin area. The aim of this study was to perform a sensitivity analysis of the rainfall time resolution on the results of a distributed hydrologic model in a flash-flood prone basin. Within such a catchment, floods are produced by heavy rainfall events with a large convective component. A second objective of the current paper is the proposal of a methodology that improves the radar rainfall estimation at a higher spatial and temporal resolution. Composite radar data from a network of three C-band radars with 6-min temporal and 2 × 2 km2 spatial resolution were used to feed the RIBS distributed hydrological model. A modification of the Window Probability Matching Method (gauge-adjustment method) was applied to four cases of heavy rainfall to improve the observed rainfall sub-estimation by computing new Z/R relationships for both convective and stratiform reflectivities. An advection correction technique based on the cross-correlation between two consecutive images was introduced to obtain several time resolutions from 1 min to 30 min. The RIBS hydrologic model was calibrated using a probabilistic approach based on a multiobjective methodology for each time resolution. A sensitivity analysis of rainfall time resolution was conducted to find the resolution that best represents the hydrological basin behaviour.

A look to the HyMeX program (Una mirada al programa HyMeX)

The international HyMeX (Hydrological Mediterranean Experiment) program aims to improve our understanding of the water cycle in the Mediterranean, using a multidisciplinary and multiscale approach and with emphasis on extreme events. This program will improve our understanding and our predictive ability of hydrometeorological hazards including their evolution within the next century. One of the most important results of the program will be its observational campaigns, which will greatly improve the data available, leading to significant scientific results. The interest of the program for the Spanish research groups is described, as the active participation of some of them in the design and execution of the observational activities. At the same time, due to its location, Spain is key to the program, being a good observation platform. HyMeX will enrich the work of the Spanish research groups, it will improve the predictive ability of the weather services, will help us to have a better understanding of the impacts of hydrometeorological extremes on our society and will lead to better strategies for adapting to climate change.

Lidar measurement of surface melt for a temperate Alpine glacier at the seasonal and hourly scales

This study shows how a new generation of terrestrial laser scanners can be used to investigate glacier surface ablation and other elements of glacial hydrodynamics at exceptionally high spatial and temporal resolution. The study area is an Alpine valley glacier, Haut Glacier d'Arolla, Switzerland. Here we use an ultra-long-range lidar RIEGL VZ-6000 scanner, having a laser specifically designed for measurement of snow- and ice-cover surfaces. We focus on two timescales: seasonal and daily. Our results show that a near-infrared scanning laser system can provide high-precision elevation change and ablation data from long ranges, and over relatively large sections of the glacier surface. We use it to quantify spatial variations in the patterns of surface melt at the seasonal scale, as controlled by both aspect and differential debris cover. At the daily scale, we quantify the effects of ogive-related differences in ice surface debris content on spatial patterns of ablation. Daily scale measurements point to possible hydraulic jacking of the glacier associated with short-term water pressure rises. This latter demonstration shows that this type of lidar may be used to address subglacial hydrologic questions, in addition to motion and ablation measurements.