Mechanism of collision model for bedload transport

The current study introduces a new particle–particle collision model as well as a classical particle–wall collision model to simulate the process of bedload transport. Large eddy simulation and two-way coupling Lagrangian point–particle models also are applied. Flow conditions for all simulations are done based on previous experiments and four comparison cases are studied. The bedload function result obtained using this model agrees with the previous studies, which proves that the model is reasonable and effective in dealing with bedload transport in turbulent flow. Otherwise, failure to consider the particle–particle collision, Basset force, or Saffman lift force in the model lead to underestimating the bedload transport rate.     

Velocity field and tectonic strain in Southern Spain and surrounding areas derived from GPS episodic measurements

The goal of this paper is to study the velocity field and deformation parameters in Southern Spain and surrounding areas (Ibero-Maghrebian region) using GPS episodic measurements. Results are compared to those previously published as well as deformation parameters derived from seismic data. For this purpose, a geodetic GPS network of 12 stations was observed during eight field campaigns from 1998 to 2005 by the San Fernando Naval Observatory (ROA), Spain. Relative GPS velocities in the Gulf of Cadiz with respect to the stable part of Eurasia are ～4.1 mm/yr in a NW–SE to NNW–SSE direction. In the Betics, Alboran Sea and North of Morocco, velocities are ～4.4 mm/yr in a NW–SE direction, and they are ～2.3 mm/yr in a N–S direction in the eastern part of the Iberian Peninsula. These results are in agreement with the anticlockwise rotation of the African plate. GPS strain tensors are determined from the velocity model, to obtain a more realistic crustal deformation model. The Gulf of Cadiz is subjected to uniform horizontal compression in a NNW–SSE direction, with a rotation to N–S in the Alboran Sea and Northern Morocco. An extensional regime in a NW–SE direction, which rotates to W–E, is present in the Internal Betics area. In the Betic, Alboran Sea and North of Morocco regions we compare seismic deformation rates from shallow earthquakes with the determined GPS deformation rates. The comparison indicates a seismic coupling of 27%, while the remaining 73% might be generated in aseismic processes. Deformations measured in the Ibero-Maghrebian region with GPS could be interpreted in terms of either elastic loading or ductile deformation.     

Dynamic analysis of train–bridge system subjected to non‐uniform seismic excitations

Based on the theory of dynamic wheel–rail interactions, a dynamic model of coupled train–bridge system subjected to earthquakes is established, in which the non‐uniform characteristics of the seismic wave input from different foundations are considered. The bridge model is based on the modal comprehension analysis technique. Each vehicle is modelled with 31 degrees of freedom. The seismic loads are imposed on the bridge by using the influence matrix and exerted on the vehicles through the dynamic wheel–rail interaction relationships. The normal wheel–rail interaction is tackled by using the Hertzian contact theory, and the tangent wheel–rail interaction by the Kalker linear theory and the Shen–Hedrick–Elkins theory. A computer code is developed. A case study is performed to a continuous bridge on the planned Beijing–Shanghai high‐speed railway in China. Through input of typical seismic waves with different propagation velocities to the train–bridge system, the histories of the train running through the bridge are simulated and the dynamic responses of the bridge and the vehicles are calculated. The influences of train speed and seismic wave propagation velocity on the dynamic responses of the bridge–vehicle system are studied. The critical train speeds are proposed for running safety on high‐speed railway bridges under earthquakes of various intensities. Copyright © 2006 John Wiley & Sons, Ltd.     

A numerical model for the dynamic analysis of inclined pile groups

The paper presents a numerical model for the dynamic analysis of pile groups with inclined piles in horizontally layered soil deposits. Piles are modelled with Euler–Bernoulli beams, while the soil is supposed to be constituted by independent infinite viscoelastic horizontal layers. The pile–soil–pile interaction as well as the hysteretic and geometric damping is taken into account by means of two‐dimensional elastodynamic Green's functions. Piles cap is considered by introducing a rigid constraint; the condensation of the problem permits a consistent derivation of both the dynamic impedance matrix of the soil–foundation system and the foundation input motion. These quantities are those used to perform inertial soil–structure interaction analyses in the framework of the substructure approach. Furthermore, the model allows evaluating the kinematic stress resultants in piles resulting from waves propagating in the soil deposit, taking into account the pile–soil–pile interactions. The model validation is carried out by performing accuracy analyses and comparing results in terms of dynamic impedance functions, kinematic response parameters and pile stress resultants, with those furnished by 3D refined finite element models. To this purpose, classical elastodynamic solutions are adopted to define the soil–pile interaction problem. The model results in low computational demands without significant loss of precision, compared with more rigorous approaches or refined finite element models. Copyright © 2015 John Wiley & Sons, Ltd.     

Dynamic response of an elastic-viscoplastic system in modal co-ordinates

A dynamic analysis of elastic–viscoplastic systems, incorporating the modal co-ordinate transformation technique, is presented. The formulation results in uncoupled incremental equations of motion with respect to the modal co-ordinates. The elastic–viscoplastic model adopted allows the analysis not to involve yielding regions and loading/unloading processes. An implicit Runge–Kutta scheme together with the Newton–Raphson method are used to solve the non-linear constitutive equations. Stability and accuracy of the numerical solution are improved by utilizing a local time step sub-incrementing procedure. Applications of the analyses to multi-storey shear buildings show that good results can be obtained for the maximum displacement response by including only a few lower modes in the computation, but the prediction of the ductility factor response tends to underestimate the peak values when too few modes are used. In addition, stable and valid results can be obtained even with a sizable time step increment.     

Channel and floodplain response to recent abrupt climate change: The tyne basin,Northern England

This paper examines the timing, nature and magnitude of river response in upland, piedmont and lowland reaches of the Tyne basin, northern England, to high-frequency (20–30 year) changes in climate and flood regime since 1700 AD. Over this period fluvial activity has been characterized by alternating phases of river-bed incision and stability coinciding with non-random, decadal-scale fluctuations in flood frequency and hydroclimate that appear to be linked to changes in large-scale upper atmospheric circulation patterns. Episodes of widespread channel bed incision (1760–1799, 1875–1894, 1955–1969) result from a higher frequency of large floods (> 20 year return period) and cool, wet climate under meridional circulation regimes. Phases of more moderate floods (5–20 year return period), corresponding to zonal circulation types (1820–1874, 1920–1954), are characterized by enhanced lateral reworking and sediment transfer in upper reaches of the catchment, and channel narrowing and infilling downstream. Rates of fluvial activity are reduced in intermediate periods (1800–1819, 1895–1919) with no dominant circulation regime associated with lower flood frequency and magnitude. The results of this study provide a valuable guide for forecasting probable drainage basin and channel response to future climate change.     

Discrete modelling of vertical track–soil coupling for vehicle–track dynamics

This paper presents a coupled lumped mass model (CLM model) for the vertical dynamic coupling of railway track through the soil. The well-known Winkler model and its extensions are analysed and fitted on the result obtained numerically with a finite–infinite element model in order to validate the approach in a preliminary step. A mass–spring–damper system with frequency independent parameters is then proposed for the interaction between the foundations, representing the contact area of the track with the soil. The frequency range of track–soil coupling is typically under 100 Hz. Analytical expressions are derived for calibrating the system model with homogeneous and layered half-spaces. Numerical examples are derived, with emphasis on soil stiffness and layering. The dynamic analysis of a track on various foundation models is compared with a complete track–soil model, showing that the proposed CLM model captures the dynamic interaction of the track with the soil and is reliable to predict the vertical track deflection and the reaction forces acting on the soil surface.     

Ionospheric slab thickness and its seasonal variations observed by GPS

The ionospheric slab thickness, the ratio of the total electron content (TEC) to the F2-layer peak electron density (NmF2), is closely related to the shape of the ionospheric electron density profile Ne (h) and the TEC. Therefore, the ionospheric slab thickness is a significant parameter representative of the ionosphere. In this paper, the continuous GPS observations in South Korea are firstly used to study the equivalent slab thickness (EST) and its seasonal variability. The averaged diurnal medians of December–January–February (DJF), March–April–May (MAM), June–July–August (JJA) and September–October–November (SON) in 2003 have been considered to represent the winter, spring, summer and autumn seasons, respectively. The results show that the systematic diurnal changes of TEC, NmF2 and EST significantly appeared in each season and the higher values of TEC and NmF2 are observed during the equinoxes (semiannual anomaly) as well as in the mid-daytime of each season. The EST is significantly smaller in winter than in summer, but with a consistent variation pattern. During 14–16 LT in daytime, the larger EST values are observed in spring and autumn, while the smaller ones are in summer and winter. The peaks of EST diurnal variation are around 10–18 LT which are probably caused by the action of the thermospheric wind and the plasmapheric flow into the F2-region.     

Experimental study of near-shore pile-to-pile interaction

This paper presents the results of lateral impact load field tests carried out on a system of three steel pipe piles vibro-driven into soft clay in a near-shore marine environment, with the aim of evaluating the pile–soil–pile dynamic interaction. Piles are arranged in an “L” shaped horizontal layout and are instrumented with accelerometers at their free heads. The obtained results show the complex dynamic behaviour at very small strain of the vibrating soil–water–piles system. The role of different type of waves in the pile to pile interaction is investigated by analyzing the results in the time and frequency domains and by means of a time–frequency analysis. The effects of the pile spacing and input direction on these interaction mechanisms are also presented. Finally, important dynamic parameters of the soil, such as the velocities of the shear waves and surface waves (Scholte waves) of the upper soil are directly estimated from the time delays between signals recorded at the pile heads.     

A model for the 3D kinematic interaction analysis of pile groups in layered soils

The paper presents a numerical model for the analysis of the soil–structure kinematic interaction of single piles and pile groups embedded in layered soil deposits during seismic actions. A finite element model is considered for the pile group and the soil is assumed to be a Winkler‐type medium. The pile–soil–pile interaction and the radiation problem are accounted for by means of elastodynamic Green's functions. Condensation of the problem permits a consistent and straightforward derivation of both the impedance functions and the foundation input motion, which are necessary to perform the inertial soil–structure interaction analyses. The model proposed allows calculating the internal forces induced by soil–pile and pile‐to‐pile interactions. Comparisons with data available in literature are made to study the convergence and validate the model. An application to a realistic pile foundation is given to demonstrate the potential of the model to catch the dynamic behaviour of the soil–foundation system and the stress resultants in each pile. Copyright © 2009 John Wiley & Sons, Ltd.     

Zircon fission‐track and U–Pb double dating using femtosecond laser ablation–inductively coupled plasma–mass spectrometry: A technical note

We present a new LA–ICP–MS system for zircon fission‐track (FT) and U–Pb double dating, whereby a femtosecond laser combined with galvanometric optics simultaneously ablates multiple spots to measure average surface U contents. The U contents of zircon measured by LA–ICP–MS and standardized with the NIST SRM610 glass are comparable to those measured by the induced FT method, and have smaller analytical errors. LA–ICP–MS FT dating of seven zircon samples including three IUGS age standards is as accurate as the external detector method, but can give a higher‐precision age depending on the counting statistics of the U content measurement. Double dating of the IUGS age standards gives FT and U–Pb ages that are in agreement. A chip of the Nancy 91500 zircon has a homogeneous U content of 84 ppm, suggesting the possibility of using this zircon as a matrix‐matched U‐standard for FT dating. When using the Nancy 91500 zircon as a U‐standard, a zeta calibration value of 42–43 year cm² for LA–ICP–MS FT dating is obtained. While this value is strictly valid only for the particular session, it can serve as a reference for other studies.     

Drought hotspot analysis and risk assessment using probabilistic drought monitoring and severity–duration–frequency analysis

Drought hotspot identification requires continuous drought monitoring and spatial risk assessment. The present study analysed drought events in the agriculture‐dominated mid‐Mahanadi River Basin in Odisha, India, using crop water stress as a drought indicator. This drought index incorporated different factors that affect crop water deficit such as the cropping pattern, soil characteristics, and surface soil moisture. The drought monitoring framework utilized a relevance vector machine model‐based classification that provided the uncertainty associated with drought categorization. Using the proposed framework, drought hotspots are identified in the study region and compared with indices based on precipitation and soil moisture. Further, a bivariate copula is employed to model the agricultural drought characteristics and develop the drought severity–duration–frequency (S–D–F) relationships. The drought hotspot maps and S–D–F curves are developed for different locations in the region. These provided useful information on the site‐specific drought patterns and the characteristics of the devastating droughts of 2002 and 2012, characterized by an average drought duration of 7 months at several locations. The site‐specific risk of short‐ and long‐term agricultural droughts are then investigated using the conditional copula. The results suggest that the conditional return periods and the S–D–F curves are valuable tools to assess the spatial variability of drought risk in the region.     

Inelastic displacement ratios for soil–structure systems allowed to uplift

The simultaneous effects of soil–structure interaction, foundation uplift and inelastic behavior of the superstructure on total displacement response of soil–structure systems are investigated. The superstructure is modeled as an equivalent single‐degree‐of‐freedom system with bilinear behavior mounted on a rigid foundation resting on distributed tensionless Winkler springs and dampers. It is well known that the behavior of soil–structure systems can be well described using a limited number of nondimensional parameters. Here, by introducing two new parameters, the concept is extended to inelastic soil–structure systems in which the foundation is allowed to uplift. An extensive parametric study is conducted for a wide range of the key parameters through nonlinear time history analyses. It is shown that while uplifting soil–structure systems experience excessive displacements, in comparison with systems that are not allowed to uplift, ductility demand in the superstructure generally decreases owing to foundation uplift. A new inelastic displacement ratio (IDR) is proposed in conjunction with a nonlinear static analysis of uplifting soil–structure systems. Simplified expressions are also provided to estimate the proposed IDR. Copyright © 2014 John Wiley & Sons, Ltd.     

Transformation of spatial and perturbation derivatives of travel time at a curved interface between two arbitrary media

We consider the partial derivatives of travel time with respect to both spatial coordinates and perturbation parameters. These derivatives are very important in studying wave propagation and have already found various applications in smooth media without interfaces. In order to extend the applications to media composed of layers and blocks, we derive the explicit equations for transforming these travel–time derivatives of arbitrary orders at a general smooth curved interface between two arbitrary media. The equations are applicable to both real–valued and complex–valued travel time. The equations are expressed in terms of a general Hamiltonian function and are applicable to the transformation of travel–time derivatives in both isotropic and anisotropic media. The interface is specified by an implicit equation. No local coordinates are needed for the transformation.     

Nonlinear dynamic foundation and frame structure response observed in geotechnical centrifuge experiments

Soil–foundation–structure interaction (SFSI) and structure–soil–structure interaction (SSSI) influence the seismic response of a structure. Yet, consideration of nonlinear SFSI and SSSI in design practice is lacking. In this paper data from two centrifuge tests are examined. During each test, inelastic models of (1) a low-rise frame with shallowly embedded footings and (2) a mid-rise frame with a large basement are subjected to earthquake motions. In the first test, the structures are separated. In the second test, the structures are placed next to each other. Results show that the presence of the deep basement affects the moment–rotation behavior of the adjacent shallow footings, stiffening the response in the direction of loading towards the basement. This can be attributed to the additional restraint provided by the basement. Although the presence of the basement stiffens the response, it also limits the permanent displacements of the footing, which in turn limits physical damage to the superstructure. These results suggest that in addition to considering nonlinear SFSI effects, SSSI should be considered in the design of closely clustered structures.     

Alluvial architecture in headwater streams with special emphasis on step–pool topography

Alluvial mountain streams exhibit a range of channel forms: pool–riffle, plane bed, step–pool and cascades. Previous work suggested that these forms exist within discrete, and progressively steeper slope classes. Measurements conducted at over 100 sites in west‐central and central Idaho confirm that slope steepens progressively as one moves from pool–riffle, to plane bed, to step–pool, and finally to cascades. Median slope for pool–riffle topography is 0·0060, for plane beds 0·013, for step–pools 0·044, and for cascades 0·068. There is substantial overlap in the slopes associated with these channel forms. Pool–riffle topography was found at slopes between 0·0010 and 0·015, plane beds between 0·0010 and 0·035, step–pools between 0·015 and 0·134, and cascades between 0·050 and 0·12. Step–pools are particularly striking features in headwater streams. They are characterized by alternating steep and gentle channel segments. The steep segments (step risers) are transverse accumulations of boulder and cobbles, while the gentle segments (pools) contain finer material. Step wavelength is best correlated to step height which is in turn best correlated to the median particle size found on step risers. This result differs from past studies that have reported channel slope to be the dominant control on step wavelength. The presumed geometry and Froude number associated with the features under formative conditions are consistent with the existence field for antidunes and by extension with the hypothesis that step–pools are formed by antidunes. Copyright © 2000 John Wiley & Sons, Ltd.     

Tectonic setting of the recent damaging seismic series in the Southeastern Betic Cordillera,Spain

In this work we analyze the tectonic setting of the recent damaging seismic series occurred in the Internal Zones of the eastern Betic Cordillera (SE Spain) and surrounding areas, the tectonic region where took place the 11th May 2011 Mw 5.2 Lorca earthquake. We revisit and make a synthesis of the seven largest and damaging seismic series occurred from 1984 to 2011. We analyze their seismotectonic setting, and their geological sources under the light of recent advances in the knowledge on active faults, neotectonics, seismotectonics and stress regime, with special attention focused on the Lorca Earthquake. These seismic series are characterized by two types of focal mechanisms, produced mainly by two sets of active faults, NNW–SSE to NNE–SSW small (no larger than 20–30 km) extensional faults with some strike slip component, and E–W to NE–SW large strike slip faults (more than 50 km long) with some compressional component (oblique slip faults). The normal fault earthquakes related to the smaller faults are dominant in the interior of large crustal tectonic blocks that are bounded by the large E–W to NE–SW strike-slip faults. The strike slip earthquakes are associated to the reactivation of segments or intersegment regions of the large E–W to NE–SW faults bounding those crustal tectonic blocks. Most of the seismic series studied in this work can be interpreted as part of the background seismicity that occurs within the crustal blocks that are strained under a transpressional regime driven by the major strike slip shear corridors bounding the blocks. The seismotectonic analysis and the phenomenology of the studied series indicate that it is usual the occurrence of damaging compound earthquakes of M  \(\sim \)  5.0 associated with triggering processes driven by coseismic stress transfer. These processes mainly occur in the seismic series generated by NNW–SSE to NNE–SSW faults. These mechanical interaction processes may induce a higher frequency of occurrence of this kind of earthquakes than considered in traditional probabilistic seismic hazard assessments and it should be taken into account in future seismic hazard assessments.     

Comparison of the anisotropic-ray-theory rays and anisotropic common S-wave rays with the SH and SV reference rays in a velocity model with a split intersection singularity

We describe the behaviour of the anisotropic–ray–theory S–wave rays in a velocity model with a split intersection singularity. The anisotropic–ray–theory S–wave rays crossing the split intersection singularity are smoothly but very sharply bent. While the initial–value rays can be safely traced by solving Hamilton’s equations of rays, it is often impossible to determine the coefficients of the equations of geodesic deviation (paraxial ray equations, dynamic ray tracing equations) and to solve them numerically. As a result, we often know neither the matrix of geometrical spreading, nor the phase shift due to caustics. We demonstrate the abrupt changes of the geometrical spreading and wavefront curvature of the fast anisotropic–ray–theory S wave. We also demonstrate the formation of caustics and wavefront triplication of the slow anisotropic–ray–theory S wave.Since the actual S waves propagate approximately along the SH and SV reference rays in this velocity model, we compare the anisotropic–ray–theory S–wave rays with the SH and SV reference rays. Since the coupling ray theory is usually calculated along the anisotropic common S–wave rays, we also compare the anisotropic common S–wave rays with the SH and SV reference rays.     

Damping reduction factors and code‐based design equation for structures using semi‐active viscous dampers

This study uses a semi‐active viscous damper with three different control laws to reshape the structural hysteresis loop and mitigate structural response, referred to as 1–4, 1–3 and 2–4 devices, respectively. The 1–4 control law provides damping in all four quadrants of the force‐displacement graph (it behaves like a standard viscous damper), the 1–3 control law provides resisting forces only in the first and third quadrants, and the 2–4 control law provides damping in the second and fourth quadrants. This paper first outlines the linear single degree of freedom structural performance when the three types of semi‐active viscous dampers are applied. The results show that simultaneous reduction in both displacement and base‐shear demand is only available with the semi‐active 2–4 device. To enable guidelines for adding a 2–4 device into the design procedure, damping reduction factors (RF_ξs) are developed, as they play an important role and provide a means of linking devices to design procedures. Three methods are presented to obtain RF_ξ and equivalent viscous damping of a structure with a 2–4 semi‐active viscous damper. In the first method, the relationship between RF_ξ and the damping of a semi‐active structure can be obtained by calculating the area under the force‐deformation diagram. The second and third method modified the Eurocode8 formula of RF_ξ and smoothed results from analysis, respectively. Finally, a simple method is proposed to incorporate the design or retrofit of structures with simple, robust and reliable 2–4 semi‐active viscous dampers using standard design approaches. Copyright © 2016 John Wiley & Sons, Ltd.     

Observational evidence of deep convection over Indonesian sector in relation with major stratospheric warming events of 2003–04 and 2005–06

The major sudden stratospheric warming (SSW) events of 2003–04 and 2005–06 are considered to investigate changes in equatorial convection due to circulation changes associated with the SSW events. It is observed that the SSW events are accompanied by a considerable decrease in Outgoing Longwave Radiation (OLR), a proxy for tropical convection, over equatorial latitudes (15°N–15°S) in the Indonesian sector (90°E–150°E). However, unlike noted by earlier observations, the zonal mean OLR does not show any notable relationship with the SSW events. It can be explained from the latitude–longitude map of potential vorticity (PV) at 100 hPa, which shows a tongue of high PV emanating from high latitudes towards equator and converges in the longitude band of 90°E–150°E on the day of peak warming at 1 hPa in the case of 2003–04 and 10 hPa in the case of 2005–06. The latitude-height map of Eliassen–Palm (EP) vector and its divergence show convergence of EP flux in the upper troposphere at latitudes even lower than 20°N on these days. Further, vertical winds computed from the convergence of momentum flux are upward indicating convective activity at low-latitudes and downward at mid-latitudes.