Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

A comparison of wetness indices for the prediction of observed connected saturated areas under contrasting conditions

For lack of other widely available spatial information, topography is often used to predict water fluxes and water quality in mesoscale watersheds. Such data have however proven to be misleading in many environments where large and flat valley bottoms and/or highly conducive soil covers determine water storage and water transport mechanisms. Also, the focus is generally on the prediction of saturation areas regardless of whether they are connected to the catchment hydrographic network or rather present in isolated topographic depressions. Here soil information was coupled with terrain data towards the targeted prediction of connected saturated areas. The focus was on the 30 km² Girnock catchment (Cairngorm Mountains, northeast Scotland) and its 3 km² sub‐catchment, Bruntland Burn in which seven field surveys were done to capture actual maps of connected saturated areas in both dry and humid conditions. The 1 km² resolution UK Hydrology of Soil Types (HOST) classification was used to extract relevant, spatially variable, soil parameters. Results show that connected saturated areas were fairly well predicted by wetness indices but only in wet conditions when they covered more than 30% of the whole catchment area. Geomorphic indices including information on terrain shape, steepness, aspect, soil texture and soil depth showed potential but generally performed poorly. Indices based on soil and topographic data did not have more predictive power than those based on topographic information only: this was attributed to the coarse resolution of the HOST classification. Nevertheless, analyses provided interesting insights into the scale‐dependent water storage and transport mechanisms in both study catchments. Copyright © 2013 John Wiley & Sons, Ltd.     

Conceptual modelling to assess how the interplay of hydrological connectivity,catchment storage and tracer dynamics controls nonstationary water age estimates

Although catchment storage is an intrinsic control on the rainfall–runoff response of streams, direct measurement remains a major challenge. Coupled models that integrate long‐term hydrometric and isotope tracer data are useful tools that can provide insights into the dynamics of catchment storage and the volumes of water involved. In this study, we use a tracer‐aided hydrological model to characterize catchment storage as a dynamic control on system function related to streamflow generation, which also allows direct estimation of the nonstationarity of water ages. We show that in a wet Scottish upland catchment dominated by runoff generation from riparian peats (histosols) with high water storage, nonstationarity in water age distributions is only clearly detectable during more extreme wet and dry periods. This is explained by the frequency and longevity of hydrological connectivity and the associated relative importance of flow paths contributing younger or older waters to the stream. Generally, these saturated riparian soils represent large mixing zones that buffer the time variance of water age and integrate catchment‐scale partial mixing processes. Although storage simulations depend on model performance, which is influenced by input variability and the degree of isotopic damping in the stream, a longer‐term storage analysis of this model indicates a system that is only sensitive to more extreme hydroclimatic variability. Copyright © 2014 John Wiley & Sons, Ltd.     

Testing a spatially distributed tracer‐aided runoff model in a snow‐influenced catchment: Effects of multicriteria calibration on streamwater ages

Integrating stable isotope tracers into rainfall‐runoff models allows investigation of water partitioning and direct estimation of travel times and water ages. Tracer data have valuable information content that can be used to constrain models and, in integration with hydrometric observations, test the conceptualization of catchment processes in model structure and parameterization. There is great potential in using tracer‐aided modelling in snow‐influenced catchments to improve understanding of these catchments' dynamics and sensitivity to environmental change. We used the spatially distributed tracer‐aided rainfall‐runoff (STARR) model to simulate the interactions between water storage, flux, and isotope dynamics in a snow‐influenced, long‐term monitored catchment in Ontario, Canada. Multiple realizations of the model were achieved using a combination of single and multiple objectives as calibration targets. Although good simulations of hydrometric targets such as discharge and snow water equivalent could be achieved by local calibration alone, adequate capture of the stream isotope dynamics was predicated on the inclusion of isotope data in the calibration. Parameter sensitivity was highest, and most local, for single calibration targets. With multiple calibration targets, key sensitive parameters were still identifiable in snow and runoff generation routines. Water ages derived from flux tracking subroutines in the model indicated a catchment where runoff is dominated by younger waters, particularly during spring snowmelt. However, resulting water ages were most sensitive to the partitioning of runoff sources from soil and groundwater sources, which was most realistically achieved when isotopes were included in the calibration. Given the paucity of studies where hydrological models explicitly incorporate tracers in snow‐influenced regions, this study using STARR is an important contribution to satisfactorily simulating snowpack dynamics and runoff generation processes, while simultaneously capturing stable isotope variability in snow‐influenced catchments.     

West African palaeosynoptic patterns at the last glacial maximum

Summary The composite method and an atmospheric moisture budget is applied to study the present synoptic summer situation in the West African Sahel. This allows the information of the synoptic scale flow systems at the last glacial maximum to be obtained. Relying on the results of general circulation models and of local geological findings the palaeosynoptic situation was found to consist of the same constituents as today's given pattern. In both time slices the strength of the midtropospheric African Easterly Jet determines the intensity and frequency of the Sahelian synoptic disturbances: the Easterly Waves and the Squall Lines. The glacial rainfall intensity of the Squall Lines, the most effective rain-bearing system of the Sahel, is calculated to be 2/3 of the modern one. Additionally, only half the number of Squall Lines propagated across the West African continent per summer month. Thus the glacial rainfall amount of the Sahel was roughly 30% of the value observed today.With 6 Figures     

Spatial organization of groundwater dynamics and streamflow response from different hydropedological units in a montane catchment

Groundwater dynamics play an important role in runoff generation and hydrologic connectivity between hillslopes and streams. We monitored a network of 14 shallow groundwater (GW) wells in a 3.2 km² experimental catchment in the Scottish Highlands. Wells were placed in three contrasting landscape units with different hydropedological characteristics and different topographic positions relative to the stream network, encompassing a catena sequence from freely draining podzols on steeper hillslopes to increasingly thick peats (histosols) in the valley bottom riparian zone. GW dynamics were characterized by statistical analyses of water table fluctuations, estimation of variabilities in lag times and hysteresis response in relation to streamflow. The three landscape units had distinct storage–discharge relationships and threshold responses with a certain GW level above which lateral flow dominates. Steeper hillslopes with freely draining podzols were characterized by GW fluctuations of around 150 cm in the underlying drift. GW usually showed peak response up to several hours after stream flow. During persistent wet periods the water table remained in the soil profile for short spells and connected shallow flow paths in the near surface horizons to the lower hillslopes. In the peaty gleys in the lower foot slopes, GW was characterized by a water table generally within 20 cm of the soil surface, though at some locations this could fall to 50 cm in extreme dry periods. GW responses were usually a few hours prior to the stream responses. In riparian peats, the water table was also usually less than 20 cm deep and responded several hours before the stream. These riparian peat soils remain at, or very near saturation with near‐continuous GW–surface water connectivity. In contrast, the steeper slopes remain disconnected for prolonged periods and need large recharge events to overcome storage thresholds. GW responses vary seasonally, and landscape controls on the spatial organization of GW dynamics are strongest at low flows and in small events. During wettest periods, limited storage and extensive saturation weaken such controls. This study demonstrated that montane catchments can have highly dynamic GW stores, which are important in generating both storm flows and baseflows. Copyright © 2016 John Wiley & Sons, Ltd.     

Groundwater isoscapes in a montane headwater catchment show dominance of well‐mixed storage

We conducted an integrated groundwater–surface water monitoring programme in a 3.2‐km² experimental catchment in the Scottish Highlands by sampling all springs, seepages, and wells in six, spatially extensive synoptic surveys over a 2‐year period. The catchment has been glaciated, with steep hillslopes and a flat valley bottom. There is around 70% glacial drift cover in lower areas. The solid geology, which outcrops at higher elevations, is granite and metamorphic schist. The springs and seepages generally occur at the contact between the solid geology and drift or at breaks of slopes in the valley bottom. Samples were analysed for stable isotopes, Gran alkalinity and electrical conductivity. Despite the surveys encompassing markedly different antecedent conditions, the isotopic composition of groundwater at each location exhibited limited temporal variability, resulting in a remarkable persistence of spatial patterns indicating well‐mixed shallow, groundwater stores. Moreover, line‐conditioned excess values derived from the isotope data indicated no evidence of fractionation affecting the groundwater, which suggests that most recharge occurs in winter. The alkalinity and electrical conductivity of groundwater reflected geological differences in the catchment, being highest where more weatherable calcareous rocks outcrop at higher altitudes in the catchment. Springs draining these areas also had the most variable isotope composition, which indicated that they have shorter residence times than the drift covered part of the catchment. The study showed that even in geologically heterogeneous upland catchments, groundwater can be characterized by a consistent isotopic composition, reflecting rapid mixing in the recharge zone. Our work, thus, emphasizes the critical role of groundwater in upland catchments and provides tracer data that can help constrain quantitative groundwater models.     

Does the incorporation of process conceptualization and tracer data improve the structure and performance of a simple rainfall‐runoff model in a Scottish mesoscale catchment?

A geomorphological instantaneous unit hydrograph (GIUH) rainfall‐runoff model was applied in a 31 km² montane catchment in Scotland. Modelling was based on flow path length distributions derived from a digital terrain model (DTM). The model was applied in two ways; a single landscape unit response based on the DTM alone, and a two‐landscape unit response, which incorporated the distribution of saturated areas derived from field‐validated geographic information system (GIS) analysis based on a DTM and soil maps. This was to test the hypothesis that incorporation of process‐information would enhance the model performance. The model was applied with limited multiple event calibration to produce parameter sets which could be applied to a spectrum of events with contrasting characteristics and antecedent conditions. Gran alkalinity was used as a tracer to provide an additional objective measure for assessing model performance. The models captured the hydrological response dynamics of the catchment reasonably well. In general, the single landscape unit approach produced the best individual model performance statistics, though the two‐landscape unit approach provided a range of models, which bracketed the storm hydrograph response more realistically. There was a tendency to over‐predict the rising limb of the hydrograph, underestimate large storm event peaks and anticipate the hydrograph recession too rapidly. Most of these limitations could be explained by the simplistic assumptions embedded within the GIUH approach. The modelling also gave feasible predictions of stream water chemistry, though these could not be used as a basis for model rejection. Nevertheless, the study suggested that the approach has potential for prediction of hydrological response in ungauged montane headwater basins. Copyright © 2007 John Wiley & Sons, Ltd.