一种新的统计回归模型及其建模方案

根据经典门限自回归模型的基本思想,引入半截多项式变换,导出了一种新的统计回归模型,并提出了相应的一整套建模方案。这种模型和建模方案的特点是:（1）解决了统计回归分析中逐段线性化模型的检验问题;（2）在确定统计回归模型中各待估参数和变量——包括门限变量、门限值、阶数、时滞和回归系数时显得十分方便、快捷。试验结果表明,根据此方案建立的统计回归模型具有较高的拟合和预报精度,同时具有良好的稳定性。     

Some Insights into Topographic,Elastic and Self-gravitation Interaction in Modelling Ground Deformation and Gravity Changes in Active Volcanic Areas

Surface displacements and gravity changes due to volcanic sources are influenced by medium properties. We investigate topographic, elastic and self-gravitation interaction in order to outline the major factors that are significant in data modelling. While elastic-gravitational models can provide a suitable approximation to problems of volcanic loading in areas where topographic relief is negligible, for prominent volcanoes the rough topography could affect deformation and gravity changes to a greater extent than self-gravitation. This fact requires the selection, depending on local relief, of a suitable model for use in the interpretation of surface precursors of volcanic activity. We use the three-dimensional Indirect Boundary Element Method to examine the effects of topography on deformation and gravity changes in models of magma chamber inflation/deflation. Topography has a significant effect on predicted surface deformation and gravity changes. Both the magnitude and pattern of the geodetic signals are significantly different compared to half-space solutions. Thus, failure to account for topographic effects in areas of prominent relief can bias the estimate of volcanic source parameters, since the magnitude and pattern of deformation and gravity changes depend on such effects.     

Absorption properties of urban/suburban aerosols in China

The broadband diffuse radiation method is improved to retrieve the aerosol refractive index imaginary part （AIP） and broadband （400-1000 nm mean） single scattering albedo （SSA）. In this method, four sets of SSA selection criteria are proposed for quality control. The method is used to retrieve AIP, SSA and absorptive optical thickness （AbOT） from routine hourly-exposed pyrheliometer and paranometer measurements over 11 sites （meteorological observatories） in China during 1998-2003. Apart from one suburban site （Ejin Qi）, the other urban sites are all located around big or medium cities. As shown in the retrieval results, annual mean SSA during 1998-2003 changes from 0.941 （Wuhan） to 0.849 （Lanzhou）, and AIP from 0.0054 to 0.0203. The 11-site average annual mean SSA and AIP are 0.898 and 0.0119, respectively. SSA during winter is smaller for most sites. There is an evidently positive correlation between SSA and aerosol optical thickness （AOT） for all sites. There is also a positive correlation between SSA and relative humidity for most sites, but a negative correlation for a few sites, such as Kashi and ǚrǚmqi in Northwest China.     

Climatic change and Chinese population growth dynamics over the last millennium

The climate–population relationship has long been conceived. Although the topic has been repeatedly investigated, most of the related works are Eurocentric or qualitative. Consequently, the relationship between climate and population remains ambiguous. In this study, fine-grained temperature reconstructions and historical population data sets have been employed to statistically test a hypothesized relationship between temperature change and population growth (i.e., cooling associated with below average population growth) in China over the past millennium. The important results were: (1) Long-term temperature change significantly determined the population growth dynamics of China. However, spatial variation existed, whilst population growth in Central China was shown to be responsive to both long- and short-term temperature changes; in marginal areas, population growth was only sensitive to short-term temperature fluctuations. (2) Temporally, the temperature–population relationship was obscured in some periods, which was attributable to the factors of drought and social buffers. In summary, a temperature–population relationship was mediated by geographic factors, the aridity threshold, and social factors. Given the upcoming threat posed by climate change to human societies, this study seeks to improve our knowledge and understanding of the climate–society relationship.     

Scaling of Wall-Normal Turbulence Intensity and Vertical Eddy Structures in the Atmospheric Surface Layer

Adequate high-quality data on three-dimensional velocities in the atmospheric surface layer (height \(\delta \)) were acquired in the field at the Qingtu Lake Observation Array. The measurement range occupies nearly the entire logarithmic layer from approximately \(0.006\delta \)–\(0.2\delta \). The turbulence intensity and eddy structures of the velocity fluctuations in the logarithmic region were primarily analyzed, and their variations in the z (wall-normal) direction were revealed. The primary finding was that the turbulent intensity of wall-normal velocity fluctuations exhibits a sharp upswing in the logarithmic region, which differs from classic scaling law and laboratory results. The upswing of the wall-normal turbulence intensity in the logarithmic region is deemed to be linear based on an ensemble of 20 sets of data. In addition, the wall-normal extent of the correlated structures and wall-normal spectra were compared to low Reynolds number results in the laboratory.     

UV attenuation in the cloudy atmosphere

Ultraviolet (UV) energy absorption plays a very important role in the Earth–atmosphere system. Based on observational data for Beijing, we suggest that some atmospheric constituents utilize or transfer UV energy in chemical and photochemical (C&P) reactions, in addition to those which absorb UV energy directly. These constituents are primarily volatile organic compounds (VOCs) emitted from both vegetative and anthropogenic sources. The total UV energy loss in the cloudy atmosphere for Beijing in 1990 was 78.9 Wm⁻². This attenuation was caused by ozone (48.3 Wm⁻²), other compounds in the atmosphere (26.6 Wm⁻²) and a scattering factor (4.0 Wm⁻²). Our results for a cloudy atmosphere in the Beijing area show that the absorption due to these other compounds occurs largely through the mediation of water vapor. This fraction of energy loss has not been fully accounted for in previous models. Observations and previous models results suggest that 1) a cloudy atmosphere absorbs 25∼30 Wm⁻² more solar shortwave radiation than models predict; and 2) aerosols can significantly decrease the downward mean UV-visible radiation and the absorbed solar radiation at the surface by up to 28 and 23 Wm⁻², respectively. Thus, quantitative study of UV and visible absorption by atmospheric constituents involved in homogeneous and heterogeneous C&P reactions is important for atmospheric models.     

Estimation of Surface-Layer Scaling Parameters in the Unstable Boundary Layer: Implications for Orographic Flow Speed-Up

A method is proposed for estimating the surface-layer depth \((z_s)\) and the friction velocity \((u_*)\) as a function of stability (here quantified by the Obukhov length, L) over the complete range of unstable flow regimes. This method extends that developed previously for stable conditions by Argaín et al. (Boundary-Layer Meteorol 130:15–28, 2009), but uses a qualitatively different approach. The method is specifically used to calculate the fractional speed-up \((\varDelta S)\) in flow over a ridge, although it is suitable for more general boundary-layer applications. The behaviour of \(z_s \left( L\right) \) and \(u_*\left( L\right) \) as a function of L is indirectly assessed via calculation of \(\varDelta S\left( L\right) \) using the linear model of Hunt et al. (Q J R Meteorol Soc 29:16–26, 1988) and its comparison with the field measurements reported in Coppin et al. (Boundary-Layer Meteorol 69:173–199, 1994) and with numerical simulations carried out using a non-linear numerical model, FLEX. The behaviour of \(\varDelta S\) estimated from the linear model is clearly improved when \(u_*\) is calculated using the method proposed here, confirming the importance of accounting for the dependences of \(z_s\left( L \right) \) and \(u_*\left( L \right) \) on L to better represent processes in the unstable boundary layer.     

Characteristics of carbonyl compounds in ambient air of Shanghai,China

The levels of carbonyl compounds in Shanghai ambient air were measured in five periods from January 2007 to October 2007 (covering winter, high-air-pollution days, spring, summer and autumn). A total of 114 samples were collected and eighteen carbonyls were identified. Formaldehyde, acetaldehyde and acetone were the most abundant carbonyls and their mean concentrations of 19.40 ± 12.00, 15.92 ± 12.07 and 11.86 ± 7.04 μg m⁻³ respectively, in the daytime for five sampling periods. Formaldehyde and acetaldehyde showed similar diurnal profiles with peak mixing ratios in the morning and early afternoon during the daytime. Their mean concentrations were highest in summer and lowest in winter. Acetone showed reversed seasonal variation. The high molecular weight (HMW, ≥C5) carbonyls also showed obvious diurnal variations with higher concentrations in the daytime in summer and autumn, while they were all not detected in winter. Formaldehyde and acetaldehyde played an important role in removing OH radicals in the atmosphere, but the contribution of acetone was below 1%. The carbonyls levels in high-air-pollution days were reported. More carbonyl species with higher concentrations were found in high-air-pollution days than in spring. These carbonyls were transported with other pollutants from north and northwest in March 27 to April 2, 2007 and then mixed with local sources. Comparing with Beijing and Guangzhou, the concentrations of formaldehyde and acetaldehyde in Shanghai were the highest, which indicated that the air pollution in Shanghai was even worse than expected.     

Adaptation of Pressure Based CFD Solvers for Mesoscale Atmospheric Problems

General purpose Computational Fluid Dynamics (CFD) solvers are frequently used in small-scale urban pollution dispersion simulations without a large extent of ver- tical flow. Vertical flow, however, plays an important role in the formation of local breezes, such as urban heat island induced breezes that have great significance in the ventilation of large cities. The effects of atmospheric stratification, anelasticity and Coriolis force must be taken into account in such simulations. We introduce a general method for adapting pressure based CFD solvers to atmospheric flow simulations in order to take advantage of their high flexibility in geometrical modelling and meshing. Compressibility and thermal stratification effects are taken into account by utilizing a novel system of transformations of the field variables and by adding consequential source terms to the model equations of incompressible flow. Phenomena involving mesoscale to microscale coupled effects can be analyzed without model nesting, applying only local grid refinement of an arbitrary level. Elements of the method are validated against an analytical solution, results of a reference calculation, and a laboratory scale urban heat island circulation experiment. The new approach can be applied with benefits to several areas of application. Inclusion of the moisture transport phenomena and the surface energy balance are important further steps towards the practical application of the method.     

CHARACTERISTICS OF LARGE PULSES OF RADIATED FIELDS FROM CLOUD FLASH DISCHARGES

Observations show that the sequences of microsecond-scale pulses are produced in cloud flashdischarge processes.In the present work we report on measurements of larger bipolar pulses(LBP)of the radiation fields which are classified into four types according to their waveshapes andpolarities of initial peak.Either positive or negative pulse contains two kinds of waveshapes:onewith a smooth rise and the other with 2 or 3 fast-rising pulse superimposed on the initial half cycle.The pulse waveshape characteristics in both time and frequency domains are as follows.Among2908 samples there are 45% of the positive polarity and 55% of the negative polarity.The positivepulses tend to have distinctly smaller values than negative pulses in rise time,half width and zero-crossing time of the initial peak.The peaks of positive pulse spectra tend to occur at a higherfrequency than those of negative pulses.Despite differences above,positive and negative pulseshave a similar total duration.