A document by Mike Kirkby. Click on the document to view its contents.

A knowledge of the exchanges of energy and water over the terrestrial surface is the first step to understand the ecohydrological mechanisms, particularly in water-limited ecosystems in the dryland environments. However, patterns of energy exchange and evapotranspiration (ET) are not well understood in the oasis-desert ecotone, which plays an important role in protecting oasis against the threat of desertification in northwestern China’s arid regions. Here the continuous measurements of surface energy fluxes were made using eddy covariance in conjunction with auxiliary measurements for two years (2014-2015) at a shrubland within an oasis-desert ecotone in the arid regions, northwestern China. Statistical analysis on 30-min time scale indicates that about 50% of daytime net radiation (Rn) over the shrubland is dissipated as H on average, which peaks in spring; one third Rn is consumed by soil heat flux (G). Only 9% of Rn was consumed for latent heat flux (λE), which peaks in summer (21% in 2014 and 16% in 2015), corresponding to the season with highest rainfall among all seasons. Daily mean ET is about 1 mm·d−1 during growing season of the shrub species. The rapid and transient increase in ET occurs following a rainfall event. A switch in surface soil moisture from 0.04 to 0.11 m3·m−3 causes an increase in Rn by about 11% and λE by 151% at the shrubland, respectively. Accumulated annual ET were 195 and 181 mm in 2014 and 2015, respectively, exceeding the corresponding P by about 87 and 77 mm, indicating that groundwater may be another important source of water for ET over the shrubland aside from P. These results provide valuable insight into the mechanisms of sustaining energy and water balance at the ecotone, and then produce some management guidelines for allocating water resources and protecting vegetation.

Desert pavements are critical for maintaining ecological stability and promoting near-surface hydrological cycle in arid regions. However, few studies have reported the desert pavements on ecological on fluvial fans. Although desert pavement surfaces appear to be barren and flat, we found that the surfaces were featured by mosaic pattern of desert pavement (DP) and bare ground (BG). In this study, we investigated the effects of mosaic DP on water infiltration and vegetation distribution at six sites (i.e. one on the hillside and five in the sectors of fluvial fans) along a southwest belt transect on the fluvial fans in the Northern Linze County, in the middle of Hexi Corridor. The results showed that significant differences of Mosaic DP between hillside and sectors of fans were found in pavement thickness, thickness of vesicular horizon (Av thickness), particle composition and bulk density, rather than soil moisture content (SMC), gravel coverage and surface gravel size. The mosaic DP can inhibit water infiltration by pavement layer, where the sorptivity (S), initial infiltration rate (iint) and steady- state infiltration rate (isat) and infiltration time (T) averaged 1.30 cm/min-0.5, 5.03 cm/min, 0.23 cm/min, and 12.76 min respectively. If pavement layer was scalped, the S, iint and isat increased by 0.75 cm/min-0.5, 2.90 cm/min and 0.13 cm/min, respectively, and the T was shortened by 5.34 min. Water infiltration was mainly controlled by the pavement layer thickness (+), Av thickness (-), surface gravel coverage (-), and fine earth (+) and fine gravel (-) of pavement layer. Mosaic DP grew less shrubs than mosaic BG where distributed plenty of herbs. It can be concluded that desert pavements can keep vegetation stability by self-regulating rainfall. This study would deepen our understanding of the eco-hydrological cycle of pavement landscape in arid regions.

Complex networks of both natural and engineered flow paths control the hydrology of streams in major cities through spatio-temporal variations in connection and disconnection of water sources. We used spatially extensive and temporally intensive sampling of water stable isotopes to disentangle the hydrological sources of the heavily urbanized Panke catchment (≈ 220 km²) in the north of Berlin, Germany. The isotopic data enabled us to partition stream water sources across the catchment using a Bayesian mixing analysis. The upper part of the catchment streamflow here is dominated by groundwater from gravel aquifers underlying surrounding agricultural land. In dry summer periods, streamflow becomes intermittent; possibly as a result of local groundwater abstractions. Urban storm drainage is also an important part of runoff generation, dominating the responses to precipitation events. Although this dramatically changes the isotopic composition of the stream, it only accounts for 10-15% of annual streamflow. Moving downstream, subtle changes in sources and isotope signatures occur as catchment characteristic vary and the stream is affected by different tributary inflows. However, effluent from a wastewater treatment plant (WWTP) serving 700,000 people dominates the stream in the lower catchment where urbanisation effects are more dramatic. The associated increase in sealed surfaces downstream also reduces the relative contribution of groundwater to streamflow. The volume and isotopic composition of storm runoff is again dominated by urban drainage. As a result, only about 10% of annual runoff in the lower catchment comes from urban storm drains. The study shows the potential of stable water isotopes as inexpensive tracers in urban catchments that can provide a more integrated understanding of the complex hydrology of major cities. This offers an important evidence base for guiding the plans to develop and re-develop urban catchments to protect, restore and enhance the ecological and amenity value of these important resources.

Suspended sediment delivery and deposition in proglacial lakes is generally sensitive to a wide range of hydrometeorologic and geomorphic controls. High discharge conditions are of particular importance in many glaciolacustrine records, with individual floods potentially recorded as distinctive turbidites. We used an extensive network of surface sediment cores and hydroclimatic monitoring data to analyze recent flood turbidites and associated sediment transfer controls over instrumental periods at Eklutna Lake, western Chugach Mountains, Alaska. Close to a decade of fluvial data from primary catchment tributaries show a dominating influence of discharge on sediment delivery, with various interconnections with other related hydroclimatic controls. Multivariate fluvial models highlight and help quantify some complexities in sediment transfer, including intra-annual variations, meteorological controls, and the influence of subcatchment glacierization. Sediments deposited in Eklutna Lake during the last half century are discontinuously varved and contain multiple distinctive turbidites. Over a 30-year period of stratigraphic calibration, we correlate the four thickest flood turbidites (1989, 1995, 2006, 2012) to specific regional storms. The studied turbidites correlate with late-summer and early-autumn rainstorms with a magnitude of relatively instantaneous sedimentation 3 to 15 times greater than annual background accumulation. Our network of sediment core data captured the broad extent and sediment variability among the study turbidites and background sediment yield. Within-lake spatial modelling of deposition quantifies variable rates of downlake thinning and sediment focusing effects, and highlights especially large differences between the thickest flood turbidites and background sedimentation. This we primarily relate to strongly contrasting dispersion processes controlled by inflow current strength and turbidity. Sediment delivery is of interest for this catchment because of reservoir and water supply operations. Furthermore, although smaller floods may not be consistently represented, the lake likely contains a valuable proxy record of regional flooding proximal to major population centers of south-central Alaska including Anchorage.

The Ardières catchment, (150 km²), in Beaujolais (France), belongs to the first European catchments where surface water contamination by pesticides has been highlighted in the late 1980s. Research on this site mainly aims at better understanding organic pesticides and trace elements hydrological pathways to watercourses, and at evaluating subsequent contamination and impact on aquatic ecosystems. Landscape as well as instream processes are studied in order to highlight catchment vulnerability to contaminants and possible mitigation actions. A consistent hydrological and chemical monitoring of the Ardières River and one of its tributaries, the Morcille River, has been taking place since 2002. It was supplemented by biological measurements on aquatic micro-organisms and macroinvertebrates, more particularly after 2005. The results show the importance of long-term study to account for the kinetics of contaminant transfer in hydro-biogeochemical systems. Physico-chemical, ecological and ecotoxicological measurements all showed spatial and temporal variability in water quality and a gradient of impact on community structures and ecological functions as a function of the pressure of human activity. Results allowed to develop indicators of toxic impacts and resilience of communities and provided avenues for action to improve water quality.

Quantifying species-specific tree transpiration across watershed zones is important for estimating watershed evapotranspiration (ET) and predicting drought effects on vegetation. The objectives of this study are to 1) assess sap flux density (Js) and tree-level transpiration (Ts) across three contrasting zones (riparian buffer, mid-hillslope, and upland-hillslope), 2) determine how species-specific Js responds to vapor pressure deficit (VPD), and 3) compare watershed-level transpiration (Tw) derived from each zone. We measured Js and Ts in eight tree species in the three zones in a 12-ha forested watershed. In the dry year of 2015, loblolly pine (Pinus taeda), Virginia pine (Pinus virginiana), and sweetgum (Liquidambar styraciflua) Js rates were significantly higher in the buffer when compared to the other two zones. In contrast, Js in tulip poplar (Liriodendron tulipifera) and red maple (Acer rubrum) were significantly lower in the buffer than in the mid-hillslope. Daily Ts varied by zone and ranged from 10 to 93 liters in the dry year and 9 to 122 liters in the wet year. Js responded nonlinearly to VPD in all trees and zones. Annual Tw based on scaled-Js data was 447 mm, 377 mm, and 340 mm for the buffer, mid-hillslope, and upland-hillslope, respectively. We conclude that large spatial variability in Js and scaled Tw were driven by differences in soil moisture at each zone and forest composition. Consequently, spatial heterogeneity of vegetation and soil moisture must be considered when accurately quantifying watershed level ET.

Quantifying and understanding catchment sediment yields is crucial both from a scientific and environmental management perspective. To deepen the understanding of land use impacts and climate change on sediment load, we explore mechanisms of the suspended sediment yield formation in the Northern Caucasus during the Anthropocene. We examine how sediment flux of various river basins with different land-use/landcover and glacier cover changes during the 1925-2018 period. Our analysis is based on observed mean annual suspended sediment discharges (SSD, kg·s−1) and annual fluxes (SSL, t·yr−1) from 33 Roshydromet gauging stations (Russia). SSL series have been analyzed to detect statistically significant changes during the 1925-2018 period. The occurrence of abrupt change points in SSD was investigated using cumulative sum (CUSUM) charts. We found that SSL has decreased by −1.81% per year on average at most gauges. However, the decline was not linear. Several transition years are expected in the region: increasing trends from the 1950s and decreasing trends from 1988-1994. Correlation analyses showed that variation in SSL trend values is mainly explained by gauging station altitude, differences in land use (i.e., the fraction of cropland), and catchment area. Nonetheless, more accurate quantifications of SSL trend values and more refined characterizations of the catchments regarding (historical) land use, soil types/lithology, weather conditions, and topography may reveal other tendencies.

Tropical montane cloud forests (TMCF) receive additional (‘occult’) inputs of water from fog and wind-driven rain. Together with the concomitant reduction in evaporative losses, this typically leads to high soil moisture levels (often approaching saturation) that are likely to promote rapid subsurface flow via macropores. Although TMCF make up an estimated 6.6% of all remaining montane tropical forest and occur mostly in steep headwater areas that are protected in the expectation of reduced downstream flooding, TMCF hillslope hydrological functioning has rarely been studied. To better understand the hydrology of a TMCF on heterogeneously layered volcanic ash soils under very high net rainfall conditions (up to 6,740 mm yr-1), we examined temporal and spatial soil water dynamics and their contribution to subsurface runoff and stormflow for a year (1 July 2003 – 30 June 2004) in a small headwater catchment on the Atlantic (windward) slope near Monteverde, NW Costa Rica. Particular attention was paid to partitioning of water fluxes into lateral subsurface flow and vertical percolation. The presence of a gravelly volcanic layer at shallow depth of very high hydraulic conductivity intercalated between two layers of lower conductivity, controlled both surface infiltration and delayed vertical water movement deeper into the soil profile. Soil water fluxes during rainfall were dominated by rapid lateral flow in the gravelly layer, particularly at high soil water contents. In turn, this lateral subsurface flow controlled the magnitude and timing of stormflow from the catchment. Stormflow amount increased rapidly once topsoil moisture exceeded a threshold value of ~0.58 cm3 cm-3. Responses were not affected appreciably by rainfall intensity because soil hydraulic conductivities across the profile largely exceeded prevailing rainfall intensities.

Bushfires are common in Australia and can cause vegetation loss and affect hydrological processes such as interception, evapotranspiration, soil water storage and streamflow. This study investigates bushfire impacts on catchment mean annual streamflow for 14 Australian catchments that have been severely impacted by the 2009 Victoria bushfire, the second worst bushfire disaster in Australia. A statistical approach based on sensitivity coefficients was used for quantifying the climate variability impacts on streamflow and the time-trend analysis method was used to estimate the annual streamflow changes due to bushfire respectively. Our results show that bushfire has caused a noticeable increase in mean annual streamflow in majority of burnt catchments for an immediate post bushfire period (2009-2015), when the bushfire impact on streamflow is evidently larger than the climate change impact. Furthermore, the bushfire impact on mean annual streamflow linearly increases with the burnt percentage area. These results provide strong evidence for evaluating large-scale bushfire impact on streamflow at small to medium sized catchments, and guidance for process-based hydrological models for simulating bushfire impacts on hydrological processes for the immediate period after bushfire.

The assessment of climate change impacts on water resources and flood risk is typically underpinned by hydrological models calibrated and selected based on observed streamflow records. Yet, changes in climate are rarely accounted for when selecting hydrological models, which compromises their ability to robustly represent future changes in catchment hydrology. In this paper, we test a simple framework for selecting an ensemble of calibrated hydrological model structures in catchments where changing climatic conditions have been observed. We start by considering 78 model structures produced using the FUSE modular modelling framework and rely on a Pareto scheme to select model structures maximizing model efficiency in both wet and dry periods. The application of this approach in three case study basins in Peru enables the identification of structures with good robustness, but also good performance according to hydrological signatures not used for model selection. We also highlight that some model structures that perform well according to traditional efficiency metrics have low performance in contrasting climates or suspicious internal states and fluxes. Importantly, the model selection approach followed here helps to reduce the spread in precipitation elasticities and temperature sensitivities, providing a clearer picture of future hydrological changes. Overall, this work demonstrates the potential of using contrasting climatic conditions in a multi-objective framework to produce robust and credible simulations, and to constrain structural uncertainties in hydrological projections.

The δ2H and δ18O composition of 77 precipitation samples collected between January 2014 and April 2019 from two sites across the Guadalquivir Basin, SW Spain, were analyzed. Sampling was performed within intervals of at least 14 days basis if rain occurred but frequently intervals were longer according to the rainfall incidence. Precipitation weighted averages and local meteoric water lines are presented for use in hydrological applications. δ2H and δ18O values show a remarkably high variability attributed to the influence of Atlantic and Mediterranean vapor sources. Precipitation weighted average d-excess values of 12.5-13 ‰ confirm the Western Mediterranean influence. Temperature and amount effects were found to be weak but significant influence of secondary evaporation for single rainfall events during summer was identified by reduced d-excess values and enriched isotopic signatures plotting close or below the Global Meteoric Water line (GMWL). Isotopic signatures of both sites are very similar in general and any temperature related urban effects of Seville city compared to the rural site Doñana could not be identified with the present data

The main objective of the study is to evaluate the roles of climate change and human factors on runoff, baseflow, and hydrological drought characteristics at a watershed scale. The novelty of the study is to assess separately the cascading, indirect, accumulative effects of climate change and human factors on hydrological drought, i.e. runoff and baseflow. This involved analyzing change points to divide the available hydrometeorological data into a baseline and a perturbed period. We applied two hydrological models, SWAT and HBV-light, and two nonparametric climate elasticity of runoff to identify the contribution of climate change and human factors in influencing runoff and baseflow processes. The hydrological models were used to simulate naturalized runoff and baseflow during the perturbed period. The temporal variation in the characteristics of the baseflow regime is expressed as baseflow index. Drought indices, standardized runoff index and standardized baseflow index were used as hydrological drought indicators. A significant change in runoff reduction in the Kamienna watershed began in 1982, suggesting that human factors play a dominant role in influencing runoff. In addition, we found that an increase in baseflow and a decrease in hydrological drought events in the 2010s are a positive long-term effect of human factors such as construction of dams in the watershed. Finally, analyses of changes in land cover dynamics in the watershed over the past four decades using satellite imagery are used to confirm the presence of human interventions.

Despite the importance of headwater catchments for western United States’ water supply, these regions are often poorly understood, particularly with respect to quantitative understanding of evapotranspiration (ET) fluxes. Heterogeneity of land cover, physiography, and atmospheric patterns in these high-elevation regions lead to difficulty in developing spatially-distributed characterization of ET. As the largest terrestrial water flux behind precipitation, ET represents a significant fraction of the water budget for any watershed. Likewise, groundwater is the largest available freshwater store and has been shown to play a large role in the water balance, even in headwater systems. Using an eddy covariance tower in the East River Catchment, a Colorado River headwaters basin, this study estimates water and energy fluxes in high-elevation, complex systems to better constrain ET estimates and calculate overall water and energy budgets, including losses from groundwater. The eddy covariance method is used to estimate ET from years 2017 through 2019 at a saturated, riparian end-member site. Owing to complexities in near surface atmospheric structure such as stable boundary layers over snowpack and shallow terrain driven flow from surrounding landscape features, energy flux and ET estimates were limited to the warm season when energy closure residuals from the eddy-covariance system were reliably less than 30 %, a threshold commonly used in eddy covariance energy flux estimation. The resulting ET estimations are useful for constraining water budget estimates at this energy-limited site, which uses groundwater for up to 84 % of ET in the summer months. We also compared East River ET magnitudes and seasonality to two other flux towers (Niwot Ridge, CO and Valles Caldera, NM), located in the Rocky Mountains. This data is useful for constraining ET estimates in similar end-member locations across the East River Catchment. Our results show that groundwater-fed ET is a significant component of the water balance and groundwater may supply riparian ET even during low-snow years.

Delays in forest recovery from terrestrial acidification combined with climate change is leading Acadian Forest ecosystems into new territory. Kejimkujik Calibrated Catchments (KCC) Study Program was established in an around Kejimkujik National Park and Historic Site (KNPHS) in Southwest Nova Scotia (SWNS) in the late 1970s to increase our understanding of the impacts of acid precipitation on relatively pristine ecosystems. KCC now have one of the longest continuously monitored water chemistry records in North America, with data collection beginning in 1980. Its infrastructure includes three gauged streams, twelve forest inventory plots, an atmospheric deposition monitoring station, and three streams with continuous water quality monitoring and regular lab analysis of stream chemistry, and recent LiDAR coverage. The KCC fits into a wider network of monitored lakes. Data collected at the KCC form a key datapoint in comparisons of catchment response to terrestrial acidification in the context of a warming climate, due to their high and increasing DOC levels, highly dilute waters, lowland topography and extensive wetlands. KCC are also emerging as an important source of information for species at risk protection as SWNS was declared one of the 11 national priority places for biodiversity protection.

Agricultural land use can increase runoff and erosion leading to detrimental downstream impacts. This paper examines the impact of agricultural land cover change on runoff and sediment generation at event scales using a model-based approach. SHETRAN, a physically based, spatially distributed model, was applied in two southwest England catchments to represent: (a) changes in the land cover (cropland extent and spatial arrangement), (b) changes in crop type, and (c) use of riparian buffer strips. A total of 84 simulated events within a 4-year period were used to quantify flow and sediment generation impacts. We found past changes in land cover resulted in significant differences in sediment yield (p < 0.05). Linear regression showed an increase in flow and sediment yield proportional to increases in arable crop area (p < 0.001). The spatial arrangement of cropped fields and riparian buffer strips produced no significant differences in event flow or sediment yield (p > 0.05). However, buffer strip scenarios compared with the base run showed sediment load reductions in specific events, up to 20% and 15% for woodland and grass riparian buffers, respectively. When comparing crop types with and without the use of post-harvest cover crops, we observed non-significant differences (both Qv and Sy). However, large reductions in modelled sediment yields occurred for some events (e.g., up to 60% for winter cereals, 50% for maize and 74% for spring cereals). For these scenarios, examination of rainfall event magnitude emphasized the importance of ground cover in mitigating soil erosion for maize and spring cereals, but not for winter cereals. Our findings indicate that significant changes in sediment delivery at the event scale are associated with cropland extent and crop types, depending on rainfall magnitude, but not on the spatial arrangement of cropped fields or the use of riparian buffer strips.

Flow regimes are critical for determining physical and biological processes in rivers, and their classification and regionalization traditionally seeks to link patterns of flow to physiographic, climate and other information. There are many approaches to, and rationales for, catchment classification, with those focused on streamflow often seeking to relate a particular response characteristic to a physical property or climatic driver. Rationales include such topics as Prediction in Ungauged Basins (PUB), helping with experimental approaches, and providing guidance for model selection in poorly understood hydrological systems. While scale and time are important considerations for classification, the Annual Daily Hydrograph (ADH) is a first-order easily visualized integrated expression of catchment function, and over many years is a distinct hydrological signature. In this study, we use t-SNE, a state-of-the-art technique of dimensionality reduction, to classify 17110 ADHs for 304 reference catchments in mountainous Western North America. t-SNE is chosen over other conventional methods of dimensionality reduction (e.g. PCA) as it presents greater separability of ADHs, which are projected on a 2D map where the similarities are evaluated according to their map distance. We then utilize a Deep Learning encoder to upgrade the non-parametric t-SNE to a parametric approach, enhancing its capability to address ‘unseen’ samples. Results showed that t-SNE was an effective classifier as it successfully clustered ADHs of similar flow regimes on the 2D map. In addition, many compact clusters on the 2D map in the coastal Pacific Northwest suggest information redundancy in the local hydrometric network. The t-SNE map provides an intuitive way to visualize the similarity of high-dimensional data of ADHs, groups catchments with like characteristics, and avoids the reliance on subjective hydrometric indicators.

Gravel augmentation has become common practice to mitigate the effects of decline in upstream sediment supply in gravel-bed rivers. The success of such rehabilitation schemes relies partly on the monitoring strategy and efforts. When long-term monitoring is lacking, some aspects of rehabilitation initiatives suffer more than others, such as insights into functions and functionalities of the river system. Despite temperature being a fundamental parameter determining the general health of river ecosystems, a limited number of studies have tested whether gravel-augmentation can aid restoring thermal functions. With the help of airborne thermal infrared (TIR) imagery, this paper explores the potential positive feedbacks through the monitoring of gravel augmentation actions, of different magnitude, taken on 3 rivers of the Rhône basin in France. A specific trajectory-based Before-After-Control-Impact (BACI) framework using simple indicators, combined with a TIR-based Control-Impact strategy, was designed to assess the success of thermal function restorations based on dynamic fuzzy references. Results indicate that restoring forms is not sufficient to restore thermal functions. The control-impact strategy shows limitations in the sense that two neighbouring reaches can display similar planform characteristics but different thermal functions; what is observed in a control reach should not necessarily be expected following rehabilitation. When assessing thermal processes, a before-after strategy is needed to either serve as a target or help define an adequate target in accordance with changes in the catchment and channel adjustments and responsiveness. We therefore recommend a trajectory-based BACI assessment to identify current biogeophysical conditions within which rehabilitation can be assessed. From a technical perspective, airborne TIR proved to be useful to rapidly map surface temperature over dozens of kilometres at high resolution, and can be advocated as a powerful tool to monitor and diagnose thermal functions of gravel-bed rivers. With an increasing number of rehabilitation schemes, and increasing pressure of global changes on rivers, we suggest that monitoring of water temperature, even with simple but well-designed sampling strategies, becomes a routine part of river rehabilitation projects.