Spatial patch structure and adaptive strategy for desert shrub of Reaumuria soongorica in arid ecosystem of the Heihe River Basin

In many arid ecosystems, vegetation frequently occurs in high-cover patches interspersed in a matrix of low plant cover. However, theoretical explanations for shrub patch pattern dynamics along climate gradients remain unclear on a large scale. This context aimed to assess the variance of the Reaumuria soongorica patch structure along the precipitation gradient and the factors that affect patch structure formation in the middle and lower Heihe River Basin (HRB). Field investigations on vegetation patterns and heterogeneity in soil properties were conducted during 2014 and 2015. The results showed that patch height, size and plant-to-patch distance were smaller in high precipitation habitats than in low precipitation sites. Climate, soil and vegetation explained 82.5% of the variance in patch structure. Spatially, R. soongorica shifted from a clumped to a random pattern on the landscape towards the MAP gradient, and heterogeneity in the surface soil properties (the ratio of biological soil crust (BSC) to bare gravels (BG)) determined the R. soongorica population distribution pattern in the middle and lower HRB. A conceptual model, which integrated water availability and plant facilitation and competition effects, was revealed that R. soongorica changed from a flexible water use strategy in high precipitation regions to a consistent water use strategy in low precipitation areas. Our study provides a comprehensive quantification of the variance in shrub patch structure along a precipitation gradient and may improve our understanding of vegetation pattern dynamics in the Gobi Desert under future climate change.

     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

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.     

Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

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).     

Flexural deformation controls on Late Quaternary sediment dispersal in the Garo-Rajmahal Gap,NW Bengal Basin

Subsurface deformation is a driver for river path selection when deformation rates become comparable to the autogenic mobility rate of rivers. Here we combine geomorphology, soil and sediment facies analyses, and geophysical data of the Late Quaternary sediments of the central Garo-Rajmahal Gap in Northwest Bengal to link subsurface deformation with surface processes. We show variable sedimentation characteristics, from slow rates (<0.8 mm/year) in the Tista megafan at the foot of the Himalaya to nondeposition at the exposed surface of the Barind Tract to the south, enabling the development of mature soils. Combined subsidence in the Tista fan and uplift of the Barind Tract are consistent with a N-S flexural response of the Indian plate to loading of the Himalaya Mountains given a low value of elastic thickness (15–25 km). Provenance analysis based on bulk strontium concentration suggests a dispersal of sediment consistent with this flexural deformation—in particular the abandonment of the Barind Tract by a Pleistocene Brahmaputra River and the current extents of the Tista megafan lobes. Overall, these results highlight the control by deeply rooted deformation patterns on the routing of sediment by large rivers in foreland settings.     

Jurassic igneous rocks of the central Sanandaj–Sirjan zone (Iran) mark a propagating continental rift,not a magmatic arc

Jurassic igneous bodies of the Sanandaj–Sirjan zone (SaSZ) in SW Iran are generally considered as a magmatic arc but critical evaluation of modern geochronology, geochemistry and radiogenic isotopes challenges this conclusion. There is no evidence for sustained igneous activity along the ~1,200 km long SaSZ, as expected for a convergent plate margin; instead activity was brief at most sites and propagated NW at ~20 mm/a. Jurassic igneous rocks define a bimodal suite of gabbro‐diorite and granite. Chemical and isotopic compositions of mafic rocks indicate subcontinental lithospheric mantle sources that mostly lacked subduction‐related modifications. The arc‐like features of S‐type granites reflect massive involvement of Cadomian crust and younger sediments to generate felsic melts in response to mafic intrusions. We conclude that Jurassic SaSZ igneous activity occurred in a continental rift, not an arc. SaSZ igneous rocks do not indicate that subduction along the SW margin of Eurasia began in Jurassic time.     

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.