Suspended sediment transport along an idealised tidal embayment: settling lag, residual transport and the interpretation of tidal signals

We present semi-analytical solutions for suspended sediment concentration (SSC) and residual sediment transport in a simple mathematical model of a short tidal embayment. These solutions allow us to investigate in some detail the characteristic tidal and semi-tidal variation of SSC and the processes leading to residual sediment transport, including settling and scour lags, the roles of ‘local’ and ‘advective’ contributions, and the presence of internally or externally generated overtides. By interpreting the transport mechanisms in terms of the classic conceptual models of settling lag we clarify how these models may be expressed in mathematical terms. Our results suggest that settling lag is usually a more important process than scour lag, and that a local model which neglects advection may predict the direction of net sediment transport incorrectly. Finally, we discuss our results in the context of other transport processes and morphodynamic feedback.     

Data assimilation (4D-VAR) to forecast flood in shallow-waters with sediment erosion

In this paper, the four-dimensional variational data assimilation technique (4D-VAR) is presented as a tool to forecast floods. Our study is limited to purely hydrological flows and supposes that the weather, here a big rain, has been already forecasted by meteorological services. The technique consists in minimizing, in the sense of Lagrange, the cost function: a measure of the difference between calculated data and available observations, here the water level. This is done under constraints that are the equations of the physical model. In our case, we modified the shallow-water equations to include a simplified sediment transport model. The steepest descent algorithm is then used to find the minimum. This is made possible because we can compute analytically the gradient of the cost function by using the adjoint equations of the model. As an application of the 4D-VAR technique, the overflowing of the Chicoutimi River at the Chute-Garneau dam, during the 1996 flood, is investigated. It is found that the 4D-VAR method reduces the error in the water height forecast even when the erosion model is not activated. In terms of Lyapunov exponents, we estimate the predictability horizon of such an event to be about half-an-hour after a big rain. However, this limit of predictability can be increased by using more observations or by using a finer computational grid.     

Signatures in flowing fluid electric conductivity logs

Flowing fluid electric conductivity logging provides a means to determine hydrologic properties of fractures, fracture zones, or other permeable layers intersecting a borehole in saturated rock. The method involves analyzing the time-evolution of fluid electric conductivity (FEC) logs obtained while the well is being pumped and yields information on the location, hydraulic transmissivity, and salinity of permeable layers. The original analysis method was restricted to the case in which flows from the permeable layers or fractures were directed into the borehole (inflow). Recently, the method was adapted to permit treatment of both inflow and outflow, including analysis of natural regional flow in the permeable layer. A numerical model simulates flow and transport in the wellbore during flowing FEC logging, and fracture properties are determined by optimizing the match between simulation results and observed FEC logs. This can be a laborious trial-and-error procedure, especially when both inflow and outflow points are present. Improved analyses methods are needed. One possible tactic would be to develop an automated inverse method, but this paper takes a more elementary approach and focuses on identifying the signatures that various inflow and outflow features create in flowing FEC logs. The physical insight obtained provides a basis for more efficient analysis of these logs, both for the present trial and error approach and for a potential future automated inverse approach. Inflow points produce distinctive signatures in the FEC logs themselves, enabling the determination of location, inflow rate, and ion concentration. Identifying outflow locations and flow rates typically requires a more complicated integral method, which is also presented in this paper.     

Transport rate of calcareous sand in unidirectional flow

This paper examines the transport of calcareous sand in unidirectional flow and its prediction through existing sediment transport models. A flume experiment of four sand samples collected on Oahu, Hawaii, provides 29 sets of sediment transport data in the bed-form and suspended transport stages. The measured transport data are compared with direct predictions from four energy-based transport models developed for siliceous particles. Corrections for the grain-size, fall velocity, and critical velocity of calcareous sand based on recent research are applied to the models and the results are compared with the direct calculations and measured data. The comparison illustrates the important role particle shape plays in the transport of calcareous sand. All four sediment transport models give consistent predictions and good agreement with the majority of the measured data. Two of the models respond positively to the corrections in both the bed-form and suspended transport stages indicating that such an approach may provide an interim solution for the transport of calcareous sand.     

Influence of vegetation on aeolian sand transport rate from a backshore to a foredune at Hasaki, Japan

The influence of vegetation on aeolian sediment transport rate in the region from a backshore to a foredune was investigated at the Hasaki Coast in Japan, where an onshore wind was predominant and the creeping beach grasses Carex kobomugi and Calystegia soldanella were major species. The comparison of cross-shore distributions of the cross-shore component of aeolian sand transport rate with and without vegetation, which were estimated on the basis of the beach profile changes and a mass conservation equation, showed that the creeping grasses influenced the aeolian sand transport rate. The landward aeolian sand transport rate rapidly decreased landward from the seaward limit of vegetation when the grasses grew. The aeolian sand transport rate reduced by 95% with a vegetation cover of 28%. On the other hand, when the grasses were absent, the landward aeolian sand transport rate did not decrease near the seaward vegetation limit, but near the foot of the foredune.     

Modeling the influence of settling velocity on cohesive sediment transport in Tanshui River estuary

Settling velocities of suspended cohesive sediment in estuaries vary over a range of several orders in magnitude. Variations in the suspended sediment concentration are often considered as the principal cause. Turbulence and the suspended sediment concentration, as well as other factors such as salinity, dissolved organic substances, flocculation ability, and the rate of floc growth affect setting velocities. A laterally–averaged finite difference model for hydrodynamics and cohesive sediment transport is developed and applied in the Tanshui River estuary, Taiwan. The model has been calibrated and verified with water surface elevation, longitudinal velocity, salinity, and cohesive sediment measured. The overall performance of the model is in qualitative agreement with the available data. The model is used to investigate the influence of settling velocity on cohesive sediment transport dynamics. The simulation indicates that the turbidity maximum zone is near Kuan–Du. When settling velocities increase the surface cohesive sediment concentration at Kuan–Du station trends to decrease and bottom cohesive sediment concentration increases. Both surface and bottom cohesive sediment concentrations decrease at Taipei Bridge and Pa–Ling Bridge. This implies that suspended sediment advected seaward and deposited. There is consequently a net seaward flux of suspended sediment near surface, and a net landward flux near the bed.     

Assimilation and Fractional Crystallization Controlled by Transport Process of Crustal Melt: Implications from an Alkali Basalt-Dacite Suite from Rishiri Volcano, Japan

Mechanisms of fractional crystallization with simultaneous crustalassimilation (AFC) are examined for the Kutsugata and Tanetomilavas, an alkali basalt–dacite suite erupted sequentiallyfrom Rishiri Volcano, northern Japan. The major element variationswithin the suite can be explained by boundary layer fractionation;that is, mixing of a magma in the main part of the magma bodywith a fractionated interstitial melt transported from the mushyboundary layer at the floor. Systematic variations in SiO₂ correlatewith variations in the Pb, Sr and Nd isotopic compositions ofthe lavas. The geochemical variations of the lavas are explainedby a constant and relatively low ratio of assimilated mass tocrystallized mass (‘r value’). In the magma chamberin which the Kutsugata and Tanetomi magmas evolved, a strongthermal gradient was present and it is suggested that the marginalpart of the reservoir was completely solidified. The assimilantwas transported by crack flow from the partially fused floorcrust to the partially crystallized floor mush zone throughfractures in the solidified margin, formed mainly by thermalstresses resulting from cooling of the solidified margin andheating of the crust. The crustal melt was then mixed with thefractionated interstitial melt in the mushy zone, and the mixedmelt was further transported by compositional convection tothe main magma, causing its geochemical evolution to be characteristicof AFC. The volume flux of the assimilant from the crust tothe magma chamber is suggested to have decreased progressivelywith time (proportional to t^–1/2), and was about 3 x 10^–2m/year at t = 10 years and 1 x 10^–2 m/year at t = 100years. It has been commonly considered that the heat balancebetween magmas and the surrounding crust controls the couplingof assimilation and fractional crystallization processes (i.e.absolute value of r). However, it is inferred from this studythat the ratio of assimilated mass to crystallized mass canbe controlled by the transport process of the assimilant fromthe crust to magma chambers. KEY WORDS: assimilation and fractional crystallization; mass balance model; magma chamber; melt transport; Pb isotope     

Hillslope erosion measurement—a simple approach to a complex process

The measurement of hillslope erosion can be a difficult, costly and time‐consuming activity. Many techniques are available, ranging from using environmental tracers, to LiDAR. Erosion measurements using erosion pins are assessed and compared with regional scale erosion data, hillslope data obtained using ¹³⁷Cs and erosion modelling results. The pins produced erosion rates which are within the range determined using ¹³⁷Cs and model data but above that of regional denudation rates. Our findings demonstrate that inexpensive erosion pins can provide reliable data on hillslope erosion. © 2015 Commonwealth of Australia. Hydrological Processes © 2015 John Wiley & Sons Ltd.     

Episodic granite accumulation and extraction from the mid‐crust

The processes of long‐range granitic magma transfer from mid‐ and lower crustal anatectic zones to upper crustal pluton emplacement sites remain controversial in the literature. This is partly because feeder networks that could have accommodated this large‐scale magma transport remain elusive in the field. Existing granite ascent models are based largely on numerical and theoretical studies that seek to demonstrate the viability of fracture‐controlled magma transport through dykes or self‐propagating hydrofractures. In most cases, the models present very little supporting field evidence, such as sufficiently voluminous near‐ or within‐source magma accumulations, to support their basic premises. We document large (deca‐ to hectometre‐scale), steeply dipping and largely homogeneous granite lenses in suprasolidus (~5 kbar, ~750 °C) mid‐crustal rocks in the Damara Belt in Namibia. The lenses are surrounded by and connected to shallowly dipping networks of stromatic leucogranites in the well‐layered gneisses of the deeply incised Husab Gorge. The outcrops define a four‐stage process from (i) the initial formation and growth of large, subvertical magma‐filled lenses as extension fractures developed at high angles to the subhorizontal regional extension in relatively competent wall‐rock layers. This stage is followed by (ii) the simultaneous lateral inflation and (iii) subcritical vertical growth of the lenses to a critical length that (iv) promotes fracture destabilization, buoyancy‐driven upward fracture mobilization and, consequently, vertical magma transport. These field observations are compared with existing numerical models and are used to constrain, by referring to the dimensions of the largest preserved inflated leucogranite lens, an estimate of the minimum fracture length (~100 m) and volume (~2.4 × 10⁵ m³) required to initiate buoyancy‐driven brittle fracture propagation in this particular mid‐crustal section. The critical values and field relationships compare favourably with theoretical models of magma ascent along vertical self‐propagating hydrofractures which close at their tails during propagation. This process leaves behind subtle wake‐like structures and thin leucogranite trails that mark the path of magma ascent. Reutilization of such conduits by repeated inflation and drainage is consistent with the episodic accumulation and removal of magma from the mid‐crust and is reflected in the sheeted nature of many upper crustal granitoid plutons.