Tree‐ring records unveil long‐term influence of the Pacific Decadal Oscillation on snowpack dynamics in the Stikine River basin,northern British Columbia

Long hydroclimate records are essential elements for the assessment and management of changing freshwater resources. These records are especially important in transboundary watersheds where international cooperation is required in the joint planning and management process of shared basins. Dendrochronological techniques were used to develop a multicentury record of April 1 snow water equivalent (SWE) for the Stikine River basin in northern British Columbia, Canada, from moisture‐sensitive white spruce (Picea glauca) tree rings. Explaining 43% of the instrumental SWE variability, to our knowledge, this research represents the first attempt to develop long‐term snowpack reconstructions in northern British Columbia. The results indicated that 15 extreme low April 1 SWE events occurred from 1789 to the beginning of the instrumental record in 1974. The reconstruction record also shows that the occurrence of hydrological extremes in the Stikine River basin is characterized by persistent below‐average periods in SWE consistent with phase shifts of the Pacific Decadal Oscillation (PDO). Spectral analyses indicate a very distinct in‐phase (positive) relationship between the multidecadal frequencies of variability (~40 years) extracted from the SWE tree‐ring reconstruction and other reconstructed winter and spring PDO indices. Comparison of the reconstructed SWE record with other tree‐ring‐derived PDO proxy records shows coherence at multidecadal frequencies of variability. The research has significant implications for regional watershed management by highlighting the hydrological response of the Stikine River basin to prior climate changes.     

Regional variation of recession flow power‐law exponent

Recession flows of a basin provide valuable information about its storage–discharge relationship as during recession periods discharge occurs due to depletion of storage. Storage–discharge analysis is generally performed by plotting ?dQ/dt against Q , where Q is discharge at time t . For most real world catchments, ?dQ/dt versus Q show a power‐law relationship of the type: ?dQ/dt = kQ^α . Because the coefficient k varies across recession events significantly, the exponent α needs to be computed separately for individual recession events. The median α can then be considered as the representative α for the basin. The question that arises here is what are the basin characteristics that influence the value of α ? Studies based on a small number of basins (up to 50 basins) reveal that α has good relationship with several basin characteristics. However, whether such a relationship is universal remains an important question, because a universal relationship would allow prediction of the value of α for any ungauged basin. To test this hypothesis, here, we study data collected from a relatively large number of basins (358 basins) in USA and examine the influence of 35 different physio‐climatic characteristics on α . We divide the basins into 2 groups based on their longitudes and test the relationship between α and basin characteristics separately for the two groups. The results indicate that α is not identically influenced by different basin characteristics for the two datasets. This may suggest that the power‐law exponent α of a region is determined by the way local physio‐climatic forces have shaped the landscape.     

Assessment of climate change impacts on the hydrology of the Peruvian Amazon–Andes basin

In this article, we propose an investigation of the modifications of the hydrological response of two Peruvian Amazonas–Andes basins in relationship with the modifications of the precipitation and evapotranspiration rates inferred by the IPCC. These two basins integrate around 10% of the total area of the Amazonian basin. These estimations are based on the application of two monthly hydrological models, GR2M and MWB3, and the climatic projections come from BCM2, CSMK3 and MIHR models for A1B and B1 emission scenarios (SCE A1B and SCE B1). Projections are approximated by two simple scenarios (anomalies and horizon) and annual rainfall rates, evapotranspiration rates and discharge were estimated for the 2020s (2008–2040), 2050s (2041–2070) and 2080s (2071–2099). Annual discharge shows increasing trend over Requena basin (Ucayali river), Puerto Inca basin (Pachitea river), Tambo basin (Tambo river) and Mejorada basin (Mantaro river) while discharge shows decreasing trend over the Chazuta basin (Huallaga river), the Maldonadillo basin (Urubamba river) and the Pisac basin (Vilcanota river). Monthly discharge at the outlet of Puerto Inca, Tambo and Mejorada basins shows increasing trends for all seasons. Trends to decrease are estimated in autumn discharge over the Requena basin and spring discharge over Pisac basin as well as summer and autumn discharges over both the Chazuta and the Maldonadillo basins. Copyright © 2011 John Wiley & Sons, Ltd.     

Climate and hydrological changes in the Ob River Basin during 1936–2017

Changes in the hydrological regimes of Arctic rivers could affect the thermohaline circulation of the Arctic Ocean. In this study, we analysed spatiotemporal variations in temperature and precipitation in the Ob River Basin regions during 1936–2017 based on data from the Global Precipitation Climatology Center. Changes in discharge and response to climate change were examined based on monthly observed data during the same period. It is indicated the Ob River Basin experienced significant overall rapid warming and wetting (increased precipitation) in the study period, with average rates of 0.20°C (10 year⁻¹) and 5.3 mm (10 year⁻¹), respectively. The annual spatial variations of temperature and precipitation showed different scales in different regions. The discharge in spring and winter significantly increased at a rate of 384.1 and 173.1 m³/s (10 year⁻¹), respectively. Hydrograph separation indicated infiltration and supported that deep flow paths increased the contribution of groundwater to base flow. Meanwhile, the variation of the ratio of Q_max/Q_min suggested that the basin storage and the mechanism of discharge generation have significantly changed. The hydrological processes were influenced by changes of permafrost in a certain in the Ob River Basin. An increase in the recession coefficient (RC) implies that the permafrost degradation in the basin due to climate warming affected hydrological processes in winter. Permafrost degradation affected the Q_max/Q_min more significantly in the warm season than RC due to the enhanced infiltration that converted more surface water into groundwater in the cold season. The impact of precipitation on discharge, including surface flow and base flow, was more significant than temperature at the annual and seasonal scales in the Ob River Basin. The base flow was more obviously influenced by temperature than surface flow. The results of this study are significant for analyses of the basin water budget and freshwater input to the Arctic Ocean.     

Snow distribution from SSM/I and its relationships to the hydroclimatology of the Mackenzie River Basin,Canada

The spatial and temporal distribution of snow cover extent (SCE) and snow water equivalent (SWE) play vital roles in the hydrology of northern watersheds. We apply remotely sensed Special Sensor Microwave Imager (SSM/I) data from 1988 to 2007 to explore the relationships between snow distribution and the hydroclimatology of the Mackenzie River Basin (MRB) of Canada and its major sub-basins. The Environment Canada (EC) algorithm is adopted to retrieve the SWE from SSM/I data. Moderate Resolution Imaging Spectroradiometer (MODIS) 8-day maximum snow cover extent products (MOD10A2) are used to estimate the different thresholds of retrieved SWE from SSM/I to classify the land cover as snow or no snow for various sub-basins in the MRB. The sub-basins have varying topography and hence different thresholds that range from 10 mm to 30 mm SWE. The accuracy of snow cover mapping based on the combination of several thresholds for the different sub-basins reaches ≈ 90%. The northern basins are found to have stronger linear relationships between the date on which snow cover fraction (SCF) reaches 50% or when SWE reaches 50% and mean air temperatures, than the southern basins. Correlation coefficients between SCF, SWE, and hydroclimatological variables show the new SCF products from SSM/I perform better than SWE from SSM/I to analyze the relationships with the regional hydroclimatology. Statistical models relating SCF and SWE to runoff indicate that the SCF and SWE from EC algorithms are able to predict the discharge in the early snow ablation seasons in northern watersheds.     

The effect of El Niño Southern Oscillation cycles on the decadal scale suspended sediment behavior of a coastal dry‐summer subtropical catchment

Rivers display temporal dependence in suspended sediment–water discharge relationships. Although most work has focused on multi‐decadal trends, river sediment behavior often displays sub‐decadal scale fluctuations that have received little attention. The objectives of this study were to identify inter‐annual to decadal scale fluctuations in the suspended sediment–discharge relationship of a dry‐summer subtropical river, infer the mechanisms behind these fluctuations, and examine the role of El Niño Southern Oscillation climate cycles. The Salinas River (California) is a moderate sized (11 000 km²), coastal dry‐summer subtropical catchment with a mean discharge (Q_mean) of 11.6 m³ s^?1. This watershed is located at the northern most extent of the Pacific coastal North America region that experiences increased storm frequency during El Niño years. Event to inter‐annual scale suspended sediment behavior in this system was known to be influenced by antecedent hydrologic conditions, whereby previous hydrologic activity regulates the suspended sediment concentration–water discharge relationship. Fine and sand suspended sediment in the lower Salinas River exhibited persistent, decadal scale periods of positive and negative discharge corrected concentrations. The decadal scale variability in suspended sediment behavior was influenced by inter‐annual to decadal scale fluctuations in hydrologic characteristics, including: elapsed time since small (~0.1 × Q_mean), and moderate (~10 × Q_mean) threshold discharge values, the number of preceding days that low/no flow occurred, and annual water yield. El Niño climatic activity was found to have little effect on decadal‐scale fluctuations in the fine suspended sediment–discharge relationship due to low or no effect on the frequency of moderate to low discharge magnitudes, annual precipitation, and water yield. However, sand concentrations generally increased in El Niño years due to the increased frequency of moderate to high magnitude discharge events, which generally increase sand supply. Copyright © 2014 John Wiley & Sons, Ltd.     

Water temperature dynamics in High Arctic river basins

Despite the high sensitivity of polar regions to climate change and the strong influence of temperature upon ecosystem processes, contemporary understanding of water temperature dynamics in Arctic river systems is limited. This research gap was addressed by exploring high‐resolution water column thermal regimes for glacier‐fed and non‐glacial rivers at eight sites across Svalbard during the 2010 melt season. Mean water column temperatures in glacier‐fed rivers (0.3–3.2 °C) were lowest and least variable near the glacier terminus but increased downstream (0.7–2.3 °C km^–1). Non‐glacial rivers, where discharge was sourced primarily from snowmelt runoff, were warmer (mean: 2.9–5.7 °C) and more variable, indicating increased water residence times in shallow alluvial zones and increased potential for atmospheric influence. Mean summer water temperature and the magnitude of daily thermal variation were similar to those of some Alaskan Arctic rivers but low at all sites when compared with alpine glacierized environments at lower latitudes. Thermal regimes were correlated strongly (p < 0.01) with incoming short‐wave radiation, air temperature, and river discharge. Principal drivers of thermal variability were inferred to be (i) water source (i.e. glacier melt, snowmelt, groundwater); (ii) exposure time to the atmosphere; (iii) prevailing meteorological conditions; (iv) river discharge; (v) runoff interaction with permafrost and buried ice; and (vi) basin‐specific geomorphological features (e.g. channel morphology). These results provide insight into the potential changes in high‐latitude river systems in the context of projected warming in polar regions. We hypothesize that warmer and more variable temperature regimes may prevail in the future as the proportion of bulk discharge sourced from glacial meltwater declines and rivers undergo a progressive shift towards snow water and groundwater sources. Importantly, such changes could have implications for aquatic species diversity and abundance and influence rates of ecosystem functioning in high‐latitude river systems. Copyright © 2012 John Wiley & Sons, Ltd.     

Biogeochemistry of dissolved carbon,major, and trace elements during spring flood periods on the Ob River

Detailed knowledge of the flood period of Arctic rivers remains one of the few factors impeding rigorous prediction of the effect of climate change on carbon and related element fluxes from the land to the Arctic Ocean. In order to test the temporal and spatial variability of element concentration in the Ob River (western Siberia) water during flood period and to quantify the contribution of spring flood period to the annual element export, we sampled the main channel year round in 2014–2017 for dissolved C, major, and trace element concentrations. We revealed high stability (approximately ≤10% relative variation) of dissolved C, major, and trace element concentrations in the Ob River during spring flood period over a 1‐km section of the river channel and over 3 days continuous monitoring (3‐hr frequency). We identified two groups of elements with contrasting relationship to discharge: (a) DIC and soluble elements (Cl, SO₄, Li, B, Na, Mg, Ca, P, V, Cr, Mn, As, Rb, Sr, Mo, Ba, W, and U) negatively correlated (p < 0.05) with discharge and exhibited minimal concentrations during spring flood and autumn high flow and (b) DOC and particle‐reactive elements (Al, Fe, Ti, Y, Zr, Nb, Cs, REEs, Hf, Tl, Pb, and Th), some nutrients (K), and metalloids (Ge, Sb, and Te), positively correlated (p < 0.05) with discharge and showed the highest concentrations during spring flood. We attribute the decreased concentration of soluble elements with discharge to dilution by groundwater feeding and increased concentration of DOC and particle‐reactive metals with discharge to leaching from surface soil, plant litter, and suspended particles. Overall, the present study provides first‐order assessment of fluxes of major and trace elements in the middle course of the Ob River, reveals their high temporal and spatial stability, and characterizes the mechanism of river water chemical composition acquisition.     

Increasing discharge from the Mackenzie River system to the Arctic Ocean

The Mackenzie River, Canada's longest and largest river system, provides the greatest Western Hemisphere discharge to the Arctic Ocean. Recent reports of declining flows have prompted concern because (1) this influences Arctic Ocean salinity, stratification and polar ice; (2) a major tributary, the Peace River, has large hydroelectric projects, and further dams are proposed; and (3) the system includes the extensive and biodiverse Peace–Athabasca, Slave and Mackenzie deltas. To assess hydrological trends over the past century that could reflect climate change, we analysed historic patterns of river discharges. We expanded the data series by infilling for short gaps, calculating annual discharges from early summer‐only records (typical r² > 0.9), coordinating data from sequential hydrometric gauges (requiring r² > 0.8) and advancing the data to 2013. For trend detection, Pearson correlation provided similar outcomes to non‐parametric Kendall's τ and Spearman's ρ tests. There was no overall pattern for annual flows of the most southerly Athabasca River (1913–2013), while the adjacent, regulated Peace River displayed increasing flows (1916–2013, p < 0.05). These rivers combine to form the Slave River, which did not display an overall trend (1917–2013). The more northerly, free‐flowing Liard River is the largest tributary and displayed increasing annual flows (1944–2013, p < 0.01, ~3.5% per decade) because of increasing winter, spring, and summer flows, and annual maximum and minimum flows also increased. Following from the tributary contributions, the Mackenzie River flows gradually increased (Fort Simpson 1939–2013, p < 0.05, ~1.5% per decade), but the interannual patterns for the Liard and other rivers were correlated with the Pacific Decadal Oscillation, complicating the pattern. This conclusion of increasing river flows to the Arctic Ocean contrasts with some prior reports, based on shorter time series. The observed flow increase is consistent with increasing discharges of the large Eurasian Arctic drainages, suggesting a common northern response to climate change. Analyses of historic trends are strengthened with lengthening records, and with the Pacific Decadal Oscillation influence, we recommend century‐long records for northern rivers. Copyright © 2016 John Wiley & Sons, Ltd.     

Theoretical link between rainfall and flood magnitude

The rainfall–run‐off convolution integral is analytically solved for several models for the elementary hydrograph. These solutions can be combined with available rainfall frequency analyses to predict flood flows along streams for different recurrence intervals, using no free parameters for gauged streams and one estimable parameter for ungauged streams. Extreme discharge magnitudes at gauged sites can be typically estimated within a factor of two of actual records, using no historical data on extreme flows. The flow predictions reproduce several important characteristics of the flood phenomenon, such as the slope of the regression line between observed extreme flows and basin area on the conventional logQ versus logA plot. Importantly, for the models and data sets investigated, the storm duration of greatest significance to flooding was found to approximate the intrinsic transport timescale of the particular watershed, which increases with basin size. Thus, storms that deliver extraordinary amounts of rainfall over a particular time interval will most greatly activate basins whose time constants approximately equal that interval. This theoretical finding is supported by examination of the regional hydrological response to the massive storms of September 14, 2008, and April 28–30, 2017, which caused extraordinary record flooding of basins of about 5–100 km² and 500–4,000 km², respectively, but produced few records in basins that were larger or smaller than those ranges.     

Measuring water availability with limited ground data: assessing the feasibility of an entirely remote‐sensing‐based hydrologic budget of the Rufiji Basin,Tanzania, using TRMM,GRACE, MODIS,SRB, and AIRS

This study explores the feasibility of an entirely satellite remote sensing (RS)‐based hydrologic budget model for a ground data‐constrained basin, the Rufiji basin in Tanzania, from the balance of runoff (Q), precipitation (P), storage change (ΔS), and evapotranspiration (ET). P was determined from the Tropical Rainfall Measuring Mission, ΔS from the Gravity Recovery and Climate Experiment, and ET from the Moderate Resolution Imaging Spectroradiometer, the surface radiation budget, and the Atmosphere Infrared Radiation Sounder. Q was estimated as a residual of the water balance and tested against measured Q for a sub‐basin of the Rufiji (the Usangu basin) where ground measurements were available (R² = 0.58, slope = 1.9, root mean square error = 29 mm/month, bias = 14%). We also tested a geographical information system (GIS)‐driven (ArcCN‐runoff) runoff model (R² = 0.64, slope = 0.43, root mean square error = 39 mm/month). We conducted an error propagation analysis from each of the model's hydrologic components (P, ET, and ΔS). We find that the RS‐based model amplitude is most sensitive to ET and slightly less so to P, whereas the model's seasonal trends are most sensitive to ?S. Although RS–GIS‐driven models are becoming increasingly used, our results indicate that long‐term water resource assessment policy and management may be more appropriate than ‘instantaneous’ or short‐term water resource assessment. However, our analyses help develop a series of tools and techniques to progress our understanding of RS–GIS in water resource management of data‐constrained basins at the level of a water resource manager. Copyright © 2012 John Wiley & Sons, Ltd.     

Simulating pan‐Arctic runoff with a macro‐scale terrestrial water balance model

A terrestrial hydrological model, developed to simulate the high‐latitude water cycle, is described, along with comparisons with observed data across the pan‐Arctic drainage basin. Gridded fields of plant rooting depth, soil characteristics (texture, organic content), vegetation, and daily time series of precipitation and air temperature provide the primary inputs used to derive simulated runoff at a grid resolution of 25 km across the pan‐Arctic. The pan‐Arctic water balance model (P/WBM) includes a simple scheme for simulating daily changes in soil frozen and liquid water amounts, with the thaw–freeze model (TFM) driven by air temperature, modelled soil moisture content, and physiographic data. Climate time series (precipitation and air temperature) are from the National Centers for Environmental Prediction (NCEP) reanalysis project for the period 1980–2001. P/WBM‐generated maximum summer active‐layer thickness estimates differ from a set of observed data by an average of 12 cm at 27 sites in Alaska, with many of the differences within the variability (1σ) seen in field samples. Simulated long‐term annual runoffs are in the range 100 to 400 mm year^?1. The highest runoffs are found across northeastern Canada, southern Alaska, and Norway, and lower estimates are noted along the highest latitudes of the terrestrial Arctic in North America and Asia. Good agreement exists between simulated and observed long‐term seasonal (winter, spring, summer–fall) runoff to the ten Arctic sea basins (r = 0·84). Model water budgets are most sensitive to changes in precipitation and air temperature, whereas less affect is noted when other model parameters are altered. Increasing daily precipitation by 25% amplifies annual runoff by 50 to 80% for the largest Arctic drainage basins. Ignoring soil ice by eliminating the TFM sub‐model leads to runoffs that are 7 to 27% lower than the control run. The results of these model sensitivity experiments, along with other uncertainties in both observed validation data and model inputs, emphasize the need to develop improved spatial data sets of key geophysical quantities (particularly climate time series) to estimate terrestrial Arctic hydrological budgets better. Copyright © 2003 John Wiley & Sons, Ltd.     

Intercomparison of snow water equivalent observations in the Northern Great Plains

In the Northern Great Plains, melting snow is a primary driver of spring flooding, but limited knowledge of the magnitude and spatial distribution of snow water equivalent (SWE) hampers flood forecasting. Passive microwave remote sensing has the potential to enhance operational river flow forecasting but is not routinely incorporated in operational flood forecasting. We compare satellite passive microwave estimates from the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR‐E) to the National Oceanic and Atmospheric Administration Office of Water Prediction (OWP) airborne gamma radiation snow survey and U.S. Army Corps of Engineers (USACE) ground snow survey SWE estimates in the Northern Great Plains from 2002 to 2011. AMSR‐E SWE estimates compare favourably with USACE SWE measurements in the low relief, low vegetation study area (mean difference = ?3.8 mm, root mean squared difference [RMSD] = 34.7 mm), but less so with OWP airborne gamma SWE estimates (mean difference = ?9.5 mm, RMSD = 42.7 mm). An error simulation suggests that up to half of the error in the former comparison is potentially due to subpixel scale SWE variability, limiting the maximum achievable RMSD between ground and satellite SWE to approximately 26–33 mm in the Northern Great Plains. The OWP gamma versus AMSR‐E SWE comparison yields larger error than the point‐scale USACE versus AMSR‐E comparison, despite a larger measurement footprint (5–7 km² vs. a few square centimetres, respectively), suggesting that there are unshared errors between the USACE and OWP gamma SWE data.     

Modelling particle residence times in agricultural river basins using a sediment budget model and fallout radionuclide tracers

Contemporary patterns in river basin sediment dynamics have been widely investigated but the timescales associated with current sediment delivery processes have received much less attention. Furthermore, no studies have quantified the effect of recent land use change on the residence or travel times of sediment transported through river basins. Such information is crucial for understanding contemporary river basin function and responses to natural and anthropogenic disturbances or management interventions. To address this need, we adopt a process‐based modelling approach to quantify changes in spatial patterns and residence times of suspended sediment in response to recent agricultural land cover change. The sediment budget model SedNet was coupled with a mass balance model of particle residence times based on atmospheric and fluvial fluxes of three fallout radionuclide tracers (⁷Be, excess ²¹⁰Pb and ¹³⁷Cs). Mean annual fluxes of suspended sediment were simulated in seven river basins (38–920 km²) in south‐west England for three land cover surveys (1990, 2000 and 2007). Suspended sediment flux increased across the basins from 0.5–15 to 1.4–37 kt y^‐1 in response to increasing arable land area between consecutive surveys. The residence time model divided basins into slow (upper surface soil) and rapid (river channel and connected hillslope sediment source area) transport compartments. Estimated theoretical residence times in the slow compartment decreased from 13–48 to 5.6–14 ky with the increase in basin sediment exports. In contrast, the short residence times for the rapid compartment increased from 185–256 to 260–368 d as the modelled connected source area expanded with increasing sediment supply from more arable land. The increase in sediment residence time was considered to correspond to longer sediment travel distances linked to larger connected source areas. This novel coupled modelling approach provides unique insight into river basin responses to recent environmental change not otherwise available from conventional measurement techniques. © 2014 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Spatial interpolation of snow water equivalency using surface observations and remotely sensed images of snow‐covered area

As demand for water continues to escalate in the western Unites States, so does the need for accurate monitoring of the snowpack in mountainous areas. In this study, we describe a simple methodology for generating gridded‐estimates of snow water equivalency (SWE) using both surface observations of SWE and remotely sensed estimates of snow‐covered area (SCA). Multiple regression was used to quantify the relationship between physiographic variables (elevation, slope, aspect, clear‐sky solar radiation, etc.) and SWE as measured at a number of sites in a mountainous basin in south‐central Idaho (Big Wood River Basin). The elevation of the snowline, obtained from the SCA estimates, was used to constrain the predicted SWE values. The results from the analysis are encouraging and compare well to those found in previous studies, which often utilized more sophisticated spatial interpolation techniques. Cross‐validation results indicate that the spatial interpolation method produces accurate SWE estimates [mean R² = 0·82, mean mean absolute error (MAE) = 4·34 cm, mean root mean squared error (RMSE) = 5·29 cm]. The basin examined in this study is typical of many mid‐elevation mountainous basins throughout the western United States, in terms of the distribution of topographic variables, as well as the number and characteristics of sites at which the necessary ground data are available. Thus, there is high potential for this methodology to be successfully applied to other mountainous basins. Copyright © 2010 John Wiley & Sons, Ltd.     

Factors controlling the chemical evolution of travertine‐depositing rivers of the Barkly karst,northern Australia

Groundwaters feeding travertine‐depositing rivers of the northeastern segment of the Barkly karst (NW Queensland, Australia) are of comparable chemical composition, allowing a detailed investigation of how the rate of downstream chemical evolution varies from river to river. The discharge, pH, temperature, conductivity and major‐ion concentrations of five rivers were determined by standard field and laboratory techniques. The results show that each river experiences similar patterns of downstream chemical evolution, with CO₂ outgassing driving the waters to high levels of calcite supersaturation, which in turn leads to widespread calcium carbonate deposition. However, the rate at which the waters evolve, measured as the loss of CaCO₃ per kilometre, varies from river to river, and depends primarily upon discharge at the time of sampling and stream gradient. For example, Louie Creek (Q = 0·11 m³ s^?1) and Carl Creek (Q = 0·50 m³ s^?1) have identical stream gradients, but the loss of CaCO₃ per kilometre for Louie Creek is twice that of Carl Creek. The Gregory River (Q = 3·07 m³ s^?1), O'Shanassy River (Q = 0·57 m³ s^?1) and Lawn Hill Creek (Q = 0·72 m³ s^?1) have very similar gradients, but the rate of hydrochemical evolution of the Gregory River is significantly less than either of the other two systems. The results have major implications for travertine deposition: the stream reach required for waters to evolve to critical levels of calcite supersaturation will, all others things being equal, increase with increasing discharge, and the length of reach over which travertine is deposited will also increase with increasing discharge. This implies that fossil travertine deposits preserved well downstream of modern deposition limits are likely to have been formed under higher discharge regimes. Copyright © 2002 John Wiley & Sons, Ltd.     

Decoupling the effects of deforestation and climate variability in the Tapajós river basin in the Brazilian Amazon

Tropical river basins are experiencing major hydrological alterations as a result of climate variability and deforestation. These drivers of flow changes are often difficult to isolate in large basins based on either observations or experiments; however, combining these methods with numerical models can help identify the contribution of climate and deforestation to hydrological alterations. This paper presents a study carried out in the Tapaj?s River (Brazil), a 477,000 km² basin in South‐eastern Amazonia, in which we analysed the role of annual land cover change on daily river flows. Analysis of observed spatial and temporal trends in rainfall, forest cover, and river flow metrics for 1976 to 2008 indicates a significant shortening of the wet season and reduction in river flows through most of the basin despite no significant trend in annual precipitation. Coincident with seasonal trends over the past 4 decades, over 35% of the original forest (140,000 out of 400,000 km²) was cleared. In order to determine the effects of land clearing and rainfall variability to trends in river flows, we conducted hindcast simulations with ED2 + R, a terrestrial biosphere model incorporating fine scale ecosystem heterogeneity arising from annual land‐use change and linked to a flow routing scheme. The simulations indicated basin‐wide increases in dry season flows caused by land cover transitions beginning in the early 1990s when forest cover dropped to 80% of its original extent. Simulations of historical potential vegetation in the absence of land cover transitions indicate that reduction in rainfall during the dry season (mean of ?9 mm per month) would have had an opposite and larger magnitude effect than deforestation (maximum of +4 mm/month), leading to the overall net negative trend in river flows. In light of the expected increase in future climate variability and water infrastructure development in the Amazon and other tropical basins, this study presents an approach for analysing how multiple drivers of change are altering regional hydrology and water resources management.     

Study of dynamic behaviour of recession curves

Since Brutsaert and Neiber (1977), recession curves are widely used to analyse subsurface systems of river basins by expressing ? dQ/dt as a function of Q, which typically take a power law form: ? dQ/dt = kQ^α, where Q is the discharge at a basin outlet at time t. Traditionally recession flows are modelled by single reservoir models that assume a unique relationship between ? dQ/dt and Q for a basin. However, recent observations indicate that ? dQ/dt ‐ Q relationship of a basin varies greatly across recession events, indicating the limitation of such models. In this study, the dynamic relationship between ? dQ/dt and Q of a basin is investigated through the geomorphological recession flow model which models recession flows by considering the temporal evolution of its active drainage network (the part of the stream network of the basin draining water at time t). Two primary factors responsible for the dynamic relationship are identified: (i) degree of aquifer recharge (ii) spatial variation of rainfall. Degree of aquifer recharge, which is likely to be controlled by (effective) rainfall patterns, influences the power law coefficient, k. It is found that k has correlation with past average streamflow, which confirms the notion that dynamic ? dQ/dt ‐ Q relationship is caused by the degree of aquifer recharge. Spatial variation of rainfall is found to have control on both the exponent, α, and the power law coefficient, k. It is noticed that that even with same α and k, recession curves can be different, possibly due to their different (recession) peak values. This may also happen due to spatial variation of rainfall. Copyright © 2012 John Wiley & Sons, Ltd.     

Spatial patterns of suspended sediment yields in a humid tropical watershed in Costa Rica

An Erratum has been published for this article in Hydrological Processes 16(5) 2002, 1130–1131. Humid tropical regions are often characterized by extreme variability of fluvial processes. The Rio Terraba drains the largest river basin, covering 4767 km², in Costa Rica. Mean annual rainfall is 3139±419_sd mm and mean annual discharge is 2168±492_sd mm (1971–88). Loss of forest cover, high rainfall erosivity and geomorphologic instability all have led to considerable degradation of soil and water resources at local to basin scales. Parametric and non‐parametric statistical methods were used to estimate sediment yields. In the Terraba basin, sediment yields per unit area increase from the headwaters to the basin mouth, and the trend is generally robust towards choice of methods (parametric and LOESS) used. This is in contrast to a general view that deposition typically exceeds sediment delivery with increase in basin size. The specific sediment yield increases from 112±11·4_sd t km^?2 year^?1 (at 317·9 km² on a major headwater tributary) to 404±141·7_sd t km^?2 year^?1 (at 4766·7 km²) at the basin mouth (1971–92). The analyses of relationships between sediment yields and basin parameters for the Terraba sub‐basins and for a total of 29 basins all over Costa Rica indicate a strong land use effect related to intensive agriculture besides hydro‐climatology. The best explanation for the observed pattern in the Terraba basin is a combined spatial pattern of land use and rainfall erosivity. These were integrated in a soil erosion index that is related to the observed patterns of sediment yield. Estimated sediment delivery ratios increase with basin area. Intensive agriculture in lower‐lying alluvial fans exposed to highly erosive rainfall contributes a large part of the sediment load. The higher elevation regions, although steep in slope, largely remain under forest, pasture, or tree‐crops. High rainfall erosivity (>7400 MJ mm ha^?1 h^?1 year ^?1) is associated with land uses that provide inadequate soil protection. It is also associated with steep, unstable slopes near the basin mouth. Improvements in land use and soil management in the lower‐lying regions exposed to highly erosive rainfall are recommended, and are especially important to basins in which sediment delivery ratio increases downstream with increasing basin area. Copyright © 2001 John Wiley & Sons, Ltd.