Hydrochemical Dynamic Characteristics and Evolution of Underground Brine in the Mahai Salt Lake of the Qaidam Basin Qinghai-Tibet Plateau

The mineral rock salts present in the Mahai Salt Lake of the Qaidam basin exhibit high solubilities in water. In addition, the multicomponent underground brine exhibits a high salinity and is easily precipitated. In the natural state, brine transport in the brine layer is extremely slow, and the brine is in a relatively stable chemical equilibrium state with the rock salt media. However, during mining, both the seepage and the chemical fields fluctuate significantly, thereby disrupting the equilibrium and leading to variations in the chemical composition and dynamic characteristics of the brine. Therefore, we selected underground brine from the Mahai Salt Lake, collecting a total of 183 brine samples over three stages of mining（i.e., the early stage of underground brine extraction, the initial stage of mining, and the later stage of mining）. Using a range of analytical techniques, the chemical dynamics of the underground brine water and its evolution were systematically studied. We found that evaporation and enrichment were the main mechanisms of underground brine evolution in the Mahai Salt Lake, with cation exchange and mineral dissolution/precipitation being key factors in determining the dynamic characteristics and evolution of the brine.

Dynamic evolution of shear rate-dependent behavior of rock discontinuity under shearing condition

Rock blocks sliding along discontinuities can cause serious disasters,such as landslides,earthquakes,or rock bursts.The shear rate-dependent behavior is a typical time-dependent behavior of a rock discontinuity,and it is closely related to the stability of a rock block.To further study the shear rate-dependent behavior of rock discontinuities,shear tests with alternating shear rates(SASRs)were conducted on rock discontinuities with various surface morphologies.The dynamic evolution of the shear rate dependency was studied in detail based on the shear test results,and three stages were identified with respect to the shear stress and shear deformation states.The test results revealed that dynamic changes in shear stiffness and the energy storage abilities of the rock discontinuities occurred in relation to the shear rate-dependent behavior of crack growth,which increased with an increase in normal stress and/or the joint roughness coefficient.The stage of decreasing shear stiffness corresponded to a stage of noticeable shear rate-dependency,and the shear rate was found to have no influence on the initial crack stress.

Spatiotemporal evolution of ecosystem service value and topographic gradient effect in the Da-Xiao Liangshan Mountains in Sichuan Province,China

The Da-Xiao Liangshan mountains are critical ecological function areas and essential ecological barriers in the upper reaches of the Yangtze River in China.This study selected a total of six periods of land use land cover(LULC)data from 1995 to 2020,to estimate the ecosystem service value(ESV)and analyzed its spatiotemporal evolution and topographic gradient divergence.The results showed that:(1)The ESV increased by 1.1 billion yuan,with an increase rate of 1.47%from 1995 to 2020.Two time periods,2005–2010 and 2015–2020,showed more significant increases than other periods.(2)The elevation and slope of mountainous areas determine the type of land use and further influence the spatial pattern of ESV.(3)Although woodland and grassland are the main land use types of the study area(more than 90%),the hydrological regulation function of the water area partially compensated for the impact of the encroachment of the built-up area on the ESV of grassland.(4)The spatial distribution of ESVs showed an inverted V-shaped characteristic as the topographic gradient increased,with the dominant position being the 5th topographic gradient zone.Finally,this study provided relevant recommendations for ecosystem protection and optimization.The findings of this study clarified the influence of topographical factors on the spatial differentiation of ESV and provided novel insights into ecosystem protection.

Numerical simulation of wheel wear evolution for heavy haul railway

The prediction of the wheel wear is a fundamental problem in heavy haul railway. A numerical methodology is introduced to simulate the wheel wear evolution of heavy haul freight car. The methodology includes the spatial coupling dynamics of vehicle and track, the three-dimensional rolling contact analysis of wheel-rail, the Specht's material wear model, and the strategy for reproducing the actual operation conditions of railway. The freight vehicle is treated as a full 3D rigid multi-body model. Every component is built detailedly and various contact interactions between parts are accurately simulated, taking into account the real clearances. The wheel-rail rolling contact calculation is carried out based on Hertz's theory and Kalker's FASTSIM algorithm. The track model is built based on field measurements. The material loss due to wear is evaluated according to the Specht's model in which the wear coefficient varies with the wear intensity. In order to exactly reproduce the actual operating conditions of railway,dynamic simulations are performed separately for all possible track conditions and running velocities in each iterative step.Dimensionless weight coefficients are introduced that determine the ratios of different cases and are obtained through site survey. For the wheel profile updating, an adaptive step strategy based on the wear depth is introduced, which can effectively improve the reliability and stability of numerical calculation. At last, the wear evolution laws are studied by the numerical model for different wheels of heavy haul freight vehicle running in curves. The results show that the wear of the front wheelset is more serious than that of the rear wheelset for one bogie, and the difference is more obvious for the outer wheels. The wear of the outer wheels is severer than that of the inner wheels. The wear of outer wheels mainly distributes near the flange and the root; while the wear of inner wheels mainly distributes around the nominal rolling circle. For the outer wheel of front wheelset of each bogie, the development of wear is gradually concentrated on the flange and the developing speed increases continually with the increase of traveled distance.

Dynamic evolution patterns of the degree of ground resistivity anisotropy and the sets mogenic process

Following a new train of thinking, this paper has explored first the potential information in the ground resistivitydata observed by the existing geoelectric observation system, investigated and proposed a new dimensionlessgeoelectric precursor factor, the degree of ground resistivity anisotropy, S, and studied the characteristics of dynamic evolution pattern of S during the seismogenic process. The results show that, during the seismogenic process, the degree of ground resistivity anisotropy (S) displays h process of 'normal' → 'abnormal strengthening(amplitude, range)' → 'abnormal weakening'→ 'earthquake occurrence'→ 'normal'. The earthquake wouldoccur at the time when the S value has entered the late stage of strengthening and turns to weaken and in the gradient belt on the margin ofS anomaly region. The dynamic evolution pattern ofS reflects the changes of the tectonicstress field during the seismogenic process. Therefore, it would be possible to trace the process of earthquake generation and occurrence from the dynamic evolution pattern ofS so as to service eaJ'thquake prediction.

A semi-analytic model for the study of 1/1 resonant dynamics of the planar elliptic restricted co-orbital problem

Mean motion resonances(MMRs)are widespread in our Solar System.Moreover,resonant dynamics has always been an essential topic in planetary research.Recently,the research about exoplanets and the potential Planet Nine with large eccentricity has given rise to our interests in the secular dynamics inside MMRs in the elliptic model.In this paper,we study the fixed points of the averaged Hamiltonian and the long-term stable regions of the 1/1 resonance(or co-orbital motion)in the elliptic restricted three-body problem(ERTBP)systematically.Numerical integrations prove those test particles trapped in"apsidal co-rotation",where both the resonant angleφ_(res)and the secular angleΔω(or apsidal longitude differences)librate simultaneously,always survive the long-term simulations.Furthermore,utilizing a semianalytical method based on the adiabatic approach,three families of long-term fixed points of the averaged Hamiltonian of the planar ERTBP inside the 1/1 resonance have been found.We call them QS-points,Hpoints,and T-points here,whose values of the(φ_(res),Δω)are(0°,180°),(180°,0°),and(±60°,±60°),respectively.All the fixed points of the averaged Hamiltonian of the co-orbital motion in the ERTBP are presented in the e-e’plane(’represents the elements of the planet in this paper).We find that QS-points and T-points always exist for the arbitrary eccentricity of a planet,while H-points only exist for the cases of low e’and very high e.Furthermore,we measure the libration width in terms of eccentricity,Δe,around these stable equilibrium points in the e-Δωphase-space portraits.The"apsidal co-rotation"around all the stable equilibrium points is presented in the e-e’plane.All these results are effectively confirmed by numerical experiments.The long-term stable zones around these periodic orbits in the e-e’plane are significant for the research of the co-orbital motion in the ERTBP.Above all,these practical approaches that we proposed can also be used to study the secular dynamics of other MMRs.

Solving chemical dynamic optimization problems with ranking-based differential evolution algorithms

Dynamic optimization problems(DOPs) described by differential equations are often encountered in chemical engineering. Deterministic techniques based on mathematic programming become invalid when the models are non-differentiable or explicit mathematical descriptions do not exist. Recently, evolutionary algorithms are gaining popularity for DOPs as they can be used as robust alternatives when the deterministic techniques are invalid. In this article, a technology named ranking-based mutation operator(RMO) is presented to enhance the previous differential evolution(DE) algorithms to solve DOPs using control vector parameterization. In the RMO, better individuals have higher probabilities to produce offspring, which is helpful for the performance enhancement of DE algorithms. Three DE-RMO algorithms are designed by incorporating the RMO. The three DE-RMO algorithms and their three original DE algorithms are applied to solve four constrained DOPs from the literature. Our simulation results indicate that DE-RMO algorithms exhibit better performance than previous non-ranking DE algorithms and other four evolutionary algorithms.

Microstructural evolution of a superaustenitic stainless steel during a two-step deformation process

Single-and two-step hot compression experiments were carried out on 16Cr25Ni6Mo superaustenitic stainless steel in the temperature range from 950 to 1150°C and at a strain rate of 0.1 s^（-1）. In the two-step tests, the first pass was interrupted at a strain of 0.2; after an interpass time of 5, 20, 40, 60, or 80 s, the test was resumed. The progress of dynamic recrystallization at the interruption strain was less than 10%. The static softening in the interpass period increased with increasing deformation temperature and increasing interpass time. The static recrystallization was found to be responsible for fast static softening in the temperature range from 950 to 1050°C. However, the gentle static softening at 1100 and 1150°C was attributed to the combination of static and metadynamic recrystallizations. The correlation between calculated fractional softening and microstructural observations showed that approximately 30% of interpass softening could be attributed to the static recovery. The microstructural observations illustrated the formation of fine recrystallized grains at the grain boundaries at longer interpass time. The Avrami kinetics equation was used to establish a relationship between the fractional softening and the interpass period. The activation energy for static softening was determined as 276 kJ/mol.

Sand barrier morphological evolution based on time series remote sensing images:a case study of Anhaiao,Pingtan

The morphological evolution of the sand barrier in the Anhaiao coastal zone of Pingtan from 1996 to 2018 was studied.Tidal correction was used to refine the location of the coastline.A standard deviation ellipse method was applied to further analyze the movement of the barrier head with the axis and rotation angle.A natural neighbor interpolation(NNI)method was carried out to calculate the terrain of the intertidal area,and the erosion and deposition characteristics were illustrated based on the terrain.The results showed that the northern part of the sand barrier facing the lagoon area was deposited over the whole studied period,while erosion has always occurred in the southern part of the sand barrier facing the open sea.The erosion and deposition were slightly different on both sides of the barrier head due to hydrodynamic turbulence.The middle sand barrier moved 102.60 m away from its original location in 1996,and the end of the barrier moved 65.45 m.The head of the sand barrier continued moving 379 m to the northwest.Consequently,the preliminary morphological evolution of the sand barrier corresponding to the distance and direction of movement was detected.

Features and physical proces of the dynamic evolution pattern of ground resistivity precursor front

In order to quantitatively describe the geoelectric precursor anomaly in the short-impending process of earthquakes, a new geoelectric precursor index — (monthly) relative change rate of ground resistivity, R ρ (t) , is designed. Using this index and choosing the internationally accepted ground resistivity data before the Tangshan M =7.8 earthquake of July 28, 1976, the features of dynamic evolution pattern of R ρ(t) are studied. The results show that: ① about 10～9 months before earthquake, the ground resistivity in a certain range around the epicentral region begins to display the anomaly of accelerating descent, and the rate of descent is higher in the epicentral region than in surrounding areas; ② with the shortening of countdown before earthquake, the R ρ(t) value in epicentral region increases gradually (ground resistivity value decreases at an increasing rate); ③ the R ρ(t) value has the epicentral area as a center and its contour lines propagate towards surrounding areas with the shortening of countdown before earthquake; ④after the R ρ(t) value in epicentral region has descended at increasing rate to reach an extremity [ R ρ(t) = (7.0], it turns to descend at decreasing rate (2～3 months) and earthquake occurs when it accelerates again. Meanwhile, earthquake occurs when the contour lines of R ρ(t) stop propagating towards surrounding areas and turn to shrink back (2～3 months later). Its physical process can be explained by the″ swollen hypothesis″ of Prof. Fu and the theory of ″Slip-weakening and rockmass instability″ of Mei, Niu, et al ..

Time Evolution,Dynamical Quantum Fluctuation and High-Order Squeezing Feature in Polariton System——Ⅰ

We have set up a new reduced model Hamiltonian for the polariton system, in which the nonlinear interaction contains the rotating term k l （ a ＋ b ＋ ab＋） and the attractive two-mode squeezed coupling - k2 （ a ＋ b＋ ＋ ab ） . The dynamical evolution of this system has been solved and the nonclassical features relevant to the second-order and high-order squeezing have been obtained in an analytical form. For the first time, in contrast to the existing result, we have confirmed for the phonon field that the attractive two-mode squeezed interaction will not only result in the second-order and high-order squeezing in X-component with the time evolution, but also in time average. Furthermore, the phenomena of collapse and revival of inversion will occur as well in the time evolution of the average number of photon and phonon, as also in the second-order and high-order squeezing of photon field, particularly, in the high-order squeezing of phonon field.

Dynamic aggregation evolution of competitive societies of cooperative and noncooperative agents

We propose an evolution model of cooperative agent and noncooperative agent aggregates to investigate the dynamic evolution behaviors of the system and the effects of the competing microscopic reactions on the dynamic evolution. In this model, each cooperative agent and noncooperative agent are endowed with integer values of cooperative spirits and nonco- operative spirits, respectively. The cooperative spirits of a cooperative agent aggregate and the noncooperative spirits of a noncooperative agent aggregate change via four competing microscopic reaction schemes： the win-win reaction between two cooperative agents, the lose-lose reaction between two noncooperative agents, the win-lose reaction between a coop- erative agent and a noncooperative agent （equivalent to the migration of spirits from cooperative agents to noncooperative agents）, and the cooperative agent catalyzed decline of noncooperative spirits. Based on the generalized Smoluchowski＇s rate equation approach, we investigate the dynamic evolution behaviors such as the total cooperative spirits of all coop- erative agents and the total noncooperative spirits of all noncooperative agents. The effects of the three main groups of competition on the dynamic evolution are revealed. These include： （i） the competition between the lose-lose reaction and the win-lose reaction, which gives rise to respectively the decrease and increase in the noncooperative agent spirits; （ii） the competition between the win-win reaction and the win-lose reaction, which gives rise to respectively the increase and decrease in the cooperative agent spirits; （iii） the competition between the win-lose reaction and the catalyzed-decline reaction, which gives rise to respectively the increase and decrease in the noncooperative agent spirits.

The Investigation of the Dynamical Evolution of Extrasolar Three-planetary System GJ 3138

This article is devoted to studying the dynamical evolution and orbital stability of compact extrasolar threeplanetary system GJ 3138. In this system, all semimajor axes are less than 0.7 au. The modeling of planetary motion is performed using the averaged semi-analytical motion theory of the second order in planetary masses,which the authors construct. Unknown and known with errors orbital elements vary in allowable limits to obtain a set of initial conditions. Each of these initial conditions is applied for the modeling of planetary motion. The assumption about the stability of observed planetary systems allows to eliminate the initial conditions leading to excessive growth of the orbital eccentricities and inclinations and to identify those under which these orbital elements conserve moderate values over the whole modeling interval. Thus, it becomes possible to limit the range of possible values of unknown orbital elements and determine their most probable values in terms of stability.

Stellar Abundance and Galactic Chemical Evolution through LAMOST Spectroscopic Survey

A project of a spectroscopic survey of Galactic structure and evolution with a Large sky Area Multi-Object fiber Spectroscopic Telescope （LAMOST） is presented. The spectroscopic survey consists of two observational modes for various targets in our Galaxy. One is a major survey of the Milky Way aimed at a systematic study of the stellar abundance and Galactic chemical evolution through low resolution （R = 1000 2000） spectroscopy. Another is a follow-up observation with medium resolution （R = 10000） spectrographs aimed at detailed studies of the selected stars with different chemical composition, kinematics and dynamics.

Regional difference and dynamic evolution of development quality of power industry in China

To achieve the goals of carbon peaking and carbon neutrality and maintain high-quality economic growth,China is currently striving to improve the quality of development of its power sector.In this regard,revealing the regional differences and evolutionary trends in the development quality of China’power sector has a high value to inspire the next improvement direction toward how to integrate regional power recourses to an overall optimization level.Motived by this purpose,this paper uses the entropy method to evaluate the com‐prehensive and subsystem indices of the development quality of the power industry,and reveals their re‐gional differences and evolutionary trends with the help of the Dagum Gini coefficient and Kernel density es‐timation methods.The findings show that:There are obvious regional differences in the development quality of China’s power industry,and the differences are steadily declining in all regions except the West.Regional differences are mainly derived from inter-regional differences,with the largest inter-regional differences in the East-Northeast region.Intra-regional differences show a distribution pattern of East>West>North‐east>Center.

A Dynamic Evolution Model of Airline Networks

Empirical data show that most of the degree distribution of airline networks assume a double power law. In this work, firstly, we assume cities as sites, flight between two cities as an edge between two sites, and build a dynamic evolution model for airline networks by improving the BA model, in which the conception of attractiveness plays a decisive role in the course of evolution of the networks. To this end, we discuss whether the attractiveness depends on the site label s or not separately, finally we obtain analytic degree distribution. As a result, if the attractiveness of a site is independent of the degree distribution of sites, which will follow the double power law, otherwise, it will be scale-free. Moreover, degree distribution depends on the parameters of the models, and some parameters aye more sensitive than others.

Dynamical Evolution of an Effective Two-Level System with PT Symmetry

We investigate the dynamics of parity-and time-reversal（PT） symmetric two-energy-level atoms in the presence of two optical and one radio-frequency fields. The strength and relative phase of fields can drive the system from the unbroken to the broken PT symmetric regions. Compared with the Hermitian model, Rabi-type oscillation is still observed, and the oscillation characteristics are also adjusted by the strength and relative phase in the region of the unbroken symmetry. At the exception point, the oscillation breaks down. To better understand the underlying properties we study the effective Bloch dynamics and find that the emergence of the PT components of the fixed points is the feature of the PT symmetry breaking and the projections in the x–y plane can be controlled with high flexibility compared with the standard two-level system with the PT symmetry. It helps to study the dynamic behavior of the complex PT symmetric model.

Collaboration strategy for software dynamic evolution of multi-agent system

As the ability of a single agent is limited while information and resources in multi-agent systems are distributed, cooperation is necessary for agents to accomplish a complex task. In the open and changeable environment on the Internet, it is of great significance to research a system flexible and capable in dynamic evolution that can find a collaboration method for agents which can be used in dynamic evolution process. With such a method, agents accomplish tasks for an overall target and at the same time, the collaborative relationship of agents can be adjusted with the change of environment. A method of task decomposition and collaboration of agents by improved contract net protocol is introduced. Finally, analysis on the result of the experiments is performed to verify the improved contract net protocol can greatly increase the efficiency of communication and collaboration in multi-agent system.

Scale effect of coastal landscape pattern stability and driving forces:a case study of Guangdong Province,China

The long-term dynamic evolution and underlying mechanisms of coastal landscape pattern stability,driven by strong anthropogenic interference and consequently climate change,are topics of major interest in national and international scientific research.Guangdong Province,located in southeastern China,has been undergoing rapid urbanization over several decades.In this study,we quantitatively determined the scale threshold characteristics of coastal landscape pattern stability in Guangdong Province,from the dual perspective of spatial heterogeneity and spatial autocorrelation.An analysis of the spatiotemporal evolution of the coastal landscape was conducted after the optical scale was determined.Then,we applied the geodetector statistical method to quantitatively explore the mechanisms underlying coastal landscape pattern stability.Based on the inflection point of landscape metrics and the maximum value of the MoranⅠindex,the optimal scale for analyzing coastal landscape pattern stability in Guangdong Province was 240 m×240 m.Within the past several decades,coastal landscape pattern stability increased slightly and then decreased,with a turning point around 2005.The most significant variations in coastal landscape pattern stability were observed in the transition zone of rural-urban expansion.A q-statistics analysis showed that the explanatory power of paired factors was greater than that of a single driving factor;the paired factors with the greatest impact on coastal landscape pattern stability in Guangdong Province were the change in gross industrial output and change in average annual precipitation from 2010 to 2015,based on a q value of 0.604.These results will contribute to future efforts to achieve sustainable coastal development and provide a scientific basis and technical support for the rational planning and utilization of resources in large estuarine areas,including marine disaster prevention and seawall ecological restoration.

Capture Efficiency Analysis in the Circular Restricted Three-body Problem

Temporary capture efficiency is studied in the framework of the circular restricted three-body problem in two steps.First, a non-uniform distribution of test particles around the secondary's orbit is obtained by fully accounting the secondary's gravitational influence. Second, the capture efficiency is computed based on the non-uniform distribution. Several factors influencing the result are discussed. By studying the capture efficiency in the circular restricted three-body problem of different mass ratios, a power-law relation between the capture efficiency(p) and the mass ratio(μ) is established, which is given by p ≈ 0.27 × μ^(0.53), within the range of 3.0035 ×10^(-6)≤ μ ≤ 3.0034 × 10^(-5). Taking the Sun–Earth system as an example, the influence from the orbit eccentricity of the secondary on the non-uniform distribution and the capture efficiency is studied. Our studies find that the secondary's orbit eccentricity has a negative influence on the capture efficiency.