Estimating the properties of weathered bedrock and pilerock interaction from the geological strength index

The effect of variable rock mass properties on pile-rock interaction poses a great challenge to the design of stabilizing piles and numerical analysis of pile-rock interaction. The paper presents a novel method to estimate the properties of weathered bedrock, which can be applied to routine design of landslide-stabilizing piles for collivial landslides. The Ercengyan landslide located in the Three Gorges Reservoir, China, is the area of interest for this study. A geological investigation and triaxial tests were conducted to estimate the basic parameters, including Geological Strength Index(GSI), uniaxial compressive strength σ_(ci) and Hoek-Brown constant m_i of intact bedrock in the study area. Hoek-Brown criterion was used to estimate mechanical properties of the weathered rock, including elastic modulus E_m, cohesion c, friction angle Φ, and normal ultimate lateral resistance p_(max). A parametric study was performed to evaluate the effect of parameterizations of GSI, σ_(ci) and m_i on the bedrock properties and p-y curves. The estimated rock mass properties were used with PLAXIS 2D software to simulate pile-rock interaction. Effect of GSI on stress at the pile-rock interface and in the rock, pile bending moment, pile shear force, and p-y curve were analysed.     

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

Analytical study of the swing-by maneuver in an elliptical system

The objective of the present paper is to derive a set of analytical equations that describe a swing-by maneuver realized in a system of primaries that are in elliptical orbits. The goal is to calculate the variations of energy, velocity and angular momentum as a function of the usual basic parameters that describe the swing-by maneuver, as done before for the case of circular orbits. In elliptical orbits the velocity of the secondary body is no longer constant, as in the circular case, but it varies with the position of the secondary body in its orbit. As a consequence, the variations of energy, velocity and angular momentum become functions of the magnitude and the angle between the velocity vector of the secondary body and the line connecting the primaries. The “patched-conics” approach is used to obtain these equations. The configurations that result in maximum gains and losses of energy for the spacecraft are shown next, and a comparison between the results obtained using the analytical equations and numerical simulations are made to validate the method developed here.     

Fine-resolution analysis of exoplanetary distributions by wavelets: hints of an overshooting iceline accumulation

We investigate 1D exoplanetary distributions using a novel analysis algorithm based on the continuous wavelet transform. The analysis pipeline includes an estimation of the wavelet transform of the probability density function (p.d.f.) without pre-binning, use of optimized wavelets, a rigorous significance testing of the patterns revealed in the p.d.f., and an optimized minimum-noise reconstruction of the p.d.f. via matching pursuit iterations.In the distribution of orbital periods, \(P\), our analysis revealed a narrow subfamily of exoplanets within the broad family of “warm Jupiters”, or massive giants with \(P\gtrsim 300~\mbox{d}\), which are often deemed to be related with the iceline accumulation in a protoplanetary disk. We detected a p.d.f. pattern that represents an upturn followed by an overshooting peak spanning \(P\sim 300\mbox{--}600~\mbox{d}\), right beyond the “period valley”. It is separated from the other planets by p.d.f. concavities from both sides. It has at least 2-sigma significance.In the distribution of planet radii, \(R\), and using the California Kepler Survey sample properly cleaned, we confirm the hints of a bimodality with two peaks about \(R=1.3R_{\oplus }\) and \(R=2.4R_{ \oplus }\), and the “evaporation valley” between them. However, we obtain just a modest significance for this pattern, 2-sigma only at the best. Besides, our follow-up application of the Hartigan and Hartigan dip test for unimodality returns 3 per cent false alarm probability (merely 2.2-sigma significance), contrary to 0.14 per cent (or 3.2-sigma), as claimed by Fulton et al. (2017).     

Understanding the Role of Mass-Unloading in a Filament Eruption

We describe a partial filament eruption on 11 December 2011 that demonstrates that the inclusion of mass is an important next step for understanding solar eruptions. Observations from the Solar Terrestrial Relations Observatory-Behind (STEREO-B) and the Solar Dynamics Observatory (SDO) spacecraft were used to remove line-of-sight projection effects in filament motion and correlate the effect of plasma dynamics with the evolution of the filament height. Flux cancellation and nearby flux emergence are shown to have played a role in increasing the height of the filament prior to eruption. The two viewpoints allow the quantitative estimation of a large mass-unloading, the subsequent radial expansion, and the eruption of the filament to be investigated. A 1.8 to 4.1 lower-limit ratio between gravitational and magnetic-tension forces was found. We therefore conclude that following the loss-of-equilibrium of the flux-rope, the radial expansion of the flux-rope was restrained by the filamentary material until 70% of the mass had evacuated the structure through mass-unloading.     

Influence of periodic orbits on the formation of giant planetary systems

The late-stage formation of giant planetary systems is rich in interesting dynamical mechanisms. Previous simulations of three giant planets initially on quasi-circular and quasi-coplanar orbits in the gas disc have shown that highly mutually inclined configurations can be formed, despite the strong eccentricity and inclination damping exerted by the disc. Much attention has been directed to inclination-type resonance, asking for large eccentricities to be acquired during the migration of the planets. Here we show that inclination excitation is also present at small to moderate eccentricities in two-planet systems that have previously experienced an ejection or a merging and are close to resonant commensurabilities at the end of the gas phase. We perform a dynamical analysis of these planetary systems, guided by the computation of planar families of periodic orbits and the bifurcation of families of spatial periodic orbits. We show that inclination excitation at small to moderate eccentricities can be produced by (temporary) capture in inclination-type resonance and the possible proximity of the non-coplanar systems to spatial periodic orbits contributes to maintaining their mutual inclination over long periods of time.     

Assessment of alternative adsorption models and global sensitivity analysis to characterize hexavalent chromium loss from soil to surface runoff

We investigate our ability to assess transfer of hexavalent chromium, Cr(VI), from the soil to surface runoff by considering the effect of coupling diverse adsorption models with a two‐layer solute transfer model. Our analyses are grounded on a set of two experiments associated with soils characterized by diverse particle size distributions. Our study is motivated by the observation that Cr(VI) is receiving much attention for the assessment of environmental risks due to its high solubility, mobility, and toxicological significance. Adsorption of Cr(VI) is considered to be at equilibrium in the mixing layer under our experimental conditions. Four adsorption models, that is, the Langmuir, Freundlich, Temkin, and linear models, constitute our set of alternative (competing) mathematical formulations. Experimental results reveal that the soil samples characterized by the finest grain sizes are associated with the highest release of Cr(VI) to runoff. We compare the relative abilities of the four models to interpret experimental results through maximum likelihood model calibration and four model identification criteria (i.e., the Akaike information criteria [AIC and AIC_C] and the Bayesian and Kashyap information criteria). Our study results enable us to rank the tested models on the basis of a set of posterior weights assigned to each of them. A classical variance‐based global sensitivity analysis is then performed to assess the relative importance of the uncertain parameters associated with each of the models considered, within subregions of the parameter space. In this context, the modelling strategy resulting from coupling the Langmuir isotherm with a two‐layer solute transfer model is then evaluated as the most skilful for the overall interpretation of both sets of experiments. Our results document that (a) the depth of the mixing layer is the most influential factor for all models tested, with the exception of the Freundlich isotherm, and (b) the total sensitivity of the adsorption parameters varies in time, with a trend to increase as time progresses for all of the models. These results suggest that adsorption has a significant effect on the uncertainty associated with the release of Cr(VI) from the soil to the surface runoff component.     

Study of the flood control scheduling scheme for the Three Gorges Reservoir in a catastrophic flood

The Three Gorges Project is the world's largest water conservancy project. According to the design standards for the 1,000‐year flood, flood diversion areas in the Jingjiang reach of the Yangtze River must be utilized to ensure the safety of the Jingjiang area and the city of Wuhan. However, once these areas are used, the economic and life loss in these areas may be very great. Therefore, it is vital to reduce this loss by developing a scheme that reduces the use of the flood diversion areas through flood regulation by the Three Gorges Reservoir (TGR), under the premise of ensuring the safety of the Three Gorges Dam. For a 1,000‐year flood on the basis of a highly destructive flood in 1954, this paper evaluates scheduling schemes in which flood diversion areas are or are not used. The schemes are simulated based on 2.5‐m resolution reservoir topography and an optimized model of dynamic capacity flood regulation. The simulation results show the following. (a) In accord with the normal flood‐control regulation discharge, the maximum water level above the dam should be not more than 175 m, which ensures the safety of the dam and reservoir area. However, it is necessary to utilize the flood diversion areas within the Jingjiang area, and flood discharge can reach 2.81 billion m³. (b) In the case of relying on the TGR to impound floodwaters independently rather than using the flood diversion areas, the maximum water level above the dam reaches 177.35 m, which is less than the flood check level of 180.4 m to ensure the safety of the Three Gorges Dam. The average increase of the TGR water level in the Chongqing area is not more than 0.11 m, which indicates no significant effect on the upstream reservoir area. Comparing the various scheduling schemes, when the flood diversion areas are not used, it is believed that the TGR can execute safe flood control for a 1,000‐year flood, thereby greatly reducing flood damage.     

A near-inertial current event in the homogeneous deep layer of the northern Sea of Japan during winter

Under strong surface wind forcing during winter, direct current observations in the northern Sea of Japan show the existence of strong near-inertial currents in the deep water that is characterized by the extremely homogeneous vertical structures of temperature and salinity. However, the mechanism generating internal waves in the deep water of the northern Sea of Japan has not been well understood. In this study, to clarify the dynamical link between the surface wind forcing and near-inertial currents in the deep water of the northern Sea of Japan, we drive a general circulation model taking into account realistic wind stress, ocean bottom and land topography. In the northern Sea of Japan, the numerical results show that vertically coherent horizontal currents with a speed of ~ 0.05 m s^?1 are excited throughout the homogeneous deep water. A two-layer model successfully reproduces the pattern of the horizontal current velocities shown by the general circulation model, indicating that internal waves emanate westward from the northwestern coast of Japan through coastal adjustment to the strong wind forcing event and, while propagating into the ocean interior, they excite evanescent near-inertial response throughout the lower layer below the interface.