Estimation of the triggering effect energy in relation to model sample failure

Based on the analysis of acoustic emission of model samples under the action of a discretely varied axial load, the minimum energy level of the triggering effect e _tr is estimated as an increment in the accumulated (maximum) potential energy Π that leads to failure. The value of the ratio e _tr/Π determined at this stage of the process averages K _tr = 1 × 10^?8 for samples with various physicomechanical properties. The ratio K _tr is also calculated for levels of the accumulated energy Π _i corresponding to loads k _p of 0.3–1.0. The obtained results are prerequisite for controllable regulation of the energy impact level and deformation regime of a model sample.     

Efficiency of electron acceleration in the Earth’s magnetosphere by whistler mode waves

The efficiency of energetic electron cyclotron acceleration in the Earth’s magnetosphere in different regimes of electron resonant interaction with parallel propagating whistler mode waves of variable frequency, specifically, with chorus ELF-VLF emissions, is considered. The regime of stochastic acceleration, typical of the interaction between particles and noise-like emissions, and particle acceleration in the regime of nonlinear trapping by a quasimonochromatic wave field are discussed. The specific feature of the latter regime consists in its non-diffuse character, i.e., the definite sign of the energy variation depending on the frequency variation in the wave packet. The trapped electron energy becomes higher if frequency increases within an element, which is typical of chorus emissions. For the parameters typical of chorus emissions (the amplitude of a wave magnetic field B _～ = 10² nT, the initial frequency ω ～ 0.3ω _H , and the frequency variation &;Dω ～ 0.15ω _H , where ω _H is the electron gyrofrequency), the energy increase during one act of such an interaction at L = 4?5 exceeds the rms variation in the energy of untrapped electron (during stochastic acceleration) by one-two orders of magnitude. The estimates indicate that a considerable fraction (several tens of percent) of the chorus element energy can be absorbed by electrons accelerated in the trapping regime during a single hop.     

Depositional processes and gas pore pressure in pyroclastic flows: an experimental perspective

The depositional processes and gas pore pressure in pyroclastic flows are investigated through scaled experiments on transient, initially fluidized granular flows. The flow structure consists of a sliding head whose basal velocity decreases backwards from the front velocity (U _f) until onset of deposition occurs, which marks transition to the flow body where the basal deposit grows continuously. The flows propagate in a fluid-inertial regime despite formation of the deposit. Their head generates underpressure proportional to U _f ² whereas their body generates overpressure whose values suggest that pore pressure diffuses during emplacement. Complementary experiments on defluidizing static columns prove that the concept of pore pressure diffusion is relevant for gas-particle mixtures and allow characterization of the diffusion timescale (t _d) as a function of the material properties. Initial material expansion increases the diffusion time compared with the nonexpanded state, suggesting that pore pressure is self-generated during compaction. Application to pyroclastic flows gives minimum diffusion timescales of seconds to tens of minutes, depending principally on the flow height and permeability. This study also helps to reconcile the concepts of en masse and progressive deposition of pyroclastic flow units or discrete pulses. Onset of deposition, whose causes deserve further investigation, is the most critical parameter for determining the structure of the deposits. Even if sedimentation is fundamentally continuous, it is proposed that late onset of deposition and rapid aggradation in relatively thin flows can generate deposits that are almost snapshots of the flow structure. In this context, deposition can be considered as occurring en masse, though not strictly instantaneously.     

Effects of relative motion between gas and liquid on 1-dimensional steady flow in silicic volcanic conduits: 2. Origin of diversity of eruption styles

We investigate the origin of diversity of eruption styles in silicic volcanoes on the basis of a 1-dimensional steady conduit flow model that considers vertical relative motion between gas and liquid (i.e., vertical gas escape). The relationship between the assemblage of steady-state solutions in the conduit flow model and magma properties or geological conditions is expressed by a regime map in the parameter space of the ratio of liquid-wall friction force to liquid–gas interaction force (non-dimensional number ε), and a normalized conduit length Λ. The regime map developed in the companion paper shows that when ε is smaller than a critical value ε_cr, a solution of explosive eruption exists for a wide range of Λ, whereas an effusive solution exists only when Λ ~ 1. On the other hand, when ε > ε_cr, an effusive solution exists for a wide range of Λ. Diversity of eruption styles observed in nature is explained by the change in ε accompanied by the change in magma viscosity during magma ascent. As magma ascends, the magma viscosity increases because of gas exsolution and crystallization, leading to the increase in ε. For the viscosity of hydrous silicic magma at magma chamber, ε is estimated to be smaller than ε_cr, indicating that an explosive solution exists for wide ranges of geological parameters. When magma flow rate is small, the viscosity of silicic magma drastically increases because of extensive crystallization at a shallow level in the conduit. In this case, ε can be greater than ε_cr; as a result, a stable effusive solution co-exists with an explosive solution.     

Oscillations of the Earth with periods of 9–57 min in the background seismic process and the energy flux direction in the range of the free oscillation 0 S 2

Retrospective data of monitoring under conditions of low seismic activity are used to identify free oscillations of the Earth, including the fundamental mode, the oscillation with a central period of 54 min (₀ S ₂ ^m multiplet), split into five lines with azimuthal numbers m = ?2, ?1, 0, 1, 2. It is shown that some lines of this oscillation are also recorded in atmospheric pressure variation spectra and group in ensembles of observations around frequencies predicted by the ₀ S ₂ ^m splitting theory. This phenomenon is discovered in data recorded both synchronously and in different time intervals. A causal relationship involved in the oscillation under study is determined on the basis of the examination of the direction of the acoustic energy flux. The energy flux in the region of the ₀ S ₂ ^m multiplet is shown to be directed from the Earth toward the atmosphere. This suggests that deep processes in the Earth are capable of exciting its upper shells.     

Ab initio study of the composition dependence of the pressure-induced spin crossover in perovskite (Mg1 − x,Fex)SiO3

We present ab initio calculations of the zero-temperature iron high- to low-spin crossover in (Mg_1 ? xFe_x)SiO₃ perovskite at pressures relevant to Earth's lower mantle. Equations of state are fit for a range of compositions and used to predict the Fe spin transition pressure and associated changes in volume and bulk modulus. We predict a dramatic decrease in transition pressure as Fe concentration increases. This trend is contrary to that seen in ferropericlase, and suggests the energetics for spin crossover is highly dependent on the structural environment of Fe. Both Local Density Approximation (LDA) and Generalized Gradient Approximation (GGA) exchange-correlation methods are used, and both methods reproduce the same compositional trends. However, GGA gives a significantly higher transition pressure than LDA. The spin transition is made easier by the decreasing spin-flip energy with pressure but is also driven by the change in volume from high to low spin. Volume trends show that high-spin Fe²⁺ is larger than Mg²⁺ even under pressure, but low-spin Fe²⁺ is smaller at ambient conditions and approximately the same size as Mg²⁺ under high pressure, indicating that low-spin Fe²⁺ is less compressible than high-spin Fe²⁺. We find large changes between high- and low-spin in the slope of volume with Fe concentration. Although these changes are small in absolute magnitude for small Fe content, they are still important when measured per Fe and could be relevant for calculating partitioning coefficients in the lower mantle.     

T-phase Observations in Northern California: Acoustic to Seismic Coupling at a Weakly Elastic Boundary

—?Plans for a hydroacoustic network intended to monitor compliance with the CTBT call for the inclusion of five T-phase stations situated at optimal locations for the detection of seismic phases converted from ocean-borne T phases. We examine factors affecting the sensitivity of land-based stations to the seismic T phase. The acoustic to seismic coupling phenomenon is described by upslope propagation of an acoustic ray impinging at a sloping elastic wedge. We examine acoustic to seismic coupling characteristics for two cases; the first in which the shear velocity of the bottom is greater than the compressional velocity of the fluid (i.e., v _p > v _s > v _w ), the second is a weakly elastic solid in which v _s << v _w < v _p . The former is representative of velocities in solid rock, which might be encountered at volcanic islands; the latter is representative of marine sediments. For the case where v _s > v _w , we show that acoustic energy couples primarily to shear wave energy, except at very high slope angles. We show that the weakly elastic solid (i.e., v _s << v _w ) behaves nearly like a fluid bottom, with acoustic energy coupling to both P and S waves even at low slope angles.¶We examine converted T-wave arrivals at northern California seismic stations for two event clusters; one a series of earthquakes near the Hawaiian Islands, the other a series of nuclear tests conducted near the Tuamoto archipelago. Each cluster yielded characteristic arrivals at each station which were consistent from event to event within a cluster, but differed between clusters. The seismic T-phases consisted of both P- and S-wave arrivals, consistent with the conversion of acoustic to seismic energy at a gently sloping sediment-covered seafloor. In general, the amplitudes of the seismic T phases were highest for stations nearest the continental slope, where seafloor slopes are greatest, however noise levels decrease rapidly with increasing distance from the coastline, so that T-wave arrivals were observable at distances reaching several hundred kilometers from the coast. Signal-to-noise levels at the seismic stations are lower over the entire frequency spectrum than at the Pt. Sur hydrophone nearby, and decrease more rapidly with increasing frequency, particularly for stations furthest from the continental slope.     

Ion microprobe studies of water in silicate melts: Temperature-dependent water diffusion in obsidian

The temperature dependence of water diffusivity in rhyolite melts over the range 650–950°C and [P_T(H₂O] = 700 bars is evaluated from water concentration-distance profiles measured in glass with an ion microprobe. Diffusivities are exponentially dependent on concentration over this temperature range and vary from about 10^?8 cm²/s at 650°C to about 10^?7 cm²/s at 950°C at 2 wt.% water. Water solubility also varies with temperature at a rate of ?0.14 wt. per 100°C increase. The avtivation energy (E_a) appears to be constant at 19 ± 1kal/mole for 1, 2,and 3 wt.% H₂O. Comparison of these data with results for cation diffusion indicates that this value is a minimum E_a for diffusion of any species in a rhyolite melt.Compensation plots of log₁₀D₀ (the frequency factor) versus E_a indicate that hydrous rhyolite melts follow the same trend as anhydrous basalts. D₀ increases for H₂O and Ca²⁺ [1] as E_a decreases. This suggests that these molecules may diffuse by different mechanisms than do monovalent cations, and that hydration of the melt affects diffusion of Ca²⁺ and H₂O differently than it does monovalent cation diffusion. The results imply that dramatic increases in cation diffusivities by hydration [1] may occur with additions of less than 1 wt.% H₂O.     

Hot springs mediate spatial exchange of heat and mass in the enclosed sediment domain: A stability perspective

In many natural environments, such as in underwater hot springs and hydrothermal vents, thermal gradients are accompanied with changes in the concentration of chemical compounds transported to the seawater, causing the so-called double-diffusive, mixed convection. To study the physical scenarios in such systems, a vertical channel filled with a porous medium saturated with saline water is considered. The motion in the sediment-filled channel is induced by two buoyancy forces and an external pressure gradient, similar to the situation in a vent with an upward flow direction. The fluid flow has been modeled by an extended Darcy model, and the flow instability mechanisms have been studied numerically. The linear stability analysis is performed considering a wide range of Darcy number (Da = 10⁻⁵ -10⁻⁸). The instability boundary curve showed three distinct dynamic regimes: (i) Rayleigh-Taylor (R-T), (ii) log-log non-linear variation, and (iii) log-log linear variation. The domain of different regimes were sensitive to external pressure gradient as well as permeability. Similar to cross-diffusive natural convection in pure viscous fluids, a linear relationship between logarithmic absolute values of critical thermal Rayleigh number (∣Ra_T∣) and solute Rayleigh number (Ra_C) is found in the third regime. Based on the permeability, for any solute Rayleigh number (Ra_C), there existed a minimum value of Reynolds number (Re), below which R-T type of instability appeared. Above this minimum value, the instability was due to two buoyancy forces, known as buoyant instability. Simulations of secondary flow via energy analysis demonstrated the development of complex dynamics at the critical state in all three regimes characterized by transition of multi to uni-cellular structures and vice verse.     

A new approach to assessing the water footprint of hydroelectric power based on allocation of water footprints among reservoir ecosystem services

Hydroelectric power is an important energy source to meet the growing demand for energy, and large amounts of water are consumed to generate this energy. Previous studies often assumed that the water footprint of hydroelectric power equaled the reservoir’s water footprint, but failed to allocate the reservoir water footprint among the many beneficiaries; dealing with this allocation remains a challenge. In this study, we developed a new approach to quantify the water footprint of hydroelectric power (WF_h) by separating it from the reservoir water footprint (WF) using an allocation coefficient (η_h) based on the ratio of the benefits from hydroelectric power to the total ecosystem service benefits. We used this approach in a case study of the Three Gorges Reservoir, the world’s largest reservoir, which provides multiple ecosystem services. We found large differences between the WF_h and the water footprint of per unit of hydroelectric production (PWF_h) calculated using η_h and those calculated without this factor. From 2003 to 2012, η_h decreased sharply (from 0.76 in 2005 to 0.41 in 2012), which was due to the fact that large increases in the value of non-energy ecosystem services, and particularly flood control. In 2009, flood control replaced hydroelectricity as the largest ecosystem service of water from the Three Gorges Reservoir. Using our approach, WF_h and PWF_h averaged 331.0 × 10⁶ m³ and 1.5 m³ GJ⁻¹, respectively. However, these values would almost double without allocating water footprints among different reservoir ecosystem services. Thus, previous studies have overestimated the WF_h and PWF_h of reservoirs, especially for reservoirs that serve multiple purposes. Thus, the allocation coefficient should not be ignored when calculating the WF of a product or service.     

Thermal diffusivity of rhyolitic glasses and melts: effects of temperature, crystals and dissolved water

Thermal diffusivity (D) was measured using laser-flash analysis on pristine and remelted obsidian samples from Mono Craters, California. These high-silica rhyolites contain between 0.013 and 1.10?wt% H₂O and 0 to 2?vol% crystallites. At room temperature, D _glass varies from 0.63 to 0.68?mm²?s^?1, with more crystalline samples having higher D. As T increases, D _glass decreases, approaching a constant value of ??0.55?mm²?s^?1 near 700?K. The glass data are fit with a simple model as an exponential function of temperature and a linear function of crystallinity. Dissolved water contents up to 1.1?wt% have no statistically significant effect on the thermal diffusivity of the glass. Upon crossing the glass transition, D decreases rapidly near ??1,000?K for the hydrous melts and ??1,200?K for anhydrous melts. Rhyolitic melts have a D _melt of ??0.51?mm²?s^?1. Thermal conductivity (k?=?D·??·C _P) of rhyolitic glass and melt increases slightly with T because heat capacity (C _P) increases with T more strongly than density (??) and D decrease. The thermal conductivity of rhyolitic melts is ??1.5?W?m^?1?K^?1, and should vary little over the likely range of magmatic temperatures and water contents. These values of D and k are similar to those of major crustal rock types and granitic protoliths at magmatic temperatures, suggesting that changes in thermal properties accompanying partial melting of the crust should be relatively minor. Numerical models of shallow rhyolite intrusions indicate that the key difference in thermal history between bodies that quench to obsidian, and those that crystallize, results from the release of latent heat of crystallization. Latent heat release enables bodies that crystallize to remain at high temperatures for much longer times and cool more slowly than glassy bodies. The time to solidification is similar in both cases, however, because solidification requires cooling through the glass transition in the first case, and cooling only to the solidus in the second.     

Effects of relative motion between gas and liquid on 1-dimensional steady flow in silicic volcanic conduits: 1. An analytical method

We investigate the effects of vertical relative motion between gas and liquid on eruption styles by formulating a model for 1-dimensional steady flow in volcanic conduits. As magma ascends and decompresses, volatiles exsolve and volume fraction of gas increases. As a result, magma fragmentation occurs and the flow changes from bubbly flow to gas-pyroclast flow. In our model, a transitional region (‘permeable flow region’) is introduced between the bubbly flow region and the gas-pyroclast flow region. In this region, both the gas and the liquid are continuous phases, allowing the efficient vertical escape of gas through the permeable structure. We describe the features of conduit flow with relative motion of gas and liquid using non-dimensional numbers α, γ and ε. The parameter α represents the ratio of effects of wall friction to gravitational load, and is proportional to magma flow rate. The parameter γ represents the degree of decompression for the gas-pyroclast flow to reach the sound velocity at α = 1, and is proportional to r_c²/μ for given magma temperature and initial volatile content, where r_c is conduit radius and μ is liquid viscosity. The parameter ε is defined as the ratio of liquid–wall friction force to liquid–gas interaction force in the permeable flow region, and represents the efficiency of gas escape from magma. The values of γ and ε are determined only by magma properties and geological conditions such as liquid viscosity, magma permeability and conduit radius. We formulate a 1-dimensional steady-state conduit flow model to find non-dimensional magma flow rate α as a function of magma properties and geological conditions (e.g., γ and ε) under given boundary conditions. When the relative motion is taken into account with the assumption that magma fragmentation occurs when the gas volume fraction reaches some critical values, the pressure at the fragmentation level (P_f) decreases as the magma flow rate (α) decreases or the efficiency of gas escape (ε) increases, because gas escape suppresses the increase in the gas volume fraction accompanying magma ascent. When ε is so large that P_f is below the atmospheric pressure (P_a), the flow reaches the vent before fragmentation at low α. On the other hand, when ε is so small that P_f is greater than P_a, the flow reaches the vent after fragmentation at high α. These steady-state solutions of the flow at low and high α correspond to effusive and explosive eruptions, respectively. We present a graphical method to systematically find α. On the basis of the graphical method, a simple regime map showing the relationship between the assemblage of the solutions of conduit flow and the magma properties or the geological conditions is obtained.     

A parametric representation of Kamchatka seismicity over time

This paper presents a set of seismicity parameters that are estimated at the Kamchatka Branch of the Geophysical Service, Russian Academy of Sciences based on the regional catalog data with the purpose of routine monitoring of the current seismic situation in the region. The focus is on the identification of changes in the seismic regime (seismic quiescences and seismicity increases) in earth volumes adjacent to the maturing rupture zone of a large earthquake. The techniques we use include estimation of the seismicity level for the region using the SOUS’09 scale; calculation of the variations in the slope of the recurrence relation, identification of statistically significant anomalies in the slope using the Z test, and calculation of the seismic activity A ₁₀; monitoring the RTL parameter and variations in the area of seismogenic ruptures; using the Z test to detect areas of statistically significant decreases in the rate of seismicity; and identification of earthquake clusters. We furnish examples of such anomalies in these seismicity parameters prior to large earthquakes in Kamchatka.     

A bed‐load transport model for rough turbulent open‐channel flows on plane beds

Data from flume studies are used to develop a model for predicting bed‐load transport rates in rough turbulent two‐dimensional open‐channel flows moving well sorted non‐cohesive sediments over plane mobile beds. The object is not to predict transport rates in natural channel flows but rather to provide a standard against which measured bed‐load transport rates influenced by factors such as bed forms, bed armouring, or limited sediment availability may be compared in order to assess the impact of these factors on bed‐load transport rates. The model is based on a revised version of Bagnold's basic energy equation i_bs_b = e_bω, where i_b is the immersed bed‐load transport rate, ω is flow power per unit area, e_b is the efficiency coefficient, and s_b is the stress coefficient defined as the ratio of the tangential bed shear stress caused by grain collisions and fluid drag to the immersed weight of the bed load. Expressions are developed for s_b and e_b in terms of G, a normalized measure of sediment transport stage, and these expressions are substituted into the revised energy equation to obtain the bed‐load transport equation i_b = ω G ^3·4. This equation applies regardless of the mode of bed‐load transport (i.e. saltation or sheet flow) and reduces to i_b = ω where G approaches 1 in the sheet‐flow regime. That i_b = ω does not mean that all the available power is dissipated in transporting the bed load. Rather, it reflects the fact that i_b is a transport rate that must be multiplied by s_b to become a work rate before it can be compared with ω. It follows that the proportion of ω that is dissipated in the transport of bed load is i_bs_b/ω, which is approximately 0·6 when i_b = ω. It is suggested that this remarkably high transport efficiency is achieved in sheet flow (1) because the ratio of grain‐to‐grain to grain‐to‐bed collisions increases with bed shear stress, and (2) because on average much more momentum is lost in a grain‐to‐bed collision than in a grain‐to‐grain one. Copyright © 2006 John Wiley & Sons, Ltd.     

Shear viscosity of rhyolite-vapor emulsions at magmatic temperatures by concentric cylinder rheometry

The viscosity of natural rhyolitic melt from Lipari, Aeolian Islands and melt-bubble emulsions (30–50 vol% porosity) generated from Lipari rhyolite have been measured in a concentric cylinder rheometer at temperatures and shear rates in the range 925–1150°C and 10⁻³–10^−1.2 s⁻¹, respectively, in order to better understand the dependence of emulsion shear viscosity on temperature and shear rate in natural systems. Bubble-free melt exhibits Newtonian–Arrhenian behavior in the temperature range 950–1150°C with an activation energy of 395±30 kJ/mol; the shear viscosity is given by log η_m=−8.320+20624/T. Suspensions were prepared from natural rhyolite glass to which small amounts of Na₂SO₄ were added as a ‘foaming agent’. Reasonably homogeneous magmatic mixtures with an approximate log-normal distribution of bubbles were generated by this technique. Suspension viscosity varied from 10^6.1 to 10^8.37 Pa s and systematically correlates with temperature and porosity in the shear stress range (10^4.26–10^5.46 Pa) of the experiments. The viscosity of melt-bubble emulsions is described in terms of the relative viscosity, η_r=η_e/η_m where η_e is the emulsion viscosity and η_m is the viscosity of melt of the same composition and temperature. The dependence of relative viscosity on porosity for magmatic emulsions depends on the magnitude of the capillary number Ca≡G/(σr_b⁻¹η_m⁻¹), the ratio of viscous forces acting to deform bubbles to interfacial forces resisting bubble deformation. For inviscid bubbles in magmatic flows three regimes may be identified. For Ca<0.1, bubbles are nearly spherical and relative viscosity is an increasing function of porosity. For dilute systems, η_r=1+φ given by the classical result of Taylor [Proc. R. Soc. London A 138 (1932) 41–48]. For Ca in the range 0.1<Ca<10, emulsions behave as power law fluids and the relative viscosity depends on shear rate (or Ca) as well as porosity. At high Ca (Ca>10) an asymptotic regime is reached in which relative viscosity decreases with increasing porosity and is independent of Ca. Our experiments were carried out for 30<Ca<925 in order to quantify the maximal effect of bubbles in reducing the viscosity of magmatic emulsions relative to single-phase melt at identical conditions of shear rate and temperature. The viscosity of a 50 vol% emulsion is a factor of five smaller than that of melt alone. Rheometric measurements obtained in this study are useful in constraining models of magma transport and volcanic eruption mechanics relevant to transport of volatile-saturated magma in the crust and upper mantle.     

Regression analysis of reported earthquake precursors. I. Presentation of data

Around 700 reported precursors of about 350 earthquakes, including the negative observations, have been compiled in 11 categories with 31 subdivisions. The data base is subjected to an initial sorting and screening by imposing three restrictions on the ranges of main shock magnitude (M≥4.0), precursory time (t≤20 years), and the epicentral distance of observation points (X _m≤4.10^0.3M ). Of the 31 subcategories of precursory phenomena, 18 with 9 data points or more are independently studied by regressing their precursory times against magnitude. The preliminary results tend to classify the precursors into three groups:

1. The precursors which show weak or no correlation between time and the magnitude of the eventual main shock. Examples of this group are foreshocks and precursory tilt.
2. The precursors which show clear scaling with magnitude. These include seismic velocity ratio (V _p/V_s), travel time delay, duration of seismic quiescence, and, to some degree, the variation ofb-value, and anomalous seismicity.
3. The precursors which display clustering of precursory times around a mean value, which differs for different precursors from a few hours to a few years. Examples include the conductivity rate, geoelectric current and potential, strain, water well level, geochemical anomalies, change of focal mechanism, and the enhancement of seismicity reported only for larger earthquakes. Some of the precursors in this category, such as leveling changes and the occurrence of microseismicity, show bimodal patterns of precursory times and may partially be coseismic.

In addition, each category with a sufficient number of reported estimates of distance and signal amplitude is subjected to multiple linear regression. The usefulness of these regressions at this stage appears to be limited to specifying which of the parameters shows a more significant correlation. Standard deviations of residuals of precursory time against magnitude are generally reduced when observation distance enters as a second independent variable. The effect is more pronounced for water well level and conductivity rate changes. While a substantial portion of the data seem to suffer from personal bias, hence should be regarded as noise, the observations of a number of strain sensitive phenomena such as strain, water well level, and conductivity rate changes, appear to be internally more consistent. For instance, their precursory times suggest a scaling relationship with the strain energy surface density associated with the main shock. The scaling is not identical for all three phenomena so that they may constitute the imminent, short- and intermediate-term manifestation of the same process, i.e. strain loading, respectively.     

The effect of wave-induced turbulence on intertidal mudflats: Impact of boat traffic and wind

Semi-diurnal and fortnightly surveys were carried out to quantify the effects of wind- and navigation-induced high-energy events on bed sediments above intertidal mudflats. The mudflats are located in the upper fluvial part (Oissel mudflat) and at the mouth (Vasière Nord mudflat) of the macrotidal Seine estuary. Instantaneous flow velocities and mudflat bed elevation were measured at a high frequency and high resolution with an acoustic doppler velocimeter (ADV) and an ALTUS altimeter, respectively. Suspended particulate matter concentrations were estimated by calibrating the ADV acoustic backscattered intensity with bed sediments collected at the study sites. Turbulent bed shear stress values were estimated by the turbulent kinetic energy method, using velocity variances filtered from the wave contribution. Wave shear stress and maximum wave–current shear stress values were calculated with the wave–current interaction (WCI) model, which is based on the bed roughness length, wave orbital velocities and the wave period (T_S). In the fluvial part of the estuary, boat passages occurred unevenly during the surveys and were characterized by long waves (T_S>50 s) induced by the drawdown effect and by short boat-waves (T_S<10 s). Boat waves generated large bottom shear stress values of 0.5 N m⁻² for 2–5 min periods and, in burst of several seconds, larger bottom shear stress values up to 1 N m⁻². At the mouth of the estuary, west south-west wind events generated short waves (T_S<10 s) of H_S values ranging from 0.1 to 0.3 m. In shallow-water environment (water depth <1.5 m), these waves produced bottom shear stress values between 1 and 2 N m⁻². Wave–current shear stress values are one order of magnitude larger than the current-induced shear stress and indicate that navigation and wind are the dominant hydrodynamic forcing parameters above the two mudflats. Bed elevation and SPM concentration time series showed that these high energy events induced erosion processes of up to several centimetres. Critical erosion shear stress (τ_ce) values were determined from the SPM concentration and bed elevation measurements. Rough τ_ce values were found above 0.2 N m⁻² for the Oissel mudflat and about 1 N m⁻² for the Vasière Nord mudflat.     

Numerical modeling of bedload and suspended load contributions to morphological evolution of the Seine Estuary (France)

This numerical modeling study (i) assesses the influence of the sediment erosion process on the sediment dynamics and subsequent morphological changes of a mixed-sediment environment, the macrotidal Seine estuary, when non-cohesive particles are dominant within bed mixtures (non-cohesive regime), and (ii) investigates respective contributions of bedload and suspended load in these dynamics. A three dimensional (3D) process-based morphodynamic model was set up and run under realistic forcings (including tide, waves, wind, and river discharge) during a 1-year period. Applying erosion homogeneously to bed sediment in the non-cohesive regime, i.e., average erosion parameters in the erosion law (especially the erodibility parameter, E₀), leads to higher resuspension of fine sediment due to the presence of coarser fractions within mixtures, compared to the case of an independent treatment of erosion for each sediment class. This results in more pronounced horizontal sediment flux (two-fold increase for sand, +30% for mud) and erosion/deposition patterns (up to a two-fold increase in erosion over shoals, generally associated with some coarsening of bed sediment). Compared to observed bathymetric changes, more relevant erosion/deposition patterns are derived from the model when independent resuspension fluxes are considered in the non-cohesive regime. These results suggest that this kind of approach may be more relevant when local grain-size distributions become heterogeneous and multimodal for non-cohesive particles. Bedload transport appears to be a non-dominant but significant contributor to the sediment dynamics of the Seine Estuary mouth. The residual bedload flux represents, on average, between 17 and 38% of the suspended sand flux, its contribution generally increasing when bed sediment becomes coarser (can become dominant at specific locations). The average orientation of residual fluxes and erosion/deposition patterns caused by bedload generally follow those resulting from suspended sediment dynamics. Sediment mass budgets cumulated over the simulated year reveal a relative contribution of bedload to total mass budgets around 25% over large erosion areas of shoals, which can even become higher in sedimentation zones. However, bedload-induced dynamics can locally differ from the dynamics related to suspended load, resulting in specific residual transport, erosion/deposition patterns, and changes in seabed nature.     

Effective energy loss per electron-ion pair in proton aurora

Effective energy loss per electron-ion pair produced, <xi>(E ₀), as a function of a particle’s initial energy has been obtained for proton transport in the atmosphere. The influence of some transport parameters on the shape of <xi>(E ₀) has been studied. Comparisons with the case of electron transport and with other results were made. It has been shown that: 1. for E ₀>1 keV, <xi>(E ₀) varies within the range 30–36 eV; 2. as E ₀ increases the value of <xi>(E ₀) tries to attain an asymptotic value that is the same as for electrons (\approx35 eV); 3. <xi>(E ₀) strongly depends on the average energy of secondary electrons, but the energy distribution of secondary electrons is not as important. The range of possible changes in <xi>(E ₀) associated with discrepancies in cross sections has been obtained.