Turbulence flux measurement above the overstory of a subtropical Pinus plantation over the hilly region in southeastern China

Continuous turbulence flux measurement using the eddy covariance (EC) technique was made from January 1 to December 31 in 2003 at two and three canopy heights of a subtropical Pinus plantation on the red earth hilly region in southeastern China. To be able to make sure that the measured turbulence flux will equal the net ecosystem/atmosphere exchange, the quality of the data has to be assessed. Three criteria were investigated here, including the power spectra and cospectra analyses, flux variance similarity (integral turbulence test) and energy balance closure. The spectral analyses suggested that above-canopy power spectral slopes for all velocity components and scalars such as CO2, H2O and air temperature followed the expected -2/3 power law in the inertial subrange, and their cospectral slopes were close to -4/3 power law in the inertial subrange. The important contribution of large-scale motions to energy and mass transfer above the canopy at higher measurement level was also confirmed by the spectral analyses. The eddy covariance systems have the ability to resolve fluctuations associated with small-scale eddies and did not induce an obvious underestimation of the measured turbulence flux. The Monin-Obukhov similarity functions for the normalized standard deviation of vertical wind speed and air temperature were well-defined functions of atmospheric stability at two heights above the forest canopy, which indicated that turbulence flux measurements made at two heights were within the surface layer. Nocturnal flux underestimation and departures of this normalized standard deviation of vertical wind speed similarity function from that expected from Monin-Obukhov theory were a function of friction velocity. Thus, an optimal criterion of friction velocity was determined to be greater than 0.2-0.3 m s-1 for nocturnal fluxes so that the eddy covariance flux measurement was under high turbulent mixing conditions. Energy balance closure reached about 72%-81% at the studied site, which was comparable to the 10%-30% of energy imbalance reported in the literature. However, the energy balance closure could be only used as a useful reference criterion.     

Resistive instabilities in rapidly rotating fluids: Linear theory of the convective modes

Abstract

This paper analyzes the linear stability of a rapidly-rotating, stratified sheet pinch in a gravitational field, g, perpendicular to the sheet. The sheet pinch is a layer (O ? z ? d) of inviscid, Boussinesq fluid of electrical conductivity σ, magnetic permeability μ, and almost uniform density ρ _o; z is height. The prevailing magnetic field. B _o(z), is horizontal at each z level, but varies in direction with z. The angular velocity, Ω, is vertical and large (Ω ? V_A/d, where V_A = B₀√(μρ₀) is the Alfvén velocity). The Elsasser number, Λ = σB² ₀/2Ωρ₀, measures σ. A (modified) Rayleigh number, R = gβd²/ρ₀V² _A, measures the buoyancy force, where β is the imposed density gradient, antiparallel to g. A Prandtl number, P_K = μσK, measures the diffusivity, k, of density differences.     

On conditions for the unique inversion of seismic travel-time curves

It is shown that when the travel-time curve of a refracted wave from a surface source is known and at least one of the following conditions is satisfied, i.e. when

1. the travel-time curve of a wave reflected from a horizontal interface lying below the deepest low velocity layer is known, or
2. the travel-time curve of a wave from a deep source situated below the deepest low velocity layer is known, or
3. the measureH(u)=mes {z∶z≥0,v ^?1 (z)≥u} is analytical in some segment [c, d], where \(0< c< d< \infty , c< a_n , H(a_n ) = \bar z_n ,\bar z_n\) is the depth of the lower end of the deepest low velocity layer and in the interval [c, ∞) an analytical functionH(u)) exists which providesH(u)≡H(u)) ifu∈[c, d], then (1) velocityv(z) outside the low velocity layers and (2) the measureH _k (u)=mes {z∶z∈L _k,v ^?1 (z)≥u} for each low velocity layerL _k,k=1, 2, ..., n, are defined unambiguously.

     

Observations of bedforms on a dissipative macrotidal beach

Field measurements of wave ripples and megaripples were made with a Sand Ripple Profiler in the surf and shoaling zones of a sandy macrotidal dissipative beach at Perranporth, UK in depths 1–6 m and significant wave heights up to 2.2 m. A frequency domain partitioning approach allowed quantification of height (η), length (λ) and migration rate of ripples and megaripples. Wave ripples with heights up to 2 cm and wavelengths ～20 cm developed in low orbital velocity conditions (u _m?<?0.65 m/s) with mobility number ψ?<?25. Wave ripple heights decreased with increasing orbital velocity and were flattened when mean currents were >0.1 m/s. Wave ripples were superimposed on top of megaripples (η?=?10 cm, λ?=?1 m) and contributed up to 35 % of the total bed roughness. Large megaripples with heights up to 30 cm and lengths 1–1.8 m developed when the orbital velocity was 0.5–0.8 m/s, corresponding to mobility numbers 25–50. Megaripple heights and wavelengths increased with orbital velocity but reduced when mean current strengths were >0.15 m/s. Wave ripple and megaripple migrations were generally onshore directed in the shoaling and surf zones. Onshore ripple migration rates increased with onshore-directed (+ve) incident wave skewness. The onshore migration rate reduced as offshore-directed mean flows (undertow) increased in strength and reached zero when the offshore-directed mean flow was >0.15 m/s. The migration pattern was therefore linked to cross-shore position relative to the surf zone, controlled by competition between onshore-directed velocity skewness and offshore-directed mean flow.     

岩盐断层带摩擦滑动的速度依赖性转换及其地震学意义

为了深入理解断层带摩擦滑动速度依赖性转换及其机制,利用双轴摩擦实验对干燥及含水条件下岩盐断层带摩擦的速度依赖性进行了实验研究,并观测了摩擦滑动过程中的声发射,分析了断层带的微观结构.实验结果表明,干燥岩盐断层带在0.1~100 μm/s的速度范围内表现为速度弱化,增大σ₂会使断层带向速度强化转变;含水条件下岩盐断层带在1~100 μm/s的速度范围内表现为速度弱化,而在0.1~0.01 μm/s的速度范围内表现为速度强化,速度依赖性转换出现在0.1~1 μm/s,其中断层表现为振荡或应力释放时间较长的黏滑事件;岩盐断层带在干燥条件下表现出很强的声发射活动,每个黏滑均对应一丛声发射事件,而在含水条件下一次黏滑只对应一个声发射事件.显微观察表明,局部化的脆性破裂是速度弱化域的主要变形机制,分布式的碎裂流动是干燥岩盐断层带在速度强化域的变形机制,颗粒边界迁移以及压溶作用的塑性变形是含水条件下岩盐断层带在速度强化域的主要变形机制,而脆性破裂和塑性变形共同控制着速度依赖性转换域断层带的变形.水的存在促进岩盐发生塑性变形,进而导致断层带从速度弱化向速度强化转换.上述结果有助于理解断层带上地震活动的特征和慢地震的机制.     

戈壁地区夏季碳和水热循环及其湍流特征

本文利用2006年夏季大气边界层观测资料,采用涡动相关法分析了我国西北地区戈壁下垫面碳收支及水热循环的规律和特征,并分析了大气湍流特征.结果表明:夏季白天的CO₂湍流通量呈逆输送特征,即CO₂白天向下输送,夜间向上输送,平均数值为-0.199 mg·m^-2·s^-1,整体上表现为碳汇;戈壁地区湿度小,其数值受水平来流的影响较大,日变化特征不明显;温度的归一化标准差与稳定度参数的关系满足Monin-Obukhov Similarity (MOS)理论;温度和CO₂的能谱相似;互谱uc与uθ,wc与wθ相似;水汽和CO₂的输送主要受水平方向湍流的影响.     

A model for the simulation of turbulent and eddy diffusion processes at heights of 0–65 km

A generalized turbulent diffusion model has been developed which evaluates the time rate of growth of a simulated cloud of particles released into a turbulent (i.e. diffusive) atmosphere. The general model, in the form of second-order differential equations, computes the three-dimensional size of the cloud as a function of time. Parameters which influence the cloud growth, and which are accounted for in the model equations, are: (1) length scales and velocity magnitudes of the diffusive field, (2) rate of viscous dissipation , (3) vertical stability as characterized by the relative adiabatic lapse rate (1/T)(g/C _p +T/z), and (4) vertical shear in the mean horizontal winds , and , for a given height and of spatial extent equal to that of the diffusing cloud. Sample results for near ground level and for upper stratospheric heights are given. For the atmospheric boundary layer case, the diffusive field is microscale turbulence. In the upper stratospheric case it is considered to be a field of highly interactive and dispersive gravity waves.     

新加坡地应力测量

2008年和2009年, 通过和新加坡宜康地质技术有限责任公司合作, 在新加坡18个浅钻孔中进行了水压致裂地应力测量. 通过测量获得了新加坡地壳浅部的现今应力状态, 为新加坡拟建地下管缆隧道设计施工提供了地应力依据, 也为世界应力环境数据库提供了参考. 测试结果表明, 在38.00~56.95 m深度域内最大水平主应力σ_h值为3.45~9.64 MPa, 最小水平主应力σ_h为2.49~6.28 MPa. 水平主应力和用上覆岩层重量计算的垂直应力σv之间的关系为σ_h>σ_h>σv,最大、最小水平主应力平均值分别是垂直应力平均值的5.13倍和3.23倍, 表明新加坡应力场以水平应力为主导. 初步分析认为, 新加坡陆域浅部构造应力场的量值大致呈南、北两端偏低, 中间高的态势, 在北端最大与最小水平主应力值都小于它们各自的平均值, 而南端则与平均值相当. 最大水平主应力优势方向为NE—NEE向, 该方向与前人在其邻区获得的构造应力场方向和由震源机制解反映的主压应力方向较为一致. 新加坡现今构造应力场状态与印度板块、澳大利亚板块、菲律宾板块与欧亚板块的汇聚作用有着密切的联系. 根据库仑摩擦滑动准则, 断层面摩擦系数取0.6~1.0, 分析了研究区内断层产生摩擦滑动的可能性.     

An investigation of the dependence of structural parameters of the middle atmosphere on the intensity of photochemical and dynamical processes

Using a numerical model for temperature and neutral and ion composition behaviour at middle atmospheric heights, an analysis has been made of the dependence of atmospheric structural parameters on temperature, solar activity, and on turbulent transfer intensity. For mesospheric heights, an inverse dependence of the nitric oxide density on the temperature has been found. It is thus possible to explain experimentally obtained temperature variations over a cyclc of solar activity at mesospheric and lower thermospheric heights. Numerical simulation results indicate that the temperature in the height range 75–120 km depends considerably on both the absolute values of turbulent transfer coefficients and their vertical gradients.     

Physical mechanism of strong negative storm effects in the daytime ionospheric F2 region observed with EISCAT

A self-consistent method for daytime F-region modelling was applied to EISCAT observations during two periods comprising the very disturbed days 3 April 1992 and 10 April 1990. The observed strong N_e decrease at F2-layer heights originated from different physical mechanisms in the two cases. The negative F2-layer storm effect with an N_mF2 decrease by a factor of 6.4 on 3 April 1992 was produced by enhanced electric fields (E 85 mV/m) and strong downward plasma drifts, but without any noticeable changes in thermos-pheric parameters. The increase of the O⁺ + N₂ reaction rate resulted in a strong enrichment of the ionosphere with molecular ions even at F2-layer heights. The enhanced electric field produced a wide mid-latitude daytime trough on 03 April 1992 not usually observed during similar polarization jet events. The other strong negative storm effect on 10 April 1990 with a complete disappearance of the F2-layer maximum at the usual heights was attributed mainly to changes in neutral composition and temperature. A small value for the shape parameter S in the neutral temperature profile and a low neutral temperature at 120 km indicate strong cooling of the lower thermosphere. We propose that this cooling is due to increased nitric oxide concentration usually observed at these heights during geomagnetic storms.     

Effects of different gap filling methods and land surface energy balance closure on annual net ecosystem exchange in a semiarid area of China

Based on eddy covariance measurements over two kinds of land surfaces(a degraded grassland and a maize cropland)in a semiarid area of China in 2005 and 2008,the effects of different gap filling methods,energy balance closure and friction velocity threshold(u*)on annual net ecosystem exchange(NEE)were analyzed.Six gap filling methods,including mean diurnal variation(MDV),marginal distribution sampling(MDS),and nonlinear regressions method,were investigated by generating secondary datasets with four different artificial gap lengths(ranging in length from single half-hours to 12 consecutive days).The MDS generally showed a good overall performance especially for long gaps,with an annual sum bias error less than 5 g C m-2 yr-1.There was a large positive annual sum bias error for nonlinear regressions,indicating an overestimate on net ecosystem respiration.The offset in the annual sum NEE for four nonlinear regressions was from 8.0 to 30.8 g C m-2 yr-1.As soil water content was a limiting factor in the semiarid area,the nonlinear regressions considering both soil temperature and soil water content as controlling variables had a better performance than others.The performance of MDV was better in daytime than in nighttime,with an annual sum bias error falling between-2.6 and-13.4 g C m-2 yr-1.Overall,the accuracy of the gap filling method was dependent on the type of the land surface,gap length,and the time of day when the data gap occurred.The energy balance ratio for the two ecosystems was nearly 80%.Turbulent intensity had a large impact on energy balance ratio.Low energy balance ratio was observed under low friction velocity during the night.When there was a large fetch distance in a wind direction,a low energy balance ratio was caused by mismatch of the footprints between the available energy and turbulent fluxes.The effect of energy balance correction on CO2 flux was evaluated by assuming the imbalance caused by the underestimation of sensible heat flux and latent heat flux.The results showed an average increase of 10 g C m-2 yr-1 for annual NEE in both ecosystems with an energy balance correction.On the other hand,the u*threshold also have a large impact on annual sum NEE.Net carbon emission increased 37.5 g C m-2 yr-1 as u*threshold increased from 0.1 to 0.2 m s-1,indicating a large impact of imposing u*threshold on net ecosystem carbon exchange.     

插值公式、相关系数和采样间隔对GRACE Follow-On 星间加速度精度的影响

本文基于星间加速度法开展了插值公式、相关系数和采样间隔对GRACE Follow-On星间加速度精度影响的研究. 模拟结果表明:1)适当增加数值微分公式的插值点数可有效提高插值精度. 基于9点Newton插值公式,星间加速度的插值误差为4.401×10^-13 m·s^-2,分别基于7点、5点和3点插值公式,插值误差增加了1.192倍、6.912倍和274.029倍. 2)适当增大相关系数可有效降低星间加速度的误差. 基于相关系数0.99,星间加速度方差为3.777×10^-24 m²·s^-4,分别基于相关系数0.90、0.70、0.50和0.00,方差增加了9.780倍、22.404倍、26.217倍和26.820倍. 3)随着采样间隔增大,星间加速度方差逐渐降低,但卫星观测值的空间分辨率也同时降低,因此合理选取采样间隔有利于地球重力场精度的提高. 4)基于9点Newton插值公式、相关系数(K波段测量系统星间距离和星间速度0.85、GPS轨道位置和轨道速度0.95、星载加速度计非保守力0.90)和采样间隔10 s,利用预处理共轭梯度迭代法,精确和快速反演了120阶GRACE Follow-On地球重力场,在120阶处累计大地水准面精度为4.602×10^-4 m.     

The α-Effect and Current Helicity for Fast Sheared Rotators

We explore the f -effect and the small-scale current helicity, , for the case of weakly compressible magnetically driven turbulence that is subjected to the differential rotation. No restriction is applied to the amplitude of angular velocity, i.e., the derivations presented are valid for an arbitrary Coriolis number, z * = 2 z cor , though the differential rotation itself is assumed to be weak. The expressions obtained are used to explore the possible distributions of f -effect and h c in convection zones (CZ) of the solar-type stars. Generally, our theory gives f { { > 0 in the northern hemisphere of the Sun and the opposite case in the southern hemisphere. In most cases the h c has the opposite sign to f { { . However, we show that in the depth of CZ where the influence of rotation upon turbulence (associated with z *) and the radial shear of angular velocity are strong, the distribution of f { { might be drastically different from a classical cos è -dependence, where è is colatitude. It is shown that f { { has a negative sign at the bottom and below of CZ at mid latitudes. There, the distribution of h c is also different from cos è , but it does not change its sign with the depth. Further, we briefly consider these quantities in the disk geometry. The application of the developed theory to dynamos in the accretion disk is more restrictive because they usually have a strong differential rotation, | ‘ log z / ‘ log r | > 1.     

西昌昔格达组黏土岩动力特性试验研究

在大量动三轴试验的基础上,研究了不同固结条件(σ3c=100,200,400 kPa,kc=1.0,1.5,2.0)下西昌昔格达组黏土岩的动应力-动应变关系、动弹性模量、阻尼比及动强度特性。研究结果表明:西昌昔格达组黏土岩在一定条件下的动应力-动应变关系符合幂函数模型。动应力随着固结围压σ3c或固结主压力比kc增大而增大。在相同围压下,达到相同动应变时,偏压固结状态比均压状态所需要的动应力大。当其他条件相同时,动弹性模量随固结压力或固结主压力比增大而增大,随着动应变的增加而减小并趋于稳定。最大动弹性模量Ed0与σ3c/Pa、kc均有良好的幂函数关系,且对不同固结应力状态条件可以归一。阻尼比随动应变增大而增大,但增大幅度随动应变增大而迅速减小,阻尼比随围压或固结主压力比增大而减小,但它们对阻尼比影响较小。不同固结应力状态下(λ/λmax)与(1-Ed/Ed0)有良好的幂函数关系。西昌昔格达组饱和黏土岩动剪应力τd随固结围压σ3c或固结主压力比kc增大而增大,随振动次数增加而减小。τd/σm均随固结围压σ3c增大而减小,随固结压力比kc增大而增大。动剪应力τd与平均压力σm、固结主压力比kc有较好的线性关系。     

Airborne microorganisms: Their occurrence,distribution and dependence on environmental factors — Especially on organic compounds of air-pollution

For a period of approximately six months the concentrations of airborne microorganisms were determined by the uptake of their day-characteristic at local places differing by the extent of air pollution. By simultaneously recording the meteorological parameters as temperature, relative humidity, direction and velocity of wind and intensity of solar radiation, it was intended to find out the relationship between these factors and the viability of airborne microorganisms. The present data show an oscillation of the concentration counts over a wide range and seem not yet to depend on the local situation of air pollution.In a further series of experiments the atmospheric aerosol <10 m was collected on air-filters. The decreased substances were extracted from the filter material by aqueous and organic solutants and tested for their antimicrobial activities. In this context the seasonal dependence of the microbicidal activity of organic compounds of aerosol is discussed.     

Hydromagnetic spin-up under conditions of the earth's core-mantle boundary

Summary The disturbances of the velocity and magnetic fields close to the Earth's core-mantle boundary, caused by sudden irregular changes in the Earth's rotation, are investigated. The problem leads to the investigation of the structure of the Ekman-Hartman hydromagnetic boundary layer, the magnetic diffusive region and the currentless region. Precise Laplacean inversions of the images of all disturbances in the Earth's core-mantle system are obtained for the limiting case of a zero magnetic Prandtl number, =0. The disturbance of the velocity in the direction of the axis of rotation (Ekman suction) in the currentless region has the nature of inertial oscillations with a frequency of 2. Additional disturbances (with respect to the case of =0) of the velocity in the azimuthal and radial directions, particularly for the EHL and MDR region, are determined for the case when 0< « 1. The disturbance in the velocity again has the character of inertial oscillations with the frequency 2, being exponentially damped in EHL asexp (–2²t) and in MDR asexp (–²t).     

基于GRACE卫星重力数据确定地球重力场模型WHU-GM-05

基于卫星轨道运动的能量积分方程,可导出利用卫星跟踪卫星数据求解地球重力场的实用公式.本文在Jekeli给出的公式基础上导出了基于能量守恒方程利用两颗低-低卫星跟踪的扰动位差求解重力位系数的严密关系式.基于两颗GRACE卫星的观测数据,采用本文导出的严密能量积分方法求解得到120阶的GRACE地球重力场模型,命名为WHU-GM-05;将WHU-GM-05模型与国际上同类重力场模型EIGEN-GRACE系列和GGM02S分别在阶方差和大地水准面高等方面作了比较,并与美国和中国的部分地区GPS水准观测值进行了精度分析.结果表明基于本文推导的严密双星能量守恒方程得到的WHU-GM-05重力场模型精度与国际上同类重力场模型的精度相当.     

沙尘暴电效应的实验观测研究

利用国内大型风沙物理风洞实验模拟沙尘暴电现象,研究风沙起电机理,结果表 明,不同风速下不同沙粒会产生不同极性的电场强度和电位效应,风沙电随风速增大而增强 ,且随沙粒度增大而减小. 在沙漠区的16m,8m,4m和1m高度上观测到27次不同沙尘暴天气 过程的电场和风速随时间变化. 结果表明, 在晴天4个高度上的电场均为小正电场值,电场 随高度降低而减小,最大电场强度在5kV/m以下,日风速变化对各层电场起伏没有较大影响 . 有沙尘天气,各高度上的电场强度随风速变化而变化. 16m高度上电场均为负值,平均值 为 -20kV/m;中层8m 电场一般为较高正电场值,达到10～40kV/m,与16m高度上电场呈反相 关;下层1m 电场值变化一般很小,在1kV/m以下. 在强沙尘暴天气4个高度上的电场均为负 电值,电场值随高度降低而减小,16m高度上最大平均电场强度达到-200kV/m以上,瞬时值 超过 -2500kV/m,与晴天电场矢量相反.     

白云母对岩盐断层带摩擦速度依赖性影响的实验研究

为更好地理解层状硅酸盐对断层强度、滑动速度依赖性及地震活动特征的影响,利用双轴摩擦实验对含白云母岩盐断层带在干燥及含水条件下摩擦的速度依赖性进行了实验研究,并观测了摩擦滑动过程中的声发射,分析了断层带的微观结构.实验结果表明,干燥条件下含白云母岩盐断层带在0.1 ～ 100μm/s的速度范围内表现为黏滑和速度弱化,增大σ2会使断层带从速度弱化向速度强化转化,速度依赖性转换出现在0.1 μm/s,其中断层滑动表现为稳滑或应力释放时间较长的黏滑事件;含水条件下含白云母岩盐断层带在0.05 ～0.01μm/s的速度范围内表现为速度强化,0.1 ～10μm/s的速度范围内表现为速度弱化,50～100μm/s的速度范围内又转换为速度强化行为.含白云母岩盐断层带在干燥条件下一次黏滑伴随一个或一丛声发射事件,而在含水条件下与稳滑相对应,滑动过程中并未记录到声发射事件.显微结构观察表明,速度弱化域的主要变形机制是岩盐颗粒的脆性破裂和局部化的滑动;干燥条件下,速度强化域的主要变形机制是岩盐颗粒的均匀破裂;含水条件下2个速度强化域对应不同的微观机制,高速域的速度强化受控于岩盐颗粒在白云母相互连结形成的网状结构上的滑动及其均匀碎裂作用,而低速域的速度强化还受岩盐的压溶作用控制.通过与岩盐断层带摩擦实验结果对比可知,白云母的存在对于燥岩盐断层带摩擦滑动方式和速度依赖性没有显著影响,而在含水条件下白云母的存在使得岩盐断层带滑动趋于稳定.实验结果为分析含层状硅酸盐断层的强度和稳定性提供了依据.此外,在速度依赖性转换域上观察到的应力缓慢释放的现象进一步证实了在岩盐断层带摩擦滑动过程中观察到的现象,这对慢地震机制研究具有参考意义.