Depth–energy and depth–force relationships in open channel flows: Analytical findings

In the present work the depth–specific energy relationship and the depth–total force relationship in open channel flows of wide rectangular cross-section are analytically inverted.     

Analytical inversion of specific energy–depth relationship in channels with parabolic cross-sections

Abstract

The specific energy–depth relationship in open channels with parabolic cross-sections is analytically inverted. Two nondimensional expressions of the specific energy are considered, depending on the prescribed quantity (specific energy or alternate depth). The inversion of such functions consists of finding the roots of cubic and quartic equations. By solving a quartic equation for a given discharge and for each value of the specific energy, a subcritical depth and a supercritical depth are found analytically. In this case, two acceptable roots are recognized and two other roots are discarded, on the basis of their physical meaning. Moreover, by solving a cubic equation for a given discharge and for each value of the depth, the other corresponding alternate depth is found analytically. Then, one acceptable root is recognized and two other complex conjugate roots are discarded. Finally, different examples are presented, to show the efficiency of the proposed solutions. Such analytical solutions can be easily used in natural rivers and parabolic channels.

Citation Vatankhah, A. R. & Valiani, A. (2011) Analytical inversion of specific energy–depth relationship in channels with parabolic cross-sections. Hydrol. Sci. J. 56(5), 834–840.     

Temporal variation in community–environment relationships and stream classifications in benthic diatoms: Implications for bioassessment

Diatom based biomonitoring tools are proven to be a practical way to indicate stream conditions, but only little is known how stable diatom inferred stream classifications are in time. We studied annual variation in diatom community–environment relationships, community structure and diatom indices (Index of Pollution Sensitivity, IPS and Trophic Diatom Index, TDI) during three consecutive years (2010–2012) in four drainage basins distributed in Finland. We also used a cluster analysis to examine if stream classifications resulted in distinct and temporally stable community types. We found only subtle temporal variation in community–environment relationships, nutrient concentrations (N, P) and conductivity consistently being the main factors structuring communities. According to Mantel tests, diatom communities resembled each other significantly at different years, and the values of IPS and TDI indices remained relatively stable in the basins through time. The stream classification based on diatoms also resulted in temporally stable and statistically distinct community types. We thus suggest that sampling of diatoms, e.g., in every three years seems to be a reliable procedure to assess biological water quality. Generally, to choose the correct metrics to assess water quality is essential, and both biological and physicochemical factors should be considered.     

Advected fields in maps: II. Dynamo action in the stretch–fold–shear map

The growth of magnetic field is considered in the stretch–fold–shear map in the limit of weak diffusion. Numerical results are given for insulating, perfectly conducting and periodic boundary conditions. The resulting eigenvalue branches and magnetic fields are related to eigenvalue branches for perfect dynamo action, obtained for zero diffusion using a complex variable formulation.

The effect of diffusion on these perfect dynamo modes depends on their structure, growth rate and the diffusive boundary conditions employed. In some cases, the effect of diffusion is a small perturbation, giving a correction going to zero in the limit of weak diffusion, with a scaling exponent given analytically. In other cases weak diffusion can entirely destroy a perfect dynamo branch. Diffusive boundary layers can also generate entirely new branches.

These different cases are elucidated, and within the framework of the asymptotic approximations used (which do not constitute a rigorous proof), it is seen that for all three boundary conditions employed, the stretch–fold–shear map is a fast dynamo.     

Comparing flow regime,channel hydraulics,and biological communities to infer flow–ecology relationships in the Mara River of Kenya and Tanzania

Abstract

Equatorial rivers of East Africa exhibit unusually complex seasonal and inter-annual flow regimes, and aquatic and adjacent terrestrial organisms have adapted to cope with this flow variability. This study examined the annual flow regime over the past 40 years for three gauging stations on the Mara River in Kenya and Tanzania, which is of international importance because it is the only perennial river traversing the Mara-Serengeti ecoregion. Select environmental flow components were quantified and converted to ecologically relevant hydraulic variables. Vegetation, macroinvertebrates, and fish were collected and identified at target study sites during low and high flows. The results were compared with available knowledge of the life histories and flow sensitivities of the riverine communities to infer flow–ecology relationships. Management implications are discussed, including the need to preserve a dynamic environmental flow regime to protect ecosystems in the region. The results for the Mara may serve as a useful model for river basins of the wider equatorial East Africa region.

Editor Z.W. Kundzewicz; Guest editor M. Acreman     

Analytical solutions for sequentially coupled one-dimensional reactive transport problems – Part II: Special cases,implementation and testing

This is Part-II of a two-part article that presents analytical solutions to multi-species reactive transport equations coupled through sorption and sequential first-order reactions. In Part-I, we provide the mathematical derivations and in this article we discuss the computational techniques for implementing these solutions. We adopt these techniques to develop a general computer code and use it to verify the solutions. We also simplify the general solutions for various special-case transport scenarios involving zero initial condition, identical retardation factors and zero advection. In addition to this, we derive specialized solution expressions for zero dispersion and steady-state conditions. Whereever possible, we compare these special-case solutions against previously published analytical solutions to establish the validity of the new solution. Finally, we test the new solution against other published analytical and semi-analytical solutions using a set of example problems.     

Runoff-generated debris flows: Observation of initiation conditions and erosion–deposition dynamics along the channel at Cancia (eastern Italian Alps)

In the Dolomitic region, abundant coarse hillslope sediment is commonly found at the toe of rocky cliffs. Ephemeral channels originate where lower permeability bedrock surfaces concentrate surface runoff. Debris flows initiate along such channels following intense rainfall and determine the progressive erosion and deepening of the channels. Sediment recharge mechanisms include rock fall, dry ravel processes and channel-bank failures. Here we document debris flow activity that took place in an active debris flow basin during the year 2015. The Cancia basin is located on the southwestern slope of Mount Antelao (3264 m a.s.l.) in the dolomitic region of the eastern Italian Alps. The 2.5 km² basin is incised in dolomitic limestone rocks. The data consist of repeated topographic surveys, distributed rainfall measurements, time-lapse (2 s) videos of two events and pore pressure measurements in the channel bed. During July and August 2015, two debris flow events occurred, following similarly intense rainstorms. We compared rainfall data to existing rainfall triggering thresholds and simulated the hydrological response of the headwater catchment with a distributed model in order to estimate the total and peak water discharge. Our data clearly illustrate how debris entrainment along the channel is the main contributor to the overall mobilized volume and that erosion is dominant when the channel slope exceeds 16°. Further downstream, sediment accumulation and depletion occurred alternately for the two successive events, indicating that sediment availability along the channel also influences the flow behaviour along the prevailing-transport reach. The comparison between monitoring data, topographical analysis and hydrological simulation allows the estimation of the average solid concentration of the two events and suggests that debris availability has a significant influence on the debris flow volume. © 2020 John Wiley & Sons, Ltd.     

The near-field method: a modified equivalent linear method for dynamic soil–structure interaction analysis. Part II: verification and example application

Usually for modeling of soil in a direct soil–structure interaction (SSI) problem, the equivalent linear soil properties are used. However, this approach is not valid in the vicinity of a foundation, where the soil experiences large strains and a high level of nonlinearity because of structural vibrations. The near-field method was developed and described in a companion paper to overcome this limitation. This method considers the effects of large strains and suggests a shear modulus and a damping ratio further modified in the near-field of a foundation. Validity and performance of this approach are evaluated, application examples are explained and the results of a parametric study about the role of soil and structure parameters in the extent of SSI effects on the nonlinear seismic response of structures are presented in this paper. One real existing and five, ten, fifteen and twenty story moment-resisting frame steel buildings with two different site conditions corresponding to firm and soft soils are considered and the responses obtained from the near-field method are compared with the recorded and rigorous responses. Moreover, various SSI modeling techniques are employed to investigate the accuracy and performance of each approach. The results show that the near-field method is a simple yet accurate enough approach for analysis of direct SSI problems.     

Space–time relationships between volcanic associations of different alkalinities: The Belogolovskii Massif,Sredinnyi Range,Kamchatka. Part II. Geochemistry of volcanic rocks and magma sources

Data are presented relating to volcanic series in the Belogolovskii Massif, Sredinnyi Range, Kamchatka. We discuss new geochronologic data, the distributions of rare elements and platinum elements in the rocks, and list the isotope characteristics of volcanic series with normal and moderate alkalinities. We show that the Late Pliocene to Early Pleistocene rocks that belong to the moderate alkaline series of the Belogolovskii volcanic massif are different from rocks in the normally alkaline series of the Late Miocene to Middle Pliocene volcanogenic basement in having higher concentrations of the HFSE and LILE components. We propose a model for the generation of moderate alkaline magmas involving a heterogeneous depleted and a heterogeneous enriched source of material. According to the isotope data, one of these sources may be the subducted oceanic lithosphere of the Pacific and the Commander-Islands type, while the other source was recycled material of the Indian MORB type.     

Empirical statistical characterization and regionalization of amplitude–duration–frequency curves for extreme peak flows in the Lake Victoria Basin,East Africa

Abstract

This paper focuses on a regionalization attempt to partly solve data limitation problems in statistical analysis of high flows to derive discharge–duration–frequency (QDF) relationships. The analysis is based on 24 selected catchments in the Lake Victoria Basin (LVB) in East Africa. Characteristics of the theoretical QDF relationships were parameterized to capture their slopes of extreme value distributions (evd), tail behaviour and scaling measures. To enable QDF estimates to be obtained for ungauged catchments, interdependence relationships between the QDF parameters were identified, and regional regression models were developed to explain the regional difference in these parameters from physiographic characteristics. In validation of the regression models, from the lowest (5 years) to the highest (25 years) return periods considered, the percentage bias in the QDF estimates ranged from –2% for the 5-year return period to 27% for 25-year return period.

Editor D. Koutsoyiannis     

Pre-explosive conduit conditions of the 1997 Vulcanian explosions at Soufrière Hills Volcano,Montserrat: II. Overpressure and depth distributions

A type example of Vulcanian eruptive dynamics is the series of 88 explosions that occurred between August and October 1997 at Soufrière Hills volcano on Montserrat Island. These explosions are interpreted to be caused by the pressurization of a conduit by a shallow highly crystalline and degassed magma plug. We test such an interpretation by combining the pressures and porosities of the pre-explosive magma column proposed by Burgisser et al. (2010, doi:10.1016/j.jvolgeores.2010.04.008) into a physical model that reconstructs a depth-referenced density profile of the column for four mechanisms of pressure buildup. Each mechanism yields a different overpressure profile: 1) gas accumulation, 2) conduit wall elasticity, 3) microlite crystallization, and 4) magma flowage. For the first three mechanisms, the three-part vertical layering of the conduit prior to explosion was spatially distributed as a dense cap atop the conduit with a thickness of a few tens of meters, a transition zone of 400–700 m with heterogeneous vesicularities, and, at greater depth, a more homogeneous, low-porosity zone that brings the total column length to ~ 3.5 km. A shorter column can be obtained with mechanism 4: a dense cap of less than a few meters, a heterogeneous zone of 200–500 m, and a total column length as low as 2.5 km. Inflation/deflation cycles linked to a periodic overpressure source offer a dataset that we use to constrain the four overpressure mechanisms. Magma flowage is sufficient to cause periodic edifice deformation through semi-rigid conduit walls and build overpressures able to trigger explosions. Gas accumulation below a shallow plug is also able to build such overpressures and can occur regardless of magma flowage. The concurrence of these three mechanisms offers the highest likelihood of building overpressures leading to the 1997 explosion series. We also explore the consequences of sudden (eruptive) overpressure release on our magmatic columns to assess the role of syn-explosive vesiculation and pre-fragmentation column expansion. We find that large shallow overpressures and efficient syn-explosive vesiculation cause the most dramatic pre-fragmentation expansion. This leads us to depict two end-member pictures of a Vulcanian explosion. The first case corresponds to the widely accepted view that the downward motion of a fragmentation front controls column evacuation. In the second case, syn-explosive column expansion just after overpressure release brings foamed-up magma up towards an essentially stationary and shallow fragmentation front.     

Lava channel roofing,overflows, breaches and switching: insights from the 2008–2009 eruption of Mt. Etna

During long-lived basaltic eruptions, overflows from lava channels and breaching of channel levées are important processes in the development of extensive 'a'ā lava flow-fields. Short-lived breaches result in inundation of areas adjacent to the main channel. However, if a breach remains open, lava supply to the original flow front is significantly reduced, and flow-field widening is favoured over lengthening. The development of channel breaches and overflows can therefore exert strong control over the overall flow-field development, but the processes that determine their location and frequency are currently poorly understood. During the final month of the 2008–2009 eruption of Mt. Etna, Sicily, a remote time-lapse camera was deployed to monitor events in a proximal region of a small ephemeral lava flow. For over a period of ~10 h, the flow underwent changes in surface elevation and velocity, repeated overflows of varying vigour and the construction of a channel roof (a required prelude to lava tube formation). Quantitative interpretation of the image sequence was facilitated by a 3D model of the scene constructed using structure-from-motion computer vision techniques. As surface activity waned during the roofing process, overflow sites retreated up the flow towards the vent, and eventually, a new flow was initiated. Our observations and measurements indicate that flow surface stagnation and flow inflation propagated up-flow at an effective rate of ~6 m h⁻¹, and that these processes, rather than effusion rate variations, were ultimately responsible for the most vigorous overflow events. We discuss evidence for similar controls during levée breaching and channel switching events on much larger flows on Etna, such as during the 2001 eruption.     

Some relationships between VAN Seismic Electric Signals (SES) and earthquake parameters: Study for the 1995–2007 period

A 12-year period experimental data, from 1 January 1995, to 20 August 2007, have been examined for possible correlations between Seismic Electric Signals (SES) of the VAN method and source parameters of the corresponding earthquakes in Western Greece. During that period 13 earthquakes of magnitude M _W ≥ 5 with epicenters in the area 19–24°E, 36–41°N and available the CMT solutions (Centroid Moment Tensor focal mechanism solutions) have been found to be preceded by a SES each time at one of the three VAN stations of IOA, PIR or PAT. The results of IOA and PIR stations are compared to those reported by Uyeda et al. 1999, for the previous period 1983–1994. The IOA station seems to be sensitive to earthquakes with thrust type mechanisms being mainly consistent to its past behavior. The PIR station, which is moved from its previous position by few kilometers, exhibits some changes. It detects now both strike slip and thrust type earthquakes and becomes sensitive to new areas indicating probably the strong dependence of the SES station detection ability upon its location.     

Incentives for field hydrology and data sharing: collaboration and compensation: reply to “A need for incentivizing field hydrology,especially in an era of open data”*

We thank Allen and Berghuijs for continuing the discussion on field hydrology and data sharing and discuss two incentives to promote data collection and sharing in hydrological sciences: a collaborative attitude and additional funding to make data publicly available.     

The construction of a reliable absolute chronology for the last two millennia in an anthropogenically disturbed peat bog: Limitations and advantages of using a radio-isotopic proxy and age–depth modelling

The main aim of this paper is to present the pitfalls connected with the construction of reliable chronologies for anthropogenically disturbed peatlands over the last two millennia based on ²¹⁰Pb and ¹⁴C dating, i.e. the period of the strongest human impact on these ecosystems. The following hypotheses have been formulated: i) parts of peatlands suspected to be affected by peat extraction may possess traces of mechanical disturbances undetectable using different analyses based on biota proxy; ii) failure to consider information included in radionuclide date inversions may contribute to the establishment of misleading chronologies. To test these hypotheses, different scenarios of chronology based on high resolution ²¹⁰Pb and ¹⁴C dating from a peat core retrieved from the Puścizna Krauszowska bog (southern Poland) have been analysed. Nowadays, this mire is intensively exploited by humans; however, it contains remains of dome considered undisturbed, from which the core presented in this paper was collected. The set of dates revealed the presence of marked ¹⁴C date inversions (mechanical disturbances) which, if inappropriately interpreted before the age–depth modelling process, may lead to the establishment of misleading chronologies, and thus an incorrect interpretation of biota proxy records, e.g. pollen. Those sections of peat profiles with prominent age inversions and/or strong discrepancies between the peat accumulation rate and bulk density should be rejected from age–depth modelling, even if interpretable chronologies can be obtained.     

Ambient noise energy bursts observation and modeling: Trapping of harmonic structure-soil induced–waves in a topmost sedimentary layer

We study the nature of energy bursts that appeared in the frequency range 3–5 Hz in ambient seismic noise recorded in the Grenoble basin (French Alps) during a seismological array experiment. A close agreement is found between the identified azimuths of such noise bursts with the location of an industrial chimney. In-situ measurements of the chimney dynamic characteristics show a coincidence between the frequency of the first harmonic mode of the chimney and the fundamental frequency of a thin surficial layer that overlay the deep sediment fill. The interaction between the chimney and the surficial layer is then numerically simulated using simple impedance models and two soil profiles. Simulations exhibit a satisfactory agreement with observations and suggest that energy bursts result of inertial structure-soil interaction favored by resonance effects between the first harmonic mode of the structure and the fundamental frequency of the topmost layer.     

Displacement of large-scale open solar magnetic fields from the zone of active longitudes and the heliospheric storm of November 3–10, 2004: 2. “Explosion” of singularity and dynamics of sunspot formation and energy release

A more detailed scenario of one stage (August–November 2004) of the quasibiennial MHD process “Origination ... and dissipation of the four-sector structure of the solar magnetic field” during the decline phase of cycle 23 has been constructed. It has been indicated that the following working hypothesis on the propagation of an MHD disturbance westward (in the direction of solar rotation) and eastward (toward the zone of active longitudes) with the displacement of the large-scale open solar magnetic field (LOSMF) from this zone can be constructed based on LOSMF model representations and data on sunspot formation, flares, active filaments, and coronal ejections as well as on the estimated contribution of sporadic energy release to the flare luminosity and kinetic energy of ejections: (1) The “explosion” of the LOSMF singularity and the formation in the explosion zone of an anemone active region (AR), which produced the satellite sunspot formation that continued west and east of the “anemone,” represented a powerful and energy-intensive source of MHD processes at this stage. (2) This resulted in the origination of two “governing” large-scale MHD processes, which regulated various usual manifestations of solar activity: the fast LOSMF along the neutral line in the solar atmosphere, strongly affecting the zone of active longitudes, and the slow LOSMF in the outer layers of the convection zone. The fronts of these processes were identified by powerful (about 10³¹ erg) coronal ejections. (3) The collision of a wave reflected from the zone of active longitudes with the eastern front of the hydromagnetic impulse of the convection zone resulted in an increase in LOSMF magnetic fluxes, origination of an active sector boundary in the zone of active longitudes, shear-convergent motions, and generation and destabilization of the flare-productive AR 10696 responsible for the heliospheric storm of November 3–10, 2004.