Fluvial sediment supply and pioneer woody seedlings as a control on bar‐surface topography

Plants influence river channel topography, but our understanding of the interaction among plants, flow, and sediment is limited, especially when sediment supply is variable. Using laboratory experiments in a recirculating flume with live seedlings in a mobile sand bed, we demonstrate how varying the balance between sediment supply and transport capacity shifts the relationship between plants and bar‐surface topography. Each experimental trial contrasted two sediment conditions, in which initially supply was maintained in equilibrium with transport via sediment recirculation, followed by sediment deficit, in which transport capacity exceeded supply, which was set to zero. For both sediment balances, the topographic response was sensitive to plant size, with larger plants inducing greater aggradation relative to a baseline condition. During sediment equilibrium, the positive relationship between plant size and topographic change also depended on species morphology (multi‐stemmed shrubs versus single‐stemmed plants). Plant morphology effects disappeared when the sediment balance shifted to a deficit, but the presence of plants had a greater impact on the magnitude of change compared to the topographic response under sediment equilibrium. Our results suggest that the interactions among sediment supply, plants, and topography may be strongest on rivers with a balance in sediment supply and transport capacity. Because of the large variability in fluvial sediment supply resulting from natural and anthropogenic influences, these interactions will differ spatially (e.g. longitudinally through a watershed) and at different temporal scales, from single flood events to longer time periods. Copyright © 2016 John Wiley & Sons, Ltd.     

Bar dynamics and bifurcation evolution in a modelled braided sand‐bed river

Morphodynamics in sand‐bed braided rivers are associated with simultaneous evolution of mid‐channel bars and channels on the braidplain. Bifurcations around mid‐channel bars are key elements that divide discharge and sediment. This, in turn, may control the evolution of connected branches, with effects propagating to both upstream and downstream bifurcations. Recent works on bifurcation stability and development hypothesize major roles of secondary flow and gradient advantage. However, this has not been tested for channel networks within a fully developed dynamic braided river. A reason for this is a lack of detailed measurements with sufficient temporal and spatial length, covering multiple bifurcations. Therefore we used a physics‐based numerical model to generate a dataset of bathymetry, flow and sediment transport of an 80 km river reach with self‐formed braid bars and bifurcations. The study shows that bar dissection due to local transverse water surface gradients is the dominant bifurcation initiation mechanism, although conversion of unit bars into compound bars dominates in the initial stage of a braided river. Several bifurcation closure mechanisms are equally important. Furthermore, the study showed that nodal point relations for bifurcations are unable to predict short‐term bifurcation evolution in a braided river. This is explained by occurrence of nonlinear processes and non‐uniformity within the branches, in particular migrating bars and larger‐scale backwater‐effects, which are not included in the nodal point relations. Planform morphology, on the other hand, has predictive capacity: bifurcation angle asymmetry and bar‐tail limb shape are indicators for near‐future bifurcation evolution. Remote sensing data has predictive value, for which we developed a conceptual model for interactions between bars, bifurcations and channels in the network. We conducted a preliminary test of the conceptual model on satellite images of the Brahmaputra. Copyright © 2015 John Wiley & Sons, Ltd.     

Spatio‐temporal sedimentation patterns in beaver ponds along the Chevral river,Ardennes, Belgium

With the recovery of the European beaver (Castor fiber) and their capacity to engineer fluvial landscapes, questions arise as to how they influence sediment transport, including the spatio‐temporal trends and patterns of sedimentation in beaver ponds. The Chevral river (Ardennes, Belgium) contains two beaver dam sequences, which appeared in 2004. Volumes of sediment deposited behind the dams were measured, and grain‐size distribution patterns were determined. Flow discharges and sediment fluxes were measured at the inflow and outflow of each dam sequence. Between 2004 and 2011, 1710.1 m³ of sediment was deposited behind the beaver dams, with an average sediment thickness of 25.1 cm. The thickness of the sediment layer was significantly (p < 0.001) related to the area of the beaver ponds. Along the stream, beaver pond sediment thickness displayed a sinusoidal deposition pattern, in which ponds with thick sediment layers were preceded by a series of ponds with thinner sediment layers. A downstream textural coarsening in the dam sequences was also observed, probably because of dam failures subsequent to surges. Differences in sediment flux between the inflow and outflow at the beaver pond sequence were related to the river hydrograph, with deposition taking place during the rising limbs and slight erosion during the falling limbs. The 7‐year‐old sequences have filtered 190.19 ton of sediment out of the Chevral river, which is of the same order of magnitude as the 374.4 ton measured in pond deposits, with the difference between the values corresponding to beaver excavations (60.24 ton), inflow from small tributaries, and runoff from the valley flanks. Hydrogeomorphic effects of C. fiber and Castor canadensis activity are similar in magnitude. The detailed analysis of sedimentation in beaver pond sequences confirms the potential of beavers to contribute to river and wetland restoration, and catchment management. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatial and temporal variations in bank erosion on sand‐bed streams in the seasonally wet tropics of northern Australia

Bank erosion rates and processes across a range of spatial scales are poorly understood in most environments, especially in the seasonally wet tropics of northern Australia where sediment yields are among global minima. A total of 177 erosion pins was installed at 45 sites on four sand‐bed streams (Tributaries North and Central, East Tributary and Ngarradj) in the Ngarradj catchment in the Alligator Rivers Region. Bank erosion was measured for up to 3·5 years (start of 1998/99 wet season to end of 2001/02 wet season) at three spatial scales, namely a discontinuous gully (0·6 km²) that was initiated by erosion of a grass swale between 1975 and 1981, a small continuous channel (2·5 km²) on an alluvial fan that was formed by incision of a formerly discontinuous channel between 1964 and 1978, and three medium‐sized, continuous channels (8·5–43·6 km²) with riparian vegetation. The bank erosion measurements during a period of average to above‐average rainfall established that substantial bank erosion occurred during the wet season on the two smaller channels by rapid lateral migration (Tributary Central) and by erosion of gully sidewalls due to a combination of within‐gully flows and overland flow plunging over the sidewalls (Tributary North). Minor bank erosion also occurred during the dry season by faunal activity, by desiccation and loss of cohesion of the sandy bank sediments and by dry flow processes. The larger channels with riparian vegetation (East Tributary and Ngarradj) did not generate significant amounts of sediment by bank erosion. Deposition (i.e. negative pin values) was locally significant at all scales. Bank profile form and channel planform exert a strong control on erosion rates during the wet season but not during the dry season. Copyright © 2006 Commonwealth Government of Australia.     

Hydraulic and geomorphic processes in an overbank flood along a meandering,gravel‐bed river: implications for chute formation

Hydraulic interactions between rivers and floodplains produce off‐channel chutes, the presence of which influences the routing of water and sediment and thus the planform evolution of meandering rivers. Detailed studies of the hydrologic exchanges between channels and floodplains are usually conducted in laboratory facilities, and studies documenting chute development are generally limited to qualitative observations. In this study, we use a reconstructed, gravel‐bedded, meandering river as a field laboratory for studying these mechanisms at a realistic scale. Using an integrated field and modeling approach, we quantified the flow exchanges between the river channel and its floodplain during an overbank flood, and identified locations where flow had the capacity to erode floodplain chutes. Hydraulic measurements and modeling indicated high rates of flow exchange between the channel and floodplain, with flow rapidly decelerating as water was decanted from the channel onto the floodplain due to the frictional drag provided by substrate and vegetation. Peak shear stresses were greatest downstream of the maxima in bend curvature, along the concave bank, where terrestrial LiDAR scans indicate initial floodplain chute formation. A second chute has developed across the convex bank of a meander bend, in a location where sediment accretion, point bar development and plant colonization have created divergent flow paths between the main channel and floodplain. In both cases, the off‐channel chutes are evolving slowly during infrequent floods due to the coarse nature of the floodplain, though rapid chute formation would be more likely in finer‐grained floodplains. The controls on chute formation at these locations include the flood magnitude, river curvature, floodplain gradient, erodibility of the floodplain sediment, and the flow resistance provided by riparian vegetation. Copyright © 2015 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.     

Comparison of fractional bed‐material load ­computation methods in sand‐bed channels

This paper presents a comparison of different fractional bed‐material load computation methods in sand‐bed channels. These methods include the direct computation by size fraction approach of Einstein, Laursen and Toffaleti; the bed material fraction (BMF) approach using equations of Engelund and Hansen, Ackers and White, Yang, and Karim's modified BMF method; and the transport capacity fraction (TCF) approach using the transport capacity distribution functions of Karim and Kennedy, Li, and Wu and Molinas. Over 150 sets of flume and field data in the sand size range containing a total of 1007 data points are used to evaluate and compare the accuracy of these methods. Statistical analysis and graphical comparison are utilized to demonstrate the performance and variations in different methods. Overall, the Einstein method underpredicts the transport rate for finer sizes and overpredicts for the coarser sizes, while the other methods overestimate the finer fractions and underestimate the coarser fractions. The Wu and Molinas method, which was developed to account for these deficiencies, is shown to significantly improve fractional bed‐material load computations. The Karim and Kennedy method is also found to be applicable. Copyright © 2000 John Wiley & Sons, Ltd.     

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.     

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.     

Morphological effects of different channel‐forming discharges in a gravel‐bed river

The study analyses the morphological response of a gravel‐bed river to discharges of different magnitude (from moderate events that occur several times a year to a 12‐year flood) and so defines the range of formative discharges for single morphological units (channels, bars, islands) and a range of magnitude of morphological activity from the threshold discharges for gravel transport and minor bar modification up to flows causing major morphological changes. The study was conducted on the Tagliamento River, a large gravel‐bed river in north‐eastern Italy, using two different methods, analysis of aerial photographs and field observation of painted gravel particles. The available photographs (five flights from August 1997 to November 2002) and the two commissioned flights (June 2006 and April 2007) do not define periods with a single flood event, but the intervals are short enough (11 to 22 months) to have a limited number of flood events in each case. The fieldwork, which involved cross‐section survey, grain‐size analysis and observation of painted sediments, complemented the aerial surveys by allowing analysis of channel response to single flood events. Substantial morphological changes (e.g. bank erosion of several tens of metres up to more than 100 m) associated with flood events with a recurrence interval between 1·1 year and 12 years have been documented. Multiple forming discharges were defined based on the activity of different morphological units. Discharges equal to 20–50% of the bankfull discharge are formative for the channels, whereas the bankfull discharge (1·1 year flood in this case of the Tagliamento River) is formative for low bars. Larger floods, but still relatively frequent (with a recurrence interval less than five years), are required for full gravel transport on high bars and significant morphological changes of islands. Copyright © 2009 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.     

Field studies on factors affecting very fine and ultra fine particulate organic matter deposition in low‐gradient sand‐bed streams

The knowledge on particle deposition in streams is mainly based on investigations in mountain streams. No data exist from low‐gradient sand‐bed streams that largely differ in the morphological and hydraulic factors proposed to affect deposition. To identify physical control on particle deposition in low‐gradient streams, we assessed deposition of very fine and ultra fine organic particulate matter in 18 sand‐bed stream reaches. We added particles derived from lake sediment and assessed the mean transport distance S_P and the deposition velocity v_dep. Additionally, reach hydraulics were estimated by injections of a conservative solute tracer (NaCl). Among the low‐gradient streams, particle deposition kinetics were variable but similar to deposition in mountain streams. S_P was solely related to the flow velocity. This relation was confirmed when comprising published data on deposition of fine organic particles. An association between particle deposition and transient storage factors was insignificant. We found significance of the transient storage to S_P only for repeated measures within a single reach, when flow velocity and benthic conditions were nearly constant. Measured v_dep/v_fall ratios were much larger than unity in most reaches. Evidence from this relation suggests that the vertical transport of very fine and ultra fine organic particulate matter through the water column was caused mainly by vertical mixing. Copyright © 2006 John Wiley & Sons, Ltd.     

Assessment of intermediate fine sediment storage in a braided river reach (southern French Prealps)

This paper presents a field investigation on river channel storage of fine sediments in an unglaciated braided river, the Bès River, located in a mountainous region in the southern French Prealps. Braided rivers transport a very large quantity of bedload and suspended sediment load because they are generally located in the vicinity of highly erosive hillslopes. Consequently, these rivers play an important role because they supply and control the sediment load of the entire downstream fluvial network. Field measurements and aerial photograph analyses were considered together to evaluate the variability of fine sediment quantity stored in a 2·5‐km‐long river reach. This study found very large quantities of fine sediment stored in this reach: 1100 t per unit depth (1 dm). Given that this reach accounts for 17% of the braided channel surface area of the river basin, the quantities of fine sediment stored in the river network were found to be approximately 80% of the mean annual suspended sediment yields (SSYs) (66 200 t year^?1), comparable to the SSYs at the flood event scale: from 1000 t to 12 000 t depending on the flood event magnitude. These results could explain the clockwise hysteretic relationships between suspended sediment concentrations and discharges for 80% of floods. This pattern is associated with the rapid availability of the fine sediments stored in the river channel. This study shows the need to focus on not only the mechanisms of fine sediment production from hillslope erosion but also the spatiotemporal dynamics of fine sediment transfer in braided rivers. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Spatial and temporal factors controlling short‐term sedimentation in a salt and freshwater tidal marsh,Scheldt estuary,Belgium, SW Netherlands

During a one‐year period temporal and spatial variations in suspended sediment concentration (SSC) and deposition were studied on a salt and freshwater tidal marsh in the Scheldt estuary (Belgium, SW Netherlands) using automatic water sampling stations and sediment traps. Temporal variations were found to be controlled by tidal inundation. The initial SSC, measured above the marsh surface at the beginning of inundation events, increases linearly with inundation height at high tide. In accordance with this an exponential relationship is observed between inundation time and sedimentation rates, measured over 25 spring–neap cycles. In addition both SSC and sedimentation rates are higher during winter than during summer for the same inundation height or time. Although spatial differences in vegetation characteristics are large between and within the studied salt and freshwater marsh, they do not affect the spatial sedimentation pattern. Sedimentation rates however strongly decrease with increasing (1) surface elevation, (2) distance from the nearest creek or marsh edge and (3) distance from the marsh edge measured along the nearest creek. Based on these three morphometric parameters, the spatio‐temporal sedimentation pattern can be modelled very well using a single multiple regression model for both the salt and freshwater marsh. A method is presented to compute two‐dimensional sedimentation patterns, based on spatial implementation of this regression model. Copyright © 2003 John Wiley & Sons, Ltd.     

The influence of flow discharge variations on the morphodynamics of a diffluence–confluence unit on a large river

Bifurcations are key geomorphological nodes in anabranching and braided fluvial channels, controlling local bed morphology, the routing of sediment and water, and ultimately defining the stability of their associated diffluence–confluence unit. Recently, numerical modelling of bifurcations has focused on the relationship between flow conditions and the partitioning of sediment between the bifurcate channels. Herein, we report on field observations spanning September 2013 to July 2014 of the three‐dimensional flow structure, bed morphological change and partitioning of both flow discharge and suspended sediment through a large diffluence–confluence unit on the Mekong River, Cambodia, across a range of flow stages (from 13 500 to 27 000 m³ s^?1). Analysis of discharge and sediment load throughout the diffluence–confluence unit reveals that during the highest flows (Q = 27 000 m³ s^?1), the downstream island complex is a net sink of sediment (losing 2600 ± 2000 kg s^?1 between the diffluence and confluence), whereas during the rising limb (Q = 19 500 m³ s^?1) and falling limb flows (Q = 13 500 m³ s^?1) the sediment balance is in quasi‐equilibrium. We show that the discharge asymmetry of the bifurcation varies with discharge and highlight that the influence of upstream curvature‐induced water surface slope and bed morphological change may be first‐order controls on bifurcation configuration. Comparison of our field data to existing bifurcation stability diagrams reveals that during lower (rising and falling limb) flow the bifurcation may be classified as unstable, yet transitions to a stable condition at high flows. However, over the long term (1959–2013) aerial imagery reveals the diffluence–confluence unit to be fairly stable. We propose, therefore, that the long‐term stability of the bifurcation, as well as the larger channel planform and morphology of the diffluence–confluence unit, may be controlled by the dominant sediment transport regime of the system. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Bedload yields for sand‐bed streams in the Ngarradj Creek catchment,Northern Territory,Australia

Bedload yields were calculated by 39 methods at the East Tributary gauge, nine methods at Upper Swift Creek gauge and 11 methods at Swift Creek gauge in the Ngarradj Creek catchment in northern Australia. These methods involved combining various significant bedload rating curves determined for a measured bedload data set for a 4‐year period with either the hourly or daily hydrographs or flow duration curves for the same period, 1 September 1998 to 31 August 2002. Bedload ratings were both statistically significant (ρ ≤ 0.05) and explained at least 60% of the variance in bedload flux. Bias corrections were used with all methods based on log₁₀‐transformed ratings. Estimated mean annual bedload yields varied by three orders of magnitude at the East Tributary gauge and by two orders of magnitude at Upper Swift Creek and Swift Creek gauges. Hourly discharges usually produced higher estimated yields than daily discharges. The bedload rating‐flow duration curve technique overestimates yields and bias correction methods always produce even higher yields. Ratings using both immersed bedload weight and adjusted immersed bedload weight always under‐predicted yields because they contain an implicit threshold of motion condition that is at least four times greater than that predicted by Bagnold's threshold equation. Such a result questions the applicability of Bagnold's threshold equation to the Ngarradj Creek catchment. The best estimates of mean annual bedload yield at East Tributary, Upper Swift Creek and Swift Creek gauges are 600 ± 170 (SE), 1065 ± 150 and 1795 ± 270 t/year, respectively. © 2015 Commonwealth of Australia. Hydrological Processes © 2015 John Wiley & Sons Ltd.