Coupling of the Calculated Freezing and Thawing Front Parameterization in the Earth System Model CAS-ESM

The soil freezing and thawing process affects soil physical properties,such as heat conductivity,heat capacity,and hydraulic conductivity in frozen ground regions,and further affects the processes of soil energy,hydrology,and carbon and nitrogen cycles.In this study,the calculation of freezing and thawing front parameterization was implemented into the earth system model of the Chinese Academy of Sciences(CAS-ESM)and its land component,the Common Land Model(CoLM),to investigate the dynamic change of freezing and thawing fronts and their effects.Our results showed that the developed models could reproduce the soil freezing and thawing process and the dynamic change of freezing and thawing fronts.The regionally averaged value of active layer thickness in the permafrost regions was 1.92 m,and the regionally averaged trend value was 0.35 cm yr–1.The regionally averaged value of maximum freezing depth in the seasonally frozen ground regions was 2.15 m,and the regionally averaged trend value was–0.48 cm yr–1.The active layer thickness increased while the maximum freezing depth decreased year by year.These results contribute to a better understanding of the freezing and thawing cycle process.

Wetting front migration model of ion-adsorption rare earth during the multi-hole unsaturated liquid injection

In the process of ion-adsorption rare earth ore leaching,the migration characteristics of the wetting front in multi-hole injection holes and the influence of wetting front intersection effect on the migration distance of wetting fronts are still unclear.Besides,wetting front migration distance and leaching time are usually required to optimize the leaching process.In this study,wetting front migration tests of ionadsorption rare earth ores during the multi-hole fluid injection(the spacing between injection holes was 10 cm,12 cm and 14 cm)and single-hole fluid injection were completed under the constant water head height.At the pre-intersection stage,the wetting front migration laws of ion-adsorption rare earth ores during the multi-hole fluid injection and single-hole fluid injection were identical.At the postintersection stage,the intersection accelerated the wetting front migration.By using the Darcy’s law,the intersection effect of wetting fronts during the multi-hole liquid injection was transformed into the water head height directly above the intersection.Finally,based on the Green-Ampt model,a wetting front migration model of ion-adsorption rare earth ores during the multi-hole unsaturated liquid injection was established.Error analysis results showed that the proposed model can accurately simulate the infiltration process under experimental conditions.The research results enrich the infiltration law and theory of ion-adsorption rare earth ores during the multi-hole liquid injection,and this study provides a scientific basis for optimizing the liquid injection well pattern parameters of ion-adsorption rare earth in situ leaching in the future.

Spatial Inhomogeneity of Atmospheric CO_(2) Concentration and Its Uncertainty in CMIP6 Earth System Models

This paper provides a systematic evaluation of the ability of 12 Earth System Models(ESMs)participating in the Coupled Model Intercomparison Project Phase 6(CMIP6)to simulate the spatial inhomogeneity of the atmospheric carbon dioxide(CO_(2))concentration.The multi-model ensemble mean(MME)can reasonably simulate the increasing trend of CO_(2) concentration from 1850 to 2014,compared with the observation data from the Scripps CO_(2) Program and CMIP6 prescribed data,and improves upon the CMIP5 MME CO_(2) concentration(which is overestimated after 1950).The growth rate of CO_(2) concentration in the northern hemisphere(NH)is higher than that in the southern hemisphere(SH),with the highest growth rate in the mid-latitudes of the NH.The MME can also reasonably simulate the seasonal amplitude of CO_(2) concentration,which is larger in the NH than in the SH and grows in amplitude after the 1950s(especially in the NH).Although the results of the MME are reasonable,there is a large spread among ESMs,and the difference between the ESMs increases with time.The MME results show that regions with relatively large CO_(2) concentrations(such as northern Russia,eastern China,Southeast Asia,the eastern United States,northern South America,and southern Africa)have greater seasonal variability and also exhibit a larger inter-model spread.Compared with CMIP5,the CMIP6 MME simulates an average spatial distribution of CO_(2) concentration that is much closer to the site observations,but the CMIP6-inter-model spread is larger.The inter-model differences of the annual means and seasonal cycles of atmospheric CO_(2) concentration are both attributed to the differences in natural sources and sinks of CO_(2) between the simulations.

Investigating the effect of the Tibetan Plateau on the ITCZ using a coupled Earth system model

本文利用耦合地球气候系统模式研究了青藏高原对热带辐合带(ITCZ)的影响.我们研究发现热带大西洋ITCZ的位置对青藏高原存在与否有明显的敏感性.与目前真实情况相比,移除青藏高原会导致北半球海面降温,南半球海面升温.这种海面温度变化在大西洋表现得尤为明显,导致热带大西洋最大海温中心向南移动,从而迫使大气对流中心向南移动,即表现为ITCZ的南移.相应地,夏季热带大气Hadley环流的上升支也发生明显南移.北(南)半球海洋变冷(变暖)这种态势要求增强跨赤道向北的大气经向热量输送,从而维持各个半球的能量平衡,而这需要ITCZ位置的南移才能实现.本文研究表明,青藏高原的存在在现今ITCZ气候态的形成中可能扮演了重要角色.

Major Modes of Short-Term Climate Variability in the Newly Developed NUIST Earth System Model(NESM)

A coupled earth system model（ESM） has been developed at the Nanjing University of Information Science and Technology（NUIST） by using version 5.3 of the European Centre Hamburg Model（ECHAM）, version 3.4 of the Nucleus for European Modelling of the Ocean（NEMO）, and version 4.1 of the Los Alamos sea ice model（CICE）. The model is referred to as NUIST ESM1（NESM1）. Comprehensive and quantitative metrics are used to assess the model＇s major modes of climate variability most relevant to subseasonal-to-interannual climate prediction. The model＇s assessment is placed in a multi-model framework. The model yields a realistic annual mean and annual cycle of equatorial SST, and a reasonably realistic precipitation climatology, but has difficulty in capturing the spring–fall asymmetry and monsoon precipitation domains. The ENSO mode is reproduced well with respect to its spatial structure, power spectrum, phase locking to the annual cycle, and spatial structures of the central Pacific（CP）-ENSO and eastern Pacific（EP）-ENSO; however, the equatorial SST variability,biennial component of ENSO, and the amplitude of CP-ENSO are overestimated. The model captures realistic intraseasonal variability patterns, the vertical-zonal structures of the first two leading predictable modes of Madden–Julian Oscillation（MJO）, and its eastward propagation; but the simulated MJO speed is significantly slower than observed. Compared with the T42 version, the high resolution version（T159） demonstrates improved simulation with respect to the climatology, interannual variance, monsoon–ENSO lead–lag correlation, spatial structures of the leading mode of the Asian–Australian monsoon rainfall variability, and the eastward propagation of the MJO.

Global Climate Internal Variability in a 2000-year Control Simulation with Community Earth System Model(CESM)

Using the low-resolution(T31, equivalent to 3.75°× 3.75°) version of the Community Earth System Model(CESM) from the National Center for Atmospheric Research(NCAR), a global climate simulation was carried out with fixed external forcing factors(1850 Common Era.(C.E.) conditions) for the past 2000 years. Based on the simulated results, spatio-temporal structures of surface air temperature, precipitation and internal variability, such as the El Nio-Southern Oscillation(ENSO), the Atlantic Multi-decadal Oscillation(AMO), the Pacific Decadal Oscillation(PDO), and the North Atlantic Oscillation(NAO), were compared with reanalysis datasets to evaluate the model performance. The results are as follows: 1) CESM showed a good performance in the long-term simulation and no significant climate drift over the past 2000 years; 2) climatological patterns of global and regional climate changes simulated by the CESM were reasonable compared with the reanalysis datasets; and 3) the CESM simulated internal natural variability of the climate system performs very well. The model not only reproduced the periodicity of ENSO, AMO and PDO events but also the 3–8 years variability of the ENSO. The spatial distribution of the CESM-simulated NAO was also similar to the observed. However, because of weaker total irradiation and greenhouse gas concentration forcing in the simulation than the present, the model performances had some differences from the observations. Generally, the CESM showed a good performance in simulating the global climate and internal natural variability of the climate system. This paves the way for other forced climate simulations for the past 2000 years by using the CESM.

Global air-sea CO_(2) exchange fl ux since 1980s: results from CMIP6 Earth System Models

The ocean could profoundly modulate the ever-increasing atmospheric CO_(2) by air-sea CO_(2) exchange process,which is also able to cause signifi cant changes of physical and biogeochemical properties in return.In this study,we assessed the long-term average and spatial-temporal variability of global air-sea CO_(2) exchange fl ux(F CO_(2))since 1980s basing on the results of 18 Coupled Model Intercomparison Project Phase 6(CMIP6)Earth System Models(ESMs).Our fi ndings indicate that the CMIP6 ESMs simulated global CO_(2) sink in recent three decades ranges from 1.80 to 2.24 Pg C/a,which is coincidence with the results of cotemporaneous observations.What’s more,the CMIP6 ESMs consistently show that the global oceanic CO_(2) sink has gradually intensifi ed since 1980s as well as the observations.This study confi rms the simulated F CO_(2) could reach agreements with the observations in the aspect of primary climatological characteristics,however,the simulation skills of CIMP6 ESMs in diverse open-sea biomes are unevenness.None of the 18 CMIP6 ESMs could reproduce the observed F CO_(2) increasement in the central-eastern tropical Pacifi c and the midlatitude Southern Ocean.Defi ciencies of some CMIP6 ESMs in reproducing the atmospheric pressure systems of the Southern Hemisphere and the El Niño-Southern Oscillation(ENSO)mode of the tropical Pacifi c are probably the major causes.

Earth System Model FGOALS-s2: Coupling a Dynamic Global Vegetation and Terrestrial Carbon Model with the Physical Climate System Model

Earth System Models （ESMs） are fundamental tools for understanding climate-carbon feedback. An ESM version of the Flexible Global Ocean-Atmosphere-Land System model （FGOALS） was recently developed within the IPCC AR5 Coupled Model Intercomparison Project Phase 5 （CMIP5） modeling framework, and we describe the development of this model through the coupling of a dynamic global vegetation and terrestrial carbon model with FGOALS-s2. The performance of the coupled model is evaluated as follows. The simulated global total terrestrial gross primary production （GPP） is 124.4 PgC yr-I and net pri- mary production （NPP） is 50.9 PgC yr-1. The entire terrestrial carbon pools contain about 2009.9 PgC, comprising 628.2 PgC and 1381.6 PgC in vegetation and soil pools, respectively. Spatially, in the tropics, the seasonal cycle of NPP and net ecosystem production （NEP） exhibits a dipole mode across the equator due to migration of the monsoon rainbelt, while the seasonal cycle is not so significant in Leaf Area Index （LAI）. In the subtropics, especially in the East Asian monsoon region, the seasonal cycle is obvious due to changes in temperature and precipitation from boreal winter to summer. Vegetation productivity in the northern mid-high latitudes is too low, possibly due to low soil moisture there. On the interannual timescale, the terrestrial ecosystem shows a strong response to ENSO. The model- simulated Nifio3.4 index and total terrestrial NEP are both characterized by a broad spectral peak in the range of 2-7 years. Further analysis indicates their correlation coefficient reaches -0.7 when NEP lags the Nifio3.4 index for about 1-2 months.

Lateral earth pressure of granular backfills on retaining walls with expanded polystyrene geofoam inclusions under limited surcharge loading

Existing studies have focused on the behavior of the retaining wall equipped with expanded polystyrene(EPS)geofoam inclusions under semi-infinite surcharge loading rather than limited surcharge loading.In this paper,the failure mode and the earth pressure acting on the rigid retaining wall with EPS geofoam inclusions and granular backfills(henceforth referred to as EPS-wall),under limited surcharge loading are investigated through two-and three-dimensional model tests.The testing results show that different from the sliding of almost all the backfill in the EPS-wall under semi-infinite surcharge loading,only an approximately triangular backfill slides in the wall under limited surcharge loading.The distribution of the lateral earth pressure on the EPS-wall under limited surcharge loading is non-linear,and the distribution changes from the increase of the wall depth to the decrease with the increase of the limited surcharge loading.An approach based on the force equilibrium of a differential element is developed to predict the lateral earth pressure behind the EPS-wall subjected to limited surcharge loading,and its performance was fully validated by the three-dimensional model tests.

Regional earth system modeling：review and future directions

耦合区域地球系统模式(RESM)现在仍处于早期的发展阶段,在本文中,我们对其近期相关进展进行了回顾。到目前为止,已经开发出了耦合的区域大气–海洋–海冰、大气–气溶胶和大气–生物圈模式,但总体它们仅在有限的区域得到应用,需要更多的工作来评估其在更多区域的可移植性。我们认为RESM发展中的未来挑战,是在大气、海洋、冰冻圈、生物圈、化学圈以外,同时将人类及其活动成分以完全相互作用的方式引入进来。

Calculation of Interaction Parameters from Immiscible Phase Diagram of Alkali Metal or Alkali Earth Metal-Halide System by Means of Subregular Solution Model

In this paper, the interaction parameters in the subregular solution model, λ1 and λ2, are regarded as a linear function of temperature, T. Therefore, the molar excess Gibbs energy of A-B binary system may be reexpressed as follows:Gm^E=xAxB[(λ11+λ12T)+(λ21+λ22T)xB]The calculation of the model parameters, λ11, λ12, λ21and λ22, was carried out numerically from the phase diagrams for 11 alkali metal-alkali halide or alkali earth metal-halide systems. In addition, artificial neural network trained by known data has been used to predict the values of these model parameters. The predicted results are in good agreement with the .calculated ones. The applicability of the subregular solution model to the alkali metal-alkali halide or alkali earth metal-halide systems were tested by comparing the available experimental composition along the boundary of miscibility gap with the calculated ones which were obtained by using genetic algorithm. The good agreement between the calculated and experimental results across the entire liquidus is valid evidence in support of the model.

A modified stochastic model for LS+AR hybrid method and its application in polar motion short-term prediction

Short-term(up to 30 days)predictions of Earth Rotation Parameters(ERPs)such as Polar Motion(PM:PMX and PMY)play an essential role in real-time applications related to high-precision reference frame conversion.Currently,least squares(LS)+auto-regressive(AR)hybrid method is one of the main techniques of PM prediction.Besides,the weighted LS+AR hybrid method performs well for PM short-term prediction.However,the corresponding covariance information of LS fitting residuals deserves further exploration in the AR model.In this study,we have derived a modified stochastic model for the LS+AR hybrid method,namely the weighted LS+weighted AR hybrid method.By using the PM data products of IERS EOP 14 C04,the numerical results indicate that for PM short-term forecasting,the proposed weighted LS+weighted AR hybrid method shows an advantage over both the LS+AR hybrid method and the weighted LS+AR hybrid method.Compared to the mean absolute errors(MAEs)of PMX/PMY sho rt-term prediction of the LS+AR hybrid method and the weighted LS+AR hybrid method,the weighted LS+weighted AR hybrid method shows average improvements of 6.61%/12.08%and 0.24%/11.65%,respectively.Besides,for the slopes of the linear regression lines fitted to the errors of each method,the growth of the prediction error of the proposed method is slower than that of the other two methods.

Simulations of dissolved oxygen concentration in CMIP5 Earth system models

The climatologies of dissolved oxygen concentration in the ocean simulated by nine Earth system models(ESMs) from the historical emission driven experiment of CMIP5(Phase 5 of the Climate Model Intercomparison Project) are quantitatively evaluated by comparing the simulated oxygen to the WOA09 observation based on common statistical metrics. At the sea surface, distribution of dissolved oxygen is well simulated by all nine ESMs due to well-simulated sea surface temperature(SST), with both globally-averaged error and root mean square error(RMSE) close to zero, and both correlation coefficients and normalized standard deviation close to 1. However, the model performance differs from each other at the intermediate depth and deep ocean where important water masses exist. At the depth of 500 to 1 000 m where the oxygen minimum zones(OMZs) exist, all ESMs show a maximum of globally-averaged error and RMSE, and a minimum of the spatial correlation coefficient. In the ocean interior, the reason for model biases is complicated, and both the meridional overturning circulation(MOC) and the particulate organic carbon flux contribute to the biases of dissolved oxygen distribution. Analysis results show the physical bias contributes more. Simulation bias of important water masses such as North Atlantic Deep Water(NADW), Antarctic Bottom Water(AABW) and North Pacific Intermediate Water(NPIW) indicated by distributions of MOCs greatly affects the distributions of oxygen in north Atlantic, Southern Ocean and north Pacific, respectively.Although the model simulations of oxygen differ greatly from each other in the ocean interior, the multi-model mean shows a better agreement with the observation.

Optimisation of Thermal Comfort of Building in a Hot and Dry Tropical Climate: A Comparative Approach between Compressed Earth/Concrete Block Envelopes

Compressed earth blocks (CEB) are an alternative to cement blocks in the construction of wall masonry. However, the optimal architectural construction methods for adequate thermal comfort for occupants in hot and arid environments are not mastered. This article evaluates the influence of architectural and constructive modes of buildings made of CEB walls and concrete block walls, to optimize and compare their thermal comfort in the hot and dry tropical climate of Ouagadougou, Burkina Faso. Two identical pilot buildings whose envelopes are made of CEB and concrete blocks were monitored for this study. The thermal models of the pilot buildings were implemented in the SketchUp software using an extension of EnergyPlus. The models were empirically validated after calibration against measured thermal data from the buildings. The models were used to do a parametric analysis for optimization of the thermal performances by simulating plaster coatings on the exterior of walls, airtight openings and natural ventilation depending on external weather conditions. The results show that the CEB building displays 7016 hours of discomfort, equivalent to 80.1% of the time, and the concrete building displays 6948 hours of discomfort, equivalent to 79.3% of the time. The optimization by modifications reduced the discomfort to 2918 and 3125 hours respectively;i.e. equivalent to only 33.3% for the CEB building and 35.7% for the concrete building. More study should evaluate thermal optimizations in buildings in real time of usage such as residential buildings commonly used by the local middle class. The use of CEB as a construction material and passive means of improving thermal comfort is a suitable ecological and economical option to replace cementitious material.

Dark Galaxies, Sun-Earth-Moon Interaction, Tunguska Event—Explained by WUM

Great experimental results and observations achieved by Astronomy in the last decades revealed new unexplainable phenomena. Astronomers have conclusive new evidence that a recently discovered “dark galaxy” is, in fact, an object the size of a galaxy, made entirely of dark matter. They found that the speed of the Earth’s rotation varies randomly each day. 115 years ago, the Tunguska Event was observed, and astronomers still do not have an explanation of It. Main results of the present article are: 1) Dark galaxies explained by the spinning of their Dark Matter Cores with the surface speed at equator less than the escape velocity. Their Rotational Fission is not happening. Extrasolar systems do not emerge;2) 21-cm Emission explained by the self-annihilation of Dark Matter particles XIONs (5.3 μeV);3) Sun-Earth-Moon Interaction explained by the influence of the Sun’s and the Moon’s magnetic field on the electrical currents of the charged Geomagma (the 660-km layer), and, as a result, the Earth’s daylength varies;4) Tunguska Event explained by a huge atmospheric explosion of the Superbolide, which was a stable Dark Matter Bubble before entering the Earth’s atmosphere.

气候模式CAS-ESM-C对1月热带太平洋流场模态的模拟评估

利用气候模式CAS-ESM-C从1922年起84年的模拟资料,对1月份热带太平洋上层流场作复EOF分解及小波分析,并与实况以及理论解析解作对比讨论,以考察模式对赤道大洋上层流场的模拟能力,得到主要结论:(1)复EOF分解前3个模态的方差贡献为53.5%、12.9%、9.5%,累积方差贡献为75.9%,累积方差贡献比实况更高。(2)第一、二模态空间场与实况相比总体相像,流场都为赤道所俘获,在俘获区内的流场均以偏纬向流为主;差异在于模拟资料分析的赤道俘获区范围较实况要大,流场的经向流分量及越赤道流也较实况明显。(3)第一、二模态实时间系数序列无线性变化趋势,而实况则有。复EOF模态年际及年代际变化与实况相同或相近;第一、二模态中3~7年的年际变化是厄尔尼诺与南方涛动(ENSO)的反映;第一模态22~23年的年代际变化受北太平洋主要气候模态北太平洋年代际振荡(PDO)对热带太平洋的影响,而第二模态13年的年代际变化是受北太平洋次要气候模态北太平洋环流振荡(NPGO)对热带太平洋的影响;第一、二模态还都有峰值16年的年代际变化,这可能与印尼穿越流有关。(4)模拟资料分析的结果具有理论解析解中流场为赤道所俘获及流场为纬向流的特点,只是解析解中因无风应力强迫,流场呈纯纬向流。(5)第一(二)模态在赤道太平洋东部(中部)有海温动力异常,并可称之为东(中)部型ENSO模态。气候模式CAS-ESM-C对热带太平洋上层流场的模拟表现较佳。

A Brief Introduction to BNU-HESM1.0 and Its Earth Surface Temperature Simulations

Integrated assessment models and coupled earth system models both have their limitations in understanding the interactions between human activity and the physical earth system. In this paper, a new human--earth system model, BNU- HESM1.0, constructed by combining the economic and climate damage components of the Dynamic Integrated Model of Climate Change and Economy to the BNU-ESM model, is introduced. The ability of BNU-HESM1.0 in simulating the global CO2 concentration and surface temperature is also evaluated. We find that, compared to observation, BNU-HESM1.0 underestimates the global CO2 concentration and its rising trend during 1965-2005, due to the uncertainty in the economic components. However, the surface temperature simulated by BNU-HESM1.0 is much closer to observation, resulting from the overestimates of surface temperature by the original BNU-ESM model. The uncertainty of BNU-ESM falls within the range of present earth system uncertainty, so it is the economic and climate damage component of BNU-HESM 1.0 that needs to be improved through further study. However, the main purpose of this paper is to introduce a new approach to investigate the complex relationship between human activity and the earth system. It is hoped that it will inspire further ideas that prove valuable in guiding human activities appropriate for a sustainable future climate.

Seepage safety monitoring model for an earth rock dam under influence of high-impact typhoons based on particle swarm optimization algorithm

Extreme hydrological events induced by typhoons in reservoir areas have presented severe challenges to the safe operation of hydraulic structures. Based on analysis of the seepage characteristics of an earth rock dam, a novel seepage safety monitoring model was constructed in this study. The nonlinear influence processes of the antecedent reservoir water level and rainfall were assumed to follow normal distributions. The particle swarm optimization （PSO） algorithm was used to optimize the model parameters so as to raise the fitting accuracy. In addition, a mutation factor was introduced to simulate the sudden increase in the piezometric level induced by short-duration heavy rainfall and the possible historical extreme reservoir water level during a typhoon. In order to verify the efficacy of this model, the earth rock dam of the Siminghu Reservoir was used as an example. The piezometric level at the SW1-2 measuring point during Typhoon Fitow in 2013 was fitted with the present model, and a corresponding theoretical expression was established. Comparison of fitting results of the piezometric level obtained from the present statistical model and traditional statistical model with monitored values during the typhoon shows that the present model has a higher fitting accuracy and can simulate the uprush feature of the seepage pressure during the typhoon perfectly.

Component Content Soft-Sensor Based on Hybrid Models in Countercurrent Rare Earth Extraction Process

In consideration of the online measurement of the component content in rare earth countercurrent extraction separation process, the soft sensor method based on hybrid modeling was proposed to measure the rare earth component content. The hybrid models were composed of the extraction equilibrium calculation model and the Radial Basis Function (RBF) Neural Network (NN) error compensation model; the parameters of compensation model were optimized by the hierarchical genetic algorithms (HGA). In addition, application experiment research of this proposed method was carried out in the rare earth separation production process of a corporation. The result shows that this method is effective and can realize online measurement for the component content of rare earth in the countercurrent extraction.

Evaluation of behaviors of earth and rockfill dams during construction and initial impounding using instrumentation data and numerical modeling

In this study,the behavior of Gavoshan dam was evaluated during construction and the first impounding.A two-dimensional（2D） numerical analysis was conducted based on a finite difference method on the largest cross-section of the dam using the results of instrument measurements and back analysis.These evaluations will be completed in the case that back analysis is carried out in order to control the degree of the accuracy and the level of confidence of the measured behavior since each of the measurements could be controlled by comparing it to the result obtained from the numerical model.Following that,by comparing the results of the numerical analysis with the measured values,it is indicated that there is a proper consistency between these two values.Moreover,it was observed that the dam performance was suitable regarding the induced pore water pressure,the pore water pressure ratio r;,settlement,induced stresses,arching degree,and hydraulic fracturing probability during the construction and initial impounding periods.The results demonstrated that the maximum settlement of the core was 238 cm at the end of construction.In the following 6 years after construction（initial impounding and exploitation period）,the accumulative settlement of the dam was 270 cm.It is clear that 88% of the total settlement of the dam took place during dam construction.The reason is that the clay core was smashed in the wet side,i.e.the optimum moisture content.Whereas the average curving ratio was 0.64 during dam construction; at the end of the initial impounding,the maximum amount of curving ratio in the upstream was 0.81,and the minimum（critical） amount in the downstream was 0.52.It was also concluded that this dam is safe in comparison with the behaviors of other similar dams in the world.