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.     

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.     

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.     

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.     

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

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

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

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

RegCM3对全球变暖背景下中国东部季风降水和雨型变化的模拟

The authors investigate possible changes of monsoon rainfall and associated seasonal (June-JulyAugust) anomaly patterns over eastern China in the late 21st century under the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emission Scenarios (SRES) A2 emission scenario as simulated by a high-resolution regional climate model (RegCM3) nested in a general circulation model (FvGCM/CCM3).Two sets of multi-decadal simulations are performed at 20-km grid spacing for present day and future climate conditions.Results show that the RegCM3 reproduces the mean rainfall distribution;however the evolution of the monsoon rain belt from South China to North China is not well simulated.Concerning the rain pattern classifications,RegCM3 overestimates the occurrence of Pattern 1 (excessive rainfall in northern China) and underestimates that of Pattern 2 (increased rainfall over the Huai River basin).Under future climate conditions,RegCM3 projects less occurrence of Pattern 1,more of Pattern 2,and little change of Pattern 3 (rainfall increase along the Yangtze River).These results indicate that there might be increased rainfall over the Huai-Yellow River area and reduced rainfall over North China in the future,while rainfall over the lower reaches of the Yangtze River basin is not modified significantly.Uncertainties exist in the present study are also discussed.     

NUMERICAL SIMULATION OF CLOUD MICROPHYSICAL CHARACTERISTICS OF LANDFALL TYPHOON KROSA

In this study,the super typhoon KROSA(2007)was simulated using a mesoscale numerical model Global and Regional Assimilation and Prediction System(GRAPES)with a two-moment mixed-phase microphysics scheme.Local rainfall observations,radar and satellite data were also used to analyze the precipitation structure and microphysical features.It was shown that low-level jets and unstable temperature stratification provided this precipitation process with favorable weather condition.Heavy rainfall centers were located in the north and east part of KROSA with the maxima of 6-hourly total rainfall during the simulation more than 100 mm.The quantities of column solid water and column liquid water were generally equivalent,indicating the important role of ice phase in precipitation formation.Results of CloudSat showed that strong convection occurred in the eyewall around the cyclonic center.According to the simulation results,heavy precipitation in the northeast part of the typhoon was mainly triggered by convective clouds,accompanied by the strongest updraft under the melting level.In the southwest part of KROSA,precipitation intensity was rather homogeneous.The ascending center occurred in high-level cold clouds,favoring the formation and growth of ice particles.     

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.     

Seasonal Cumulative Effect of Ural Blocking Episodes on the Frequent Cold events in China during the Early Winter of 2020/21

Starting in mid-November,China was hit by several cold events during the early winter of 2020/21.The lowest temperature observed at Beijing station on 7 January reached?19.6°C.In this paper,we show that the outbreak of the record-breaking extreme cold event can be attributed to a huge merging Ural blocking(UB)ridge over the Eurasian region.The sea-ice cover in the Kara and East Siberia Seas(KESS)in autumn was at its lowest value since 1979,which could have served as a precursor signal.Further analysis shows that several successive UB episodes occurred from 1 September 2020 to 10 January 2021.The persistent UB that occurred in late September/early October 2020 may have made an important contribution to the October historical minimum of sea ice in the KESS region.Our results also show that,after each UB episode in winter,significant upward propagation of wave activity occurred around 60°E,which resulted in weakening the stratospheric vortex.Meanwhile,each UB episode also caused a significant reduction in sea-ice extent in KESS and a significant weakening of the westerly jet in mid-high-latitude Eurasia.Results suggest that the Arctic vortex,which is supposed to enhance seasonally,became weaker and more unstable than the climatic mean under the seasonal cumulative effects of UB episodes,KESS warming,and long-lasting negative-phase North Atlantic Oscillation(NAO-).Those seasonal cumulative effects,combined with the impact of La Ni?a winter,led to the frequent occurrence of extreme cold events.