Hydrologic characteristics of lake‐ and stream‐side riparian wetted margins in the McMurdo Dry Valleys,Antarctica

Water is a limiting factor for life in the McMurdo Dry Valleys (MDV), Antarctica. The active layer (seasonally thawed soil overlying permafrost) accommodates dynamic hydrological and biological processes for 10–16 weeks per year. Wetted margins (visually wetted areas with high moisture content) adjacent to lakes and streams are potential locations of great importance in the MDV because of the regular presence of liquid water, compared with the rest of the landscape where liquid water is rare. At 11 plots (four adjacent to lakes, seven adjacent to streams), soil particle size distribution, soil electrical conductivity, soil water content and isotopic signature, width of the wetted margin, and active layer thaw depth were characterised to determine how these gradients influence physicochemical properties that determine microbial habitat and biogeochemical cycling. Sediments were generally coarse‐grained in wetted margins adjacent to both lakes and streams. Wetted margins ranged from 1·04 to 11·01 m in average length and were found to be longer at lakeside sites than streamside. Average thaw depths ranged from 0·12 to 0·85 m, and were found to be deepest under lake margins. Lake margins also had much higher soil electrical conductivity, steeper topographic gradients, but more gradual soil moisture gradients than stream margins. Patterns of soil water δ¹⁸O and δD distribution indicate capillary action and evaporation from wetted margins; margin pore waters generally demonstrated isotopic enrichment with distance from the shore, indicating evaporation of soil water. Lake margin pore waters were significantly more negative in D_XS (D_XS = δD‐8δ¹⁸O) than streamside pore waters, indicating a longer history of evaporation there. Differences between lake and stream margins can be explained by the more consistent availability of water to lake margins than stream margins. Differences in margin characteristics between lakes and streams have important consequences for the microbial habitat of these margins and their functional role in biogeochemical cycling at these terrestrial–aquatic interfaces. Copyright © 2009 John Wiley & Sons, Ltd.     

Micro-hole and multigrain quartz luminescence dating of Paleodeltas at Lake Fryxell,McMurdo Dry Valleys (Antarctica), and relevance for lake history

Relict (perched) lacustrine deltas around the perennially ice-covered lakes in the Taylor Valley, Antarctica, imply that these lakes were up to 40 times larger in area than at present since the last glacial maximum (LGM). These deltas have been used to constrain ice-margin positions in Taylor Valley, and the boundaries of the proposed LGM ice-damned Glacial Lake Washburn. The timing of these high lake levels has depended on ¹⁴C chronologies of algal layers within relict lacustrine deltas. To provide additional geochronometric data for the post-LGM lake-level history, we applied photon-stimulated-luminescence (PSL) sediment dating to polymineral fine silt and sand-size quartz from 7 perched-delta and 3 active-delta sites of different elevations along 3 major meltwater streams entering Lake Fryxell. Our PSL dating of 4 quartz-sand samples from core tops in the seasonal ice-free moat of Lake Fryxell (elevation ∼18 m a.s.l.) and two core-top moat samples from the seasonal moat of Lake Vanda in nearby Wright Valley establish that adequate PSL clock zeroing (by daylight) occurs in regional, modern shoreline deposits. Minimum-age micro-hole PSL results from the moats are consistently near 100 a. Minimum-age micro-hole age estimates for the deltas range from ∼50 to 100 a near the present lake level up to 13.4 ± 1.3 ka at 240 m. These are systematically younger than the comparable, reservoir-uncorrected, ¹⁴C ages that range from 7 ka (cal yr BP) to 13 ka (cal yr BP) near lake level up to 20 ka (cal yr BP) at 220–240 m elevation. Our results indicate the occurrence of a dramatic discrepancy between PSL minimum-age and ¹⁴C age estimates that is presently unresolved.     

Seasonal sediment fluxes forcing supraglacial melting on the Wright Lower Glacier,McMurdo Dry Valleys,Antarctica

Cold‐based glaciers exist in low temperature and low humidity environments in which shortwave radiation is the largest source of energy to the glacier surface and the energy budget is very sensitive to the surface albedo. Consequently, the presence of relatively low volumes of debris on glacier surfaces has a significant impact on the timing, magnitude and rate of ablation at the surface. The aim of this study is to understand how the presence of sediment on the glacier surface at the start of the melt season can affect meltwater generation and delivery on a cold‐based glacier. A combination of field measurements, energy balance modelling and chemical mixing modelling were used on the Wright Lower Glacier, McMurdo Dry Valleys, Antarctica, between October 2005 and January 2006 to address this aim. In this system, sediment was transported onto the glacier surface during the winter months (March–October) by foehn winds, which reduced surface albedo at the start of the summer melt season. The areas of the glacier on which sediment accumulated began to melt earlier than other parts of the glacier and experienced a longer melt season. Over the study period, the total ablation on the dirty surfaces was nine times greater than for clean ice. Ablation on the dirty surfaces is dominated by melting, whereas sublimation dominates the clean ice. As the sediment was unevenly distributed over the glacier surface, the variation in melt amount and timing drove the development of a cryoconite hole system. Copyright © 2012 John Wiley & Sons, Ltd.     

Quantifying the altered hydrologic connectivity of forest roads resulting from decommissioning and relocation

The decommissioning of roads is occurring in many forest environments with the aim of reducing the negative impacts of road runoff on water quality and aquatic habitat. Works associated with decommissioning are expensive so prior assessment of the outcomes of various options is merited. This paper presents a method of quantifying the degree to which a road is hydrologically connected to the stream network and thus the likely impacts of constructing a road of different configurations upon water quality. The method permits comparisons between different road network management options and is useful for assessing the likely result of decommissioning works. Emphasis is placed on quantifying the uncertainty of key performance measures. The procedures developed here are an extension of the probabilistic ‘volume to breakthrough’ model recently formulated by Australian water quality researchers and allow the quantification of road/stream connectivity without the need for extensive parameterization. To demonstrate its utility, the model was applied to an actual road decommissioning and replacement project in southeast Australia. Road areas and drainage outlets were surveyed in the field and flow paths to streams derived from a 1 metre resolution LiDAR based digital elevation model. The results demonstrate that the actual road decommissioning examined in this case was unlikely to reduce runoff to the stream network and that the overall impact of the works in conditions of design storms are likely to result in a net reduction in stream water quality. Copyright © 2007 John Wiley & Sons, Ltd.     

Spatial and temporal patterns of snow accumulation and aerial ablation across the McMurdo Dry Valleys,Antarctica

Snow in the McMurdo Dry Valleys is a potential source of moisture for subnivian soils in a cold desert ecosystem. In a water‐limited environment, enhanced soil moisture is expected to provide more favourable conditions for subnivian soil communities. In addition, snow cover insulates the underlying soil from air temperature extremes. Quantifying the spatial and temporal patterns of seasonal snow accumulation and ablation is necessary to understand these dynamics. Repeat high‐resolution imagery acquired for the 2009–2010 austral summer was used to map the seasonal distribution of snow across Taylor and Wright valleys, Southern Victorialand, Antarctica. An edge detection algorithm was used to perform an object‐based classification of snow‐covered area. Coupled with topographic parameters obtained from a 30‐m digital elevation model, unique distribution patterns were characterized for five regions within the neighbouring valleys. Time lapses of snow distribution in each region provide insight into spatially variable aerial ablation rates (change in area of landscape covered by snow) across the region. A strong coastal to interior gradient of decreasing snow‐covered area was evident for both Taylor and Wright valleys. The surrounding regions of Lake Fryxell, Lake Hoare, Lake Bonney, Lake Brownworth, and Lake Vanda exhibited losses of snow‐covered area of 9.61 km² (?93%), 1.63 km² (?72%), 1.07 km² (?97%), 2.60 km² (?82%), and 0.25 km² (?96%), respectively, as measured from peak accumulation in October to mid‐January. Differences in aerial ablation rates within and across local regions suggest that both topographic variation and regional microclimates influence the ablation of seasonal snow cover. Copyright © 2012 John Wiley & Sons, Ltd.     

Development of push moraines in deeply frozen sediment adjacent to a cold-based glacier in the McMurdo Dry Valleys,Antarctica

Previous studies of push moraines have concluded that structural architecture is controlled by a combination of stress exerted by flowing ice and the rheology of the deforming sediment. However, the relationship between deformation processes and the thermal/hydrologic conditions within the sediment remains ambiguous. Using a combination of surface mapping, logging exposures and ground penetrating radar we examine the relationships between surface morphology, material properties and development of deformation structures in deeply frozen sediment that has been deformed by Joyce Glacier, a cold-based glacier in the McMurdo Dry Valleys. The structural architecture of the push moraines is characterized by a combination of brittle and ductile deformation structures that have produced a structurally complex pattern of thrust faults, low angle listric thrusts and recumbent folds that extend 400 m beyond the glacier margin. Deformation is ductile where the ice concentration exceeds c. 65% volume and predominantly brittle where the ice concentration is less than c. 65%. The change in rheology reflects transition in behaviour from the material having a predominantly frictional character when the ice is limited to pore spaces to a non-frictional character in which strength is primarily determined by the cohesive strength of the ice. This work shows that glaciotectonic deformation can occur in deeply frozen permafrost where there is no liquid pore water. We conclude that the presence of liquid porewater is not a necessary condition for the development of glaciotectonic deformation or for the formation of push moraines. © 2019 John Wiley & Sons, Ltd.     

Spatial–temporal variability and hydrologic connectivity of runoff generation areas in a North Alabama pasture—implications for phosphorus transport

This study delineated spatially and temporally variable runoff generation areas in the Sand Mountain region pasture of North Alabama under natural rainfall conditions, and demonstrated that hydrologic connectivity is important for generating hillslope response when infiltration‐excess (IE) runoff mechanism dominates. Data from six rainfall events (13·7–32·3 mm) on an intensively instrumented pasture hillslope (0·12 ha) were analysed. Analysis of data from surface runoff sensors, tipping bucket rain gauge and HS‐flume demonstrated spatial and temporal variability in runoff generation areas. Results showed that the maximum runoff generation area, which contributed to runoff at the outlet of the hillslope, varied between 67 and 100%. Furthermore, because IE was the main runoff generation mechanism on the hillslope, the data showed that as the rainfall intensity changed during a rainfall event, the runoff generation areas expanded or contracted. During rainfall events with high‐intensity short‐ to medium‐duration, 4–8% of total rainfall was converted to runoff at the outlet. Rainfall events with medium‐ to low‐intensity, medium‐duration were found less likely to generate runoff at the outlet. In situ soil hydraulic conductivity (k) was measured across the hillslope, which confirmed its effect on hydrologic connectivity of runoff generation areas. Combined surface runoff sensor and k‐interpolated data clearly showed that during a rainfall event, lower k areas generate runoff first, and then, depending on rainfall intensity, runoff at the outlet is generated by hydrologically connected areas. It was concluded that in IE‐runoff‐dominated areas, rainfall intensity and k can explain hydrologic response. The study demonstrated that only connected areas of low k values generate surface runoff during high‐intensity rainfall events. Identification of these areas would serve as an important foundation for controlling nonpoint source pollutant transport, especially phosphorus. The best management practices can be developed and implemented to reduce transport of phosphorus from these hydrologically connected areas. Copyright © 2009 John Wiley & Sons, Ltd.     

Investigating hydrologic connectivity and its association with threshold change in runoff response in a temperate forested watershed

Hydrological studies across varied climatic and physiographic regions have observed small changes in the ‘states of wetness’; based on average soil moisture, can lead to dramatic changes in the amount of water delivered to the stream channel. This non-linear behaviour of the storm response has been attributed to a critical switching in spatial organization of shallow soil moisture and hydrologic connectivity. However, much of the analysis of the role of soil moisture organization and connectivity has been performed in small rangeland catchments. Therefore, we examined the relationship between hydrologic connectivity and runoff response within a temperate forested watershed of moderate relief. We have undertaken spatial surveys of shallow soil moisture over a sequence of storms with varying antecedent moisture conditions. We analyse each survey for evidence of hydrologic connectivity and we monitor the storm response from the catchment outlet. Our results show evidence of a non-linear response in runoff generation over small changes in measures of antecedent moisture conditions; yet, unlike the previous studies of rangeland catchments, in this forested landscape we do not observe a significant change in geostatistical hydrologic connectivity with variations in antecedent moisture conditions. These results suggest that a priori spatial patterns in shallow soil moisture in forested terrains may not always be a good predictor of critical hydrologic connectivity that leads to threshold change in runoff generation, as has been the case in rangeland catchments. Copyright © 2007 John Wiley & Sons, Ltd.     

GLUE: 20 years on

This paper reviews the use of the Generalized Likelihood Uncertainty Estimation (GLUE) methodology in the 20 years since the paper by Beven and Binley in Hydrological Processes in (1992), which is now one of the most highly cited papers in hydrology. The original conception, the on‐going controversy it has generated, the nature of different sources of uncertainty and the meaning of the GLUE prediction uncertainty bounds are discussed. The hydrological, rather than statistical, arguments about the nature of model and data errors and uncertainties that are the basis for GLUE are emphasized. The application of the Institute of Hydrology distributed model to the Gwy catchment at Plynlimon presented in the original paper is revisited, using a larger sample of models, a wider range of likelihood evaluations and new visualization techniques. It is concluded that there are good reasons to reject this model for that data set. This is a positive result in a research environment in that it requires improved models or data to be made available. In practice, there may be ethical issues of using outputs from models for which there is evidence for model rejection in decision making. Finally, some suggestions for what is needed in the next 20 years are provided. © 2013 The Authors. Hydrological Processes published by John Wiley & Sons, Ltd.     

A temporal stable isotopic (δ18O, δD,d‐excess) comparison in glacier meltwater streams,Taylor Valley,Antarctica

In this paper, we describe the importance of hyporheic dynamics within Andersen Creek and Von Guerard Stream, Taylor Valley, Antarctica, from the 2010–2011 melt season using natural tracers. Water collection started at flow onset and continued, with weekly hyporheic‐zone sampling. The water δ¹⁸O and δD values were isotopically lighter in the beginning and heavier later in the season. D‐excess measurements were used as an indicator of mixing because an evaporative signature was evident and distinguishable between 2 primary end‐members (glacier meltwater and hyporheic zone). Hyporheic‐zone influence on the channel water was variable with a strong control on streamwater chemistry, except at highest discharges. This work supports previous research indicating that Von Guerard Stream has a large, widespread hyporheic zone that varies in size with time and discharge. Andersen Creek, with a smaller hyporheic zone, displayed hyporheic‐zone solute interaction through the influence from subsurface hypersaline flow. Overall, the evolution of Taylor Valley hyporheic‐zone hydrology is described seasonally. In mid‐December, the hyporheic zone is a dynamic system exchanging with the glacier meltwater in the channel, and with diminishing flow in January, the hyporheic zone drains back into the channel flow also impacting stream chemistry. This work adds new information on the role of hyporheic zone–stream interaction in these glacier meltwater streams.     

CHOSEN: A synthesis of hydrometeorological data from intensively monitored catchments and comparative analysis of hydrologic extremes

Comparative hydrology has been hampered by limited availability of geographically extensive, intercompatible monitoring data on comprehensive water balance stores and fluxes. These limitations have, for example, restricted comprehensive assessment of multiple dimensions of wetting and drying related to climate change and hampered understanding of why widespread changes in precipitation extremes are uncorrelated with changes in streamflow extremes. Here, we address this knowledge gap and underlying data gap by developing a new data synthesis product and using that product to detect trends in the frequencies and magnitudes of a comprehensive set of hydroclimatic and hydrologic extremes. CHOSEN (Comprehensive Hydrologic Observatory Sensor Network) is a database of streamflow, soil moisture, and other hydroclimatic and hydrologic variables from 30 study areas across the United States. An accompanying data pipeline provides a reproducible, semi-automated approach for assimilating data from multiple sources, performing quality assurance and control, gap-filling and writing to a standard format. Based on the analysis of extreme events in the CHOSEN dataset, we detected hotspots, characterized by unusually large proportions of monitored variables exhibiting trends, in the Pacific Northwest, New England, Florida and Alaska. Extreme streamflow wetting and drying trends exhibited regional coherence. Drying trends in the Pacific Northwest and Southeast were often associated with trends in soil moisture and precipitation (Pacific Northwest) and evapotranspiration-related variables (Southeast). In contrast, wetting trends in the upper Midwest and the Rocky Mountains showed few univariate associations with other hydroclimatic extremes, but their latitudes and elevations suggested the importance of changing snowmelt characteristics. On the whole, observed trends are incompatible with a ‘drying-in-dry, wetting-in-wet’ paradigm for climate-induced hydrologic changes over land. Our analysis underscores the need for more extensive, longer-term observational data for soil moisture, snow and evapotranspiration.     

Variation in hydrologic connectivity as a result of microtopography explained by discharge to catchment size relationship

In hydrological terms, raised bogs are often approximated by simple models as in the acrotelm–catotelm concept. However, raised bogs are often characterized by a pronounced surface topography, causing large changes in connectivity of contributing areas on the bog. In this study, daily regression of measured discharges versus catchment areas is used to quantify the impact of surface topography on catchment connectivity within a raised bog. The resulting coefficient of determination shows the strength of the relationship between the discharge and catchment area over time under different hydrological conditions. Monitoring of discharge, water table, transmissivity, and basic weather data on a raised bog (1.9 km²) in eastern central Estonia took place from May 2008 to June 2010. Contributing areas, calculated based on the outlet's discharge volume (V _Q ) divided by the net precipitation volume ( ), of the outlet containing the central pool‐ridge system varied between 1×10^?3 and 0.7 km², suggesting significant differences in connectivity between hydrological events. Correlation between discharge and theoretical catchment size was high (R ²>0.75) when the water table was close to the surface (less than 5 cm below peat surface), and consequently, transmissivities were also high (up to 1,030m²d^?1), which led to connectivity of local storage elements, such as pools and hollows. However, a water table below this threshold resulted in large parts of the catchment being disconnected. The importance of water table depths on catchment connectivity suggests the need to reconsider the hydrological concept of raised bogs; to incorporate these shallow flow components and better understand residence time and consequently transport of solutes, such as DOC, from patterned peatlands.     

近60年鄱阳湖东部湖湾水文连通变化及其对湿地植物与候鸟的影响

鄱阳湖碟形湖生物群落分布特征很大程度上受区域与主湖区水文连通性影响。由于水位波动,东部湖湾水文连通性受鄱阳湖主湖体水位影响较为敏感。本文以东部湖湾为例,分析主湖区水位变异程度及其影响,并从连通天数和发生时间等方面定量表征东部湖湾与主湖区的水文连通性,进而讨论水文连通性变化及其对湿地植物生境和候鸟栖息地的潜在影响。结果表明,当水位在13～16 m时,东部湖湾与主湖区存在着良好的水文连通关系,湿地植物适宜生境面积最大。2003年以后,鄱阳湖主湖区水位普遍降低,低枯水位持续时间延长,IHA/RVA法分析表明主湖区水位发生了中等程度改变,整体改变度为40.2%,东部湖湾与主湖区的连通关系发生明显改变,年连通天数减少了46.2 d,变化幅度为15.9%,而非连通期发生时间提前约1个月且年内时间跨度更长。水文连通性减弱造成3 10月东部湖湾平均水位下降了0.6 m,其中9 10月减少了约1 m,使薹草（Carexsp．）和苦草（Vallisneria natans）适宜水深对应的水面面积在多数月份呈增加趋势,特别是9 10月（幅度超过40%）,但枯水提前导致沉水植物面积减少,湿生植物生物量增加。水...     

Radio‐echo layering in West Antarctica: a spreadsheet dataset

Of the various information recovered from radio‐echo sounding (RES) of polar ice sheets, internal layering is currently under‐utilized by glaciologists, due in part to a lack of available data. Here, RES layering of the West Antarctic Ice Sheet, from the 1970s RES survey of approximately 70 per cent of this ice mass, is made available in a series of spreadsheets. Three types of internal layers are evident in the dataset. The first is continuous layers that have a stratigraphic appearance and can often be traced easily for hundreds of kilometres. The second is buckled layering, which also resembles stratigraphy and can sometimes be traced over tens of kilometres (although layer identification can often be difficult). The roughness of these layers is often greater than the bed at the same wavelength. The third is highly distorted or absent layering, which is not possible to trace laterally. Despite debate concerning the origin of RES layers, they are thought by most glaciologists to represent isochronous surfaces. The pattern of internal layering is potentially of importance to glaciologists for three reasons. (1) The position of undeformed layers below the ice surface is a function of accumulation rate, ice flow and basal melting conditions. Numerical modelling (including new ‘data assimilation’ techniques) could be used to discriminate between these processes, so revealing important information about the ice sheet and its environment. (2) Buckled layers are deformed by ice flow process, and so their occurrence can be related to the flow dynamics of the ice sheet. (3) Very buckled layers are often associated with ice stream flow, which allows their location to mark the positions of past and present fast‐flowing ice. Copyright © 2005 John Wiley & Sons, Ltd.     

Quantifying nitrogen cycling beneath a meander of a low gradient,N‐impacted,agricultural stream using tracers and numerical modelling

In watersheds impacted by nitrate from agricultural fertilizers, nitrification and denitrification may be decoupled as denitrification in the hyporheic zone is not limited to naturally produced nitrate. While most hyporheic research focuses on the 1–2 m of sediment beneath the stream bed, there are a limited number of studies that quantify nitrogen (N) cycling at larger hyporheic scales (10s of metres to kms). We conducted an investigation to quantify N cycling through a single meander of a low gradient, meandering stream, draining an agricultural watershed. Chemistry (major ions and N species) and hydrologic data were collected from the stream and groundwater beneath the meander. Evidence indicates that nearly all the shallow groundwater flowing beneath the meander originates as stream water on the upgradient side of the meander, and returns to the stream on the downgradient side. We quantified the flux of water beneath the meander using a numerical model. The flux of N into and out of the meander was quantified by multiplying the concentration of the important N species (nitrate, ammonium, dissolved organic nitrogen (DON)) by the modelled water fluxes. The flux of N into the meander is dominated by nitrate, and the flux of N out of the meander is dominated by ammonium and DON. While stream nitrate varied seasonally, ammonium and DON beneath the meander were relatively constant throughout the year. When stream nitrate concentrations are high (>2 mg litre^?1), flow beneath the meander is a net sink for N as more N from nitrate in stream water is consumed than is produced as ammonium and DON. When stream nitrate concentrations are low (<2 mg litre^?1), the flux of N entering is less than exiting the meander. On an annual basis, the meander hyporheic flow serves as a net sink for N. Copyright © 2007 John Wiley & Sons, Ltd.     

Influence of the hyporheic zone on the phosphorus dynamics of a large gravel‐bed river,Garonne River,France

Phosphorus (P) concentrations in sediments and in surface and interstitial water from three gravel bars in a large river (Garonne River, southern France) were measured daily, downstream of a wastewater treatment plant for a city of 740 000 inhabitants (Toulouse). Measurements were made of vertical hydraulic gradient (VHG), total dissolved phosphorus (TDP), soluble reactive phosphorus (SRP) and total phosphorus (TP) in water and of three extractable forms of phosphorus (water extractable, NaOH extractable and H₂SO₄ extractable) in hyporheic sediments from the gravel bars. Dissolved phosphorus was the major contributor to TP (74–79%) in both interstitial and surface waters on all sampling dates, and in most cases surface water P concentrations were significantly higher than interstitial concentrations. Hyporheic sediment TP concentrations ranged between 269 and 465 µg g^?1 and were highest in fine sediment fractions. Acid‐extractable P, a non‐bioavailable form, represented at least 95% of sediment TP. A positive relationship was observed between VHG and TP in two of the gravel bars, with wells that were strongly downwelling having lower TP concentrations. These results suggest that in downwelling zones, hyporheic sediments can trap surface‐derived dissolved P, and that much of this P becomes stored in refractory particulate forms. Bioavailable P is mainly present in dissolved form and only occupies a small fraction of total P, with particulate P comprising the majority of total P. Copyright © 2009 John Wiley & Sons, Ltd.     

Modelling the hydrologic effects of dynamic land‐use change using a distributed hydrologic model and a spatial land‐use allocation model

This study develops a novel approach for modelling and examining the impacts of time–space land‐use changes on hydrological components. The approach uses an empirical land‐use change allocation model (CLUE‐s) and a distributed hydrological model (DHSVM) to examine various land‐use change scenarios in the Wu‐Tu watershed in northern Taiwan. The study also uses a generalized likelihood uncertainty estimation approach to quantify the parameter uncertainty of the distributed hydrological model. The results indicate that various land‐use policies—such as no change, dynamic change and simultaneous change—have different levels of impact on simulating the spatial distributions of hydrological components in the watershed study. Peak flow rates under simultaneous and dynamic land‐use changes are 5·71% and 2·77%, respectively, greater than the rate under the no land‐use change scenario. Using dynamic land‐use changes to assess the effect of land‐use changes on hydrological components is more practical and feasible than using simultaneous land‐use change and no land‐use change scenarios. Furthermore, land‐use change is a spatial dynamic process that can lead to significant changes in the distributions of ground water and soil moisture. The spatial distributions of land‐use changes influence hydrological processes, such as the ground water level of whole areas, particularly in the downstream watershed. Copyright © 2010 John Wiley & Sons, Ltd.     

The effects of topographic depressions on multiscale overland flow connectivity: A high‐resolution spatiotemporal pattern analysis approach based on connectivity statistics

In watershed modelling, the traditional practice of arbitrarily filling topographic depressions in digital elevation models has raised concerns. Advanced high‐resolution remote sensing techniques, including airborne scanning laser altimetry, can identify naturally occurring depressions that impact overland flow. In this study, we used an ensemble physical and statistical modelling approach, including a 2D hydraulic model and two‐point connectivity statistics, to quantify the effects of depressions on high‐resolution overland flow patterns across spatial scales and their temporal variations in single storm events. Computations for both models were implemented using graphic processing unit‐accelerated computing. The changes in connectivity statistics for overland flow patterns between airborne scanning laser altimetry‐derived digital elevation models with (original) and without (filled) depressions were used to represent the shifts of overland flow response to depressions. The results show that depressions can either decrease or increase (to a lesser degree and shorter duration) the probability that any two points (grid locations) are hydraulically connected by overland flow pathways. We used macro‐connectivity states (Φ) as a watershed‐specific indicator to describe the spatiotemporal thresholds of connectivity variability caused by depressions. Four states of Φ are identified in a studied watershed, and each state represents different magnitudes of connectivity and connectivity changes (caused by depressions). The magnitude of connectivity variability corresponds to the states of Φ, which depend on the topological relationship between depressions, the rising/recession limb, and the total rainfall amount in a storm event. In addition, spatial distributions of connectivity variability correlate with the density of depression locations and their physical structures, which cause changes in streamflow discharge magnitude. Therefore, this study suggests that depressions are “nontrivial” in watershed modelling, and their impacts on overland flow should not be neglected. Connectivity statistics at different spatial scales and time points within a watershed provide new insights for characterizing the distributed and accumulated effects of depressions on overland flow.     

Assimilation of streamflow and in situ soil moisture data into operational distributed hydrologic models: Effects of uncertainties in the data and initial model soil moisture states

We assess the potential of updating soil moisture states of a distributed hydrologic model by assimilating streamflow and in situ soil moisture data for high-resolution analysis and prediction of streamflow and soil moisture. The model used is the gridded Sacramento (SAC) and kinematic-wave routing models of the National Weather Service (NWS) Hydrology Laboratory’s Research Distributed Hydrologic Model (HL-RDHM) operating at an hourly time step. The data assimilation (DA) technique used is variational assimilation (VAR). Assimilating streamflow and soil moisture data into distributed hydrologic models is new and particularly challenging due to the large degrees of freedom associated with the inverse problem. This paper reports findings from the first phase of the research in which we assume, among others, perfectly known hydrometeorological forcing. The motivation for the simplification is to reduce the complexity of the problem in favour of improved understanding and easier interpretation even if it may compromise the goodness of the results. To assess the potential, two types of experiments, synthetic and real-world, were carried out for Eldon (ELDO2), a 795-km² headwater catchment located near the Oklahoma (OK) and Arkansas (AR) border in the U.S. The synthetic experiment assesses the upper bound of the performance of the assimilation procedure under the idealized conditions of no structural or parametric errors in the models, a full dynamic range and no microscale variability in the in situ observations of soil moisture, and perfectly known univariate statistics of the observational errors. The results show that assimilating in situ soil moisture data in addition to streamflow data significantly improves analysis and prediction of soil moisture and streamflow, and that assimilating streamflow observations at interior locations in addition to those at the outlet improves analysis and prediction of soil moisture within the drainage areas of the interior stream gauges and of streamflow at downstream cells along the channel network. To assess performance under more realistic conditions, but still under the assumption of perfectly known hydrometeorological forcing to allow comparisons with the synthetic experiment, an exploratory real-world experiment was carried out in which all other assumptions were lifted. The results show that, expectedly, assimilating interior flows in addition to outlet flow improves analysis as well as prediction of streamflow at stream gauge locations, but that assimilating in situ soil moisture data in addition to streamflow data provides little improvement in streamflow analysis and prediction though it reduces systematic biases in soil moisture simulation.