Variability of dry sediment bulk density between and within retention ponds and its impact on the calculation of sediment yields

Monitoring sediment yields from catchments is important for assessing overall denudation rates and the impact of environmental change. One of the methods used to assess sediment yield is by quantifying sedimentation rates in reservoirs, lakes or small ponds. Before reliable sediment yield values (t ha^?1 a^?1) can be computed from such sedimentation records, the measured sediment volumes need to be converted to sediment masses using representative values of the dry sediment bulk density. In textbooks, simple relations predicting dry sediment bulk density from sediment texture, time since deposition and hydrologic condition are presented. In this study, 13 small flood retention ponds in central Belgium were sampled to reveal the variability in dry sediment bulk density and to test the commonly used relations to predict dry sediment bulk density. Dry sediment bulk density varies not only between the selected ponds (0·78–1·35 t m^?3) but also within individual ponds (coefficient of variation at 95 per cent ranges from 7 to 80 per cent). The observed variability can be attributed primarily to the hydrologic condition of the retention pond and, also, to sediment texture. The existing relations are not a reliable predictor for the observed dry bulk densities, because they are primarily based on sediment texture. Thus, when using volumetric sedimentation data from small ponds with varying hydrologic condition to predict sediment yield, existing relations predicting dry sediment bulk density cannot be applied. Instead, frequent and dense sampling of sediments is necessary to calculate a representative value of the dry sediment bulk density. Copyright © 2001 John Wiley & Sons, Ltd.     

DuMu: DUNE for multi-{phase, component, scale, physics, …} flow and transport in porous media

DuMu^x is a free and open-source simulator for flow and transport processes in porous media, based on the Distributed and Unified Numerics Environment (DUNE). Its main intention is to provide a sustainable and consistent framework for the implementation and application of model concepts, constitutive relations, discretizations, and solvers. The paper provides an overview of DuMu^x with the focus on software-related aspects. Selected examples highlight the multi-scale and the parallel capabilities.     

Factors controlling the morphology and volume (V)–length (L) relations of permanent gullies in the northern Ethiopian Highlands

Small‐scale aerial photographs and high‐resolution satellite images, available for Ethiopia since the second half of the twentieth century as for most countries, allow only the length of gullies to be determined. Understanding the development of gully volumes therefore requires that empirical relations between gully volume (V) and length (L) are established in the field. So far, such V–L relations have been proposed for a limited number of gullies/environments and were especially developed for ephemeral gullies. In this study, V–L relations were established for permanent gullies in northern Ethiopia, having a total length of 152 km. In order to take the regional variability in environmental characteristics into account, factors that control gully cross‐sectional morphology were studied from 811 cross‐sections. This indicated that the lithology and the presence of check dams or low‐active channels were the most important controls of gully cross‐sectional shape and size. Cross‐sectional size could be fairly well predicted by their drainage area. The V–L relation for the complete dataset was V = 0 · 562 L ^1·381 (n = 33, r² = 0 · 94, with 34 · 9% of the network having check dams and/or being low‐active). Producing such relations for the different lithologies and percentages of the gully network having check dams and/or being low‐active allows historical gully development from historical remote sensing data to be assessed. In addition, gully volume was also related to its catchments area (A) and catchment slope gradient (S_c). This study demonstrates that V–L and V–A × S_c relations can be very suitable for planners to assess gully volume, but that the establishment of such relations is necessarily region‐specific. Copyright © 2013 John Wiley & Sons, Ltd.     

Converting Radial Distortions Due to Geopotential Models to Geopotential Stokes Coefficients

A methodology for improving geopotential models has been developed. Theoretical relations have been derived converting coefficients in harmonic expansions for radial distortions due to geopotential models into geopotential Stokes coefficients. Terms of the order of 10 ^–10 in magnitude have been retained.     

Attenuation of Coda Waves in the Saurashtra Region, Gujarat (India)

The attenuation characteristics based on coda waves of two areas—Jamnagar and Junagarh of Saurashtra, Gujarat (India)—have been investigated in the present study. The frequency dependent relationships have been developed for both the areas using single back scattering model. The broadband waveforms of the vertical components of 33 earthquakes (Mw 1.5–3.5) recorded at six stations of the Jamnagar area, and broadband waveforms of 68 earthquakes (Mw 1.6–5) recorded at five stations of the Junagarh area have been used for the analysis. The estimated relations for the Junagarh area are: Q _c?=?(158?±?5)f^(0.99±0.04) (lapse time : 20?s), Q _c?=?(170?±?4.4)f^(0.97±0.02) (lapse time : 30?s) and Q _c?=?(229?±?6.6)f^(0.94±0.03) (lapse time : 40?s) and for the Jamnagar area are: Q _c?=?(178?±?3)f^(0.95±0.05) (lapse time : 20?s), Q _c?=?(224?±?6)f^(0.98±0.06) (lapse time : 30?s) and Q _c?=?(282?±?7)f^(0.91±0.03) (lapse time : 40?s). These are the first estimates for the areas under consideration. The Junagarh area appears to be more attenuative as compared to the Jamnagar area. The increase in Q _c values with lapse time found here for both the areas show the depth dependence of Q _c as longer lapse time windows will sample larger area. The rate of decay of attenuation (Q ^?1) with frequency for the relations obtained here is found to be comparable with those of other regions of the world though the absolute values differ. A comparison of the coda-Q estimated for the Saurashtra region with those of the nearby Kachchh region shows that the Saurashtra region is less heterogeneous. The obtained relations are expected to be useful for the estimation of source parameters of the earthquakes in the Saurashtra region of Gujarat where no such relations were available earlier. These relations are also important for the simulation of earthquake strong ground motions in the region.     

LONG-TERM PERSPECTIVE ON THE NATURE OF THE AIR–WATER TEMPERATURE RELATIONSHIP: A CASE STUDY

Monthly mean data from a 90 year period relating to a small catchment (142.4 km²) in north-central Austria were used to provide a long-term perspective on the nature of the air–water temperature relationship. Annual mean values of air and water temperature were related in a relatively insensitive and scattered way (r² < 55%, b < 0.65), whereas the relationship for monthly mean values was closer and steeper (r² > 95%, b > 0.65). A separate regression equation was needed to describe the behaviour of monthly mean water temperatures as the air temperatures fell below freezing. Analysis of air–water temperature regressions for individual months revealed a series of relations which were generally more scattered and less, but variously, sensitive than the ensemble relationship of monthly mean values. Monthly mean water temperatures could be predicted from the ensemble air–water temperature relationship and from the relations for individual months with root mean square errors of > 1.0 and < 0.8°C, respectively. Segmentation of air–water temperature regressions according to air temperatures above and below freezing did not significantly improve water temperature prediction. Hysteresis in, and the relatively low slope of, the air–water temperature relationships in the study catchment reflected the importance of snowmelt in the flow regime. © 1997 by John Wiley & Sons Ltd.     

Phosphorus export dynamics and hydrobiogeochemical controls across gradients of scale,topography and human impact

Concentration‐discharge (c‐Q) plots are routinely used as an integrated signal of watershed response to infer solute sources and travel pathways. However, the interpretation of c‐Q data can be difficult unless these data are fitted using statistical models. Such models are frequently applied for geogenic solutes, but it is unclear to what extent they might aid in the investigation of nutrient export patterns, particularly for total dissolved phosphorus (TDP) which is a critical driver of downstream eutrophication problems. The goal of the present study was therefore to statistically model c‐Q relations (where c is TDP concentrations) in a set of contrasting watersheds in the Northern Great Plains—ranging in size from 0.2 to 1000+ km²—to assess the controls of landscape properties on TDP transport dynamics. Six statistical models were fitted to c‐Q data, notably (a) one linear model, (b) one model assuming that c‐Q relations are driven by the mixing of end‐member waters from different landscape locations (i.e., hydrograph separation), (c) one model relying on a biogeochemical stationarity hypothesis (i.e., power law), (d) one model hypothesizing that c‐Q relations change as a function of the solute subsurface contact time (i.e., hyperbolic model), and (e) two models assuming that solute fluxes are mostly dependent on reaction rates (i.e., chemical models). Model performance ranged from mediocre (R² < 0.2) to very good (R² > 0.9), but the hydrograph separation model seemed most universal. No watershed was found to exhibit chemostatic behaviour, but many showed signs of dilution or enrichment behaviour. A tendency toward a multi‐model fit and better model performance was observed for watersheds with moderate slope and higher effective drainage area. The relatively poor model performance obtained outside these conditions illustrates the likely importance of controls on TDP concentrations in the region that are independent of flow dynamics.     

Heavy element affinities in Apollo 17 samples

²⁰⁴Pb, Bi, Tl and Zn in samples from Apollo 17 exhibit relationships not found in samples from other sites.²⁰⁴Pb, Tl and Zn in residues remaining after dilute acid leaching are correlated with one another. Orange soil 74220, which is enriched in²⁰⁴Pb, Tl and Zn, is included in these relationships. In addition the submicron metallic phase generally associated with agglutinate formation is correlated with all three of these elements; this relationship has already been reported for²⁰⁴Pb in other samples. Thus, orange soil and agglutinates appear to be involved in concentrating heavy volatile metals. A process other than mixing is required to account for this. As a consequence of the isolation of the landing site by the surrounding massifs, local supply and recycling of volatile trace elements in soils may account for some of the interelement relations.     

Long-range Earthquake Forecasting with Every Earthquake a Precursor According to Scale

Scaling relations previously derived from examples of the precursory scale increase before major earthquakes show that the precursor is a long-term predictor of the time, magnitude, and location of the major earthquake. These relations are here taken as the basis of a stochastic forecasting model in which every earthquake is regarded as a precursor. The problem of identifying those earthquakes that are actually precursory is thus set aside, at the cost of limiting the strength of the resulting forecast. The contribution of an individual earthquake to the future distribution of hazard in time, magnitude and location is on a scale determined, through the scaling relations, by its magnitude. Provision is made for a contribution to be affected by other earthquakes close in time and location, e.g., an aftershock may be given low weight. Using the New Zealand catalogue, the model has been fitted to the forecasting of shallow earthquakes exceeding magnitude 5.75 over the period 1965–2000. It fits the data much better than a baseline Poisson model with a location distribution based on proximity to the epicenters of past earthquakes. Further, the model has been applied, with unchanged parameters, to the California region over the period 1975–2001. There also, it performs much better than the baseline model fitted to the same region over the period 1951–1974; the likelihood ratio is 10¹⁵ in favor of the present model. These results lend credence to the precursory scale increase phenomenon, and show that the scaling relations are pervasive in earthquake catalogues. The forecasting model provides a new baseline model against which future refinements, and other proposed models, can be tested. It may also prove to be useful in practice. Its applicability to other regions has still to be established.     

Iron oxidation state in lower mantle mineral assemblages: II. Inclusions in diamonds from Kankan, Guinea

Inclusions of ferropericlase and former (Mg,Fe)(Si,Al)O₃ perovskite in diamonds from Kankan, Guinea believed to originate in the lower mantle were studied using Mössbauer spectroscopy to determine Fe³⁺/ΣFe. Fe³⁺ concentration in the (Mg,Fe)(Si,Al)O₃ inclusion is consistent with empirical relations relating Fe³⁺/ΣFe to Al concentration, supporting the inference that it crystallised in the perovskite structure at lower mantle conditions. In ferropericlase there is a nearly linear variation of trivalent cation abundance with monovalent cation abundance, suggesting a substitution of the form Na_0.5M_0.5³⁺O (M=Fe³⁺, Cr³⁺, Al³⁺). Excess positive charge is likely balanced by cation vacancies, where their abundance is observed to increase with increasing iron concentration, consistent with high-pressure experiments. The abundance of cation vacancies is related to oxygen fugacity, where ferropericlase inclusions from Kankan and São Luiz (Brazil) are inferred to have formed at conditions more oxidising than Fe-(Mg,Fe)O equilibrium, but more reducing than Re-ReO₂ equilibrium. Fe²⁺/Mg partition coefficients between (Mg,Fe)(Si,Al)O₃ and ferropericlase were calculated for inclusions co-existing in the same diamond using Mössbauer data and empirical relations based on high-pressure experimental work. Most values are consistent with high-pressure experiments, suggesting that these inclusions equilibrated at lower mantle conditions. The measured ferropericlase Fe³⁺ concentrations are consistent with diamond formation in a region of redox gradients, possibly arising from the subduction of oxidised material into reduced lower mantle. Reduction of carbonate to form ferropericlase and diamond is consistent with a slight shift of Kankan δ¹³C values to isotopically heavy compositions compared to the worldwide dataset, and could supply the oxygen necessary to satisfy the high Fe³⁺ concentration in (Mg,Fe)(Si,Al)O₃ perovskite, as well as account for the high proportion of ferropericlase in the lower mantle paragenesis. The heterogeneity of lower mantle diamond sources indicates that the composition of lower mantle diamonds do not necessarily reflect those of the bulk mantle.     

The Taylor Creek Rhyolite of New Mexico: a rapidly emplaced field of lava domes and flows

The Tertiary Taylor Creek Rhyolite of southwest New Mexico comprises at least 20 lava domes and flows. Each of the lavas was erupted from its own vent, and the vents are distributed throughout a 20 km by 50 km area. The volume of the rhyolite and genetically associated pyroclastic deposits is at least 100 km³ (denserock equivalent). The rhyolite contains 15%–35% quartz, sanidine, plagioclase, ±biotite, ±hornblende phenocrysts. Quartz and sanidine account for about 98% of the phenocrysts and are present in roughly equal amounts. With rare exceptions, the groundmass consists of intergrowths of fine-grained silica and alkali feldspar. Whole-rock major-element composition varies little, and the rhyolite is metaluminous to weakly peraluminous; mean SiO₂ content is about 77.5±0.3%. Similarly, major-element compositions of the two feldsparphenocryst species also are nearly constant. However, whole-rock concentrations of some trace-elements vary as much as several hundred percent. Initial radiometric age determinations, all K–Ar and fission track, suggest that the rhyolite lava field grew during a period of at least 2 m.y. Subsequent ⁴⁰Ar/³⁹Ar ages indicate that the period of growth was no more than 100 000 years. The time-space-composition relations thus suggest that the Taylor Creek Rhyolite was erupted from a single magma reservoir whose average width was at least 30 km, comparable in size to several penecontemporaneous nearby calderas. However, this rhyolite apparently is not related to a caldera structure. Possibly, the Taylor Creek Phyolite magma body never became sufficiently volatile rich to produce a large-volume pyroclastic eruption and associated caldera collapse, but instead leaked repeatedly to feed many relatively small domes and flows.The new ⁴⁰Ar/³⁹Ar ages do not resolve preexisting unknown relative-age relations among the domes and flows of the lava field. Nonetheless, the indicated geologically brief period during which Taylor Creek Rhyolite magma was erupted imposes useful constraints for future evaluation of possible models for petrogenesis and the origin of trace-element characteristics of the system.     

The use of mixing models and chemical geothermometers for estimating underground temperatures in geothermal systems

Application of various chemical geothermometers and mixing models indicate underground temperatures of 260°C, 280°C and 265°C in the Geysir, Hveravellir and Landmannalaugar geothermal fields in Iceland, respectively. Mixing of the hot water with cold water occurs in the upflow zones of all these geothermal systems. Linear relations between chloride, boron and δ¹⁸O constitute the main evidence for mixing, which is further substantiated by chloride, silica and sulphate relations in the Geysir and Hveravellir fields.A new carbonate-silica mixing model is proposed which is useful in distinguishing boiled and non-boiled geothermal waters. This model can also be used to estimate underground temperatures using data from warm springs. This model, as well as the chloride-enthalpy model and the Na-Li, and CO₂-gas geothermometers, invariably yield similar results as the quartz geothermometer sometimes also does. By contrast, the Na-K and the Na-K-Ca geothermometers yield low values in the case of boiling hot springs, largely due to loss of potassium from solution in the upflow. The results of these geothermometers are unreliable for mixed waters due to leaching subsequent to mixing.     

Iron oxidation state in lower mantle mineral assemblages: I. Empirical relations derived from high-pressure experiments

The oxidation state of iron can significantly influence the physical and chemical properties of lower mantle minerals. To improve methods for estimation of Fe³⁺/∑Fe, synthetic assemblages of (Mg,Fe)(Si,Al)O₃ perovskite and (Mg,Fe)O ferropericlase were synthesised from oxide starting mixtures in Re or Fe capsules at 26 GPa and 1650-1850 °C using a multianvil press. (Mg,Fe)(Si,Al)O₃ majorite was also present in some of the run products. Both electron energy loss spectra (EELS) and Mössbauer spectra were measured for each run product, and a robust fitting method was developed for Mössbauer spectra using EELS results as a standard that enabled Fe³⁺/∑Fe of (Mg,Fe)(Si,Al)O₃ perovskite to be determined from Mössbauer spectra of multiphase assemblages. There is a close to linear variation between Fe³⁺/∑Fe and Al concentration in (Mg,Fe)(Si,Al)O₃ perovskite, independent of oxygen fugacity. The concentration of Fe³⁺ in (Mg,Fe)O increases with increasing iron concentration along curves of constant oxygen fugacity, where higher oxygen fugacity stabilises greater Fe³⁺ concentrations. Fe²⁺/Mg partition coefficients calculated from chemical composition data corrected for measured Fe³⁺/∑Fe showed values nearly identical within experimental error for all samples, and independent of Al concentration and oxygen fugacity. Simple empirical relations were derived to calculate Fe³⁺/∑Fe in (Mg,Fe)(Si,Al)O₃ perovskite and (Mg,Fe)O ferropericlase samples for which no Mössbauer or EELS data were available, and tested by applying them to calculation of Fe²⁺/Mg partition coefficients from literature data for (Mg,Fe)(Si,Al)O₃ perovskite-(Mg,Fe)O assemblages where only total iron concentrations had been measured. Results showed Fe²⁺/Mg partition coefficients that were equal to existing values within experimental error, hence confirming the validity of the empirical relations.     

Radon anomalies and volcanic eruptions

A well-documented case of ²²²Rn anomaly preceding the eruptive activity of Karymsky volcano (Kamchatka) was recently reported in the literature. Stimulated by this example, we have attempted to utilize the available data on radon emanation from rocks, its solubility, and its circulation in waters to discuss how a ²²²Rn anomaly can be produced by magma approaching the surface. It is shown that the most likely process of radon release is the flushing of gases through pore fluids. Heating of extensive fracture surfaces by high-temperature gases may also be important. In order to survive in detectable amounts after moving distances greater than a few meters, radon must be transported by fast-moving fluids, such as those rising toward the surface through a fracture or in the ascending limbs of fluids convecting in very porous rocks under high temperature gradients. The pattern observed at Karymsky volcano is interpreted on the basis of these relations.     

Surfactant and irrigation effects on wettable soils: runoff,erosion, and water retention responses

Surfactants are chemical compounds that can change the contact angle of a water drop on solid surfaces and are commonly used to increase infiltration into water repellent soil. Since production fields with water repellent soil often contain areas of wettable soil, surfactants applied to such fields worldwide will likely be applied to wettable soil, with unknown consequences for irrigation‐induced erosion, runoff, or soil water relations. We evaluated surfactant and simulated sprinkler irrigation effects on these responses for three wettable, Pacific Northwest soils, Latahco and Rad silt loams, and Quincy sand. Along with an untreated control, we studied three surfactants: an alkyl polyglycoside (APG) in solution at a concentration of 18 g active ingredient (AI) kg^?1, a block copolymer at 26 g kg^?1, and a blend of the two at 43 g kg^?1. From 2005 to 2009 in the laboratory, each surfactant was sprayed at a rate of 46·8 l ha^?1 onto each soil packed by tamping into 1·2‐ by 1·5‐m steel boxes. Thereafter, each treated soil was irrigated twice at 88 mm h^?1 with surfactant‐free well water. After each irrigation, runoff and sediment loss were measured and soil samples were collected. While measured properties differed among soils and irrigations, surfactants had no effect on runoff, sediment loss, splash loss, or tension infiltration, compared to the control. Across all soils, however, the APG increased volumetric water contents by about 3% (significant at p≤0·08) at matric potentials from 0 to ? 20 kPa compared to the control. With a decrease in the liquid–solid contact angle on treated soil surfaces, surfactant‐free water appeared able to enter, and be retained in pores with diameters ≥ 15 µm. All told, surfactants applied at economic rates to these wettable Pacific Northwest soils posed little risk of increasing either runoff or erosion or harming soil water relations. Moreover, by increasing water retention at high potentials, surfactants applied to wettable soils may allow water containing pesticides or other agricultural chemicals to better penetrate soil pores, thereby increasing the efficacy of the co‐applied materials. Copyright © 2010 John Wiley & Sons, Ltd.     

Application of CCA for study on modern lake diatoms and environment in the Tibetan Plateau

The relations between lake surface sediment diatoms and water environmental variables were revealed effectively by use of a new multivariate canonical correspondence analysis (CCA) based on 45 lakes in the Tibetan Plateau. Water depth, conductivity, Cl^?, Mg²⁺, K⁺ and pH, identified from 12 contemporary water environmental variables, can significantly and independently explain the diatom distributions (p<0.05). The first two axes (λ1=0.34, λ2=0.27) capture 16.1% of the variance in the species data, and account for 57.4% of the variance in diatom-environment relationship. The deletion of redundant environmental variables and unusual samples do not influence the explanation to diatom data. The final CCA result indicates that the water depth and the salinity are the two important environmental gradients and influence the diatom distribution in the plateau lakes. The water depth correlates with axis 1, while conductivity, Cl^?, Mg²⁺ and K⁺, indicating the direction of salinity changes, correlate with both of the first two axes. The definition of diatom-environment model may provide a basis for further quantitative inference on diatom-environment transfer function.     

Spatio‐temporal patterns of nutrient concentrations and export in a north‐eastern European lowland catchment

In 2002–2004 we undertook six sampling campaigns during representative hydrological stages in a 901 km² Estonian lowland catchment to quantify the spatial and seasonal variability of in‐stream dissolved inorganic nitrogen (DIN) and dissolved reactive phosphorus (DRP) concentrations and to identify the influence of land cover and landscape structure. Using a synoptic approach we mapped concentrations in all stream orders. Using linear regression, the relations between the share of agricultural land and log‐transformed in‐stream concentrations were explored. Both the share of agricultural land in the entire ‘area of influence’ upstream from a sampling location, as well as the share in a 150‐m buffer around the stream were used as linear regression input variables. Log‐transformed DIN and DRP concentration variability was highest for lower order streams, while it averaged out in higher order streams during all seasons. Between‐season variation in export can mainly be attributed to discharge variation. In extremely dry periods, there are no significant relations between land cover/structure and in‐stream ln(DIN) concentrations and only weak relations for ln(DRP) concentrations. In other seasons, the share of agricultural land in the upstream area can explain concentrations in higher order streams better than in lower order streams. The prediction of ln(DIN) concentrations in lower order streams can be improved by using the share of agricultural land in a 150‐m buffer as an input variable. This indicates that hydrological connectivity must be taken into account for lower order streams, while land cover shares are enough to explain concentrations for higher order streams. Copyright © 2009 John Wiley & Sons, Ltd.     

Thorium isotopes （22STh, 230Th, 232Th） and applications in reconstructing the Yangtze and Yellow River floods

In the past decades, the floods of the Yangtze and Yellow River introduced unexpected changes of the ecological community and sedimentary dynamics in the East China Sea （ECS）. To reconstruct the flood events in the ECS, 228Th, 230Th and 232Th have been examined in a sediment core. The specific activities of three thorium isotopes have good positive relations with fine fractions （〈63 μm）, indicating that Th activity concentrations heavily depend upon the sediment grain size. The size-normalized activities of 228Th, 23-Th and 232Th showed significant variations. Coincidences between the higher Th activities and historical floods of the Yangtze and Yellow River demonstrated that size-normalized Th recorded the two rivers＇ flood events. The activity ratios of thorium isotopes, i.e. 230Th/232Th and 228Th/232Th, also showed similar patterns to the historical river floods. In three periods （1740s, 1840-1860s and 1930-1960s）, characterized by frequent floods, the thorium activity ratios were fairly low and close to the Yangtze and Yellow River estuary sediments, coinciding with the less oceanic 228Th and 230Th contributions during the flooding periods. Accordingly, these results support the size-normalized Th activity and thorium ratios as proxies of the river floods in coastal seas.     

Shaping and maintaining a medium-sized main channel in the Lower Yellow River

This paper studies relations between bankfull discharge,lateral cross section variation and the incoming flow and sediment condition in the Lower Yellow River using measured data from 1950 to 2003.Since 1950 the bankfull discharge has obviously decreased and the ratio of channel width to flow depth has increased.The critical annual average incoming sediment coefficient（defined as the ratio of sediment concentration to discharge） and discharge at the Huayuankou station are approximately 0.012 and 1,850 m³s^-1,respectively,for no accumulative deposition occurring in the reach from Huayuankou to Lijin.On this basis,a mathematical model is used to study the scale of the main channel in the Lower Yellow River and its corresponding bankfull discharge under possible incoming flow and sediment conditions in the near future.The main factors influencing the scale of the main channel are analyzed,and measures to shape and maintain a medium-sized channel are discussed.The results show the effect of various water and sediment combinations released from the Xiaolangdi Reservoir on the shaping of the main channel and suggest that under recent incoming flow and sediment conditions,it is possible to shape and maintain a medium-sized channel with a bankfull discharge of approximate 4,000 m³ s^-1.     

Experimental study of the velocity of density currents in convergent and divergent channels

The head velocity of the density current in the convergent and divergent channel is a key parameter for evaluating the extent to which suspended material travels, and for determining the type and distribution of sediment in the water body. This study experimentally evaluated the effects of the reach degree of convergence and divergence on the head velocity of the density current. Experiments were conducted in the flume with 6.0 m long, 0.72 m width and 0.6 m height. The head velocity was measured at three convergent degrees （-8°; -12°; -26°）, at three divergent degrees （8°; 12°; 26°） and two slopes （0.009, 0.016） for various discharges. The measured head velocity of the density current is compared with the head velocity of the density current in the constant cross section channel. Based on non-dimensional and statistical analysis, relations as linear multiple regression are offered for predicting head velocity of the density current in the convergent, divergent and constant cross section channel. Also the results of this research show that for the same slope and discharge, the head velocity of the density current in the convergent and divergent channel are greater and less than the head velocity of the constant cross section, respectively.