AN INTEGRATION METHOD WITH FITTING CUBIC SPLINE FUNCTIONS TO A NUMERICAL MODEL OF 2ND-ORDER SPACE-TIME DIFFERENTIAL REMAINDER——FOR AN IDEAL GLOBAL SIMULATION CASE WITH PRIMITIVE ATMOSPHERIC EQUATIONS

In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.     

Classification of Northeast China Cold Vortex Activity Paths in Early Summer Based on K-means Clustering and Their Climate Impact

The classification of the Northeast China Cold Vortex(NCCV)activity paths is an important way to analyze its characteristics in detail.Based on the daily precipitation data of the northeastern China(NEC)region,and the atmospheric circulation field and temperature field data of ERA-Interim for every six hours,the NCCV processes during the early summer(June)seasons from 1979 to 2018 were objectively identified.Then,the NCCV processes were classified using a machine learning method(k-means)according to the characteristic parameters of the activity path information.The rationality of the classification results was verified from two aspects,as follows:(1)the atmospheric circulation configuration of the NCCV on various paths;and(2)its influences on the climate conditions in the NEC.The obtained results showed that the activity paths of the NCCV could be divided into four types according to such characteristics as the generation origin,movement direction,and movement velocity of the NCCV.These included the generation-eastward movement type in the east of the Mongolia Plateau(eastward movement type or type A);generation-southeast longdistance movement type in the upstream of the Lena River(southeast long-distance movement type or type B);generationeastward less-movement type near Lake Baikal(eastward less-movement type or type C);and the generation-southward less-movement type in eastern Siberia(southward less-movement type or type D).There were obvious differences observed in the atmospheric circulation configuration and the climate impact of the NCCV on the four above-mentioned types of paths,which indicated that the classification results were reasonable.     

Phase Two of the Integrative Monsoon Frontal Rainfall Experiment (IMFRE-II) over the Middle and Lower Reaches of the Yangtze River in 2020

Phase Two of the Integrative Monsoon Frontal Rainfall Experiment(IMFRE-II)was conducted over the middle and lower reaches of the Yangtze River during the period 16 June to 19 July 2020.This paper provides a brief overview of the IMFRE-II field campaign,including the multiple ground-based remote sensors,aircraft probes,and their corresponding measurements during the 2020 mei-yu period,as well as how to use these numerous datasets to answer scientific questions.The highlights of IMFRE-II are:(1)to the best of our knowledge,IMFRE-II is the first field campaign in China to use ground-based,airborne,and spaceborne platforms to conduct comprehensive observations over the middle and lower reaches of the Yangtze River;and(2)seven aircraft flights were successfully carried out,and the spectra of ice particles,cloud droplets,and raindrops at different altitudes were obtained.These in-situ measurements will provide a“cloud truth”to validate the ground-based and satellite-retrieved cloud and precipitation properties and quantitatively estimate their retrieval uncertainties.They are also crucial for the development of a warm(and/or cold)rain conceptual model in order to better understand the cloud-to-rain conversion and accretion processes in mei-yu precipitation events.Through an integrative analysis of ground-based,aircraft,and satellite observations and model simulations,we can significantly improve our cloud and precipitation retrieval algorithms,investigate the microphysical properties of cloud and precipitation,understand in-depth the formation and dissipation mechanisms of mei-yu frontal systems,and improve cloud microphysics parameterization schemes and model simulations.     

THE IMPACTS OF MADDEN-JULIAN OSCILLATION ON SPRING RAINFALL IN EAST CHINA

Phase composite analyses are conducted to investigate the possible effect of the Madden–Julian oscillation(MJO)on the spring rainfall anomalies in East China by using the Real-time Multivariate MJO(RMM)index from Australian Meteorological Bureau.The results show that the rainfall anomalies over the mid-and lower-valley of Yangtze River are positive when the MJO shifts eastward to the mid-and eastern-Indian Ocean,and anomalous precipitation over South China are positive when the MJO moves further eastward to the maritime continent,whereas spring rainfall anomalies over East China are negative in the other MJO episodes.The MJO impacts on the precipitation over East China result from the changes in large-scale atmospheric circulation as well as vorticity and water vapor transportation in the mid-and lower-troposphere.     

CHARACTERISTICS OF TROPICAL CYCLONES IMPACTING A NUCLEAR POWER PLANT IN SOUTHERN FUJIAN AND ESTIMATE OF MINIMUM CENTRAL PRESSURE

Based on tropical cyclone (TC) data for the period 1949 to 2008 and following the Gumbe-I method, Pearson-Ⅲ method and determinacy method, this article estimates the possible minimum central pressure of TCs affecting southern Fujian where a nuclear power will be located. Results show that the observed minimum central pressure of TCs agrees well with the results determined with the methods above and there is little difference between them (the minimum central pressure is 867.4 hPa and 868.1 hPa, respectively, in a 1,000-yr return period). Established with the theory of atmospheric dynamics, the determinacy method yields a result of 867.28 hPa/1000 years, only a little smaller than the result of the probability method. Because of randomicity in parameter adjustment with the Pearson-Ⅲ method whereas the determinacy method is theoretically solid and its estimates are the smallest of the three methods, it is therefore reasonable, for security and conservative concerns, to adopt the result determined with the determinacy method as the possible maximum intensity of TC (with the central pressure being 867.28 hPa in a 1,000-yr return period).     

A Comparison of Two Bulk Microphysics Parameterizations for the Study of Aerosol Impacts on an Idealized Supercell

Idealized supercell storms are simulated with two aerosol-aware bulk microphysics schemes(BMSs),the Thompson and the Chen-Liu-Reisner(CLR),using the Weather Research and Forecast(WRF)model.The objective of this study is to investigate the parameterizations of aerosol effects on cloud and precipitation characteristics and assess the necessity of introducing aerosols into a weather prediction model at fine grid resolution.The results show that aerosols play a decisive role in the composition of clouds in terms of the mixing ratios and number concentrations of liquid and ice hydrometeors in an intense supercell storm.The storm consists of a large amount of cloud water and snow in the polluted environment,but a large amount of rainwater and graupel instead in the clean environment.The total precipitation and rain intensity are suppressed in the CLR scheme more than in the Thompson scheme in the first three hours of storm simulations.The critical processes explaining the differences are the auto-conversion rate in the warm-rain process at the beginning of storm intensification and the low-level cooling induced by large ice hydrometeors.The cloud condensation nuclei(CCN)activation and auto-conversion processes of the two schemes exhibit considerable differences,indicating the inherent uncertainty of the parameterized aerosol effects among different BMSs.Beyond the aerosol effects,the fall speed characteristics of graupel in the two schemes play an important role in the storm dynamics and precipitation via low-level cooling.The rapid intensification of storms simulated with the Thompson scheme is attributed to the production of hail-like graupel.     

Verification of Subseasonal-to-Seasonal Forecasts for Major Stratospheric Sudden Warmings in Northern Winter from 1998/99 to 2012/13

This study reports verification results of hindcast data of four systems in the subseasonal-to-seasonal(S2S)prediction project for major stratospheric sudden warmings(MSSWs)in northern winter from 1998/99 to 2012/13.This report deals with average features across all MSSWs,and possible differences between two MSSW types(vortex displacement and split types).Results for the average features show that stratospheric forecast verifications,when further averaged among the four systems,are judged to be successful for lead times around 10 d or shorter.All systems are skillful for lead times around 5 d,whereas the results vary among the systems for longer lead times.A comparison between the MSSW types overall suggests larger forecast errors or lower skill for MSSWs of the vortex split type,although the differences do not have strong statistical significance for almost all cases.This limitation is likely to at least partly reflect the small sample size of the MSSWs available.     

Hydro-climatic Characteristics of Yarlung Zangbo River Basin since the Last Glacial Maximum

Global climate changes significantly impact the water condition of big rivers in glacierized high mountains. However,there is a lack of studies on hydrological changes within river basins caused by climate changes over a geological timescale due to the impossibility of direct observations. In this study, we examine the hydro-climatic variation of the Yarlung Zangbo River Basin in the Tibet Plateau since the Last Glacial Maximum(LGM) by combining δ18 O proxy records in Indian and Omani caves with the simulated Indian summer monsoon, surface temperature, precipitation, evapotranspiration and runoff via the Community Climate System Model and the reconstructed glacier coverage via the Parallel Ice Sheet Model. The mean river runoff was kept at a low level of 145 billion cubic meters per year until an abrupt increase at a rate of 8.7 million cubic meters per year in the B?lling-Aller?d interval(BA). The annual runoff reached a maximum of 250 billion cubic meters in the early Holocene and then reduced to the current value of 180 billion cubic meters at a rate of 6.4 million cubic meters per year. The low runoff in the LGM and Heinrich Stadial 1(HS1) is likely attributed to such a small contribution of precipitation to runoff and the large glacier cover. The percentage of precipitation to runoff was only 20%during the LGM and HS1. Comparison of glacier area among different periods indicates that the fastest deglaciation occurred during the late HS1, when nearly 60% of glacier area disappeared in the middle reach, 50% in the upper reach,and 30% in the lower reach. The rapid deglaciation and increasing runoff between the late HS1 and BA may have accelerated widespread ice-dam breaches and led to extreme outburst flood events. Combining local geological proxy records and regional simulations could be a useful approach for the study of paleo-hydrologic variations in big river basins.     

南京三千公尺高空之风向与天气之预测

晚近日本籐原笑平(Fujiwhora)博士于地球物理杂志发表「根据三千公尺高空等压线,以预测天气之一例证」一文,谓日本最近用三千公尺高空之等压线,作每日天气之预测,已得相当成就。氏之经验法则谓自九月以迄五月,日本太平洋沿岸,三千公尺高空之等压线,来自西南者,可形去致雨,而来自西北者,则可期晴明。此种倾向颇为显著。     

Impact of Assimilation of Radiosonde and UAV Observations from the Southern Ocean in the Polar WRF Model

Weather forecasting in the Southern Ocean and Antarctica is a challenge above all due to the rarity of observations to be assimilated in numerical weather prediction(NWP)models.As observations are expensive and logistically challenging,it is important to evaluate the benefit that additional observations could bring to NWP.Atmospheric soundings applying unmanned aerial vehicles(UAVs)have a large potential to supplement conventional radiosonde sounding observations.Here,we applied UAV and radiosonde sounding observations from an RV Polarstern cruise in the ice-covered Weddell Sea in austral winter 2013 to evaluate the impact of their assimilation in the Polar version of the Weather Research and Forecasting(Polar WRF)model.Our experiments revealed small to moderate impacts of radiosonde and UAV data assimilation.In any case,the assimilation of sounding data from both radiosondes and UAVs improved the analyses of air temperature,wind speed,and humidity at the observation site for most of the time.Further,the impact on the results of 5-day-long Polar WRF experiments was often felt over distances of at least 300 km from the observation site.All experiments succeeded in capturing the main features of the evolution of near-surface variables,but the effects of data assimilation varied between different cases.Due to the limited vertical extent of the UAV observations,the impact of their assimilation was limited to the lowermost 1?2-km layer,and assimilation of radiosonde data was more beneficial for modeled sea level pressure and near-surface wind speed.     

A COMPARATIVE ANALYSIS OF ATMOSPHERIC AND OCEANIC CONDITIONS BEFORE THE OCCURRENCE OF TWO TYPES OF EL NIO EVENTS

The atmospheric and oceanic conditions before the onset of EP El Ni?o and CP El Ni?o in nearly 30 years are compared and analyzed by using 850 hPa wind, 20℃ isotherm depth, sea surface temperature and the Wheeler and Hendon index. The results are as follows: In the western equatorial Pacific, the occurrence of the anomalously strong westerly winds of the EP El Ni?o is earlier than that of the CP El Ni?o. Its intensity is far stronger than that of the CP El Ni?o. Two months before the El Ni?o, the anomaly westerly winds of the EP El Ni?o have extended to the eastern Pacific region, while the westerly wind anomaly of the CP El Ni?o can only extend to the west of the dateline three months before the El Ni?o and later stay there. Unlike the EP El Ni?o, the CP El Ni?o is always associated with easterly wind anomaly in the eastern equatorial Pacific before its onset. The thermocline depth anomaly of the EP El Ni?o can significantly move eastward and deepen. In addition, we also find that the evolution of thermocline is ahead of the development of the sea surface temperature for the EP El Ni?o. The strong MJO activity of the EP El Ni?o in the western and central Pacific is earlier than that of the CP El Ni?o. Measured by the standard deviation of the zonal wind square, the intensity of MJO activity of the EP El Ni?o is significantly greater than that of the CP El Ni?o before the onset of El Ni?o.     

ANALYSIS OF “BIMODAL PATTERN” STORM ACTIVITY CHARACTERISTICS OF THE BAY OF BENGAL

Storms that occur at the Bay of Bengal (BoB) are of a bimodal pattern, which is different from that of the other sea areas. By using the NCEP, SST and JTWC data, the causes of the bimodal pattern storm activity of the BoB are diagnosed and analyzed in this paper. The result shows that the seasonal variation of general atmosphere circulation in East Asia has a regulating and controlling impact on the BoB storm activity, and the “bimodal period” of the storm activity corresponds exactly to the seasonal conversion period of atmospheric circulation. The minor wind speed of shear spring and autumn contributed to the storm, which was a crucial factor for the generation and occurrence of the “bimodal pattern” storm activity in the BoB. The analysis on sea surface temperature (SST) shows that the SSTs of all the year around in the BoB area meet the conditions required for the generation of tropical cyclones (TCs). However, the SSTs in the central area of the bay are higher than that of the surrounding areas in spring and autumn, which facilitates the occurrence of a “two-peak” storm activity pattern. The genesis potential index (GPI) quantifies and reflects the environmental conditions for the generation of the BoB storms. For GPI, the intense low-level vortex disturbance in the troposphere and high-humidity atmosphere are the sufficient conditions for storms, while large maximum wind velocity of the ground vortex radius and small vertical wind shear are the necessary conditions of storms.     

Dominant cloud microphysical processes of a torrential rainfall event in Sichuan,China

High-resolution numerical simulation data of a rainstorm triggering debris flow in Sichuan Province of China simulated by the Weather Research and Forecasting (WRF) Model were used to study the dominant cloud microphysical processes of the torrential rainfall. The results showed that: (1) In the strong precipitation period, particle sizes of all hydrometeors increased, and mean-mass diameters of graupel increased the most significantly, as compared with those in the weak precipitation period; (2) The terminal velocity of raindrops was the strongest among all hydrometeors, followed by graupel’s, which was much smaller than that of raindrops. Differences between various hydrometeors’ terminal velocities in the strong precipitation period were larger than those in the weak precipitation period, which favored relative motion, collection interaction and transformation between the particles. Absolute terminal velocity values of raindrops and graupel were significantly greater than those of air upward velocity, and the stronger the precipitation was, the greater the differences between them were; (3) The orders of magnitudes of the various hydrometeors’ sources and sinks in the strong precipitation period were larger than those in the weak precipitation period, causing a difference in the intensity of precipitation. Water vapor, cloud water, raindrops, graupel and their exchange processes played a major role in the production of the torrential rainfall, and there were two main processes via which raindrops were generated: abundant water vapor condensed into cloud water and, on the one hand, accretion of cloud water by rain water formed rain water, while on the other hand, accretion of cloud water by graupel formed graupel, and then the melting of graupel formed rain water.     

An assessment of the predictability of the East Asian Subtropical Westerly Jet based on TIGGE data

The predictability of the position, spatial coverage and intensity of the East Asian subtropical westerly jet(EASWJ) in the summers of 2010 to 2012 was examined for ensemble prediction systems(EPSs) from four representative TIGGE centers,including the ECMWF, the NCEP, the CMA, and the JMA. Results showed that each EPS predicted all EASWJ properties well, while the levels of skill of all EPSs declined as the lead time extended. Overall, improvements from the control to the ensemble mean forecasts for predicting the EASWJ were apparent. For the deterministic forecasts of all EPSs, the prediction of the average axis was better than the prediction of the spatial coverage and intensity of the EASWJ. ECMWF performed best, with a lead of approximately 0.5–1 day in predictability over the second-best EPS for all EASWJ properties throughout the forecast range. For probabilistic forecasts, differences in skills among the different EPSs were more evident in the earlier part of the forecast for the EASWJ axis and spatial coverage, while they departed obviously throughout the forecast range for the intensity. ECMWF led JMA by about 0.5–1 day for the EASWJ axis, and by about 1–2 days for the spatial coverage and intensity at almost all lead times. The largest lead of ECMWF over the relatively worse EPSs, such as NCEP and CMA, was approximately 3–4 days for all EASWJ properties. In summary, ECMWF showed the highest level of skill for predicting the EASWJ, followed by JMA.     

IMPACT OF ATMOSPHERIC LOW-FREQUENCY OSCILLATION ON THE IMMIGRATION OF NILAPARVATA LUGENS(STL) IN HUNAN AND JIANGXI PROVINCES

Various features of the atmospheric environment affect the number of migratory insects, besides their initial population. However, little is known about the impact of atmospheric low-frequency oscillation(10 to 90 days) on insect migration. A case study was conducted to ascertain the influence of low-frequency atmospheric oscillation on the immigration of brown planthopper, Nilaparvata lugens(Stl), in Hunan and Jiangxi provinces. The results showed the following:(1) The number of immigrating N. lugens from April to June of 2007 through 2016 mainly exhibited a periodic oscillation of 10 to 20 days.(2) The 10-20 d low-frequency number of immigrating N. lugens was significantly correlated with a low-frequency wind field and a geopotential height field at 850 h Pa.(3) During the peak phase of immigration, southwest or south winds served as a driving force and carried N. lugens populations northward, and when in the back of the trough and the front of the ridge, the downward airflow created a favorable condition for N. lugens to land in the study area. In conclusion, the northward migration of N. lugens was influenced by a low-frequency atmospheric circulation based on the analysis of dynamics. This study was the first research connecting atmospheric low-frequency oscillation to insect migration.     

中国冰冻圈变化的适应研究：进展与展望

冰冻圈变化的适应研究是冰冻圈科学领域的新兴研究方向,是当今自然科学与社会科学交叉融合跨学科集成研究的典型代表。起步于2007年的中国冰冻圈变化适应研究,经历了早期的探索,研究重点由评价脆弱性发展为量化冰冻圈变化的影响,形成以影响/风险—脆弱性—适应全链条的完善的研究体系,研究方法突破传统的指标体系赋权法的不足,初步实现了定量化,有机结合影响/风险、脆弱性、适应三方面的研究结果,使冰冻圈变化的适应措施由偏重宏观性、普适性开始转向更有针对性。未来中国冰冻圈变化的适应研究应拓展、完善和深化现有的理论体系,构建冰冻圈与社会经济耦合模型,科学量化冰冻圈全要素变化的影响,建立不同利益相关者与科学家共同参与的研究新模式,科学有效应对与适应冰冻圈变化及其影响。     

淮河流域水文极值预测模型研究

为探索气候变化影响下水文极值的非平稳性和预测方法,建立了水文极值非平稳广义极值（GEV）分布的统计预测模型。利用1952-2010年淮河上游流域累计面雨量和流量年最大值资料、同期500 hPa环流特征量资料以及17个CMIP5模式对环流特征量的模拟结果,筛选出对水文极值影响显著的年平均北半球极涡强度指数作为GEV分布参数的预测因子。分析了在RCP2.6、RCP4.5和RCP8.5情景下2006-2050年淮河上游流域水文极值对气候变化的响应。结果表明,10年以下与10年以上重现期的水文极值在非平稳过程中呈现前者下降而后者上升的相反变化趋势;多模型预测的集合平均在未来情景中均呈现上升趋势,情景排放量越大增幅越大,重现期越长增幅也越大。与极值的常态相比,极值的极端态更易受气候变化影响。     

A political ecology of maladaptation: Insights from a Gramscian theory of the State

This article analyzes a political process in the aftermath of a disaster and explains why and how it led to maladaptation. Grounding Gramsci’s theory of the State on a case of post-disaster response to a fatal mudslide in the city of Sarno in Italy, this research argues that, under certain conditions, civil society and the ruling classes may coalesce to produce policies that are maladaptive. We unpack the mechanisms through which consent was reproduced in Sarno, and show how the claims of civil society were articulated and fused with the hegemonic goals of capital circulation and economic growth, reaffirming a view of government as only a provider of safety. A Gramscian treatment of the State as a process, and not as a thing, highlights that the main barrier to adaptation is not the lack of techno-managerial solutions. It is the lack of political struggle around the social reconfiguration of the logic and functions of the State.     

ECMWF极端天气指数在华东台风降水预报中的应用和检验

选取2015—2018年影响华东地区的13个台风个例,分析降水极端天气指数EFI(Extreme Forecast Index)和SOT(Shift of Tails)与台风降水之间的统计关系。结果表明：EFI和SOT与降水气候百分位之间存在明显的正相关关系。EFI和SOT越大,强降水发生概率越高。随着预报时效的增加,EFI和SOT指数对暴雨和大暴雨的预报效果逐渐变差。对于短期(72 h以内的时效),EFI预报技巧优于SOT,而随着预报时效的延长,SOT的预报技巧逐渐接近并超过EFI。以TS评分最大为标准兼顾合理的预报偏差,得到两种极端天气指数不同预报时效、不同等级暴雨的预报阈值。总体而言,事件越极端,EFI和SOT的预报阈值越大,对于暴雨和大暴雨,EFI指数的预报阈值随着预报时效的延长有减小趋势,而SOT的预报阈值基本保持不变。在台风极端降水预报中,EFI和SOT可以作为EC定量降水预报的补充,有助于减少强降水的漏报,并提早发出预警信息。     

Recent Advances in Predictability Studies in China （1999-2002）

Since the last International Union of Geodesy and Geophysics (IUGG) General Assembly (1999), the predictability studies in China have made further progress during the period of 1999-2002. Firstly, three predictability sub-problems in numerical weather and climate prediction are classified, which are concerned with the maximum predictability time, the maximum prediction error, and the maximum allowable initial error, and then they are reduced into three nonlinear optimization problems. Secondly, the concepts of the nonlinear singular vector (NSV) and conditional nonlinear optimal perturbation (CNOP) are proposed,which have been utilized to study the predictability of numerical weather and climate prediction. The results suggest that the nonlinear characteristics of the motions of atmosphere and oceans can be revealedby NSV and CNOP. Thirdly, attention has also been paid to the relations between the predictability and spatial-temporal scale, and between the model predictability and the machine precision, of which the investigations disclose the importance of the spatial-temporal scale and machine precision in the study of predictability. Also the cell-to-cell mapping is adopted to analyze globally the predictability of climate,which could provide a new subject to the research workers. Furthermore, the predictability of the summer rainfall in China is investigated by using the method of correlation coefficients. The results demonstrate that the predictability of summer rainfall is different in different areas of China. Analysis of variance, which is one of the statistical methods applicable to the study of predictability, is also used to study the potential predictability of monthly mean temperature in China, of which the conclusion is that the monthly mean temperature over China is potentially predictable at a statistical significance level of 0.10. In addition,in the analysis of the predictability of the T106 objective analysis/forecasting field, the variance and the correlation coefficient are calculated to explore the distribution characteristics of the mean-square errors.Finally, the predictability of short-term climate prediction is investigated by using statistical methods or numerical simulation methods. It is demonstrated that the predictability of short-term climate in China depends not only on the region of China being investigated, but also on the time scale and the atmospheric internal dynamical process.