显式自由表面模型时间分裂格式的一个积分方案

讨论了有关显式自由表面模型的时间分裂格式中的分裂误差。在传统的积分方案中，由于采用蛙跳格式，斜压动量方程具有三层的时间离散结构，但是，正压动量方程在相应的时间层上只有二层的离散结构。因此，正压方程与斜压方程的时间离散结构是不相容的。分析结果表明，这是导致分裂误差的原因。本提出一种能保持时间离散结构相容的正压和斜压方程的时间积分方案。该积分方案的分裂误差很小因此特别适合于需要长期积分的气候研究。该方案在北太平洋的一个自由表面模型中得到了应用。     

Ekman动量近似下中间边界层模式中的风场结构

发展了一个准三维的、中等复杂的边界层动力学模式，该模式包含了EKman动量近似下的惯性加速度和Blackadar的非线性湍流粘性系数，它进一步改进了Tan和Wu(1993）提出的边界层理论模型。该模型在数值计算复杂性上与经典Ekman模式相类似，但由于包含了Ekman动量近似下的惯性项，使得该模式比传统Ekman模式更近于实际过程。中详细地比较了该模式与其他简化边界层模式在动力学上的差异，结果表明：在经典的Ekman模式中，由于忽略了流动的惯性项作用，导致在气旋性切变气流（反气旋性切变气流）中风速和边界层顶部的垂直速度的高估（低估），而在半地转边界层模式中，由于高估了流动惯性项的作用，结果与经典Ekman模式相反。同样，该模式可以应用于斜压边界层，对于Ekman动量下的斜压边界层风场同时具有经典斜压边界层和Ekman动量近似边界层的特征。     

改进的切线性模式对一个边界层模式变分资料同化的改善

大气、海洋中边界层物理过程的强非线性会对切线性模式近似（尤其是当积分时间长时）有较大的影响，从而给相应的4维变分资料同化问题的求解造成困难。本构造了Mellor-Yamada湍流封闭模式（level 2.5)的一个改进的切线性模式，相比通常的切线性模式和简化的切线性模式可以提高对非线性扰动的逼近。利用这个改进的切线性模式的伴随模式分别进行了1到7天的变分资料同化试验，得到了满意的结果。而用通常的伴随模式和简化的伴随模式都无法得到令人满意的结果。     

Observations Of Nocturnal Drainage Flow In A Shallow Gully

The formation of cold air drainage flows in a shallow gully is studied during CASES-99 (Cooperative Atmosphere-Surface Exchange Study). Fast and slow response wind and temperature measurements were obtained on an instrumented 10-m tower located in the gully and from a network of thermistors and two-dimensional sonic anemometers, situated across the gully. Gully flow formed on clear nights even with significant synoptic flow. Large variations in surface temperature developed within an hour after sunset and in situ cooling was the dominant factor in wind sheltered locations. The depth of the drainage flow and the height of the down-gully wind speed maximum were found to be largest when the external wind speed above the gully flow is less than 2 m s^-1. The shallow drainage current is restricted to a depth of a few metres, and is deepest when the stratification is stronger and the external flow is weaker. During the night the drainage flow breaks down, sometimes on several occasions, due to intermittent turbulence and downward fluxes of heat and momentum. The near surface temperature may increase by 6 ° C in less than 30 min due to the vertical convergence of downward heat flux. The mixing events are related to acceleration of the flow above the gully flow and decreased Richardson number. These warming events also lead to warming of the near surface soil and reduction of the upward soil heat flux. To examine the relative importance of different physical mechanisms that could contribute to the rapid warming, and to characterize the turbulence generated during the intermittent turbulent periods, the sensible heat budget is analyzed and the behaviour of different turbulent parameters is discussed.     

Atmospheric Surface-Layer Processes During the Total Solar Eclipse Of 11 August 1999

The central path of the total solar eclipse (TSE) of 11 August 1999 mostly passedthrough a region of active monsoon in India, with the eclipse ending around localsunset. Measurements in the surface layer (SL) were carried out close to the eclipseaxis at Akola (20°42' N, 77°2' E) in central India. The globalradiation flux reaching the surface vanishes around totality at 1803:24 (LT), followedby a small recovery before again dipping to zero at sunset. The temperatures in the SL, and subsoil at 50-mm depth, show a local minimum with a lag of about 10 min after the second contact, whereas the lag appears to vanish when the temperature series is detrended. The SL exhibits near-neutral, though generally stable, conditions from about 1500 hr itself due to heavy cloud cover followed by the eclipse-induced cooling of the surface. The wind component across the eclipse axis vanishes at totality, the wind vector aligning with the azimuth of the traversing moon shadow. The deceleration of the mean flow can be due to the combined effect of the colder surface and downward heat flux, where the locally altered horizontal temperature gradients may cause, as in this instance, the cessation of the cross flow.The oscillations in temperature and wind that show significant peaks, around 90–100 min as well as the semi-period of the eclipse near 60 min, persist for several hours past the eclipse event. A fresh analysis of the published data on the TSE of 7 March 1970 shows spectral peaks in the temperature nearly coincident with those already reported from the surface pressure records. The oscillations in the SL fields during the two TSE events are very similar implying that the source mechanism is also likely to be the same in both the cases.     

The Production of Methanol by Flowering Plants and the Global Cycle of Methanol

Methanol has been recognised as an important constituent of the background atmosphere, but little is known about its overall cycle in the biosphere/atmosphere system. A model is proposed for the production and emission to the atmosphere of methanol by flowering plants based on plant structure and metabolic properties, particularly the demethylation of pectin in the primary cell walls. This model provides a framework to extend seven sets of measurements of methanol emission rates to the global terrestrial biosphere. A global rate of release of methanol from plants to the atmosphere of 100 Tg y^–1 is calculated. A separate model of the global cycle of methanol is constructed involving emissions from plant growth and decay, atmospheric and oceanic chemical production, biomass burning and industrial production. Removal processes occur through hydroxyl radical attack in the atmosphere, in clouds and oceans, and wet and dry deposition. The model successfully reproduces the methanol concentrations in the continental boundary-layer and the free atmosphere, including the inter-hemispheric gradient in the free atmosphere. The model demonstrates a new concept in global biogeochemistry, the coupling of plant cell growth with the global atmospheric concentration of methanol. The model indicates that the ocean provides a storage reservoir capable of holding at least 66 times more methanol than the atmosphere. The ocean surface layer reservoir essentially buffers the atmospheric concentration of methanol, providing a physically based smoothing mechanism with a time constant of the order of one year.     

Aerosol climatology for the planetary boundary layer derived from regular lidar measurements

Regular aerosol extinction and backscatter measurements using a UV Raman Lidar have been performed for almost 3 years in Hamburg in the frame of the German Lidar Network. A set of 92 aerosol extinction and 164 aerosol backscatter profiles has been used for statistical investigations. Mean values and variances of the aerosol extinction and backscatter in the boundary layer have been calculated. Large fluctuations during the whole year have been found. The measured aerosol extinction over Hamburg shows a seasonal cycle with highest values in early fall and a second less prominent peak in spring.An analysis of the data using back trajectories showed a dependence of the aerosol extinction on the origin of the air mass. The residence time of the air mass over industrialized areas was found to be an important parameter for the measured aerosol extinction at Hamburg. However, only a small part of the total variability could be explained by the air mass origin.For 75 cases of aerosol extinction measurements under cloud-free conditions, the aerosol backscatter profile and therefore, the lidar ratio as a function of altitude could be determined. Winter measurements of the lidar ratio are often close to model results for maritime aerosol, the summer measurements are close to the model results for urban or continental aerosols.The high quality of the data has been proven by intercomparisons with other lidar systems and with star photometer measurements of the aerosol optical depth during the Lindenberg Aerosol Characterization Experiment (LACE'98) field campaign.     

Formation of a shelfbreak front by freshwater discharge

The circulation and transport of freshwater generated by an idealized buoyant source is studied using a three-dimensional primitive equation model. Freshwater enters the continental shelf, turns anticyclonically and moves downstream in the direction of Kelvin wave propagation. In the region close to the source, the flow reaches an equilibrium in the bottom boundary layer so that freshwater does not spread offshore any further. This offshore equilibrium distance increases as we move downstream until the freshwater is able to feel the presence of the shelfbreak. A shelfbreak front forms and the shelfbreak prevents any further offshore spreading of freshwater in the bottom boundary layer.Two complimentary mechanisms are responsible for the slow cross-shelf migration of freshwater and subsequent trapping of shelfbreak fronts: bottom stress and topographic changes. The shelfbreak creates an active, dynamic process preventing leakage from the continental shelf region to the slope region. However, the dynamical process that traps the front to the shelfbreak is still unclear.The location of the shelfbreak front depends on four dimensionless parameters: scaled inlet volume transport, scaled breadth, scaled “diffusivity” and scaled shelf width. We develop empirical relations for predicting the location of the frontal bottom intersection, given these parameters.     

可控深度电法对水平柱体和垂直界面的计算和实验

球状地球体可控深度电法的效应系数可近似地全面表达可控深度电法和各种效应。通过点源二度体的电算模拟实验，证明了该方法在水平柱状体和地下存在电阻率垂直界面的条件下，也是有效的。     

静压预制桩工程性状研究──以河南某电厂技改工程为例

以河南某电厂技改工程为例 ,根据现场检测结果 ,探讨了影响静压预制桩沉桩深度成因 ,分析了静压桩承载力的分布规律。研究结果证明静压桩的沉桩深度和承载力与地基土的性状密切相关。对于河南某电厂技改工程场地 ,尽管桩端进入砂较小深度 ,却可以在桩端平面附近良好土层的作用下提供较大的端阻力。