To Diversify or to Assimilate—Speculating on Multicultural Education

Multicultural education(MCED) is an issue not appropriately treated in China from the perspective of its adaptiveness and legitimacy.This paper tries to clarify the concept in a light of three different scholastic schools.To achieve diversity or assimilation as the ultimate purpose of MCED is also tentatively discussed.

Winter wheat yield estimation based on assimilated Sentinel-2 images with the CERES-Wheat model

Assimilating Sentinel-2 images with the CERES-Wheat model can improve the precision of winter wheat yield estimates at a regional scale. To verify this method, we applied the ensemble Kalman filter(EnKF) to assimilate the leaf area index(LAI) derived from Sentinel-2 data and simulated by the CERES-Wheat model. From this, we obtained the assimilated daily LAI during the growth stage of winter wheat across three counties located in the southeast of the Loess Plateau in China: Xiangfen, Xinjiang, and Wenxi. We assigned LAI weights at different growth stages by comparing the improved analytic hierarchy method, the entropy method, and the normalized combination weighting method, and constructed a yield estimation model with the measurements to accurately estimate the yield of winter wheat. We found that the changes of assimilated LAI during the growth stage of winter wheat strongly agreed with the simulated LAI. With the correction of the derived LAI from the Sentinel-2 images, the LAI from the green-up stage to the heading–filling stage was enhanced, while the LAI decrease from the milking stage was slowed down, which was more in line with the actual changes of LAI for winter wheat. We also compared the simulated and derived LAI and found the assimilated LAI had reduced the root mean square error(RMSE) by 0.43 and 0.29 m^(2) m^(–2), respectively, based on the measured LAI. The assimilation improved the estimation accuracy of the LAI time series. The highest determination coefficient(R2) was 0.8627 and the lowest RMSE was 472.92 kg ha^(–1) in the regression of the yields estimated by the normalized weighted assimilated LAI method and measurements. The relative error of the estimated yield of winter wheat in the study counties was less than 1%, suggesting that Sentinel-2 data with high spatial-temporal resolution can be assimilated with the CERES-Wheat model to obtain more accurate regional yield estimates.

Numerical simulation of scatterometer assimilated wind and ocean wave in eastern China seas and adjacent waters

Using the latest version of Mesoscale Modeling System （MMSv3）, we assimilated wind data from the scatterometer and built a model to assimilate the wind field over eastern China seas and adjacent waters and applied the wave model WAVEWATCH-Ⅲ to test the sea area with assimilative wind and blended wind of QSCAT and NCEP as driving forces. High precision and resolution numerical wave results were obtained. Analysis indicated that if we replace the model wind result with the blended wind, better sea surface wind results and wave results could be obtained.

THE APPLICATION OF ASSIMILATED AIRCRAFT DATA IN SIMULATING A HEAVY RAIN OVER SOUTH CHINA IN JUNE 2005

Regular and irregular observational data are used to analyze and simulate a torrential rain over the south of China on 18 - 24 June 2005. Since the regular data cannot depict the rainfall system fully, GRAPES model is used to simulate this process. Different data are assimilated for 12 hours by its simulating system and different analysis data are obtained. In order to analyze how well the model forecast has been improved with the addition of assimilated aircraft data, these different analysis data are used as the first-guess data to conduct two control numerical simulation tests. From these tests, it is proved that be model that adds aircraft assimilation data can simulate the main region of precipitation, which is more consistent with the observed precipitation than the model that does not, and that the accuracy rate is also improved. These numerical simulation tests not only show that it is necessary and capable to improve the modeling of this torrential rain process by using aircraft data, but also lays the foundation for forecasting heavy rains in the south of China based on aircraft data.

Simulated Indonesian Throughflow in Makassar Strait across the SODA3 products

The Indonesian Throughflow(ITF), which connects the tropical Pacific and Indian oceans, plays important roles in the inter-ocean water exchange and regional or even global climate variability. The Makassar Strait is the main inflow passage of the ITF, carrying about 77% of the total ITF volume transport. In this study, we analyze the simulated ITF in the Makassar Strait in the Simple Ocean Data Assimilation version 3(SODA3) datasets. A total of nine ensemble members of the SODA3 datasets, of which are driven by different surface forcings and bulk formulas, and with or without data assimilation, are used in this study. The annual mean water transports(i.e.,volume, heat and freshwater) are related to the combination of surface forcing and bulk formula, as well as whether data assimilation is employed. The phases of the seasonal and interannual variability in water transports cross the Makassar Strait, are basically consistent with each other among the SODA3 ensemble members. The interannual variability in Makassar Strait volume and heat transports are significantly correlated with El Ni?oSouthern Oscillation(ENSO) at time lags of-6 to 7 months. There is no statistically significant correlation between the freshwater transport and the ENSO. The Makassar Strait water transports are not significantly correlated with the Indian Ocean Dipole(IOD), which may attribute to model deficiency in simulating the propagation of semiannual Kelvin waves from the Indian Ocean to the Makassar Strait.

Construction and Application of a Regional Kilometer-Scale Carbon Source and Sink Assimilation Inversion System(CCMVS-R)

CO_(2)is one of the most important greenhouse gases(GHGs)in the earth’s atmosphere.Since the industrial era,anthropogenic activities have emitted excessive quantities of GHGs into the atmosphere,resulting in climate warming since the 1950s and leading to an increased frequency of extreme weather and climate events.In 2020,China committed to striving for carbon neutrality by 2060.This commitment and China’s consequent actions will result in significant changes in global and regional anthropogenic carbon emissions and therefore require timely,comprehensive,and objective monitoring and verification support(MVS)systems.The MVS approach relies on the top-down assimilation and inversion of atmospheric CO_(2)concentrations,as recommended by the Intergovernmental Panel on Climate Change(IPCC)Inventory Guidelines in 2019.However,the regional high-resolution assimilation and inversion method is still in its initial stage of development.Here,we have constructed an inverse system for carbon sources and sinks at the kilometer level by coupling proper orthogonal decomposition(POD)with four-dimensional variational(4DVar)data assimilation based on the weather research and forecasting-greenhouse gas(WRF-GHG)model.Our China Carbon Monito ring and Verification Support at the Regional level(CCMVS-R)system can continuously assimilate information on atmospheric CO_(2)and other related information and realize the inversion of regional and local anthropogenic carbon emissions and natural terrestrial ecosystem carbon exchange.Atmospheric CO_(2)data were collected from six ground-based monito ring sites in Shanxi Province,China to verify the inversion effect of regio nal anthropogenic carbon emissions by setting ideal and real experiments using a two-layer nesting method(at 27 and 9 km).The uncertainty of the simulated atmospheric CO_(2)decreased significantly,with a root-mean-square error of CO_(2)concentration values between the ideal value and the simulated after assimilation was close to 0.The total anthropogenic carbon emissions in Shanxi Province in 2019 from the assimilated inversions were approximately 28.6%(17%-38%)higher than the mean of five emission inventories using the bottomup method,showing that the top-down CCMVS-R system can obtain more comprehensive information on anthropogenic carbon emissions.

Shoot Apex Demand Determines Assimilate and Nutrients Partitioning and Nutrient-uptake Rate in Tobacco Plants

Our previous experiment revealed that apex-removed plants have larger root systems but a lower K＋-uptake rates than intact tobacco plants. Since the apex is not only a center of growth and metabolism, but also an important place of auxin synthesis and export, the aims of this study were to distinguish whether the apex demand or auxin synthesized in the apex regulates assimilate and nutrients partitioning within plant, and to explain the reason for the lower K＋-uptake rate of the apex-removed plant. In comparison with the control plant, covering the shoot apex with a black transparent plastic bag reduced net increases in dry matter and nutrients; however, the distribution of the dry matter and nutrients between shoot and roots and nutrient-uptake rates were not changed. Removal of the shoot apex shifted the dry mass and nutrients distributions to roots, and reduced the rate of nutrient uptake. Application of 1-naphthylacetic acid （NAA） could partly replace the role of the removed apex, stimulated assimilate and nutrient deposition into the treated tissue, and enhanced the reduced plasma membrane ATPase activity of roots to the control level. However, treatment of the apex-removed plants with NAA could not rescue the reduced nutrient uptake rate and the shifted assimilates and nutrients partitioning caused by excision of the apex. Higher nutrient uptake rate of the intact plants could not be explained by root growth parameters, such as total root surface area and number of root tips. The results from the present study indicate that strong apex demand determined assimilates and nutrients partitioning and nutrient-uptake rate in tobacco （Nicotiana tabacum） plants.

Assimilation of GOES-R Geostationary Lightning Mapper Flash Extent Density Data in GSI 3DVar, EnKF, and Hybrid En3DVar for the Analysis and Short-Term Forecast of a Supercell Storm Case

Capabilities to assimilate Geostationary Operational Environmental Satellite “R-series ”(GOES-R) Geostationary Lightning Mapper(GLM) flash extent density(FED) data within the operational Gridpoint Statistical Interpolation ensemble Kalman filter(GSI-EnKF) framework were previously developed and tested with a mesoscale convective system(MCS) case. In this study, such capabilities are further developed to assimilate GOES GLM FED data within the GSI ensemble-variational(EnVar) hybrid data assimilation(DA) framework. The results of assimilating the GLM FED data using 3DVar, and pure En3DVar(PEn3DVar, using 100% ensemble covariance and no static covariance) are compared with those of EnKF/DfEnKF for a supercell storm case. The focus of this study is to validate the correctness and evaluate the performance of the new implementation rather than comparing the performance of FED DA among different DA schemes. Only the results of 3DVar and pEn3DVar are examined and compared with EnKF/DfEnKF. Assimilation of a single FED observation shows that the magnitude and horizontal extent of the analysis increments from PEn3DVar are generally larger than from EnKF, which is mainly caused by using different localization strategies in EnFK/DfEnKF and PEn3DVar as well as the integration limits of the graupel mass in the observation operator. Overall, the forecast performance of PEn3DVar is comparable to EnKF/DfEnKF, suggesting correct implementation.

Comparison of Adaptive Simulation Observation Experiments of the Heavy Rainfall in South China and Sichuan Basin

This study examines the effectiveness of adaptive observation experiments using the ensemble transformation sensitivity(ETS) method to improve precipitation forecasts during heavy rainfall events in South China and the Sichuan Basin. High-resolution numerical models are employed to simulate adaptive observations. By identifying the sensitive areas of key weather system positions 42 hours before heavy rainfall events, the adaptive observations improve the prediction of jet streams, strong winds, and shear lines, which are essential for accurate heavy rainfall forecasting. This improvement is reflected in both the precipitation structure and location accuracy within the verification region. In South China, targeted observations enhance rainfall predictions by improving water vapor transport. In the Sichuan Basin, adaptive observations refine water vapor transport and adjust vortex dynamics. This research highlights the importance of accurately predicting shear lines and jet streams for forecasting heavy rainfall in these areas. Overall, this study found that adaptive observation enhances the precipitation forecast skills of the structure and location for heavy rainfall in South China and the Sichuan Basin, emphasizing their potential utility in operational numerical weather prediction.

The effect of glutathione on glucosinolate biosynthesis through the sulfur assimilation pathway in pakchoi associated with the growth conditions

Glucosinolates(GSLs) are a group of nitrogen-and sulfur-containing secondary metabolites, synthesized primarily in members of the Brassicaceae family, that play an important role in food flavor, plant antimicrobial activity, resistance to insect attack, stress tolerance, and human anti-cancer effects. As a sulfur-containing compound, glutathione has a strong connection with GSLs biosynthesis as a sulfur donor or redox system, and exists in reduced(glutathione;GSH) and oxidized(glutathione disulfide;GSSG) forms. However, the mechanism of GSH regulating GSLs biosynthesis remainds unclear. Hence, the exogenous therapy to pakchoi under normal growth condition and sulfur deficiency condition were conducted in this work to explore the relevant mechanism. The results showed that exogenous application of buthionine sulfoximine, an inhibitor of GSH synthesis, decreased the transcript levels of GSLs synthesis-related genes and transcription factors, as well as sulfur assimilation-related genes under the normal growth condition. Application of exogenous GSH inhibited the expression of GSLs synthesis-and sulfur assimilation-related genes under the normal condition, while the GSLs biosynthesis and the sulfur assimilation pathway were activated by exogenous application of GSH when the content of GSH in vivo of plants decreased owing to sulfur deficiency. Moreover,exogenous application of GSSG increased the transcript levels of GSLs synthesis-and sulfur assimilation-related genes under the normal growth condition and under sulfur deficiency. The present work provides new insights into the molecular mechanisms of GSLs biosynthesis underlying glutathione regulation.

Application of the finite analytic numerical method to a flowdependent variational data assimilation

An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection term,the discrete method needs to be chosen very carefully.The finite analytic method is an alternative scheme to solve the advection-diffusion equation.As a combination of analytical and numerical methods,it not only has high calculation accuracy but also holds the characteristic of the auto upwind.To demonstrate its ability,the one-dimensional steady and unsteady advection-diffusion equation numerical examples are respectively solved by the finite analytic method.The more widely used upwind difference method is used as a control approach.The result indicates that the finite analytic method has higher accuracy than the upwind difference method.For the two-dimensional case,the finite analytic method still has a better performance.In the three-dimensional variational assimilation experiment,the finite analytic method can effectively improve analysis field accuracy,and its effect is significantly better than the upwind difference and the central difference method.Moreover,it is still a more effective solution method in the strong flow region where the advective-diffusion filter performs most prominently.

Optimal Assimilation of Microwave Upper-Level Sounding Data in CMA-GFS

Various approaches have been proposed to minimize the upper-level systematic biases in global numerical weather prediction(NWP)models by using satellite upper-air sounding channels as anchors.However,since the China Meteorological Administration Global Forecast System(CMA-GFS)has a model top near 0.1 hPa(60 km),the upper-level temperature bias may exceed 4 K near 1 hPa and further extend to 5 hPa.In this study,channels 12–14 of the Advanced Microwave Sounding Unit A(AMSU-A)onboard five satellites of NOAA and METOP,whose weighting function peaks range from 10 to 2 hPa are all used as anchor observations in CMA-GFS.It is shown that the new“Anchor”approach can effectively reduce the biases near the model top and their downward propagation in three-month assimilation cycles.The bias growth rate of simulated upper-level channel observations is reduced to±0.001 K d^(–1),compared to–0.03 K d^(–1)derived from the current dynamic correction scheme.The relatively stable bias significantly improves the upper-level analysis field and leads to better global medium-range forecasts up to 10 days with significant reductions in the temperature and geopotential forecast error above 10 hPa.

Assimilating FY-4A AGRI Radiances with a Channel-Sensitive Cloud Detection Scheme for the Analysis and Forecasting of Multiple Typhoons

This paper presents an attempt at assimilating clear-sky FY-4A Advanced Geosynchronous Radiation Imager(AGRI)radiances from two water vapor channels for the prediction of three landfalling typhoon events over the West Pacific Ocean using the 3DVar data assimilation(DA)method along with the WRF model.A channel-sensitive cloud detection scheme based on the particle filter(PF)algorithm is developed and examined against a cloud detection scheme using the multivariate and minimum residual(MMR)algorithm and another traditional cloud mask–dependent cloud detection scheme.Results show that both channel-sensitive cloud detection schemes are effective,while the PF scheme is able to reserve more pixels than the MMR scheme for the same channel.In general,the added value of AGRI radiances is confirmed when comparing with the control experiment without AGRI radiances.Moreover,it is found that the analysis fields of the PF experiment are mostly improved in terms of better depicting the typhoon,including the temperature,moisture,and dynamical conditions.The typhoon track forecast skill is improved with AGRI radiance DA,which could be explained by better simulating the upper trough.The impact of assimilating AGRI radiances on typhoon intensity forecasts is small.On the other hand,improved rainfall forecasts from AGRI DA experiments are found along with reduced errors for both the thermodynamic and moisture fields,albeit the improvements are limited.

Added Benefit of the Early-Morning-Orbit Satellite Fengyun-3E on the Global Microwave Sounding of the Three-Orbit Constellation

The three-orbit constellation can comprehensively increase the spatial coverage of polar-orbiting satellites,but the polar-orbiting satellites currently in operation are only mid-morning-orbit and afternoon-orbit satellites.Fengyun-3E(FY-3E)was launched successfully on 5 July 2021 in China.As an early-morning-orbit satellite,FY-3E can help form a complete three-orbit observation system together with the mid-morning and afternoon satellites in the current mainstream operational system.In this study,we investigate the added benefit of FY-3E microwave sounding observations to the midmorning-orbit Meteorological Operational satellite-B(Met Op-B)and afternoon-orbit Fengyun-3D(FY-3D)microwave observations in the Chinese Meteorological Administration global forecast system(CMA-GFS).The results show that the additional FY-3E microwave temperature sounder-3(MWTS-3)and microwave humidity sounder-2(MWHS-2)data can increase the global coverage of microwave temperature and humidity sounding data by 14.8% and 10.6%,respectively.It enables the CMA-GFS to achieve nearly 100% global coverage of microwave-sounding observations at each analysis time.Furthermore,after effective quality control and bias correction,the global biases and standard deviations of the differences between observations and model simulations are also reduced.Based on the Advanced Microwave Sounding Unit A and the Microwave Humidity Sounder onboard Met Op-B,and the MWTS-2 and MWHS-2 onboard FY-3D,adding the microwave sounding data of FY-3E can further reduce the errors of analysis results and improve the global prediction skills of CMA-GFS,especially for the southern-hemisphere forecasts within 96 hours,all of which are significant at the 95% confidence level.

A Cooperated Imperialist Competitive Algorithm for Unrelated Parallel Batch Machine Scheduling Problem

This study focuses on the scheduling problem of unrelated parallel batch processing machines(BPM)with release times,a scenario derived from the moulding process in a foundry.In this process,a batch is initially formed,placed in a sandbox,and then the sandbox is positioned on a BPM formoulding.The complexity of the scheduling problem increases due to the consideration of BPM capacity and sandbox volume.To minimize the makespan,a new cooperated imperialist competitive algorithm(CICA)is introduced.In CICA,the number of empires is not a parameter,and four empires aremaintained throughout the search process.Two types of assimilations are achieved:The strongest and weakest empires cooperate in their assimilation,while the remaining two empires,having a close normalization total cost,combine in their assimilation.A new form of imperialist competition is proposed to prevent insufficient competition,and the unique features of the problem are effectively utilized.Computational experiments are conducted across several instances,and a significant amount of experimental results show that the newstrategies of CICAare effective,indicating promising advantages for the considered BPMscheduling problems.

Scientific Advances and Weather Services of the China Meteorological Administration’s National Forecasting Systems during the Beijing 2022 Winter Olympics

Since the Beijing 2022 Winter Olympics was the first Winter Olympics in history held in continental winter monsoon climate conditions across complex terrain areas,there is a deficiency of relevant research,operational techniques,and experience.This made providing meteorological services for this event particularly challenging.The China Meteorological Administration(CMA)Earth System Modeling and Prediction Centre,achieved breakthroughs in research on short-and medium-term deterministic and ensemble numerical predictions.Several key technologies crucial for precise winter weather services during the Winter Olympics were developed.A comprehensive framework,known as the Operational System for High-Precision Weather Forecasting for the Winter Olympics,was established.Some of these advancements represent the highest level of capabilities currently available in China.The meteorological service provided to the Beijing 2022 Games also exceeded previous Winter Olympic Games in both variety and quality.This included achievements such as the“100-meter level,minute level”downscaled spatiotemporal resolution and forecasts spanning 1 to 15 days.Around 30 new technologies and over 60 kinds of products that align with the requirements of the Winter Olympics Organizing Committee were developed,and many of these techniques have since been integrated into the CMA’s operational national forecasting systems.These accomplishments were facilitated by a dedicated weather forecasting and research initiative,in conjunction with the preexisting real-time operational forecasting systems of the CMA.This program represents one of the five subprograms of the WMO’s high-impact weather forecasting demonstration project(SMART2022),and continues to play an important role in their Regional Association(RA)II Research Development Project(Hangzhou RDP).Therefore,the research accomplishments and meteorological service experiences from this program will be carried forward into forthcoming highimpact weather forecasting activities.This article provides an overview and assessment of this program and the operational national forecasting systems.

A multi-scale second-order autoregressive recursive filter approach for the sea ice concentration analysis

To effectively extract multi-scale information from observation data and improve computational efficiency,a multi-scale second-order autoregressive recursive filter(MSRF)method is designed.The second-order autoregressive filter used in this study has been attempted to replace the traditional first-order recursive filter used in spatial multi-scale recursive filter(SMRF)method.The experimental results indicate that the MSRF scheme successfully extracts various scale information resolved by observations.Moreover,compared with the SMRF scheme,the MSRF scheme improves computational accuracy and efficiency to some extent.The MSRF scheme can not only propagate to a longer distance without the attenuation of innovation,but also reduce the mean absolute deviation between the reconstructed sea ice concentration results and observations reduced by about 3.2%compared to the SMRF scheme.On the other hand,compared with traditional first-order recursive filters using in the SMRF scheme that multiple filters are executed,the MSRF scheme only needs to perform two filter processes in one iteration,greatly improving filtering efficiency.In the two-dimensional experiment of sea ice concentration,the calculation time of the MSRF scheme is only 1/7 of that of SMRF scheme.This means that the MSRF scheme can achieve better performance with less computational cost,which is of great significance for further application in real-time ocean or sea ice data assimilation systems in the future.

Bridging element-free Galerkin and pluri-Gaussian simulation for geological uncertainty estimation in an ensemble smoother data assimilation framework

The facies distribution of a reservoir is one of the biggest concerns for geologists,geophysicists,reservoir modelers,and reservoir engineers due to its high importance in the setting of any reliable decisionmaking/optimization of field development planning.The approach for parameterizing the facies distribution as a random variable comes naturally through using the probability fields.Since the prior probability fields of facies come either from a seismic inversion or from other sources of geologic information,they are not conditioned to the data observed from the cores extracted from the wells.This paper presents a regularized element-free Galerkin(R-EFG)method for conditioning facies probability fields to facies observation.The conditioned probability fields respect all the conditions of the probability theory(i.e.all the values are between 0 and 1,and the sum of all fields is a uniform field of 1).This property achieves by an optimization procedure under equality and inequality constraints with the gradient projection method.The conditioned probability fields are further used as the input in the adaptive pluri-Gaussian simulation(APS)methodology and coupled with the ensemble smoother with multiple data assimilation(ES-MDA)for estimation and uncertainty quantification of the facies distribution.The history-matching of the facies models shows a good estimation and uncertainty quantification of facies distribution,a good data match and prediction capabilities.

Lateral root elongation in maize is related to auxin synthesis and transportation mediated by N metabolism under a mixed NO_(3)^(–) and NH_(4)^(+) supply

A mixed nitrate (NO_(3)^(–)) and ammonium (NH_(4)^(+)) supply can promote root growth in maize (Zea mays),however,the changes in root morphology and the related physiological mechanism under different N forms are still unclear.Here,maize seedlings were grown hydroponically with three N supplied in three different forms (NO_(3)^(–)only,75/25 NO_(3)^(–)/NH_(4)^(+)and NH_(4)^(+)only).Compared with sole NO_(3)^(–)or NH_(4)^(+),the mixed N supply increased the total root length of maize but did not affect the number of axial roots.The main reason was the increased total lateral root length,while the average lateral root (LR) length in each axle was only slightly increased.In addition,the average LR density of 2nd whorl crown root under mixed N was also increased.Compared with sole nitrate,mixed N could improve the N metabolism of roots (such as the N influx rate,nitrate reductase (NR) and glutamine synthase (GS)enzyme activities and total amino content of the roots).Experiments with exogenously added NR and GS inhibitors suggested that the increase in the average LR length under mixed N was related to the process of N assimilation,and whether the NR mediated NO synthesis participates in this process needs further exploration.Meanwhile,an investigation of the changes in root-shoot ratio and carbon (C) concentration showed that C transportation from shoots to roots may not be the key factor in mediating lateral root elongation,and the changes in the sugar concentration in roots further proved this conclusion.Furthermore,the synthesis and transportation of auxin in axial roots may play a key role in lateral root elongation,in which the expression of ZmPIN1B and ZmPIN9 may be involved in this pathway.This study preliminarily clarified the changes in root morphology and explored the possible physiological mechanism under a mixed N supply in maize,which may provide some theoretical basis for the cultivation of crop varieties with high N efficiency.

Forty-year investigation of wave power in energetic region of Persian Gulf in Iranian territorial waters by using short-term and new long-term stability assessment parameters

The wave power in high potential area of the northern Persian Gulf(near to Iranian coastal areas)is assessed by taking into account the temporal and spatial distributions of wave power for a period of forty years.For this purpose,assimilated wind data of European Centre for Medium-Range Weather Forecasting Interim Reanalysis(ERA-Interim),and hydrography data of General Bathymetric Chat of the oceans(GEBCO)are used as SWAN model.Seven locations are investigated in the study area by considering the amount of coefficient of variation,the amount of average annual power,and the short-term and a new long-term(decadal variability index)power stability assessment parameters.The results showed more stability in the eastern parts of the study area and concluded that a narrow line between the point which is in the middle and another point which is in the eastern middle part of the study area may be the best locations for more investigation and the feasibility study for energy converter farms.Also,it is found that the middle part of the study region with about 2.5 kW/m power is the most energetic area.It is concluded that the dominant direction of wave power distribution in all points is the northeast,and this dominant direction has not changed during the forty-year period.It is observed that the mean annual energy increases with a slight slope in the total 40 a,and this increasing trend is more obvious in the fourth decade.Although it is observed that the wave power of the second decade has the most stability and the least variation,the wave power in the fourth decade has the most variation.Moreover,the results showed that the study region’s wave power increase by approximately a mean change rate of 0.027 kW/(m·a);and the maximum change rate of wave power was in the northwest part and the minimum change rate of wave power was in the southeast part which were about 0.036 kW/(m·a)and 0.014 kW/(m·a),respectively.