A new era of flood control strategies from the perspective of managing the 2020 Yangtze River flood

Flood control of the Yangtze River is an important part of China’s national water security.In July 2020,due to continuous heavy rainfall,the water levels along the middle-lower reaches of the Yangtze River and major lakes constantly exceeded the warning levels,in which Taihu Lake exceeded its highest safety water level and some stations of Poyang Lake reached their highest water levels in its history.In August 2020,another huge flood occurred in the Minjiang River and the Jialing River in the upper Yangtze River,and some areas of Chongqing Municipality and other cities along the rivers were inundated,resulting in great pressure on flood control and high disaster losses.The 2020 Yangtze River flood has received extensive media coverage and raised concerns on the roles of the Three Gorges Dam and other large reservoirs in flood control.Here we analyze the changes in the pattern of the Yangtze River flood control by comparing the strategies to tackle the three heavy floods occurring in 1954,1998,and 2020.We propose that the overall strategy of the Yangtze River flood control in the new era should adhere to the principle of"Integration of storage and drainage over the entire Yangtze River Basin,with draining floods downstream as the first priority"by using both engineering and non-engineering measures.On the basis of embankments,the engineering measures should use the Three Gorges Dam and other large reservoirs as the major regulatory means,promote the construction of key flood detention areas,keep the floodways clear,and maintain the ecosystem services of wetlands and shoals.In terms of non-engineering measures,we should strengthen adaptive flood risk management under climate change,standardize the use of lands in flood detention areas,give space to floods,and promote the implementation of flood risk maps and flood insurance policies.The ultimate goal of this new flood control system is to enhance the adaptability to frequent floods and increase the resilience to extreme flood disasters.     

The soil Microbial Carbon Pump as a new concept for terrestrial carbon sequestration

Soil is a huge terrestrial carbon pool, which has higher carbon storage than the sum of atmospheric and terrestrial vegetation carbon. Small fluctuations in soil carbon pool can affect regional carbon flux and global climate change. As soil organic carbon plays key roles in soil carbon storage and sequestration, studying its composition, sources and stability mechanism is a key to deeply understand the functions of terrestrial ecosystem and how it will respond to climate changes. The recently-proposed concept of soil Microbial Carbon Pump(MCP) emphasizes the importance of soil microbial anabolism and its contributions to soil carbon formation and stabilization, which can be applied for elucidating the source, formation and sequestration of soil organic carbon. This article elaborates MCP-mediated soil carbon sequestration mechanism and its influencing factors, as well as representative scientific questions we may explore with the soil MCP conceptual framework.     

Geophysical evidence for thrusting of crustal materials from orogenic belts over both sides of the Yangtze Block and its geological significance

Based on deep geophysical detections, we have reconstructed the crustal structure from the eastern margin of the Tibetan Plateau to the Jiangnan-Xuefeng orogenic belt. The results suggest that the Yangtze Block was overthrusted by crustal materials in its NW direction from the eastern Tibetan Plateau but in its SE direction from the Jiangnan orogen. These overthrusting effects control the crustal structure from the western Sichuan to the western area of the Jiangnan orogen-Xuefeng orogenic belt. The eastward extruded materials from the eastern Tibetan Plateau were blocked by the rigid basement in the Sichuan Basin, where upper-middle crust was overthrusted whereas the lower crust was underthrusted beneath the Sichuan Basin. The underthrusted unit was absorbed by crustal folding, shortening and thickening in the Yangtze Block, forming the Xiongpo and Longquan Mountains tectonic belts and resulting in the NW-directed thrusting of the Pujiang-Chengdu-Deyang fault, and the western hillsiden fault in the Longquan Mountain. These results provide resolution to the controversy where the eastward extrusion material from the Qinghai-Tibet Plateau had gone. Overall, that Yangtze Block was subjected to thrusting of the crustal materials from the orogenic belts over its both sides. This finding has implications for the study of the intracontinental orogenic mechanism in South China, the reconstruction of tectonic evolutionary history and the kinematics processes during the lateral extrusion of the Tibet Plateau.     

Variations of widespread extreme cold and warm days in winter over China and their possible causes

The two leading modes of winter surface air temperature(SAT) over China during 1961–2017 are a spatially consistent pattern and a north-south dipole pattern. Based on the two leading modes, the characteristics of the extreme cold and warm days in the two patterns, defined by the standard deviation larger than 1.28 or smaller than-1.28 in the time series of the two leading modes, are analyzed. With the increase of winter SAT during 1961–2017, the number of spatially consistent extreme cold days decreased and their occurrence was restricted to late December to early January, whereas the number of spatially consistent extreme warm days increased significantly in January and February. Global warming is associated with an increase in the spatially consistent extreme warm days and a decrease in spatially consistent extreme cold days, but has little relation to the sum of extreme cold and warm days of either the spatially consistent or north-south dipole pattern. The Siberian High(SH) is the main factor controlling the sum of spatially consistent extreme warm and cold days. The strong(weak) SH before(after) the1990 s corresponds to an increase(decrease) in the sum of the spatially consistent extreme warm and cold days. The occurrences of extreme south-cold-north-warm and extreme south-warm-north-cold days are related to the north-south difference of the SH.When the center of the SH is in mid-high latitudes, the extreme south-warm-north-cold(south-cold-north-warm) days occur more(less) often. During the winters of 1961–2017, the total number of extreme cold and warm days of the north-south dipole pattern changes negligibly. The North Atlantic meridional overturning circulation(AMOC) may be the main factor affecting the sum of the extreme cold and warm days of the two types of SAT pattern in China.     

Existence,morphology and structure of the Yellow Sea Warm Current Branch approaching waters offshore Qingdao,China

In this study, we focused on full-region cruise survey data, near-bottom continuous mooring observations and sea surface wind products from the western South Yellow Sea in winter; after ensuring the data reliability and accuracy, we processed and analyzed the data. Image resolution experiments were carried out to determine the lowest recognition resolutions for all image types, which represent the resolution characteristics of the data. The existence of a warm water tongue originating from the Yellow Sea Warm Current(YSWC) that approached waters offshore Qingdao was confirmed. For the first time, a high salinity water tongue, corresponding to the warm water tongue, was described and found to be more representative of the YSWC branch path. This warm tongue is a sign of the branch originating from the YSWC, which we defined as the Yellow Sea Warm Current Branch approaching waters offshore Qingdao(YSWC-QDB). The pattern of the warm and salty water tongues showed remarkable rear, branching middle, shrinking neck and expanding top regions. These patterns showed a temporal feature of the tongues, and were the result of multi-temporal branches in front of the YSWC main section as well as the YSWC-QDB crossing the southwestward path of the extension of the North Shandong Coastal Current flowing along the southeast coast of the Shandong Peninsula(NSCC-SESE). Analysis using mooring data at a sensitive and representative station also showed the existence of the YSWC-QDB. It is a probabilistic event that manifests as a northwestward flow that decreases gradually from the bottom to the surface in the early cold air transit stage and consistent in the whole water column profile in the later stage. It varies quasi-periodically with weather processes. It also transports some of the YSWC water stored in the entrance area of the Bohai and Yellow seas under winter wind conditions to the western South Yellow Sea as a compensatory current. This current, caused by northerly winds, especially northwest winds and obstruction of the NSCC-SESE, was present, and strong water reduction and compensation caused significant residual sea level oscillations. The compensatory current, if caused by strong northwest wind,began to appear when its direction was opposite to the wind direction. In addition, confirmation of the YSWC-QDB provides an oceanographic basis for the short cooling time and rapid warming in the Qingdao area in winter. This research provides a basis for further studies of the YSWC-QDB at high spatial and temporal resolutions using large sea surface datasets. For monsoon basin dynamics, this study can also be extended to the whole Bohai and Yellow seas and closed or semi-closed basins on the continental margin.     

Review of Mercury's dynamic magnetosphere:Post-MESSENGER era and comparative magnetospheres

This review paper summarizes the research of Mercury’s magnetosphere in the Post-MESSENGER era and compares its dynamics to those in other planetary magnetospheres,especially to those in Earth’s magnetosphere.This review starts by introducing the planet Mercury,including its interplanetary environment,magnetosphere,exosphere,and conducting core.The frequent and intense magnetic reconnection on the dayside magnetopause,which is represented by the flux transfer event"shower",is reviewed on how they depend on magnetosheath plasma β and magnetic shear angle across the magnetopause,followed by how it contributes to the flux circulation and magnetosphere-surface-exosphere coupling.In the next,Mercury’s magnetosphere under extreme solar events,including the core induction and the reconnection erosion on the dayside magnetosphere,as well as the responses of the nightside magnetosphere,are reviewed.Then,the dawn-dusk properties of the plasma sheet,including the features of the ions,the structure of the current sheet,and the dynamics of magnetic reconnection,are summarized.The last topic is devoted to the particle energization in Mercury’s magnetosphere,which includes the energization of the Kelvin-Helmholtz waves on the magnetopause boundaries,reconnection-generated magnetic structures,and the cross-tail electric field.In each chapter,the last section discusses the open questions related to each topic,which can be considered by the simulations and the future spacecraft mission.We end this paper by summarizing the future Bepi Colombo opportunities,which is a joint mission of ESA and JAXA and is en route to Mercury.     

Detailed mapping of the surface rupture of the 12 February 2014 Yutian Ms7.3 earthquake,Altyn Tagh fault,Xinjiang, China

Large-scale detailed mapping plays a key role in revealing the rupture characteristics and mechanisms of strong earthquakes.Relatively few studies have been performed on the surface ruptures of large earthquakes in central and western Tibet due to its remote nature and high elevation.Based on high-resolution unmanned aerial vehicle(UAV)photography,we mapped the coseismic surface rupture of the 2014 Yutian M_s7.3 earthquake.Along the western Altyn Tagh fault system,the earthquake produced~37 km of surface rupture along the South Xor Kol fault(southern section S1),Xor Kol fault(central section S2)and Ashikule fault(northern section S3).Section S1 has a 16-km-long surface rupture with an average sinistral offset of 52±25 cm and a maximum offset of~90 cm,while section S3 has a 14.2-km-long surface rupture with an average sinistral offset of 36±21 cm and a maximum offset of~84 cm.A compilation of 5308 cracks yields an average crack width along the southern section of 85±71 cm and a maximum width of~700 cm;the average width along the central section is 39±21 cm,and the maximum width is 243 cm;and the average width along the northern section is 61±44 cm with a maximum of~340 cm.In addition,the average cumulative opening across rupture zone is 3.4±2.9 m along the southern section,with a maximum of~17 m;4.3±3.6 m along the central section,with a maximum of~13 m;and 1.7±1.6 m along the northern section,with a maximum of~6 m.Evidently,the average crack width and cumulative opening decrease towards bends and steps along the fault.A global synthesis of surface rupture distributions corresponding to strike-slip earthquakes indicates that the rupture zone is wider near the complex parts of fault geometries(such as bends,steps and fault bifurcations)than along straight sections,suggesting that the fault geometry has an obvious control on the surface rupture width.The widespread cracks at the intersection between the Xor Kol and South Xor Kol faults may indicate that an extensional regime is more likely to produce distributed offfault deformation,which provides an observational constraint for the numerical simulation of dynamic rupture on a fault.In addition to coseismic surface rupture,the Yutian earthquake also produced a large number of gravity-driven slides on alluvial fans with gentle slopes.The friction efficiency of the water-bearing salt layer beneath fans could decrease the sliding threshold and trigger instability under surface shaking.These distributed deformations and gravity-driven slides reflect the coupling between the rupture propagation and fault geometry and indicate that the rupture may have propagated in two directions along the Ashikule fault after passing through a step.Therefore,the investigation of coseismic surface rupture provides important observational constraints on the dynamic rupture process.     

Coupled and splitting bedload sediment transport models based on a modified flux-wave approach

Numerical modeling of free-surface flow over a mobile bed with predominantly bedload sediment transport can be done by solving the shallow water and Exner equations using coupled and splitting approaches.The coupled method uses a coupling of the governing equations at the same time step leading to a non-conservative solution.The splitting method solves the Exner and the shallow water equations in a separate manner,and is only capable of modeling weak free-surface and bedload interactions.In the current study,an extended version of a Godunov-type wave propagation algorithm is presented for modeling of morphodynamic systems using both coupled and splitting approaches.In the introduced coupled method the entire morphodynamic system is solved in the form of a conservation law.For the splitting technique,a new wave Riemann decomposition is defined which enables the scheme to be utilized for mild and strong interactions.To consider the bedload sediment discharge within the Exner equation,the Smart and Meyer-Peter&Müller formulae are used.It was found that the coupled solution gives accurate predictions for all investigated flow regimes including propagation over a dry-state using a Courant-Friedrichs-Lewy(CFL)number equal to 0.6.Furthermore,the splitting method was able to model all flow regimes with a lower CFL number of 0.3.     

Numerical simulations of shake-table experiment for dynamic soil-pile-structure interaction in liquefiable soils

A shake-table experiment on pile foundations in liquefi able soils composed of liquefi able sand and overlying soft clay is studied. A three-dimensional(3D) effective stress fi nite element(FE) analysis is employed to simulate the experiment. A recently developed multi-surface elasto-plastic constitutive model and a fully coupled dynamic inelastic FE formulation(u-p) are used to model the liquefaction behavior of the sand. The soil domains are discretized using a solid-fl uid fully coupled(u-p) 20-8 noded brick element. The pile is simulated using beam-column elements. Upon careful calibration, very good agreement is obtained between the computed and the measured dynamic behavior of the ground and the pile. A parametric analysis is also conducted on the model to investigate the effect of pile-pinning, pile diameter, pile stiffness, ground inclination angle, superstructure mass and pile head restraints on the ground improvement. It is found that the pile foundation has a noticeable pinning effect that reduces the lateral soil displacement. It is observed that a larger pile diameter and fi xed pile head restraints contribute to decreasing the lateral pile deformation; however, a higher ground inclination angle tends to increase the lateral pile head displacements and pile stiffness, and superstructure mass seems to effectively infl uence the lateral pile displacements.     

Recent progress of seismic research on tall buildings in China Mainland

To study the seismic behavior of high strength concrete fi lled double-tube(CFDT) columns,each consisting of an external square steel tube and an internal circular steel tube,quasi-static tests on eight CFDT column specimens were conducted.The test variables included the width-to-thickness ratio(β1) and the area ratio(β2) of the square steel tube,the wall thickness of the circular steel tube,and the axial force(or the axial force ratio) applied to the CFDT columns.The test results indicate that for CFDT columns with a square steel tube with β1 of 50.1 and 24.5,local buckling of the specimen was found at a drift ratio of 1/150 and 1/50,respectively.The lateral force-displacement hysteretic loops of all specimens were plump and stable.Reducing the width-to-thickness ratio of the square steel tube,increasing its area ratio,or increasing the wall thickness of the internal circular steel tube,led to an increased fl exural strength and deformation capacity of the specimens.Increasing the design value of the axial force ratio from 0.8 to 1.0 may increase the fl exural strength of the specimens,while it may also decrease the ultimate deformation capacity of the specimen with β1 of 50.1.     

Natural gas leakage of Mizhi gas reservoir in Ordos Basin, recorded by natural gas fluid inclusion

Abundant natural gas inclusions were found in calcite veins filled in fractures of Central Fault Belt across the centre of Ordos Basin. Time of the calcite veins and characteristics of natural gas fluid inclusion were investigated by means of dating of thermolum luminescence (TL) and analyzing stable isotope of fluid inclusion. Results show that natural gas inclusion formed at 130―140℃ with salinity of 5.5 wt%―6.0 wt% NaCl. It indicates that natural gas inclusion is a kind of thermal hydrocarbon fluid formed within the basin. Method of opening inclusion by heating was used to analyze composition of fluid inclusion online, of which the maximal hydrocarbon gas content of fluid inclusion contained in veins is 2.4219 m3/t rock and the maximal C1/Σci ratio is 91%. Laser Raman spectroscopy (LRS) was used to analyze chemistry of individual fluid inclusion in which the maximal hydrocarbon gas content is 91.6% compared with little inorganic composition. Isotope analysis results of calcite veins show that they were deposited in fresh water, in which the δ13CPDB of calcite veins is from -5.75‰ to 15.23‰ andδ18OSMOW of calcite veins is from 21.33‰ to 21.67‰. Isotope results show thatδ13C1 PDB of natural gas fluid inclusion is from -21.36‰ to -29.06‰ and δDSMOW of that is from -70.89‰ to -111.03‰. It indicates that the gas of fluid inclusion formed from coal source rocks and it is the same as that of natural gas of Mizhi gas reservoir. Results of TL dating show that time of calcite vein is (32.4±3.42)×104 a, which is thought to be formation time of gas inclusion. It indicated that natural gas inclusion contained in calcite veins recorded natural gas leakage from Mizhi gas reservoir through the Central Fault Belt due to Himalayan tectonic movement.     

SP-3400气相色谱仪与GC-2400便携式色谱仪对比观测

气相色谱法具有分离效率高、分析速度快、分析灵敏度高等特点,基于SP-3400气相色谱仪及GC-2400便携式色谱仪的工作原理、工作方式、工作条件等,对庐江地震台2套色谱仪观测结果进行对比,分析存在的误差,通过H2高值异常现象,探讨色谱仪运行及维护工作。初步分析结果表明:2套色谱仪定量分析精度均较高,观测数据相关性较好,因不同仪器生产厂家所提供的色谱工作站软件不同,部分测项观测结果存在一定差异。     

Natural gas origins of large and medium-scale gas fields in China sedimentary basins

China sedimentary basins present abundant natural gas resource thanks to its unique geological settings.Marine highly-matured hydrocarbon source rocks,widespread coal-measure strata and low temperature Quaternary saline strata,etc.,indicate the wide foreground of China natural gas resources. Up to now,most of the petroliferous basins have been discovered to have wholesale natural gas accumulation from Precambrian,Paleozoic,Mesozoic to Cenozoic in the east,the central,the west and the coast of China.These large and medium-scale gas reservoirs are mainly composed of hydrocarbon gas with big dry coefficient,tiny non-hydrocarbon,wide carbon isotope distribution and varying origin types,the hydrocarbon gas includes coal-formed gas,oil-formed gas,biogenic gas and inorganic gas, etc.Coal-formed gas is the main type of China natural gas resources,in particular several explored large-scale gas fields(>100 billion cubic meter)of Kela 2,Sulige and Daniudi,etc.,they all belong to coal-formed gas fields or the gas fields consisting mostly of coal-formed gas.Oil-formed gas is also abundant in China marine basins,for example marine natural gas of Sichuan Basin generated from crude oil cracking gas.Primary and secondary biogenic gas fields were discovered respectively in the Qaidam Basin and Western Slope of Songliao Basin.In addition,inorganic gases are mainly distributed in the eastern China,in particular the Songliao Basin with abundant carbon dioxide accumulation,indicating that the eastern China present large exploration potential of inorganic gas.     

Geological characteristics of large gas provinces and large gas fields in China

Based on the examination of the global researches on oil and gas provinces and large gas fields and the analysis of the features,attributes and distribution of large gas provinces and gas fields,this paper puts forward three indicators of determining large oil and gas provinces:spatial indicator,reservoir-forming indicator and resource indicator.It classifies the gas accumulated areas and large gas provinces in China and analyzes the controlling factors on the distribution of large gas provinces and large gas fields:the lateral distribution is mainly controlled by high-energy sedimentary facies and constructive diagenetic facies,palaeo-highs and their periclinal zones,deep faults,etc,and the vertical distribution is mainly controlled by unconformities,series of evaporates and deep low-velocity highly-conductive beds,etc.It also reveals the main geological characteristics of large gas provinces and large gas fields in China.Large gas fields in four-type basins have their own characteristics and onland large gas fields are dominantly developed in foreland basins and craton basins;there are three types of gas sources,of which,coal is the main source with high gas generating intensity and varying origins;reservoir rocks of the large gas fields(provinces)are of various types and dominated generally by low-middle permeability and porosity pore-type reservoirs;structural traps and litho-stratigraphic traps coexist in Chinese large gas fields and form dense high abundance and large-area low and middle-abundance large gas fields;most of the large gas fields have late hydrocarbon-generation peaks and reservoir formation,and experienced the process of multiple-stage charging and late finalization; large gas provinces(fields)have good sealing and preservation conditions,and evaporates seals are largely developed in large and extra-large gas fields.This paper intends to shed light on the exploration and development of large gas fields(provinces)through analyzing their geological characteristics.     

Variation in noble gas isotopic composition of gas samples from the Aegean arc, Greece

In contrast to most other arcs with oceanic plate subduction, the Aegean arc is characterized by continent–continent subduction. Noble gas abundances and isotopic compositions of 45 gas samples have been determined from 6 volcanoes along the arc, 2 islands in the back-arc region and 7 sites in the surrounding areas. The ³He/⁴He ratios of the samples ranged from 0.027R_A to 6.2R_A (R_A denotes the atmospheric ³He/⁴He ratio of 1.4×10⁻⁶), demonstrating that even the maximum ³He/⁴He ratio in the region is significantly lower than the maximum ratios of most oceanic subduction systems, which are equal to the MORB value of 8±1 R_A. Regional variations in the ³He/⁴He ratio were observed both along and across the arc. The maximum ³He/⁴He ratio was obtained from Nisyros volcano located in the eastern end of the arc, and the ratio decreased westward possibly reflecting the difference in potential degree of crustal assimilation or the present magmatic activity in each volcano. Across the volcanic arc, the ³He/⁴He ratio decreased with an increasing distance from the arc front, reaching a low ratio of 0.063R_A in Macedonia, which suggested a major contribution of radiogenic helium derived from the continental crust. At Nisyros, a temporal increase in ³He/⁴He ratio due to ascending subsurface magma was observed after the seismic crisis of 1995–1998 and mantle neon was possibly detected. The maximum ³He/⁴He ratio (6.2R_A) in the Aegean region, which is significantly lower than the MORB value, is not probably due to crustal assimilation at shallow depth or addition of slab-derived helium to MORB-like mantle wedge, but inherent characteristics of the subcontinental lithospheric mantle (SCLM) beneath the Aegean arc.     

油气勘探中的地震速度泡沫(Fizz)现象

地下岩石孔隙中含有天然气会造成地震波传播速度的下降,但是这种含气量与速度下降量之间的关系一直没有被很清楚地认识到,因此造成很多勘探和开发上的决策失误.本文通过对Wood方程和Gassmann－Biot理论分析,定量地描述了岩石孔隙中天然气含量与波速之间的关系,叙述了由于极少量气造成很大速度变化的“速度泡沫”（Fizz gas）现象以及对这种速度假象进行识别的方法,并通过对中国某油田XD区块的实例研究验证了这种现象的存在及识别.     

Characteristics of shallow gas hydrate in Okhotsk Sea

Multidisciplinary field investigations were carried out in Okhotsk Sea by R/V Akademik M.A. Lavrentyev (LV) of the Russian Academy of Sciences (RAS) in May 2006, supported by funding agencies from Ko- rea, Russia, Japan and China. Geophysical data including echo-sounder, bottom profile, side-scan- sonar, and gravity core sample were obtained aimed to understand the characteristics and formation mechanism of shallow gas hydrates. Based on the geophysical data, we found that the methane flare detected by echo-sounder was the evidence of free gas in the sediment, while the dome structure de- tected by side-scan sonar and bottom profile was the root of gas venting. Gas hydrate retrieved from core on top of the dome structure which was interbedded as thin lamination or lenses with thickness varying from a few millimeters to 3 cm. Gas hydrate content in hydrate-bearing intervals visually amounted to 5%―30% of the sediment volume. This paper argued that gases in the sediment core were not all from gas hydrate decomposition during the gravity core lifting process, free gases must existed in the gas hydrate stability zone, and tectonic structure like dome structure in this paper was free gas central, gas hydrate formed only when gases over-saturated in this gas central, away from these struc- tures, gas hydrate could not form due to low gas concentration.     

基于气象色谱仪产出数据的数据库设计

为了更好的管理BF-2002色谱工作站从气象色谱仪SP3400采集的数据,设计了气象色谱仪数据库,实现了对色谱仪产出的原始数据,谱图数据及气样分析结果等数据的统一管理,方便了观测数据的使用和分析研究。本文就数据库设计及创建过程中所应用的技术进行了简要阐述。     

用于数字化气体观测的脱气-集气装置研制与应用

介绍用于地震地下流体台站数字化气氡、气汞等气体观测中的脱气-集气装置的研究.该装置采用流体溅射的负压脱气原理,可以有效地将自流井地下水中的气体脱析出来,在脱气过程中,可实现对装置的喷射嘴堵塞等清理的功能.     

Natural gas diffusion model and diffusion computation in well Cai25 Bashan Group oil and gas reservoir

Natural gas diffusion through the cap rock is mainly by means of dissolving in water, so its concentration can be replaced by solubility, which varies with temperature, pressure and salinity in strata. Under certain geological conditions the maximal solubility is definite, so the diffusion computation can be handled approximately by stable state equation. Furthermore, on the basis of the restoration of the paleo-buried history, the diffusion is calculated with the dynamic method, and the result is very close to the real diffusion value in the geological history.