Channel response to an extreme flood and sediment pulse in a mixed bedrock and gravel‐bed river

We exploit a natural experiment caused by an extreme flood (~500 year recurrence interval) and sediment pulse derived from more than 2500 concurrent landslides to explore the influence of valley‐scale geomorphic controls on sediment slug evolution and the impact of sediment pulse passage and slug deposition and dispersion on channel stability and channel form. Sediment slug movement is a crucial process that shapes gravel‐bed rivers and alluvial valleys and is an important mechanism of downstream bed material transport. Further, increased bed material transport rates during slug deposition can trigger channel responses including increases in lateral mobility, channel width, and alluvial bar dominance. Pre‐ and post‐flood LiDAR and aerial photographs bracketing the 2007 flood on the Chehalis River in south‐western Washington State, USA, document the channel response with high spatial and temporal definition. The sediment slug behaved as a Gilbert Wave, with both channel aggradation and sequestration of large volumes of material in floodplains of headwaters' reaches and reaches where confined valleys enter into broad alluvial valleys. Differences between the valley form of two separate sub‐basins impacted by the pulse highlight the important role channel and channel‐floodplain connectivity play in governing downstream movement of sediment slug material. Finally, channel response to the extreme flood and sediment pulse illustrate the connection between bed material transport and channel form. Specifically, the channel widened, lateral channel mobility increased, and the proportion of the active channel covered by bars increased in all reaches in the study area. The response scaled tightly with the relative amount of bed material sediment transport through individual reaches, indicating that the amount of morphological change caused by the flood was conditioned by the simultaneous introduction of a sediment pulse to the channel network. Copyright © 2015 John Wiley & Sons, Ltd.     

The role of channel morphology on the mobility and dispersion of bed sediment in a small gravel‐bed stream

This study uses a unique 10‐year tracer dataset from a small gravel‐bed stream to examine bed mobility and sediment dispersion over long timescales and at a range of spatial scales. Seasonal tracer data that captured multiple mobilizing events was examined, while the effects of morphology on bed mobility and sediment dispersion were captured at three spatial scales: within morphological units (unit scale), between morphological units (reach scale) and between reaches with different channel morphologies (channel scale). This was achieved by analyzing both reach‐average mobility and travel distance data, as well as the development of ‘mobility maps’ that capture the spatial variability in tracer mobility within the channel. The tracer data suggest that sediment transport in East Creek remains near critical the majority of the time, with only rare large events resulting in high mobility rates and grain travel distances large enough to move sediment past dominant bedforms. While a variable capturing both the magnitude and frequency of flow events within a season yielded a better predictor to sediment mobility and dispersion than peak discharge alone, the distribution of events of different magnitude within the season played a large role in determining tracer mobility rates and travel distances. The effects of morphology differed depending on the analysis scale, demonstrating the importance of scale, and therefore study design, when examining the effect of morphology on sediment transport. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluating fine sediment mobilization and storage in a gravel‐bed river using controlled reservoir releases

Two controlled flow events were generated by releasing water from a reservoir into the Olewiger Bach, located near Trier, Germany. This controlled release of near bank‐full flows allowed an investigation of the fine sediment (<63 µm) mobilized from channel storage. Both a winter (November) and a summer (June) release event were generated, each having very different antecedent flow conditions. The characteristics of the release hydrographs and the associated sediment transport indicated a reverse hysteresis with more mass, but smaller grain sizes, moving on the falling limb. Fine sediment stored to a depth of 10 cm in the gravels decreased following the release events, indicating the dynamic nature and importance of channel‐stored sediments as source materials during high flow events. Sediment traps, filled with clean natural gravel, were buried in riffles before the release of the reservoir water and the total mass of fine sediment collected by the traps was measured following the events. Twice the mass of fine sediment was retained by the gravel traps compared with the natural gravels, which may be due to their altered porosity. Although the amount of fine sediment collected by the traps was not significantly related to measures of gravel structure, it was found to be significantly correlated to measures of local flow velocity and Froude number. A portion of the traps were fitted with lids to restrict surface exchange of water and sediment. These collected the highest amounts of event‐mobilized sediments, indicating that inter‐gravel lateral flows, not just surface infiltration of sediments, are important in replenishing and redistributing the channel‐stored fines. These findings regarding the magnitude and direction of fine sediment movement in gravel beds are significant in both a geomorphic and a biological context. Copyright © 2006 John Wiley & Sons, Ltd.     

The spatial and temporal patterns of aggradation in a temperate,upland, gravel‐bed river

Intensive field monitoring of a reach of upland gravel‐bed river illustrates the temporal and spatial variability of in‐channel sedimentation. Over the six‐year monitoring period, the mean bed level in the channel has risen by 0·17 m with a maximum bed level rise of 0·5 m noted at one location over a five month winter period. These rapid levels of aggradation have a profound impact on the number and duration of overbank flows with flood frequency increasing on average 2·6 times and overbank flow time increasing by 12·8 hours. This work raises the profile of coarse sediment transfer in the design and operation of river management, specifically engineering schemes. It emphasizes the need for the implementation of strategic monitoring programmes before engineering work occurs to identify zones where aggradation is likely to be problematic. Exploration of the sediment supply and transfer system can explain patterns of channel sedimentation. The complex spatial, seasonal and annual variability in sediment supply and transfer raise uncertainties into the system's response to potential changes in climate and land‐use. Thus, there is a demand for schemes that monitor coarse sediment transfer and channel response. Copyright © 2009 John Wiley & Sons, Ltd.     

Effect of changes in fine‐grained matrix on bedload sediment transport in a gravel‐bed river

While clay and silt matrices of gravel‐bed rivers have received attention from ecologists concerned variously with the deteriorating environments of benthic and hyporheic organisms, their impact on sediment entrainment and transport has been explored less. A recent increase of such a matrix in the bed of Nahal Eshtemoa, an ephemeral river of the northern Negev, has more than doubled the boundary shear stress needed to initiate bedload, from 7 N m^‐2 (τ^* = 0.027) during the flash floods of 1991–2001 to 15 N m^‐2 (τ^* = 0.059) during those of 2008–2009. The relation between bedload flux and boundary shear stress continues to be well‐defined, but it is displaced. The matrix now contains a significant amount of silt and clay size material. The reasons for the increased entrainment threshold of bedload are explored. Large‐scale laser scanning of the dry bed reveals a reduction in grain‐scale morphological roughness, while artificial in situ tests of matrix integrity indicate considerable cohesion. The implications for adopting bed material sampling strategies that account for matrix development are assessed. Copyright © 2012 John Wiley & Sons, Ltd.     

Comparison of hydraulic geometry between sand‐ and gravel‐bed rivers in relation to channel pattern discrimination

A comparison has been made between the hydraulic geometry of sand‐ and gravel‐bed rivers, based on data from alluvial rivers around the world. The results indicate a signi?cant difference in hydraulic geometry among sand‐ and gravel‐bed rivers with different channel patterns. On this basis, some diagrams for discrimination of meandering and braided channel patterns have been established. The relationships between channel width and water discharge, between channel depth and water discharge, between width–depth ratio and water discharge and between channel slope and water discharge can all be used for channel pattern discrimination. The relationship between channel width and channel depth can also be used for channel pattern discrimination. However, the accuracy of these relationships for channel pattern discrimination varies, and the depth–discharge relationship is a better discriminator of pattern type than the classic slope–discharge function. The cause for this difference has been explained qualitatively. To predict the development of channel patterns under different natural conditions, the pattern discriminator should be searched on the basis of independent or at least semi‐independent variables. The relationship between stream power and bed material grain size can be used to discriminate channel patterns, which shows a better result than the discriminator using the slope–discharge relationship. Copyright © 2004 John Wiley & Sons, Ltd.     

Insights into bedload transport processes of a large regulated gravel‐bed river

A comprehensive monitoring programme focusing on bedload transport behaviour was conducted at a large gravel‐bed river. Innovative monitoring strategies were developed during five years of preconstruction observations accompanying a restoration project. A bedload basket sampler was used to perform 55 cross‐sectional measurements, which cover the entire water discharge spectrum from a 200‐year flood event in 2013 to a rare low flow event. The monitoring activities provide essential knowledge regarding bedload transport processes in large rivers. We have identified the initiation of motion under low flow conditions and a decrease in the rate of bedload discharge with increasing water discharge around bankfull conditions. Bedload flux strongly increases again during high flood events when the entire inundation area is flooded. No bedload hysteresis was observed. The effective discharge for bedload transport was determined to be near mean flow conditions, which is therefore at a lower flow discharge than expected. A numerical sediment transport model was able to reproduce the measured sediment transport patterns. The unique dataset enables the characterisation of bedload transport patterns in a large and regulated gravel‐bed river, evaluation of modern river engineering measures on the Danube, and, as a pilot project has recently been under construction, is able to address ongoing river bed incision, unsatisfactory ecological conditions for the adjacent national park and insufficient water depths for inland navigation. Copyright © 2017 John Wiley & Sons, Ltd.     

Modeling the effects of pulsed versus chronic sand inputs on salmonid spawning habitat in a low‐gradient gravel‐bed river

It is widely recognized that high supplies of fine sediment, largely sand, can negatively impact the aquatic habitat quality of gravel‐bed rivers, but effects of the style of input (chronic vs. pulsed) have not been examined quantitatively. We hypothesize that a continuous (i.e. chronic) supply of sand will be more detrimental to the quality of aquatic habitat than an instantaneous sand pulse equal to the integrated volume of the chronic supply. We investigate this issue by applying a two‐dimensional numerical model to a 1 km long reach of prime salmonid spawning habitat in central Idaho. Results show that in both supply scenarios, sand moves through the study reach as bed load, and that both the movement and depth of sand on the streambed mirrors the hydrograph of this snowmelt‐dominated river. Predictions indicate greater and more persistent mortality of salmonid embryos under chronic supplies than pulse inputs, supporting our hypothesis. However, predicted mortality varies both with salmonid species and location of spawning. We found that the greatest impacts occur closer to the location of the sand input under both supply scenarios. Results also suggest that reach‐scale morphology may modulate the impact of sand loads, and that under conditions of high sand loading climate‐related increases in flow magnitude could increase embryo mortality through sand deposition, rather than streambed scour. Copyright © 2013 John Wiley & Sons, Ltd.     

Life of a bifurcation in a gravel‐bed braided river

Channel bifurcation is a key element in braided rivers, determining the water and sediment distribution and hence controlling the morphological evolution. Recent theoretical and experimental findings, as well as field observations, showed that bifurcations in gravel‐bed braided rivers are often asymmetrical and highly unstable. In this paper field data are presented on a bifurcation in the Tagliamento River, northeast Italy. The planform configuration of the bifurcation and its temporal evolution was monitored by an automatic digital camera during a series of seven floods with different magnitudes. This remote sensing technique allowed a high temporal resolution (pictures were acquired every hour) that was proved to be essential in a highly dynamic system as the one considered here. Digitized maps of the channels provided information on the location of the bifurcation, the width of the anabranches, the angle between them, along with the occurrence and migration of sediment bars. Data were acquired at two different water levels, giving the possibility to compare low and high flow conditions. The monitored bifurcation is largely unstable and shows sudden changes in the water distribution, mainly driven by the bar migrating in the upstream channel and entering the distributaries. A relationship between width asymmetry and flood magnitude was observed, confirming previous analyses. Moreover, recent theoretical findings were applied, in order to test the possibility to estimate general trends in bifurcation evolution. The analysis pointed out the relevance of a correct assessment of the characteristic temporal scales, as the bifurcation evolves on a timescale similar to that of bar migration and flood duration. Understanding the interactions between these processes is therefore crucial in order to increase the ability to model and predict the morphological evolution of a braided network. Copyright © 2012 John Wiley & Sons, Ltd.     

A two‐dimensional morphodynamic model of gravel‐bed river with floodplain vegetation

Artificially straight river channels tend to be unstable, and ultimately develop into river meanders through bank erosion and point‐bar deposition. In this paper account is taken of the effects of riparian and floodplain vegetation on bank strength, floodplain flow resistance, shear stress partitioning, and bedload transport. This is incorporated into an existing 2D hydrodynamic‐morphological model. By applying the new model to an initially straight and single‐threaded channel, the way that its planform and cross‐sectional geometry evolve for different hydraulic and floodplain vegetation conditions is demonstrated. The results show the formation and upstream migration of gravel bars, confluence scouring and the development of meandering and braiding channel patterns. In cases where the channel becomes unstable, the instability grows out of bar formation. The resulting braiding patterns are similar to analytical results. The formation of a transition configuration requires a strong influence from vegetation. Copyright © 2010 John Wiley & Sons, Ltd.     

Particle path length distributions in meandering gravel‐bed streams: results from physical models

In gravel‐bed rivers with well‐de?ned pool–bar morphology, the path length of transported bed particles must be, at least during ‘channel‐forming’ ?ows, equal to the length scale of the morphology. This is the basis for some methods for estimating bed material transport rates. However, previous data, especially from ?eld tests, are often strongly positively skewed with mean much shorter than the pool–bar spacing. One possible explanation is that positively skewed distributions occur only in channels lacking distinct pool–bar topography or only at lower discharges in pool–bar channels. A series of ?ume experiments using ?uorescent tracers was used to measure path length distributions in low‐sinuosity meandering channels to assess the relation with channel morphology and ?ow conditions. At channel‐forming ?ows, 55 to 75 per cent of the tracer grains were deposited on the ?rst point bar downstream of the point of tracer input, with 15 per cent passing beyond the ?rst bar. Path length distributions are symmetrical with mean equal to the pool–bar spacing and can be described with a Cauchy distribution. In some cases there was a secondary mode close to the point of tracer introduction; this bimodal distribution ?ts a combined gamma–Cauchy distribution. Only when discharge was reduced below the channel‐forming ?ow were frequency distributions unimodal and positively skewed with no relation to the pool–bar spacing. Thus, path length distributions become more symmetrical, and mean path length increases to coincide with pool–bar spacing, as ?ow approaches channel‐forming conditions. This is a substantial modi?cation of existing models of particle transfer in gravel‐bed rivers. Copyright © 2003 John Wiley & Sons, Ltd.     

Transport and storage of bed material in a gravel‐bed channel during episodes of aggradation and degradation: a field and flume study

The dynamics of sediment transport capacity in gravel‐bed rivers is critical to understanding the formation and preservation of fluvial landforms and formulating sediment‐routing models in drainage systems. We examine transport‐storage relations during cycles of aggradation and degradation by augmenting observations of three events of channel aggradation and degradation in Cuneo Creek, a steep (3%) gravel‐bed channel in northern California, with measurements from a series of flume runs modeling those events. An armored, single‐thread channel was formed before feed rates were increased in each aggradation run. Output rates increased as the channel became finer and later widened, steepened, and braided. After feed rates were cut, output rates remained high or increased in early stages of degradation as the incising channel remained fine‐grained, and later decreased as armoring intensified. If equilibrium was not reached before sediment feed rate was cut, then a rapid transition from a braided channel to a single‐thread channel caused output rates for a given storage volume to be higher during degradation than during aggradation. Variations in channel morphology, and surface bed texture during runs that modeled the three cycles of aggradation and degradation were similar to those observed in Cuneo Creek and provide confidence in interpretations of the history of change: Cuneo Creek aggraded rapidly as it widened, shallowed, and braided, then degraded rapidly before armoring stabilized the channel. Such morphology‐driven changes in transport capacity may explain the formation of flood terraces in proximal channels. Transport‐storage relations can be expected to vary between aggradation and degradation and be influenced by channel conditions at the onset of changes in sediment supply. Published in 2011. This article is a US Government work and is in the public domain in the USA.     

Persistence of the surface texture of a gravel‐bed river during a large flood

We present herein clear field evidence for the persistence of a coarse surface layer in a gravel‐bed river during flows capable of transporting all grain sizes present on the channel bed. Detailed field measurements of channel topography and bed surface grain size were made in a gravel‐bed reach of the Colorado River prior to a flood in 2003. Runoff produced during the 2003 snowmelt was far above average, resulting in a sustained period of high flow with a peak discharge of 27 m³/s (170% of normal peak flow); all available grain sizes within the study reach were mobilized in this period of time. During the 2003 peak flow, the river avulsed immediately upstream of the study reach, thereby abandoning approximately one half kilometer of the former channel. The abandonment was rapid (probably within a few hours), leaving the bed texture essentially frozen in place at the peak of the flood. All locations sampled prior to the flood were resampled following the stream abandonment. In response to the high flow, the surface median grain size (D_50s) coarsened slightly in the outer part of the bend while remaining nearly constant along the inner part of the bend, resulting in an overall increase from 18 to 21 mm for the study reach. Thus, the coarse bed surface texture persisted despite shear stresses throughout the bend that were well above the critical entrainment value. This may be explained because the response of the bed texture to increases in flow strength depends primarily upon the continued availability of the various grain size percentiles in the supply, which in this case was essentially unlimited for all sizes present in the channel. Copyright © 2007 John Wiley & Sons, Ltd.     

A coupled sediment routing and lateral migration model for gravel‐bed rivers

We describe additions made to a multi‐size sediment routing model enabling it to simulate width adjustment simultaneously alongside bed aggradation/incision and fining/coarsening. The model is intended for use in single thread gravel‐bed rivers over annual to decadal timescales and for reach lengths of 1–10 km. It uses a split‐channel approach with separate calculations of flow and sediment transport in the left and right sides of the channel. Bank erosion is treated as a function of excess shear stress with bank accretion occurring when shear stress falls below a second, low, threshold. A curvature function redistributes shear stress to either side of the channel. We illustrate the model through applications to a 5·6‐km reach of the upper River Wharfe in northern England. The sediment routing component with default parameter values gives excellent agreement with field data on downstream fining and down‐reach reduction in bedload flux, and the width‐adjustment components with approximate calibration to match maximum observed rates of bank shifting give plausible patterns of local change. The approach may be useful for exploring interactions between sediment delivery, river management and channel change in upland settings. Copyright © 2011 John Wiley & Sons, Ltd.     

Fine‐grained bed sediment storage within the main channel systems of the Frome and Piddle catchments,Dorset, UK

Lowland permeable catchments in the UK are particularly prone to sedimentation problems, on account of the increased fine sediment loadings generated by recent land‐use change and their stable seasonal hydrological regimes, which are frequently depleted by groundwater abstraction. Fine‐grained sediment storage on the bed of the main channel systems of the Frome (437 km²) and Piddle (183 km²) catchments, Dorset, UK, has been examined at 29 sites using a sediment remobilization technique. Measurements encompassed the period February 2003–July 2004. At individual sites in the Frome, average values ranged between 410 and 2630 g m^?2, with an overall mean of 918 g m^?2. In the Piddle, the average values for individual sites varied between 260 and 4340 g m^?2, with an overall mean of 1580 g m^?2. Temporal variations in fine bed sediment storage at each site were appreciable, with the coefficients of variation ranging between 43 and 155% in the Frome and between 33 and 160% in the Piddle. Average reach‐scale specific bed sediment storage increased markedly downstream along each main stem from 2 to 29 t km^?1 (Frome) and from 4 to 19 t km^?1 (Piddle). Total fine sediment storage on the channel bed of the Frome varied between 479 t (5 t km^?1) and 1694 t (17 t km^?1), with a mean of 795 t (7 t km^?1), compared with between 371 t (5 t km^?1) and 1238 t (14 t km^?1) with a mean of 730 t (9 t km^?1) in the Piddle. During the study period, fine bed sediment storage was typically equivalent to 18% (Frome) and 57% (Piddle) of the mean annual suspended sediment flux at the study catchment outlets. Copyright © 2006 John Wiley & Sons, Ltd.     

Benefits of design modifications to the Birkbeck bedload sampler illustrated by flash‐floods in an ephemeral gravel‐bed channel

The benefits of three simple modifications to the design of a Birkbeck bedload slot‐sampling system that has been continuously operating in Nahal Eshtemoa, Israel, since the early 1990s are demonstrated. The modifications include the deployment of a removable slot cover which delays the accumulation of sediment, so allowing sampling at late stages of a flood and, in conjunction with other samplers, extending the period of sampling during a flood wave; inclusion of a slot the size of which is adjustable so that that the probability of sampling the largest clast sizes in transit as bedload can be increased post‐installation, once knowledge is gained about the bedload grain‐size distribution; and a sampler side‐wall door that allows stratification and textural changes within the accumulated bedload to be identified, so promoting intelligent sampling of the deposit for grain‐size determination. Results from seven flash‐floods are presented and discussed, with recommendations for bedload monitoring, particularly in rivers where sediment flux is high and dynamic sediment records are inevitably short because of instrumental limitations. Copyright © 2006 John Wiley & Sons, Ltd.     

Response of river channel patterns to the spatially varying suspended sediment concentrations

River channel pattern transformation is dealt with in a broad background of suspended sediment concentration, varying from low, medium, high to hyperconcentration. Based on data from about 100 alluvial rivers in China, suspended sediment transport rate has been plotted against mean annual water discharge, showing that all points can be divided into four belts by three straight lines, as stable braided pattern, meandering pattern with ordinary sediment concentrations, wandering braided pattern, and meandering pattern with hyperconcentration of sediment. This picture of channel pattern transformation can be well explained by the law of the water flow's energy expenditure varying with its sediment concentration. The energy expenditure increases with sediment concentration, reaching a maximum, then declines. Rivers falling in different ranges of sediment concentrations adjust their own energy expenditure in different manners, leading to occurrence of different channel patterns. Project supported by the National Natural Science Foundation of China (Grant No. 49671011).     

Flume‐ and field‐based evaluation of a time‐integrated suspended sediment sampler for the analysis of sediment properties

Suspended sediment has been identified as a vector for nutrient and contaminant transport in the fluvial environment. A time‐integrated sampler (the Phillips sampler), which emerged over a decade ago as a cost‐effective tool for in situ suspended sediment collection, is increasingly being used to collect samples for the analysis of sediment properties such as particle size composition, and nutrient and contaminant concentrations. This study evaluates the sampler under both flume and field conditions for efficiency in the mass and grain size of the suspended sediment collected. The sampler was tested in a flume using both kaolinite and sediment samples (sieved to < 180 µm) collected from the Quesnel River, British Columbia, Canada. In the kaolinite trails, the sampler preferentially collected coarser grain sizes compared to the original sediment, probably due to finer sediment remaining in suspension and therefore passing through the sampler, and also possibly due to flocculation of the kaolinite upon introduction to the flume. Conversely, the sampler collected river sediment that was finer than the original sediment, probably due to some settling of coarser sediment observed at the bottom of the flume. Once allowance was made for these operational issues associated with the flume, maximum sediment mass efficiency for kaolinite and river sediment was 43% and 87%, respectively. Sediment collected by the time‐integrated sampler during field deployment and adjacent channel bed sediment were also compared. The sampler collected sediment with a representative grain size distribution. However, there were differences in the geochemical (arsenic and selenium) concentrations of channel bed sediment and sediment collected by the Phillips sampler which may be a function of differences in the behavior of geochemical elements associated with the two types of sediment. This work suggests that further research is needed to evaluate the role of the Phillips sampler in collecting sediment for contaminant and nutrient analysis. Copyright © 2014 John Wiley & Sons, Ltd.     

Challenges to the representation of suspended sediment transfer using a depth‐averaged flux

The sediment saturation recovery process (i.e. the adaptation of suspended sediment concentration [SSC] to local forcing) is the main feature of the non‐equilibrium suspended sediment transport (SST) frequently occurring in fluvial, estuarine and coastal waters. In order to quantitatively describe this phenomenon, a series solution is analytically derived, including the evolution of both vertical SSC profile and near‐bed sediment flux (NBSF), and is verified by net erosion and net deposition experiments, respectively. The results suggest that the sediment saturation recovery process involves vertically varying fluxes that are not represented correctly by depth‐averaging. Consequently, a vertical two‐dimensional (2D) combined scheme is established and applied respectively in to a dredged trench and to a sand wave feature to demonstrate this argument. By analyzing the variations of the calculated depth‐averaged SSC and NBSF we reveal that the equilibrium state presented by the sediment carrying capacity (SCC) form of the NBSF, which is usually applied in depth‐integrated SST models, lags behind the actual dynamic bed equilibrium state. Moreover, the key factor α, the so‐called saturation recovery coefficient within this form, is not only a function of local Rouse number but also is influenced by the local SSC profile. Finally, a three‐dimensional (3D) non‐orthogonal curvilinear body‐fitted SST model is developed and validated in the Yangtze estuary, China, combined with the in situ hourly hydrographic data from August 14–15, 2007 during spring tide in the wet season. Model results confirm that the vertically varying sediment saturation recovery process, the discrepancies between the actual and SCC form of NBSF and non‐constant value of α are significant in actual real geomorphic cases. The quantitative morphological change resulting from variations in environmental conditions may not be correctly represented by uncorrected depth‐integrated SST models if they do not treat the effects of vertical motion on the sediment saturation recovery process. Copyright © 2016 John Wiley & Sons, Ltd.