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.     

EFFECTS OF PACIFIC SSTA ON SUMMER PRECIPITATION OVER EASTERN CHINA, PART Ⅰ: OBSERVATIONAL ANALYSIS

With the methods of REOF (Rotated Empirical Orthogonal Function), the summer precipitation from 43 stations over eastern China for the 1901 - 2000 period was examined. The results show that South China and Southwest China, the middle and lower reaches of Changjiang River, North China and the southwestern of Northeast China are the three main areas of summer rainfall anomaly. Furthermore, correlation analysis is used in three time series of three mostly summer rainfall modes and four seasonal Pacific SSTA (Sea Surface Temperature Anomaly), and the results suggest that the Pacific SSTA which notably causes the summer rainfall anomaly over eastern China are the SSTA of the preceding winter over Kuroshio region of Northwest Pacific, SSTA of the preceding spring in the eastern and central equatorial Pacific, and SSTA of the current summer in the central region of middle latitude. The relationship between summer precipitation over eastern China and SSTA of Pacific key regions was further verified by SVD (Singular Value Decomposition) analysis.The composite analysis was used to analyze the features of atmospheric general circulation in the years of positive and negative precipitation anomaly. Its results were used to serve as the base of numerical simulation analysis.     

EFFECTS OF PACIFIC SSTA ON SUMMER PRECIPITATION OVER EASTERN CHINA, PART II: NUMERICAL SIMULATIONS

Based on an observational analysis, seven numerical experiments are designed to study the impacts of Pacific SSTA on summer precipitation over eastern China and relevant physical mechanism by NCAR CCM3. The numerical simulation results show that preceding winter SSTA in the Kuroshio region leads to summer precipitation anomaly over the Yangtze River valleys by modifying atmospheric general circulation over eastern Asia and middle-high latitude. West Pacific subtropical high is notably affected by preceding spring SSTA over the middle and east of Equator Pacific; SSTA of the central region of middle latitude in the corresponding period causes the summer rainfall anomaly over eastern China so as to trigger the atmospheric Eurasia-Pacific teleconnection pattern.     

EFFECTS OF PACIFIC SSTA ON SUMMER PRECIPITATION OVER EASTERN CHINA, PART Ⅱ: NUMERICAL SIMULATIONS

Based on an observational analysis, seven numerical experiments are designed to study the impacts of Pacific SSTA on summer precipitation over eastem China and relevant physical mechanism by NCAR CCM3. The numerical simulation results show that preceding winter SSTA in the Kuroshio region leads to summer precipitation anomaly over the Yangtze River valleys by modifying atmospheric general circulation over eastern Asia and middle-high latitude. West Pacific subtropical high is notably affected by preceding spring SSTA over the middle and east of Equator Pacific; SSTA of the central region of middle latitude in the corresponding period causes the summer rainfall anomaly over eastern China so as to trigger the atmospheric Eurasia-Pacific teleconnection pattern.     

EFFECTS OF PACIFIC SSTA ON SUMMER PRECIPITATION OVER EASTERN CHINA, PART Ⅱ: NUMERICAL SIMULATIONS

Based on an observational analysis, seven numerical experiments are designed to study the impacts of Pacific SSTA on summer precipitation over eastem China and relevant physical mechanism by NCAR CCM3. The numerical simulation results show that preceding winter SSTA in the Kuroshio region leads to summer precipitation anomaly over the Yangtze River valleys by modifying atmospheric general circulation over eastern Asia and middle-high latitude. West Pacific subtropical high is notably affected by preceding spring SSTA over the middle and east of Equator Pacific; SSTA of the central region of middle latitude in the corresponding period causes the summer rainfall anomaly over eastern China so as to trigger the atmospheric Eurasia-Pacific teleconnection pattern.     

春季热带南大西洋海表温度异常与夏季亚澳季风区降水异常的联系

利用1979—2019年Hadley中心的海表温度资料、GPCP的降水资料以及NCEP-DOE的再分析资料等,分析了北半球春季热带南大西洋海表温度异常与北半球夏季亚澳季风区降水异常的联系。研究表明,北半球春季热带南大西洋海表温度异常与随后夏季热带西太平洋到南海（澳大利亚东侧海域到热带东印度洋）地区的降水异常为显著负相关（正相关）关系。北半球春季热带南大西洋的海表温度正异常可以引起热带大西洋和热带太平洋间的异常垂直环流,其中异常上升支（下沉支）位于热带大西洋（热带中太平洋）。热带中太平洋的异常下沉气流和低层辐散气流引起热带中西太平洋低层的异常东风,后者有利于热带中东太平洋海表温度出现负异常。通过Bjerknes正反馈机制,热带中东太平洋海表温度异常从北半球春季到夏季得到发展。热带中东太平洋海表温度负异常激发的Rossby波使得北半球夏季热带西太平洋低层出现一对异常反气旋。此时,850 hPa上热带西太平洋到海洋性大陆地区为显著的异常东风,有利于热带西太平洋到南海（澳大利亚东侧海域到热带东印度洋）地区出现异常的水汽辐散（辐合）,导致该地区降水减少（增加）。     

夏半年热带太平洋中部型海温异常与热带印度洋海盆模对同期中国东部降水的共同影响

利用1961—2013年NCEP/NCAR再分析资料和Had ISST月平均海表面温度资料,分析了夏半年热带太平洋中部型海温异常与热带印度洋海盆模(Indian Ocean Basin M ode,IOBM)的特征,并研究了不同位相配置时二者对同期中国东部气候的共同影响。结果表明:1)太平洋中部型海温异常指数与印度洋海盆模指数几乎相互独立。太平洋中部型海温异常与IOBM同位相变化(记为PPNN事件)和反位相变化(记为PNNP事件)时,热带印太地区海温异常分别呈三级型和偶极型分布。2)不同位相配置对中国东部地区降水异常的影响及其影响机制存在显著差异:当发生PPNN事件时,水汽从海洋性大陆(Maritime Continent,MC)地区向江淮流域输送;热带海温异常引起大气产生Gill型响应,维持了中国东部的环流异常;M C地区通过经向三圈异常垂直环流引起江淮流域降水异常增多。发生PNNP事件时,Gill型环流响应中心西移,长江流域降水偏少,水汽辐散;同时MC地区对流层低层准定常Rossby波能传播也有利于长江流域扰动的维持。这些结果对深刻认识中国东部地区夏半年降水异常成因和印度洋/太平洋海温异常不同分布的作用具有重要意义。     

中国东部夏季分区降水对海温异常响应特征的研究

用刘征宇最新发展的广义平衡反馈方法结合经验正交分析和旋转经验正交分析法,研究中国东部6个分区夏季降水异常对各海盆海温异常(SSTA)模态的响应特征,探讨了东北区降水异常对SSTA的响应机制。结果表明:各分区降水异常对SSTA的响应特征有明显差异。东北区降水异常对热带太平洋SSTA的回应显著,当SSTA表现为类似El Niño模时,该区夏季降水增加;江淮区降水异常同时受到中纬度以及热带SSTA的强迫,对热带太平洋SSTA的类似La Niña Modoki模、北大西洋SSTA三极型模等回应显著;西南地区的降水异常主要与中纬度SSTA有密切关联。对流层高度场对热带太平洋SSTA类似El Niño模的直接回应在热带太平洋上空为显著正异常。通过东亚—太平洋型遥相关波列使东北地区上空高度场出现负异常,并在其西侧形成负异常中心,对流层低层为气旋性环流异常。这一环流回应与东北地区夏季降水偏多时的环流特征相近。     

河北夏季降水异常的变化及其与环流和海温异常的可能联系

利用河北省及周边73个台站1961-2010年夏季降水资料、NCEP／NCAR再分析资料、NOAA太平洋海温资料,应用SVD等多种统计方法探讨了河北夏季降水异常变化及其与环流和同期及前期太平洋海表温度异常（SSTA）的可能联系.研究结果表明：河北省夏季降水异常存在全省一致和东北-西南相反两种模态;夏季欧亚地区500 hPa高度场表现为从西到东的负、正、负、正的遥相关型时,河北降水偏多,最明显的区域主要位于河北省的东部;夏季东北部多雨、其他地区（保定、邢台和衡水）干旱与北太平洋SSTA的PDO正位相、赤道中东太平洋SSTA的El Niño位相具有较好的耦合关系;最大协方差分析（MCA）的诊断显示,前期冬季北太平洋SSTA类似PDO正位相且赤道中东太平洋SSTA出现El Niño型态,对河北省夏季多雨（尤其是东北部地区）的预测具有较好的指示意义.夏季850 hPa上在我国东部地区表现出明显的偏北（南）风距平,不（有）利于南方暖湿水汽向北的输送,是造成河北夏季降水异常少（多）的主要因素.     

ENSO对中国夏季降水可预测性变化的研究

众多研究表明,ENSO对东亚夏季风尤其是中国夏季降水存在很大影响,已成为中国夏季降水首要的预测因子。传统的预测模型认为,当前期ENSO为暖位相状态时,夏季中国主要雨带位置偏南,长江流域降水偏多;反之,当前期ENSO为冷位相状态时,夏季中国主要雨带位置偏北,长江流域降水偏少。基于1951—2003年中国160站月平均降水资料和同时段的NOAA ERSST海表温度资料,讨论了中国夏季降水和前冬Nino3区海温关系的年代际变化。分析结果显示,近20年来二者相关性已大大衰减。作为中国夏季降水的主要预测指标,ENSO的指示意义也相应减弱。在1951—1974年,依据前冬Nino3区SSTA预测夏季降水符号准确率在67%以上的站数有43站,但在1980—2003年,同样准确率的站数只有15站。在前一个研究时段,这43站呈区域性分布于东北地区、黄河和长江流域,但后一个研究时段内的15站分布分散,不利于区域性预测。相关分析结果表明,在20世纪70年代中期之前,当前冬赤道东太平洋海温偏高时,华北和江南南部的多数测站夏季降水偏多,淮河流域降水偏少,同时梅雨开始偏晚。反之,当前冬赤道东太平洋海温偏低时,华北和江南南部夏季降水易偏少,淮河流域降水则偏多,同时梅雨开始偏早。但在20世纪80年代之后,上述对应关系较难成立。因此,在汛期预测业务中参考ENSO的作用时必须充分考虑年代际背景的差异。     

多因子和多尺度合成中国夏季降水预测模型及预报试验

根据青藏高原60个站平均的月积雪深度、热带太平洋Nino 3区月海温和中国160个站月降水量等资料,用小波变换和相关分析,分析了1958～1998年秋冬季青藏高原异常雪盖与El Nino-南方涛动(ENSO)的关系、多时间尺度变化的特征及其与中国夏季降水的相关型式.并取青藏高原积雪和Nino 3区海温的年际变化、年代际变化和线性趋势三种不同时间尺度的小波分量作为预报因子,对我国夏季降水距平作线性回归,建立了相应的预测模型.最后,利用1999～2002年的独立资料进行了预报试验,并在2003年和2004年应用于实际预报.研究表明,青藏高原雪盖与ENSO这两个物理因子彼此具有一定的独立性.它们都是多时间尺度现象,并与中国夏季降水有较好的关系.在不同时间尺度上不仅有不同的相关型式,而且相对贡献也有变化.回归预测模型的拟合情况和预报试验表明,综合考虑前期秋冬季青藏高原雪盖和ENSO这两个物理因子的年际变化、年代际变化和线性趋势作为预报因子建立的预测我国夏季降水距平分布的模型,有一定的预报能力.