Directional interpolation infinite element for dynamic problems in saturated porous media

In this study, a directional interpolation infinite element suited to a saturated porous medium is presented to account for dynamic problems with semi-infi     

A multiscale work-analysis approach for geotechnical structures

This paper presents a second-order work analysis in application to geotechnical problems by using a novel effective multiscale approach. To abandon complicated equations involved in conventional phenomenological models, this multiscale approach employs a micromechanically-based formulation, in which only four parameters are involved. The multiscale approach makes it possible a coupling of the finite element method (FEM) and the micromechanically-based model. The FEM is used to solve the boundary value problem (BVP) while the micromechanically-based model is utilized at the Gauss point of the FEM. Then, the multiscale approach is used to simulate a three-dimensional triaxial test and a plain-strain footing. On the basis of the simulations, material instabilities are analyzed at both mesoscale and global scale. The second-order work criterion is then used to analyze the numerical results. It opens a road to interpret and understand the micromechanisms hiding behind the occurrence of failure in geotechnical issues.     

A machine learning approach for predicting computational intensity and domain decomposition in parallel geoprocessing

ABSTRACT

High performance computing is required for fast geoprocessing of geospatial big data. Using spatial domains to represent computational intensity (CIT) and domain decomposition for parallelism are prominent strategies when designing parallel geoprocessing applications. Traditional domain decomposition is limited in evaluating the computational intensity, which often results in load imbalance and poor parallel performance. From the data science perspective, machine learning from Artificial Intelligence (AI) shows promise for better CIT evaluation. This paper proposes a machine learning approach for predicting computational intensity, followed by an optimized domain decomposition, which divides the spatial domain into balanced subdivisions based on the predicted CIT to achieve better parallel performance. The approach provides a reference framework on how various machine learning methods including feature selection and model training can be used in predicting computational intensity and optimizing parallel geoprocessing against different cases. Some comparative experiments between the approach and traditional methods were performed using the two cases, DEM generation from point clouds and spatial intersection on vector data. The results not only demonstrate the advantage of the approach, but also provide hints on how traditional GIS computation can be improved by the AI machine learning.     

ECOLOGICAL SERIES OF SOIL ANIMALS IN DARLIDAI MOUNTAIN

The ecological series of soil animals under the broad-leaved and pine mixed forest in Darlidai Mountainwas studied. Seven sample plots were selected according to different altitude gradients, which belong to different vegeta-tion types. By investigating and analyzing soil animals in every sample plot it is found that there are 45 groups and 1956individuals, which axe involved in 3 phylums, 7 classes, 16 orders, respectively. The altitude is a key factor which af-fects ecological series of soil animals. Both the groups and individuals of soil animals increase with altitude increasingunder certain conditions, which contrastes with ordinary cases, resulting from special micro-climate in studied area. Thegroups and individuls of soil animals are the most under the broad-leaved and pine forest on the top of the mountain, andthe least under Picea-Abies forest in the foot of the mountain.     

Isotopic evidence for holocene climatic changes in Bosten Lake, Southern Xinjiang, China

Based on the δ¹³C and δ¹⁸O fluctuation of lacustrine carbonate, CaCO₃ content and sporo-pollen data, a palaeoclimatic history of Bosten Lake during the Holocene has been outlined, several stages of climatic changes are divided, and the following result es are obtained: (1) Palaeoclimatic changes revealed by carbonate isotope around Bosten Lake are basically identical with that revealed by other geological records in Xinjiang. Environmental changes presented apparent Westlies Style model: during cold period, relative humidity increased, δ¹⁸O, δ¹³C and CaCO₃ appeared low; but in warm periods, the dry regime aggravated. (2) The temperature reflected by δ¹⁸O exist evident features being increase in the late period during the Holocene. Together with the δ¹³C, pollen and CaCO₃ analyses, several cold and warm phases which are of broad regional significance can be identified. The warm peaks occurred at about 11.0 ka B.P., 9.4 ka B.P., 7.5 ka B.P., 5.0 ka B.P., 3.0 ka B.P. and 2.0 ka B.P.; the cold peaks at 11.5 ka B.P., 10.5 ka B.P., 8.8 ka B.P., 5.5 ka B.P., 3.3 ka B.P., 2.2 ka B.P. and 1.5 ka B.P.. (3) Several climatic events with the nature of “abrupt climatic changes” are revealed in the periods of 11.0 ka B.P. −10.5 ka B.P., 9.4 ka B.P. −8.8 ka B.P., 5.5 ka B.P. −5.0 ka B.P. and 2.0 ka B.P. −1.5 ka B.P.. (4) The results show that carbonate isotopic record of lacustrine sediment in arid area is very sensitive to climatic changes, and may be play a very important role in understanding the features and mechanism of palaeoclimatic changes.     

富氏谱分析法模拟地壳垂直运动的地学依据及其策应力模型

应用地壳波浪与镶嵌构造学说对富氏谱分析法提取地壳垂直形变信息的科学性做了地质学意义上的阐释 ,并提出了根据多期形变资料提取特定波段上构造策应力的数学模型     

One zone modeling of irregular variability of stellar convective envelope

One zone modeling of the irregular variability of red super-giants is intended with regard to the nonlinear coupling of finite amplitude pulsation with convection. The nonlocal mixing length is employed for the evaluation of the convective flux, the turbulent pressure and the turbulent power of temperature fluctuations. The radial pulsation and the Boussinesq convection are assumed for simplicity. The one zone is defined as the layer having the entropy maximum and the minimum at the bottom and at the top, respectively. The quasi-adiabatic approximation is consistent with this definition in fixing the zone to the same mass range. The spatial derivatives are evaluated under the assumption of homologous changes with the equilibrium homologous parameters. Then, a set of 6 simultaneous first order nonlinear ordinary differential equations are obtained as the one zone representation of the irregular variability of the convective envelope.     

The β Cephei instability strip

We carried out a series of linear stability analyses of the radial and low-degree non-radial p modes for stellar models with initial masses of     . The stellar models were computed by using convective overshoot distance     , 0.25 and 0.40  H _P. Our numerical results show that the β Cephei instability strip forms a horn-shaped region pointing upwards near the main sequence on the Hertzsprung–Russell diagram (HRD). The lower part of the instability strip for the radial modes join the zero-age main-sequence (ZAMS) at     , while the top of the instability strip extends up to     . The instability strip for the non-radial modes is even wider. The overall instability strip is dominated by the radial and non-radial fundamental modes. The first overtone (the radial-order index     is also pulsationally unstable. We have shown that the β Cephei stability is almost independent of the overshoot parameter d _over used for the stellar models, while it depends critically on the metal abundance. With decreasing metal abundance, the instability region shrinks and eventually disappears for     .     

Optical Flash of GRB 990123: Constraints on the Physical Parameters of the Reverse Shock

The optical flash accompanying GRB 990123 is believed to be powered by the reverse shock of a thin shell. With the best-fit physical parameters for GRB 990123 and the assumption that the parameters in the optical flash are the same as in the afterglow, we show that: 1) the shell is thick rather than thin, and we have provided the light curve for the thick shell case which coincides with the observation; 2) the theoretical peak flux of the optical flash accounts for only 3×10~-4 of the observed. In order to remove this discrepancy, the physical parameters, the electron energy and magnetic ratios, εe and εB, should be 0.61 and 0.39, which are very different from their values for the late afterglow.     

Silicon limitation on primary production and its destiny in Jiaozhou Bay, China——Ⅳ：Study on cross-bay transect from estuary to ocean

The authors analyzed the data collected in the Ecological Station Jiaozhou Bay from May 1991 to November 1994, including 12 seasonal investigations, to determine the characteristics, dynamic cycles and variation trends of the silicate in the bay. The results indicated that the rivers around Jiaozhou Bay provided abundant supply of silicate to the bay. The silicate concentration there depended on river flow variation. The horizontal variation of silicate concentration on the transect showed that the silicate concentration decreased with distance from shorelines. The vertical variation of it showed that silicate sank and deposited on the sea bottom by phytoplankton uptake and death, and zooplankton excretion. In this way, silicon would endlessly be transferred from terrestrial sources to the sea bottom. The silicon took up by phytoplankton and by other biogeochemical processes led to insufficient silicon supply for phytoplankton growth. In this paper, a 2D dynamic model of river flow versus silicate concentration was established by which silicate concentrations of 0.028–0.062 μmol/L in seawater was yielded by inputting certain seasonal unit river flows (m³/s), or in other words, the silicate supply rate; and when the unit river flow was set to zero, meaning no river input, the silicate concentrations were between 0.05–0.69 μmol/L in the bay. In terms of the silicate supply rate, Jiaozhou Bay was divided into three parts. The division shows a given river flow could generate several different silicon levels in corresponding regions, so as to the silicon-limitation levels to the phytoplankton in these regions. Another dynamic model of river flow versus primary production was set up by which the phytoplankton primary production of 5.21–15.55 (mgC/m²·d)/(m³/s) were obtained in our case at unit river flow values via silicate concentration or primary production conversion rate. Similarly, the values of primary production of 121.98–195.33 (mgC/m²·d) were achieved at zero unit river flow condition. A primary production conversion rate reflects the sensitivity to silicon depletion so as to different phytoplankton primary production and silicon requirements by different phytoplankton assemblages in different marine areas. In addition, the authors differentiated two equations (Eqs. 1 and 2) in the models to obtain the river flow variation that determines the silicate concentration variation, and in turn, the variation of primary production. These results proved further that nutrient silicon is a limiting factor for phytoplankton growth. This study was funded by NSFC (No. 40036010), and the Director's Fund of the Beihai Sea Monitoring Center, the State Oceanic Administration.