Propagation of slip along frictional surfaces

The propagation of slip along a pre-existing frictional plane is formulated for the faults containing interstitial fluid. When normal and shear stresses satisfy the effective frictional law, a frictional sliding occurs stably or unstably, depending on the inhomogeneities of the surface. For a slip to rapidly sweep the surface, the pore pressure of fluid must exceed the critical value that is related to the physical or geometrical irregularities of the surface. If this condition fails, stable sliding is expected, analogous to a seismic fault creep. This prediction makes the role of water in seismic faulting more clear.     

Analytical solution of transient flow in a sloping soil layer with recharge

Abstract

An analytical solution of planar flow in a sloping soil layer described by the linearized extended Boussinesq equation is presented. The solution consists of the sum of steady-state and transient-series solutions, the latter in a separation-of-variables form, and can satisfy an arbitrary initial condition via collocation; this feature reduces the number of series terms, making the solution efficient. Key parameter is the dimensionless linearization depth η _o (R), R being the dimensionless recharge. The variable η _o (R), not the slope, characterizes the flow as kinematic or diffusive, and R ≈ 0.2 demarcates the two regimes. The transient series converges rapidly for large η _o (large R, near-diffusive flow) and slowly as η _o → 0 (kinematic flow). The quasi-steady (QS) state method of Verhoest & Troch is also analysed and it is shown that the QS depth profiles approximate the transient ones well, only if Δt exceeds a system-dependent transition time between flow states (possibly >>1 day). In an application example for a 30-day recharge series, the QS solution differs from the transient one by as much as 20% (RMSE = 15%), does not track recharge changes as well and fails to conserve mass.     

Boundary layer models of ocean floor spreading

Summary The equations of conservations of momentum and energy scaled with the characteristic values of the mantle indicate the presence of the upper boundary layer to produce the estimated rate of the ocean floor spreading by convection and the importance of the frictional heating. The depth of the upper boundary layer can be estimated from the balance of the viscous force with the horizontal pressure gradient at the sea floor. It is of the orders of 100 km and becomes deeper for the Pacific than for the Atlantic Ocean and also with frictional heating than without it. The frictional heating increases the surface heat flow of the heat conduction by ten to twenty percent for the Pacific Ocean but only by a few percent for the Atlantic Ocean. The similarity solutions are determined for the temperature and horizontal velocity in the upper boundary layer. These solutions are expressed in power series of the variabley x ^–n , wherex, y, andn are horizontal and vertical coordinates and numerical constant, respectively. Both temperature and horizontal velocity within the boundary layer are higher for the Pacific than for the Atlantic Ocean. When a larger viscosity is applied, it causes the increase of horizontal velocity below the surface because of the surface boundary conditions of the finite velocity and of vanishment of the velocity shear. The higher horizontal velocity generates higher temperature because it advects hotter material from the mid-ocean ridge site. The direct effect of frictional heating on the temperature distribution of the similarity solution is almost negligible, since the shear zone is deep and near the lower boundary of the upper boundary layer. In the similarity solution, the surface heat flow which is increased by the frictional heating is given as the boundary value. The effect of the frictional heating is important below the mid-ocean ridge.     

Explicit data-driven models for prediction of pressure fluctuations occur during turbulent flows on sloping channels

Stochastic Environmental Research and Risk Assessment - Pressure fluctuations are among the favorite topics for hydraulic engineers due to their critical role in the design and safe operation of...     

Unsteady boundary layer flows induced by a wavy plate

Summary Unsteady boundary layer flows generated in an incompressible, homogeneous, nonrotating viscous fluid bounded by a rigid wavy plate are studied theoretically. The Laplace transform method is employed to obtain exact solutions of the unsteady boundary layer equations in a wavy plate configuration. The structures of the unsteady velocity distribution and the associated boundary layers are determined explicitly and several particular solutions are recovered as special cases of this analysis. The physical interpretation of the mathematical results are examined.     

Two-dimensional scattering of SH waves in a soil layer underlain with a sloping bedrock

A two-dimensional analysis is applied to examine the effect that a sloping bedrock half-space has on the amplification of an anti-plane shear wave. The direct boundary integral equation method is used for the two-dimensional analysis. The particular soil–rock configuration investigated includes a homogeneous soil layer underlain by a sloping rock half-space. The rock half-space dips for a horizontal distance L and then becomes horizontal so that the overlying soil layer has a thickness H that remains constant from this point to infinity. The materials in the soil–rock configuration are considered viscoelastic except in the rock half-space below soil layer thickness H, which is considered elastic. This limitation in damping is due to the correction used for the truncation of the half-space boundary. Four cases are used to study the relationship between rock slope and surface displacement, vertical, 1:2, 1:4, 1:8. Surface displacements are determined for each of these cases for half-space incidence angles of 90, 75, and 60°. To allow for applicability to a wide range of problems, results are determined as a function of dimensionless parameters. In addition, solutions from a one-dimensional analysis are compared with the results of the two-dimensional analysis to establish limits outside of which a one-dimensional analysis suffices.     

Structure of numerically simulated katabatic and anabatic flows along steep slopes

Direct numerical simulation (DNS) is applied to investigate properties of katabatic and anabatic flows along thermally perturbed (in terms of surface buoyancy flux) sloping surfaces in the absence of rotation. Numerical experiments are conducted for homogeneous surface forcings over infinite planar slopes. The simulated flows are the turbulent analogs of the Prandtl (1942) one-dimensional laminar slope flow. The simulated flows achieve quasi-steady periodic regimes at large times, with turbulent fluctuations being modified by persistent low-frequency oscillatory motions with frequency equal to the product of the ambient buoyancy frequency and the sine of the slope angle. These oscillatory wave-type motions result from interactions between turbulence and ambient stable stratification despite the temporal constancy of the surface buoyant forcing. The structure of the mean-flow fields and turbulence statistics in simulated slope flows is analyzed. An integral dynamic similarity constraint for steady slope/wall flows forced by surface buoyancy flux is derived and quantitatively verified against the DNS data.     

An asperity model to simulate rupture along heterogeneous fault surfaces

A model has been developed to simulate the statistical and mechanical nature of rupture on a heterogeneous strike-slip fault. The model is based on the progressive failure of circular asperities of varying sizes and strengths along a fault plane subjected to a constant far-field shear displacement rate. The basis of the model is a deformation and stress intensity factory solution for a single circular asperity under a unidirectional shear stress. The individual asperities are unified through the fault stiffness and the far-field stress and displacement. During fault deformation asperities can fail and reheal, resulting in changes in the local stresses in the asperities, stress drops, and changes in the stiffness of the fault. Depending on how the stress is redistributed following asperity failure and on the strenghts of the neighboring asperities an earthquake event can be the failure of one or more asperities. Following an earthquake event seismic source parameters such as the stress drop, energy change, and moment magnitude are calculated. Results from the model show a very realistic pattern of earthquake rupture, with reasonable source parameters, the proper magnitude-frequency behavior, and the development of characteristic earthquakes. Also the progression ofb-values in the model gives some insight into the phenomenon of self-organized criticality.     

Boundary layer characteristics and turbulent exchange mechanisms in highly complex terrain

The Mesoscale Alpine Programme’s Riviera project investigated the turbulence structure and related exchange processes in an Alpine valley by combining a detailed experimental campaign with high-resolution numerical modelling. The present contribution reviews published material on the Riviera Valley’s boundary layer structure and discusses new material on the near-surface turbulence structure. The general conclusion of the project is that despite the large spatial variability of turbulence characteristics and the crucial influence of topography at all scales, the physical processes can accurately be understood and modelled. Nevertheless, many of the “text book characteristics” like the interaction between the valley and slope wind systems or the erosion of the nocturnal valley inversion need reconsideration, at least for small non-ideal valleys like the Riviera Valley. The project has identified new areas of research such as post-processing methods for turbulence variables in complex terrain and new approaches for the surface energy balance when advection is non-negligible. The exchange of moisture and heat between the valley atmosphere and the free troposphere is dominated by local “secondary” circulations due to the curvature of the valley axis. Because many curved valleys exist, and operational models still have rather poor resolution, parameterization of these processes may be required.     

Experimental study of the Ekman layer instability in steady or oscillating flows

To investigate the instabilities of steady and oscillating Ekman layers, an 8 m×2 m horizontal plate was moved at controlled speed in homogeneous water at rest in solid body rotation in the “Coriolis” 13 m diameter rotating tank. For a steady Ekman layer two distinct wave types were found, in agreement with previous experimental or numerical studies. Type I was stationary, was oriented positively with respect to the flow direction and had a wavelength of about 10 times the Ekman layer thickness. Type II was oriented negatively with respect to the flow direction and had a wavelength which was more than 20 times the Ekman layer thickness and a phase-speed between 0.3 and 0.5 the forcing interior velocity. The growth rates of both type I and type II waves for various Reynolds numbers Re (computed with the Ekman layer thickness) were estimated and their Re-variations qualitatively agree with previous numerical results. For an oscillating Ekman layer, experimental results depended strongly on Ro_t, the temporal Rossby number: only when Ro_t<1 was it possible to observe either type I or type II instabilities. Moreover, for all Ro_t and average to high Re, there was a noticeable upward turbulent transport occurring during each cycle between the flow maximum and the flow reversal. Such an upward turbulent transport is consistent with observations in the English Channel where maximum upward benthic movements and maximum turbidity were recorded at the flow reversal, hence Ekman layer instabilities and transition to turbulence are likely to occur in shallow tidal seas where they may be relevant for sediment resuspension and transport as well as for some biological processes.     

Accuracy of kinematic wave and diffusion wave approximations for space-independent flows on infiltrating surfaces

Error equations for the kinematic wave and diffusion wave approximations were derived under simplified conditions for space-independent flows occurring on infiltrating planes or channels. These equations specify error as a function of time in the flow hydrograph. The kinematic wave, diffusion wave and dynamic wave solutions were parameterized through a dimensionless parameter γ which is dependent on the initial conditions. This parameter reflects the effect of initial flow depth, channel bed slope, lateral inflow and channel roughness when the initial condition is non-vanishing; it reflects the effect of bed slope, channel roughness, lateral inflow and infiltration when the initial condition is vanishing. The error equations were found to be the Riccati equation.     

Parameterization of the logarithmic layer of double-averaged streamwise velocity profiles in gravel-bed river flows

The logarithmic layer of double-averaged (in time and space) streamwise velocity profiles obtained from field measurements made in the Swiss rivers, Venoge and Chamberonne is parameterized and discussed. Velocity measurements were made using a 3D Acoustic Doppler Velocity Profiler. Both riverbeds are hydraulically rough, composed of coarse gravel, with relative submergences (h/D₅₀) of 5.25 and 5.96, respectively. From the observations, the flow may be divided into three different layers: a roughness layer near the bed, an equivalent logarithmic layer and a surface or outer layer. It was found that a logarithmic law can describe the double-averaged profiles in the layer 0.30 < z/h < 0.75. The parameterization of the logarithmic law is discussed. Special emphasis is given to the geometric parameters roughness and zero-displacement heights and to the equivalent von Karman constant.     

On similarity solution of an unsteady laminar boundary layer along a flat plate

Summary A unified analysis has been made to obtain all possible similarity solutions of the steady and unsteady, forced flow, inside a boundary layer along a flat plate. Though previously, attempts were made to obtain similarity solutions of a steady boundary layer flow neglecting viscous dissipation term in the energy conservation equation but the treatments were not complete. Here we have taken account of the viscous dissipation term. In the steady case it has been shown that for a similarity solution of both velocity and temperature, there should be a relation between the undisturbed flow outside the boundary layer and the temperature of the plate. It has been shown that the similarity solution exists in the unsteady case if we neglet the viscous dissipation term in the energy equation.     

Boundary layer dissipation and the nonlinear interaction of equatorial baroclinic and barotropic rossby waves

Two-layer equatorial primitive equations for the free troposphere in the presence of a thin atmospheric boundary layer and thermal dissipation are developed here. An asymptotic theory for the resonant nonlinear interaction of long equatorial baroclinic and barotropic Rossby waves is derived in the presence of such dissipation. In this model, a self-consistent asymptotic derivation establishes that boundary layer flows are generated by meridional pressure gradients in the lower troposphere and give rise to degenerate equatorial Ekman friction. That is to say, the asymptotic model has the property that the dissipation matrix has one eigenvalue which is nearly zero: therefore the dynamics rapidly dissipates flows with pressure at the base of the troposphere and creates barotropic/baroclinic spin up/spin down. The simplified asymptotic equations for the amplitudes of the dissipative equatorial barotropic and baroclinic waves are studied by linear theory and integrated numerically. The results indicate that although the dissipation slightly weakens the tropics to midlatitude connection, strong localized wave packets are nonetheless able to exchange energy between barotropic and baroclinic waves on intraseasonal timescales in the presence of baroclinic mean shear. Interesting dissipation balanced wave-mean flow states are discovered through numerical simulations. In general, the boundary layer dissipation is very efficient for flows in which the barotropic and baroclinic components are of the same sign at the base of the free troposphere whereas the boundary layer dissipation is less efficient for flows whose barotropic and baroclinic components are of opposite sign at the base of the free troposphere.     

安远盆地边缘断层高密度电法探测及活动特征

结合钻孔资料,采用高密度电法对安远盆地边缘的毛毛山南麓活动断层(F1)和金强河-毛毛山北麓活动断裂(F3)进行了探测.结果表明,F1断层错断了盆地区内晚更新世晚期至全新世早期的Ⅲ级阶地;受拉分构造的影响,该断层最新活动显示了向安远拉分盆地方向迁移的阶梯状递减的张性正断裂性质.F3断层在走向上呈锯齿状,错断了山前Ⅰ级冲洪积阶地及全新统至上更新统地层,断距大,基岩破碎带宽,新活动主要表现为断层附近原地连续沉降.这两条断层的发育特征及展布位置的确定,为该区工程稳定性评价及建设提供了依据.     

Influence of soil properties on active layer thermal propagation along the western Antarctic Peninsula

The rate of energy transfer through soils is an important factor governing the active layer (seasonal thaw layer) in polar regions. Energy is transferred through conductive and convective means, which are primarily influenced by the bulk density and water content of soils. With global temperatures changing, it becomes important to understand how soil properties influence heat transfer and active layer depths in climatically sensitive regions, such as the Antarctic Peninsula. In this study we analyzed conductive energy transfer through several soil types on Amsler Island and Cierva Point in the central region of the western Antarctic Peninsula. Active layer temperatures on Amsler Island were monitored every three hours using iButton thermistors installed at regular depth intervals down to 2 m. Soil textures were loamy to sandy with water contents between 5 and 27%. Freezing and thawing transmission rates for all soils ranged from 1.4 to 6.9 cm/day. Thermal transmission rates were fastest in sandy soils with low water contents, indicating that the large, interconnected pores of the sandy soils facilitated the quick movement of heat with water flow through the soil profile. Snow accumulation differences also played a significant role on winter thermal propagation by providing a thermal barrier between the ground surface and atmosphere. Although there was a wide range in thermal transmission among the soils, active layer depths had little variation (7.8–9.7 m). This consistency derives from the greater dependence of very thick active layers on long‐term climatic conditions rather than on soil properties. The presence of thick moss significantly slowed thermal transmission and decreased active layer thicknesses. These effects primarily are due to the high heat capacity of water and air retained within the moss, slowing thermal transmission rates, acting as a thermal buffer between atmospheric conditions and the underlying soils. Copyright © 2016 John Wiley & Sons, Ltd.     

Instability of flows near the onset of convection in a rotating layer with stress-free horizontal boundaries

We investigate instability of convective flows of simple structure (rolls, standing and travelling waves) in a rotating layer with stress-free horizontal boundaries near the onset of convection. We show that the flows are always unstable to perturbations, which are linear combinations of large-scale modes and short-scale modes, whose wave numbers are close to those of the perturbed flows. Depending on asymptotic relations of small parameters α (the difference between the wave number of perturbed flows and the critical wave number for the onset of convection) and ε (ε² being the overcriticality and the perturbed flow amplitude being O(ε)), either small-angle or Eckhaus instability is prevailing. In the case of small-angle instability for rolls the largest growth rate scales as ε^8/5, in agreement with results of Cox and Matthews (Cox, S.M. and Matthews, P.C., Instability of rotating convection. J. Fluid. Mech., 2000, 403, 153–172) obtained for rolls with k = k _c . For waves, the largest growth rate is of the order ε^4/3. In the case of Eckhaus instability the growth rate is of the order of α².