On the correlation of nonlinear variables containing secular trend variations： numerical experiments

Due to global warming, the general circulation, underlying surfaces characteristics, and geophysical and meteorological elements all show evident secular trends. This paper points out that when calculating the correlation of two variables containing their own obvious secular trends, the interannual correlation characteristics between the two variables may be distorted （overestimated or underestimated）. Numerical experiments in this paper show that if two variables have opposite secular trends, the correlation coefficient between the two variables is reduced （the positive correlation is underestimated, or the negative correlation is overestimated）; and if the two variables have the same sign of secular trends, the correlation coefficient between the two variables is increased （the positive correlation is overestimated, or the negative correlation is underestimated）. Numerical experiments also suggest that the effect of secular trends on the interannual correlation of the two variables is interchangeable, that is to say, as long as the values of the two trends are not changed, the two variables interchange their positions, and the effect of the secular trends on the interannual correlation coefficient of the two variables remains the same. If the two variables have the same-（opposite-） sign trends, the effect of secular trends on the interannal correlation coefficient is more （less） distinctive. A meteorological example is given.

NUMERICAL EXPERIMENTS ON MULTI-LEVEL STATISTICAL ESTIMATION OF DYNAMIC BALANCE CONSTRAINTS IN GRAPES-3DVAR

This paper further explores the estimating and expressing of dynamic balance constraints using statistical methods in GRAPES-3DVAR(Version GM). Unlike the single-level scheme which only considers the coupling between mass and wind at one level, the multi-level scheme considers the coupling between their vertical profiles and calculates the balanced mass field at each layer using the rotational wind at all model levels. A reformed ridge regression method is used in the new scheme to avoid the multicollinearity problem and reduce the noises caused by unbalanced mesoscale disturbances. The results of numerical experiments show that the new scheme can get more reasonable vertical mass field, reduce the magnitude of the adjustment by the initialization, and improve the potential temperature analysis performance. Furthermore, the results of forecast verification in January(winter) and July(summer) both confirm that the new scheme can significantly improve the temperature forecast accuracy and bring slight positive effects to the pressure and wind forecast.

Numerical Simulation and Analysis of Storm Surges Under Different Extreme Weather Event and Typhoon Experiments in the South Yellow Sea

In this study,a coupled tide-surge-wave model was developed and applied to the South Yellow Sea.The coupled model simulated the evolution of storm surges and waves caused by extreme weather events,such as tropical cyclones,cold waves,extratropical cyclones coupled with a cold wave,and tropical cyclones coupled with a cold wave.The modeled surge level and significant wave height matched the measured data well.Simulation results of the typhoon with different intensities revealed that the radius to the maximum wind speed of a typhoon with 1.5 times wind speed decreased,and its influence range was farther away from the Jiangsu coastal region;moreover,the impact on surge levels was weakened.Thereafter,eight hypothetical typhoons based on Typhoon Chan-hom were designed to investigate the effects of varying typhoon tracks on the extreme value and spatial distribution of storm surges in the offshore area of Jiangsu Province.The typhoon along path 2 mainly affected the Rudong coast,and the topography of the Rudong coast was conducive to the increase in surge level.Therefore,the typhoon along path 2 induced the largest surge level,which reached up to 2.91 m in the radial sand ridge area.The maximum surge levels in the Haizhou Bay area and the middle straight coastline area reached up to 2.37 and 2.08 m,respectively.In terms of typhoons active in offshore areas,the radial sand ridge area was most likely to be threatened by typhoon-induced storm surges.

Preliminary Results from Numerical Experiments of a Heavy Rain Process with PENN STATE / NCAR MM5

PENN STATE/NCAR MM5 is used to simulate precipitation of the heavy rain process during 12-13 July 1994.The effects of different PBL parameterizations,resolvable scale moisture schemes and cumulus parameterization on the process rainfall simulation are investigated.By comparing the results of hydrostatic and nor.hydrostatic experiments,the nonhydrostatic impact upon precipitation is also examined.It is found in this study thai PENN STATE/NCAR MM5 has advantage not only in theory but also in simulating results with real data.In MM5.however,the selection of physical processes,especially water-cycle process,is very important and crucial to precipitation forecast of the case.It is concluded that the model with Grell (1993) scheme for convection and condensation method for resolvable scale precipitation captured the rainstorm during 12-13 July 1994 in Beijing area more successfully

NUMERICAL EXPERIMENTS FOR THE EFFECTS OF TWO MODEL INITIALIZATION SCHEMES ON RAINFALL FORECAST IN THE 2008 FLOODING SEASON

In this paper, based on heavy rain numerical forecast model AREM(Advanced Regional Eta Model), two different initialization schemes, LAPS and GRAPES-3DVAR, are used to run assimilation experiments of AREM-LAPS and AREM-3DVAR with the same data source(NCEP forecast field, surface data and radio-soundings) during the period from 21 May to 30 July 2008 to investigate the effect of the two initialization schemes on the rainfall simulation. The result suggests that:(1) the forecast TS score by the AREM-LAPS is higher than that by the AREM-3DVAR for rainfall in different areas, at different valid time and with different intensity, especially for the heavy rain, rainstorm and extremely heavy rain;(2) the AREM-3DVAR can generally simulate the average rainfall distribution, but the forecast area is smaller and rainfall intensity is weaker than the observation, while the AREM-LAPS significantly improves the forecast;(3) the AREM-LAPS gives a better forecast for the south-north shift of rainfall bands and the rainfall intensity variation than the AREM-3DVAR;(4) the AREM-LAPS can give a better reproduction for the daily change in the mean-rainfall-rate of the main rain band, and rainfall intensity changes in the eastern part of Southwest China, the coastal area in South China, the middle-lower valleys of Yangtze river, the Valleys of Huaihe river, and Shandong peninsula, with the rainfall intensity roughly close to the observation, while the rainfall intensity simulated by the AREM-3DVAR is clearly weaker than the observation, especially in the eastern part of Southwest China; and(5) the comparison verification between the AREM-LAPS and AREM-3DVAR for more than 10 typical rainfall processes in the summer of 2008 indicates that the AREM-LAPS gives a much better forecast than AREM-3DVAR in rain-band area, rainfall location and intensity, and in particular, the rainfall intensity forecast is improved obviously.

Experiments in numerical modelling of the Pacific sea surface temperature anomalies

By using the atmosphere-ocean coupled model (CGCM) which is composed of a 2-level global atmospheric general circulation model and a 4-layer Pacific oceanic general circulation model developed in the Institute of Atmospheric Physics of Chinese Academy of Sciences, and two model climatological fields got from the two independent models' numerical integrations respectively, the Pacific sea surface temperature anomalies (SSTA) from 1988 to 1989 are simulated in this paper with observed atmospheric general circulation data and sea surface temperature fields as initial conditions and monthly coupling scheme. In order to remove systematic biases of the model climatological fields, interaction variables between atmosphere and ocean are also corrected simultaneously. The experiments show that the simulation results can be improved effectively if these interaction variables are corrected in spite of the fact that there always exist systematic biases in independent numerical simulations of atmospheric part and oceanic part within CGCM. The basic characteristics of the observed Pacific SSTA in September and October 1988 have been simulated by using the correction scheme, such as the negative SSTA domain in the whole E-quatorial Pacific east to 150°E and the positive SSTA domain in the Western Pacific, the northern subtropical Pacific and nearly the whole Southern Pacific. Further numerical simulations show that the model can simulate not only the SSTA in the Pacific and its seasonal variations but also its interannual changes (for example, La Nino event in the Equatorial Pacific terminated after May 1989) to a certain degree. Furthermore, some problems existing in experiment processes and what we shoud do in the following stage are also discussed and analysed in this paper.

Some numerical experiments of the dynamics of the wintertime circulation in the Northern Huanghai Sea and the Bohai Sea

-Some numerical experiments have been developed to explain dynamically the mechanism of the formation of the wintertime circulation in the Northern Huanghai Sea and the Bohai Sea. From the results of the numerical experiments with Coriolis parameter having positive and negative values respectively, it is found that the current pattern characterized by flowing into the Bohai Sea in the northern part of the Bohai Straits and flowing out of the Bohai Sea in the southern part of the Bohai Straits is controlled basically by the bottom topography with 'the deepness on the northern side and shallowness on the southern side'in the Bohai Straits, and is not related to Coriolis parameter having positive and negative values. It also shows that the Huanghai Warm Current regarded as the inllow condition passing through the in-ports at the water boundary of the numerical experiments only influences the region near the water boundary and the eastern region of the Northern Huanghai Sea. In addition, the mechanism of the formation of several mesoscale eddies and the influence of the bottom topography on the wintertime circulation in the basin have been discussed.

Numerical experiments on the evolution in coronal magnetic configurations including a filament in response to the change in the photosphere

We investigate equilibrium height of a flux rope, and its internal equilib- rium in a realistic plasma environment by carrying out numerical simulations of the evolution of systems including a current-carrying flux rope. We find that the equilib- rium height of a flux rope is approximately described by a power-law function of the relative strength of the background field. Our simulations indicate that the flux rope can escape more easily from a weaker background field. This further confirms that a catastrophe in the magnetic configuration of interest can be triggered by a decrease in strength of the background field. Our results show that it takes some time to reach internal equilibrium depending on the initial state of the flux rope. The plasma flow inside the flux rope due to the adjustment for the internal equilibrium of the flux rope remains small and does not last very long when the initial state of the flux rope com- mences from the stable branch of the theoretical equilibrium curve. This work also confirms the influence of the initial radius of the flux rope in its evolution; the results indicate that a flux rope with a larger initial radius erupts more easily. In addition, by using a realistic plasma environment and a much higher resolution in our simulations, we notice some different characteristics compared to previous studies in Forbes.

EXPERIMENTS ON ASSIMILATION OF INITIAL VALUES IN NUMERICAL PREDICTION OF A WARM-SECTOR PRECIPITATION IN SOUTH CHINA

In order to understand the impact of initial conditions upon prediction accuracy of short-term forecast and nowcast of precipitation in South China, four experiments i.e. a control, an assimilation of conventional sounding and surface data, testing with nudging rainwater data and the assimilation of radar-derived radial wind, are respectively conducted to simulate a case of warm-sector heavy rainfall that occurred over South China, by using the GRAPES_MESO model. The results show that (1) assimilating conventional surface and sounding observations helps improve the 24-h rainfall forecast in both the area and order of magnitude; (2) nudging rainwater contributes to a significant improvement of nowcast, and (3) the assimilation of radar-derived radial winds distinctly improves the 24-h rainfall forecast in both the area and order of magnitude. These results serve as significant technical reference for the study on short-term forecast and nowcast of precipitation over South China in the future.

A Numerical Study on the Effects of Orography During the Landfall Process of Typhoon Haitang

The landfall process of typhoon Haitang which affected East China seriously was simulated by using the MM5 model and the track, intensity, precipitation and structure of typhoon were successfully reproduced. Then through the sensitive test, the effects of terrain were analyzed Results show that the irregular track during the period of typhoon passing through Taiwan and later landfalling at Fujian was in relation to the occurring and developing of orthographic impressed depression. The amount of rainfall was enhanced more than one time and the strength of typhoon was weakened 4 to 5 hPa. It is found that the effect of terrain on the structure of typhoon is limited at low level and is backward in space compared with the one at high level. In addition, the phenomenon that the equivalent temperature in the typhoon＇s moving direction inclines to the west on the eye of landfall may be concerned with the terrain.

A Numerical Study on Effects of Land-Surface Heterogeneity from' Combined Approach' on Atmospheric ProcessPart II: Coupling-Model Simulations

Two land surface schemes, one the standard Biosphere / Atmosphere Transfer Scheme Version ie (BOZ) and the other B1Z based on B0Z and heterogeneously-treated by' combined approach' , were co 'pled to the meso-scale model MM4, respectively. Through the calculations of equations from the companion paper, parameters representing land surface heterogeneity and suitable for the coupling models were found out. Three cases were simulated for heavy rainfalls during 36 hours, and the sensitivity of short-term weather modeling to the land surface heterogeneity was tested. Through the analysis of the simulations of the three heavy rainfalls, it was demonstrated that BIZ, compared with BOZ, could more realistically reflect the features of the land surface heterogeneity, therefore could more realistically reproduce the circulation and precipitation amount in the heavy rainfall processes of the three cases. This shows that even short-term weather is sensitive to the land surface heterogeneity, which is more obvious with time passing, and whose influence is more pronounced in the lower layer and gradually extends to the middle and upper layer. Through the analysis of these simulations with BlZ, it is suggested that the bulk effect of smaller-scale fluxes (i.e., the momentum, water vapor and sensible heat fluxes) near the s ig nificantly-heterogeneous land surface is to change the larger-scale (i.e., meso-scale) circulation, and then to influence the development of the low-level jets and precipitation. And also, the complexity of the land-atmosphere interaction was shown in these simulations.

A Numerical Experiment Study for Effects of the GrasslandDesertification on Summer Drought in North China

In this paper, the summer climate of 1991 in North China is simulated by using the high-resolution regional climate model (RegCM2) and the effects of the grassland desertification on summer drought in the central and the northern parts of North China as well as Mongolia are studied. It shows that the regional climate model essentially catches the characteristics on distribution and seasonal variation of the precipitation that keep good agreement with the observation. The desertification makes precipitation in the central part of North China during its flood period decrease obviously in July. The border of the precipitation or the soil moisture reduction in the desertification region extends about one latitude southeastward and beyond the southeast edge of the desertification. Thus, vegetation in the border region approaches desertification further. However, there appears evident difference of variation of precipitation over the whole desertification region. The grassland desertification greatly changes the transfers of fluxes between land and atmosphere. The secondary circulation or secondary circulation cells in the desertification region are excited and as a result moisture transport is changed. The variation of flux transfers between land and atmosphere as well as the vertical motion of atmosphere is closely related to that of precipitation.

A Numerical Study of the Interaction between Two Simultaneous Storms:Goni and Morakot in September 2009

Significant anomalous tracks were observed when the severe tropical storm Goni （0907） and typhoon Morakot （0908） in September 2009 were evaluated in short distances. The relationship between the two is regarded as a case of binary interaction. Based on an MM5 model （fifth=generation mesoscale model of Pennsylvania State University-National Center for Atmospheric Research）, in this study a series of sensitivity experiments were designed to determine the binary interaction between them. The sensitivity of the storm characteristics to the binary interaction was demonstrated through modeling experiments with different TC intensities and sizes based on the bogus vortices initialization. Furthermore, the contributions of large-scale environmental flow and the effects of interaction between the motions of the cyclones were distinguished by numerical experiments using only one of the TC vortices. Results from these experiments show that Morakot （0908） had a greater impact on the motion of Goni （0907）, whereas Goni （0907） had a relatively limited impact on Morakot （0908）. At the upper level, the northeasterly jet flow in the third quadrant of Morakot （0908） enhanced the upper-level divergence of Goni （0907） and had an important role in maintaining and increasing Goni＇s （0907） intensity. And at the lower level, Morakot （0908）, with strong convergence and ascending airflow, made a stable transport channel of southwesterly warm and wet flow, thus supporting the lower-level water vapor convergence of Goni （0907）. Goni （0907）, which was located upriver of the southwesterly flow, intercepted part of the water vapor transportation in the southwesterly flow, causing the water vapor convergence to strengthen while the water vapor convergence of Morakot （0908） weakened.

NUMERICAL SIMULATION OF THE RELATIONSHIP BETWEEN THE MAINTENANCE AND INCREASE IN HEAVY RAINFALL OF THE LANDING TROPICAL STORM BILIS AND MOISTURE TRANSPORT FROM LOWER LATITUDES

The NCEP/NCAR reanalysis, Japan Meteorological Agency(JMA) tropical cyclone tracks and intensive surface observations are used to diagnose the features of moisture transport of tropical storm Bilis(No. 0604), which is simulated by the WRF(weather research and forecasting) mesoscale numerical model. It is shown that the Bilis was linked with the moisture channel in the lower latitudes after its landing. Meanwhile, the cross-equatorial flows over 80°-100°E and Somali were active and brought abundant water vapor into the tropical storm, facilitating the maintenance of the landing storm with intensified heavy rainfall along its path. The simulation suggested that the decreased water vapor from lower latitudes prevents the maintenance of Bilis and the development of rainfall. While the cross-equatorial flows over 80°-100°E and Somali were in favor of keeping the cyclonic circulation over land. If the moisture supply fro m the Somali jet stream was reduced, the strength and area of heavy rainfall in tropical cyclone would be remarkably weakened. Consequently, the decreased water vapor from lower latitudes can remarkably suppress the deep convection in tropical storm, then Bilis was damped without the persistent energy support and the rainfall was diminished accordingly.

A Numerical Study on Effects of Land-Surface Heterogeneity from "Combined Approach" on Atmospheric Process Part I: Principle and Method

A method based on Giorgi (1997a, 1997b) and referred to as ’ combined approach’, which is a combi-nation of mosaic approach and analytical-statistical-dynamical approach, is proposed. Compared with those of other approaches, the main advantage of the combined approach is that it not only can represent both interpatch and intrapatch variability, but also cost less computational time when the land surface heterogeneity is considered. Because the independent variable of probability density function (PDF) is ex-tended to the single valued function of basic meteorological characteristic quantities, which is much more universal, the analytical expressions of the characteristic quantities (e.g., drag coefficient, snow coverage, leaf surface aerodynamical resistance) affected by roughness length are derived , when the roughness length(and / or the zero plane displacement) heterogeneity has been mainly taken into account with the approach. On the basis of the rule which the PDF parameters should follow, we choose a function y of the roughness length z 0 as the PDF independent variable, and set different values of the two parameters width ratio αn and height ratio γ of PDF (here a linear, symmetric PDF is applied) for sensitivity experiments, from which some conclusions can be drawn, e.g., relevant characteristic terms show different sensitivities to the heterogeneous characteristic (i.e., roughness length), which suggests that we should consider the heterogeneities of the more sensitive terms in our model instead of the heterogeneities of the rest, and which also implies that when the land surface scheme is coupled into the global or regional atmospheric model, sensitivity tests against the distribution of the heterogeneous characteristic are very necessary; when the parameter αn is close to zero, little heterogeneity is represented, and αn differs with cases, which have an upper limit of about 0.6; in the reasonable range of αn, a peak-like distribution of roughness length can be depicted by a small value of γ, etc.. Key words Representation of land surface heterogeneity - “ Combined approach” - Numerical experiment This work was supported by the National Sciences Foundation of China, Grant No.49875005 and the State Key Project (973) G19990434 (03).

ENSEMBLE PREDICTION EXPERIMENTS OF TRACKS OF TROPICAL CYCLONES BY USING MULTIPLE CUMULUS PARAMETERIZATION SCHEMES

Ensemble prediction experiments of the tracks of eight tropical cyclones occurring between 2004-2006 over the western Pacific have been performed by using MM5 with five cumulus parameterization schemes. The results show that the predictions of the tracks of the tropical cyclones are sensitive to the selection of cumulus parameterization schemes. Each scheme has its own advantage and disadvantage, and the predications without cumulus parameterization schemes are not the worst, sometimes even better than the others. And all of the three ensemble methods improve the predictions of the tracks significantly, among which the ensemble method without parameterization schemes, the Grell, Betts-Miller and Kain-Fritsch schemes are the best.

Numerical experiment rock fragmentation by combined dynamic and static loads under dual-cutter head

This paper puts forward a new rock fragmentation loading method of dual-cutter head combined dynamic and static loads. By applying the numerical simulation software - RFPA2D, we have done numerical experiment about the sihstone＇ s crushing effect by dynamic load on single cutter head without confining pressure, dynamic load on single cut- ter head with confining pressure 10 MPa and different dual-cutter heads spacing by combined dynamic and static loads with confining pressure 10 MPa. Experimental results show that the confining pressure can obviously affect the rock frag- mentation effect. Combined dynamic and static loads can greatly improve the rock fragmentation effect. There exists an optimal spacing of dual-cutter head that can make the rock fragmentation achieve the desired effect. Through analyzing the acoustic emission accumulative energy and quantity, the authors make a conclusion that the optimum spacing is 30 mm.

NUMERICAL STUDY ON THE FORMATION OF THE SOUTH CHINA SEA WARM CURRENT II. BAROCLINIC CASE

In this part, Levitus’ climatological temperature and salinity are incorporated in the numerical model developed in Part I. Diagnostic and prognostic experiment on the thermohaline circulation were conducted. The smooth Levitus’ data do not include any information on the South China Sea Warm Current (SCSWC), so it is not in the model produced diagnostic thermohaline circulation. Although the SCSWC does not appear in the wind driven circulation in the barotropic case, it appears in the prognostic wind driven circulation in the baroclinic case. This implies that the differing circulation patterns between barotropic case and baroclinic case are due to the stratification. The prognostic thermohaline circulation with wind stress and inflow/outflow transports at open boundaries are also discussed. Coupling of density and dynamic forces makes the circulation pattern more complicated. Even though the stratification is not always a direct cause of the formation of the SCSWC, it is at least an indirect cause.

NUMERICAL ANALYSIS ON THE ROLE OF TROPICAL STORM NAMTHEUN IN THE UNUSUAL TRACKS OF THREE TROPICAL CYCLONES IN 2010

Tropical cyclones(TCs) Lionrock,Kompasu,and Namtheun were formed successively within 40 hours in 2010.Over the next several days afterwards,these TCs exhibited unusual movements which made operational prediction difficult.Verifications are performed on the forecasts of the tracks of these TCs with six operational models,including three global and three regional models.Results showed that the trends of TC tracks could be reproduced by these models,whereas trajectory turning points and landfall locations were not simulated effectively.The special track of Lionrock should be associated with its direct interaction with Namtheun,according to a conceptual model of binary TC interaction.By contrast,the relation between Kompasu and Namtheun satisfied the criteria for a semi-direct interaction.Numerical experiments based on the Global and Regional Assimilation and Prediction System-Tropical Cyclone forecast Model(GRAPES-TCM) further confirmed the effect of Namtheun on the unusual tracks of Lionrock and Kompasu.Finally,the physical mechanism of binary TC interaction was preliminarily proposed.

A Preliminary Study on 1D Numerical Experiment of Water Debris Flow in Gully

In order to improve and enhance the numerical modeling methods and its application on debris flow problems,a preliminary study has been proposed in accordance with the corrected water-sediment numerical model on the premise of revised resistance and sediment capacity equations.Firstly,an overview the recent re- search achievements on numerical simulation of debris flow has been conducted,the results shown that a gener- al numerical model for debris flow can not be existed at all because the complex rheol...