OVERLAND FLOW AND SURFACE RUNOFF / Hydrological Science and Engineering

Cape Verde, off the coast of Senegal in western Africa, is a volcanic archipelago where soil and water conservation techniques play an important role in the overall subsistence of half a million inhabitants. In fact, the step slopes in the more agricultural islands due to it's volcanic origin, together with semi-arid and arid environments (the country is located in the Sahelian region), characterized by a very irregular wet season, with high intensity rainfall events, make life tough. The hard conditions lead during the first half of the XX century to frequent cycles of drought with severe implications on the local populations, with impressive numbers of deaths by famine, and a decrease of the number of local inhabitants by more than halve in some islands. Maintain the soil in place and the water inside the soil was there after a mater of survival, and the CapeVerdians implemented over the last half century a number of soil and water conservation techniques that cover all the landscape. In this work, we monitored a number of slope soil and water conservation techniques, such as terraces, half moons, live barriers, etc, together with two cultural strategies, used to plant corn and beans on one side and peanuts on the other, with a semi-quantitative methodology, to evaluate their effectiveness. A discussion is given on the costs and effectiveness of the techniques to reduce overland flow production and therefore erosion, and to promote rainfall infiltration.

Fracture network characterization using hydrological and geophysical data

Les décisions de gestion des eaux souterraines doivent souvent être justiffées par des modèles quantitatifs d'aquifères qui tiennent compte de l'hétérogénéité des propriétés hydrauliques. Les aquifères fracturés sont parmi les plus hétérogènes et très difficiles à étudier. Dans ceux-ci, les fractures connectées, d'ouverture millimètrique, peuvent agir comme conducteurs hydrauliques et donc créer des écoulements très localisés. Le manque général d'informations sur la distribution spatiale des fractures limite la possibilité de construire des modèles quantitatifs de flux et de transport. Les données qui conditionnent les modèles sont généralement spatialement limitées, bruitées et elles ne représentent que des mesures indirectes de propriétés physiques. Ces limitations aux données peuvent être en partie surmontées en combinant différents types de données, telles que les données hydrologiques et de radar à pénétration de sol plus commun ément appelé géoradar. L'utilisation du géoradar en forage est un outil prometteur pour identiffer les fractures individuelles jusqu'à quelques dizaines de mètres dans la formation. Dans cette thèse, je développe des approches pour combiner le géoradar avec les données hydrologiques affn d'améliorer la caractérisation des aquifères fracturés. Des investigations hydrologiques intensives ont déjà été réalisées à partir de trois forage adjacents dans un aquifère cristallin en Bretagne (France). Néanmoins, la dimension des fractures et la géométrie 3-D des fractures conductives restaient mal connue. Affn d'améliorer la caractérisation du réseau de fractures je propose dans un premier temps un traitement géoradar avancé qui permet l'imagerie des fractures individuellement. Les résultats montrent que les fractures perméables précédemment identiffées dans les forages peuvent être caractérisées géométriquement loin du forage et que les fractures qui ne croisent pas les forages peuvent aussi être identiffées. Les résultats d'une deuxième étude montrent que les données géoradar peuvent suivre le transport d'un traceur salin. Ainsi, les fractures qui font partie du réseau conductif et connecté qui dominent l'écoulement et le transport local sont identiffées. C'est la première fois que le transport d'un traceur salin a pu être imagé sur une dizaines de mètres dans des fractures individuelles. Une troisième étude conffrme ces résultats par des expériences répétées et des essais de traçage supplémentaires dans différentes parties du réseau local. En outre, la combinaison des données de surveillance hydrologique et géoradar fournit la preuve que les variations temporelles d'amplitude des signaux géoradar peuvent nous informer sur les changements relatifs de concentrations de traceurs dans la formation. Par conséquent, les données géoradar et hydrologiques sont complémentaires. Je propose ensuite une approche d'inversion stochastique pour générer des modèles 3-D de fractures discrètes qui sont conditionnés à toutes les données disponibles en respectant leurs incertitudes. La génération stochastique des modèles conditionnés par géoradar est capable de reproduire les connexions hydrauliques observées et leur contribution aux écoulements. L'ensemble des modèles conditionnés fournit des estimations quantitatives des dimensions et de l'organisation spatiale des fractures hydrauliquement importantes. Cette thèse montre clairement que l'imagerie géoradar est un outil utile pour caractériser les fractures. La combinaison de mesures géoradar avec des données hydrologiques permet de conditionner avec succès le réseau de fractures et de fournir des modèles quantitatifs. Les approches présentées peuvent être appliquées dans d'autres types de formations rocheuses fracturées où la roche est électriquement résistive.

Integration of high-resolution geophysical data for the characterization of alluvial aquifers : new approaches and implications for hydrological predictions

Abstract Accurate characterization of the spatial distribution of hydrological properties in heterogeneous aquifers at a range of scales is a key prerequisite for reliable modeling of subsurface contaminant transport, and is essential for designing effective and cost-efficient groundwater management and remediation strategies. To this end, high-resolution geophysical methods have shown significant potential to bridge a critical gap in subsurface resolution and coverage between traditional hydrological measurement techniques such as borehole log/core analyses and tracer or pumping tests. An important and still largely unresolved issue, however, is how to best quantitatively integrate geophysical data into a characterization study in order to estimate the spatial distribution of one or more pertinent hydrological parameters, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first develop a strategy for the assimilation of several types of hydrogeophysical data having varying degrees of resolution, subsurface coverage, and sensitivity to the hydrologic parameter of interest. In this regard a novel simulated annealing (SA)-based conditional simulation approach was developed and then tested in its ability to generate realizations of porosity given crosshole ground-penetrating radar (GPR) and neutron porosity log data. This was done successfully for both synthetic and field data sets. A subsequent issue that needed to be addressed involved assessing the potential benefits and implications of the resulting porosity realizations in terms of groundwater flow and contaminant transport. This was investigated synthetically assuming first that the relationship between porosity and hydraulic conductivity was well-defined. Then, the relationship was itself investigated in the context of a calibration procedure using hypothetical tracer test data. Essentially, the relationship best predicting the observed tracer test measurements was determined given the geophysically derived porosity structure. Both of these investigations showed that the SA-based approach, in general, allows much more reliable hydrological predictions than other more elementary techniques considered. Further, the developed calibration procedure was seen to be very effective, even at the scale of tomographic resolution, for predictions of transport. This also held true at locations within the aquifer where only geophysical data were available. This is significant because the acquisition of hydrological tracer test measurements is clearly more complicated and expensive than the acquisition of geophysical measurements. Although the above methodologies were tested using porosity logs and GPR data, the findings are expected to remain valid for a large number of pertinent combinations of geophysical and borehole log data of comparable resolution and sensitivity to the hydrological target parameter. Moreover, the obtained results allow us to have confidence for future developments in integration methodologies for geophysical and hydrological data to improve the 3-D estimation of hydrological properties.

Determination of extreme floods using a distributed hydrological model

Hydrological models developed for extreme precipitation of PMP type are difficult to calibrate because of the scarcity of available data for these events. This article presents the process and results of calibration for a distributed hydrological model at fine scale developed for the estimation of probable maximal floods in the case of a PMP. This calibration is done on two Swiss catchments for two events of summer storms. The calculation done is concentrated on the estimation of the parameters of the model, divided in two parts. The first is necessary for the computation of flow speeds while the second is required for the determination of the initial and final infiltration capacities for each terrain type. The results, validated with the Nash equation show a good correlation between the simulated and observed flows. We also apply this model on two Romanian catchments, showing the river network and estimated flow.

Calibration of high-resolution geophysical data with tracer test measurements to improve hydrological predictions

Relationships between porosity and hydraulic conductivity tend to be strongly scale- and site-dependent and are thus very difficult to establish. As a result, hydraulic conductivity distributions inferred from geophysically derived porosity models must be calibrated using some measurement of aquifer response. This type of calibration is potentially very valuable as it may allow for transport predictions within the considered hydrological unit at locations where only geophysical measurements are available, thus reducing the number of well tests required and thereby the costs of management and remediation. Here, we explore this concept through a series of numerical experiments. Considering the case of porosity characterization in saturated heterogeneous aquifers using crosshole ground-penetrating radar and borehole porosity log data, we use tracer test measurements to calibrate a relationship between porosity and hydraulic conductivity that allows the best prediction of the observed hydrological behavior. To examine the validity and effectiveness of the obtained relationship, we examine its performance at alternate locations not used in the calibration procedure. Our results indicate that this methodology allows us to obtain remarkably reliable hydrological predictions throughout the considered hydrological unit based on the geophysical data only. This was also found to be the case when significant uncertainty was considered in the underlying relationship between porosity and hydraulic conductivity.

Kleptoparasitism, disturbance and predation of yellow-legged gulls on Audouin's gulls in three colonies of the western Mediterranean

The impact of yellow-legged gulls on Audouin´s gulls was studied by means of observations from a blind, at the Ebro Delta, the Chafarinas Islands and the Columbretes Islands colonies, during different stages of Audouin´s Gulls breeding cycle. The rates of predation (upon eggs and chicks) and kleptoparasitism (aerial and on courtship and chick-feeding regurgitates) were recorded to evaluate this impact. Kleptoparasitism and predation rates recorded at the three colonies may be considered low when compared with similar studies dealing with other gull species. It seems that interactions did not represent a threat to the population dynamics of Audouin´s Gulls in any of the three colonies. The highest rates of aerial kleptoparasitism, courtship feeding and chick mortality were recorded at the Columbretes Islands, and the lowest at the Ebro Delta. The ratio of the number of yellow-legged gulls to the number of Audouin´s gulls seems to influence the rate of disturbances. Food availability was also likely to play an important role in the number of interactions, which increased when food was in shorter supply. The study at the Ebro Delta and Chafarinas Islands was restricted to some subcolonies and hence results may not be necessarily representative of the entire colonies. The existence of some yellow-legged gull individuals specialised as predators was recorded. Finally, we discuss the suitability of several conservation measures commonly applied in colonies where both species breed syntopically.

Numerical analysis of the impact of charcoal production on soil hydrological behavior, runoff response and erosion susceptibility

The impact of charcoal production on soil hydraulic properties, runoff response and erosion susceptibility were studied in both field and simulation experiments. Core and composite samples, from 12 randomly selected sites within the catchment of Kotokosu were taken from the 0-10 cm layer of a charcoal site soil (CSS) and adjacent field soils (AFS). These samples were used to determine saturated hydraulic conductivity (Ksat), bulk density, total porosity, soil texture and color. Infiltration, surface albedo and soil surface temperature were also measured in both CSS and AFS. Measured properties were used as entries in a rainfall runoff simulation experiment on a smooth (5 % slope) plot of 25 x 25 m grids with 10 cm resolutions. Typical rainfall intensities of the study watershed (high, moderate and low) were applied to five different combinations of Ks distributions that could be expected in this landscape. The results showed significantly (p < 0.01) higher flow characteristics of the soil under charcoal kilns (increase of 88 %). Infiltration was enhanced and runoff volume reduced significantly. The results showed runoff reduction of about 37 and 18 %, and runoff coefficient ranging from 0.47-0.75 and 0.04-0.39 or simulation based on high (200 mm h-1) and moderate (100 mm h-1) rainfall events over the CSS and AFS areas, respectively. Other potential impacts of charcoal production on watershed hydrology were described. The results presented, together with watershed measurements, when available, are expected to enhance understanding of the hydrological responses of ecosystems to indiscriminate charcoal production and related activities in this region.

Hydrological processes obtained on the plot scale under four simulated rainfall tests during the cycle of different crop systems

The cropping system influences the interception of water by plants, water storage in depressions on the soil surface, water infiltration into the soil and runoff. The aim of this study was to quantify some hydrological processes under no tillage cropping systems at the edge of a slope, in 2009 and 2010, in a Humic Dystrudept soil, with the following treatments: corn, soybeans, and common beans alone; and intercropped corn and common bean. Treatments consisted of four simulated rainfall tests at different times, with a planned intensity of 64 mm h-1 and 90 min duration. The first test was applied 18 days after sowing, and the others at 39, 75 and 120 days after the first test. Different times of the simulated rainfall and stages of the crop cycle affected soil water content prior to the rain, and the time runoff began and its peak flow and, thus, the surface hydrological processes. The depth of the runoff and the depth of the water intercepted by the crop + soil infiltration + soil surface storage were affected by the crop systems and the rainfall applied at different times. The corn crop was the most effective treatment for controlling runoff, with a water loss ratio of 0.38, equivalent to 75 % of the water loss ratio exhibited by common bean (0.51), the least effective treatment in relation to the others. Total water loss by runoff decreased linearly with an increase in the time that runoff began, regardless of the treatment; however, soil water content on the gravimetric basis increased linearly from the beginning to the end of the rainfall.

Anthropogenic disturbance as a driver of microspatial and microhabitat segregation of cytotypes of Centaurea stoebe and cytotype interactions in secondary contact zones.

BACKGROUND AND AIMS: In a mixed-ploidy population, strong frequency-dependent mating will lead to the elimination of the less common cytotype, unless prezygotic barriers enhance assortative mating. However, such barriers favouring cytotype coexistence have only rarely been explored. Here, an assessment is made of the mechanisms involved in formation of mixed-ploidy populations and coexistence of diploid plants and their closely related allotetraploid derivates from the Centaurea stoebe complex (Asteraceae). METHODS: An investigation was made of microspatial and microhabitat distribution, life-history and fitness traits, flowering phenology, genetic relatedness of cytotypes and intercytotype gene flow (cpDNA and microsatellites) in six mixed-ploidy populations in Central Europe. KEY RESULTS: Diploids and tetraploids were genetically differentiated, thus corroborating the secondary origin of contact zones. The cytotypes were spatially segregated at all sites studied, with tetraploids colonizing preferentially drier and open microhabitats created by human-induced disturbances. Conversely, they were rare in more natural microsites and microsites with denser vegetation despite their superior persistence ability (polycarpic life cycle). The seed set of tetraploid plants was strongly influenced by their frequency in mixed-ploidy populations. Triploid hybrids originated from bidirectional hybridizations were extremely rare and almost completely sterile, indicating a strong postzygotic barrier between cytotypes. CONCLUSIONS: The findings suggest that tetraploids are later immigrants into already established diploid populations and that anthropogenic activities creating open niches favouring propagule introductions were the major factor shaping the non-random distribution and habitat segregation of cytotypes at fine spatial scale. Establishment and spread of tetraploids was further facilitated by their superior persistence through the perennial life cycle. The results highlight the importance of non-adaptive spatio-temporal processes in explaining microhabitat and microspatial segregation of cytotypes.

Using spatially distributed parameters and multi-response objective functions to solve parameterization of complex applications of semi-distributed hydrological models

Application of semi-distributed hydrological models to large, heterogeneous watersheds deals with several problems. On one hand, the spatial and temporal variability in catchment features should be adequately represented in the model parameterization, while maintaining the model complexity in an acceptable level to take advantage of state-of-the-art calibration techniques. On the other hand, model complexity enhances uncertainty in adjusted model parameter values, therefore increasing uncertainty in the water routing across the watershed. This is critical for water quality applications, where not only streamflow, but also a reliable estimation of the surface versus subsurface contributions to the runoff is needed. In this study, we show how a regularized inversion procedure combined with a multiobjective function calibration strategy successfully solves the parameterization of a complex application of a water quality-oriented hydrological model. The final value of several optimized parameters showed significant and consistentdifferences across geological and landscape features. Although the number of optimized parameters was significantly increased by the spatial and temporal discretization of adjustable parameters, the uncertainty in water routing results remained at reasonable values. In addition, a stepwise numerical analysis showed that the effects on calibration performance due to inclusion of different data types in the objective function could be inextricably linked. Thus caution should be taken when adding or removing data from an aggregated objective function.

Disturbance, age and size structure in stands of Pinus uncinata Ram

Las estructuras de tamaño y edad de dos bosques monoespecificos subalpinos de Pinus uncinata en el Pirineo, fueron estudiadas e interpretadas introduciendo una clasificación sencilla de los árboles según la exposición de su copa a la radiación solar directa. Ambas poblaciones muestran distribuciones de diámetro de J invertida irregulares, cuya Interpretación es incierta. Por el contrario, la estructura de edad en ambos bosques muestra picos visibles, los cuales, si se tienen en cuenta la exposición de la copa, indican que la regeneración se produce en claros del bosque. Se analizaron las series de anillos de crecimiento de los árboles de cada sitio para reconstruir su historia reciente de perturbaciones Identificando los incrementos bruscos del crecimiento. El régimen de perturbaciones en los dos bosques es diferente, y en ambos casos muestra un efecto madelador análogo sobre la estructura de edad. Se concluye que P. uncinata se comporta como una especie heliófila que puede formar bosques maduros o multicohorte bajo determinados regímenes de perturbaciones suaves, gracias a su capacidad para iniciar la regeneración en claros pequeños. Este proceso tiene lugar en un lapso de tiempo corto después de la incidencia de las perturbaciones.

Benthic macroinvertebrates assemblages and macrohabitat connectivity in Mediterranean-climate streams of northern California

Abstract. Drought leads to a loss of longitudinal and lateral hydrologic connectivity, which causes direct or indirect changes in stream ecosystem properties. Changes in macrohabitat availability from a rifflepool sequence to isolated pools are among the most conspicuous consequences of connectivity loss. Macroinvertebrate assemblages were compared among 3 distinct stream macrohabitats (riffles [R], pools connected to riffles [Pc], disconnected pools [Pd]) of 19 Mediterranean-climate sites in northern California to examine the influence of loss of habitat resulting from drought disturbance. At the time of sampling, 10 sites were perennial and included R and Pc macrohabitats, whereas 9 sites were intermittent and included only Pd macrohabitats. Taxa richness was more variable in Pd, and taxa richness was significantly lower in Pd than in Pc but not R. These results suggested a decline in richness between Pc and Pd that might be associated with loss of connectivity. Lower Ephemeroptera, Plecoptera, and Trichoptera (EPT) richness relative to Odonata, Coleoptera, and Heteroptera (OCH) richness was observed for Pd than R and Pc macrohabitats. Family composition was more similar between R and Pc than between R or Pc and Pd macrohabitats. This similarity may be associated with greater connectivity between R and Pc macrohabitats. Correspondence analysis indicated that macroinvertebrate composition changed along a gradient from R to Pc and Pd that was related to a perennialintermittent gradient across sites. High variability among macroinvertebrate assemblages in Pd could have been related to variability in the duration of intermittency. In cluster analysis, macroinvertebrate assemblages were grouped by macrohabitat first and then by site, suggesting that the macrohabitat filter had a greater influence on macroinvertebrate assemblages than did local site characteristics. Few taxa were found exclusively in Pc, and this macrohabitat shared numerous taxa with R and Pd, indicating that Pc may act as a bridge between R and Pd during drought. Drought is regarded as a ramp disturbance, but our results suggest that the response of macroinvertebrate assemblages to the loss of hydrological connectivity among macrohabitats is gradual, at least in Mediterranean-climate streams where drying is gradual. However, the changes may be more dramatic in arid and semiarid streams or in Mediterranean-climate streams if drying is rapid.

Integration of local-scale hydrological and regional-scale geophysical based on a nonlinear Bayesian sequential simulation approach

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed a downscaling procedure based on a non-linear Bayesian sequential simulation approach. The basic objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity, which is available throughout the model space. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariate kernel density function. This method is then applied to the stochastic integration of low-resolution, re- gional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this downscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the downscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.