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.     

Influence of the NAO on Wintertime Surface Air Temperature over East Asia: Multidecadal Variability and Decadal Prediction

In this paper,we investigate the influence of the winter NAO on the multidecadal variability of winter East Asian surface air temperature(EASAT)and EASAT decadal prediction.The observational analysis shows that the winter EASAT and East Asian minimum SAT(EAmSAT)display strong in-phase fluctuations and a significant 60-80-year multidecadal variability,apart from a long-term warming trend.The winter EASAT experienced a decreasing trend in the last two decades,which is consistent with the occurrence of extremely cold events in East Asia winters in recent years.The winter NAO leads the detrended winter EASAT by 12-18 years with the greatest significant positive correlation at the lead time of 15 years.Further analysis shows that ENSO may affect winter EASAT interannual variability,but does not affect the robust lead relationship between the winter NAO and EASAT.We present the coupled oceanic-atmospheric bridge(COAB)mechanism of the NAO influences on winter EASAT multidecadal variability through its accumulated delayed effect of~15 years on the Atlantic Multidecadal Oscillation(AMO)and Africa-Asia multidecadal teleconnection(AAMT)pattern.An NAO-based linear model for predicting winter decadal EASAT is constructed on the principle of the COAB mechanism,with good hindcast performance.The winter EASAT for 2020-34 is predicted to keep on fluctuating downward until~2025,implying a high probability of occurrence of extremely cold events in coming winters in East Asia,followed by a sudden turn towards sharp warming.The predicted 2020/21 winter EASAT is almost the same as the 2019/20 winter.     

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.     

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.     

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.     

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.     

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.     

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

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

Extreme Cold Events from East Asia to North America in Winter 2020/21: Comparisons,Causes, and Future Implications

Three striking and impactful extreme cold weather events successively occurred across East Asia and North America during the mid-winter of 2020/21.These events open a new window to detect possible underlying physical processes.The analysis here indicates that the occurrences of the three events resulted from integrated effects of a concurrence of anomalous thermal conditions in three oceans and interactive Arctic-lower latitude atmospheric circulation processes,which were linked and influenced by one major sudden stratospheric warming(SSW).The North Atlantic warm blob initiated an increased poleward transient eddy heat flux,reducing the Barents-Kara seas sea ice over a warmed ocean and disrupting the stratospheric polar vortex(SPV)to induce the major SSW.The Rossby wave trains excited by the North Atlantic warm blob and the tropical Pacific La Nina interacted with the Arctic tropospheric circulation anomalies or the tropospheric polar vortex to provide dynamic settings,steering cold polar air outbreaks.The long memory of the retreated sea ice with the underlying warm ocean and the amplified tropospheric blocking highs from the midlatitudes to the Arctic intermittently fueled the increased transient eddy heat flux to sustain the SSW over a long time period.The displaced or split SPV centers associated with the SSW played crucial roles in substantially intensifying the tropospheric circulation anomalies and moving the jet stream to the far south to cause cold air outbreaks to a rarely observed extreme state.The results have significant implications for increasing prediction skill and improving policy decision making to enhance resilience in“One Health,One Future”.     

自动气象站出现故障编发报文的几种方法

1引言伊春自动气象站是国家基本站(现改为一级站),每天担负8次天气报,资料参加全球交换。当计算机出现故障时,天气报是否能在规定的时间内准确的编发出是至关重要的。出现故障时,在短时间内准确的手工编发出报文难度很大,容易出现错情或造成迟报等重大错情。为了避免这种情况的发生,总结了以下几种方法进行编发报,降低了出错的几率。     

Causes of the Extreme Hot Midsummer in Central and South China during 2017:Role of the Western Tropical Pacific Warming

This study investigates why an extreme hot midsummer occurred in Central and South China(CSC) during 2017. It is shown that the western North Pacific subtropical high(WNPSH) was abnormally intensified and westward-extending,resulting in anomalous high pressure and consequent extreme heat over CSC. The abnormal WNPSH was favored by the warming of the western tropical Pacific(WTP), which was unrelated to ENSO and manifested its own individual effect.The WTP warming enhanced the convection in-situ and led to anomalous high pressure over CSC via a local meridional circulation. The influence of the WTP was confirmed by CAM4 model experiments. A comparison between the 2017 midsummer and 2010 midsummer(with a stronger WNPSH but weaker extreme heat) indicated that the influence of the WNPSH on extreme heat can be modulated by the associated precipitation in the northwestern flank.The role of the WTP was verified by regression analyses on the interannual variation of the WTP sea surface temperature anomaly(SSTA). On the other hand, the WTP has undergone prominent warming during the past few decades, resulting from decadal to long-term changes and favoring extreme warm conditions. Through a mechanism similar to the interannual variation, the decadal to long-term changes have reinforced the influence of WTP warming on the temperature over CSC,contributing to the more frequent hot midsummers recently. It is estimated that more than 50% of the temperature anomaly over CSC in the 2017 midsummer was due to the WTP warming, and 40% was related to the decadal to long-term changes of the WTP SSTA.     

The Distribution and Uncertainty Quantification of Wind Profile in the Stochastic General Ekman Momentum Approximation Model

The general Ekman momentum approximation boundary-layer model(GEM) can be effectively used to describe the physical processes of the boundary layer. However, eddy viscosity, which is an approximated value, can lead to uncertainty in the solutions. In this paper, stochastic eddy viscosity is taken into consideration in the GEM, and generalized polynomial chaos is used to quantify the uncertainty. The goal of uncertainty quantification is to investigate the effects of uncertainty in the eddy viscosity on the model and to subsequently provide reliable distribution of simulation results. The performances of the stochastic eddy viscosity and generalized polynomial chaos method are validated based on three different types of eddy viscosities, and the results are compared based on the Monte Carlo method. The results indicate that the generalized polynomial chaos method can be accurately and efficiently used in uncertainty quantification for the GEM with stochastic eddy viscosity.     

Stratospheric Ozone-induced Cloud Radiative Effects on Antarctic Sea Ice

Recent studies demonstrate that the Antarctic Ozone Hole has important influences on Antarctic sea ice.While most of these works have focused on effects associated with atmospheric and oceanic dynamic processes caused by stratospheric ozone changes,here we show that stratospheric ozone-induced cloud radiative effects also play important roles in causing changes in Antarctic sea ice.Our simulations demonstrate that the recovery of the Antarctic Ozone Hole causes decreases in clouds over Southern Hemisphere(SH)high latitudes and increases in clouds over the SH extratropics.The decrease in clouds leads to a reduction in downward infrared radiation,especially in austral autumn.This results in cooling of the Southern Ocean surface and increasing Antarctic sea ice.Surface cooling also involves ice-albedo feedback.Increasing sea ice reflects solar radiation and causes further cooling and more increases in Antarctic sea ice.     

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.     

关于三次观测站02时地面温度计算的探讨

在气象观测规范107页和新规范111页规定无自记仪器的项目,三次观测站02时地面温度应用(当日地面最低温度 前一日20时地面温度)÷2求得。在以往的工作中,发现在特殊的天气条件下,利用规范规定计算出的02时地面温度低于当日夜间地面最低温度,有些还低的很多,通过对乌伊岭气象站和铁力气象站历年资料的统计,发现这种情况普遍存在。     

识别视程障碍现象中的几个问题

1引言关于轻雾、霾、浮尘、烟幕等天气现象的识别,多年积累了一定的经验。但是在遇到某些复杂情况时应当如何辨别和记录,仍有不少争议。主要有以下3点:轻雾和霾是否可利用相对湿度大小作为区别的主要依据;轻雾和霾能否互相转换,交替出现;轻雾、霾和烟雾等能否同时存在,同时记录。本地区虽同时出现次数不多,但观测员对出现时的正确识别要重视。     

Performance of the Wind Farm Parameterization Scheme Coupled with the Weather Research and Forecasting Model under Multiple Resolution Regimes for Simulating an Onshore Wind Farm

We use the Wind Farm Parameterization(WFP) scheme coupled with the Weather Research and Forecasting model under multiple resolution regimes to simulate turbulent wake dynamics generated by a real onshore wind farm and their influence at the local meteorological scale. The model outputs are compared with earlier modeling and observation studies. It is found that higher vertical and horizontal resolutions have great impacts on the simulated wake flow dynamics. The corresponding wind speed deficit and turbulent kinetic energy results match well with previous studies. In addition, the effect of horizontal resolution on near-surface meteorology is significantly higher than that of vertical resolution. The wake flow field extends from the start of the wind farm to downstream within 10 km, where the wind speed deficit may exceed 4%. For a height of 150 m or at a distance of about 25 km downstream, the wind speed deficit is around 2%. This indicates that, at a distance of more than 25 km downstream, the impact of the wind turbines can be ignored. Analysis of near-surface meteorology indicates a night and early morning warming near the surface, and increase in near-surface water vapor mixing ratio with decreasing surface sensible and latent heat fluxes. During daytime, a slight cooling near the surface and decrease in the near-surface water vapor mixing ratio with increasing surface sensible and latent heat fluxes is noticed over the wind farm area.     

Increases in Anthropogenic Heat Release from Energy Consumption Lead to More Frequent Extreme Heat Events in Urban Cities

With economic development and rapid urbanization,increases in Gross Domestic Product and population in fastgrowing cities since the turn of the 21st Century have led to increases in energy consumption.Anthropogenic heat flux released to the near-surface atmosphere has led to changes in urban thermal environments and severe extreme temperature events.To investigate the effects of energy consumption on urban extreme temperature events,including extreme heat and cold events,a dynamic representation scheme of anthropogenic heat release(AHR)was implemented in the Advanced Research version of the Weather Research and Forecasting(WRF)model,and AHR data were developed based on energy consumption and population density in a case study of Beijing,China.Two simulations during 1999?2017 were then conducted using the developed WRF model with 3-km resolution with and without the AHR scheme.It was shown that the mean temperature increased with the increase in AHR,and more frequent extreme heat events were produced,with an annual increase of 0.02?0.19 days,as well as less frequent extreme cold events,with an annual decrease of 0.26?0.56 days,based on seven extreme temperature indices in the city center.AHR increased the sensible heat flux and led to surface energy budget changes,strengthening the dynamic processes in the atmospheric boundary layer that reduce AHR heating efficiency more in summer than in winter.In addition,it was concluded that suitable energy management might help to mitigate the impact of extreme temperature events in different seasons.