Including spatial distribution in a data‐driven rainfall‐runoff model to improve reservoir inflow forecasting in Taiwan

Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.     

Mid‐Miocene initiation of E‐W extension and recoupling of the Himalaya

The Tibetan plateau is host to numerous ~N‐S striking graben that have accommodated E‐W directed extension. The development of these structures has been interpreted to reflect a variety of different geological processes including plateau collapse, oroclinal bending or mid‐to‐lower crustal flow. New ⁴⁰Ar/³⁹Ar thermochronology and quartz c‐axis data from the Thakkhola graben of west‐central Nepal show that E‐W extension was ongoing at least locally by the early Miocene (ca. 17 Ma). Our new, and previously published chronologic information on the initiation of graben across the orogen shows that they typically developed immediately after cessation of the South Tibetan detachment system, a structural network that facilitated differential southward movement of the upper and middle crust. We interpret this fundamental switch in orogen kinematics to reflect recoupling of the middle and upper Himalayan crust such that the subsequent widespread flow of the mid‐to‐lower crust out of the system to the east forced brittle accommodation in the upper crust.     

Late Holocene environmental conditions in Coronation Gulf,southwestern Canadian Arctic Archipelago: evidence from dinoflagellate cysts,other non‐pollen palynomorphs,and pollen

Boxcore 99LSSL‐001 (68.095° N, 114.186° W; 211 m water depth) from Coronation Gulf represents the first decadal‐scale marine palynology and late Holocene sediment record for the southwestern part of the Northwest Passage. The record was studied for organic‐walled microfossils (dinoflagellate cysts, non‐pollen palynomorphs), pollen, terrestrial spores, and sediment characteristics. ²¹⁰Pb, ¹³⁷Cs, and three accelerator mass spectrometry ¹⁴C dates constrain the chronology. Three prominent palaeoenvironmental zones were identified. During the interval AD 1470–1680 (Zone I), the climate was warmer and wetter than at present, and environmental conditions were more favourable to biological activity and northward boreal forest migration, with reduced sea‐ice and a longer open‐water (growing) season. The interval AD 1680–1940 (Zone II) records sea‐ice increase, and generally cool, polar conditions during the Little Ice Age. During AD 1940–2000 (Zone III), organic microfossils indicate an extended open‐water season and decreased sea‐ice, with suggested amelioration surpassing that of Zone I. Although more marine studies are needed to place this record into an appropriate context, the succession from ameliorated (Zone I) to cooler, sea‐ice influenced conditions (Zone II) and finally to 20th‐century warming (Zone III) corresponds well with several terrestrial climatic records from the neighbouring mainland and Victoria Island, and with lower‐resolution marine records to the west. Copyright © 2011 John Wiley & Sons, Ltd.     

Synthesis of soil‐hydraulic properties and infiltration timescales in wildfire‐affected soils

We collected soil‐hydraulic property data from the literature for wildfire‐affected soils, ash, and unburned soils. These data were used to calculate metrics and timescales of hydrologic response related to infiltration and surface runoff generation. Sorptivity (S) and wetting front potential (Ψ_f) were significantly different (lower) in burned soils compared with unburned soils, whereas field‐saturated hydraulic conductivity (K_fs) was not significantly different. The magnitude and duration of the influence of capillarity during infiltration was greatly reduced in burned soils, causing faster ponding times in response to rainfall. Ash had large values of S and K_fs but moderate values of Ψ_f, compared with unburned and burned soils, indicating ash has long ponding times in response to rainfall. The ratio of S²/K_fs was nearly constant (~100 mm) for unburned soils but more variable in burned soils, suggesting that unburned soils have a balance between gravity and capillarity contributions to infiltration that may depend on soil organic matter, whereas in burned soils the gravity contribution to infiltration is greater. Changes in S and K_fs in burned soils act synergistically to reduce infiltration and accelerate and amplify surface runoff generation. Synthesis of these findings identifies three key areas for future research. First, short timescales of capillary influences on infiltration indicate the need for better measurements of infiltration at times less than 1 min to accurately characterize S in burned soils. Second, using parameter values, such as Ψ_f, from unburned areas could produce substantial errors in hydrologic modeling when used without adjustment for wildfire effects, causing parameter compensation and resulting underestimation of K_fs. Third, more thorough measurement campaigns that capture soil‐structural changes, organic matter impacts, quantitative water repellency trends, and soil‐water content along with soil‐hydraulic properties could drive the development of better techniques for numerically simulating infiltration in burned areas.     

Watershed structural influences on the distributions of stream network water and solute travel times under baseflow conditions

Watershed structure influences the timing, magnitude, and spatial location of water and solute entry to stream networks. In turn, stream reach transport velocities and stream network geometry (travel distances) further influence the timing of export from watersheds. Here, we examine how watershed and stream network organization can affect travel times of water from delivery to the stream network to arrival at the watershed outlet. We analysed watershed structure and network geometry and quantified the relationship between stream discharge and solute velocity across six study watersheds (11.4 to 62.8 km²) located in the Sawtooth Mountains of central Idaho, USA. Based on these analyses, we developed stream network travel time functions for each watershed. We found that watershed structure, stream network geometry, and the variable magnitude of inputs across the network can have a pronounced affect on water travel distances and velocities within a stream network. Accordingly, a sample taken at the watershed outlet is composed of water and solutes sourced from across the watershed that experienced a range of travel times in the stream network. We suggest that understanding and quantifying stream network travel time distributions are valuable for deconvolving signals observed at watershed outlets into their spatial and temporal sources, and separating terrestrial and in‐channel hydrological, biogeochemical, and ecological influences on in‐stream observations. Copyright © 2016 John Wiley & Sons, Ltd.     

Linkages between spatio‐temporal patterns of environmental factors and distribution of plant assemblages across a boreal peatland complex

Here we examine the arrangement of plant species across an oligotrophic bog/poor fen peatland complex in the North American boreal plain and the relationships of these species to their physical and chemical environment. A semi‐uniform spatial sampling approach was utilized to describe the species assemblages, pore‐water chemistry and physical condition of 100 plots throughout a single peatland complex. Regardless of sharing the same ground cover of Sphagnum mosses, the remaining species separated into four distinct assemblages, each with unique indicators. These species groups along with associated chemical and physical factors are organized into four ecosites: bog, margin (edge) and two poor fen ecosites. The plant assemblages of this peatland have a complex relationship with numerous gradients, both physical and chemical, including depth to water table, shade, pH, nutrient and base cation. Rather than being homogenous across the landscape, most environmental variables exhibit distinct spatial patterns and do so in relationship to the plant assemblages, forming spatially distinct ecosites across the complex. Base cation concentrations play a smaller role than previously thought in differentiating these ecosites, and in addition to shade and depth to water table, nitrogen in the form of dissolved organic nitrogen was highly related to the placement of these ecosites. Many significant chemical factors appear related to evaporative water loss within the peatland complex, and these chemical factors are used to differentiate the ecosites. However, the mediation of evaporative water loss is due largely to self‐generated responses of the plant assemblages related to shade through plant morphology and peat acrotelm development related to depth to water table. We conclude that plant species and associated environmental gradients act together to form spatially distinct ecosites. The distribution of these ecosites within this large, environmentally complex peatland is largely controlled by differing self‐generated responses along the hydrotopographical gradient of differential water loss.     

Erodibility of low‐compaction steep‐sloped reclaimed surface mine lands in the southern Appalachian region,USA

The use of loose spoils on steep slopes for surface coal mining reclamation sites has been promoted by the US Department of Interior, Office of Surface Mining for the establishment of native forest, as prescribed by the Forest Reclamation Approach (FRA). Although low‐compaction spoils improve tree survival and growth, erodibility on steep slopes was suspected to increase. This study quantified a combined K_C factor (combining the effects of the soil erodibility K factor and cover management C) for low compaction, steep‐sloped (>20°) reclaimed mine lands in the Appalachian region, USA. The combined K_C factor was used because standard Unit Plot conditions required to separate these factors, per Revised Universal Soil Loss Equation (RUSLE) experimental protocols, were not followed explicitly. Three active coal mining sites in the Appalachian region of East Tennessee, each containing four replicate field plots, were monitored for rainfall and sediment yields during a 14‐month period beginning June 2009. Average cumulative erosivity for the study sites during the monitoring period was measured as 5248.9 MJ·mm·ha^?1·h^?1. The K_C ranged between 0.001 and 0.05 t·ha·h·ha^?1·MJ^?1·mm^?1, with the highest values occurring immediately following reclamation site construction as rills developed (June – August 2009). The K_C for two study sites with about an 18–20 mm spoil D₈₄ were above 0.01 t·ha·h·ha^?1·MJ^?1·mm^?1 during rill development, and below 0.003 t·ha·h·ha^?1·MJ^?1·mm^?1 after August 2009 for the post‐rill development period. The K_C values for one site with a 40 mm spoil D₈₄ were never above 0.008 t·ha·h·ha^?1·MJ^?1·mm^?1 and also on average were lower, being more similar to the other two sites after the rill development period. Based on an initial K_C factor (K_e) measured during the first few storm events, the average C factor (C_e) was estimated as 0.58 for the rill development period and 0.13 for the post‐rill development period. It appears that larger size fractions of spoils influence K_C and C_e factors on low‐compaction steep slopes reclamation sites. Copyright © 2013 John Wiley & Sons, Ltd.