Attenuation and Shock Waves in Linear Hereditary Viscoelastic Media; Strick–Mainardi,Jeffreys–Lomnitz–Strick and Andrade Creep Compliances

Dispersion, attenuation and wavefronts in a class of linear viscoelastic media proposed by Strick and Mainardi (Geophys J R Astr Soc 69:415–429, 1982) and a related class of models due to Lomnitz, Jeffreys and Strick are studied by a new method due to the author. Unlike the previously studied explicit models of relaxation modulus or creep compliance, these two classes support propagation of discontinuities. Due to an extension made by Strick, either of these two classes of models comprise both viscoelastic solids and fluids. We also discuss the Andrade viscoelastic media. The Andrade media do not support discontinuity waves and exhibit the pedestal effect.     

Regional debris-flow and debris-flood frequency–magnitude relationships

Construction of frequency–magnitude (F–M) relationships of debris floods and debris flows is challenging because of few direct observations, discontinuous event occurrence, loss of field evidence, the difficulty of accessing the sediment archive and the challenge of finding suitable statistical methods to analyse the dataset. Consultants often face budget limitations that prohibit application of the full gamut of absolute dating methods, stratigraphic analysis and analytical tools necessary to fully resolve the F–M legacy. In some cases, F–M curves are needed for watersheds without local information, or where obtaining this information is prohibitively expensive. For such watersheds, the F–M relationship may be estimated where several F–M curves have already been assembled in a specific region. Individual F–M curves are normalized by fan area or fan volume, then stratified by process type and geomorphic activity level. This paper describes the development of regional F–M curves for debris flows in southwestern British Columbia and debris flows and debris floods in the Bow River valley near Canmore, Alberta. We apply the regional relationships to other cases in Canada and the United States and demonstrate that the method can be globalized. The regional approach is compared to cases where detailed F–M relationships have been established by other means. Strong negative deviations from the regional debris-flow or debris-flood magnitude trends could signal inherent watershed stability, while strong positive deviations could signal extraordinary landslide processes, or suggest that the fan may be largely of paraglacial origin. We highlight some of these outlying cases and develop a method whereby the regional curves can be meaningfully adjusted, or reliance can be placed on lower or upper confidence bounds of the F–M curves. We caution against the indiscriminate use of the regionally based F–M curves, especially in watersheds where multiple geomorphic processes are active. © 2020 John Wiley & Sons, Ltd.     

Modelling Gravitational Instabilities: Slab Break–off and Rayleigh–Taylor Diapirism

A non-standard new code to solve multiphase viscous thermo–mechanical problems applied to geophysics is presented. Two numerical methodologies employed in the code are described: A level set technique to track the position of the materials and an enrichment of the solution to allow the strain rate to be discontinuous across the interface. These techniques have low computational cost and can be used in standard desktop PCs. Examples of phase tracking with level set are presented in two and three dimensions to study slab detachment in subduction processes and Rayleigh–Taylor instabilities, respectively. The modelling of slab detachment processes includes realistic rheology with viscosity depending on temperature, pressure and strain rate; shear and adiabatic heating mechanisms; density including mineral phase changes and varying thermal conductivity. Detachment models show a first prolonged period of thermal diffusion until a fast necking of the subducting slab results in the break–off. The influence of several numerical and physical parameters on the detachment process is analyzed: The shear heating exerts a major influence accelerating the detachment process, reducing the onset time to one half and lubricating the sinking of the detached slab. The adiabatic heating term acts as a thermal stabilizer. If the mantle temperature follows an adiabatic gradient, neglecting this heating term must be included, otherwise all temperature contrasts are overestimated. As expected, the phase change at 410 km depth (olivine–spinel transition) facilitates the detachment process due to the increase in negative buoyancy. Finally, simple plume simulations are used to show how the presented numerical methodologies can be extended to three dimensions.     

Vertical–vertical controlled-source electromagnetic instrumentation and acquisition

Vertical–vertical controlled-source electromagnetic is an alternative to other techniques for providing three-dimensional resistivity images of the subsurface. It utilizes a large and powerful vertical dipole transmitter and arrays of E-field receivers with vertical and horizontal dipole sensors. The necessary instrumentation and acquisition procedures which differ strongly from other controlled-source electromagnetic methods are described in the paper.     

Pipe–soil interaction and pipeline performance under strike–slip fault movements

The performance of pipelines subjected to permanent strike–slip fault movement is investigated by combining detailed numerical simulations and closed-form solutions. First a closed-form solution for the force–displacement relationship of a buried pipeline subjected to tension is presented for pipelines of finite and infinite lengths. Subsequently the solution is used in the form of nonlinear springs at the two ends of the pipeline in a refined finite element model, allowing an efficient nonlinear analysis of the pipe–soil system at large strike–slip fault movements. The analysis accounts for large strains, inelastic material behavior of the pipeline and the surrounding soil, as well as contact and friction conditions on the soil–pipe interface. The numerical models consider infinite and finite length of the pipeline corresponding to various angles β between the pipeline axis and the normal to the fault plane. Using the proposed closed-form nonlinear force–displacement relationship for buried pipelines of finite and infinite length, axial strains are in excellent agreement with results obtained from detailed finite element models that employ beam elements and distributed springs along the pipeline length. Appropriate performance criteria of the steel pipeline are adopted and monitored throughout the analysis. It is shown that the end conditions of the pipeline have a significant influence on pipeline performance. For a strike–slip fault normal to the pipeline axis, local buckling occurs at relatively small fault displacements. As the angle between the fault normal and the pipeline axis increases, local buckling can be avoided due to longitudinal stretching, but the pipeline may fail due to excessive axial tensile strains or cross sectional flattening. Finally a simplified analytical model introduced elsewhere, is enhanced to account for end effects and illustrates the formation of local buckling for relative small values of crossing angle.     

Using GPS–GLONASS–GALILEO data and IRI modeling for ionospheric calibration of radio telescopes and radio interferometers

VHF and HF radio signals are widely used to observe the Sun and pulsars. Nowadays, large low-frequency radio astronomical arrays (LOFAR, 30–240 MHz; MIRA, 80–300 MHz) are being constructed to record radiation of pulsars at the maximum distance. registration of the solar radio emission intensity at fixed frequencies and in the spectral VHF band is very important along with other methods of monitoring of coronal mass ejections. Interpreting radio astronomical data is known to be necessary to take into account possible distortions of these signals in the Earth ionosphere. However, in contrast to modern navigation systems (Global Position System (GPS), GLObal NAvigation Satellite System (GLONASS), GALILEO), in which a very accurate reconstruction of ionosphere parameters is a built-in function, in present-day radio astronomy a retrieve of ionosphere transfer characteristics has not been appropriately worked out. This collides with increasing requirements to accuracy of the analysis of a radio emission amplitude profile and to the angular and polarizing resolution of radio telescopes of new generation (LOFAR, SKA, etc.). We have developed a method and software to calculate the ionosphere rotation measure (RM) and dispersion measure (DM). We used the ionosphere model IRI-2001, magnetic field model IGRF-10, and the ionosphere total electron content values obtained from GPS measurements. The obtained values of DM and RM were recalculated into characteristics of the phase delay, Faraday amplitude modulation, and polarization changes. We calculated ones for different levels of geomagnetic activity as well as different angular positions of radio sources. Our main idea is to use a signal of navigation satellites (GPS, GLONASS, GALILEO) as a testing signal from a “reference” source located at minimal angle distance from a source studied. Our project allows development of methods and systems of ADAPTIVE RADIO ASTRONOMY, adaptive to the non-uniform and non-stationary ionosphere, by analogy with known systems of adaptive optics intended to adapt optical telescopes to varying conditions of the optically non-uniform and non-stationary troposphere.     

Mixed fluid sources involved in diamond growth constrained by Sr–Nd–Pb–C–N isotopes and trace elements

Sub-micrometer inclusions in diamonds carry high-density fluids (HDF) from which the host diamonds have precipitated. The chemistry of these fluids is our best opportunity of characterizing the diamond-forming environment. The trace element patterns of diamond fluids vary within a limited range and are similar to those of carbonatitic/kimberlitic melts that originate from beneath the lithospheric mantle. A convecting mantle origin for the fluid is also implied by C isotopic compositions and by a preliminary Sr isotopic study (Akagi, T., Masuda, A., 1988. Isotopic and elemental evidence for a relationship between kimberlite and Zaire cubic diamonds. Nature 336, 665–667.). Nevertheless, the major element chemistry of HDFs is very different from that of kimberlites and carbonatites, varying widely and being characterized by extreme K enrichment (up to ～ 39 wt.% on a water and carbonate free basis) and high volatile contents. The broad spectrum of major element compositions in diamond-forming fluids has been related to fluid–rock interaction and to immiscibility processes.Elemental signatures can be easily modified by a variety of mantle processes whereas radiogenic isotopes give a clear fingerprint of the time-integrated evolution of the fluid source region. Here we present the results of the first multi radiogenic-isotope (Sr, Nd, Pb) and trace element study on fluid-rich diamonds, implemented using a newly developed off-line laser sampling technique. The data are combined with N and C isotope analysis of the diamond matrix to better understand the possible sources of fluid involved in the formation of these diamonds. Sr isotope ratios vary significantly within single diamonds. The highly varied but unsupported Sr isotope ratios cannot be explained by immiscibility processes or fluid-mineral elemental fractionations occurring at the time of diamond growth. Our results demonstrate the clear involvement of a mixed fluid, with one component originating from ancient incompatible element-enriched parts of the lithospheric mantle while the trigger for releasing this fluid source was probably carbonatitic/kimberlitic melts derived from greater depths. We suggest that phlogopite mica was an integral part of the enriched lithospheric fluid source and that breakdown of this mica releases K and radiogenic Sr into a fluid phase. The resulting fluids operate as a major metasomatic agent in the sub-continental lithospheric mantle as reflected by the isotopic composition and trace element patterns of G10 garnets.     

Dynamic response of structures subjected to pounding and structure–soil–structure interaction

During strong earthquakes, adjacent structures with non-sufficient clear distances collide with each other. In addition to such a pounding, cross interaction of adjacent structures through soil can exchange the vibration energy between buildings and make the problem even more complex. In this paper, effects of both of the mentioned phenomena on the inelastic response of selected steel structures are studied. Number of stories varied between 3 and 12 and different clear distances up to the seismic codes prescribed value are considered. The pounding element is modeled within Opensees. A coupled model of springs and dashpots is utilized for through-the-soil interaction of the adjacent structures, for two types of soft soils. The pounding force, relative displacements of stories, story shears, and plastic hinge rotations are compared for different conditions as the maximum responses averaged between seven consistent earthquakes. As a result, simultaneous effects of pounding and structure–soil–structure interaction are discussed.     

Formation and Suppression of Strike–Slip Fault Systems

Strike–slip faults are a defining feature of plate tectonics, yet many aspects of their development and evolution remain unresolved. For intact materials and/or regions, a standard sequence of shear development is predicted from physical models and field studies, commencing with the formation of Riedel shears and culminating with the development of a throughgoing fault. However, for materials and/or regions that contain crustal heterogeneities (normal and/or thrust faults, joints, etc.) that predate shear deformation, kinematic evolution of strike–slip faulting is poorly constrained. We present a new plane-stress finite-strain physical analog model developed to investigate primary deformation zone evolution in simple shear, pure strike–slip fault systems in which faults or joints are present before shear initiation. Experimental results suggest that preexisting mechanical discontinuities (faults and/or joints) have a marked effect on the geometry of such systems, causing deflection, lateral distribution, and suppression of shears. A lower limit is placed on shear offset necessary to produce a throughgoing fault in systems containing preexisting structures. Fault zone development observed in these experiments provides new insight for kinematic interpretation of structural data from strike–slip fault zones on Earth, Venus, and other terrestrial bodies.     

Role of wave–mean flow interaction in sun–climate connections: Historical overview and some new interpretations and results

Quasi-decadal variations in solar irradiance – termed the 11-year solar cycle (SC) – have been linked to variations in a variety of atmospheric circulation features, including the polar vortex, the Brewer–Dobson circulation, and the quasi-biennial oscillation. These features share an underlying commonality: they are all rooted in wave–mean flow interaction. The purpose of this paper is to provide a historical overview of the connection between the SC and wave–mean flow interaction and to propose a more complete theoretical framework for solar modulated wave–mean flow interaction that includes both zonal-mean and zonally asymmetric ozone as intermediaries for communicating variations in solar spectral irradiance to the climate system. We solve a quasi-geostrophic model using the WKB formalism to highlight the physics connecting the SC to planetary wave-drag. Numerical results show the importance of the zonally asymmetric ozone field in mediating the effects of solar variability to the wave-driven circulation in the middle atmosphere.     

Quasi-geostrophic shallow-water vortex–patch equilibria and their stability

We examine the equilibrium form, properties, stability and nonlinear evolution of steadily-rotating simply-connected vortex patches in the single-layer quasi-geostrophic model of geophysical fluid dynamics. This model, valid for rotating shallow-water flow in the limit of small Rossby and Froude numbers, has an intrinsic length scale L _D called the “Rossby deformation length” relating the strength of the stratification to that of the background rotation. Here, we generate steadily-rotating vortex equilibria for a wide range of γ?=?L/L _D , where L is the typical horizontal length scale of the vortex. We vary both γ (over the range 0.02?≤?γ?≤?10) and the vortex aspect ratio λ (over the range 0?<?λ?<?1). We find two modes of instability arising at sufficiently small aspect ratio λ?<?λ _c (γ): an asymmetric (dominantly wave 3) mode at small γ (or large L _D ) and a symmetric (dominantly wave 4) mode at large γ (or small L _D ). At marginal stability, the asymmetric mode dominates for γ???3, while the symmetric mode dominates for γ???3. The nonlinear evolution of weakly-perturbed unstable equilibria results in major structural changes, in most cases producing two dominant vortex patches and thin, quasi-passive filaments. Overall, the nonlinear evolution can be classified into three principal types: (1) vacillations for a limited range of aspect ratios λ when 5?≤?γ?≤?6, (2) filamentation and a single-dominant vortex for γ???1, and (3) vortex splitting – asymmetric for 1???γ???4 and symmetric for γ???4.     

Seismic soil–pile–structure kinematic and inertial interaction—A new beam approach

The main purpose of this study is to investigate the accuracy of an advanced beam model for the soil–pile–structure kinematic and inertial interaction and demonstrate its efficiency and advantages compared to other commonly used beam or solid models. Within this context, a Beam on Nonlinear Winkler Foundation model is adopted based on the Boundary Element Method (BEM), accounting for the effects induced by geometrical nonlinearity, rotary inertia and shear deformation, employing the concept of shear deformation coefficients. The soil nonlinearity is taken into consideration by means of a hybrid spring configuration consisting of a nonlinear (p–y) spring connected in series to an elastic spring–damper model. The nonlinear spring captures the near-field plastification of the soil while the spring–damper system (Kelvin–Voigt element) represents the far-field viscoelastic character of the soil. An extensive case study is carried out on a pile-column–deck system of a bridge, found in two cohesive layers of sharply different stiffness and subjected to various earthquake excitations, providing insight to several phenomena. The results of the proposed model are compared with those obtained from a Beam-FE solution as well as from a rigorous fully three-dimensional (3-D) continuum FE scheme.     

Integrated Surface and Subsurface Hydrological Modeling with Snowmelt and Pore Water Freeze–Thaw

For the simulation of winter hydrological processes a gap in the availability of flow models existed: one either had the choice between (1) physically-based and fully-integrated, but computationally very intensive, or (2) simplified and compartamentalized, but computationally less expensive, simulators. To bridge this gap, we here present the integration of a computationally efficient representation of winter hydrological processes (snowfall, snow accumulation, snowmelt, pore water freeze–thaw) in a fully-integrated surface water-groundwater flow model. This allows the efficient simulation of catchment-scale hydrological processes in locations significantly influenced by winter processes. Snow accumulation and snowmelt are based on the degree-day method and pore water freeze–thaw is calculated with a vertical heat conduction approach. This representation of winter hydrological processes is integrated into the fully-coupled surface water-groundwater flow model HydroGeoSphere. A benchmark for pore water freeze–thaw as well as two illustrative examples are provided.     

Nonlinear seismic soil–pile–structure interactions: Shaking table tests and FEM analyses

In this paper, a soil–pile–structure model is tested on a shaking table subject to both a sinusoidal wave and the acceleration time history of the scaled 1940 El Centro earthquake. A medium-size river sand is compacted into a 1.7-m-high laminar rectangular tank to form a loose fill with a relative density of 15%. A single-storey steel structure of 2.54 ton is placed on a concrete pile cap, which is connected to the four end-bearing piles. A very distinct pounding phenomenon between soil and pile is observed; and, the acceleration response of the pile cap can be three times larger than that of the structural response. The pounding is due to the development of a gap separation between soil and pile, and the extraordinary large inertia force suffered at the top of the pile also induces cracking in the pile. To explain this observed phenomenon, nonlinear finite element method (FEM) analyses with a nonlinear gap element have been carried out. The spikes in the acceleration response of the pile cap caused by pounding can be modeled adequately by the FEM analyses. The present results suggest that one of the probable causes of pile damages is due to seismic pounding between the laterally compressed soil and the pile near the pile cap level.     

Asymmetry of Sunspot Distribution in Solar Cycles 12–24 and the Gnevyshev–Ohl Rule

Geomagnetism and Aeronomy - The nonaxisymmetric component of the sunspot distribution (longitudinal asymmetry) is considered based on the Greenwich–USAF/NOAA data for 1874–2016. Vector...     

Amplitude–frequency characteristics of ion–cyclotron and whistler-mode waves from Van Allen Probes data

Using two-hour (from 2300 UT January 25, 2013 to 0100 UT January 26, 2013) measurement data from Van Allen Probes on fluxes of energetic particles, cold plasma density, and magnetic field magnitude, we have calculated the local growth rate of electromagnetic ion–cyclotron and whistler-mode waves for field-aligned propagation. The results of these calculations have been compared with wave spectra observed by the same Van Allen Probe spacecraft. The time intervals when the calculated wave increments are sufficiently large, and the frequency ranges corresponding to the enhancement peak agree with the frequency–time characteristics of observed electromagnetic waves. We have analyzed the influence of variations in the density and ionic composition of cold plasma, fluxes of energetic particles, and their pitch-angle distribution on the wave generation. The ducted propagation of waves plays an important role in their generation during the given event. The chorus VLF emissions observed in this event cannot be explained by kinetic cyclotron instability, and their generation requires much sharper changes (“steps”) for velocity distributions than those measured by energetic particle detectors on Van Allen Probes satellites.     

Saturated area dynamics and streamflow generation from coupled surface–subsurface simulations and field observations

A distributed physically-based model describing coupled surface–subsurface flows is applied to an instrumented catchment to investigate the links between runoff generation processes and the dynamics of saturated areas. The spatial characterization of the system is obtained through geophysical measurements and in situ observations. The model is able to reproduce the dynamics of the system through the calibration of only few parameters with a clear physical interpretation, providing a solid basis for our numerical investigations. Such investigations demonstrate the important control exerted by surface topography on the time evolution of saturated area patterns, mainly mediated by topographic curvature, that dictates both the dominant streamflow generation process at the local scale and the connection-disconnection dynamics of saturated areas. The relation between hillslope water storage and streamflow, Q = f(V), is shown to be highly hysteretical and dependent on the mean saturation of the catchment: higher degrees of saturation tend to yield one-to-one relationships between streamflow and water storage. On the contrary, streamflow-water storage relations are importantly affected by the specific configuration of saturated areas connected to the outlet when the system is far from complete saturation. This observation contradicts common assumptions of a one-to-one relationship Q = f(V) often used to justify widely observed power-law Q vs. dQ/dt recession curves. Furthermore, even when Q = f(V) becomes unique at high degrees of saturation, no power-law form emerged in our runs, speculatively because of the small size of the catchment formed by a single incision and the corresponding hillslope.     

On heliophysical and geophysical phenomena during October–November 2003

A study of the geomagnetic storm of November 20–21, 2003, is presented using Solar X-rays data, solar wind parameters and magnetic index, Dst. The results suggest that very large X class flares may not cause very intense geomagnetic storms, as flares of M importance would do. Furthermore, the results suggest that the solar wind structure that was responsible for this storm is of the shock-driver gas configuration in which the sheath is the most geoeffective element. Presently it is shown that an intense storm can be driven by two successive southward B _z structures without a resultant “double dip” at the boundary of these structures within the corresponding interval of the main phase. Furthermore, this study confirms earlier results that show that pressure enhancement does not cause the direct injection of new particles into the ring current region; rather it causes a local adiabatic energization of the particles already within the ring current region.